We use data from the Evolution and Assembly of GaLaxies in their Environment (EAGLE) cosmological simulation to study properties of galaxies in the cosmic web. Galaxies become more redder and form stars at a lower rate relative to their counterparts further away from the cylindrical axis of the large-scale filaments. These trends are particularly strong for galaxies with $M_*/M_{\odot}\lesssim10^{10}$. We also find that at distances $<0.5$ Mpc from the spine of the filaments, the median gas and stellar mass fraction in filament galaxies rises sharply with decreasing distance from the spine of the filament. These results, together with matching trends in the SFR/$M_*$ and the $g-r$ colour of filament galaxies suggest that (i) the intrafilamentary gas condenses into the filament galaxies thus fuelling star formation in them, and (ii) increased number density of galaxies closer to the central axis of the filament enhances the rate of gravitational interactions among filament galaxies closer to the spine.

We study evolution of perturbations in dark matter and dark energy for spherical collapse. We study Tachyon models of dark energy using the approach outlined in Rajvanshi and Bagla (2018). We work with models that are allowed by current observations. We find that as with Quintessence models allowed by observations, dark energy perturbations do not affect evolution of perturbations in dark matter in a significant manner. Perturbations in dark energy remain small for such models. In order to compare different classes of dark energy models, we use reconstruction of potentials and study dark matter and dark energy perturbations for same expansion history in the two models. We find that dark matter perturbations carry no imprint of the class of dark energy models for the same expansion history: this is significant in that we can work with any convenient model to study clustering of dark matter. We find that the evolution of dark energy perturbations carries an imprint of the class of models in that these grow differently in Tachyon models and Quintessence models. However, the difference between these diminishes for (1+w)<<1 and hence prospects for differentiating between models using characteristics of perturbations are limited in our Universe.

We reported results of our study on non-linear spherical perturbations in quintessence models of dark energy. In the process of some follow up studies we discovered that a scaling factor in the code used for numerical calculations that should have been set to unity was set to a large value. Thus the scale of perturbations was much larger than intended, and for the larger scales the amplitude of dark matter perturbations was much higher than realistic. We provide corrected results here in this erratum. We find that there is no change in the perturbations for dark matter. The amplitude of perturbations in dark energy is much smaller than presented in the paper. Same holds true for spatial variation in the equation of state parameter.

The number of strong lens systems is expected to increase significantly in ongoing and upcoming surveys. With an increase in the total number of such systems we expect to discover many configurations that correspond to unstable caustics. In such cases, the instability can be used to our advantage for constraining the lens model. We have implemented algorithms for detection of different types of singularities in gravitational lensing. We test our approach on a variety of lens models and then go on to apply it to the inferred mass distribution for Abell 697 as an example application. We propose to represent lenses using A3-lines and singular points (A4 and D4) in the image plane. We propose this as a compact representation of complex lens systems that can capture all the details in a single snapshot.

Recent numerical simulations have suggested the probability of a fraction of the primordial stars to be ejected from the cluster of their origin. We explore the possibility that some of these can remain on the main sequence till the present epoch. We develop a semi-analytical model guided by results of cosmological simulations to study the mass accretion by these protostars as a function of the original stellar mass, and other parameters such as angular momentum and gravitational drag due to ambient gas. We also explore whether some of the protostars remain sufficiently low-mass and long-lived to survive to the present day. This requires that the protostars are ejected from the star forming region while their mass is less than $0.8\,M_{\odot}$. Assuming the protostars gain mass via the spherical Bondi-Hoyle accretion from the ambient medium, we show that Pop III protostars that initially form within a certain range of mass and are ejected with velocity larger than the escape velocity may survive to the present day on the main sequence. Thus, they may even be found in our Milky Way or its satellites. Our calculations also reveal that protostars that do not get ejected from the parent gas cloud accrete large amount of gas. We predict that these can become massive enough to be progenitors of black holes.

Recent numerical simulations have suggested the probability of a fraction of the primordial stars to be ejected from the cluster of their origin. We explore the possibility that some of these can remain on the main sequence till the present epoch. We develop a semi-analytical model guided by results of cosmological simulations to study the mass accretion by these protostars as a function of the original stellar mass, and other parameters such as angular momentum and gravitational drag due to ambient gas. We also explore whether some of the protostars remain sufficiently low-mass and long-lived to survive to the present day. This requires that the protostars are ejected from the star forming region while their mass is less than $0.8\,M_{\odot}$. Assuming the protostars gain mass via the spherical Bondi-Hoyle accretion from the ambient medium, we show that Pop III protostars that initially form within a certain range of mass and are ejected with velocity larger than the escape velocity may survive to the present day on the main sequence. Thus, they may even be found in our Milky Way or its satellites. Our calculations also reveal that protostars that do not get ejected from the parent gas cloud accrete large amount of gas. We predict that these can become massive enough to be progenitors of black holes.

We discuss the gravitational lensing of gravitational wave signals from coalescing binaries. We delineate the regime where wave effects are significant from the regime where geometric limit can be used. Further, we focus on the effect of micro-lensing and the combined effect of strong lensing and micro-lensing. We find that micro-lensing combined with strong lensing can introduce time varying phase shift in the signal and hence can lead to detectable differences in the signal observed for different images produced by strong lensing. This, coupled with the coarse localization of signal source in the sky for gravitational wave detections, can make it difficult to identify the common origin of signal corresponding to different images and use observables like time delay. In case we can reliably identify corresponding images, micro-lensing of individual images can be used as a tool to constrain properties of micro-lenses. Sources of gravitational waves can undergo microlensing due to lenses in the disk/halo of the Galaxy, or due to lenses in an intervening galaxy even in absence of strong lensing. In general the probability for this is small with one exception: Extragalactic sources of gravitational waves that lie in the galactic plane are highly likely to be micro-lensed. Wave effects are extremely important for such cases. In case of detections of such sources with low SNR, the uncertainty of occurrence of microlensing or otherwise introduces an additional uncertainty in the parameters of the source.

Dynamical models for dark energy are an alternative to the cosmological constant. It is important to investigate properties of perturbations in these models and go beyond the smooth FRLW cosmology. This allows us to distinguish different dark energy models with the same expansion history. For this, one often needs the potential for a particular expansion history. We study how such potentials can be reconstructed obtaining closed formulae for potential or reducing the problem to quadrature. We consider three classes of models here: tachyons, quintessence, and, interacting dark energy. We present results for constant w and the CPL parameterization. The method given here can be generalized to any arbitrary form of w(z).

We have used the upgraded Giant Metrewave Radio Telescope to carry out a deep (117 on-source hours) L-band observation of the Extended Groth Strip, to measure the average neutral hydrogen (HI) mass and median star formation rate (SFR) of star-forming galaxies, as well as the cosmic HI mass density, at $0.2 < z < 0.4$. This was done by stacking the HI 21cm emission and the rest-frame 1.4 GHz radio continuum from 445 blue star-forming galaxies with $\rm M_B \leq -17$ at $z_{\rm mean} \approx 0.34$. The stacked HI 21cm emission signal is detected at $\approx 7\sigma$ significance, implying an average HI mass of $\rm \langle M_{HI} \rangle = (4.93 \pm 0.70) \times 10^9 \: M_{\odot}$. We also stacked the rest-frame 1.4 GHz radio continuum emission of the same galaxies, to obtain a median SFR of $(0.54 \pm 0.06) \: {\rm M}_\odot$ yr$^{-1}$; this implies an average atomic gas depletion time scale of $\rm \langle \Delta t_{HI}\rangle \approx$ 9 Gyr, consistent with values in star-forming galaxies in the local Universe. This indicates that the star-formation efficiency does not change significantly over the redshift range $0 - 0.4$. We used the detection of the stacked HI 21cm emission signal to infer the normalized cosmic HI mass density $(\rm \rho_{HI}/\rho_{c,0})$ in star-forming galaxies at $z \approx 0.34$. Assuming the local relation between HI mass and absolute B-magnitude, we obtain $\rm \rho_{HI}/\rho_{c,0} = (4.81 \pm 0.75) \times 10^{-4}$, implying no significant evolution in $\rm \rho_{HI}/\rho_{c,0}$ from $z \approx 0.4$ to the present epoch.

We use an analytical approach to study ram pressure stripping with simple models for discs and halo gas distribution to study the phenomena in cluster, group and galaxy halos. We also study variations with galaxy properties and redshift. In each case we model the worst case scenario (i.e., maximum effect due to ram pressure). We show that there is little variation in the worst case scenario with redshift. We find that gas discs in galaxies with a higher spin parameter get stripped sooner than galaxies with a smaller spin parameter. Galaxies in cluster halos get stripped of gas more effciently as compared to group and galaxy halos: this is due to the higher infall speed and a higher density of gas in the ICM due to a greater retention of baryons. We comment on the limitations of our model and situations where a signiffcant amount of gas may be retained in galaxy disc and also give an illustration for the same. Lastly, we discuss implications for star formation in galaxies as these fall into halos.

We study the angular clustering of point sources in The GMRT (Giant Meter Wave Telescope) Sky Survey (TGSS). The survey at 150 MHz with delta > -53.5 degrees has a sky coverage of 3.6 pi steradians, i.e., 90% of the whole sky. We created subsamples by applying different total flux thresholds limit (S >> 5 sigma) for good completeness and measured the angular correlation function omega(theta) of point sources at large scales ( >= 1 degree). We find that the amplitude of angular clustering is higher for brighter subsamples, this indicates that higher threshold flux samples are hosted by massive halos and cluster strongly: this conclusions is based on the assumption that the redshift distribution of sources does not change with flux and this is supported by models of radio sources. We compare our results with other low-frequency studies of clustering of point sources and verify that the amplitude of clustering varies with the flux limit. We quantify this variation as a power law dependence of the amplitude of correlation function with the flux limit. This dependence can be used to estimate foreground contamination due to clustering of point sources for low frequency HI intensity mapping surveys for studying the epoch of reionisation.

Observations show that the expansion of the Universe is accelerating. This requires that the dominant constituent of matter in the Universe has some unusual properties like negative pressure. This exotic component has been given the name dark energy. We work with the simplest model of dark energy, the cosmological constant introduced by Einstein. We study the evolution of spherical over-densities in such a model and show that there is a minimum over-density required for collapse: perturbations with a smaller amplitude do not collapse. This threshold is interesting as even perturbations with a positive over-density and negative energy do not collapse in finite time. Further, we show that perturbations with an amplitude larger than, but comparable to the threshold value, take a very long time to collapse. We compare the solutions with the case when dark energy is absent.

Diffuse Galactic emission at low frequencies is a major contaminant for studies of redshifted $21$ cm line studies. Removal of these foregrounds is essential for exploiting the signal from neutral hydrogen at high redshifts. Analysis of foregrounds and its characteristics is thus of utmost importance. It is customary to test efficacy of foreground removal techniques using simulated foregrounds. Most simulations assume that the distribution of the foreground signal is a Gaussian random field. In this work we test this assumption by computing the binned bispectrum for the all-sky $408$ MHz map. This is done by applying different brightness temperature ($T$) thresholds in order to assess whether the cooler parts of the sky have different characteristics. We find that regions with a low brightness temperature $T < 25$ K indeed have smaller departures from a Gaussian distribution. Therefore, these regions of the sky are ideal for future H{\sc i} intensity mapping surveys.

Observations have confirmed the accelerated expansion of the universe. The accelerated expansion can be modelled by invoking a cosmological constant or a dynamical model of dark energy. A key difference between these models is that the equation of state parameter $w$ for dark energy differs from $-1$ in dynamical dark energy (DDE) models. Further, the equation of state parameter is not constant for a general DDE model. Such differences can be probed using the variation of scale factor with time by measuring distances. Another significant difference between the cosmological constant and DDE models is that the latter must cluster. Linear perturbation analysis indicates that perturbations in quintessence models of dark energy do not grow to have a significant amplitude at small length scales. In this paper we study the response of quintessence dark energy to non-linear perturbations in dark matter. We use a fully relativistic model for spherically symmetric perturbations. In this study we focus on thawing models. We find that in response to non-linear perturbations in dark matter, dark energy perturbations grow at a faster rate than expected in linear perturbation theory. We find that dark energy perturbation remains localised and does not diffuse out to larger scales. The dominant drivers of the evolution of dark energy perturbations are the local Hubble flow and a supression of gradients of the scalar field. We also find that the equation of state parameter $w$ changes in response to perturbations in dark matter such that it also becomes a function of position. The variation of $w$ in space is correlated with density contrast for matter. Variation of $w$ and perturbations in dark energy are more pronounced in response to large scale perturbations in matter while the dependence on the amplitude of matter perturbations is much weaker.

Distribution of cold gas in the post-reionization era provides an important link between distribution of galaxies and the process of star formation. Redshifted 21 cm radiation from the Hyperfine transition of neutral Hydrogen allows us to probe the neutral component of cold gas, most of which is to be found in the interstellar medium of galaxies. Existing and upcoming radio telescopes can probe the large scale distribution of neutral Hydrogen via HI intensity mapping. In this paper we use an estimate of the HI power spectrum derived using an ansatz to compute the expected signal from the large scale HI distribution at z ~ 3. We find that the scale dependence of bias at small scales makes a significant difference to the expected signal even at large angular scales. We compare the predicted signal strength with the sensitivity of radio telescopes that can observe such radiation and calculate the observation time required for detecting neutral Hydrogen at these redshifts. We find that OWFA (Ooty Wide Field Array) offers the best possibility to detect neutral Hydrogen at z ~ 3 before the SKA (Square Kilometer Array) becomes operational. We find that the OWFA should be able to make a 3 sigma or a more significant detection in 2000 hours of observations at several angular scales. Calculations done using the Fisher matrix approach indicate that a 5 sigma detection of the binned HI power spectrum via measurement of the amplitude of the HI power spectrum is possible in 1000 hours (Sarkar, Bharadwaj and Ali, 2017).

We investigate the relative importance of two galactic outflow suppression mechanisms : a) Cosmological infall of the intergalactic gas onto the galaxy, and b) the existence of a hot circumgalactic medium (CGM). Considering only radial motion, the infall reduces the speed of outflowing gas and even halts the outflow, depending on the mass and redshift of the galaxy. For star forming galaxies there exists an upper mass limit beyond which outflows are suppressed by the gravitational field of the galaxy. We find that infall can reduce this upper mass limit approximately by a factor of two (independent of the redshift). Massive galaxies ($\gtrsim 10^{12} M_{\odot}$) host large reservoir of hot, diffuse CGM around the central part of the galaxy. The CGM acts as a barrier between the infalling and outflowing gas and provides an additional source of outflow suppression. We find that at low redshifts ($z\lesssim3.5$), the CGM is more effective than the infall in suppressing the outflows. Together, these two processes give a mass range in which galaxies are unable to have effective outflows. We also discuss the impact of outflow suppression on the enrichment history of the galaxy and its environment.

We derive an estimate of the rate of formation of dark matter halos per unit volume as a function of the halo mass and redshift of formation. Analytical estimates of the number density of dark matter halos are useful in modeling several cosmological phenomena. We use the excursion set formalism for computing the formation rate of dark matter halos. We use an approach that allows us to differentiate between major and minor mergers, as this is a pertinent issue for semi-analytic models of galaxy formation. We compute the formation rate for the Press-Schechter and the Sheth-Tormen mass function. We show that the formation rate computed in this manner is positive at all scales. We comment on the Sasaki formalism where negative halo formation rates are obtained. Our estimates compare very well with N-Body simulations for a variety of models. We also discuss the halo survival probability and the formation redshift distributions using our method.

We model the distribution of neutral Hydrogen (HI hereafter) in the post-reionization universe. This model uses gravity only N-Body simulations and an ansatz to assign HI to dark matter haloes that is consistent with observational constraints and theoretical models. We resolve the smallest haloes that are likely to host HI in the simulations, care is also taken to ensure that any errors due to the finite size of the simulation box are small. We then compute the smoothed one point probability distribution function and the power spectrum of fluctuations in HI. This is compared with other predictions that have been made using different techniques. We highlight the significantly high bias for the HI distribution at small scales. This aspect has not been discussed before. We then discuss the prospects for detection with the MWA, GMRT and the hypothetical MWA5000. The MWA5000 can detect visibility correlations at large angular scales at all redshifts in the post-reionization era. The GMRT can detect visibility correlations at lower redshifts, specifically there is a strong case for a survey at z=1.3. We also discuss prospects for direct detection of rare peaks in the HI distribution using the GMRT. We show that direct detection should be possible with an integration time that is comparable to, or even less than, the time required for a statistical detection. Specifically, it is possible to make a statistical detection of the HI distribution by measuring the visibility correlation, and, direct detection of rare peaks in the HI distribution at z = 1.3 with the GMRT in less than 1000 hours of observations.

Simulated maps of the HI distribution in the post-reionization era are used to study the prospects for detection with existing and upcoming radio telescopes. We consider detection in the redshifted radiation from the hyperfine transition with a rest frame frequency of 1420 MHz. Possibility of a statistical detection using visibility correlations is discussed. We show that the MWA (Murchison Widefield Array) and the GMRT (Giant Meterwave Radio Telescope) can potentially detect signal from the HI distribution at high redshifts. MWA can detect visibility correlations at large angular scales at all redshifts accessible to it in the post-reionization era. The GMRT can detect visibility correlations at lower redshifts, specifically there is a strong case for a survey at z = 1.3. We also discuss prospects for direct detection of rare peaks in the HI distribution using the GMRT. We show that direct detection should be possible with integration time that is comparable to, or even less than, the time required for a statistical detection. Specifically, it is possible to make a statistical detection of the HI distribution by measuring the visibility correlation, and, direct detection of rare peaks in the HI distribution using the GMRT in less than 1000 hours of observations.

In the multi-fractal analysis of large scale matter distribution, the scale of transition to homogeneity is defined as the scale above which the fractal dimension of underlying point distribution is equal to the ambient dimension of the space in which points are distributed. With finite sized weakly clustered distribution of tracers obtained from galaxy redshift surveys it is difficult to achieve this equality. Recently we have defined the scale of homogeneity to be the scale above which the deviation of fractal dimension from the ambient dimension becomes smaller than the statistical dispersion. In this paper we use the relation between the fractal dimensions and the correlation function to compute the dispersion for any given model in the limit of weak clustering amplitude. We compare the deviation and dispersion for the LCDM model and discuss the implication of this comparison for the expected scale of homogeneity in the concordant model of cosmology. We estimate the upper limit to the scale of homogeneity to be close to 260 Mpc/h for the LCDM model. Actual estimates of the scale of homogeneity should be smaller than this as we have considered only statistical contribution to the dispersion in fractal dimension and we have ignored cosmic variance and contributions due to survey geometry and the selection function. We find that as long as non linear correction are insignificant, scale of homogeneity as defined above does not change with epoch. The scale of homogeneity depends very weakly on the choice of tracer of the density field. Thus the suggested definition of the scale of homogeneity is fairly robust.

We study the observational constraints of CMB temperature and polarization anisotropies on models of dark energy, with special focus on models with variation in properties of dark energy with time. We demonstrate that the key constraint from CMB observations arises from the location of acoustic peaks. An additional constraint arises from the limits on Omega_nr from the relative amplitudes of acoustic peaks. Further, we show that the distance to the last scattering surface is not how the CMB observations constrain the combination of parameters for models of dark energy. We also use constraints from Supernova observations and show that unlike the Gold and Silver samples, the SNLS sample prefers a region of parameter space that has a significant overlap with the region preferred by the CMB observations. This is a verification of a conjecture made by us in an earlier work. We discuss combined constraints from WMAP5 and SNLS observations. We find that models with w ~ -1 are preferred for models with a constant equation of state parameters. In case of models with a time varying dark energy, we show that constraints on evolution of dark energy density are almost independent of the type of variation assumed for the equation of state parameter. This makes it easy to get approximate constraints from CMB observations on arbitrary models of dark energy. Constraints on models with a time varying dark energy are predominantly due to CMB observations, with Supernova constraints playing only a marginal role.

Press-Schechter theory gives a simple, approximate functional form of the mass function of dark matter haloes. Sheth and Tormen (ST) refined this mass function to give an improved analytical fit to results of N-body simulations. These forms of the halo mass function are universal (independent of cosmology and power spectrum) when scaled in suitable variables. Using large suites of LCDM N-body simulations, studies in the last few years have shown that this universality is only approximate. We explore whether some of the deviations from universality can be attributed to the power spectrum by computing the mass function in N-body simulations of various scale-free models in an Einstein-de Sitter cosmology. This choice of cosmology does not introduce any scale into the problem. These models have the advantage of being self-similar, hence stringent checks can be imposed while running these simulations. This set of numerical experiments is designed to isolate any power spectrum dependent departures from universality of mass functions. We show explicitly that the best fit ST parameters have a clear dependence on power spectrum. Our results also indicate that an improved analytical theory with more parameters is required in order to provide better fits to the mass function.

Cosmological N-Body simulations are used for a variety of applications. Indeed progress in the study of large scale structures and galaxy formation would have been very limited without this tool. For nearly twenty years the limitations imposed by computing power forced simulators to ignore some of the basic requirements for modeling gravitational instability. One of the limitations of most cosmological codes has been the use of a force softening length that is much smaller than the typical inter-particle separation. This leads to departures from collisionless evolution that is desired in these simulations. We propose a particle based method with an adaptive resolution where the force softening length is reduced in high density regions while ensuring that it remains well above the local inter-particle separation. The method, called the Adaptive TreePM, is based on the TreePM code. We present the mathematical model and an implementation of this code, and demonstrate that the results converge over a range of options for parameters introduced in generalizing the code from the TreePM code. We explicitly demonstrate collisionless evolution in collapse of an oblique plane wave. We compare the code with the fixed resolution TreePM code and also an implementation that mimics adaptive mesh refinement methods and comment on the agreement, and disagreements in the results. We find that in most respects the ATreePM code performs at least as well as the fixed resolution TreePM in highly over-dense regions, from clustering and number density of haloes, to internal dynamics of haloes. We also show that the adaptive code is faster than the corresponding high resolution TreePM code.

We explore the prospects of using the hyperfine transition of 3He+ as a probe of the intergalactic medium. The emission signal from ionized regions during reionization is expected to be anti-correlated with 21cm maps. The predicted emission signal from Lyman-alpha blobs at lower redshifts is detectable with future radio observatories.

The epoch of reionization and formation of first stars are inter-linked topics that are of considerable interest. We use a simplified approach for studying formation of stars in collapsed haloes, and the resulting ionization of the inter-galactic medium (IGM). We consider a set of LCDM models allowed by observations of CMB temperature and polarization anisotropies for this study. We constrain parameters related to star formation with the help of observations. We constrain subsets of these parameters independently by using the observed metallicity of the inter-galactic medium at z ~ 5 and the requirement that the Thomson scattering optical depth due to an ionized IGM as determined for the model from CMB observations be reproduced. We consider a range of initial metalicities for star forming gas, and some variations of the initial mass function of stars. We find that a "normal" initial mass function (IMF) may satisfy these two constraints with a raised efficiency of star formation as compared to that seen in the local universe. Observations require a significant fraction of metals to escape from haloes to the IGM. We can also place constraints on the ratio of escape fraction for metals and ionizing photons, we find that this ratio is of order unity for most models. This highlights the importance of using the constraints arising from enrichment of the inter-galactic medium. (Abridged)

We discuss the performance characteristics of using the modification of the tree code suggested by Barnes \citep{1990JCoPh..87..161B} in the context of the TreePM code. The optimisation involves identifying groups of particles and using only one tree walk to compute force for all the particles in the group. This modification has been in use in our implementation of the TreePM code for some time, and has also been used by others in codes that make use of tree structures. In this paper, we present the first detailed study of the performance characteristics of this optimisation. We show that the modification, if tuned properly can speed up the TreePM code by a significant amount. We also combine this modification with the use of individual time steps and indicate how to combine these two schemes in an optimal fashion. We find that the combination is at least a factor of two faster than the modified TreePM without individual time steps. Overall performance is often faster by a larger factor, as the scheme of groups optimises use of cache for large simulations.

N-Body simulations are an important tool in the study of formation of large scale structures. Much of the progress in understanding the physics of galaxy clustering and comparison with observations would not have been possible without N-Body simulations. Given the importance of this tool, it is essential to understand its limitations as ignoring these can easily lead to interesting but unreliable results. In this paper we study the limitations due to the finite size of the simulation volume. In an earlier work we proposed a formalism for estimating the effects of a finite box-size on physical quantities and applied it to estimate the effect on the amplitude of clustering, mass function. Here, we extend the same analysis and estimate the effect on skewness and kurtosis in the perturbative regime. We also test the analytical predictions from the earlier work as well as those presented in this paper. We find good agreement between the analytical models and simulations for the two point correlation function and skewness. We also discuss the effect of a finite box size on relative velocity statistics and find the effects for these quantities scale in a manner that retains the dependence on the averaged correlation function.

We study the interplay of clumping at small scales with the collapse and relaxation of perturbations at larger scales using N-Body simulations. We quantify the effect of collapsed haloes on perturbations at larger scales using two point correlation function, moments of counts in cells and mass function. The purpose of the study is twofold and the primary aim is to quantify the role played by collapsed low mass haloes in the evolution of perturbations at large scales, this is in view of the strong effect seen when the large scale perturbation is highly symmetric. Another reason for this study is to ask whether features or a cutoff in the initial power spectrum can be detected using measures of clustering at scales that are already non-linear. The final aim is to understand the effect of ignoring perturbations at scales smaller than the resolution of N-Body simulations. We find that these effects are ignorable if the scale of non-linearity is larger than the average inter-particle separation in simulations. Features in in the initial power spectrum can be detected easily if the scale of these features is in the linear regime, detecting such features becomes difficult as the relevant scales become non-linear. We find no effect of features in initial power spectra at small scales on the evolved power spectra at large scales. We may conclude that in general, the effect on evolution of perturbations at large scales of clumping on small scales is very small and may be ignored in most situations.

Homogeneity and isotropy of the universe at sufficiently large scales is a fundamental premise on which modern cosmology is based. Fractal dimensions of matter distribution is a parameter that can be used to test the hypothesis of homogeneity. In this method, galaxies are used as tracers of the distribution of matter and samples derived from various galaxy redshift surveys have been used to determine the scale of homogeneity in the Universe. Ideally, for homogeneity, the distribution should be a mono-fractal with the fractal dimension equal to the ambient dimension. While this ideal definition is true for infinitely large point sets, this may not be realised as in practice, we have only a finite point set. The correct benchmark for realistic data sets is a homogeneous distribution of a finite number of points and this should be used in place of the mathematically defined fractal dimension for infinite number of points (D) as a requirement for approach towards homogeneity. We derive the expected fractal dimension for a homogeneous distribution of a finite number of points. We show that for sufficiently large data sets the expected fractal dimension approaches D in absence of clustering. It is also important to take the weak, but non-zero amplitude of clustering at very large scales into account. In this paper we also compute the expected fractal dimension for a finite point set that is weakly clustered. Clustering introduces departures in the Fractal dimensions from D and in most situations the departures are small if the amplitude of clustering is small. Features in the two point correlation function, like those introduced by Baryon Acoustic Oscillations (BAO) can lead to non-trivial variations in the Fractal dimensions where the amplitude of clustering and deviations from D are no longer related in a monotonic manner.

N-Body simulations are a very important tool in the study of formation of large scale structures. Much of the progress in understanding the physics of galaxy formation and comparison with observations would not have been possible without N-Body simulations. Given the importance of this tool, it is essential to understand its limitations as ignoring these can easily lead to interesting but unreliable results. In this paper we study the limitations due to the finite size of the simulation volume. We explicitly construct the correction term arising due to a finite box size and study its generic features for clustering of matter and also on mass functions. We show that the correction to mass function is maximum near the scale of non-linearity, as a corollary we show that the correction to the number density of haloes of a given mass changes sign at this scale; the number of haloes at small masses is over estimated in simulations. This over estimate results from a delay in mergers that lead to formation of more massive haloes. The same technique can be used to study corrections to other physical quantities. The corrections are typically small if the scale of non-linearity is much smaller than the box-size. However, there are some cases of physical interest in which the relative correction term is of order unity even though a simulation box much larger than the scale of non-linearity is used. Within the context of the concordance model, our analysis suggests that it is very difficult for present day simulations to resolve mass scales smaller than 100 solar masses accurately and the level of difficulty increases as we go to even smaller masses, though this constraint does not apply to multi-scale simulations.

We study the non-Gaussian tail of the probability distribution function of density in cosmological N-Body simulations for a variety of initial conditions. We compare the behaviour of the non-Gaussian tail in the real space with that in the redshift space. The form of the PDF in redshift space is of great significance as galaxy surveys probe this and not the real space analogue predicted using theoretical models. We model the non-Gaussian tail using the halo model. In the weakly non-linear regime the moments of counts in cells in the redshift space approach the values expected from perturbation theory for moments in real space. We show that redshift space distortions in the non-linear regime dominate over signatures of initial conditions or the cosmological background. We illustrate this using Skewness and higher moments of counts in cells, as well as using the form of the non-Gaussian tail of the distribution function. We find that at scales smaller than the scale of non-linearity the differences in Skewness, etc. for different models are very small compared to the corresponding differences in real space. We show that bias also leads to smaller values of higher moments, but the redshift space distortions are typically the dominant effect.

The dark energy component of the universe is often interpreted either in terms of a cosmological constant or as a scalar field. A generic feature of the scalar field models is that the equation of state parameter w= P/rho for the dark energy need not satisfy w=-1 and, in general, it can be a function of time. Using the Markov chain Monte Carlo method we perform a critical analysis of the cosmological parameter space, allowing for a varying w. We use constraints on w(z) from the observations of high redshift supernovae (SN), the WMAP observations of CMB anisotropies and abundance of rich clusters of galaxies. For models with a constant w, the LCDM model is allowed with a probability of about 6% by the SN observations while it is allowed with a probability of 98.9% by WMAP observations. The LCDM model is allowed even within the context of models with variable w: WMAP observations allow it with a probability of 99.1% whereas SN data allows it with 23% probability. The SN data, on its own, favors phantom like equation of state (w<-1) and high values for Omega_NR. It does not distinguish between constant w (with w<-1) models and those with varying w(z) in a statistically significant manner. The SN data allows a very wide range for variation of dark energy density, e.g., a variation by factor ten in the dark energy density between z=0 and z=1 is allowed at 95% confidence level. WMAP observations provide a better constraint and the corresponding allowed variation is less than a factor of three. Allowing for variation in w has an impact on the values for other cosmological parameters in that the allowed range often becomes larger. (Abridged)

We investigate non-linear scaling relations for two-dimensional gravitational collapse in an expanding background using a 2D TreePM code and study the strongly non-linear regime ($\bar\xi \leq 200$) for power law models. Evolution of these models is found to be scale-invariant in all our simulations. We find that the stable clustering limit is not reached, but there is a model independent non-linear scaling relation in the asymptotic regime. This confirms results from an earlier study which only probed the mildly non-linear regime ($\bar\xi \leq 40$). The correlation function in the extremely non-linear regime is a less steep function of scale than reported in earlier studies. We show that this is due to coherent transverse motions in massive haloes. We also study density profiles and find that the scatter in the inner and outer slopes is large and that there is no single universal profile that fits all cases. We find that the difference in typical density profiles for different models is smaller than expected from similarity solutions for halo profiles and transverse motions induced by substructure are a likely reason for this difference being small.

We study the interplay of clumping at small scales with the collapse and relaxation of perturbations at much larger scales. We present results of our analysis when the large scale perturbation is modelled as a plane wave. We find that in absence of substructure, collapse leads to formation of a pancake with multi-stream regions. Dynamical relaxation of plane wave is faster in presence of substructure. Scattering of substructures and the resulting enhancement of transverse motions of haloes in the multi-stream region lead to a thinner pancake. In turn, collapse of the plane wave leads to formation of more massive collapsed haloes as compared to the collapse of substructure in absence of the plane wave. The formation of more massive haloes happens without any increase in the total mass in collapsed haloes. A comparison with the Burgers' equation approach in absence of any substructure suggests that the preferred value of effective viscosity depends primarily on the number of streams in a region.

N-Body simulations are a very important tool in the study of formation of large scale structures. Much of the progress in understanding the physics of high redshift universe and comparison with observations would not have been possible without N-Body simulations. Given the importance of this tool, it is essential to understand its limitations as ignoring the limitations can easily lead to interesting but unreliable results. In this paper we study the limitations arising out of the finite size of simulation volume. This finite size implies that modes larger than the size of the simulation volume are ignored and a truncated power spectrum is simulated. If the simulation volume is large enough then the mass in collapsed haloes expected from the full power spectrum and from the truncated power spectrum should match. We propose a quantitative measure based on this approach that allows us to compute the minimum box size for an N-Body simulation. We find that the required box size for simulations of LCDM model at high redshifts is much larger than is typically used. We can also use this approach to quantify the effect of perturbations at large scales for power law models and we find that if we fix the scale of non-linearity, the required box size becomes very large as the index becomes small. The appropriate box size computed using this approach is also an appropriate choice for the transition scale when tools like MAP (Tormen and Bertschinger, 1996) that add the contribution of the missing power are used.

In this review we discuss Cosmological N-Body codes with a special emphasis on Particle Mesh codes. We present the mathematical model for each component of N-Body codes. We compare alternative methods for computing each quantity by calculating errors for each of the components. We suggest an optimum set of components that can be combined reduce overall errors in N-Body codes.

The conceptual difficulties associated with a cosmological constant have led to the investigation of alternative models in which the equation of state parameter, $w=p/\rho$, of the dark energy evolves with time. We show that combining the supernova type Ia observations {\it with the constraints from WMAP observations} restricts large variation of $\rho(z)$ at low redshifts. The combination of these two observational constraints is stronger than either one. The results are completely consistent with the cosmological constant as the source of dark energy.

Cosmological N-Body simulations have become an essential tool for studying formation of large scale structure. These simulations are computationally challenging even though the available computing power gets better every year. A number of efficient algorithms have been developed to run large simulations with better dynamic range and resolution. We discuss key algorithms in this review, focusing on techniques used and their efficacy. N-Body simulations solve a model that is an approximation of the physical model to be simulated, we discuss limitations arising from this approximation and techniques employed for solving equations. Apart from simulating models of structure formation, N-Body simulations have also been used to study aspects of gravitational clustering. Simulating formation of galaxies requires us to take many physical process into account; we review numerical implementations of key processes.

We present an algorithm for parallelising the TreePM code. We use both functional and domain decompositions. Functional decomposition is used to separate the computation of long range and short range forces, as well as the task of coordinating communications between different components. Short range force calculation is time consuming and benefits from the use of domain decomposition. We have tested the code on a Linux cluster. We get a speedup of 31.4 for 128^3 particle simulation on 33 processors; speedup being better for larger simulations. The time taken for one time step per particle is 6.5 micro seconds for a 256^3 particle simulation on 65 processors, thus a simulation that runs for 4000 time steps takes 5 days on this cluster.

We present a detailed study of cosmological effects of homogeneous tachyon matter coexisting with non-relativistic matter and radiation, concentrating on the inverse square potential and the exponential potential for the tachyonic scalar field. A distinguishing feature of these models (compared to other cosmological models) is that the matter density parameter and the density parameter for tachyons remain comparable even in the matter dominated phase. For the exponential potential, the solutions have an accelerating phase, followed by a phase with a(t) ~ t^{2/3} as t goes to infinity. This eliminates the future event horizon present in LCDM models and is an attractive feature from the string theory perspective. A comparison with supernova Ia data shows that for both the potentials there exists a range of models in which the universe undergoes an accelerated expansion at low redshifts and are also consistent with requirements of structure formation. They do require fine tuning of parameters but not any more than in the case of LCDM or quintessence models.

We describe the TreePM method for carrying out large N-Body simulations to study formation and evolution of the large scale structure in the Universe. This method is a combination of Barnes and Hut tree code and Particle-Mesh code. It combines the automatic inclusion of periodic boundary conditions of PM simulations with the high resolution of tree codes. This is done by splitting the gravitational force into a short range and a long range component. We describe the splitting of force between these two parts. We outline the key differences between TreePM and some other N-Body methods.

We estimate the signal from large scale structure at high redshifts in redshifted 21cm line. We focus on z ~ 3 and the LCDM cosmology. We assume that neutral hydrogen is to be found only in galaxies, and normalise the total content to the density parameter of neutral hydrogen in damped Lyman alpha absorption systems (DLAS). We find that the rms fluctuations due to the large scale distribution of galaxies is very small and cannot be observed at angular scales probed by present day telescopes. We have used the sensitivity of the Giant meter-wave Radio Telescope (GMRT) for comparison. We find that observations stretching over 1000hours will be required for a 3 sigma detection of such objects.

We present a detailed analysis of the error budget for the TreePM method for doing cosmological N-Body simulations. It is shown that the choice of filter for splitting the inverse square force into short and long range components suggested in Bagla (2002) is close to optimum. We show that the error in the long range component of the force contributes very little to the total error in force. Errors introduced by the tree approximation for the short range force are different from those for the inverse square force, and these errors dominate the total error in force. We calculate the distribution function for error in force for clustered and unclustered particle distributions. This gives an idea of the error in realistic situations for different choices of parameters of the TreePM algorithm. We test the code by simulating a few power law models and checking for scale invariance.

We use large volume, high resolution, N-body simulations of 3 different $\Lambda$CDM models, with different clustering strengths, to generate dark matter halo merging histories. Over the reliable range of halo masses, roughly galaxy groups to rich clusters of galaxies, we quantify the number density of major mergers for two different time intervals and compare with analytic predictions based on the extended Press-Schechter (1974) theory.

We study the observed evolution of galaxy clustering as a function of redshift. We find that the clustering of galaxies decreases as we go from observations of the local Universe to $z \sim 2$. On the other hand, clustering of the Lyman break galaxies at $z \sim 3$ is very strong, comparable to the clustering of present day galaxies. However there are three major factors to take into account while comparing clustering measurements coming from various surveys: the so-called ``scale-dependence'' effect, due to measurements being made at different scales; the ``type-selection'' effect introduced by the fact that different galaxy surveys select different populations which do not have the same clustering amplitudes; and the Malmquist bias. We correct for the first two effects and discuss the implications of Malmquist bias on the interpretation of the data at different $z$. Then we compare the observed galaxy clustering with models for the evolution of clustering in some fixed cosmologies. Correcting for the scale-dependence effect significantly reduces the discrepancies amongst different measurements. We interpret the observed clustering signal at high redshift as coming from objects which are highly biased with respect to the underlying distribution of mass; this is not the case for $z\simlt 2$ where measurements are compatible with the assumption of a much lower biasing level which only shows a weak dependence on $z$. Finally we propose a model-independent test that can be used to place a lower limit on the density parameter $\Omega_0$.

We study the evolution of correlation function of dark matter halos in the CDM class of models. We show that the halo correlation function does not evolve in proportion with the correlation function of the underlying mass distribution. Earliest halos to collapse, which correspond to rare peaks in the density field, cluster very strongly. The amplitude of halo correlation function decreases from its initial, large, value. This decrease continues till the average peaks have collapsed, after which, the amplitude grows at a slow rate. This behaviour is shown to be generic and the epoch of minimum amplitude depends only on the rms fluctuations in mass at the relevant scale and, to a much smaller extent, on the slope of the power spectrum at that scale. We discuss the relevance of this result for interpretation of observations of galaxy and quasar clustering.

Recent observations show a large concentration of galaxies at high redshift. At first sight strong clustering of galaxies at high redshifts seems to be in contradiction with the models of structure formation. In this paper we show that such structures are a manifestation of the strong clustering of rare peaks in the density field. We compute the frequency of occurrence of such large concentrations of galaxies in some models of structure formation.

Large inhomogeneities in neutral hydrogen in the universe can be detected at redshifts $z \leq 10$ using the redshifted 21cm line emission from atomic hydrogen. This paper reviews the expected evolution of neutral hydrogen and presents estimates for future surveys of HI at $z \approx 3$. We also discuss the possibility of detecting neutral hydrogen at higher redshifts.

A gamma-ray burst (GRB) releases an amount of energy similar to that of a supernova explosion, which combined with its rapid variability suggests an origin related to neutron stars or black holes. Since these compact stellar remnants form from the most massive stars not long after their birth, gamma-ray bursts should trace the star formation rate in the Universe; we show that the GRB flux distribution is consistent with this. Because of the strong evolution of the star formation rate with redshift, it follows that the dimmest known bursts have z ~ 6, much above the value usually quoted and beyond the most distant quasars. This explains the absence of bright galaxies in well-studied gamma-ray burst error boxes. The increased distances imply a peak luminosity of 8.3e51 erg/s and a rate density of 0.025 per million years per galaxy. These values are 20 times higher and 150 times lower, respectively, than are implied by fits with non-evolving GRB rates. This means that GRBs are either caused by a much rarer phenomenon than mergers of binary neutron stars, or their gamma-ray emission is often invisible to us due to beaming. Precise burst locations from optical transients will discriminate between the various models for GRBs from stellar deaths, because the distance between progenitor birth place and burst varies greatly among them. The dimmest GRBs are then the most distant known objects, and may probe the Universe at an age when the first stars were forming.

It is known that radial collapse around density peaks can explain the key features of evolution of correlation function in gravitational clustering in three dimensions. The same model also makes specific predictions for two dimensions. In this paper we test these predictions in two dimensions with the help of N-Body simulations. We find that there is no stable clustering in the extremely non-linear regime, but a nonlinear scaling relation does exist and can be used to relate the linear and the non-linear correlation function. In the intermediate regime, the simulations agree with the model.

We show that the galaxy correlation function does not evolve in proportion with the correlation function of the underlying mass distribution. Earliest galaxies cluster very strongly and the amplitude of the galaxy correlation function decreases from this large value. This continues till the average peaks have collapsed, after which, the galaxy correlation function does not evolve very strongly.

Large inhomogeneities in neutral hydrogen in the universe can be detected at redshifts $z \leq 10$ using the redshifted 21 cm line emission. We use cosmological N-Body simulations for dark matter and a simple model for baryonic collapse to estimate the signal expected from structures like proto-clusters of galaxies at high redshifts. We study: (i) the standard CDM model, (ii) a modified CDM model with less power at small scales, and (iii) a $\Lambda+$CDM model in a universe with $\Omega_0 + \Omega_\Lambda = 1$. We show that it should be possible for the next generation radio telescopes to detect such structures at the redshift 3.34 with an integration of about 100 hours. We also discuss possible schemes for enhancing signal to noise ratio to detect proto-condensates at high redshifts.

We investigate the transfer of power between different scales and coupling of modes during non-linear evolution of gravitational clustering in an expanding universe. We start with a power spectrum of density fluctuations that is exponentially damped outside a narrow range of scales and use numerical simulations to study evolution of this power spectrum. Non-Linear effects generate power at other scales with most power flowing from larger to smaller scales. The ``cascade'' of power leads to equipartition of energy at smaller scales, implying a power spectrum with index $n\approx -1$. We find that such a spectrum is produced in the range $1 < \delta < 200$ for density contrast $\delta$. This result continues to hold even when small scale power is added to the initial power spectrum. Semi-analytic models for gravitational clustering suggest a tendency for the effective index to move towards a critical index $n_c\approx -1$ in this range. For n<n_c, power in this range grows faster than linear rate, while if n>n_c, it grows at a slower rate - thereby changing the index closer to n_c. At scales larger than the narrow range of scales with initial power, a k^4 tail is produced. We demonstrate that non-linear small scales do not effect the growth of perturbations at larger scales.

Thanks to new technology of observations and fresh inputs from particle physics, cosmology has advanced on both observational and theoretical fronts. It is therefore opportune that we take stock of the cosmological situation today and examine the observational and theoretical constraints as they are now. The bottom line in this review is that despite the availability of the cosmological constant as an extra parameter for flat Friedmann models, the allowed parameter space for such models is very small. The observations that we consider here include the ages of globular clusters, measurement of Hubble's constant, abundance of rich clusters of galaxies, fraction of mass contributed by baryons in rich clusters and abundance of high redshift objects.

In an expanding universe, velocity field and gravitational force field are proportional to each other in the linear regime. Neither of these quantities evolve in time and these can be scaled suitably so that the constant of proportionality is unity and velocity and force field are equal. The Zeldovich approximation extends this feature beyond the linear regime, until formation of pancakes. Nonlinear clustering which takes place {\it after} the breakdown of Zeldovich approximation, breaks this relation and the mismatch between these two vectors increases as the evolution proceeds. We suggest that the difference of these two vectors could form the basis for a powerful, new, statistical indicator of nonlinear clustering. We define an indicator called velocity contrast, study its behaviour using N-Body simulations and show that it can be used effectively to delineate the regions where nonlinear clustering has taken place. We discuss several features of this statistical indicator and provide simple analytic models to understand its behaviour. Particles with velocity contrast higher than a threshold have a correlation function which is biased with respect to the original sample. This bias factor is scale dependent and tends to unity at large scales.

Alignment of velocity and acceleration before shell crossing, and later misalignment are used to define velocity contrast, an indicator of dynamical state of matter undergoing gravitational collapse. We use this to study bias in clustering properties of dynamically nonlinear mass.

We study the evolution of the power spectrum of gravitational potential during the nonlinear clustering in an $\Omega=1$ matter dominated phase. N-body simulations suggest that the potential does not evolve in time even in the quasilinear phase for an $n=-1$ power spectrum. For $n=-2$, the potential evolves in the quasilinear phase but not in the extreme nonlinear phase. Becauase of these facts, the evolution of the gravitational potential in spectra like CDM is less than what would have been expected naively. We discuss a class of CDM like models in which such an interesting conspiracy between the amplitude and local index occurs.

: From the epoch of recombination $(z\approx 10^3)$ till today, the typical density contrasts have grown by a factor of about $10^6$ in a Friedmann universe with $\Omega=1$. However, during the same epoch the typical gravitational potential has grown only by a factor of order unity. We present theoretical arguments explaining the origin of this approximate constancy of gravitational potential. This fact can be exploited to provide a new, powerful, approximation scheme to study the formation of nonlinear structures in the universe. The essential idea of this method is to evolve the initial distribution of particles using a gravitational potential frozen in time (Frozen Potential Approximation). This approximation provides valuable insight into understanding various features of nonlinear evolution; for example, it provides a simple explanation as to why pancakes remain thin during the evolution even in the absence of any artificial, adhesion-like, damping terms. We compare the trajectories of particles in various approximations. We also discuss a few applications of the frozen potential approximation.