# Changes in the halo formation rates due to features in the primordial spectrum

## Abstract

Features in the primordial scalar power spectrum provide a possible roadway to describe the outliers at the low multipoles in the WMAP data. Apart from the CMB angular power spectrum, these features can also alter the matter power spectrum and, thereby, the formation of the large scale structure. Carrying out a complete numerical analysis, we investigate the effects of primordial features on the formation rates of the halos. We consider a few different inflationary models that lead to features in the scalar power spectrum and an improved fit to the CMB data, and analyze the corresponding imprints on the formation of halos. Performing a Markov Chain Monte Carlo analysis with the WMAP seven year data and the SDSS halo power spectrum from LRG DR7 for the models of our interest, we arrive at the parameter space of the models allowed by the data. We illustrate that, inflationary potentials, such as the quadratic potential with sinusoidal modulations and the axion monodromy model, which generate certain repeated, oscillatory features in the inflationary perturbation spectrum, do not induce a substantial difference in the number density of halos at their best fit values, when compared with, say, a nearly scale invariant spectrum as is generated by the standard quadratic potential. However, we find that the number density and the formation rates of halos change by about – for halo masses ranging over –, for potential parameters that lie within - around the best fit values arrived at from the aforesaid joint constraints. We briefly discuss the implications of our results.

###### Contents:

## 1 Introduction

Eversince its original proposal more than three decades ago, inflation has continued to be the most efficient scenario to overcome the horizon and the flatness problems associated with the conventional, hot big bang model. Additionally, inflation also provides an efficient mechanism for sowing the seeds of the density perturbations in the early universe [1, 2]. Over the last few years, observations of the Cosmic Microwave Background (CMB) by missions such as the Wilkinson’s Microwave Anisotropy Probe (WMAP) [3, 4] and the Atacama Cosmology Telescope (ACT) [5] have provided us with observational bounds on the parameters that describe the power spectra of the primordial perturbations. Currently, apart from the strong constraints on the primordial scalar amplitude and the spectral index, there also exists an upper bound on the tensor-to-scalar ratio. Clearly, the primordial power spectra generated by inflation must be consistent with these observational constraints.

The simplest mechanism to achieve inflation consists of a canonical scalar field (referred to as the inflaton) which rolls down slowly towards the minima of its potential. The perturbations originating from such a scalar field produces an almost scale independent power spectrum [1, 2]. Such a primordial spectrum, along with a suitable choice of parameters to describe the background CDM model, is found to lead to a CMB angular power spectrum which fits the available observational data quite well. In fact, despite the constantly improving bounds, one finds that there exist many inflationary models that remain consistent with the observations.

However, there exist a few points (notably near the multipoles of , and ) in the WMAP data [3, 4], which lie outside the cosmic variance associated with the best fit theoretical curve corresponding to the CMB angular power spectra generated by slow roll inflation. Though such outliers are few in number, demanding a better fit to these data points can, obviously, lead to a more favored model of inflation. Several attempts have been made to construct inflationary potentials that lead to wiggles in the primordial power spectrum (which, in turn, cause corresponding oscillations in the CMB angular power spectrum), thereby resulting in an improved fit to the data when compared to the featureless power spectra produced by slow roll inflation [6, 7, 8, 9]. The wiggles can be generated if the inflaton goes through a step or a bump in the potential. The presence of a step or a bump leads to a period of fast roll, which then leaves its imprints on the slow roll parameters and the power spectrum [10, 11]. Apart from constructing specific models that are more favored by the data, there has also been a few efforts towards a model independent reconstruction of the primordial scalar power spectrum from the observed CMB angular power spectrum [12]. Importantly, such reconstructions too seem to support the presence of features in the primordial power spectrum. Although the reconstruction efforts suggest that the features in the scalar power spectrum are typically localized over certain scales, it has been found that, oscillatory potentials which generate features on all scales also provide a much better fit to the CMB data than the nearly scale invariant power spectrum [13, 14, 15]. It may be worth mentioning here that inflationary models leading to features also gain importance due to the fact that they can generate large levels of non-Gaussianity as is possibly indicated by the WMAP data (in this context, see, for example, Refs. [16, 17] and the references therein).

Evidently, any significant change in the primordial power spectrum will correspondingly modify the matter power spectrum evaluated today. Hence, the features in the scalar perturbation spectrum that we mentioned above will affect the matter power spectrum which, in turn, will leave its signatures on the formation rates of the dark matter halos. The effects of steps or oscillations in the inflaton potential on the predicted number of halos was recently investigated based on the perturbation spectrum evaluated using the slow roll approximation [18]. However, the slow roll parameters can often turn large in such inflationary potentials and, as a result, the power spectrum evaluated in the slow roll approximation can differ considerably from the actual power spectrum. Moreover, with the ever increasing quality of the large scale structure observations, in particular, the halo power spectrum constructed using the Luminous Red Galaxies (LRG) from the seventh data release (DR7) of the Sloan Digital Sky Survey (SDSS) [19, 20], it would be a worthwhile exercise to actually compare the models with the available data to arrive at additional constraints (i.e. apart from those obtained from the CMB data) on the primordial features.

In this paper, we utilize a recently developed and accurate Fortran 90 code to exactly evaluate the inflationary perturbation spectrum, and thereby the corresponding matter power spectrum and the number of halos formed. We compute the number densities and the formation rates of halos in a couple of inflationary models that are known to lead to certain features in the power spectrum and also an improved fit to the CMB data. Our goal will be to estimate the maximum possible change in the number density and the formation rate of halos in these models, when compared with, say, inflation driven by the simplest quadratic potential. In order to arrive at the parameter space of interest, we shall compare the models with the CMB as well as the large scale structure data. Specifically, we shall make use of the WMAP seven year (WMAP-) data [4] and the LRG halo power spectrum data from SDSS DR7 [20]. We find that the power spectra with features corresponding to the best fit values of the inflationary parameters do not typically lead to substantial deviations in the formation rates of halos. However, we find that, in certain models that we consider, the potential parameters that lie within - of the best fit values obtained from the joint constraints of the WMAP and the SDSS data can lead to a – change in the number of halos formed for halo masses ranging over –.

This paper is organized as follows. In the following section, we shall quickly describe the inflationary potentials of interest, which lead to specific features in the primordial scalar power spectrum. We shall also outline the method that we adopt to compare the models with the data. In Section 3, we shall outline the formalism to arrive at the matter power spectra and the halo formation rates from the inflationary power spectra. In Section 4, we shall provide a few essential details concerning the numerical procedures that we follow. We shall discuss the results in Section 5, and we shall close with a few concluding remarks regarding the wider implications of our results in Section 6.

Note that we shall assume the background cosmological model to be the standard, spatially flat, CDM model.

## 2 Primordial spectra with features

In this section, we shall quickly sketch a few essential points concerning the scalar power spectrum generated during inflation. We shall also discuss the inflationary models of our interest and the scalar power spectra produced by them. Further, we shall outline the methods that we adopt to compare the models with the data.

### 2.1 Essentials

We shall focus on the simplest case of inflation driven by the canonical scalar field, say, . In a spatially flat Friedmann universe, given a potential , the scalar field satisfies the differential equation

(1) |

where the overdot represents differentiation with respect to the cosmic time, is the Hubble parameter, with being the scale factor, while . The evolution of the scalar field is usually characterized by a hierarchy of the so-called slow roll parameters. The first slow roll parameter is set to be , and the higher order slow roll parameters are defined in terms of the first as follows: , where and denotes the number of e-folds.

The scalar perturbations induced by the quantum fluctuations of the inflaton can be described by the curvature perturbation, say, , whose Fourier modes satisfy the following differential equation:

(2) |

where the overprime denotes differentiation with respect to the conformal time coordinate, and . Upon solving this differential equation with suitable initial conditions for modes of cosmological interest, one arrives at the scalar power spectrum (see Section 4 for the necessary details), defined as

(3) |

with being evaluated at sufficiently late times. Smooth potentials, in general, permit slowly rolling fields (i.e. wherein the parameters are much smaller than unity), and such backgrounds generate nearly scale invariant power spectra. For instance, the archetypical quadratic potential, viz. , leads to the required COBE amplitude (for a suitable value of ) and an almost scale invariant spectrum with a scalar spectral index of , in strong conformity with the recent CMB data [3, 4].

### 2.2 Models of interest

In contrast to smooth potentials, potentials that contain either non-trivial forms or sharp changes in their slopes lead to deviations from slow roll. Departures from slow roll affect the amplitude of the modes that leave the Hubble radius during this period resulting in specific features in the inflationary scalar power spectrum, with the shape being determined by the type of deviation from slow roll (in this context, see, for instance, Ref. [21]). Interestingly, certain features in the scalar power spectrum are known to allow an improved fit to the CMB data than the conventional nearly scale invariant spectrum, as is generated by slow roll inflation [6, 7, 8, 9, 13, 14, 15].

In this work, we shall consider two types of inflationary models that have been shown to fit the data better than the power law primordial spectrum. The first type of model that we shall consider contains a step, which is typically introduced in the conventional quadratic potential at a certain location, say, , as follows [6, 9]:

(4) |

Evidently, and denote the strength and the width of the step, respectively. The field experiences a short period of fast roll as it crosses , leading to a brief burst of oscillations in the scalar power spectrum. A Markov Chain Monte Carlo (MCMC) sampling of the parameter space of the above potential points to the fact that the presence of the step aids in fitting the outliers in the WMAP data near the multipoles of and . In fact, at the expense of the three additional parameters, viz. , and , the model results in an improvement in the effective least squares parameter by about , when compared with the quadratic potential without the step. In Figure 1, we have plotted the behavior of the first two slow roll parameters in the model.

Over the last few years, two types of oscillatory inflationary potentials have drawn a considerable amount attention [13, 14, 15, 17]. The first of these consists of the above-mentioned quadratic potential that is modulated by sinusoidal oscillations as follows:

(5) |

The second is the so-called axion monodromy model that has its motivations in string theory, and is described by the potential

(6) |

In both these cases, the slow roll parameters oscillate and, in fact, continue to exhibit such a behavior until the termination of inflation (cf. Figure 1). All modes of interest are affected by these oscillations, which leave a repeated pattern that extends over all scales in the power spectrum. We shall illustrate the power spectra resulting in these models later in Section 5 (see Figure 3), wherein we shall present the results of our analysis. We would like to mention here that recent analysis indicate that such power spectra lead to an improved fit to the WMAP as well as the small scale ACT data. In the case of the quadratic potential modulated by sinusoidal oscillations, one finds that the least squares parameter reduces by about unity, whereas the monodromy model is found to lead to a much better fit to the data with an improvement in of up to , when compared with a more conventional, nearly scale invariant spectrum [14, 15].

### 2.3 Comparison with the WMAP and the SDSS data

We shall compare the models with the CMB as well as the large scale structure data. We have worked with the WMAP- data [4] and the halo power spectrum data arrived at from the LRG in SDSS DR7 [20]. We have made use of the large scale structure data to ensure that the parameter values we eventually work with to obtain the formation rate of the halos are consistent with the observed matter power spectrum.

We have made use of the publicly available codes, viz. the cosmological Boltzmann code CAMB [22, 23] and the Monte Carlo code COSMOMC [24, 25], to compute the CMB angular and the matter power spectra, and compare them with the data, respectively. As we had mentioned, we shall assume the background model to be the spatially flat CDM and we shall work with the priors on the background parameters as listed in Table 1.

Background | Lower | Upper |
---|---|---|

parameter | limit | limit |

We should note here that the step model (4), the quadratic potential with the sinusoidal modulation (5) and the axion monodromy model (6), all require four parameters to describe the potential completely. As far as the priors on these parameters are concerned, for the case of the quadratic potential with the step, we have worked with the same priors as we had worked with before [9]. In the case of the quadratic potential with superimposed sinusoidal modulations and the axion monodromy model, for the primary parameters and , we work with the same priors that we had considered in a recent analysis wherein we had compared the models with the CMB data [15]. However, while comparing with the WMAP as well as the SDSS data, we have widened the priors of the parameters and for both the oscillatory potentials as listed in Table 2. We should add that we have allowed the phase parameter to vary from to as before.

Model | Potential | Lower | Upper |
---|---|---|---|

parameter | limit | limit | |

Quadratic potential | |||

with sinusoidal modulation | |||

Axion monodromy | |||

model |

The motivations for these choice of priors are two fold. Firstly, while determining the priors on the parameters and , we have ensured that the resulting inflationary power spectra do not differ considerably from the nearly scale invariant spectrum for the faster convergence of the Markov Chains. Secondly, we have chosen the priors on and such that the scalar field does not get trapped by the oscillations in the potential.

We should mention that, unlike in our earlier efforts [9, 15], we have not taken the effects of the tensor perturbations into account, as the corresponding effects are negligible. Moreover, in the case of potentials with oscillations, which lead to fine features in the inflationary scalar power spectrum, we actually need to modify CAMB in order to ensure that the CMB angular power spectrum is evaluated at every multipole and compare them with the data [14, 15, 26]. But, we have not implemented this point here since we are only interested in the marginalized probabilities of the potential parameters. These probabilities shall indicate the extent to which deviations from a nearly scale invariant spectrum is allowed by the data, and the corresponding effects on the formation of dark matter halos. We should point out that we have not taken into account the non-linear effects on the matter power spectrum [27], but have included the SZ effect and the effects due to gravitational lensing in our analysis. Finally, we shall set the Gelman and Rubin parameter to be for convergence in all the cases.

## 3 From the primordial spectrum to the formation rate of halos

In this section, we shall quickly outline the standard formalism to arrive at the formation rate of halos from the primordial power spectrum.

### 3.1 The matter power spectrum

Given a primordial power spectrum, say, [cf. Eq. (3)], the matter power spectrum at the redshift is usually written as (see, for example, Refs. [28, 29])

(7) |

a quantity that can be conveniently expressed in units of for the modes of cosmological interest. In the above expression, is the CDM transfer function, the quantity denotes the linear growth factor of the total matter perturbation, is the non-relativistic density parameter, while and denote the scale factor and the Hubble parameter today.

If we define , then, one finds that, in the spatially flat CDM model, the quantity satisfies the differential equation [3, 30]

(8) |

where , with denoting the dimensionless density parameter associated with the cosmological constant today. In this work, we shall solve the above differential equation with suitable initial conditions to obtain the growth factor . Utilizing CAMB [22, 23] to determine the corresponding transfer function and, upon using the primordial spectrum obtained numerically in the inflationary model of interest, we shall eventually arrive at the matter power spectrum (see Section 4 on numerical methods for further details).

### 3.2 Mass functions and the halo formation rates

To arrive at the formation rate of dark matter halos, we shall first require the number density of collapsed halos with mass in the range of and in a comoving volume element. This number density, say, , is defined in terms of the root mean square fluctuation in mass through the so-called mass function as follows [31]:

(9) |

where is the mean density of non-relativistic matter in the universe. Following the convention (see, for example, Ref. [28]), we shall define the root mean square fluctuation in mass at the scale to be

(10) |

where denotes the dimensionless matter power spectrum, while is the Fourier transform of the window function that is introduced to smooth out the density perturbation. We shall work with the commonly used spherical top hat window function, whose Fourier transform is given by

(11) |

corresponding to the volume . Note that the halo mass within the window of radius is given by .

We shall make use of the Sheth-Tormen mass function to evaluate the number density of halos [32]. In contrast to the more conventional Press-Schechter mass function [33], it has been found that the Sheth-Tormen mass function fits the data from the -body simulations better. Actually, the Sheth-Tormen mass function is a generalization of the original Press-Schechter formalism for spherical collapse to the case of ellipsoidal collapse. The Sheth-Tormen mass function is defined in terms of two additional parameters and (when compared to the Press-Schechter case) as follows:

(12) |

where , with being the threshold linear overdensity for collapse. The Press-Schechter mass function corresponds to , and . However, upon comparing with the -body simulation data, the best fit values for and are found out to be and , respectively. The value of can then be arrived at from the normalization condition on , viz. that the integral of over all is unity, which leads to .

The number density of halos associated with the above Sheth-Tormen mass function is then given by

(13) |

The corresponding formation rates of the halos can be easily obtained to be [34]

(14) |

Note that the quantity proves to be negative, since the growth factor decreases as the redshift increases. As a result, it is known that the above formation rate of halos can become negative for some mass scales (i.e. when ), which in practice can not occur. Therefore, to avoid this issue and simultaneously illustrate the effects of features, we shall only plot the ratio of the formation rates in the inflationary models leading to features and the conventional, smooth, quadratic potential.

## 4 Details of the numerical methods

The slow roll approximation allows the background as well as the perturbations to be evaluated analytically during inflation. As we pointed out before, in the earlier work [18], it was the slow roll approximation that was made use of in order to arrive at the inflationary perturbation spectrum . However, we find that, in the inflationary models of interest, the slow roll parameters turn sufficiently large (cf. Figure 1) implying a breakdown of the slow roll approximation. As a result, we resort to numerical methods to compute the primordial perturbation spectrum. In fact, we utilize a Fortran 90 code that has been recently developed by us, which makes use of a Bulirsch-Stoer algorithm along with an adaptive step size control routine [35] to accurately and efficiently solve the equations governing the background and the perturbations, in order to arrive at the inflationary scalar power spectrum (see Refs. [9, 15] for further details in this regard).

Having computed the primordial power spectrum, we arrive at the matter power spectrum using the transfer function and the growth factor. As we had remarked earlier, we obtain the transfer function from CAMB, and we evaluate the growth factor by solving the differential equation (8). It should be mentioned here that the initial conditions are chosen such that is a constant and equal to unity in the early matter dominated epoch, i.e. at a sufficiently high redshift of, say, [3, 36].

After having obtained the matter power spectrum, we calculate the variance using equation (10). The integral can be evaluated numerically with the simplest of algorithms, provided the power spectrum proves to be smooth and devoid of any features. In contrast, when these exist features such as repeated oscillations, certain care is required, and we have made use of an adaptive integration routine to compute the integral involved [37]. We have carried out integral from a suitably small mode (such as ) up to a mode where the window function cuts off the integrand. Finally, we obtain the quantity by numerical differentiation. We should stress here that, keeping in mind the presence of oscillations in the power spectra, we have computed and with care and high accuracy. We should also add that we have cross checked our result by fitting the numerical values of to the Chebyshev polynomials and calculating the corresponding derivative from the polynomial (in this context, see Ref. [36]).

## 5 Results

In this section, we shall present the results of our comparison of the models of our interest with the CMB and the large scale structure data. We shall also discuss the effects of primordial spectra with features on the formation of halos.

### 5.1 Joint constraints from the WMAP and the SDSS data

In Table 3 below, we have tabulated the best fit values of the background and the potential parameters obtained from the MCMC analysis using the WMAP- and the SDSS LRG DR7 data. We have also listed the effective least squared parameter in each of the cases. For the case of the quadratic potential with and without the step, we have arrived at results similar to what we have obtained in an earlier work [9]. Also, as one would expect, we find that the background parameters are better constrained with the inclusion of the additional SDSS data [38, 39]. Moreover, it is obvious from Table 3 that the axion monodromy model does not lead to the same extent of improvement in the fit as has been obtained before (in this context, see Refs. [14, 15]). This arises due to the fact that, unlike in the earlier analysis, we have not evaluated the CMB angular power at each multipole, but have worked with the inbuilt effective sampling and interpolation routine in CAMB. However, we should stress that this does not affect our conclusions since our focus here lies on the maximum change in the formation of halos. Therefore, we are more interested in the allowed regions of the parameter space of rather in arriving at the precise best fit point. Also, importantly, as we shall discuss in the following subsection, for violent oscillations in the primordial power spectrum (when one requires computing the CMB angular power spectrum at each multipole explicitly and accurately), the percentage change in the number density of halos proves to be negligible in the observable mass bins.

Model | Quadratic | Quadratic step | Quadratic sine | Axion monodromy |
---|---|---|---|---|

0.0222 | 0.0221 | 0.0216 | 0.0225 | |

0.1162 | 0.1159 | 0.1168 | 0.1154 | |

1.038 | 1.039 | 1.036 | 1.039 | |

0.0824 | 0.0875 | 0.0836 | 0.0856 | |

-0.6545 | -0.6406 | -0.6448 | 0.9649 | |

- | ||||

- | 14.664 | - | - | |

- | - | - | ||

- | - | -2.576 | -7.61 | |

- | - | 2.208 | -1.178 | |

7515.57 | 7507.3 | 7515.12 | 7509.56 |

In Figure 2, we have illustrated the one dimensional likelihood on the parameter in the case of the quadratic potential with a step, and it is clear that the location of the step is highly constrained by data. The step affects the number density of halos over only highly localized mass scales. We have also plotted the marginalized two dimensional constraints on the parameters and for the cases of the two oscillatory potentials. It is noteworthy that the constraints are strikingly similar. In fact, the roughly triangular shape of the contours can also be understood. As the parameter decreases, the resulting oscillations in the potential and, therefore, in the inflationary perturbation spectrum turn too frequent, and the data constrains the amplitude to a smaller region.

In Figure 3, we have plotted the best fit scalar power spectrum for quadratic potential with and without the step and the two oscillatory potentials. Further below, in Figure 4, we have plotted the matter power spectrum evaluated today corresponding to the different inflationary power spectra in the previous figure. In the inset of the figure, we have highlighted the baryon acoustic oscillations and the halo power spectrum data from SDSS LRG DR7. We should add here that the theoretical best fit curves are unable to fit the data well after due to the fact that we have not taken the non-linear effects into account in arriving at the matter power spectrum.

### 5.2 Effects of features on the number density and the formation rate of halos

In this subsection, we shall discuss the effects of the features on the number density and the formation rates of halos in the different inflationary models of our interest. In order to highlight the effects purely due to the primordial features, we have frozen the values of the background cosmological parameters, viz. , , and the dimensionless baryon density parameter at the values arrived at upon comparing the smooth quadratic potential with the WMAP and SDSS data, as listed in Table 3. But, we have made use of the best fit values for the potential parameters to compute the inflationary scalar power spectrum and from thereon the matter power spectrum and the number density of halos. In Figure 5, we have plotted the percentage of change in the formation rate of halos in the Press-Schechter formalism and the number density of halos in the Sheth-Tormen formalism for different models with respect to the quadratic potential. In the case of the model with the step, the change in the number density due to the step (corresponding to parameter values within 2- around the best fit values) occurs at very high mass halos () and hence lies outside our region of interest. Due to this reason, we have only presented the results in the case of the models with oscillatory terms in the potential.

In order to arrive at the maximum possible change in the number density of halos when compared to the conventional nearly scale invariant primordial spectrum, for the models with oscillations in the potential, we have chosen values for the parameters and that they lie within (actually, at) - from the best fit values. We have chosen the parameters in such a way that they create the largest deviation from the nearly scale invariant power spectra that are allowed by the CMB and observations of the large scale structure. In Figure 5, apart from the results for the best fit values, we have plotted the number density and the formation rates of halos for the cases wherein [] is set to () and () (as indicated by the blue points in Figure 2) for the quadratic potential with sinusoidal modulations and the axion monodromy model, respectively. In arriving at these plots, we have fixed the parameters and at their best fit values as shown in the Table 3, since these parameters do not play a role in altering the features in the spectrum. We have also chosen the value of to be the best fit value for both the models.

It is evident from the figure that, for the best fit values of the parameters, the change in the number density is completely negligible (). However, we find that, for the case of the quadratic potential with sinusoidal modulation, the numbers can change by as much as for values of the potential parameters and that lie within -. It should also be highlighted that the monodromy model does not seem to lead to the same extent of change in the number density and the rate of formation of halos, despite the fact that it produces fine oscillations in the primordial as well as the linear matter power spectra (cf. Figures 3 and 4). Actually, while the unbinned number density does indicate a – change, we find that, the change proves to be smaller when we bin the numbers in logarithmic mass bins, i.e. , of . Evidently, binning seems to average out the rapid oscillations, resulting in a smaller extent of change in the numbers.

## 6 Discussion

In this work, we have investigated the effects of primordial features on the matter power spectrum as well as the number of halos formed and their rate of formation. Similar work in this context [18] had suggested that a small change in the parameters describing the inflaton potential would lead to a drastic change in the number of halos formed. The earlier work was based on the inflationary perturbation spectrum that was arrived at based on the slow roll approximation. In contrast, we have carried a complete and accurate numerical analysis. Further, we have made use of the Sheth-Tormen mass function (instead of the older Press-Schechter one) which is known to fit the data from the -body simulations better. We have included the baryon acoustic oscillations in our analysis to have a more realistic comparison. Moreover, to arrive at the parameter space of interest, all the potentials considered in our work have been constrained by an MCMC analysis (using COSMOMC) against the WMAP-7 and SDSS LRG DR7 datasets. We find that, the best fit values for the potential parameters (with the background parameters kept fixed) lead to hardly any change in the number of halos formed when compared to the conventional quadratic potential that generates a nearly scale invariant primordial spectrum. However, partly consistent with the earlier result, we find that values for the potential parameters that lie within - of the best fit values indeed lead to a reasonable change in the number of halos formed and in their formation rates. For instance, we find that, with superimposed sinusoidal modulations, the quadratic potential leads to as much as a – change in the halo number density and the rate of formation. Needless to mention, the step of comparing the models against the data is crucial as this imposes real bounds on the extent of changes in the numbers involved. It is worthwhile to note that the inclusion of SDSS data reduces the maximum change in number density to about , when compared to the case wherein one works with the parameters constrained by the WMAP data alone.

We would like to close this paper with the following remarks. As we had pointed out before, while comparing with the SDSS data, we have not taken into account the non-linear effects on the matter power spectrum. It is for this reason that the theoretical curve had not fit the observational data well on small scales (cf. Figure 4). Clearly, a more complete analysis should involve modeling of the non-linear effects and their inclusion in evaluating the matter power spectrum [27]. For instance, it will be interesting to compare the results on the number of halos formed in numerical simulations, evolved from primordial spectra with features, with the small scale data.

## Acknowledgments

The author wishes to thank Shiv Sethi and L. Sriramkumar for discussions as well as comments on the manuscript. Computational work for this study has been carried out using the cluster computing facilities at Harish-Chandra Research Institute, Allahabad, India (http://cluster.hri.res.in/).

### Footnotes

- Current address: Asia Pacific Center for Theoretical Physics, Pohang, Gyeongbuk 790-784, Korea. E-mail: dhiraj@apctp.org

### References

- E. W. Kolb and M. S. Turner, The Early Universe (Addison-Wesley, Redwood City, California, 1990); S. Dodelson, Modern Cosmology (Academic Press, San Diego, U.S.A., 2003); V. F. Mukhanov, Physical Foundations of Cosmology (Cambridge University Press, Cambridge, England, 2005); S. Weinberg, Cosmology (Oxford University Press, Oxford, England, 2008); R. Durrer, The Cosmic Microwave Background (Cambridge University Press, Cambridge, England, 2008); D. H. Lyth and A. R. Liddle, The Primordial Density Perturbation (Cambridge University Press, Cambridge, England, 2009); P. Peter, J-P. Uzan and J. Brujic, Primordial Cosmology (Oxford University Press, Oxford, England, 2009).
- H. Kodama and M. Sasaki, Prog. Theor. Phys. Suppl. 78, 1 (1984); V. F. Mukhanov, H. A. Feldman and R. H. Brandenberger, Phys. Rep. 215, 203 (1992); J. E. Lidsey, A. Liddle, E. W. Kolb, E. J. Copeland, T. Barreiro and M. Abney, Rev. Mod. Phys. 69, 373 (1997); D. H. Lyth and A. Riotto, Phys. Rep. 314, 1 (1999); A. Riotto, arXiv:hep-ph/0210162; J. Martin, arXiv:hep-th/0406011; B. Bassett, S. Tsujikawa and D. Wands, Rev. Mod. Phys. 78, 537 (2006); W. H. Kinney, arXiv:0902.1529 [astro-ph.CO]; L. Sriramkumar, Curr. Sci. 97, 868 (2009); D. Baumann, arXiv:0907.5424v1 [hep-th].
- J. Dunkley et al., Astrophys. J. Suppl. 180, 306 (2009); E. Komatsu et al., Astrophys. J. Suppl. 180, 330 (2009).
- D. Larson et al., Astrophys. J. Suppl. 192, 16 (2011); E. Komatsu et al., Astrophys. J. Suppl. 192, 18 (2011).
- J. Dunkley et al., Astrophys. J. 739, 52 (2011).
- J. A. Adams, B. Cresswell and R. Easther, Phys. Rev. D 64, 123514 (2001); L. Covi, J. Hamann, A. Melchiorri, A. Slosar and I. Sorbera, Phys. Rev. D 74, 083509 (2006); J. Hamann, L. Covi, A. Melchiorri and A. Slosar, Phys. Rev. D 76, 023503 (2007); M. J. Mortonson, C. Dvorkin, H. V. Peiris and W. Hu, Phys. Rev. D 79, 103519 (2009).
- M. Joy, V. Sahni and A. A. Starobinsky, Phys. Rev. D 77, 023514 (2008); M. Joy, A. Shafieloo, V. Sahni and A. A. Starobinsky, JCAP 0906, 028 (2009).
- R. K. Jain, P. Chingangbam, J.-O. Gong, L. Sriramkumar and T. Souradeep, JCAP 0901, 009 (2009); R. K. Jain, P. Chingangbam, L. Sriramkumar and T. Souradeep, Phys. Rev. D 82, 023509 (2010).
- D. K. Hazra, M. Aich, R. K. Jain, L. Sriramkumar and T. Souradeep, JCAP 1010, 008 (2010).
- A. A. Starobinsky, Sov. Phys. JETP Lett. 55, 489 (1992).
- C. Dvorkin and W. Hu, Phys. Rev. D 81, 023518 (2010); W. Hu, Phys. Rev. D 84, 027303 (2011).
- S. Hannestad, Phys. Rev. D 63, 043009 (2001); S. L. Bridle, A. M. Lewis, J. Weller and G. Efstathiou, Mon. Not. Roy. Astron. Soc. 342, L72 (2003); P. Mukherjee and Y. Wang, Astrophys. J. 599, 1 (2003); S. Hannestad, JCAP 0404, 002 (2004); A. Shafieloo and T. Souradeep, Phys. Rev. D 70, 043523 (2004); D. Tocchini-Valentini, Y. Hoffman and J. Silk, Mon. Not. Roy. Astron. Soc. 367, 1095 (2006); A. Shafieloo, T. Souradeep, P. Manimaran, P. K. Panigrahi and R. Rangarajan, Phys. Rev. D 75, 123502 (2007); A. Shafieloo and T. Souradeep, Phys. Rev. D 78, 023511 (2008); R. Nagata and J. Yokoyama, Phys. Rev. D 79, 043010 (2009); G. Nicholson and C. R. Contaldi, JCAP 0907, 011 (2009).
- C. Pahud, M. Kamionkowski and A. R. Liddle, Phys. Rev. D 79, 083503 (2009).
- R. Flauger, L. McAllister, E. Pajer, A. Westphal and G. Xu, JCAP 1006, 009 (2010).
- M. Aich, D. K. Hazra, L. Sriramkumar and T. Souradeep, arXiv:1106.2798v2 [astro-ph.CO].
- J. Maldacena, JHEP 0305, 013 (2003); D. Seery and J. E. Lidsey, JCAP 0506, 003 (2005); X. Chen, Adv. Astron. 2010, 638979 (2010).
- X. Chen, R. Easther, E. A. Lim, JCAP 0804, 010 (2008); R. Flauger and E. Pajer, JCAP 1101, 017 (2011); F. Arroja, A. E. Romano and M. Sasaki, Phys. Rev. D 84, 123503 (2011); J. Martin and L. Sriramkumar, JCAP 1201, 008 (2012); D. K. Hazra, L. Sriramkumar and J. Martin, arXiv:1201.0926v1 [astro-ph.CO]; F. Arroja and M. Sasaki, JCAP 1208, 012 (2012).
- L. F. S. Rodrigues and R. Opher, Phys. Rev. D 82, 023501 (2010).
- See, for instance, http://www.sdss.org/
- K. Abazajian et. al. Astrophys. J. Suppl. 182, 543-558 (2009)
- R. K. Jain, P. Chingangbam and L. Sriramkumar, JCAP 10, 003 (2007).
- See, http://camb.info/.
- A. Lewis, A. Challinor and A. Lasenby, Astrophys. J. 538, 473 (2000).
- See, http://cosmologist.info/cosmomc/.
- A. Lewis and S. Bridle, Phys. Rev. D 66, 103511 (2002).
- Z. Huang, JCAP 1206, 012 (2012).
- R. E. Smith et al., Mon. Not. Roy. Astron. Soc. 341, 1311 (2003); M. Tegmark et al., Phys. Rev. D 74, 123507 (2006); W. J. Percival et al., Astrophys. J. 657, 645 (2007).
- H. Mo, F. v. d. Bosch and S. White, Galaxy Formation and Evolution (Cambridge University Press, Cambridge, England, 2010).
- M. Takada, E. Komatsu and T. Futamase, Phys. Rev. D 73, 083520 (2006).
- L. Wang and P. J. Steinhardt, Astrophys. J. 508, 483 (1998); E. V. Linder and A. Jenkins, Mon. Not. Roy. Astron. Soc. 346, 573 (2003).
- A. Jenkins, C. S. Frenk, S. D. M. White, J. M. Colberg, S. Cole, A. E. Evrard, H. M. P. Couchman and N. Yoshida, Mon. Not. Roy. Astron. Soc. 321, 372 (2001).
- R. K. Sheth, H. J. Mo and G. Tormen, Mon. Not. Roy. Astron. Soc. 323, 1 (2001); R. K. Sheth and G. Tormen, Mon. Not. Roy. Astron. Soc. 329, 61 (2002).
- W. H. Press and P. Schechter, Astrophys. J. 187, 425 (1974).
- S. Sasaki, Publ. Astron. Soc. Jap. 46, 427 (1994); E. Ripamonti, Mon. Not. Roy. Astron. Soc. 376, 709 (2007); S. Mitra, G. Kulkarni, J. S. Bagla and J. K. Yadav, Bull. Astron. Soc. Ind. 39, 1,(2011).
- W. H. Press, S. A. Teukolsky, W. T. Vetterling and B. P. Flannery, Numerical Recipes in Fortran 90, Second edition (Cambridge University Press, Cambridge, England, 1996).
- See, http://www.mpa-garching.mpg.de/komatsu/CRL/index.html.
- See, for instance, http://www.nag.com/ and http://www.netlib.org/.
- F. Finelli, J. Hamann, S. M. Leach and J. Lesgourgues, JCAP 04 011 (2010).
- M. Benetti, S. Pandolfi, M. Lattanzi, M. Martinelli and A. Melchiorri, arXiv:1210.3562 [astro-ph.CO].