Showing 1–34 of 34 results for author: Khandai, N

Role of ionizing background on the statistics of metal absorbers in hydrodynamical simulations

Authors: Sukanya Mallik, Raghunathan Srianand, Soumak Maitra, Prakash Gaikwad, Nishikanta Khandai

Abstract: We study the statistical properties of O VI, C IV, and Ne VIII absorbers at low-$z$ (i.e., $z<0.5$) using Sherwood simulations with "WIND" only and "WIND+AGN" feedback and Massive black simulation that incorporates both "WIND" i.e. outflows driven by stellar feedback and AGN feedbacks. For each simulation, by considering a wide range of metagalactic ionizing UV background (UVB), we show the statis… ▽ More We study the statistical properties of O VI, C IV, and Ne VIII absorbers at low-$z$ (i.e., $z<0.5$) using Sherwood simulations with "WIND" only and "WIND+AGN" feedback and Massive black simulation that incorporates both "WIND" i.e. outflows driven by stellar feedback and AGN feedbacks. For each simulation, by considering a wide range of metagalactic ionizing UV background (UVB), we show the statistical properties such as distribution functions of column density ($N$), $b$-paramerer and velocity spread ($ΔV_{90}$), the relationship between $N$ and $b$-parameter and the fraction of Lya absorbers showing detectable metal lines as a function of $N$(H I) are influenced by the UVB used. This is because UVB changes the range in density, temperature, and metallicity of gas contributing to a given absorption line. For simulations considered here, we show the difference in some of the predicted distributions between different simulations is similar to the one obtained by varying the UVB for a given simulation. Most of the observed properties of O VI absorbers are roughly matched by Sherwood simulation with "WIND+AGN" feedback when using the UVB with a lower O VI ionization rate. However, this simulation fails to produce observed distributions of C IV and fraction of H I absorbers with detectable metals. Therefore, in order to constrain different feedback processes and/or UVBs, using observed properties of H I and metal ions, it is important to perform simultaneous analysis of various observable parameters. △ Less

Submitted 22 May, 2023; v1 submitted 25 January, 2023; originally announced January 2023.

Comments: 21 pages, 11 figures, 3 tables, accepted in MNRAS

Halo mass function in scale invariant models

Authors: Swati Gavas, J S Bagla, Nishikanta Khandai, Girish Kulkarni

Abstract: Sheth-Tormen mass function has been widely used to quantify the abundance of dark matter halos. It is a significant improvement over the Press-Schechter mass function as it uses ellipsoidal collapse in place of spherical collapse. Both of these mass functions can be written in a form that is universal, i.e., independent of cosmology and power spectrum when scaled in suitable variables. However, co… ▽ More Sheth-Tormen mass function has been widely used to quantify the abundance of dark matter halos. It is a significant improvement over the Press-Schechter mass function as it uses ellipsoidal collapse in place of spherical collapse. Both of these mass functions can be written in a form that is universal, i.e., independent of cosmology and power spectrum when scaled in suitable variables. However, cosmological simulations have shown that this universality is approximate. In this paper, we investigate the power spectrum dependence of halo mass function through a suite of dark-matter-only N-body simulations of seven power-law models in an Einstein-de Sitter cosmology. This choice of cosmology and a power-law power spectrum ensures the self-similar evolution of dark matter distribution, allowing us to isolate the power spectrum dependence of mass function. We find that the mass function shows a clear non-universality. We present fits for the parameters of the Sheth-Tormen mass function for a range of power-law power-spectrum indices. We find a mild evolution in the overall shape of the mass function with the epoch. Finally, we extend our result to LCDM cosmology. We show that the Sheth-Tormen mass function with parameter values derived from a matched power-law EdS cosmology provides a better fit to the LCDM mass function than the standard Sheth-Tormen mass function. Our results indicate that an improved analytical theory is required to provide better fits to the mass function. △ Less

Submitted 3 April, 2023; v1 submitted 31 October, 2022; originally announced October 2022.

Comments: 12 pages, 10 figures. This is a much expanded and upgraded version of 0908.2702. Accepted for Publication in the MNRAS

The Dark Matter Halos of HI Selected Galaxies

Authors: Saili Dutta, Nishikanta Khandai, Sandeep Rana

Abstract: We present the neutral hydrogen mass ($M_{\text{HI}}$) function (HIMF) and velocity width ($w_{50}$) function (HIWF) based on a sample of 7857 galaxies from the 40% data release of the ALFALFA survey ($α.40$). The low mass (velocity width) end of the HIMF (HIWF) is dominated by the blue population of galaxies whereas the red population dominates the HIMF (HIWF) at the high mass (velocity width) en… ▽ More We present the neutral hydrogen mass ($M_{\text{HI}}$) function (HIMF) and velocity width ($w_{50}$) function (HIWF) based on a sample of 7857 galaxies from the 40% data release of the ALFALFA survey ($α.40$). The low mass (velocity width) end of the HIMF (HIWF) is dominated by the blue population of galaxies whereas the red population dominates the HIMF (HIWF) at the high mass (velocity width) end. We use a deconvolution method to estimate the HI rotational velocity ($V_{\text{rot}}$) functions (HIVF) from the HIWF for the total, red, and blue samples. The HIWF and HIVF for the red and blue samples are well separated at the knee of the function compared to their HIMFs. We then use recent stacking results from the ALFALFA survey to constrain the halo mass ($M_{\text{h}}$) function of HI-selected galaxies. This allows us to obtain various scaling relations between $M_{\text{HI}} - w_{50} - V_{\text{rot}} - M_{\text{h}}$, which we present. The $M_{\text{HI}} - M_{\text{h}}$ relation has a steep slope ~2.10 at small masses and flattens to ~0.34 at masses larger than a transition halo mass, $\log_{10}(M_{\text{ht}} h^2_{70}/M_{\odot})=10.62$. Our scaling relation is robust and consistent with a volume-limited sample of $α.40$. The $M_{\text{HI}} - M_{\text{h}}$ relation is qualitatively similar to the $M_{\text{star}} - M_{\text{h}}$ relation but the transition halo mass is smaller by ~1.4 dex compared to that of the $M_{\text{star}} - M_{\text{h}}$ relation. Our results suggest that baryonic processes like heating and feedback in larger mass halos suppress HI gas on a shorter time scale compared to star-formation. △ Less

Submitted 9 December, 2021; v1 submitted 6 August, 2021; originally announced August 2021.

Comments: 17 pages, 11 figures, 2 tables, Accepted for publication in MNRAS

Redshift space three-point correlation function of IGM at $z<0.48$

Authors: Soumak Maitra, Raghunathan Srianand, Prakash Gaikwad, Nishikanta Khandai

Abstract: The Ly$α$ forest decomposed into Voigt profile components allow us to study clustering properties of the intergalactic medium and its dependence on various physical quantities. Here, we report the first detections of probability excess of low-z (i.e z<0.48) Ly$α$ absorber triplets over a scale of $r_\parallel\leq8$pMpc with a maximum amplitude of $8.76^{+1.96}_{-1.65}$ at a longitudinal separation… ▽ More The Ly$α$ forest decomposed into Voigt profile components allow us to study clustering properties of the intergalactic medium and its dependence on various physical quantities. Here, we report the first detections of probability excess of low-z (i.e z<0.48) Ly$α$ absorber triplets over a scale of $r_\parallel\leq8$pMpc with a maximum amplitude of $8.76^{+1.96}_{-1.65}$ at a longitudinal separation of 1-2pMpc. We measure non-zero three-point correlation ($ζ=4.76^{+1.98}_{-1.67}$) only at this scale with reduced three-point correlation value of Q=$0.95^{+0.39}_{-0.38}$. The measured $ζ$ shows an increasing trend with increasing HI column density ($N_{HI}$) while Q does not show any $N_{HI}$ dependence. About 88% of the triplets contributing to $ζ$ (at z$\le0.2$) have nearby galaxies (whose distribution is known to be complete for 0.1L$_*$ at z<0.1 and for L$_*$ at z~0.25 and within 20' to the quasar sightlines) within a velocity separation of 500$kms^{-1}$ and a median impact parameter of 405pkpc. The measured impact parameters are consistent with appreciable number of triplets at z$\le0.2$ not originating from individual galaxies but tracing the underlying galaxy distribution. Frequency of occurrence of high-b absorbers in triplets (~85%) is a factor~3 higher than that found among the full sample (~32%). Using four different cosmological simulations, we quantify the effect of peculiar velocities, feedback effects and show that most of the observed trends are broadly reproduced. However, $ζ$ at small scales ($r_\parallel<1$pMpc) and b-dependence of $ζ$ in simulations are found inconsistent with the observations. This could either be related to the fact that none of these simulations reproduce the observed b-distribution and $N_{HI}$ distribution for $N_{HI}>10^{14}$cm$^{-2}$ self-consistently or to the widespread of signal-to-noise ratio in the observed data. △ Less

Submitted 13 July, 2021; v1 submitted 10 December, 2020; originally announced December 2020.

Comments: 23 pages, 14 figures. Submitted to MNRAS

Simulated X-ray Emission in Galaxy Clusters with Feedback from Active Galactic Nuclei

Authors: Rudrani Kar Chowdhury, Soumya Roy, Suchetana Chatterjee, Nishikanta Khandai, Craig L. Sarazin, Tiziana Di Matteo

Abstract: To investigate the effect of feedback from active galactic nuclei (AGN) on their surrounding medium, we study the diffuse X-ray emission from galaxy groups and clusters by coupling the Astrophysical Plasma Emission Code (APEC) with the cosmological hydrodynamic simulation involving AGN feedback. We construct a statistical sample of synthetic Chandra X-ray photon maps to observationally characteriz… ▽ More To investigate the effect of feedback from active galactic nuclei (AGN) on their surrounding medium, we study the diffuse X-ray emission from galaxy groups and clusters by coupling the Astrophysical Plasma Emission Code (APEC) with the cosmological hydrodynamic simulation involving AGN feedback. We construct a statistical sample of synthetic Chandra X-ray photon maps to observationally characterize the effect of AGN on the ambient medium. We show that AGN are effective in displacing the hot X-ray emitting gas from the centers of groups and clusters, and that these signatures remain evident in observations of the X-ray surface brightness profiles. △ Less

Submitted 11 November, 2020; originally announced November 2020.

Comments: Accepted for publication in Astronomische Nachrichten, proceedings of the 9th International Workshop on Astronomy and Relativistic Astrophysics

The Distribution of Neutral Hydrogen in the Color-Magnitude Plane of Galaxies

Authors: Saili Dutta, Nishikanta Khandai

Abstract: We present the conditional HI (neutral hydrogen) Mass Function (HIMF) conditioned on observed optical properties, $M_{\text{r}}$ ($r$-band absolute magnitude) and $C_{\text{ur}}$ ($u-r$ color), for a sample of 7709 galaxies from ALFALFA (40% data release - $α.40$) which overlaps with a common volume in SDSS DR7. Based on the conditional HIMF we find that the luminous red, luminous blue and faint b… ▽ More We present the conditional HI (neutral hydrogen) Mass Function (HIMF) conditioned on observed optical properties, $M_{\text{r}}$ ($r$-band absolute magnitude) and $C_{\text{ur}}$ ($u-r$ color), for a sample of 7709 galaxies from ALFALFA (40% data release - $α.40$) which overlaps with a common volume in SDSS DR7. Based on the conditional HIMF we find that the luminous red, luminous blue and faint blue populations dominate the total HIMF at the high-mass end, knee and the low-mass end respectively. We use the conditional HIMF to derive the underlying distribution function of $Ω_{\text{HI}}$ (HI density parameter), $p(Ω_{\text{HI}})$, in the color-magnitude plane of galaxies. The distribution, $p(Ω_{\text{HI}})$, peaks in the blue cloud at $M_{\text{r}}^{\text{max}}=$ $-19.25, C_{\text{ur}}^{\text{max}}=1.44$ but is skewed. It has a long tail towards faint blue galaxies and luminous red galaxies. We argue that $p(Ω_{\text{HI}})$ can be used to reveal the underlying relation between cold gas, stellar mass and the star formation rate (SFR) in an unbiased way; that is the derived relation does not suffer from survey or sample selection. △ Less

Submitted 28 October, 2020; originally announced October 2020.

Comments: 5 pages, 3 figures, accepted for publication in MNRAS Letters

Cosmological Simulation of Galaxy Groups and Clusters-I: Global Effect of Feedback from Active Galactic Nuclei

Authors: Rudrani Kar Chowdhury, Suchetana Chatterjee, Anto . I. Lonappan, Nishikanta Khandai, Tiziana DiMatteo

Abstract: In this study we quantify the properties of the gas and dark matter around active galactic nuclei (AGN) in simulated galaxy groups and clusters and analyze the effect of AGN feedback on the surrounding intra-cluster (group) medium. Our results suggest downsizing of AGN luminosity with host halo mass, supporting the results obtained from clustering studies of AGN. By examining the temperature and d… ▽ More In this study we quantify the properties of the gas and dark matter around active galactic nuclei (AGN) in simulated galaxy groups and clusters and analyze the effect of AGN feedback on the surrounding intra-cluster (group) medium. Our results suggest downsizing of AGN luminosity with host halo mass, supporting the results obtained from clustering studies of AGN. By examining the temperature and density distribution of the gas in the vicinity of AGN we show that due to feedback from the central engine, the gas gets displaced from the centre of the group/cluster resulting in a reduction of the density but an enhancement of temperature. We show that these effects are pronounced at both high and low redshifts and propose new observables to study the effect of feedback in higher redshift galaxies. We also show that the average stellar mass is decreased in halos in the presence of AGN feedback confirming claims from previous studies. Our work for the first time uses a fully cosmological-hydrodynamic simulation to evaluate the global effects of AGN feedback on their host dark matter halos as well as galaxies at scales of galaxy groups and clusters. △ Less

Submitted 18 November, 2019; originally announced November 2019.

Comments: 11 pages, 8 figures, Accepted for publication at Astrophysical Journal

Journal ref: The Astrophysical Journal, Volume 889, Number 1, Year 2020

The Population of Galaxies that Contribute to The HI Mass Function

Authors: Saili Dutta, Nishikanta Khandai, Biprateep Dey

Abstract: We look at the contribution of different galaxy populations to the atomic hydrogen (HI) mass function (HIMF) and the HI density parameter, $Ω_{\text{HI}}$, in the local Universe. Our analysis is based on a sample of 7857 HI-selected galaxies selected from a volume common to the SDSS and ALFALFA surveys (40$\%$ catalog -- $α.40$). We define different populations of galaxies in the color(u-r)-magnit… ▽ More We look at the contribution of different galaxy populations to the atomic hydrogen (HI) mass function (HIMF) and the HI density parameter, $Ω_{\text{HI}}$, in the local Universe. Our analysis is based on a sample of 7857 HI-selected galaxies selected from a volume common to the SDSS and ALFALFA surveys (40$\%$ catalog -- $α.40$). We define different populations of galaxies in the color(u-r)-magnitude($M_{\text{r}}$) plane and compute the HIMF for each of them. Additionally we compute the HIMF for dark galaxies; these are undetected in SDSS and represent $\sim 2\%$ of the total sample. We find that the luminous red population dominates the total HIMF for $\log_{10}(M_{\text{HI}}h^2_{70}/M_{\odot}) \geq 10.4$. The full red population -- luminous and faint -- represents about $\sim 17\%$ of the $Ω_{\text{HI}}$ budget, while that of the dark population is $\sim 3\%$. The HIMF about the knee, $\log_{10}(M_{\text{HI}}h^2_{70}/M_{\odot}) \in [8,10.4]$, is dominated by the faint and luminous blue populations, the latter dominating at larger masses in this interval. Their total contribution to $Ω_{\text{HI}}$ is $\sim 55-70\%$, the variation depending on the definition of population. The dominant populations at the low mass end, $\log_{10}(M_{\text{HI}}h^2_{70}/M_{\odot}) \leq 8.0$ are the faint blue and faint bluer populations, the latter's dominance being sensitive to its definition. The full blue (blue--bluer luminous and faint) population represents $\sim 80\%$ of $Ω_{\text{HI}}$. A bimodal HIMF suggested by our results is however not seen since the amplitude of the HIMF of the luminous red population is small compared to that of the luminous blue population. △ Less

Submitted 25 March, 2020; v1 submitted 6 September, 2019; originally announced September 2019.

Comments: 17 pages, 10 figures, 3 tables, Accepted in MNRAS

Large-scale 3D mapping of the intergalactic medium using the Lyman Alpha Forest

Authors: Melih Ozbek, Rupert A. C. Croft, Nishikanta Khandai

Abstract: Maps of the large-scale structure of the Universe at redshifts 2-4 can be made with the Lyman-alpha forest which are complementary to low redshift galaxy surveys. We apply the Wiener interpolation method of Caucci et al. to construct three-dimensional maps from sets of Lyman-alpha forest spectra taken from cosmological hydrodynamic simulations. We mimic some current and future quasar redshift surv… ▽ More Maps of the large-scale structure of the Universe at redshifts 2-4 can be made with the Lyman-alpha forest which are complementary to low redshift galaxy surveys. We apply the Wiener interpolation method of Caucci et al. to construct three-dimensional maps from sets of Lyman-alpha forest spectra taken from cosmological hydrodynamic simulations. We mimic some current and future quasar redshift surveys (BOSS, eBOSS and MS-DESI) by choosing similar sightline densities. We use these appropriate subsets of the Lyman-alpha absorption sightlines to reconstruct the full three dimensional Lyman-alpha flux field and perform comparisons between the true and the reconstructed fields. We study global statistical properties of the intergalactic medium (IGM) maps with auto-correlation and cross-correlation analysis, slice plots, local peaks and point by point scatter. We find that both the density field and the statistical proper- ties of the IGM are recovered well enough that the resulting IGM maps can be meaningfully considered to represent large-scale maps of the Universe in agreement with Caucci et al., on larger scales and for sparser sightlines than had been tested previously. Quantitatively, for sightline parameters comparable to current and near future surveys the correlation coefficient between true and reconstructed fields is r > 0.9 on scales > 30 h^-1 Mpc. The properties of the maps are relatively insensitive to the precise form of the covariance matrix used. The final BOSS quasar Lyman-alpha forest sample will allow maps to be made with a resolution of ~ 30 h^-1 Mpc over a volume of ~ 15 h^-3 Gpc^3 between redshifts 1.9 and 2.3. △ Less

Submitted 4 February, 2016; originally announced February 2016.

Comments: 14 pages, 13 figures, 4 tables

Journal ref: MNRAS (March 11, 2016) 456 (4): 3610-3623

Luminosity function of [OII] emission-line galaxies in the MassiveBlack-II simulation

Authors: KwangHo Park, Tiziana Di Matteo, Shirley Ho, Rupert Croft, Stephen M. Wilkins, Yu Feng, Nishikanta Khandai

Abstract: We examine the luminosity function (LF) of [OII] emission-line galaxies in the high-resolution cosmological simulation MassiveBlack-II (MBII). From the spectral energy distribution of each galaxy, we select a sub-sample of star-forming galaxies at $0.06 \le z \le 3.0$ using the [OII] emission line luminosity L([OII]). We confirm that the specific star formation rate matches that in the GAMA survey… ▽ More We examine the luminosity function (LF) of [OII] emission-line galaxies in the high-resolution cosmological simulation MassiveBlack-II (MBII). From the spectral energy distribution of each galaxy, we select a sub-sample of star-forming galaxies at $0.06 \le z \le 3.0$ using the [OII] emission line luminosity L([OII]). We confirm that the specific star formation rate matches that in the GAMA survey. We show that the [OII] LF at z=1.0 from the MBII shows a good agreement with the LFs from several surveys below L([OII])=$10^{43.0}$ erg/s while the low redshifts ($z \le 0.3$) show an excess in the prediction of bright [OII] galaxies, but still displaying a good match with observations below L([OII])=$10^{41.6}$ erg/s. Based on the validity in reproducing the properties of [OII] galaxies at low redshift ($z \le 1$), we forecast the evolution of the [OII] LF at high redshift ($z \le 3$), which can be tested by upcoming surveys such as the HETDEX and DESI. The slopes of the LFs at bright and faint ends range from -3 to -2 showing minima at z=2. The slope of the bright end evolves approximately as 1/(z+1) at z=2 while the faint end evolves as ~3/(z+1) at $0.6 \le z \le 2$. In addition, a similar analysis is applied for the evolution of [OIII] LFs, which is to be explored in the forthcoming survey WFIRST-AFTA. Finally, we show that the auto-correlation function of [OII] and [OIII] emitting galaxies shows a rapid evolution from z=2 to 1. △ Less

Submitted 20 August, 2015; originally announced August 2015.

Comments: 12 pages, 13 figures, 3 tables, accepted for publication in MNRAS

Galaxy shapes and alignments in the MassiveBlack-II hydrodynamic and dark matter-only simulations

Authors: Ananth Tenneti, Rachel Mandelbaum, Tiziana Di Matteo, Alina Kiessling, Nishikanta Khandai

Abstract: We compare the shapes and intrinsic alignments of galaxies in the MassiveBlack-II cosmological hydrodynamic simulation (MBII) to those in a dark matter-only (DMO) simulation performed with the same volume (100$h^{-1}$Mpc)$^{3}$, cosmological parameters, and initial conditions. Understanding the impact of baryonic physics on galaxy shapes and alignments and their relation to the dark matter distrib… ▽ More We compare the shapes and intrinsic alignments of galaxies in the MassiveBlack-II cosmological hydrodynamic simulation (MBII) to those in a dark matter-only (DMO) simulation performed with the same volume (100$h^{-1}$Mpc)$^{3}$, cosmological parameters, and initial conditions. Understanding the impact of baryonic physics on galaxy shapes and alignments and their relation to the dark matter distribution should prove useful to map the intrinsic alignments of galaxies from hydrodynamic to dark matter-only simulations. We find that dark matter subhalos are typically rounder in MBII, and the shapes of stellar matter in low mass galaxies are more misaligned with the shapes of the dark matter of the corresponding subhalos in the DMO simulation. At $z=0.06$, the fractional difference in the mean misalignment angle between MBII and DMO simulations varies from $\sim 28 \% - 12 \%$ in the mass range $10^{10.8} - 6.0 \times 10^{14} h^{-1}M_{\odot}$. We study the dark matter halo shapes and alignments as a function of radius, and find that while galaxies in MBII are more aligned with the inner parts of their dark matter subhalos, there is no radial trend in their alignments with the corresponding subhalo in the DMO simulation. This result highlights the importance of baryonic physics in determining the alignment of the galaxy with respect to the inner parts of the halo. Finally, we compare the ellipticity-direction (ED) correlation for galaxies to that for dark matter halos, finding that it is suppressed on all scales by stellar-dark matter misalignment. In the projected shape-density correlation ($w_{δ+}$), which includes ellipticity weighting, this effect is partially canceled by the higher mean ellipticities of the stellar component, but differences of order $30-40\%$ remain on scales $> 1$ Mpc over a range of subhalo masses, with scale-dependent effects below $1$ Mpc. △ Less

Submitted 28 May, 2015; v1 submitted 12 May, 2015; originally announced May 2015.

Comments: 15 pages, 17 figures, submitted to MNRAS, reference added

Intrinsic alignments of galaxies in the MassiveBlack-II simulation: analysis of two-point statistics

Authors: Ananth Tenneti, Sukhdeep Singh, Rachel Mandelbaum, Tiziana Di Matteo, Yu Feng, Nishikanta Khandai

Abstract: The intrinsic alignment of galaxies with the large-scale density field is an important astrophysical contaminant in upcoming weak lensing surveys. We present detailed measurements of the galaxy intrinsic alignments and associated ellipticity-direction (ED) and projected shape ($w_{g+}$) correlation functions for galaxies in the cosmological hydrodynamic MassiveBlack-II (MB-II) simulation. We caref… ▽ More The intrinsic alignment of galaxies with the large-scale density field is an important astrophysical contaminant in upcoming weak lensing surveys. We present detailed measurements of the galaxy intrinsic alignments and associated ellipticity-direction (ED) and projected shape ($w_{g+}$) correlation functions for galaxies in the cosmological hydrodynamic MassiveBlack-II (MB-II) simulation. We carefully assess the effects on galaxy shapes, misalignment of the stellar component with the dark matter shape and two-point statistics of iterative weighted (by mass and luminosity) definitions of the (reduced and unreduced) inertia tensor. We find that iterative procedures must be adopted for a reliable measurement of the reduced tensor but that luminosity versus mass weighting has only negligible effects. Both ED and $w_{g+}$ correlations increase in amplitude with subhalo mass (in the range of $10^{10} - 6.0\times 10^{14}h^{-1}M_{\odot}$), with a weak redshift dependence (from $z=1$ to $z=0.06$) at fixed mass. At $z \sim 0.3$, we predict a $w_{g+}$ that is in reasonable agreement with SDSS LRG measurements and that decreases in amplitude by a factor of $\sim 5$--18 for galaxies in the LSST survey. We also compared the intrinsic alignments of centrals and satellites, with clear detection of satellite radial alignments within their host halos. Finally, we show that $w_{g+}$ (using subhalos as tracers of density) and $w_{δ+}$ (using dark matter density) predictions from the simulations agree with that of non-linear alignment models (NLA) at scales where the 2-halo term dominates in the correlations (and tabulate associated NLA fitting parameters). The 1-halo term induces a scale dependent bias at small scales which is not modeled in the NLA model. △ Less

Submitted 11 February, 2015; v1 submitted 25 September, 2014; originally announced September 2014.

Comments: 25 pages, 27 figures, revised after referee comments, accepted for publication in MNRAS

Galaxy Shapes and Intrinsic Alignments in The MassiveBlack-II Simulation

Authors: Ananth Tenneti, Rachel Mandelbaum, Tiziana Di Matteo, Yu Feng, Nishikanta Khandai

Abstract: The intrinsic alignment of galaxy shapes with the large-scale density field is a contaminant to weak lensing measurements, as well as being an interesting signature of galaxy formation and evolution (albeit one that is difficult to predict theoretically). Here we investigate the shapes and relative orientations of the stars and dark matter of halos and subhalos (central and satellite) extracted fr… ▽ More The intrinsic alignment of galaxy shapes with the large-scale density field is a contaminant to weak lensing measurements, as well as being an interesting signature of galaxy formation and evolution (albeit one that is difficult to predict theoretically). Here we investigate the shapes and relative orientations of the stars and dark matter of halos and subhalos (central and satellite) extracted from the MassiveBlack-II simulation, a state-of-the-art high resolution hydrodynamical cosmological simulation which includes stellar and AGN feedback in a volume of $(100{h^{-1}\mathrm{Mpc}})^3$. We consider redshift evolution from $z=1$ to $0.06$ and mass evolution within the range of subhalo masses, $10^{10} -6.0 \times 10^{14.0}{h^{-1}M_{\odot}}$. The shapes of the dark matter distributions are generally more round than the shapes defined by stellar matter. The projected root-mean-square (RMS) ellipticity per component for stellar matter is measured to be $e_\text{rms} = 0.28$ at $z=0.3$ for $M_{subhalo}> 10^{12.0}{h^{-1}M_{\odot}}$, which compares favourably with observational measurements. We find that the shapes of stellar and dark matter are more round for less massive subhalos and at lower redshifts. By directly measuring the relative orientation of the stellar matter and dark matter of subgroups, we find that, on average, the misalignment between the two components is larger for less massive subhalos. The mean misalignment angle varies from $\sim 30^{\circ}-10^{\circ}$ for $M \sim 10^{10} - 10^{14} {h^{-1}M_{\odot}}$ and shows a weak dependence on redshift. We also compare the misalignment angles in central and satellite subhalos at fixed subhalo mass, and find that centrals are more misaligned than satellites. We present fitting formulae for the shapes of dark and stellar matter in subhalos and also the probability distributions of misalignment angles. △ Less

Submitted 17 March, 2014; originally announced March 2014.

Comments: 18 pages, 18 figures, submitted to MNRAS

Journal ref: 2014, MNRAS, 441(1): 470-485

The MassiveBlack-II Simulation: The Evolution of Halos and Galaxies to z~0

Authors: Nishikanta Khandai, Tiziana Di Matteo, Rupert Croft, Stephen M. Wilkins, Yu Feng, Evan Tucker, Colin DeGraf, Mao-Sheng Liu

Abstract: (Abridged for arXiv)We investigate the properties of halos, galaxies and blackholes to z=0 in the high resolution hydrodynamical simulation MassiveBlack-II (MBII) which evolves a LCDM cosmology in a comoving volume Vbox=100(Mpc/h)^3. MBII is the highest resolution simulation of this size which includes a self-consistent model for star formation, black hole accretion and associated feedback. We pro… ▽ More (Abridged for arXiv)We investigate the properties of halos, galaxies and blackholes to z=0 in the high resolution hydrodynamical simulation MassiveBlack-II (MBII) which evolves a LCDM cosmology in a comoving volume Vbox=100(Mpc/h)^3. MBII is the highest resolution simulation of this size which includes a self-consistent model for star formation, black hole accretion and associated feedback. We provide a simulation browser web application which enables interactive search and tagging of halos, subhalos and their properties and publicly release our galaxy catalogs. Our analysis of the halo mass function (MF) in MBII reveals that baryons have strong effects, with changes in the halo abundance of 20-35% below the knee of the MF (Mhalo < 10^13.2 Msun/h at z=0) when compared to fits based on dark matter only simulations. We provide a fitting function for the halo MF out to redshift z=11 and discuss how the onset of non-universality in the MF limits the accuracy of our fit. We study the halo occupation distribution and clustering of galaxies, in particular the evolution and scale dependence of stochasticity and bias finding reasonable agreement with observational data. The shape of the cosmic spectral energy distribution predicted by MBII is consistent with observations, but lower in amplitude. The Galaxy Stellar Mass Function (GSMF) function is broadly consistent with observations at z>=2. At z<2, the population of passive low mass (for M*<10^9 Msun) galaxies in MBII makes the GSMF too steep compared to observations whereas at the high mass end (M*>10^11 Msun) galaxies hosting bright AGNs make significant contributions to the GSMF. The quasar bolometric luminosity function is also largely consistent with observations. We note however that more efficient AGN feedback (beyond simple thermal coupling used here) is likely necessary for the largest, rarest objects/clusters at low redshifts. △ Less

Submitted 4 February, 2014; originally announced February 2014.

Comments: 26 pages, 25 figures. Submitted to MNRAS. High-resolution version and MBII galaxy catalogs can be found at http://mbii.phys.cmu.edu/data/

Nonparametric 3D map of the IGM using the Lyman-alpha forest

Authors: Jessi Cisewski, Rupert A. C. Croft, Peter E. Freeman, Christopher R. Genovese, Nishikanta Khandai, Melih Ozbek, Larry Wasserman

Abstract: Visualizing the high-redshift Universe is difficult due to the dearth of available data; however, the Lyman-alpha forest provides a means to map the intergalactic medium at redshifts not accessible to large galaxy surveys. Large-scale structure surveys, such as the Baryon Oscillation Spectroscopic Survey (BOSS), have collected quasar (QSO) spectra that enable the reconstruction of HI density fluct… ▽ More Visualizing the high-redshift Universe is difficult due to the dearth of available data; however, the Lyman-alpha forest provides a means to map the intergalactic medium at redshifts not accessible to large galaxy surveys. Large-scale structure surveys, such as the Baryon Oscillation Spectroscopic Survey (BOSS), have collected quasar (QSO) spectra that enable the reconstruction of HI density fluctuations. The data fall on a collection of lines defined by the lines-of-sight (LOS) of the QSO, and a major issue with producing a 3D reconstruction is determining how to model the regions between the LOS. We present a method that produces a 3D map of this relatively uncharted portion of the Universe by employing local polynomial smoothing, a nonparametric methodology. The performance of the method is analyzed on simulated data that mimics the varying number of LOS expected in real data, and then is applied to a sample region selected from BOSS. Evaluation of the reconstruction is assessed by considering various features of the predicted 3D maps including visual comparison of slices, PDFs, counts of local minima and maxima, and standardized correlation functions. This 3D reconstruction allows for an initial investigation of the topology of this portion of the Universe using persistent homology. △ Less

Submitted 8 January, 2014; originally announced January 2014.

High redshift supermassive blackholes: accretion through cold flows

Authors: Yu Feng, Tiziana Di Matteo, Rupert Croft, Nishikanta Khandai

Abstract: We use zoom-in techniques to re-simulate three high-redshift (z > 5.5) halos which host 10^9 solar mass blackholes from the ~ Gpc volume, MassiveBlack cosmological hydrodynamic simulation. We examine a number of factors potentially affecting supermassive blackhole growth at high redshift in cosmological simulations. These include numerical resolution, feedback prescriptions and formulation of smoo… ▽ More We use zoom-in techniques to re-simulate three high-redshift (z > 5.5) halos which host 10^9 solar mass blackholes from the ~ Gpc volume, MassiveBlack cosmological hydrodynamic simulation. We examine a number of factors potentially affecting supermassive blackhole growth at high redshift in cosmological simulations. These include numerical resolution, feedback prescriptions and formulation of smoothed particle hydrodynamics. We find that varying the size of the region over which feedback energy is deposited directly, either for fixed number of neighbours or fixed volume makes very little difference to the accretion history of blackholes. Changing mass resolution by factors of up to 64 also does not change the blackhole growth history significantly. We find that switching from the density-entropy formulation to the pressure-entropy formulation of smoothed particle hydrodynamics slightly increases the accretion rate onto blackholes. In general numerical details appear to have small effects on the main fueling mechanism for blackholes at these high redshifts. We examine the fashion by which this occurs, finding that the insensitivity to simulation technique seems to be a hallmark of the cold flow feeding picture of these high-z supermassive blackholes. We show that the gas that participates in critical accretion phases, in these massive objects at z > 6~7 is in all cases colder, denser, and forms more coherent streams than the average gas in the halo. This is also mostly the case when the blackhole accretion is feedback regulated (z < 6), however the distinction is less prominent. For our resimulated halos, cold flows appear to be a viable mechanism for forming the most massive blackholes in the early universe, occurring naturally in LambdaCDM models of structure formation. Not requiring fine tuning of numerical parameters, they seem to be physically inevitable in these objects. △ Less

Submitted 4 December, 2013; originally announced December 2013.

Comments: 15 pages, 12 figures

Theoretical predictions for the effect of nebular emission on the broad band photometry of high-redshift galaxies

Authors: Stephen M. Wilkins, William Coulton, Joseph Caruana, Rupert Croft, Tiziana Di Matteo, Nishikanta Khandai, Yu Feng, Andrew Bunker, Holly Elbert

Abstract: By combining optical and near-IR observations from the Hubble Space Telescope with NIR photometry from the Spitzer Space Telescope it is possible to measure the rest-frame UV-optical colours of galaxies at z=4-8. The UV-optical spectral energy distribution of star formation dominated galaxies is the result of several different factors. These include the joint distribution of stellar masses, ages,… ▽ More By combining optical and near-IR observations from the Hubble Space Telescope with NIR photometry from the Spitzer Space Telescope it is possible to measure the rest-frame UV-optical colours of galaxies at z=4-8. The UV-optical spectral energy distribution of star formation dominated galaxies is the result of several different factors. These include the joint distribution of stellar masses, ages, and metallicities, and the subsequent reprocessing by dust and gas in the ISM. Using a large cosmological hydrodynamical simulation we investigate the predicted spectral energy distributions of galaxies at high-redshift with a particular emphasis on assessing the potential contribution of nebular emission. We find that the average pure stellar UV-optical colour correlates with both luminosity and redshift such that galaxies at lower-redshift and higher-luminosity are typically redder. Assuming the escape fraction of ionising photons is close to zero, the effect of nebular emission is to redden the UV-optical 1500-V_w colour by, on average, 0.4 mag at z=8 declining to 0.25 mag at z=4. Young and low-metallicity stellar populations, which typically have bluer pure stellar UV-optical colours, produce larger ionising luminosities and are thus more strongly affected by the reddening effects of nebular emission. This causes the distribution of 1500-V_w colours to narrow and the trends with luminosity and redshift to weaken. The strong effect of nebular emission leaves observed-frame colours critically sensitive to the source redshift. For example, increasing the redshift by 0.1 can result in observed frame colours changing by up to ~0.6. These predictions reinforce the need to include nebular emission when modelling the spectral energy distributions of galaxies at high-redshift and also highlight the difficultly in interpreting the observed colours of individual galaxies without precise redshifts. △ Less

Submitted 28 August, 2013; originally announced August 2013.

Comments: 13 pages, 17 figures, accepted in MNRAS

Interpreting the observed UV continuum slopes of high-redshift galaxies

Authors: Stephen M. Wilkins, Andrew Bunker, William Coulton, Rupert Croft, Tiziana Di Matteo, Nishikanta Khandai, Yu Feng

Abstract: The observed UV continuum slope of star forming galaxies is strongly affected by the presence of dust. Its observation is then a potentially valuable diagnostic of dust attenuation, particularly at high-redshift where other diagnostics are currently inaccesible. Interpreting the observed UV continuum slope in the context of dust attenuation is often achieved assuming the empirically calibrated Meu… ▽ More The observed UV continuum slope of star forming galaxies is strongly affected by the presence of dust. Its observation is then a potentially valuable diagnostic of dust attenuation, particularly at high-redshift where other diagnostics are currently inaccesible. Interpreting the observed UV continuum slope in the context of dust attenuation is often achieved assuming the empirically calibrated Meurer et al. (1999) relation. Implicit in this relation is the assumption of an intrinsic UV continuum slope ($β=-2.23$). However, results from numerical simulations suggest that the intrinsic UV continuum slopes of high-redshift star forming galaxies are bluer than this, and moreover vary with redshift. Using values of the intrinsic slope predicted by numerical models of galaxy formation combined with a Calzetti et al. (2000) reddening law we infer UV attenuations ($A_{1500}$) $0.35-0.5\,{\rm mag}$ ($A_{V}$: $0.14-0.2\,{\rm mag}$ assuming Calzetti et al. 2000) greater than simply assuming the Meurer relation. This has significant implications for the inferred amount of dust attenuation at very-high ($z\approx 7$) redshift given current observational constraints on $β$, combined with the Meurer relation, suggest dust attenuation to be virtually zero in all but the most luminous systems. △ Less

Submitted 6 February, 2013; originally announced February 2013.

Comments: 6 pages, 4 figures, accepted in MNRAS

Confronting predictions of the galaxy stellar mass function with observations at high-redshift

Authors: Stephen M. Wilkins, Tiziana Di Matteo, Rupert Croft, Nishikanta Khandai, Yu Feng, Andrew Bunker, William Coulton

Abstract: We investigate the evolution of the galaxy stellar mass function at high-redshift ($z\ge 5$) using a pair of large cosmological hydrodynamical simulations: {\em MassiveBlack} and {\em MassiveBlack-II}. By combining these simulations we can study the properties of galaxies with stellar masses greater than $10^{8}\,{\rm M_{\odot}}\,h^{-1}$ and (co-moving) number densities of… ▽ More We investigate the evolution of the galaxy stellar mass function at high-redshift ($z\ge 5$) using a pair of large cosmological hydrodynamical simulations: {\em MassiveBlack} and {\em MassiveBlack-II}. By combining these simulations we can study the properties of galaxies with stellar masses greater than $10^{8}\,{\rm M_{\odot}}\,h^{-1}$ and (co-moving) number densities of $\log_{10}(φ\, [{\rm Mpc^{-3}\,dex^{-1}}\,h^{3}])>-8$. Observational determinations of the galaxy stellar mass function at very-high redshift typically assume a relation between the observed UV luminosity and stellar mass-to-light ratio which is applied to high-redshift samples in order to estimate stellar masses. This relation can also be measured from the simulations. We do this, finding two significant differences with the usual observational assumption: it evolves strongly with redshift and has a different shape. Using this relation to make a consistent comparison between galaxy stellar mass functions we find that at $z=6$ and above the simulation predictions are in good agreement with observed data over the whole mass range. Without using the correct UV luminosity and stellar mass-to-light ratio, the discrepancy would be up to two orders of magnitude for large galaxies $>10^{10}\,{\rm M_{\odot}}\,h^{-1}$. At $z=5$, however the stellar mass function for low mass $<10^{9}\,{\rm M_{\odot}}\,h^{-1}$ galaxies is overpredicted by factors of a few, consistent with the behaviour of the UV luminosity function, and perhaps a sign that feedback in the simulation is not efficient enough for these galaxies. △ Less

Submitted 12 January, 2013; originally announced January 2013.

Comments: 7 pages, 4 figures, accepted in MNRAS

Growth and anisotropy of ionization fronts near high redshift quasars in the MassiveBlack simulation

Authors: Yu Feng, Rupert A. C. Croft, Tiziana Di Matteo, Nishikanta Khandai

Abstract: We use radiative transfer to study the growth of ionized regions around the brightest, z=8 quasars in a large cosmological hydrodynamic simulation that includes black hole growth and feedback (the MassiveBlack simulation). We find that in the presence of the quasar s the comoving HII bubble radii reach 10 Mpc/h after 20 My while with the stellar component alone the HII bubbles are smaller by at le… ▽ More We use radiative transfer to study the growth of ionized regions around the brightest, z=8 quasars in a large cosmological hydrodynamic simulation that includes black hole growth and feedback (the MassiveBlack simulation). We find that in the presence of the quasar s the comoving HII bubble radii reach 10 Mpc/h after 20 My while with the stellar component alone the HII bubbles are smaller by at least an order of magnitude. Our calculations show that several features are not captured within an analytical growth model of Stromgren spheres. The X-ray photons from hard quasar spectra drive a smooth transition from fully neutral to partially neutral in the ionization front. However the transition from partially neutral to fully ionized is significantly more complex. We measure the distance to the edge of bubbles as a function of angle and use the standard deviation of these distances as a diagnostic of the isotropy of ionized regions. We find that the overlapping of nearby ionized regions from clustered halos not only increases the anisotropy, but also is the main mechanism which allows the outer radius to grow. We therefore predict that quasar ionized bubbles at this early stage in the reionization process should be both significantly larger and more irregularly shaped than bubbles around star forming galaxies. Before the star formation rate increases and the Universe fully reionizes, quasar bubbles will form the most striking and recognizable features in 21cm maps. △ Less

Submitted 5 November, 2012; v1 submitted 31 August, 2012; originally announced August 2012.

Comments: 11 pages, 10 figures. Updated after referee report

Growth of Early Supermassive Black Holes and the High-Redshift Eddington Ratio Distribution

Authors: Colin DeGraf, Tiziana Di Matteo, Nishikanta Khandai, Rupert Croft

Abstract: Using a new large-scale (~ 0.75 Gpc)^3 hydrodynamic cosmological simulation we investigate the growth rate of supermassive black holes in the early universe (z > 4.75). Remarkably, we find a clear peak in the typical Eddington ratio at black hole masses of 4-8 * 10^7 solar masses (typically found in halos of ~7 * 10^11 to 10^12 solar masses), independent of redshift and indicative that most of BH… ▽ More Using a new large-scale (~ 0.75 Gpc)^3 hydrodynamic cosmological simulation we investigate the growth rate of supermassive black holes in the early universe (z > 4.75). Remarkably, we find a clear peak in the typical Eddington ratio at black hole masses of 4-8 * 10^7 solar masses (typically found in halos of ~7 * 10^11 to 10^12 solar masses), independent of redshift and indicative that most of BH growth occurs in the cold-flow dominated regime. Black hole growth is by and large regulated by the evolution of gas density. The typical Eddington ratio at a given mass scales simply as cosmological density (1+z)^3 and the peak is caused by the competition between increased gas density available in more massive hosts, and a decrease due to strong AGN feedback that deprives the black hole of sufficient gas to fuel further rapid growth in the high mass end. In addition to evolution in the mean Eddington ratio, we show that the distribution of Eddington ratio among both mass-selected and luminosity-selected samples is approximately log-normal. We combine these findings into a single log-normal fitting formula for the distribution of Eddington ratios as a function of (M_BH,z). This formula can be used in analytic and semi analytic models for evolving black hole populations, predicting black hole masses of observed quasars, and, in conjunction with the observed distribution of Eddington ratios, can be used to constrain the black hole mass function. △ Less

Submitted 25 January, 2012; originally announced January 2012.

Comments: 6 pages, 5 figures

The Formation of Galaxies Hosting z~6 Quasars

Authors: Nishikanta Khandai, Yu Feng, Colin DeGraf, Tiziana Di Matteo, Rupert A. C. Croft

Abstract: We investigate the formation and properties of galaxies hosting z~6 quasars, in the gigaparsec scale cosmological hydrodynamical simulation: MassiveBlack, which includes a self-consistent model for star formation, black hole accretion and associated feedback. We show that the MassiveBlack reproduces current estimates of the galaxy stellar mass function z=5, 6. We find that quasar hosts in the simu… ▽ More We investigate the formation and properties of galaxies hosting z~6 quasars, in the gigaparsec scale cosmological hydrodynamical simulation: MassiveBlack, which includes a self-consistent model for star formation, black hole accretion and associated feedback. We show that the MassiveBlack reproduces current estimates of the galaxy stellar mass function z=5, 6. We find that quasar hosts in the simulation are compact gas rich systems with high star formations rates of SFR ~ 100-1000 Msun/yr consistent with observed properties of Sloan quasar hosts in the redshift range 5.5 < z < 6.5. We show that the star-forming gas in these galaxies predominantly originates from high density cold streams which are able to penetrate the halo and grow the galaxy at the center. MassiveBlack predicts a deviation from the local Mbh-sigma and Mbh-Mstar relation implying that black holes are relatively more massive for a given stellar host at these redshifts. △ Less

Submitted 2 November, 2011; originally announced November 2011.

Comments: 12 pages, 8 figures

Dark matter halo occupation: environment and clustering

Authors: Rupert Croft, Tiziana Di Matteo, Nishikanta Khandai, Volker Springel, Anirban Jana, Jeffrey Gardner

Abstract: We use a large dark matter simulation of a LambdaCDM model to investigate the clustering and environmental dependence of the number of substructures in a halo. Focusing on redshift z=1, we find that the halo occupation distribution is sensitive at the tens of percent level to the surrounding density and to a lesser extent to asymmetry of the surrounding density distribution. We compute the autocor… ▽ More We use a large dark matter simulation of a LambdaCDM model to investigate the clustering and environmental dependence of the number of substructures in a halo. Focusing on redshift z=1, we find that the halo occupation distribution is sensitive at the tens of percent level to the surrounding density and to a lesser extent to asymmetry of the surrounding density distribution. We compute the autocorrelation function of halos as a function of occupation, building on the finding of Wechsler et al. (2006) and Gao and White (2007) that halos (at fixed mass) with more substructure are more clustered. We compute the relative bias as a function of occupation number at fixed mass, finding a strong relationship. At fixed mass, halos in the top 5% of occupation can have an autocorrelation function ~ 1.5-2 times higher than the mean. We also compute the bias as a function of halo mass, for fixed halo occupation. We find that for group and cluster sized halos, when the number of subhalos is held fixed, there is a strong anticorrelation between bias and halo mass. Such a relationship represents an additional challenge to the halo model. △ Less

Submitted 19 September, 2011; originally announced September 2011.

Comments: 13 pages, 10 figures, MNRAS submitted

Terapixel imaging of cosmological simulations

Authors: Yu Feng, Rupert A. C. Croft, Tiziana Di Matteo, Nishikanta Khandai, Randy Sargent, Illah Nourbakhsh, Paul Dille, Chris Bartley, Volker Springel, Anirban Jana, Jeffrey Gardner

Abstract: The increasing size of cosmological simulations has led to the need for new visualization techniques. We focus on Smoothed Particle Hydrodynamical (SPH) simulations run with the GADGET code and describe methods for visually accessing the entire simulation at full resolution. The simulation snapshots are rastered and processed on supercomputers into images that are ready to be accessed through a we… ▽ More The increasing size of cosmological simulations has led to the need for new visualization techniques. We focus on Smoothed Particle Hydrodynamical (SPH) simulations run with the GADGET code and describe methods for visually accessing the entire simulation at full resolution. The simulation snapshots are rastered and processed on supercomputers into images that are ready to be accessed through a web interface (GigaPan). This allows any scientist with a web-browser to interactively explore simulation datasets in both in spatial and temporal dimensions, datasets which in their native format can be hundreds of terabytes in size or more. We present two examples, the first a static terapixel image of the MassiveBlack simulation, a P-GADGET SPH simulation with 65 billion particles, and the second an interactively zoomable animation of a different simulation with more than one thousand frames, each a gigapixel in size. Both are available for public access through the GigaPan web interface. We also make our imaging software publicly available. △ Less

Submitted 6 July, 2011; originally announced July 2011.

Comments: 9 pages, 6 figures, Link to the Gigapan image: http://gigapan.org/gigapans/76215/ Link to the TimeMachine animation: http://timemachine.gigapan.org/wiki/Evolution_of_the_Universe Link to Gaepsi: http://web.phys.cmu.edu/~yfeng1/gaepsi

Early Black Holes in Cosmological Simulations: Luminosity Functions and Clustering Behaviour

Authors: Colin DeGraf, Tiziana Di Matteo, Nishikanta Khandai, Rupert Croft, Julio Lopez, Volker Springel

Abstract: We examine predictions for the quasar luminosity functions (QLF) and quasar clustering at high redshift (z > 4.75) using MassiveBlack, our new hydrodynamic cosmological simulation which includes a self-consistent model for black hole growth and feedback. We show that the model reproduces the Sloan QLF within observational constraints at z >= 5. We find that the high-z QLF is consistent with a reds… ▽ More We examine predictions for the quasar luminosity functions (QLF) and quasar clustering at high redshift (z > 4.75) using MassiveBlack, our new hydrodynamic cosmological simulation which includes a self-consistent model for black hole growth and feedback. We show that the model reproduces the Sloan QLF within observational constraints at z >= 5. We find that the high-z QLF is consistent with a redshift-independent occupation distribution of BHs among dark matter halos (which we provide) such that the evolution of the QLF follows that of the halo mass function. The sole exception is the bright-end at z=6 and 7, where BHs in high-mass halos tend to be unusually bright due to extended periods of Eddington growth caused by high density cold flows into the halo center. We further use these luminosity functions to make predictions for the number density of quasars in upcoming surveys, predicting there should be ~119+-28 (~87+-28) quasars detectable in the F125W band of the WIDE (DEEP) fields of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) from z=5-6, ~19+-7 (~18+-9) from z=6-7, and ~1.7+-1.5 (~1.5+-1.5) from z=7-8. We also investigate quasar clustering, finding that the correlation length is fully consistent with current constraints for Sloan quasars (r_0~17 h^{-1} Mpc at z=4 for quasars above m_i = 20.2), and grows slowly with redshift up to z=6 (r_0~22 h^{-1} Mpc). Finally, we note that the quasar clustering strength depends weakly on luminosity for low L_BH, but gets stronger at higher L_BH as the BHs are found in higher mass halos. △ Less

Submitted 6 July, 2011; originally announced July 2011.

Comments: 8 pages, 3 figures

Cold flows and the first quasars

Authors: Tiziana Di Matteo, Nishikanta Khandai, Colin DeGraf, Yu Feng, Rupert Croft, Julio Lopez, Volker Springel

Abstract: Observations of the most distant bright quasars imply that billion solar mass supermassive black holes (SMBH) have to be assembled within the first eight hundred million years. Under our standard galaxy formation scenario such fast growth implies large gas densities providing sustained accretion at critical or supercritical rates onto an initial black hole seed. It has been a long standing questio… ▽ More Observations of the most distant bright quasars imply that billion solar mass supermassive black holes (SMBH) have to be assembled within the first eight hundred million years. Under our standard galaxy formation scenario such fast growth implies large gas densities providing sustained accretion at critical or supercritical rates onto an initial black hole seed. It has been a long standing question whether and how such high black hole accretion rates can be achieved and sustained at the centers of early galaxies. Here we use our new cosmological hydrodynamic simulation (MassiveBlack) covering a volume (0.75 \Gpc)^3 appropriate for studying the rare first quasars to show that steady high density cold gas flows responsible for assembling the first galaxies produce the high gas densities that lead to sustained critical accretion rates and hence rapid growth commensurate with the existence of ~10^9 solar mass black holes as early as z~7. We find that under these conditions quasar feedback is not effective at stopping the cold gas from penetrating the central regions and hence cannot quench the accretion until the host galaxy reaches M_halo > 10^{12} solar masses. This cold-flow driven scenario for the formation of quasars implies that they should be ubiquitous in galaxies in the early universe and that major (proto)galaxy mergers are not a requirement for efficient fuel supply and growth, particularly for the earliest SMBHs. △ Less

Submitted 6 July, 2011; originally announced July 2011.

Comments: 6 pages, 4 figures

Detecting neutral hydrogen in emission at redshift z ~ 1

Authors: Nishikanta Khandai, Shiv K. Sethi, Tiziana Di Matteo, Rupert A. C. Croft, Volker Springel, Anirban Jana, Jeffrey P. Gardner

Abstract: We use a large N-body simulation to examine the detectability of HI in emission at redshift z ~ 1, and the constraints imposed by current observations on the neutral hydrogen mass function of galaxies at this epoch. We consider three different models for populating dark matter halos with HI, designed to encompass uncertainties at this redshift. These models are consistent with recent observations… ▽ More We use a large N-body simulation to examine the detectability of HI in emission at redshift z ~ 1, and the constraints imposed by current observations on the neutral hydrogen mass function of galaxies at this epoch. We consider three different models for populating dark matter halos with HI, designed to encompass uncertainties at this redshift. These models are consistent with recent observations of the detection of HI in emission at z ~ 0.8. Whilst detection of 21 cm emission from individual halos requires extremely long integrations with existing radio interferometers, such as the Giant Meter Radio Telescope (GMRT), we show that the stacked 21 cm signal from a large number of halos can be easily detected. However, the stacking procedure requires accurate redshifts of galaxies. We show that radio observations of the field of the DEEP2 spectroscopic galaxy redshift survey should allow detection of the HI mass function at the 5-12 sigma level in the mass range 10^(11.4) M_sun/h < M_halo < 10^(12.5)M_sun/h, with a moderate amount of observation time. Assuming a larger noise level that corresponds to an upper bound for the expected noise for the GMRT, the detection significance for the HI mass function is still at the 1.7-3 sigma level. We find that optically undetected satellite galaxies enhance the HI emission profile of the parent halo, leading to broader wings as well as a higher peak signal in the stacked profile of a large number of halos. We show that it is in principle possible to discern the contribution of undetected satellites to the total HI signal, even though cosmic variance limitation make this challenging for some of our models. △ Less

Submitted 14 December, 2010; v1 submitted 8 December, 2010; originally announced December 2010.

Comments: 14 pages, 9 figures, Submitted To MNRAS

Fractal Dimension as a measure of the scale of Homogeneity

Authors: Jaswant K. Yadav, J. S. Bagla, Nishikanta Khandai

Abstract: 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 achiev… ▽ More 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. △ Less

Submitted 15 March, 2010; v1 submitted 4 January, 2010; originally announced January 2010.

Comments: 7 pages, two figures. Accepted for publication in the MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.405:2009,2010

HI as a Probe of the Large Scale Structure in the Post-Reionization Universe: Visibility Correlations and Prospects for Detection

Authors: Nishikanta Khandai, Kanan K. Datta, J. S. Bagla

Abstract: 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 Wide… ▽ More 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. △ Less

Submitted 28 April, 2010; v1 submitted 26 August, 2009; originally announced August 2009.

Comments: This submission has been withdrawn by the authors. This manuscript has been merged with 0908.3796, and accepted for publication in the MNRAS.

HI as a Probe of the Large Scale Structure in the Post-Reionization Universe

Authors: J. S. Bagla, Nishikanta Khandai, Kanan K. Datta

Abstract: 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 t… ▽ More 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. △ Less

Submitted 28 April, 2010; v1 submitted 26 August, 2009; originally announced August 2009.

Comments: 15 pages, 11 figures. Accepted for publication in the MNRAS. This is a merged manuscript also containing material covered in 0908.3857

Journal ref: Mon.Not.Roy.Astron.Soc.407:567,2010

Mass function of haloes: scale invariant models

Authors: J. S. Bagla, Nishikanta Khandai, Girish Kulkarni

Abstract: 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 si… ▽ More 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. △ Less

Submitted 20 January, 2010; v1 submitted 19 August, 2009; originally announced August 2009.

Comments: 8 pages, four figures

The Adaptive TreePM: An Adaptive Resolution Code for Cosmological N-body Simulations

Authors: J. S. Bagla, Nishikanta Khandai

Abstract: 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 cosmo… ▽ More 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. △ Less

Submitted 27 May, 2009; v1 submitted 26 November, 2008; originally announced November 2008.

Comments: 18 pages, 11 figures. Accepted for publication in the MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.396:2211,2008

Effects of the size of cosmological N-Body simulations on physical quantities - III: Skewness

Authors: J. S. Bagla, Jayanti Prasad, Nishikanta Khandai

Abstract: 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 unreliabl… ▽ More 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. △ Less

Submitted 3 February, 2009; v1 submitted 8 April, 2008; originally announced April 2008.

Comments: 14 pages, 8 figures. Accepted for publication in the MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.395:918,2008

A Modified TreePM Code

Authors: Nishikanta Khandai, J. S. Bagla

Abstract: 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,… ▽ More 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. △ Less

Submitted 24 April, 2009; v1 submitted 21 February, 2008; originally announced February 2008.

Comments: 16 pages, 5 figures; Accepted for publication in Research In Astronomy and Astrophysics (RAA)