Source author record

Jun Koda

Jun Koda appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2025arXiv

Learning to detect continuous gravitational waves: an open data-analysis competition

We report results of a public data-analysis challenge, hosted on the open data-science platform Kaggle, to detect simulated continuous gravitational-wave signals (CWs). These are weak signals from rapidly spinning neutron stars that remain undetected despite extensive searches. The competition dataset consisted of a population of CW signals using both simulated and real LIGO detector data matching the conditions of actual CW searches. The competition attracted more than 1,000 participants to develop realistic CW search algorithms. We describe the top 10 approaches and discuss their applicability as a pre-processing step compared to standard CW-search approaches. For the competition's dataset, we find that top approaches can reduce the computing cost by 1 to 3 orders of magnitude at a false-dismissal probability comparable to standard CW searches. Additionally, the competition drove the development of new GPU-accelerated detection pipelines, which facilitated their adoption in other areas of gravitational-wave data analysis. We release the associated dataset, which constitutes the first open standardized benchmark for CW detection, to enable reproducible method comparisons and to encourage further developments toward the first detection of these elusive signals.

preprint2019arXiv

The impact of inhomogeneous subgrid clumping on cosmic reionization

Cosmic reionization was driven by the imbalance between early sources and sinks of ionizing radiation, both of which were dominated by small-scale structure and are thus usually treated in cosmological reionization simulations by subgrid modelling. The recombination rate of intergalactic hydrogen is customarily boosted by a subgrid clumping factor, ${\left<n^2\right>/\left<n\right>^2}$, which corrects for unresolved fluctuations in gas density ${n}$ on scales below the grid-spacing of coarse-grained simulations. We investigate in detail the impact of this inhomogeneous subgrid clumping on reionization and its observables, as follows: (1) Previous attempts generally underestimated the clumping factor because of insufficient mass resolution. We perform a high-resolution $N$-body simulation that resolves haloes down to the pre-reionization Jeans mass to derive the time-dependent, spatially-varying local clumping factor and a fitting formula for its correlation with local overdensity. (2) We then perform a large-scale $N$-body and radiative transfer simulation that accounts for this inhomogeneous subgrid clumping by applying this clumping factor-overdensity correlation. Boosting recombination significantly slows the expansion of ionized regions, which delays completion of reionization and suppresses 21 cm power spectra on large scales in the later stages of reionization. (3) We also consider a simplified prescription in which the globally-averaged, time-evolving clumping factor from the same high-resolution $N$-body simulation is applied uniformly to all cells in the reionization simulation, instead. Observables computed with this model agree fairly well with those from the inhomogeneous clumping model, e.g. predicting 21 cm power spectra to within 20% error, suggesting it may be a useful approximation.

preprint2016arXiv

The Detailed Science Case for the Maunakea Spectroscopic Explorer: the Composition and Dynamics of the Faint Universe

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

preprint2015arXiv

Fast and accurate mock catalogue generation for low-mass galaxies

We present an accurate and fast framework for generating mock catalogues including low-mass halos, based on an implementation of the COmoving Lagrangian Acceleration (COLA) technique. Multiple realisations of mock catalogues are crucial for analyses of large-scale structure, but conventional N-body simulations are too computationally expensive for the production of thousands of realisations. We show that COLA simulations can produce accurate mock catalogues with a moderate computation resource for low- to intermediate- mass galaxies in $10^{12} M_\odot$ haloes, both in real and redshift space. COLA simulations have accurate peculiar velocities, without systematic errors in the velocity power spectra for k < 0.15 h/Mpc, and with only 3-per-cent error for k < 0.2 h/Mpc. We use COLA with 10 time steps and a Halo Occupation Distribution to produce 600 mock galaxy catalogues of the WiggleZ Dark Energy Survey. Our parallelized code for efficient generation of accurate halo catalogues is publicly available at github.com/junkoda/cola_halo.

preprint2015arXiv

Searching for Modified Gravity: Scale and Redshift Dependent Constraints from Galaxy Peculiar Velocities

We present measurements of both scale- and time-dependent deviations from the standard gravitational field equations. These late-time modifications are introduced separately for relativistic and non-relativistic particles, by way of the parameters $G_{\rm matter}(k,z)$ and $G_{\rm light}(k,z)$ using two bins in both scale and time, with transition wavenumber $0.01$ Mpc$^{-1}$ and redshift 1. We emphasize the use of two dynamical probes to constrain this set of parameters, galaxy power spectrum multipoles and the direct peculiar velocity power spectrum, which probe fluctuations on different scales. The multipole measurements are derived from the WiggleZ and BOSS Data Release 11 CMASS galaxy redshift surveys and the velocity power spectrum is measured from the velocity sub-sample of the 6-degree Field Galaxy Survey. We combine with additional cosmological probes including baryon acoustic oscillations, Type Ia SNe, the cosmic microwave background (CMB), lensing of the CMB, and the temperature--galaxy cross-correlation. Using a Markov Chain Monte Carlo likelihood analysis, we find the inferred best-fit parameter values of $G_{\rm matter}(k,z)$ and $G_{\rm light}(k,z)$ to be consistent with the standard model at the $95\%$ confidence level. Furthermore, accounting for the Alcock-Paczynski effect, we perform joint fits for the expansion history and growth index gamma; we measure $γ= 0.665 \pm 0.0669$ ($68\%$ C.L) for a fixed expansion history, and $γ= 0.73^{+0.08}_{-0.10}$ ($68\%$ C.L) when the expansion history is allowed to deviate from $Λ$CDM. With a fixed expansion history the inferred value is consistent with GR at the $95\%$ C.L; alternatively, a $2σ$ tension is observed when the expansion history is not fixed, this tension is worsened by the combination of growth and SNe data.

preprint2015arXiv

The 6dF Galaxy Survey: Bulk Flows on $50-70 h^{-1}$ Mpc scales

We measure the bulk flow of the local Universe using the 6dF Galaxy Survey peculiar velocity sample (6dFGSv), the largest and most homogeneous peculiar velocity sample to date. 6dFGSv is a Fundamental Plane sample of $\sim10^4$ peculiar velocities covering the whole southern hemisphere for galactic latitude $|b| > 10^\circ$, out to redshift ${z=0.0537}$. We apply the `Minimum Variance' bulk flow weighting method, which allows us to make a robust measurement of the bulk flow on scales of $50$ and $70\,h^{-1}{\rm Mpc}$. We investigate and correct for potential bias due to the lognormal velocity uncertainties, and verify our method by constructing $Λ{\rm CDM}$ 6dFGSv mock catalogues incorporating the survey selection function. For a hemisphere of radius $50\,h^{-1}{\rm Mpc}$ we find a bulk flow amplitude of $U=248\pm58\,{\rm km}\,{\rm s}^{-1}$ in the direction $(l,b) = (318^\circ\pm20^\circ,40^\circ\pm13^\circ)$, and for $70\,h^{-1}{\rm Mpc}$ we find $U=243\pm58\,{\rm km}\,{\rm s}^{-1}$, in the same direction. Our measurement gives us a constraint on $σ_8$ of $1.01^{+1.07}_{-0.58}$. Our results are in agreement with other recent measurements of the direction of the bulk flow, and our measured amplitude is consistent with a $Λ{\rm CDM}$ prediction.

preprint2015arXiv

The BOSS-WiggleZ overlap region I: Baryon Acoustic Oscillations

We study the large-scale clustering of galaxies in the overlap region of the Baryon Oscillation Spectroscopic Survey (BOSS) CMASS sample and the WiggleZ Dark Energy Survey. We calculate the auto-correlation and cross-correlation functions in the overlap region of the two datasets and detect a Baryon Acoustic Oscillation (BAO) signal in each of them. The BAO measurement from the cross-correlation function represents the first such detection between two different galaxy surveys. After applying density-field reconstruction we report distance-scale measurements $D_V r_s^{\rm fid} / r_s = (1970 \pm 47, 2132 \pm 67, 2100 \pm 200)$ Mpc from CMASS, the cross-correlation and WiggleZ, respectively. We use correlated mock realizations to calculate the covariance between the three BAO constraints. The distance scales derived from the two datasets are consistent, and are also robust against switching the displacement fields used for reconstruction between the two surveys. This approach can be used to construct a correlation matrix, permitting for the first time a rigorous combination of WiggleZ and CMASS BAO measurements. Using a volume-scaling technique, our result can also be used to combine WiggleZ and future CMASS DR12 results. Finally, we use the cross-correlation function measurements to show that the relative velocity effect, a possible source of systematic uncertainty for the BAO technique, is consistent with zero for our samples.

preprint2015arXiv

The BOSS-WiggleZ overlap region II: dependence of cosmic growth on galaxy type

The anisotropic galaxy 2-point correlation function (2PCF) allows measurement of the growth of large-scale structures from the effect of peculiar velocities on the clustering pattern. We present new measurements of the auto- and cross- correlation function multipoles of 69,180 WiggleZ and 46,380 BOSS-CMASS galaxies sharing an overlapping volume of ~0.2 (Gpc/h)^3. Analysing the redshift-space distortions (RSD) of galaxy 2-point statistics for these two galaxy tracers, we test for systematic errors in the modelling depending on galaxy type and investigate potential improvements in cosmological constraints. We build a large number of mock galaxy catalogues to examine the limits of different RSD models in terms of fitting scales and galaxy type, and to study the covariance of the measurements when performing joint fits. For the galaxy data, fitting the monopole and quadrupole of the WiggleZ 2PCF on scales 24<s<80 Mpc/h produces a measurement of the normalised growth rate $fσ_8$(z=0.54)=0.409$\pm$0.059, whereas for the CMASS galaxies we found a consistent constraint of $fσ_8$(z=0.54)=0.466$\pm$0.074. When combining the measurements, accounting for the correlation between the two surveys, we obtain $fσ_8$(z=0.54)=0.413$\pm$0.054, in agreement with the LCDM-GR model of structure growth and with other survey measurements.

preprint2015arXiv

The Impact of Nonlinear Structure Formation on the Power Spectrum of Transverse Momentum Fluctuations and the Kinetic Sunyaev-Zel'dovich Effect

Cosmological transverse momentum fields, whose directions are perpendicular to Fourier wave vectors, induce temperature anisotropies in the cosmic microwave background via the kinetic Sunyaev-Zeldovich (kSZ) effect. The transverse momentum power spectrum contains the four-point function of density and velocity fields, $\langleδδv v\rangle$. In the post-reionization epoch, nonlinear effects dominate in the power spectrum. We use perturbation theory and cosmological $N$-body simulations to calculate this nonlinearity. We derive the next-to-leading order expression for the power spectrum with a particular emphasis on the connected term that has been ignored in the literature. While the contribution from the connected term on small scales ($k>0.1\,h\,\rm{Mpc}^{-1}$) is subdominant relative to the unconnected term, we find that its contribution to the kSZ power spectrum at $\ell = 3000$ at $z<6$ can be as large as ten percent of the unconnected term, which would reduce the allowed contribution from the reionization epoch ($z>6$) by twenty percent. The power spectrum of transverse momentum on large scales is expected to scale as $k^2$ as a consequence of momentum conservation. We show that both the leading and the next-to-leading order terms satisfy this scaling. In particular, we find that both of the unconnected and connected terms are necessary to reproduce $k^2$.

preprint2014arXiv

The 6dF Galaxy Velocity Survey: Cosmological constraints from the velocity power spectrum

We present scale-dependent measurements of the normalised growth rate of structure $fσ_{8}(k, z=0)$ using only the peculiar motions of galaxies. We use data from the 6-degree Field Galaxy Survey velocity sample (6dFGSv) together with a newly-compiled sample of low-redshift $(z < 0.07)$ type Ia supernovae. We constrain the growth rate in a series of $Δk \sim 0.03 h{\rm Mpc^{-1}}$ bins to $\sim35\%$ precision, including a measurement on scales $>300 h^{-1}{\rm Mpc}$, which represents one of the largest-scale growth rate measurement to date. We find no evidence for a scale dependence in the growth rate, or any statistically significant variation from the growth rate as predicted by the {\it Planck} cosmology. Bringing all the scales together, we determine the normalised growth rate at $z=0$ to $\sim15\%$ in a manner {\it independent} of galaxy bias and in excellent agreement with the constraint from the measurements of redshift-space distortions from 6dFGS. We pay particular attention to systematic errors. We point out that the intrinsic scatter present in Fundamental-Plane and Tully-Fisher relations is only Gaussian in logarithmic distance units; wrongly assuming it is Gaussian in linear (velocity) units can bias cosmological constraints. We also analytically marginalise over zero-point errors in distance indicators, validate the accuracy of all our constraints using numerical simulations, and demonstrate how to combine different (correlated) velocity surveys using a matrix `hyper-parameter' analysis. Current and forthcoming peculiar velocity surveys will allow us to understand in detail the growth of structure in the low-redshift universe, providing strong constraints on the nature of dark energy.

preprint2014arXiv

The WiggleZ Dark Energy Survey: Improved Distance Measurements to z = 1 with Reconstruction of the Baryonic Acoustic Feature

We present significant improvements in cosmic distance measurements from the WiggleZ Dark Energy Survey, achieved by applying the reconstruction of the baryonic acoustic feature technique. We show using both data and simulations that the reconstruction technique can often be effective despite patchiness of the survey, significant edge effects and shot-noise. We investigate three redshift bins in the redshift range 0.2<$z$<1, and in all three find improvement after reconstruction in the detection of the baryonic acoustic feature and its usage as a standard ruler. We measure model independent distance measures $D_{\mathrm V}(r_{\mathrm s}^\mathrm{fid}/r_{\mathrm s})$ of 1716 $\pm$ 83 Mpc, 2221 $\pm$ 101 Mpc, 2516 $\pm$ 86 Mpc (68% CL) at effective redshifts z = 0.44, 0.6, 0.73, respectively, where $D_{\mathrm V}$ is the volume-average-distance, and $r_{\mathrm s}$ is the sound horizon at the end of the baryon drag epoch. These significantly improved 4.8, 4.5 and 3.4 percent accuracy measurements are equivalent to those expected from surveys with up to 2.5 times the volume of WiggleZ. These measurements are fully consistent with cosmologies allowed by the analyses of the Planck Collaboration and the Sloan Digital Sky Survey.We provide the $D_{\mathrm V}(r_{\mathrm s}^\mathrm{fid}/r_{\mathrm s})$ posterior probability distributions and their covariances. When combining these measurements with temperature fluctuations measurements of Planck, the polarization of WMAP9, and the 6dF Galaxy Survey baryonic acoustic feature, we do not detect deviations from a flat LCDM model. Assuming this model we constrain the current expansion rate to $H_0$ = 67.15 $\pm$ 0.98 kms$^{-1}$Mpc$^{-1}$. Allowing the equation of state of dark energy to vary we obtain $w_\mathrm{DE}$ = -1.080 $\pm$ 0.135. When assuming a curved LCDM model we obtain a curvature value of $Ω_{\mathrm K}$ = -0.0043 $\pm$ 0.0047.

preprint2013arXiv

Are peculiar velocity surveys competitive as a cosmological probe?

Peculiar velocity surveys, which measure galaxy peculiar velocities directly from standard candles in addition to redshifts, can provide strong constraints on the linear growth rate of cosmological large-scale structure at low redshift. The improvement originates from the physical relationship between galaxy density and peculiar velocity, which substantially reduces cosmic variance. We present the results of Fisher matrix forecasts of correlated fields of galaxy density and velocity. Peculiar velocity can improve the growth rate constraints by about a factor of two compared to density alone, if we can use all the information for wavenumber k < 0.2 h/Mpc. Future peculiar velocity surveys, TAIPAN, and the all-sky HI surveys, WALLABY and WNSHS, can measure the growth rate, f*sigma8 to 3 per cent at z ~ 0.025. Although the velocity subsample is about an order of magnitude smaller than the redshift sample from the same survey, it improves the constraint by 40 per cent compared to the same survey without velocity measurements. Peculiar velocity surveys can also measure the growth rate as a function of wavenumber with 15-30 per cent uncertainties in bins with widths 0.01 h/Mpc in the range 0.02 h/Mpc < k < 0.1 h/Mpc, which is a large improvement over galaxy density only. Such measurements on very large scales can detect signatures of modified gravity or non-Gaussianity through scale-dependent growth rate or galaxy bias. We use N-body simulations to improve the modelling of auto- and cross-power spectra of galaxy density and peculiar velocity by introducing a new redshift-space distortion term to the velocity, which has been neglected in previous studies. The velocity power spectrum has a damping in redshift space, which is larger than that naively expected from the similar effect in the galaxy power spectrum.

preprint2013arXiv

Will Nonlinear Peculiar Velocity and Inhomogeneous Reionization Spoil 21cm Cosmology from the Epoch of Reionization?

The 21cm background from the epoch of reionization is a promising cosmological probe: line-of-sight velocity fluctuations distort redshift, so brightness fluctuations in Fourier space depend upon angle, which linear theory shows can separate cosmological from astrophysical information. Nonlinear fluctuations in ionization, density and velocity change this, however. The validity and accuracy of the separation scheme are tested here for the first time, by detailed reionization simulations. The scheme works reasonably well early in reionization (< 40% ionized), but not late (> 80% ionized).

preprint2012arXiv

Detecting the Rise and Fall of the First Stars by Their Impact on Cosmic Reionization

The intergalactic medium was reionized before redshift z~6, most likely by starlight which escaped from early galaxies. The very first stars formed when hydrogen molecules (H2) cooled gas inside the smallest galaxies, minihalos of mass between 10^5 and 10^8 solar masses. Although the very first stars began forming inside these minihalos before redshift z~40, their contribution has, to date, been ignored in large-scale simulations of this cosmic reionization. Here we report results from the first reionization simulations to include these first stars and the radiative feedback that limited their formation, in a volume large enough to follow the crucial spatial variations that influenced the process and its observability. We show that, while minihalo stars stopped far short of fully ionizing the universe, reionization began much earlier with minihalo sources than without, and was greatly extended, which boosts the intergalactic electron-scattering optical depth and the large-angle polarization fluctuations of the cosmic microwave background significantly. Although within current WMAP uncertainties, this boost should be readily detectable by Planck. If reionization ended as late as z_ov<~7, as suggested by other observations, Planck will thereby see the signature of the first stars at high redshift, currently undetectable by other probes.

preprint2012arXiv

Redshift Space Distortion of the 21cm Background from the Epoch of Reionization I: Methodology Re-examined

The peculiar velocity of the intergalactic gas responsible for the cosmic 21cm background from the epoch of reionization and beyond introduces an anisotropy in the three-dimensional power spectrum of brightness temperature fluctuations. Measurement of this anisotropy by future 21cm surveys is a promising tool for separating cosmology from 21cm astrophysics. However, previous attempts to model the signal have often neglected peculiar velocity or only approximated it crudely. This paper re-examines the effects of peculiar velocity on the 21cm signal in detail, improving upon past treatment and addressing several issues for the first time. (1) We show that properly accounting for finite optical depth eliminates the unphysical divergence of 21cm brightness temperature in overdense regions of the IGM found by previous work that employed the usual optically-thin approximation. (2) The approximation made previously to circumvent the diverging brightness temperature problem by capping velocity gradient can misestimate the power spectrum on all scales. (3) The observed power spectrum in redshift-space remains finite even in the optically-thin approximation if one properly accounts for the redshift-space distortion. However, results that take full account of finite optical depth show that this approximation is only accurate in the limit of high spin temperature. (4) The linear theory for redshift-space distortion results in ~30% error in the observationally relevant wavenumber range, at the 50% ionized epoch. (5) We describe and test two numerical schemes to calculate the 21cm signal from reionization simulations to incorporate peculiar velocity effects in the optically-thin approximation accurately. One is particle-based, the other grid-based, and while the former is most accurate, we demonstrate that the latter is computationally more efficient and can achieve sufficient accuracy. [Abridged]

preprint2012arXiv

Simulating Cosmic Reionization and the Radiation Backgrounds from the Epoch of Reionization

Large-scale reionization simulations are described which combine the results of cosmological N-body simulations that model the evolving density and velocity fields and identify the galactic halo sources, with ray-tracing radiative transfer calculations which model the nonequilibrium ionization of the intergalactic medium. These simulations have been used to predict some of the signature effects of reionization on cosmic radiation backgrounds, including the CMB, near-IR, and redshifted 21cm backgrounds. We summarize some of our recent progress in this work, and address the question of whether observations of such signature effects can be used to distinguish the relative contributions of galaxies of different masses to reionization.

preprint2011arXiv

Can 21-cm observations discriminate between high-mass and low-mass galaxies as reionization sources?

The prospect of detecting the first galaxies by observing their impact on the intergalactic medium as they reionized it during the first billion years leads us to ask whether such indirect observations are capable of diagnosing which types of galaxies were most responsible for reionization. We attempt to answer this by considering a set of large-scale radiative transfer simulations of reionization in sufficiently large volumes to make statistically meaningful predictions of observable signatures, while also directly resolving all atomically-cooling halos down to 10^8 M_solar. We focus here on predictions of the 21-cm background, to see if upcoming observations are capable of distinguishing a universe ionized primarily by high-mass halos from one in which both high-mass and low-mass halos are responsible, and to see how these results depend upon the uncertain source efficiencies. We find that 21-cm fluctuation power spectra observed by the first generation EoR/21-cm radio interferometer arrays should be able to distinguish the case of reionization by high-mass halos alone from that by both high- and low-mass halos, together. Some reionization scenarios yield very similar power spectra and rms evolution and thus can only be discriminated by their different mean reionization history and 21-cm PDF distributions. We find that the skewness of the 21-cm PDF distribution smoothed over LOFAR-like window shows a clear feature correlated with the rise of the rms due to patchiness. Measurements of the mean photoionization rates are sensitive to the average density of the regions being studied and therefore could be strongly skewed in certain cases. (abridged)

preprint2011arXiv

Gravothermal collapse of isolated self-interacting dark matter haloes: N-body simulation versus the fluid model

Self-Interacting Dark Matter (SIDM) is a collisional form of cold dark matter (CDM), originally proposed to solve problems that arose when the collisionless CDM theory of structure formation was compared with observations of galaxies on small scales. The quantitative impact of the proposed elastic collisions on structure formation has been estimated previously by Monte Carlo N-body simulations and by a conducting fluid model, with apparently diverging results. To improve this situation, we make direct comparisons between new Monte Carlo N-body simulations and solutions of the conducting fluid model, for isolated SIDM haloes of fixed mass. This allows us to separate cleanly the effects of gravothermal relaxation from those of continuous mass accretion in an expanding background universe. When these two methods were previously applied to halo formation with cosmological boundary conditions, they disagreed by an order of magnitude about the size of the scattering cross section required to solve the so-called 'cusp-core problem.' We show here, however, that the methods agree with each other within 20 per cent for isolated haloes. This suggests that the two methods are consistent, and that their disagreement for cosmological haloes is not caused by a breakdown of their validity. The isolated haloes studied here undergo gravothermal collapse. We compare the solutions calculated by these two methods for gravothermal collapse starting from several initial conditions. This allows us to calibrate the heat conduction which accounts for the effect of elastic hard-sphere scattering in the fluid model. The amount of tuning of the thermal conductivity parameters required to bring the two methods into close agreement for isolated haloes, however, is too small to explain the discrepancy found previously in the cosmological context.

preprint2010arXiv

Cosmic Structure Formation at High Redshift

We present some preliminary results from a series of extremely large, high-resolution N-body simulations of the formation of early nonlinear structures. We find that the high-z halo mass function is inconsistent with the Sheth-Tormen mass function, which tends to over-estimate the abundance of rare halos. This discrepancy is in rough agreement with previous results based on smaller simulations. We also show that the number density of minihaloes is correlated with local matter density, albeit with a significant scatter that increases with redshift, as minihaloes become increasingly rare. The average correlation is in rough agreement with a simple analytical extended Press-Schechter model, but can differ by up to factor of 2 in some regimes.