Search

Intensity interferometry is a technique developed many decades ago, that has recently enjoyed a renaissance thanks in part to advances in photodetector technology. We investigate the potential for long-baseline optical intensity interferometry to observe bright, active galactic nuclei (AGN) associated with rapidly accreting supermassive black holes. We argue that realistic telescope arrays similar in area to existing Cherenkov arrays, if equipped with modern high-precision single photon detectors, can achieve a sufficiently high signal to noise ratio not only to detect distant AGN, but also to study them in great detail. We explore the science potential of such observations by considering two examples. First, we find that intensity interferometric observations of bright nearby AGN can allow detailed studies of the central accretion disks powering the AGN, allowing reconstruction of many disk properties like the radial profile. Next, we argue that intensity interferometers can spatially resolve the broad-line regions of AGN at cosmological distances, and thereby provide a geometric determination of the angular diameter distances to those AGN when combined with reverberation mapping. Since this measurement can be performed for AGN at distances of hundreds of megaparsecs, this directly measures the Hubble expansion rate $H_0$, with a precision adequate to resolve the recent Hubble tension. Finally, we speculate on future applications that would be enabled by even larger intensity interferometer arrays.

The velocity profile of galaxies around voids is a key ingredient for redshift space distortion (RSD) measurements made using the void-galaxy correlation function. In this paper we use simulations to test whether the velocity profile of the tracers used to find the voids matches the velocity profile of the dark matter around these voids. A mismatch is expected and found in the inner part of voids, where tracers are very sparse. We discuss how this difference is caused by a selection effect where the void centre positions are correlated to the particular realization of the sparse tracers and their spatial distribution. In turn, this then affects the RSD void-galaxy correlation analysis. We show this by evaluating the Jacobian of the real to redshift space mapping using the tracer or matter velocity profile. Differences of the order of 20\% in the velocity profile translate into differences of the order of few percent in the Jacobian. This small discrepancy propagates to the monopole and quadrupole of the void-tracer correlation function, producing modifications of comparable magnitude to those from changes in $f\sigma_8$ at the level of the statistical uncertainties from current analyses., Comment: Author's Accepted Manuscript

We use observations of ultra-faint dwarf (UFD) galaxies to constrain the particle mass of ultra-light dark matter. Potential fluctuations created by wave interference in virialized "fuzzy" dark matter (FDM) halos dynamically heat stellar orbits in UFDs, some of which exhibit velocity dispersions of $\lesssim$ 3 km/s and sizes $\lesssim$ 40 pc. Using simulations of FDM halos, and existing measurements of sizes and stellar radial velocities in Segue 1 and Segue 2 UFDs, we derive a lower limit on the dark matter particle mass of $m_{fdm} > 3\times 10^{-19}$ eV at 99% confidence, marginalized over host halo circular velocity. This constraint is conservative as it is derived under the assumption that soliton heating is negligible, and that no other sources of non-FDM dynamical heating of stars operate to increase velocity dispersion. It can potentially be strengthened by future spectroscopic observations of additional stars in ultra-faint galaxies and by tightening theoretical constraints on the soliton size-halo mass relation. However, even the current conservative lower limit on the FDM mass makes this model indistinguishable from Cold Dark Matter at the scales probed by existing astronomical observations., Comment: 10 pages. Added references

Secondary halo properties beyond mass, such as the mass accretion rate (MAR), concentration, and the half mass scale, are essential in understanding the formation of large-scale structure and dark matter halos. In this paper, we study the impact of secondary halo properties on the galaxy-galaxy lensing observable, $\Delta\Sigma$. We build an emulator trained on N-body simulations to model $\Delta\Sigma$ and quantify the impact of different secondary parameters on the $\Delta\Sigma$ profile. We focus on the impact of MAR on $\Delta\Sigma$. We show that a 3$\sigma$ detection of variations in MAR at fixed halo mass could be achieved with the Hyper Suprime Cam survey in combination with a proxy for MAR with scatter $\sigma_{\Gamma_\mathrm{dyn}|\mathrm{obs}}<1.5$. We show that the full radial profile of $\Delta\Sigma$ depends on secondary properties at fixed halo mass. Consequently, an emulator that can perform full shape fitting yields better than 2 times improvement upon the constraints on MAR than only using the outer part of the halo. Finally, we highlight that miscentering and MAR impact the radial profile of $\Delta\Sigma$ in a similar fashion, implying that miscentering and MAR need to be modeled jointly for unbiased estimates of both effects. We show that present-day lensing data sets have the statistical capability to place constraints on halo MAR. Our analysis opens up new possibilities for observationally measuring the assembly history of the dark matter halos that host galaxies and clusters.

We study how tidal streams from globular clusters may be used to constrain the mass of ultra-light dark matter particles, called `fuzzy' dark matter (FDM). A general feature of FDM models is the presence of ubiquitous density fluctuations in bound, virialized dark matter structures, on the scale of the de Broglie wavelength, arising from wave interference in the evolving dark matter distribution. These time-varying fluctuations can disturb the motions of stars, leading to potentially observable signatures in cold thin tidal streams in our own Galaxy. The study of this effect has been hindered by the difficulty in simulating the FDM wavefunction in Milky Way-sized systems. We present a simple method to evolve realistic wavefunctions in nearly static potentials, that should provide an accurate estimate of this granulation effect. We quantify the impact of FDM perturbations on tidal streams, and show that initially, while stream perturbations are small in amplitude, their power spectra exhibit a sharp cutoff corresponding to the de Broglie wavelength of the FDM potential fluctuations. Eventually, when stream perturbations become nonlinear, fold caustics generically arise that lead to density fluctuations with universal behavior. This erases the signature of the de Broglie wavelength in the stream density power spectrum, but we show that it will still be possible to determine the FDM mass in this regime, by considering the fluctuations in quantities like angular momenta or actions., Comment: Comments welcome. To be submitted to JCAP

We present evidence of anisotropic galaxy assembly bias in the Baryon Oscillation Spectroscopic Survey Data Release 12 galaxy sample at a level exceeding $5\sigma$. We use measurements of the line-of-sight velocity dispersion $\sigma_\star$ and stellar mass $M_\star$ to perform a simple split into subsamples of galaxies. We show that the amplitude of the monopole and quadrupole moments of the power spectrum depend differently on $\sigma_\star$ and $M_\star$, allowing us to split the galaxy sample into subsets with matching monopoles but significantly different quadrupoles on all scales. Combining data from the LOWZ and CMASS NGC galaxy samples, we find $>5\sigma$ evidence for anisotropic bias on scales $k<0.15\,h\,\rm{Mpc}^{-1}$. We also examine splits using other observed properties. For galaxy samples split using $M_\star$ and projected size $R_0$, we find no significant evidence of anisotropic bias. Galaxy samples selected using additional properties exhibit strongly varying degrees of anisotropic assembly bias, depending on which combination of properties is used to split into subsets. This may explain why previous searches for this effect using the Fundamental Plane found inconsistent results. We conclude that any selection of a galaxy sample that depends on $\sigma_\star$ can give biased and incorrect Redshift Space Distortion measurements., Comment: 15 pages, 11 figures. Comments welcome

Cosmological neutrinos have their greatest influence in voids: these are the regions with the highest neutrino to dark matter density ratios. The marked power spectrum can be used to emphasize low density regions over high density regions, and therefore is potentially much more sensitive than the power spectrum to the effects of neutrino masses. Using 22,000 N-body simulations from the Quijote suite, we quantify the information content in the marked power spectrum of the matter field, and show that it outperforms the standard power spectrum by setting constraints improved by a factor larger than 2 on all cosmological parameters. The combination of marked and standard power spectrum allows to place a 4.3{\sigma} constraint on the minimum sum of the neutrino masses with a volume equal to 1 (Gpc/h)^3 and without CMB priors. Combinations of different marked power spectra yield a 6{\sigma} constraint within the same conditions., Comment: 5 pages, 3 figures

We study the effect of large-scale tidal fields on internal halo properties using a set of N-body simulations. We measure significant cross-correlations between large-scale tidal fields and several non-scalar halo properties: shapes, velocity dispersion, and angular momentum. Selection effects that couple to these non-scalar halo properties can produce anisotropic clustering even in real-space. We investigate the size of this effect and show that it can produce a non-zero quadrupole similar in size to the one generated by linear redshift-space distortions (RSD). Finally, we investigate the clustering properties of halos identified in redshift-space and find enormous deviations from the standard linear RSD model, again caused by anisotropic assembly bias. These effects could contaminate the values of cosmological parameters inferred from the observed redshift-space clustering of galaxies, groups, or 21cm emission from atomic hydrogen, if their selection depends on properties affected by halo assembly bias. We briefly discuss ways in which this effect can be measured in existing and future large-scale structure surveys., Comment: 9 pages, 4 figures. Comments welcome

Non-gravitational interactions between dark matter particles with strong scattering, but relatively small annihilation and dissipation, has been proposed to match various observables on cluster and group scales. In this paper, we present the results from large cosmological simulations which include the effects of different self-interaction scenarios. In particular we explore a model with the differential cross section that can depend on both the relative velocity of the interacting particles and the angle of scattering. We focus on how quantities, such as the stacked density profiles, subhalo counts and the splashback radius change as a function of different forms of self-interaction. We find that self-interactions not only affect the central region of the cluster, the effect well known from previous studies, but also significantly alter the distribution of subhalos and the density of particles out to the splashback radius. Our results suggest that current weak lensing data can already put constraints on the self-interaction cross-section that are only slightly weaker than the Bullet Cluster constraints ($\sigma/m \lesssim 2$ cm$^2/$g), and future lensing surveys should be able to tighten them even further making halo profiles on cluster scales a competitive probe for DM physics., Comment: To be submitted to JCAP

Gravitational potentials which change in time induce fluctuations in the observed cosmic microwave background (CMB) temperature. Cosmological structure moving transverse to our line of sight provides a specific example known as the moving lens effect. Here we explore how the observed CMB temperature fluctuations combined with the observed matter over-density can be used to infer the transverse velocity of cosmological structure on large scales. We show that near-future CMB surveys and galaxy surveys will have the statistical power to make a first detection of the moving lens effect, and we discuss applications for the reconstructed transverse velocity., Comment: 5 pages, 4 figures, comments welcome

We study the kinematics of galaxies within massive clusters, as a probe of the physics of star-formation quenching within clusters. Using N-body simulations, we argue that satellite kinematics provide information about galaxy infall that is complementary to the (instantaneous) spatial distribution of satellites. Comparing the simulation results with measurements of real cluster galaxies, we find evidence that the kinematics of red (quiescent) satellite galaxies are consistent with earlier infall times than that of blue (star-forming) satellites., Comment: 15 pages. To be submitted to ApJ

The splashback radius is a physical scale in dark matter halos that is set by the gravitational dynamics of recently accreted shells. We use analytical models and N-body simulations to study the dependence of splashback on dark energy and screened modified gravity theories. In modified gravity models, the transition from screened to unscreened regions typically occurs in the cluster outskirts, suggesting potentially observable signatures in the splashback feature. We investigate the location of splashback in both chameleon and Vainshtein screened models and find significant differences compared with $\Lambda$CDM predictions. We also find an interesting interplay between dynamical friction and modified gravity, providing a distinctive signature for modified gravity models in the behavior of the splashback feature as a function of galaxy luminosity., Comment: To be submitted to JCAP, 12 pages

We study halo assembly bias for cluster-sized halos. Previous work has found little evidence for correlations between large-scale bias and halo mass assembly history for simulated cluster-sized halos, in contrast to the significant correlation found between bias and concentration for halos of this mass. This difference in behavior is surprising, given that both concentration and assembly history are closely related to the same properties of the linear-density peaks that collapse to form halos. Using publicly available simulations, we show that significant assembly bias is indeed found in the most massive halos with $M\sim 10^{15}M_\odot$, using essentially any definition of halo age. For lower halo masses $M\sim 10^{14}M_\odot$, no correlation is found between bias and the commonly used age indicator $a_{0.5}$, the half-mass time. We show that this is a mere accident, and that significant assembly bias exists for other definitions of halo age, including those based on the time when the halo progenitor acquires some fraction $f$ of the ultimate mass at $z=0$. For halos with $M_{\rm vir}\sim 10^{14}M_\odot$, the sense of assembly bias changes sign at $f=0.5$. We explore the origin of this behavior, and argue that it arises because standard definitions of halo mass in halo finders do not correspond to the collapsed, virialized mass that appears in the spherical collapse model used to predict large-scale clustering. Because bias depends strongly on halo mass, these errors in mass definition can masquerade as or even obscure the assembly bias that is physically present. More physically motivated halo definitions using splashback should be free of this particular defect of standard halo finders., Comment: 7 pages, 4 figures, to be submitted to JCAP

Cosmology is intrinsically intertwined with questions in fundamental physics. The existence of non-baryonic dark matter requires new physics beyond the Standard Model of elemenatary-particle interactions and Einstein's general relativity, as does the accelerating expansion of the universe. Current tensions between various cosmological measurements may be harbingers of yet more new physics. Progress on understanding dark matter and cosmic acceleration requires long term, high-precision measurements and excellent control of systematics, demanding observational programs that are often outside the discovery/characterization mode that drives many areas of astronomy. We outline potential programs through which the Hubble Space Telescope (HST) could have a major impact on issues in fundamental physics in the coming years. To realize this impact, we suggest the introduction of a "HST Fundamental Physics" observational program that would be subject to a modified proposal and review process., Comment: Report of the HST and Fundamental Physics Working Group, constituted by STScI Director Ken Sembach to identify ways that HST could enable advances in fundamental physics

Self-Interacting Dark Matter (SIDM) has long been proposed as a solution to small scale problems posed by standard Cold Dark Matter (CDM). We use numerical simulations to study the effect of dark matter interactions on the morphology of disk galaxies falling into galaxy clusters. The effective drag force on dark matter leads to offsets of the stellar disk with respect to the surrounding halo, causing distortions in the disk. For anisotropic scattering cross-sections of 0.5 and 1.0$\,\textrm{cm}^{2}\textrm{g}^{-1}$, we show that potentially observable warps, asymmetries, and thickening of the disk occur in simulations. We discuss observational tests of SIDM with galaxy surveys and more realistic simulations needed to obtain detailed predictions., Comment: Title changed and minor modifications to text; matches version published in ApJ

(abridged) We investigate the signatures left by the cosmic neutrino background on the clustering of matter, CDM+baryons and halos in redshift-space using a set of more than 1000 N-body and hydrodynamical simulations with massless and massive neutrinos. We find that the effect neutrinos induce on the clustering of CDM+baryons in redshift-space on small scales is almost entirely due to the change in $\sigma_8$. Neutrinos imprint a characteristic signature in the quadrupole of the matter (CDM+baryons+neutrinos) field on small scales, that can be used to disentangle the effect of $\sigma_8$ and $M_\nu$. We show that the effect of neutrinos on the clustering of halos is very different, on all scales, to the one induced by $\sigma_8$. We find that the effects of neutrinos of the growth rate of CDM+baryons ranges from $\sim0.3\%$ to $2\%$ on scales $k\in[0.01, 0.5]~h{\rm Mpc}^{-1}$ for neutrinos with masses $M_\nu \leqslant 0.15$ eV. We compute the bias between the momentum of halos and the momentum of CDM+baryon and find it to be 1 on large scales for all models with massless and massive neutrinos considered. This point towards a velocity bias between halos and total matter on large scales that it is important to account for in order to extract unbiased neutrino information from velocity/momentum surveys such as kSZ observations. We show that baryonic effects can affect the clustering of matter and CDM+baryons in redshift-space by up to a few percent down to $k=0.5~h{\rm Mpc}^{-1}$. We find that hydrodynamics and astrophysical processes, as implemented in our simulations, only distort the relative effect that neutrinos induce on the anisotropic clustering of matter, CDM+baryons and halos in redshift-space by less than $1\%$. Thus, the effect of neutrinos in the fully non-linear regime can be written as a transfer function with very weak dependence on astrophysics., Comment: 19 pages, 12 figures

Mass around dark matter halos can be divided into "infalling" material and "collapsed" material that has passed through at least one pericenter. Analytical models and simulations predict a rapid drop in the halo density profile associated with the transition between these two regimes. Using data from SDSS, we explore the evidence for such a feature in the density profiles of galaxy clusters and investigate the connection between this feature and a possible phase space boundary. We first estimate the steepening of the outer galaxy density profile around clusters: the profiles show an abrupt steepening, providing evidence for truncation of the halo profile. Next, we measure the galaxy density profile around clusters using two sets of galaxies selected based on color. We find evidence of an abrupt change in the galaxy colors that coincides with the location of the steepening of the density profile. Since galaxies are likely to be quenched of star formation and turn red inside of clusters, this change in the galaxy color distribution can be interpreted as the transition from an infalling regime to a collapsed regime. We also measure this transition using a model comparison approach which has been used recently in studies of the "splashback" phenomenon, but find that this approach is not a robust way to quantify the significance of detecting a splashback-like feature. Finally, we perform measurements using an independent cluster catalog to test for potential systematic errors associated with cluster selection. We identify several avenues for future work: improved understanding of the small-scale galaxy profile, lensing measurements, identification of proxies for the halo accretion rate, and other tests. With upcoming data from the DES, KiDS and HSC surveys, we can expect significant improvements in the study of halo boundaries., Comment: 17 pages, 8 figures

We analyze the ability of galaxy and CMB lensing surveys to constrain massive neutrinos and new models of dark radiation. We present a Fisher forecast analysis for neutrino mass constraints with the LSST galaxy survey and the CMB S4 survey. A joint analysis of the three galaxy and shear 2-point functions, along with key systematics parameters and Planck priors, constrains the neutrino masses to $\sum m_\nu = 0.041\,$eV at 1-$\sigma$ level, comparable to constraints expected from Stage 4 CMB lensing. If low redshift information from upcoming spectroscopic surveys like DESI is included, the constraint becomes $\sum m_\nu = 0.032\,$eV. These constraints are derived having marginalized over the number of relativistic species ($N_{\rm eff}$), which is somewhat degenerate with the neutrino mass. We also explore the gain by combining LSST and CMB S4, that is, using the five relevant auto- and cross-correlations of the two datasets. We conclude that advances in modeling the nonlinear regime and the measurements of other parameters are required to ensure a neutrino mass detection. Using the same datasets, we explore the ability of LSST-era surveys to test "nonstandard" models with dark radiation. We find that if evidence for dark radiation is found from $N_{\rm eff}$ measurements, the mass of the dark radiation candidate can be measured at a 1-$\sigma$ level of $0.162\,$eV for fermionic dark radiation, and $0.137\,$eV for bosonic dark radiation, for $\Delta N_{\rm eff} = 0.15$. We also find that the NNaturalness model of Arkani-Hamed et al 2016, with extra light degrees of freedom, has a sub-percent effect on the power spectrum: even more ambitious surveys than the ones considered here will be needed to test such models., Comment: References added, minor corrections to abstract

We present a new method for simulating cosmologies that contain massive particles with thermal free streaming motion, such as massive neutrinos or warm/hot dark matter. This method combines particle and fluid descriptions of the thermal species to eliminate the shot noise known to plague conventional N-body simulations. We describe this method in detail, along with results for a number of test cases to validate our method, and check its range of applicability. Using this method, we demonstrate that massive neutrinos can produce a significant scale-dependence in the large-scale biasing of deep voids in the matter field. We show that this scale-dependence may be quantitatively understood using an extremely simple spherical expansion model which reproduces the behavior of the void bias for different neutrino parameters., Comment: 20 pages, to be submitted to JCAP

We present a novel method to detect the effects of dynamical friction in observed galaxy clusters. Following accretion into clusters, massive satellite galaxies will backsplash to systematically smaller radii than less massive satellites, an effect that may be detected by stacking the number density profiles of galaxies around clusters. We show that this effect may be understood using a simple toy model which reproduces the trends with halo properties observed in simulations. We search for this effect using SDSS redMaPPer clusters with richness 10-20) satellites at 99% confidence., Comment: 6 pages, to be submitted to JCAP

When dark matter halos are accreted by massive host clusters, strong gravitational tidal forces begin stripping mass from the accreted subhalos. This stripping eventually removes all mass beyond a subhalo's tidal radius, but the unbound mass remains in the vicinity of the satellite for at least a dynamical time t_dyn. The N-body subhalo study of Chamberlain et al. verified this picture and pointed out a useful observational consequence: measurements of subhalo correlations beyond the tidal radius are sensitive to the infall time, t_infall, of the subhalo onto its host. We perform this cross-correlation measurement using ~ 160,000 red satellite galaxies in SDSS redMaPPer clusters and find evidence that subhalo correlations do persist well beyond the tidal radius, suggesting that many of the observed satellites fell into their current host less than a dynamical time ago, t_infall < t_dyn. Combined with estimated dynamical times t_dyn ~ 3-5 Gyr and SED fitting results for the time at which satellites stopped forming stars, t_quench ~ 6 Gyr, we infer that for a significant fraction of the satellites, star formation quenched before those satellites entered their current hosts. The result holds for red satellites over a large range of cluster-centric distances 0.1 - 0.6 Mpc/h. We discuss the implications of this result for models of galaxy formation., Comment: 9 pages, 9 figures

We show that the projected number density profiles of SDSS photometric galaxies around galaxy clusters displays strong evidence for the splashback radius, a sharp halo edge corresponding to the location of the first orbital apocenter of satellite galaxies after their infall. We split the clusters into two subsamples with different mean projected radial distances of their members, $\langle R_{\rm mem}\rangle$, at fixed richness and redshift, and show that the sample with smaller $\langle R_{\rm mem}\rangle$ has a smaller ratio of the splashback radius to the traditional halo boundary $R_{\rm 200m}$, than the subsample with larger $\langle R_{\rm mem}\rangle$, indicative of different mass accretion rates for the two subsamples. The same cluster samples were recently used by Miyatake et al. to show that their large-scale clustering differs despite their similar weak lensing masses, demonstrating strong evidence for halo assembly bias. We expand on this result by presenting a 6.6-$\sigma$ detection of halo assembly bias using the cluster-photometric galaxy cross-correlations. Our measured splashback radii are smaller, while the strength of the assembly bias signal is stronger, than expectations from N-body simulations based on the $\Lambda$-dominated, cold dark matter structure formation model. Dynamical friction or cluster-finding systematics such as miscentering or projection effects are not likely to be the sole source of these discrepancies., Comment: 18 pages, 13 figures. Have comments, criticisms, new ideas or new tests, please email me. For Savitzky-Golay filter for data with error covariance, see http://ascl.net/code/v/1299

We study the abundance of substructure in the matter density near galaxies using ALMA Science Verification observations of the strong lensing system SDP.81. We present a method to measure the abundance of subhalos around galaxies using interferometric observations of gravitational lenses. Using simulated ALMA observations, we explore the effects of various systematics, including antenna phase errors and source priors, and show how such errors may be measured or marginalized. We apply our formalism to ALMA observations of SDP.81. We find evidence for the presence of a $M=10^{8.96\pm 0.12} M_{\odot}$ subhalo near one of the images, with a significance of $6.9\sigma$ in a joint fit to data from bands 6 and 7; the effect of the subhalo is also detected in both bands individually. We also derive constraints on the abundance of dark matter subhalos down to $M\sim 2\times 10^7 M_{\odot}$, pushing down to the mass regime of the smallest detected satellites in the Local Group, where there are significant discrepancies between the observed population of luminous galaxies and predicted dark matter subhalos. We find hints of additional substructure, warranting further study using the full SDP.81 dataset (including, for example, the spectroscopic imaging of the lensed carbon monoxide emission). We compare the results of this search to the predictions of $\Lambda$CDM halos, and find that given current uncertainties in the host halo properties of SDP.81, our measurements of substructure are consistent with theoretical expectations. Observations of larger samples of gravitational lenses with ALMA should be able to improve the constraints on the abundance of galactic substructure., Comment: 18 pages, 13 figures, Comments are welcome

The fact that the clustering of dark matter halos depends not only on their mass, but also the formation epoch, is a prominent, albeit subtle, feature of the cold dark matter structure formation theory, and is known as assembly bias. At low mass scales ($\sim 10^{12}\,h^{-1}M_\odot$), early-forming halos are predicted to be more strongly clustered than the late-forming ones. In this study we aim to robustly detect the signature of assembly bias observationally, making use of formation time indicators of central galaxies in low mass halos as a proxy for the halo formation history. Weak gravitational lensing is employed to ensure our early- and late-forming halo samples have similar masses, and are free of contamination of satellites from more massive halos. For the two formation time indicators used (resolved star formation history and current specific star formation rate), we do not find convincing evidence of assembly bias. For a pair of early- and late-forming galaxy samples with mean mass $M_{200c} \approx 9\times10^{11}\,h^{-1}M_\odot$, the relative bias is $1.00\pm 0.12$. We attribute the lack of detection to the possibilities that either the current measurements of these indicators are too noisy, or they do not correlate well with the halo formation history. Alternative proxies for the halo formation history that should perform better are suggested for future studies., Comment: ApJ, in press. 14 pages, 6 figures

The statistics of primordial curvature fluctuations are our window into the period of inflation, where these fluctuations were generated. To date, the cosmic microwave background has been the dominant source of information about these perturbations. Large scale structure is however from where drastic improvements should originate. In this paper, we explain the theoretical motivations for pursuing such measurements and the challenges that lie ahead. In particular, we discuss and identify theoretical targets regarding the measurement of primordial non-Gaussianity. We argue that when quantified in terms of the local (equilateral) template amplitude $f_{\rm NL}^{\rm loc}$ ($f_{\rm NL}^{\rm eq}$), natural target levels of sensitivity are $\Delta f_{\rm NL}^{\rm loc, eq.} \simeq 1$. We highlight that such levels are within reach of future surveys by measuring 2-, 3- and 4-point statistics of the galaxy spatial distribution. This paper summarizes a workshop held at CITA (University of Toronto) on October 23-24, 2014., Comment: 27 pages + references

Recent work has shown that density profiles in the outskirts of dark matter halos can become extremely steep over a narrow range of radius. This behavior is produced by splashback material on its first apocentric passage after accretion. We show that the location of this splashback feature may be understood quite simply, from first principles. We present a simple model, based on spherical collapse, that accurately predicts the location of splashback without any free parameters. The important quantities that determine the splashback radius are accretion rate and redshift., Comment: 8 pages

We investigate optimization strategies to measure primordial non-Gaussianity with future spectroscopic surveys. We forecast measurements coming from the 3D galaxy power spectrum and compute constraints on primordial non-Gaussianity parameters f_NL and n_NG. After studying the dependence on those parameters upon survey specifications such as redshift range, area, number density, we assume a reference mock survey and investigate the trade-off between number density and area surveyed. We then define the observational requirements to reach the detection of f_NL of order 1. Our results show that while power spectrum constraints on non-Gaussianity from future spectroscopic surveys can be competitive with current CMB limits, measurements from higher-order statistics will be needed to reach a sub unity precision in the measurements of the non-Gaussianity parameter f_NL., Comment: 12 pages, 10 figures

Satellites within simulated massive clusters are significantly spatially correlated with each other, even when those satellites are not gravitationally bound to each other. This correlation is produced by satellites that entered their hosts relatively recently, and is undetectable for satellites that have resided in their hosts for multiple dynamical timescales. Therefore, a measurement of clustering statistics of cluster satellites may be used to determine the typical accretion redshifts of those satellites into their observed hosts. We argue that such measurements may be used to determine the fraction of satellite galaxies that were quenched by their current hosts, thereby discriminating among models for quenching of star formation in satellite galaxies., Comment: 7 pages

We propose a method to measure the ellipticities of dark matter halos using the lens-shear-shear 3-point correlation function. This method is immune to effects of galaxy-halo misalignments that can potentially limit 2-point galaxy-galaxy lensing measurements of halo anisotropy. Using a simple model for the projected mass distributions of dark matter halos, we construct an ellipticity estimator that sums over all possible triangular configurations of the 3-point function. By applying our estimator to halos from N-body simulations, we find that systematic errors in the recovered ellipticity will be at the <5% fractional level. We estimate that future imaging surveys like LSST will have sufficient statistics to detect halo ellipticities using 3-point lensing., Comment: 9 pages, to be submitted to JCAP

In recent years, it has become possible to detect individual dark matter subhalos near images of strongly lensed extended background galaxies. Typically, only the most massive subhalos in the strong lensing region may be detected this way. In this work, we show that strong lenses may also be used to constrain the much more numerous population of lower mass subhalos that are too small to be detected individually. In particular, we show that the power spectrum of projected density fluctuations in galaxy halos can be measured using strong gravitational lensing. We develop the mathematical framework of power spectrum estimation, and test our method on mock observations. We use our results to determine the types of observations required to measure the substructure power spectrum with high significance. We predict that deep observations ($\sim10$ hours on a single target) with current facilities can measure this power spectrum at the $3\sigma$ level, with no apparent degeneracy with unknown clumpiness in the background source structure or fluctuations from detector noise. Upcoming ALMA measurements of strong lenses are capable of placing strong constraints on the abundance of dark matter subhalos and the underlying particle nature of dark matter., Comment: 8 pages, 6 figures, submitted to ApJ

We perform cosmological N-body simulations with non-Gaussian initial conditions generated from two independent fields. The dominant contribution to the perturbations comes from a purely Gaussian field, but we allow the second field to have local non-Gaussianity that need not be weak. This scenario allows us to adjust the relative importance of non-Gaussian contributions beyond the skewness, producing a scaling of the higher moments different from (and stronger than) the scaling in the usual single field local ansatz. We compare semi-analytic prescriptions for the non-Gaussian mass function, large scale halo bias, and stochastic bias against the simulation results. We discuss applications of this work to large scale structure measurements that can test a wider range of models for the primordial fluctuations than is usually explored., Comment: 35 pages, 14 figures, minor corrections and clarifications added to the text to match the JCAP published version

Recent observations have accumulated compelling evidence that some short gamma-ray bursts (SGRBs) are associated with the mergers of neutron star (NS) binaries. This would indicate that the SGRB event is associated with a gravitational-wave (GW) signal corresponding to the final inspiral of the compact binary. In addition, the radioactive decay of elements produced in NS binary mergers may result in transients visible in the optical and infrared with peak luminosities on hours-days timescales. Simultaneous observations of the inspiral GWs and signatures in the electromagnetic band may allow us to directly and independently determine both the luminosity distance and redshift to a binary. These standard sirens (the GW analog of standard candles) have the potential to provide an accurate measurement of the low-redshift Hubble flow. In addition, these systems are absolutely calibrated by general relativity, and therefore do not experience the same set of astrophysical systematics found in traditional standard candles, nor do the measurements rely on a distance ladder. We show that 15 observable GW and EM events should allow the Hubble constant to be measured with 5% precision using a network of detectors that includes advanced LIGO and Virgo. Measuring 30 beamed GW-SGRB events could constrain H_0 to better than 1%. When comparing to standard Gaussian likelihood analysis, we find that each event's non-Gaussian posterior in H_0 helps reduce the overall measurement errors in H_0 for an ensemble of NS binary mergers., Comment: To be submitted to ApJ, 8 pages in emulateapj format, 1 table, 2 figures

We investigate how strong lensing of dusty, star-forming galaxies by foreground galaxies can be used as a probe of dark matter halo substructure. We find that spatially resolved spectroscopy of lensed sources allows dramatic improvements to measurements of lens parameters. In particular we find that modeling of the full, three-dimensional (angular position and radial velocity) data can significantly facilitate substructure detection, increasing the sensitivity of observables to lower mass subhalos. We carry out simulations of lensed dusty sources observed by early ALMA (Cycle 1) and use a Fisher matrix analysis to study the parameter degeneracies and mass detection limits of this method. We find that, even with conservative assumptions, it is possible to detect galactic dark matter subhalos of ~ 10^8 M_{\odot} with high significance in most lensed DSFGs. Specifically, we find that in typical DSFG lenses, there is a ~ 55 % probability of detecting a substructure with M>10^8 M_{\odot} with more than 5 sigma detection significance in each lens, if the abundance of substructure is consistent with previous lensing results. The full ALMA array, with its significantly enhanced sensitivity and resolution, should improve these estimates considerably. Given the sample of ~100 lenses provided by surveys like the South Pole Telescope, our understanding of dark matter substructure in typical galaxy halos is poised to improve dramatically over the next few years., Comment: 11 pages, 6 figures, submitted to the Astrophysical Journal

We investigate the effect of supersonic relative velocities between baryons and dark matter, recently shown to arise generically at high redshift, on baryonic acoustic oscillation (BAO) measurements at low redshift. The amplitude of the relative velocity effect at low redshift is model-dependent, but can be parameterized by using an unknown bias. We find that if unaccounted, the relative velocity effect can shift the BAO peak position and bias estimates of the dark energy equation-of-state due to its non-smooth, out-of-phase oscillation structure around the BAO scale. Fortunately, the relative velocity effect can be easily modeled in constraining cosmological parameters without substantially inflating the error budget. We also demonstrate that the presence of the relative velocity effect gives rise to a unique signature in the galaxy bispectrum, which can be utilized to isolate this effect. Future dark energy surveys can accurately measure the relative velocity effect and subtract it from the power spectrum analysis to constrain dark energy models with high precision., Comment: 17 pages, 6 figures, submitted to JCAP

The inspirals and mergers of compact binaries are among the most promising events for ground-based gravitational-wave (GW) observatories. The detection of electromagnetic (EM) signals from these sources would provide complementary information to the GW signal. It is therefore important to determine the ability of gravitational-wave detectors to localize compact binaries on the sky, so that they can be matched to their EM counterparts. We use Markov Chain Monte Carlo techniques to study sky localization using networks of ground-based interferometers. Using a coherent-network analysis, we find that the Laser Interferometer Gravitational Wave Observatory (LIGO)-Virgo network can localize 50% of their ~8 sigma detected neutron star binaries to better than 50 sq.deg. with 95% confidence region. The addition of the Large Scale Cryogenic Gravitational Wave Telescope (LCGT) and LIGO-Australia improves this to 12 sq.deg.. Using a more conservative coincident detection threshold, we find that 50% of detected neutron star binaries are localized to 13 sq.deg. using the LIGO-Virgo network, and to 3 sq.deg. using the LIGO-Virgo-LCGT-LIGO-Australia network. Our findings suggest that the coordination of GW observatories and EM facilities offers great promise., Comment: 6 pages, 4 figures, 1 table, matches published version in ApJ (incorporates referee's comments)

The friends-of-friends algorithm (hereafter, FOF) is a percolation algorithm which is routinely used to identify dark matter halos from N-body simulations. We use results from percolation theory to show that the boundary of FOF halos does not correspond to a single density threshold but to a range of densities close to a critical value that depends upon the linking length parameter, b. We show that for the commonly used choice of b = 0.2, this critical density is equal to 81.62 times the mean matter density. Consequently, halos identified by the FOF algorithm enclose an average overdensity which depends on their density profile (concentration) and therefore changes with halo mass contrary to the popular belief that the average overdensity is ~180. We derive an analytical expression for the overdensity as a function of the linking length parameter b and the concentration of the halo. Results of tests carried out using simulated and actual FOF halos identified in cosmological simulations show excellent agreement with our analytical prediction. We also find that the mass of the halo that the FOF algorithm selects crucially depends upon mass resolution. We find a percolation theory motivated formula that is able to accurately correct for the dependence on number of particles for the mock realizations of spherical and triaxial Navarro-Frenk-White halos. However, we show that this correction breaks down when applied to the real cosmological FOF halos due to presence of substructures. Given that abundance of substructure depends on redshift and cosmology, we expect that the resolution effects due to substructure on the FOF mass and halo mass function will also depend on redshift and cosmology and will be difficult to correct for in general. Finally, we discuss the implications of our results for the universality of the mass function., Comment: 19 pages, 17 figures, submitted to ApJ supplement

We investigate the collapse and internal structure of dark matter halos. We consider halo formation from initially scale-free perturbations, for which gravitational collapse is self-similar. Fillmore and Goldreich (1984) and Bertschinger (1985) solved the one dimensional (i.e. spherically symmetric) case. We generalize their results by formulating the three dimensional self-similar equations. We solve the equations numerically and analyze the similarity solutions in detail, focusing on the internal density profiles of the collapsed halos. By decomposing the total density into subprofiles of particles that collapse coevally, we identify two effects as the main determinants of the internal density structure of halos: adiabatic contraction and the shape of a subprofile shortly after collapse; the latter largely reflects the triaxiality of the subprofile. We develop a simple model that describes the results of our 3D simulations. In a companion paper, we apply this model to more realistic cosmological fluctuations, and thereby explain the origin of the nearly universal (NFW-like) density profiles found in N-body simulations., Comment: corresponds to version published in ApJ

Motivated by the properties of early universe scenarios that produce observationally large local non-Gaussianity, we perform N-body simulations with non-Gaussian initial conditions from a generalized local ansatz. The bispectra are schematically of the local shape, but with scale-dependent amplitude. We find that in such cases the size of the non-Gaussian correction to the bias of small and large mass objects depends on the amplitude of non-Gaussianity roughly on the scale of the object. In addition, some forms of the generalized bispectrum alter the scale dependence of the non-Gaussian term in the bias by a fractional power of k. These features may allow significant observational constraints on the particle physics origin of any observed local non-Gaussianity, distinguishing between scenarios where a single field or multiple fields contribute to the curvature fluctuations. While analytic predictions for the non-Gaussian bias agree qualitatively with the simulations, we find numerically a stronger observational signal than expected. This suggests that a more precise understanding of halo formation is needed to fully explain the consequences of primordial non-Gaussianity, Comment: 29 pages, 5 figures; References added

The apparent presence of large core radii in Low Surface Brightness galaxies has been claimed as evidence in favor of warm dark matter. Here we show that WDM halos do not have cores that are large fractions of the halo size: typically, r_core/r_200 < 0.001. This suggests an astrophysical origin for the large cores observed in these galaxies, as has been argued by other authors., Comment: 8 pages, 9 figures

A longstanding puzzle of fundamental importance in modern cosmology has been the origin of the nearly universal density profiles of dark matter halos found in N-body simulations -- the so-called NFW profile. We show how this behavior may be understood, simply, by applying adiabatic contraction to peaks of Gaussian random fields. We argue that dynamical friction acts to reduce enormously the effect of random scatter in the properties of initial peaks, providing a key simplification. We compare our model predictions with results of the ultra-high resolution Via Lactea-II N-body simulation, and find superb agreement. We show how our model may be used to predict the distribution of halo properties like concentration. Our results suggest that many of the basic properties of halo structure may be understood using extremely simple physics., Comment: v2: updated references, submitted to ApJ

Recent work has shown that at high redshift, the relative velocity between dark matter and baryonic gas is typically supersonic. This relative velocity suppresses the formation of the earliest baryonic structures like minihalos, and the suppression is modulated on large scales. This effect imprints a characteristic shape in the clustering power spectrum of the earliest structures, with significant power on 100 Mpc scales featuring highly pronounced baryon acoustic oscillations. The amplitude of these oscillations is orders of magnitude larger at z=20 than previously expected. This characteristic signature can allow us to distinguish the effects of minihalos on intergalactic gas at times preceding and during reionization. We illustrate this effect with the example of 21 cm emission and absorption from redshifts during and before reionization. This effect can potentially allow us to probe physics on kpc scales using observations on 100 Mpc scales. We present sensitivity forecasts for FAST and Arecibo. Depending on parameters, this enhanced structure may be detectable by Arecibo at redshifts near z=15-20, and with appropriate instrumentation FAST could measure the BAO power spectrum with high precision. In principle, this effect could also pose a serious challenge for efforts to constrain dark energy using observations of the BAO feature at low redshift., Comment: 9 pages, 5 figures

Recent observations support the hypothesis that a large fraction of "short-hard" gamma-ray bursts (SHBs) are associated with compact binary inspiral. Since gravitational-wave (GW) measurements of well-localized inspiraling binaries can measure absolute source distances, simultaneous observation of a binary's GWs and SHB would allow us to independently determine both its luminosity distance and redshift. Such a "standard siren" (the GW analog of a standard candle) would provide an excellent probe of the relatively nearby universe's expansion, complementing other standard candles. In this paper, we examine binary measurement using a Markov Chain Monte Carlo technique to build the probability distributions describing measured parameters. We assume that each SHB observation gives both sky position and the time of coalescence, and we take both binary neutron stars and black hole-neutron star coalescences as plausible SHB progenitors. We examine how well parameters particularly distance) can be measured from GW observations of SHBs by a range of ground-based detector networks. We find that earlier estimates overstate how well distances can be measured, even at fairly large signal-to-noise ratio. The fundamental limitation to determining distance proves to be a degeneracy between distance and source inclination. Overcoming this limitation requires that we either break this degeneracy, or measure enough sources to broadly sample the inclination distribution. (Abridged), Comment: 19 pages, 10 figures. Accepted for publication in ApJ; this version incorporates referee's comments and criticisms

There is a vast menagerie of plausible candidates for the constituents of dark matter, both within and beyond extensions of the Standard Model of particle physics. Each of these candidates may have scattering (and other) cross section properties that are consistent with the dark matter abundance, BBN, and the most scales in the matter power spectrum; but which may have vastly different behavior at sub-galactic "cutoff" scales, below which dark matter density fluctuations are smoothed out. The only way to quantitatively measure the power spectrum behavior at sub-galactic scales at distances beyond the local universe, and indeed over cosmic time, is through probes available in multiply imaged strong gravitational lenses. Gravitational potential perturbations by dark matter substructure encode information in the observed relative magnifications, positions, and time delays in a strong lens. Each of these is sensitive to a different moment of the substructure mass function and to different effective mass ranges of the substructure. The time delay perturbations, in particular, are proving to be largely immune to the degeneracies and systematic uncertainties that have impacted exploitation of strong lenses for such studies. There is great potential for a coordinated theoretical and observational effort to enable a sophisticated exploitation of strong gravitational lenses as direct probes of dark matter properties. This opportunity motivates this white paper, and drives the need for: a) strong support of the theoretical work necessary to understand all astrophysical consequences for different dark matter candidates; and b) tailored observational campaigns, and even a fully dedicated mission, to obtain the requisite data., Comment: Science white paper submitted to the Astro2010 Decadal Cosmology & Fundamental Physics Science Frontier Panel

We derive radial mass profiles of four strong lensing selected clusters which show prominent giant arcs (Abell 1703, SDSS J1446+3032, SDSS J1531+3414, and SDSS J2111-0115), by combining detailed strong lens modeling with weak lensing shear measured from deep Subaru Suprime-cam images. Weak lensing signals are detected at high significance for all four clusters, whose redshifts range from z=0.28 to 0.64. We demonstrate that adding strong lensing information with known arc redshifts significantly improves constraints on the mass density profile, compared to those obtained from weak lensing alone. While the mass profiles are well fitted by the universal form predicted in N-body simulations of the LCDM model, all four clusters appear to be slightly more centrally concentrated (the concentration parameters c_vir \sim 8) than theoretical predictions, even after accounting for the bias toward higher concentrations inherent in lensing selected samples. Our results are consistent with previous studies which similarly detected a concentration excess, and increases the total number of clusters studied with the combined strong and weak lensing technique to ten. Combining our sample with previous work, we find that clusters with larger Einstein radii are more anomalously concentrated. We also present a detailed model of the lensing cluster A1703 with constraints from multiple image families, and find the dark matter inner density profile to be cuspy with the slope consistent with -1, in agreement with expectations., Comment: 41 pages, 9 figures, 7 tables, accepted for publication in ApJ

We investigate the origin of halo assembly bias, the dependence of halo clustering on assembly history. We relate halo assembly to peak properties measured in the Lagrangian space of the initial linear Gaussian random density field, and show how these same Lagrangian properties determine large-scale bias. We focus on the two regimes where assembly bias has been observed to be significant: at masses very large and very small compared to the nonlinear mass scale. At high masses, we show that assembly bias is expected from the statistics of the peaks of Gaussian random fluctuations, and we show that the extent of assembly bias found in N-body simulations of rare halos is in excellent agreement with our theoretical prediction. At low masses, we argue that assembly bias largely arises from a sub-population of low mass halos whose mass accretion has ceased. Due to their arrested development, these halos naturally become unbiased, in contrast to their anti-biased peers. We show that a simple toy model incorporating these effects can roughly reproduce the bias trends found in N-body simulations., Comment: 10 pages, color figures, downgraded resolution

We study the effect of primordial nongaussianity on large-scale structure, focusing upon the most massive virialized objects. Using analytic arguments and N-body simulations, we calculate the mass function and clustering of dark matter halos across a range of redshifts and levels of nongaussianity. We propose a simple fitting function for the mass function valid across the entire range of our simulations. We find pronounced effects of nongaussianity on the clustering of dark matter halos, leading to strongly scale-dependent bias. This suggests that the large-scale clustering of rare objects may provide a sensitive probe of primordial nongaussianity. We very roughly estimate that upcoming surveys can constrain nongaussianity at the level |fNL| <~ 10, competitive with forecasted constraints from the microwave background., Comment: 16 pages, color figures, revtex4. v2: added references and an equation. submitted to PRD. v3: simplified derivation, additional references

We report on the first results of an imaging survey to detect strong gravitational lensing targeting the richest clusters selected from the photometric data of the Sloan Digital Sky Survey (SDSS) with follow-up deep imaging observations from the Wisconsin Indiana Yale NOAO (WIYN) 3.5m telescope and the University of Hawaii 88-inch telescope (UH88). The clusters are selected from an area of 8000 deg^2 using the Red Cluster Sequence technique and span the redshift range 0.1 < z < 0.6, corresponding to a comoving cosmological volume of ~ 2 Gpc^3. Our imaging survey thus targets a volume more than an order of magnitude larger than any previous search. A total of 240 clusters were imaged of which 141 had sub-arcsecond image quality. Our survey has uncovered16 new lensing clusters with definite giant arcs, an additional 12 systems for which the lensing interpretation is very likely, and 9 possible lenses which contain shorter arclets or candidate arcs which are less certain and will require further observations to confirm their lensing origin. The number of new cluster lenses detected in this survey is likely > 30. Among these new systems are several of the most dramatic examples of strong gravitational lensing ever discovered with multiple bright arcs at large angular separation. These will likely become 'poster-child' gravitational lenses similar to Abell 1689 and CL0024+1654. The new lenses discovered in this survey will enable future sysetmatic studies of the statistics of strong lensing and its implications for cosmology and our structure formation paradigm., Comment: 19 pages, 7 pages of Figures, submitted to AJ. Fixed Typos

In recent years, gravitational lensing has been used as a means to detect substructure in galaxy-sized halos, via anomalous flux ratios in quadruply-imaged lenses. In addition to causing anomalous flux ratios, substructure may also perturb the positions of lensed images at observable levels. In this paper, we numerically investigate the scale of such astrometric perturbations using realistic models of substructure distributions. Substructure distributions that project clumps near the Einstein radius of the lens result in perturbations that are the least degenerate with the best-fit smooth macromodel, with residuals at the milliarcsecond scale. Degeneracies between the center of the lens potential and astrometric perturbations suggest that milliarcsecond constraints on the center of the lensing potential boost the observed astrometric perturbations by an order of magnitude compared to leaving the center of the lens as a free parameter. In addition, we discuss methods of substructure detection via astrometric perturbations that avoid full lens modeling in favor of local image observables and also discuss modeling of systems with luminous satellites to constrain the masses of those satellites., Comment: 15 pages, 13 figures, ApJ accepted version. 2 new figures, additional section on halo-to-halo variation

Observations of the gravitational radiation from well-localized, inspiraling compact object binaries can measure absolute source distances with high accuracy. When coupled with an independent determination of redshift through an electromagnetic counterpart, these standard sirens can provide an excellent probe of the expansion history of the Universe and the dark energy. Short gamma-ray bursts, if produced by merging neutron star binaries, would be standard sirens with known redshifts detectable by ground-based GW networks such as LIGO-II, Virgo, and AIGO. Depending upon the collimation of these GRBs, a single year of observation of their gravitational waves can measure the Hubble constant to about 2%. When combined with measurement of the absolute distance to the last scattering surface of the cosmic microwave background, this determines the dark energy equation of state parameter w to 9%. Similarly, supermassive binary black hole inspirals will be standard sirens detectable by LISA. Depending upon the precise redshift distribution, 100 sources could measure w at the 4% level., Comment: 8 pages, submitted to PRD