Showing 1–50 of 138 results for author: Hui, L

Massive-ish Particles from Small-ish Scales: Non-Perturbative Techniques for Cosmological Collider Physics from Large-Scale Structure Surveys

Authors: Samuel Goldstein, Oliver H. E. Philcox, J. Colin Hill, Lam Hui

Abstract: Massive particles produced during inflation impact soft limits of primordial correlators. Searching for these signatures presents an exciting opportunity to uncover the particle spectrum in the inflationary epoch. We present non-perturbative methods to constrain intermediate-mass scalars ($0\leq m/H<3/2$, where $H$ is the inflationary Hubble scale) produced during inflation, which give rise to a p… ▽ More Massive particles produced during inflation impact soft limits of primordial correlators. Searching for these signatures presents an exciting opportunity to uncover the particle spectrum in the inflationary epoch. We present non-perturbative methods to constrain intermediate-mass scalars ($0\leq m/H<3/2$, where $H$ is the inflationary Hubble scale) produced during inflation, which give rise to a power-law scaling in the squeezed primordial bispectrum. Exploiting the large-scale structure consistency relations and the separate universe approach, we derive models for the late-time squeezed matter bispectrum and collapsed matter trispectrum sourced by these fields. To validate our models, we run $N$-body simulations with the "Cosmological Collider" squeezed bispectrum for two different particle masses. Our models yield unbiased constraints on the amplitude of non-Gaussianity, $f_{\rm NL}^Δでるた$, from the squeezed bispectrum and collapsed trispectrum deep into the non-linear regime ($k_{\rm max}\approx 2~h/{\rm Mpc}$ at $z=0$). We assess the information content of these summary statistics, emphasizing the importance of sample variance cancellation in the matter sector. We also study the scale-dependent halo bias in our simulations. For mass-selected halos, the non-Gaussian bias estimated from our simulations agrees with predictions based on (i) separate universe simulations and (ii) universal mass functions. With further work, these results can be used to search for inflationary massive particle production with upcoming galaxy surveys. △ Less

Submitted 11 July, 2024; originally announced July 2024.

Comments: 25 pages, 11 figures; comments welcome

Fragmention in Gravitationally-Unstable Collapsar Disks and Sub-Solar Neutron Star Mergers

Authors: Brian D. Metzger, Lam Hui, Matteo Cantiello

Abstract: Although stable neutron stars (NS) can in principle exist down to masses Mns ~ 0.1Msun, standard models of stellar core-collapse predict a robust lower limit Mns >~ 1.2Msun, roughly commensurate with the Chandrasekhar mass Mch of the progenitor's iron core (electron fraction Ye ~ 0.5). However, this limit may be circumvented in sufficiently dense neutron-rich environments (Ye << 0.5) for which Mch… ▽ More Although stable neutron stars (NS) can in principle exist down to masses Mns ~ 0.1Msun, standard models of stellar core-collapse predict a robust lower limit Mns >~ 1.2Msun, roughly commensurate with the Chandrasekhar mass Mch of the progenitor's iron core (electron fraction Ye ~ 0.5). However, this limit may be circumvented in sufficiently dense neutron-rich environments (Ye << 0.5) for which Mch ~ Ye^2 is reduced to < Msun. Such physical conditions could arise in the black hole accretion disks formed from the collapse of rapidly-rotating stars (``collapsars''), as a result of gravitational instabilities and cooling-induced fragmentation, similar to models for planet formation in protostellar disks. We confirm that the conditions to form sub-solar mass NS (ssNS) may be marginally satisfied in the outer regions of massive neutrino-cooled collapsar disks. If the disk fragments into multiple ssNS, their subsequent coalescence offers a channel for precipitating sub-solar mass LIGO/Virgo gravitational-wave mergers that does not implicate primordial black holes. The model makes several additional predictions: (1) ~Hz frequency Doppler modulation of the ssNS-merger gravitational wave signals due to the binary's orbital motion in the disk; (2) at least one additional gravitational wave event (coincident within <~ hours), from the coalescence of the ssNS-merger remnant(s) with the central black hole; (3) an associated gamma-ray burst and supernova counterpart, the latter boosted in energy and enriched with r-process elements from the NS merger(s) embedded within the exploding stellar envelope (``kilonovae inside a supernova''). △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: 9 pages, 2 figures

Consistently constraining $f_{\rm NL}$ with the squeezed lensing bispectrum using consistency relations

Authors: Samuel Goldstein, Oliver H. E. Philcox, J. Colin Hill, Angelo Esposito, Lam Hui

Abstract: We introduce a non-perturbative method to constrain the amplitude of local-type primordial non-Gaussianity ($f_{\rm NL}$) using squeezed configurations of the CMB lensing convergence and cosmic shear bispectra. First, we use cosmological consistency relations to derive a model for the squeezed limit of angular auto- and cross-bispectra of lensing convergence fields in the presence of $f_{\rm NL}$.… ▽ More We introduce a non-perturbative method to constrain the amplitude of local-type primordial non-Gaussianity ($f_{\rm NL}$) using squeezed configurations of the CMB lensing convergence and cosmic shear bispectra. First, we use cosmological consistency relations to derive a model for the squeezed limit of angular auto- and cross-bispectra of lensing convergence fields in the presence of $f_{\rm NL}$. Using this model, we perform a Fisher forecast with specifications expected for upcoming CMB lensing measurements from the Simons Observatory and CMB-S4, as well as cosmic shear measurements from a Rubin LSST/Euclid-like experiment. Assuming a minimum multipole $\ell_{\rm min}=10$ and maximum multipole $\ell_{\rm max}=1400$, we forecast $σしぐま_{f_{\rm NL}}=175$ ($95$) for Simons Observatory (CMB-S4). Our forecasts improve considerably for an LSST/Euclid-like cosmic shear experiment with three tomographic bins and $\ell_{\rm min}=10$ and $\ell_{\rm max}=1400$ ($5000$) with $σしぐま_{f_{\rm NL}}=31$ ($16$). A joint analysis of CMB-S4 lensing and LSST/Euclid-like shear yields little gain over the shear-only forecasts; however, we show that a joint analysis could be useful if the CMB lensing convergence can be reliably reconstructed at larger angular scales than the shear field. The method presented in this work is a novel and robust technique to constrain local primordial non-Gaussianity from upcoming large-scale structure surveys that is completely independent of the galaxy field (and therefore any nuisance parameters such as $b_φふぁい$), thus complementing existing techniques to constrain $f_{\rm NL}$ using the scale-dependent halo bias. △ Less

Submitted 19 October, 2023; originally announced October 2023.

Comments: 14 pages, 5 figures, comments welcome!

Statistical analysis of the onset temperature of solar flares in 2010-2011

Authors: Douglas Félix da Silva, Li Hui, Paulo J. A. Simões, Adriana Valio, Joaquim C. E. R., Hugh S. Hudson, Paulo J. A. Simoes, Lyndsay Fletcher, Laura A. Hayes, Iain G. Hannah

Abstract: Understanding the physical processes that trigger solar flares is paramount to help with forecasting space weather and mitigating the effects on our technological infrastructure. A previously unknown phenomenon was recently identified in solar flares: the plasma temperature, derived from soft X-ray (SXR) data, at the onset of four flares, was revealed to be in the range 10-15 MK, without evidence… ▽ More Understanding the physical processes that trigger solar flares is paramount to help with forecasting space weather and mitigating the effects on our technological infrastructure. A previously unknown phenomenon was recently identified in solar flares: the plasma temperature, derived from soft X-ray (SXR) data, at the onset of four flares, was revealed to be in the range 10-15 MK, without evidence of gradual heating. To investigate how common the hot-onset phenomenon may be, we extend this investigation to solar flares of B1.2- X6.9 classes recorded by the X-ray Sensor (XRS) on-board the GOES-14 and GOES-15 satellites between 2010 and 2011. For this statistical study, we employed the same methodology as in recent work, where the pre-flare SXR flux of each flare is obtained manually, and the temperature and emission measure values are obtained by the flux ratio of the two GOES/XRS channels using the standard software. From 3224 events listed in the GOES flare catalog for 2010-2011, we have selected and analyzed 745 events for which the flare heliographic location was provided in the list, to investigate center-to-limb effects of the hot-onset phenomenon. Our results show that 559 out of 745 flares (75%) exhibit an onset temperature above 8.6 MK (the first quartile), with respective log10 of the emission measure values between 46.0 - 47.25 cm-3, indicating that small amounts of plasma are quickly heated to high temperatures. These results suggest that the hot-onset phenomenon is very common in solar flares. △ Less

Submitted 21 August, 2023; originally announced August 2023.

Comments: 6 pages,7 figures

Relativistic drag forces on black holes from scalar dark matter clouds of all sizes

Authors: Dina Traykova, Rodrigo Vicente, Katy Clough, Thomas Helfer, Emanuele Berti, Pedro G. Ferreira, Lam Hui

Abstract: We use numerical simulations of scalar field dark matter evolving on a moving black hole background to confirm the regime of validity of (semi-)analytic expressions derived from first principles for both dynamical friction and momentum accretion in the relativistic regime. We cover both small and large clouds (relative to the de Broglie wavelength of the scalars), and light and heavy particle mass… ▽ More We use numerical simulations of scalar field dark matter evolving on a moving black hole background to confirm the regime of validity of (semi-)analytic expressions derived from first principles for both dynamical friction and momentum accretion in the relativistic regime. We cover both small and large clouds (relative to the de Broglie wavelength of the scalars), and light and heavy particle masses (relative to the BH size). In the case of a small dark matter cloud, the effect of accretion is a non-negligible contribution to the total force on the black hole, even for small scalar masses. We confirm that this momentum accretion transitions between two regimes (wave- and particle-like) and we identify the mass of the scalar at which the transition between regimes occurs. △ Less

Submitted 19 February, 2024; v1 submitted 17 May, 2023; originally announced May 2023.

Comments: 11 pages, 5 figures. Minor corrections and references added to match published version

Journal ref: Phys.Rev.D 108 (2023) 12, L121502

Soft theorems for boosts and other time symmetries

Authors: Lam Hui, Austin Joyce, Ilia Komissarov, Klaas Parmentier, Luca Santoni, Sam S. C. Wong

Abstract: We derive soft theorems for theories in which time symmetries -- symmetries that involve the transformation of time, an example of which are Lorentz boosts -- are spontaneously broken. The soft theorems involve unequal-time correlation functions with the insertion of a soft Goldstone in the far past. Explicit checks are provided for several examples, including the effective theory of a relativisti… ▽ More We derive soft theorems for theories in which time symmetries -- symmetries that involve the transformation of time, an example of which are Lorentz boosts -- are spontaneously broken. The soft theorems involve unequal-time correlation functions with the insertion of a soft Goldstone in the far past. Explicit checks are provided for several examples, including the effective theory of a relativistic superfluid and the effective field theory of inflation. We discuss how in certain cases these unequal-time identities capture information at the level of observables that cannot be seen purely in terms of equal-time correlators of the field alone. We also discuss when it is possible to phrase these soft theorems as identities involving equal-time correlators. △ Less

Submitted 28 October, 2022; originally announced October 2022.

Comments: 50 pages

Squeezing $f_{\rm NL}$ out of the matter bispectrum with consistency relations

Authors: Samuel Goldstein, Angelo Esposito, Oliver H. E. Philcox, Lam Hui, J. Colin Hill, Roman Scoccimarro, Maximilian H. Abitbol

Abstract: We show how consistency relations can be used to robustly extract the amplitude of local primordial non-Gaussianity ($f_{\rm NL}$) from the squeezed limit of the matter bispectrum, well into the non-linear regime. First, we derive a non-perturbative relation between primordial non-Gaussianity and the leading term in the squeezed bispectrum, revising some results present in the literature. This rel… ▽ More We show how consistency relations can be used to robustly extract the amplitude of local primordial non-Gaussianity ($f_{\rm NL}$) from the squeezed limit of the matter bispectrum, well into the non-linear regime. First, we derive a non-perturbative relation between primordial non-Gaussianity and the leading term in the squeezed bispectrum, revising some results present in the literature. This relation is then used to successfully measure $f_{\rm NL}$ from $N$-body simulations. We discuss the dependence of our results on different scale cuts and redshifts. Specifically, the analysis is strongly dependent on the choice of the smallest soft momentum, $q_{\rm min}$, which is the most sensitive to primordial bispectrum contributions, but is largely independent of the choice of the largest hard momentum, $k_{\rm max}$, due to the non-Gaussian nature of the covariance. We also show how the constraints on $f_{\rm NL}$ improve at higher redshift, due to a reduced off-diagonal covariance. In particular, for a simulation with $f_{\rm NL} = 100$ and a volume of $(2.4 \text{ Gpc}/h)^3$, we measure $f_{\rm NL} = 98 \pm 12$ at redshift $z=0$ and $f_{\rm NL} = 97 \pm 8$ at $z=0.97$. Finally, we compare our results with a Fisher forecast, showing that the current version of the analysis is satisfactorily close to the Fisher error. We regard this as a first step towards the realistic application of consistency relations to constrain primordial non-Gaussianity using observations. △ Less

Submitted 6 January, 2023; v1 submitted 13 September, 2022; originally announced September 2022.

Comments: 17 pages, 8 figures. Minor changes. Matches version published in PRD

Journal ref: Phys. Rev. D 106, 123525 (2022)

Nonlinearities in Black Hole Ringdowns

Authors: Keefe Mitman, Macarena Lagos, Leo C. Stein, Sizheng Ma, Lam Hui, Yanbei Chen, Nils Deppe, François Hébert, Lawrence E. Kidder, Jordan Moxon, Mark A. Scheel, Saul A. Teukolsky, William Throwe, Nils L. Vu

Abstract: The gravitational wave strain emitted by a perturbed black hole (BH) ringing down is typically modeled analytically using first-order BH perturbation theory. In this Letter we show that second-order effects are necessary for modeling ringdowns from BH merger simulations. Focusing on the strain's $(\ell,m)=(4,4)$ angular harmonic, we show the presence of a quadratic effect across a range of binary… ▽ More The gravitational wave strain emitted by a perturbed black hole (BH) ringing down is typically modeled analytically using first-order BH perturbation theory. In this Letter we show that second-order effects are necessary for modeling ringdowns from BH merger simulations. Focusing on the strain's $(\ell,m)=(4,4)$ angular harmonic, we show the presence of a quadratic effect across a range of binary BH mass ratios that agrees with theoretical expectations. We find that the quadratic $(4,4)$ mode's amplitude exhibits quadratic scaling with the fundamental $(2,2)$ mode -- its parent mode. The nonlinear mode's amplitude is comparable to or even larger than that of the linear $(4,4)$ mode. Therefore, correctly modeling the ringdown of higher harmonics -- improving mode mismatches by up to 2 orders of magnitude -- requires the inclusion of nonlinear effects. △ Less

Submitted 22 February, 2023; v1 submitted 15 August, 2022; originally announced August 2022.

Comments: 6+2 pages, 4 figures, 1 table. Matches PRL version

Journal ref: Phys. Rev. Lett. 130, 081402 (2023)

Generation and propagation of nonlinear quasi-normal modes of a Schwarzschild black hole

Authors: Macarena Lagos, Lam Hui

Abstract: In the analysis of a binary black hole coalescence, it is necessary to include gravitational self-interactions in order to describe the transition of the gravitational wave signal from the merger to the ringdown stage. In this paper we study the phenomenology of the generation and propagation of nonlinearities in the ringdown of a Schwarzschild black hole, using second-order perturbation theory. F… ▽ More In the analysis of a binary black hole coalescence, it is necessary to include gravitational self-interactions in order to describe the transition of the gravitational wave signal from the merger to the ringdown stage. In this paper we study the phenomenology of the generation and propagation of nonlinearities in the ringdown of a Schwarzschild black hole, using second-order perturbation theory. Following earlier work, we show that the Green's function and its causal structure determines how both first-order and second-order perturbations are generated, and hence highlight that both of these solutions share some physical properties. In particular, we discuss the sense in which both linear and quadratic quasi-normal modes (QNMs) are generated in the vicinity of the peak of the gravitational potential barrier (loosely referred to as the light ring). Among the second-order perturbations, there are solutions with linear QNM frequencies (whose amplitudes are thus renormalized from their linear values), as well as quadratic QNM frequencies with a distinct spectrum. Moreover, we show using a WKB analysis that, in the eikonal limit, waves generated inside the light ring propagate towards the black hole horizon, and only waves generated outside propagate towards an asymptotic observer. These results might be relevant for recent discussions on the validity of perturbation theory close to the merger. Finally, we argue that even if nonlinearities are small, quadratic QNMs may be detectable and would likely be useful for improving ringdown models of higher angular harmonics and future tests of gravity. △ Less

Submitted 9 January, 2023; v1 submitted 15 August, 2022; originally announced August 2022.

Comments: Version accepted in PRD

Black hole superradiance with (dark) matter accretion

Authors: Lam Hui, Y. T. Albert Law, Luca Santoni, Guanhao Sun, Giovanni Maria Tomaselli, Enrico Trincherini

Abstract: Studies of black hole superradiance often focus on the growth of a cloud in isolation, accompanied by the spin-down of the black hole. In this paper, we consider the additional effect of the accretion of matter and angular momentum from the environment. We show that, in many cases, the black hole evolves by drifting along the superradiance threshold, in which case the evolution of its parameters c… ▽ More Studies of black hole superradiance often focus on the growth of a cloud in isolation, accompanied by the spin-down of the black hole. In this paper, we consider the additional effect of the accretion of matter and angular momentum from the environment. We show that, in many cases, the black hole evolves by drifting along the superradiance threshold, in which case the evolution of its parameters can be described analytically or semi-analytically. We quantify the conditions under which accretion can serve as a mechanism to increase the cloud-to-black hole mass ratio, beyond the standard maximum of about 10%. This occurs by a process we call over-superradiance, whereby accretion effectively feeds the superradiance cloud, by way of the black hole. We give two explicit examples: accretion from a vortex expected in wave dark matter and accretion from a baryonic disk. In the former case, we estimate the accretion rate by using an analytical fit to the asymptotic behavior of the confluent Heun function. Level transition, whereby one cloud level grows while the other shrinks, can be understood in a similar way. △ Less

Submitted 25 May, 2023; v1 submitted 12 August, 2022; originally announced August 2022.

Comments: 30+21 pages, 14 figures

Near-Zone Symmetries of Kerr Black Holes

Authors: Lam Hui, Austin Joyce, Riccardo Penco, Luca Santoni, Adam R. Solomon

Abstract: We study the near-zone symmetries of a massless scalar field on four-dimensional black hole backgrounds. We provide a geometric understanding that unifies various recently discovered symmetries as part of an SO(4,2) group. Of these, a subset are exact symmetries of the static sector and give rise to the ladder symmetries responsible for the vanishing of Love numbers. In the Kerr case, we compare d… ▽ More We study the near-zone symmetries of a massless scalar field on four-dimensional black hole backgrounds. We provide a geometric understanding that unifies various recently discovered symmetries as part of an SO(4,2) group. Of these, a subset are exact symmetries of the static sector and give rise to the ladder symmetries responsible for the vanishing of Love numbers. In the Kerr case, we compare different near-zone approximations in the literature, and focus on the implementation that retains the symmetries of the static limit. We also describe the relation to spin-1 and 2 perturbations. △ Less

Submitted 16 March, 2022; originally announced March 2022.

Comments: 4+3 pages, 1 figure

Journal ref: JHEP09(2022)049

Construction of Wave Dark Matter Halos: Numerical Algorithm and Analytical Constraints

Authors: Tomer D. Yavetz, Xinyu Li, Lam Hui

Abstract: We present a wave generalization of the classic Schwarzschild method for constructing self-consistent halos -- such a halo consists of a suitable superposition of waves instead of particle orbits, chosen to yield a desired mean density profile. As an illustration, the method is applied to spherically symmetric halos. We derive an analytic relation between the particle distribution function and the… ▽ More We present a wave generalization of the classic Schwarzschild method for constructing self-consistent halos -- such a halo consists of a suitable superposition of waves instead of particle orbits, chosen to yield a desired mean density profile. As an illustration, the method is applied to spherically symmetric halos. We derive an analytic relation between the particle distribution function and the wave superposition amplitudes, and show how it simplifies in the high energy (WKB) limit. We verify the stability of such constructed halos by numerically evolving the Schrödinger-Poisson system. The algorithm provides an efficient and accurate way to simulate the time-dependent halo substructures from wave interference. We use this method to construct halos with a variety of density profiles, all of which have a core from the ground-state wave function, though the core-halo relation need not be the standard one. △ Less

Submitted 5 January, 2023; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: 22 pages, 15 figures; published in Phys. Rev. D

Journal ref: Phys. Rev. D 105, 023512, Published 10 January 2022

Dynamical friction from scalar dark matter in the relativistic regime

Authors: Dina Traykova, Katy Clough, Thomas Helfer, Emanuele Berti, Pedro G. Ferreira, Lam Hui

Abstract: Light bosonic scalars (e.g. axions) may form clouds around black holes via superradiant instabilities, or via accretion if they form some component of the dark matter. It has been suggested that their presence may lead to a distinctive dephasing of the gravitational wave signal when a small compact object spirals into a larger black hole. Motivated by this, we study numerically the dynamical frict… ▽ More Light bosonic scalars (e.g. axions) may form clouds around black holes via superradiant instabilities, or via accretion if they form some component of the dark matter. It has been suggested that their presence may lead to a distinctive dephasing of the gravitational wave signal when a small compact object spirals into a larger black hole. Motivated by this, we study numerically the dynamical friction force on a black hole moving at relativistic velocities in a background scalar field with an asymptotically homogeneous energy density. We show that the relativistic scaling is analogous to that found for supersonic collisional fluids, assuming an approximate expression for the pressure correction which depends on the velocity and scalar mass. While we focus on a complex scalar field, our results confirm the expectation that real scalars would exert a force which oscillates between positive and negative values in time with a frequency set by the scalar mass. The complex field describes the time averaged value of this force, but in a real scalar the rapid force oscillations could in principle leave an imprint on the trajectory. The approximation we obtain can be used to inform estimates of dephasing in the final stages of an extreme mass ratio inspiral. △ Less

Submitted 27 October, 2021; v1 submitted 15 June, 2021; originally announced June 2021.

Comments: 17 pages, 13 figures, 1 table; Version as accepted in PRD - minor changes

Ladder Symmetries of Black Holes: Implications for Love Numbers and No-Hair Theorems

Authors: Lam Hui, Austin Joyce, Riccardo Penco, Luca Santoni, Adam R. Solomon

Abstract: It is well known that asymptotically flat black holes in general relativity have a vanishing static, conservative tidal response. We show that this is a result of linearly realized symmetries governing static (spin 0,1,2) perturbations around black holes. The symmetries have a geometric origin: in the scalar case, they arise from the (E)AdS isometries of a dimensionally reduced black hole spacetim… ▽ More It is well known that asymptotically flat black holes in general relativity have a vanishing static, conservative tidal response. We show that this is a result of linearly realized symmetries governing static (spin 0,1,2) perturbations around black holes. The symmetries have a geometric origin: in the scalar case, they arise from the (E)AdS isometries of a dimensionally reduced black hole spacetime. Underlying the symmetries is a ladder structure which can be used to construct the full tower of solutions, and derive their general properties: (1) solutions that decay with radius spontaneously break the symmetries, and must diverge at the horizon; (2) solutions regular at the horizon respect the symmetries, and take the form of a finite polynomial that grows with radius. Taken together, these two properties imply that static response coefficients -- and in particular Love numbers -- vanish. Moreover, property (1) is consistent with the absence of black holes with linear (perturbative) hair. We also discuss the manifestation of these symmetries in the effective point particle description of a black hole, showing explicitly that for scalar probes the worldline couplings associated with a non-trivial tidal response and scalar hair must vanish in order for the symmetries to be preserved. △ Less

Submitted 14 January, 2022; v1 submitted 3 May, 2021; originally announced May 2021.

Comments: Some sentences rephrased for clarity. Equations and conclusions unchanged

Journal ref: JCAP 01 (2022) 01, 032

Wave Dark Matter

Authors: Lam Hui

Abstract: We review the physics and phenomenology of wave dark matter: a bosonic dark matter candidate lighter than about 30 eV. Such particles have a de Broglie wavelength exceeding the average inter-particle separation in a galaxy like the Milky Way, and are well described as classical waves. We outline the particle physics motivations for them, including the QCD axion and ultra-light axion-like-particles… ▽ More We review the physics and phenomenology of wave dark matter: a bosonic dark matter candidate lighter than about 30 eV. Such particles have a de Broglie wavelength exceeding the average inter-particle separation in a galaxy like the Milky Way, and are well described as classical waves. We outline the particle physics motivations for them, including the QCD axion and ultra-light axion-like-particles such as fuzzy dark matter. The wave nature of the dark matter implies a rich phenomenology: (1) Wave interference leads to order unity density fluctuations on de Broglie scale. A manifestation is vortices where the density vanishes and around which the velocity circulates. There is one vortex ring per de Broglie volume on average. (2) For sufficiently low masses, soliton condensation occurs at centers of halos. The soliton oscillates and random walks, another manifestation of wave interference. The halo/subhalo abundance is suppressed at small masses, but the precise prediction from numerical wave simulations remains to be determined. (3) For ultra-light ~$10^{-22}$ eV dark matter, the wave interference substructures can be probed by tidal streams/gravitational lensing. The signal can be distinguished from that due to subhalos by the dependence on stream orbital radius/image separation. (4) Axion detection experiments are sensitive to interference substructures for moderately light masses. The stochastic nature of the waves affects the interpretation of experiments and motivates the measurement of correlation functions. Current constraints and open questions, covering detection experiments and cosmological/galactic/black-hole observations, are discussed. △ Less

Submitted 27 January, 2021; originally announced January 2021.

Comments: 44 pages, to appear in Annual Review of Astronomy and Astrophysics

Don't cross the streams: caustics from Fuzzy Dark Matter

Authors: Neal Dalal, Jo Bovy, Lam Hui, Xinyu Li

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

Submitted 26 November, 2020; originally announced November 2020.

Comments: Comments welcome. To be submitted to JCAP

Oscillations and Random Walk of the Soliton Core in a Fuzzy Dark Matter Halo

Authors: Xinyu Li, Lam Hui, Tomer D. Yavetz

Abstract: A Fuzzy Dark Matter (FDM) halo consists of a soliton core close to the center and an NFW-like density profile in the outer region. Previous investigations found that the soliton core exhibits temporal oscillations and random walk excursions around the halo center. Analyzing a set of numerical simulations, we show that both phenomena can be understood as the results of wave interference -- a suitab… ▽ More A Fuzzy Dark Matter (FDM) halo consists of a soliton core close to the center and an NFW-like density profile in the outer region. Previous investigations found that the soliton core exhibits temporal oscillations and random walk excursions around the halo center. Analyzing a set of numerical simulations, we show that both phenomena can be understood as the results of wave interference -- a suitable superposition of the ground (solitonic) state and excited states in a fixed potential suffices to account for the main features of these phenomena. Such an eigenmode analysis can shed light on the evolution of a satellite halo undergoing tidal disruption. As the outer halo is stripped away, reducing the amplitudes of the excited states, the ground state evolves adiabatically. This suggests diminished soliton oscillations and random walk excursions, an effect to consider in deducing constraints from stellar heating. △ Less

Submitted 23 November, 2020; originally announced November 2020.

Comments: 9 pages, 12 figures

Journal ref: Phys. Rev. D 103, 023508 (2021)

Growth of accretion driven scalar hair around Kerr black holes

Authors: Jamie Bamber, Katy Clough, Pedro G. Ferreira, Lam Hui, Macarena Lagos

Abstract: Scalar fields around compact objects are of interest for scalar-tensor theories of gravity and dark matter models consisting of a massive scalar, e.g. axions. We study the behaviour of a scalar field around a Kerr black hole with non trivial asymptotic boundary conditions - both non zero density and non zero angular momentum. Starting from an initial radially homogeneous configuration, a scalar cl… ▽ More Scalar fields around compact objects are of interest for scalar-tensor theories of gravity and dark matter models consisting of a massive scalar, e.g. axions. We study the behaviour of a scalar field around a Kerr black hole with non trivial asymptotic boundary conditions - both non zero density and non zero angular momentum. Starting from an initial radially homogeneous configuration, a scalar cloud is accreted, which asymptotes to known stationary configurations over time. We study the cloud growth for different parameters including black hole spin, scalar field mass, and the scalar field density and angular momentum far from the black hole. We characterise the transient growth of the mass and angular momentum in the cloud, and the spatial profile of the scalar around the black hole, and relate the results of fully non-linear simulations to an analytic perturbative expansion. We also highlight the potential for these accreted clouds to create monochromatic gravitational wave signals - similar to the signals from superradiant clouds, although significantly weaker in amplitude. △ Less

Submitted 11 March, 2021; v1 submitted 16 November, 2020; originally announced November 2020.

Comments: 21 pages, 24 figures

Journal ref: Phys. Rev. D 103, 044059 (2021)

Static response and Love numbers of Schwarzschild black holes

Authors: Lam Hui, Austin Joyce, Riccardo Penco, Luca Santoni, Adam R. Solomon

Abstract: We derive the quadratic action for the physical degrees of freedom of massless spin-0, spin-1, and spin-2 perturbations on a Schwarzschild--(A)dS background in arbitrary dimensions. We then use these results to compute the static response of asymptotically flat Schwarzschild black holes to external fields. Our analysis reproduces known facts about black hole Love numbers, in particular that they v… ▽ More We derive the quadratic action for the physical degrees of freedom of massless spin-0, spin-1, and spin-2 perturbations on a Schwarzschild--(A)dS background in arbitrary dimensions. We then use these results to compute the static response of asymptotically flat Schwarzschild black holes to external fields. Our analysis reproduces known facts about black hole Love numbers, in particular that they vanish for all types of perturbation in four spacetime dimensions, but also leads to new results. For instance, we find that neutral Schwarzschild black holes polarize in the presence of an electromagnetic background in any number of spacetime dimensions except four. Moreover, we calculate for the first time black hole Love numbers for vector-type gravitational perturbations in higher dimensions and find that they generically do not vanish. Along the way, we shed some light on an apparent discrepancy between previous results in the literature, and clarify some aspects of the matching between perturbative calculations of static response on a Schwarzschild background and the point-particle effective theory △ Less

Submitted 11 February, 2024; v1 submitted 1 October, 2020; originally announced October 2020.

Comments: 78 pages, 1 figure v2: minor corrections, v3: minor corrections, v4: fixed minor error in spin-2 matching

Journal ref: JCAP 04 (2021) 052

Fast Magnetic Reconnection with Turbulence in High Lundquist Number Limit

Authors: Yang Liping, Li Hui, Guo Fan, Li Xiaocan, Li Shengtai, He jiansen, Zhang Lei, Feng Xueshang

Abstract: We use extensive 3D resistive MHD simulations to study how large-scale current sheets will undergo fast reconnection in the high Lundquist number $S$ limit (above $\sim 10^4$), when the system is subject to different externally driven turbulence levels and the self-generated turbulence produced by 3D reconnection dynamics. We find that the normalized global reconnection rate $\sim 0.01 - 0.13$, we… ▽ More We use extensive 3D resistive MHD simulations to study how large-scale current sheets will undergo fast reconnection in the high Lundquist number $S$ limit (above $\sim 10^4$), when the system is subject to different externally driven turbulence levels and the self-generated turbulence produced by 3D reconnection dynamics. We find that the normalized global reconnection rate $\sim 0.01 - 0.13$, weakly dependent on $S$. Global reconnection with the classic inflow/outflow configurations is observed, and 3D flux ropes are hierarchically formed and ejected from reconnection regions. A statistical separation of the reconnected magnetic field lines follows a super-diffusive behavior, from which the rate is measured to be very similar to that obtained from the mixing of tracer populations. We find that the reconnection rate scales roughly linearly with the turbulence level during the peak of reconnection. This scaling is consistent with the turbulence properties produced by both the externally driven and self-generation processes. These results imply that large-scale thin current sheets tend to undergo rigorous reconnection. △ Less

Submitted 14 September, 2020; originally announced September 2020.

Comments: Accepted By ApJL

Vortices and waves in light dark matter

Authors: Lam Hui, Austin Joyce, Michael J. Landry, Xinyu Li

Abstract: In a galactic halo like the Milky Way, bosonic dark matter particles lighter than about $30$ eV have a de Broglie wavelength larger than the average inter-particle separation and are therefore well described as a set of classical waves. This applies to, for instance, the QCD axion as well as to lighter axion-like particles such as fuzzy dark matter. We show that the interference of waves inside a… ▽ More In a galactic halo like the Milky Way, bosonic dark matter particles lighter than about $30$ eV have a de Broglie wavelength larger than the average inter-particle separation and are therefore well described as a set of classical waves. This applies to, for instance, the QCD axion as well as to lighter axion-like particles such as fuzzy dark matter. We show that the interference of waves inside a halo inevitably leads to vortices, locations where chance destructive interference takes the density to zero. The phase of the wavefunction has non-trivial winding around these points. This can be interpreted as a non-zero velocity circulation, so that vortices are sites where the fluid velocity has a non-vanishing curl. Using analytic arguments and numerical simulations, we study the properties of vortices and show they have a number of universal features: (1) In three spatial dimensions, the generic defects take the form of vortex rings. (2) On average there is about one vortex ring per de Broglie volume and (3) generically only single winding ($\pm 1$) vortices are found in a realistic halo. (4) The density near a vortex scales as $r^2$ while the velocity goes as $1/r$, where $r$ is the distance to vortex. (5) A vortex segment moves at a velocity inversely proportional to its curvature scale so that smaller vortex rings move faster, allowing momentary motion exceeding escape velocity. We discuss observational/experimental signatures from vortices and, more broadly, wave interference. In the ultra-light regime, gravitational lensing by interference substructures leads to flux anomalies of $5-10 \%$ in strongly lensed systems. For QCD axions, vortices lead to a diminished signal in some detection experiments but not in others. We advocate the measurement of correlation functions by axion detection experiments as a way to probe and capitalize on the expected interference substructures. △ Less

Submitted 5 November, 2020; v1 submitted 2 April, 2020; originally announced April 2020.

Comments: 61 pages, 11 figures. v2: minor edits. v3: minor changes

Quasinormal modes, echoes and the causal structure of the Green's function

Authors: Lam Hui, Daniel Kabat, Sam S. C. Wong

Abstract: Quasinormal modes describe the return to equilibrium of a perturbed system, in particular the ringdown phase of a black hole merger. But as globally-defined quantities, the quasinormal spectrum can be highly sensitive to global structure, including distant small perturbations to the potential. In what sense are quasinormal modes a property of the resulting black hole? We explore this question for… ▽ More Quasinormal modes describe the return to equilibrium of a perturbed system, in particular the ringdown phase of a black hole merger. But as globally-defined quantities, the quasinormal spectrum can be highly sensitive to global structure, including distant small perturbations to the potential. In what sense are quasinormal modes a property of the resulting black hole? We explore this question for the linearized perturbation equation with two potentials having disjoint bounded support. We give a composition law for the Wronskian that determines the quasinormal frequencies of the combined system. We show that over short time scales the evolution is governed by the quasinormal frequencies of the individual potentials, while the sensitivity to global structure can be understood in terms of echoes. We introduce an echo expansion of the Green's function and show that, as expected on general grounds, at any finite time causality limits the number of echoes that can contribute. We illustrate our results with the soluble example of a pair of $δでるた$-function potentials. We explicate the causal structure of the Green's function, demonstrating under what conditions two very different quasinormal spectra give rise to very similar ringdown waveforms. △ Less

Submitted 23 September, 2019; originally announced September 2019.

Comments: 23 pages, 5 figures

A non-perturbative test of consistency relations and their violation

Authors: Angelo Esposito, Lam Hui, Roman Scoccimarro

Abstract: In this paper, we verify the large scale structure consistency relations using N-body simulations, including modes in the highly non-linear regime. These relations (pointed out by Kehagias & Riotto and Peloso & Pietroni) follow from the symmetry of the dynamics under a shift of the Newtonian potential by a constant and a linear gradient, and predict the absence of certain poles in the ratio betwee… ▽ More In this paper, we verify the large scale structure consistency relations using N-body simulations, including modes in the highly non-linear regime. These relations (pointed out by Kehagias & Riotto and Peloso & Pietroni) follow from the symmetry of the dynamics under a shift of the Newtonian potential by a constant and a linear gradient, and predict the absence of certain poles in the ratio between the (equal time) squeezed bispectrum and power spectrum. The consistency relations, as symmetry statements, are exact, but have not been previously checked beyond the perturbative regime. Our test using N-body simulations not only offers a non-perturbative check, but also serves as a warm-up exercise for applications to observational data. A number of subtleties arise when taking the squeezed limit of the bispectrum--we show how to circumvent or address them. An interesting by-product of our investigation is an explicit demonstration that the linear-gradient symmetry is unaffected by the periodic boundary condition of the simulations. Lastly, we verify using simulations that the consistency relations are violated when the initial conditions are non-gaussian (of the local fNL type). The methodology developed here paves the way for constraining primordial non-gaussianity using large scale structure data, including (numerous) highly non-linear modes that are otherwise hard to interpret and utilize. △ Less

Submitted 29 August, 2019; v1 submitted 27 May, 2019; originally announced May 2019.

Comments: 10 pages, 5 figures, 1 table

Journal ref: Phys. Rev. D 100, 043536 (2019)

Black Hole Hair from Scalar Dark Matter

Authors: Lam Hui, Daniel Kabat, Xinyu Li, Luca Santoni, Sam S. C. Wong

Abstract: We show that a black hole surrounded by scalar dark matter develops scalar hair. This is the generalization of a phenomenon pointed out by Jacobson, that a minimally coupled scalar with a non-trivial time dependence far away from the black hole would endow the black hole with hair. In our case, the time dependence arises from the oscillation of a scalar field with a non-zero mass. We systematicall… ▽ More We show that a black hole surrounded by scalar dark matter develops scalar hair. This is the generalization of a phenomenon pointed out by Jacobson, that a minimally coupled scalar with a non-trivial time dependence far away from the black hole would endow the black hole with hair. In our case, the time dependence arises from the oscillation of a scalar field with a non-zero mass. We systematically explore the scalar profile around the black hole for different scalar masses. In the small mass limit, the scalar field has a $1/r$ component at large radius $r$, consistent with Jacobson's result. In the large mass limit (with the Compton wavelength of order of the horizon or smaller), the scalar field has a $1/r^{3/4}$ profile yielding a pile-up close to the horizon, while distinctive nodes occur for intermediate masses. Thus, the dark matter profile around a black hole, while challenging to measure, contains information about the dark matter particle mass. As an application, we consider the case of the supermassive black hole at the center of M87, recently imaged by the Event Horizon Telescope. Its horizon size is roughly the Compton wavelength of a scalar particle of mass $10^{-20}$ eV. We consider the implications of the expected scalar pile-up close to the horizon, for fuzzy dark matter at a mass of $10^{-20}$ eV or below. △ Less

Submitted 5 June, 2019; v1 submitted 29 April, 2019; originally announced April 2019.

Comments: 34 pages. v2: references added, link to animation fixed. v3: minor changes, more references, version to appear in JCAP

Inflationary soft theorems revisited: A generalized consistency relation

Authors: Lam Hui, Austin Joyce, Sam S. C. Wong

Abstract: We reconsider the derivation of soft theorems associated with nonlinearly-realized symmetries in cosmology. Utilizing the path integral, we derive a generalized consistency relation that relates a squeezed $(N+1)$-point correlation function to an $N$-point function, where the relevant soft mode is at early rather than late time. This generalized (early-late-time) version has wider applicability th… ▽ More We reconsider the derivation of soft theorems associated with nonlinearly-realized symmetries in cosmology. Utilizing the path integral, we derive a generalized consistency relation that relates a squeezed $(N+1)$-point correlation function to an $N$-point function, where the relevant soft mode is at early rather than late time. This generalized (early-late-time) version has wider applicability than the standard consistency relation where all modes are evaluated at late times. We elucidate the conditions under which the latter follows from the former. A key ingredient is the physical mode condition: that the nonlinear part of the symmetry transformation must match the time dependence of the dominant, long wavelength physical mode. This is closely related to, but distinct from, the adiabatic mode condition. Our derivation sheds light on a number of otherwise puzzling features of the standard consistency relation: (1) the underlying nonlinearly-realized symmetries (such as dilation and special conformal transformation SCT) originate as residual gauge redundancies, yet the consistency relation has physical content---for instance, it can be violated; (2) the standard consistency relation is known to fail in ultra-slow-roll inflation, but since dilation and SCT remain good symmetries, there should be a replacement for the standard relation; (3) in large scale structure applications, it is known that the standard consistency relation breaks down if the long wavelength power spectrum is too blue. The early-late-time consistency relation helps address these puzzles. We introduce a toy model where explicit checks of this generalized consistency relation are simple to carry out. Our methodology can be adapted to cases where violations of the standard consistency relation involve additional light degrees of freedom beyond the inflaton. △ Less

Submitted 10 February, 2019; v1 submitted 14 November, 2018; originally announced November 2018.

Comments: 74 pages (main text: 38 pages, appendix: 30 pages), 3 figures, 2 tables

Journal ref: JCAP 1902 (2019) 060

Effective Field Theory of Black Hole Quasinormal Modes in Scalar-Tensor Theories

Authors: Gabriele Franciolini, Lam Hui, Riccardo Penco, Luca Santoni, Enrico Trincherini

Abstract: The final ringdown phase in a coalescence process is a valuable laboratory to test General Relativity and potentially constrain additional degrees of freedom in the gravitational sector. We introduce here an effective description for perturbations around spherically symmetric spacetimes in the context of scalar-tensor theories, which we apply to study quasi-normal modes for black holes with scalar… ▽ More The final ringdown phase in a coalescence process is a valuable laboratory to test General Relativity and potentially constrain additional degrees of freedom in the gravitational sector. We introduce here an effective description for perturbations around spherically symmetric spacetimes in the context of scalar-tensor theories, which we apply to study quasi-normal modes for black holes with scalar hair. We derive the equations of motion governing the dynamics of both the polar and the axial modes in terms of the coefficients of the effective theory. Assuming the deviation of the background from Schwarzschild is small, we use the WKB method to introduce the notion of "light ring expansion". This approximation is analogous to the slow-roll expansion used for inflation, and it allows us to express the quasinormal mode spectrum in terms of a small number of parameters. This work is a first step in describing, in a model independent way, how the scalar hair can affect the ringdown stage and leave signatures on the emitted gravitational wave signal. Potential signatures include the shifting of the quasi-normal spectrum, the breaking of isospectrality between polar and axial modes, and the existence of scalar radiation. △ Less

Submitted 31 January, 2019; v1 submitted 17 October, 2018; originally announced October 2018.

Comments: 52 pages, 2 figures, 8 appendices; matches published version

Parkes Pulsar Timing Array constraints on ultralight scalar-field dark matter

Authors: Nataliya K. Porayko, Xingjiang Zhu, Yuri Levin, Lam Hui, George Hobbs, Aleksandra Grudskaya, Konstantin Postnov, Matthew Bailes, N. D. Ramesh Bhat, William Coles, Shi Dai, James Dempsey, Michael J. Keith, Matthew Kerr, Michael Kramer, Paul D. Lasky, Richard N. Manchester, Stefan Osłowski, Aditya Parthasarathy, Vikram Ravi, Daniel J. Reardon, Pablo A. Rosado, Christopher J. Russell, Ryan M. Shannon, Renée Spiewak , et al. (5 additional authors not shown)

Abstract: It is widely accepted that dark matter contributes about a quarter of the critical mass-energy density in our Universe. The nature of dark matter is currently unknown, with the mass of possible constituents spanning nearly one hundred orders of magnitude. The ultralight scalar field dark matter, consisting of extremely light bosons with $m \sim 10^{-22}$ eV and often called "fuzzy" dark matter, pr… ▽ More It is widely accepted that dark matter contributes about a quarter of the critical mass-energy density in our Universe. The nature of dark matter is currently unknown, with the mass of possible constituents spanning nearly one hundred orders of magnitude. The ultralight scalar field dark matter, consisting of extremely light bosons with $m \sim 10^{-22}$ eV and often called "fuzzy" dark matter, provides intriguing solutions to some challenges at sub-Galactic scales for the standard cold dark matter model. As shown by Khmelnitsky and Rubakov, such a scalar field in the Galaxy would produce an oscillating gravitational potential with nanohertz frequencies, resulting in periodic variations in the times of arrival of radio pulses from pulsars. The Parkes Pulsar Timing Array (PPTA) has been monitoring 20 millisecond pulsars at two to three weeks intervals for more than a decade. In addition to the detection of nanohertz gravitational waves, PPTA offers the opportunity for direct searches for fuzzy dark matter in an astrophysically feasible range of masses. We analyze the latest PPTA data set which includes timing observations for 26 pulsars made between 2004 and 2016. We perform a search in this data set for evidence of ultralight dark matter in the Galaxy using Bayesian and Frequentist methods. No statistically significant detection has been made. We therefore place upper limits on the local dark matter density. Our limits, improving on previous searches by a factor of two to five, constrain the dark matter density of ultralight bosons with $m \leq 10^{-23}$ eV to be below $6\,\text{GeV}\,\text{cm}^{-3}$ with 95\% confidence in the Earth neighborhood. Finally, we discuss the prospect of probing the astrophysically favored mass range $m \gtrsim 10^{-22}$ eV with next-generation pulsar timing facilities. △ Less

Submitted 16 December, 2019; v1 submitted 7 October, 2018; originally announced October 2018.

Comments: Fixed a typo in Fig. 5; Published in Phys. Rev. D as Editor's Suggestion article, 17 pages, 7 figures, 4 tables

Journal ref: Phys. Rev. D 98, 102002 (2018)

Numerical and Perturbative Computations of the Fuzzy Dark Matter Model

Authors: Xinyu Li, Lam Hui, Greg L. Bryan

Abstract: We investigate nonlinear structure formation in the fuzzy dark matter (FDM) model using both numerical and perturbative techniques. On the numerical side, we examine the virtues and limitations of a Schrodinger-Poisson solver (wave formulation) versus a fluid dynamics solver (Madelung formulation). We also carry out a perturbative computation of the one-loop mass power spectrum. We find that (1) i… ▽ More We investigate nonlinear structure formation in the fuzzy dark matter (FDM) model using both numerical and perturbative techniques. On the numerical side, we examine the virtues and limitations of a Schrodinger-Poisson solver (wave formulation) versus a fluid dynamics solver (Madelung formulation). We also carry out a perturbative computation of the one-loop mass power spectrum. We find that (1) in many cases, the fluid dynamics solver is capable of producing the expected interference patterns, but it fails where destructive interference causes the density to vanish which generally occurs in the nonlinear regime. (2) The Schrodinger-Poisson solver works well in all test cases, but it is demanding in resolution: one must resolve the small de Broglie scale to obtain the correct dynamics on large scales. (3) We compare the mass power spectrum from perturbation theory against that from the Schrodinger-Poisson solver, and find good agreement in the mildly nonlinear regime. Compared with fluid perturbation theory, wave perturbation theory has a more limited range of validity. (4) As an application, we compare the Lyman-alpha forest flux power spectrum obtained from the Schrodinger-Poisson solver versus one from an N-body simulation (which is often used as an approximate method to make predictions for FDM). At redshift 5, the two, starting from the same initial condition, agree to better than 10 % on observationally relevant scales as long as the FDM mass exceeds $2 \times 10^{-23}$ eV. △ Less

Submitted 3 October, 2018; originally announced October 2018.

Journal ref: Phys. Rev. D 99, 063509 (2019)

The Implications of a Pre-reionization 21 cm Absorption Signal for Fuzzy Dark Matter

Authors: Adam Lidz, Lam Hui

Abstract: The EDGES experiment recently announced evidence for a broad absorption feature in the sky-averaged radio spectrum around 78 MHz, as may result from absorption in the 21 cm line by neutral hydrogen at z~15-20. If confirmed, one implication is that the spin temperature of the 21 cm line is coupled to the gas temperature by z=20. The known mechanism for accomplishing this is the Wouthuysen-Field eff… ▽ More The EDGES experiment recently announced evidence for a broad absorption feature in the sky-averaged radio spectrum around 78 MHz, as may result from absorption in the 21 cm line by neutral hydrogen at z~15-20. If confirmed, one implication is that the spin temperature of the 21 cm line is coupled to the gas temperature by z=20. The known mechanism for accomplishing this is the Wouthuysen-Field effect, whereby Lyman-alpha photons scatter in the intergalactic medium (IGM) and impact the hyperfine level populations. This suggests that early star formation had already produced a copious Lyman-alpha background by z=20, and strongly constrains models in which the linear matter power spectrum is suppressed on small-scales, since halo and star formation are delayed in such scenarios. Here we consider the case that the dark matter consists of ultra-light axions with macroscopic de Broglie wavelengths (fuzzy dark matter, FDM). We assume that star formation tracks halo formation and adopt two simple models from the current literature for the halo mass function in FDM. We further suppose that the fraction of halo baryons which form stars is less than a conservative upper limit of $f_\star \leq 0.05$, and that ~10^4 Lyman-alpha to Lyman-limit photons are produced per stellar baryon. We find that the requirement that the 21 cm spin temperature is coupled to the gas temperature by $z=20$ places a lower-limit on the FDM particle mass of $m_a \geq 5 \times 10^{-21} {\rm eV}$. The constraint is insensitive to the precise minimum mass of halos where stars form. As the global 21 cm measurements are refined, the coupling redshift could change and we quantify how the FDM constraint would be modified. A rough translation of the FDM mass bound to a thermal relic warm dark matter (WDM) mass bound is also provided. △ Less

Submitted 3 May, 2018; originally announced May 2018.

Comments: 10 pages, 3 figures

Journal ref: Phys. Rev. D 98, 023011 (2018)

Tests of Gravity Theories Using Supermassive Black Holes

Authors: Jeremy Sakstein, Bhuvnesh Jain, Jeremy S. Heyl, Lam Hui

Abstract: Scalar-tensor theories of gravity generally violate the strong equivalence principle, namely compact objects have a suppressed coupling to the scalar force, causing them to fall slower. A black hole is the extreme example where such a coupling vanishes, i.e. black hole has no scalar hair. Following earlier work, we explore observational scenarios for detecting strong equivalence principle violatio… ▽ More Scalar-tensor theories of gravity generally violate the strong equivalence principle, namely compact objects have a suppressed coupling to the scalar force, causing them to fall slower. A black hole is the extreme example where such a coupling vanishes, i.e. black hole has no scalar hair. Following earlier work, we explore observational scenarios for detecting strong equivalence principle violation, focusing on galileon gravity as an example. For galaxies in-falling towards galaxy clusters, the supermassive black hole can be offset from the galaxy center away from the direction of the cluster. Hence, well resolved images of galaxies around nearby clusters can be used to identify the displaced black hole via the star cluster bound to it. We show that this signal is accessible with imaging surveys, both ongoing ones such as the Dark Energy Survey, and future ground and space based surveys. Already, the observation of the central black hole in M~87 places new constraints on the galileon parameters, which we present here. $\mathcal{O}(1)$ matter couplings are disfavored for a large region of the parameter space. We also find a novel phenomenon whereby the black hole can escape the galaxy completely in less than one billion years. △ Less

Submitted 24 July, 2017; v1 submitted 7 April, 2017; originally announced April 2017.

Comments: Updated to match published version

Journal ref: ApJL, 2017, Vol. 844, L14

Ultralight scalars as cosmological dark matter

Authors: Lam Hui, Jeremiah P. Ostriker, Scott Tremaine, Edward Witten

Abstract: An intriguing alternative to cold dark matter (CDM) is that the dark matter is a light ( $m \sim 10^{-22}$ eV) boson having a de Broglie wavelength $λらむだ\sim 1$ kpc, often called fuzzy dark matter (FDM). We describe the arguments from particle physics that motivate FDM, review previous work on its astrophysical signatures, and analyze several unexplored aspects of its behavior. In particular, (i) FDM… ▽ More An intriguing alternative to cold dark matter (CDM) is that the dark matter is a light ( $m \sim 10^{-22}$ eV) boson having a de Broglie wavelength $λらむだ\sim 1$ kpc, often called fuzzy dark matter (FDM). We describe the arguments from particle physics that motivate FDM, review previous work on its astrophysical signatures, and analyze several unexplored aspects of its behavior. In particular, (i) FDM halos smaller than about $10^7 (m/10^{-22} {\rm eV})^{-3/2} M_\odot$ do not form. (ii) FDM halos are comprised of a core that is a stationary, minimum-energy configuration called a "soliton", surrounded by an envelope that resembles a CDM halo. (iii) The transition between soliton and envelope is determined by a relaxation process analogous to two-body relaxation in gravitating systems, which proceeds as if the halo were composed of particles with mass $\sim ρろーλらむだ^3$ where $ρろー$ is the halo density. (iv) Relaxation may have substantial effects on the stellar disk and bulge in the inner parts of disk galaxies. (v) Relaxation can produce FDM disks but an FDM disk in the solar neighborhood must have a half-thickness of at least $300 (m/10^{-22} {\rm eV})^{-2/3}$ pc. (vi) Solitonic FDM sub-halos evaporate by tunneling through the tidal radius and this limits the minimum sub-halo mass inside 30 kpc of the Milky Way to roughly $10^8 (m/10^{-22} {\rm eV})^{-3/2} M_\odot$. (vii) If the dark matter in the Fornax dwarf galaxy is composed of CDM, most of the globular clusters observed in that galaxy should have long ago spiraled to its center, and this problem is resolved if the dark matter is FDM. △ Less

Submitted 24 February, 2017; v1 submitted 26 October, 2016; originally announced October 2016.

Comments: 37 pages, 3 figures. Additional references and discussions on clusters and the soliton-to-host-mass relation. Accepted for publication in PRD

Journal ref: Phys. Rev. D 95, 043541 (2017)

Measurement of the dipole in the cross-correlation function of galaxies

Authors: Enrique Gaztanaga, Camille Bonvin, Lam Hui

Abstract: It is usually assumed that in the linear regime the two-point correlation function of galaxies contains only a monopole, quadrupole and hexadecapole. Looking at cross-correlations between different populations of galaxies, this turns out not to be the case. In particular, the cross-correlations between a bright and a faint population of galaxies contain also a dipole. In this paper we present the… ▽ More It is usually assumed that in the linear regime the two-point correlation function of galaxies contains only a monopole, quadrupole and hexadecapole. Looking at cross-correlations between different populations of galaxies, this turns out not to be the case. In particular, the cross-correlations between a bright and a faint population of galaxies contain also a dipole. In this paper we present the first attempt to measure this dipole. We discuss the four types of effects that contribute to the dipole: relativistic distortions, evolution effect, wide-angle effect and large-angle effect. We show that the first three contributions are intrinsic anti-symmetric contributions that do not depend on the choice of angle used to measure the dipole. On the other hand the large-angle effect appears only if the angle chosen to extract the dipole breaks the symmetry of the problem. We show that the relativistic distortions, the evolution effect and the wide-angle effect are too small to be detected in the LOWz and CMASS sample of the BOSS survey. On the other hand with a specific combination of angles we are able to measure the large-angle effect with high significance. We emphasise that this large-angle dipole does not contain new physical information, since it is just a geometrical combination of the monopole and the quadrupole. However this measurement, which is in excellent agreement with theoretical predictions, validates our method for extracting the dipole from the two-point correlation function and it opens the way to the detection of relativistic effects in future surveys like e.g. DESI. △ Less

Submitted 16 January, 2017; v1 submitted 12 December, 2015; originally announced December 2015.

Comments: 15 pages, 17 figures. v2: 20 pages, 17 figures. Section IIIc partly rewritten, new section IIId, new figures 16 and 17. Main results unchanged. Matches published version in JCAP

Journal ref: JCAP01(2017)032

Optimising the measurement of relativistic distortions in large-scale structure

Authors: Camille Bonvin, Lam Hui, Enrique Gaztanaga

Abstract: It has been shown recently that relativistic distortions generate a dipolar modulation in the two-point correlation function of galaxies. To measure this relativistic dipole it is necessary to cross-correlate different populations of galaxies with for example different luminosities or colours. In this paper, we construct an optimal estimator to measure the dipole with multiple populations. We show… ▽ More It has been shown recently that relativistic distortions generate a dipolar modulation in the two-point correlation function of galaxies. To measure this relativistic dipole it is necessary to cross-correlate different populations of galaxies with for example different luminosities or colours. In this paper, we construct an optimal estimator to measure the dipole with multiple populations. We show that this estimator increases the signal-to-noise of the dipole by up to 35 percent. Using 6 populations of galaxies, in a survey with halos and number densities similar to those of the millennium simulation, we forecast a cumulative signal-to-noise of 4.4. For the main galaxy sample of SDSS at low redshift z<0.2 our optimal estimator predicts a cumulative signal-to-noise of 2.4. Finally we forecast a cumulative signal-to-noise of 7.4 in the upcoming DESI survey. These forecasts indicate that with the appropriate choice of estimator the relativistic dipole should be detectable in current and future surveys. △ Less

Submitted 16 January, 2017; v1 submitted 11 December, 2015; originally announced December 2015.

Comments: 20 pages, 8 figures. v2: 19 pages, 6 figures. Explicit derivation of the cancellation of the cosmic variance added. Some references added, minor changes to the text, main results unchanged. Matches version published in JCAP

Journal ref: JCAP08(2016)021

Eötvös Experiments with Supermassive Black Holes

Authors: Asha Asvathaman, Jeremy S. Heyl, Lam Hui

Abstract: By examining the locations of central black holes in two elliptical galaxies, M\,32 and M\,87, we derive constraints on the violation of the strong equivalence principle for purely gravitational objects, i.e. black holes, of less than about two-thirds, $ηいーた_N<0.68$ from the gravitational interaction of M\,87 with its neighbours in the Virgo cluster. Although M\,32 appears to be a good candidate for… ▽ More By examining the locations of central black holes in two elliptical galaxies, M\,32 and M\,87, we derive constraints on the violation of the strong equivalence principle for purely gravitational objects, i.e. black holes, of less than about two-thirds, $ηいーた_N<0.68$ from the gravitational interaction of M\,87 with its neighbours in the Virgo cluster. Although M\,32 appears to be a good candidate for this technique, the high concentration of stars near its centre substantially weakens the constraints. On the other hand, if a central black hole is found in NGC 205 or one of the other satellite ellipticals of M\,31, substantially better constraints could be obtained. In all cases the constraints could improve dramatically with better astrometry. △ Less

Submitted 13 December, 2016; v1 submitted 25 June, 2015; originally announced June 2015.

Comments: 7 pages, 3 figures, accepted by Monthly Notices

Phenomenology of dark energy: general features of large-scale perturbations

Authors: Louis Perenon, Federico Piazza, Christian Marinoni, Lam Hui

Abstract: We present a systematic exploration of dark energy and modified gravity models containing a single scalar field non-minimally coupled to the metric. Even though the parameter space is large, by exploiting an effective field theory (EFT) formulation and by imposing simple physical constraints such as stability conditions and (sub-)luminal propagation of perturbations, we arrive at a number of gener… ▽ More We present a systematic exploration of dark energy and modified gravity models containing a single scalar field non-minimally coupled to the metric. Even though the parameter space is large, by exploiting an effective field theory (EFT) formulation and by imposing simple physical constraints such as stability conditions and (sub-)luminal propagation of perturbations, we arrive at a number of generic predictions. (1) The linear growth rate of matter density fluctuations is generally suppressed compared to $Λらむだ$CDM at intermediate redshifts ($0.5 \lesssim z \lesssim 1$), despite the introduction of an attractive long-range scalar force. This is due to the fact that, in self-accelerating models, the background gravitational coupling weakens at intermediate redshifts, over-compensating the effect of the attractive scalar force. (2) At higher redshifts, the opposite happens; we identify a period of super-growth when the linear growth rate is larger than that predicted by $Λらむだ$CDM. (3) The gravitational slip parameter $ηいーた$ - the ratio of the space part of the metric perturbation to the time part - is bounded from above. For Brans-Dicke-type theories $ηいーた$ is at most unity. For more general theories, $ηいーた$ can exceed unity at intermediate redshifts, but not more than about $1.5$ if, at the same time, the linear growth rate is to be compatible with current observational constraints. We caution against phenomenological parametrization of data that do not correspond to predictions from viable physical theories. We advocate the EFT approach as a way to constrain new physics from future large-scale-structure data. △ Less

Submitted 10 July, 2015; v1 submitted 9 June, 2015; originally announced June 2015.

Comments: 24 pages, 7 figures

Journal ref: JCAP 1511 (2015) no.11, 029

Emulating the CFHTLenS Weak Lensing data: Cosmological Constraints from moments and Minkowski functionals

Authors: Andrea Petri, Jia Liu, Zoltan Haiman, Morgan May, Lam Hui, Jan M. Kratochvil

Abstract: Weak gravitational lensing is a powerful cosmological probe, with non--Gaussian features potentially containing the majority of the information. We examine constraints on the parameter triplet $(Ωおめが_m,w,σしぐま_8)$ from non-Gaussian features of the weak lensing convergence field, including a set of moments (up to $4^{\rm th}$ order) and Minkowski functionals, using publicly available data from the 154deg… ▽ More Weak gravitational lensing is a powerful cosmological probe, with non--Gaussian features potentially containing the majority of the information. We examine constraints on the parameter triplet $(Ωおめが_m,w,σしぐま_8)$ from non-Gaussian features of the weak lensing convergence field, including a set of moments (up to $4^{\rm th}$ order) and Minkowski functionals, using publicly available data from the 154deg$^2$ CFHTLenS survey. We utilize a suite of ray--tracing N-body simulations spanning 91 points in $(Ωおめが_m,w,σしぐま_8)$ parameter space, replicating the galaxy sky positions, redshifts and shape noise in the CFHTLenS catalogs. We then build an emulator that interpolates the simulated descriptors as a function of $(Ωおめが_m,w,σしぐま_8)$, and use it to compute the likelihood function and parameter constraints. We employ a principal component analysis to reduce dimensionality and to help stabilize the constraints with respect to the number of bins used to construct each statistic. Using the full set of statistics, we find $Σしぐま_8\equivσしぐま_8(Ωおめが_m/0.27)^{0.55}=0.75\pm0.04$ (68% C.L.), in agreement with previous values. We find that constraints on the $(Ωおめが_m,σしぐま_8)$ doublet from the Minkowski functionals suffer a strong bias. However, high-order moments break the $(Ωおめが_m,σしぐま_8)$ degeneracy and provide a tight constraint on these parameters with no apparent bias. The main contribution comes from quartic moments of derivatives. △ Less

Submitted 18 May, 2015; v1 submitted 20 March, 2015; originally announced March 2015.

Comments: 17 pages, 10 figures, 4 tables

Journal ref: Phys. Rev. D, 91, 103511 (2015)

Lagrangian space consistency relation for large scale structure

Authors: Bart Horn, Lam Hui, Xiao Xiao

Abstract: Consistency relations, which relate the squeezed limit of an (N+1)-point correlation function to an N-point function, are non-perturbative symmetry statements that hold even if the associated high momentum modes are deep in the nonlinear regime and astrophysically complex. Recently, Kehagias & Riotto and Peloso & Pietroni discovered a consistency relation applicable to large scale structure. We sh… ▽ More Consistency relations, which relate the squeezed limit of an (N+1)-point correlation function to an N-point function, are non-perturbative symmetry statements that hold even if the associated high momentum modes are deep in the nonlinear regime and astrophysically complex. Recently, Kehagias & Riotto and Peloso & Pietroni discovered a consistency relation applicable to large scale structure. We show that this can be recast into a simple physical statement in Lagrangian space: that the squeezed correlation function (suitably normalized) vanishes. This holds regardless of whether the correlation observables are at the same time or not, and regardless of whether multiple-streaming is present. The simplicity of this statement suggests that an analytic understanding of large scale structure in the nonlinear regime may be particularly promising in Lagrangian space. △ Less

Submitted 24 February, 2015; originally announced February 2015.

Comments: 19 pages, no figures

Cosmology Constraints from the Weak Lensing Peak Counts and the Power Spectrum in CFHTLenS

Authors: Jia Liu, Andrea Petri, Zoltan Haiman, Lam Hui, Jan M. Kratochvil, Morgan May

Abstract: Lensing peaks have been proposed as a useful statistic, containing cosmological information from non-Gaussianities that is inaccessible from traditional two-point statistics such as the power spectrum or two-point correlation functions. Here we examine constraints on cosmological parameters from weak lensing peak counts, using the publicly available data from the 154 deg$^2$ CFHTLenS survey. We ut… ▽ More Lensing peaks have been proposed as a useful statistic, containing cosmological information from non-Gaussianities that is inaccessible from traditional two-point statistics such as the power spectrum or two-point correlation functions. Here we examine constraints on cosmological parameters from weak lensing peak counts, using the publicly available data from the 154 deg$^2$ CFHTLenS survey. We utilize a new suite of ray-tracing N-body simulations on a grid of 91 cosmological models, covering broad ranges of the three parameters $Ωおめが_m$, $σしぐま_8$, and $w$, and replicating the Galaxy sky positions, redshifts, and shape noise in the CFHTLenS observations. We then build an emulator that interpolates the power spectrum and the peak counts to an accuracy of $\leq 5\%$, and compute the likelihood in the three-dimensional parameter space ($Ωおめが_m$, $σしぐま_8$, $w$) from both observables. We find that constraints from peak counts are comparable to those from the power spectrum, and somewhat tighter when different smoothing scales are combined. Neither observable can constrain $w$ without external data. When the power spectrum and peak counts are combined, the area of the error "banana'' in the ($Ωおめが_m$, $σしぐま_8$) plane reduces by a factor of $\approx2$, compared to using the power spectrum alone. For a flat $Λらむだ$ cold dark matter model, combining both statistics, we obtain the constraint $σしぐま_8(Ωおめが_m/0.27)^{0.63}=0.85\substack{+0.03 \\ -0.03}$. △ Less

Submitted 1 October, 2015; v1 submitted 1 December, 2014; originally announced December 2014.

Comments: 15 pages, 14 figures. The manuscript has been revised and matches published version

Journal ref: Phys. Rev. D 91, 063507 (2015)

Soft-Pion Theorems for Large Scale Structure

Authors: Bart Horn, Lam Hui, Xiao Xiao

Abstract: Consistency relations -- which relate an N-point function to a squeezed (N+1)-point function -- are useful in large scale structure (LSS) because of their non-perturbative nature: they hold even if the N-point function is deep in the nonlinear regime, and even if they involve astrophysically messy galaxy observables. The non-perturbative nature of the consistency relations is guaranteed by the fac… ▽ More Consistency relations -- which relate an N-point function to a squeezed (N+1)-point function -- are useful in large scale structure (LSS) because of their non-perturbative nature: they hold even if the N-point function is deep in the nonlinear regime, and even if they involve astrophysically messy galaxy observables. The non-perturbative nature of the consistency relations is guaranteed by the fact that they are symmetry statements, in which the velocity plays the role of the soft pion. In this paper, we address two issues: (1) how to derive the relations systematically using the residual coordinate freedom in the Newtonian gauge, and relate them to known results in $ζぜーた$-gauge (often used in studies of inflation); (2) under what conditions the consistency relations are violated. In the non-relativistic limit, our derivation reproduces the Newtonian consistency relation discovered by Kehagias \& Riotto and Peloso & Pietroni. More generally, there is an infinite set of consistency relations, as is known in $ζぜーた$-gauge. There is a one-to-one correspondence between symmetries in the two gauges; in particular, the Newtonian consistency relation follows from the dilation and special conformal symmetries in $ζぜーた$-gauge. We probe the robustness of the consistency relations by studying models of galaxy dynamics and biasing. We give a systematic list of conditions under which the consistency relations are violated; violations occur if the galaxy bias is non-local in an infrared divergent way. We emphasize the relevance of the adiabatic mode condition, as distinct from symmetry considerations. As a by-product of our investigation, we discuss a simple fluid Lagrangian for LSS. △ Less

Submitted 3 June, 2014; originally announced June 2014.

Comments: 48 pages, no figures

Two-dimensional Lorentz invariance of spherically symmetric black holes

Authors: Lam Hui, Alberto Nicolis

Abstract: We show that a static spherically symmetric black hole, in a generic theory of gravity with generic matter fields, has a two-dimensional Lorentz symmetry. We show that a static spherically symmetric black hole, in a generic theory of gravity with generic matter fields, has a two-dimensional Lorentz symmetry. △ Less

Submitted 26 February, 2014; originally announced February 2014.

Comments: To appear in Phys. Rev. D

Single-Field Consistency Relations of Large Scale Structure. Part III: Test of the Equivalence Principle

Authors: Paolo Creminelli, Jérôme Gleyzes, Lam Hui, Marko Simonović, Filippo Vernizzi

Abstract: The recently derived consistency relations for Large Scale Structure do not hold if the Equivalence Principle (EP) is violated. We show it explicitly in a toy model with two fluids, one of which is coupled to a fifth force. We explore the constraints that galaxy surveys can set on EP violation looking at the squeezed limit of the 3-point function involving two populations of objects. We find that… ▽ More The recently derived consistency relations for Large Scale Structure do not hold if the Equivalence Principle (EP) is violated. We show it explicitly in a toy model with two fluids, one of which is coupled to a fifth force. We explore the constraints that galaxy surveys can set on EP violation looking at the squeezed limit of the 3-point function involving two populations of objects. We find that one can explore EP violations of order 10^{-3} - 10^{-4} on cosmological scales. Chameleon models are already very constrained by the requirement of screening within the Solar System and only a very tiny region of the parameter space can be explored with this method. We show that no violation of the consistency relations is expected in Galileon models. △ Less

Submitted 15 July, 2014; v1 submitted 20 December, 2013; originally announced December 2013.

Comments: 16 pages, 1 figure. v2: minor changes, JCAP published version

The Impact of Magnification and Size Bias on Weak Lensing Power Spectrum and Peak Statistics

Authors: Jia Liu, Zoltán Haiman, Lam Hui, Jan M. Kratochvil, Morgan May

Abstract: The weak lensing power spectrum is a powerful tool to probe cosmological parameters. Additionally, lensing peak counts contain cosmological information beyond the power spectrum. Both of these statistics can be affected by the preferential selection of source galaxies in patches of the sky with high magnification, as well as by the dilution in the source galaxy surface density in such regions. If… ▽ More The weak lensing power spectrum is a powerful tool to probe cosmological parameters. Additionally, lensing peak counts contain cosmological information beyond the power spectrum. Both of these statistics can be affected by the preferential selection of source galaxies in patches of the sky with high magnification, as well as by the dilution in the source galaxy surface density in such regions. If not accounted for, these biases introduce systematic errors for cosmological measurements. Here we quantify these systematic errors, using convergence maps from a suite of ray-tracing N-body simulations. At the cut-off magnitude m of on-going and planned major weak lensing surveys, the logarithmic slope of the cumulative number counts s = dlog[n(>m)]/dlog(m) is in the range 0.1 < s < 0.5. At s = 0.2, expected in the I band for LSST, the inferred values of Omega_m, w and sigma_8 are biased by many sigma (where sigma denotes the marginalized error) and therefore the biases will need to be carefully modeled. We also find that the parameters are biased differently in the (Omega_m, w, sigma_8) parameter space when the power spectrum and when the peak counts are used. In particular, w derived from the power spectrum is less affected than w derived from peak counts, while the opposite is true for the best-constrained combination of [sigma_8 Omega_m^gamma] (with gamma=0.62 from the power spectrum and gamma = 0.48 from peak counts). This suggests that the combination of the power spectrum and peak counts can help mitigate the impact of magnification and size biases. △ Less

Submitted 28 October, 2013; originally announced October 2013.

Comments: 14 pages, 11 figures

Journal ref: Phys. Rev. D 89, 023515 (2014), 14pp

Novel Probes of Gravity and Dark Energy

Authors: Bhuvnesh Jain, Austin Joyce, Rodger Thompson, Amol Upadhye, James Battat, Philippe Brax, Anne-Christine Davis, Claudia de Rham, Scott Dodelson, Adrienne Erickcek, Gregory Gabadadze, Wayne Hu, Lam Hui, Dragan Huterer, Marc Kamionkowski, Justin Khoury, Kazuya Koyama, Baojiu Li, Eric Linder, Fabian Schmidt, Roman Scoccimarro, Glenn Starkman, Chris Stubbs, Masahiro Takada, Andrew Tolley , et al. (7 additional authors not shown)

Abstract: The discovery of cosmic acceleration has stimulated theorists to consider dark energy or modifications to Einstein's General Relativity as possible explanations. The last decade has seen advances in theories that go beyond smooth dark energy -- modified gravity and interactions of dark energy. While the theoretical terrain is being actively explored, the generic presence of fifth forces and dark s… ▽ More The discovery of cosmic acceleration has stimulated theorists to consider dark energy or modifications to Einstein's General Relativity as possible explanations. The last decade has seen advances in theories that go beyond smooth dark energy -- modified gravity and interactions of dark energy. While the theoretical terrain is being actively explored, the generic presence of fifth forces and dark sector couplings suggests a set of distinct observational signatures. This report focuses on observations that differ from the conventional probes that map the expansion history or large-scale structure. Examples of such novel probes are: detection of scalar fields via lab experiments, tests of modified gravity using stars and galaxies in the nearby universe, comparison of lensing and dynamical masses of galaxies and clusters, and the measurements of fundamental constants at high redshift. The observational expertise involved is very broad as it spans laboratory experiments, high resolution astronomical imaging and spectroscopy and radio observations. In the coming decade, searches for these effects have the potential for discovering fundamental new physics. We discuss how the searches can be carried out using experiments that are already under way or with modest adaptations of existing telescopes or planned experiments. The accompanying paper on the Growth of Cosmic Structure describes complementary tests of gravity with observations of large-scale structure. △ Less

Submitted 24 September, 2013; v1 submitted 20 September, 2013; originally announced September 2013.

Comments: Report from the "Dark Energy and CMB" working group for the American Physical Society's Division of Particles and Fields long-term planning exercise ("Snowmass")

Cosmology with Minkowski functionals and moments of the weak lensing convergence field

Authors: Andrea Petri, Zoltán Haiman, Lam Hui, Morgan May, Jan M. Kratochvil

Abstract: We compare the efficiency of moments and Minkowski functionals (MFs) in constraining the subset of cosmological parameters (Omega_m,w,sigma_8) using simulated weak lensing convergence maps. We study an analytic perturbative expansion of the MFs in terms of the moments of the convergence field and of its spatial derivatives. We show that this perturbation series breaks down on smoothing scales belo… ▽ More We compare the efficiency of moments and Minkowski functionals (MFs) in constraining the subset of cosmological parameters (Omega_m,w,sigma_8) using simulated weak lensing convergence maps. We study an analytic perturbative expansion of the MFs in terms of the moments of the convergence field and of its spatial derivatives. We show that this perturbation series breaks down on smoothing scales below 5', while it shows a good degree of convergence on larger scales (15'). Most of the cosmological distinguishing power is lost when the maps are smoothed on these larger scales. We also show that, on scales comparable to 1', where the perturbation series does not converge, cosmological constraints obtained from the MFs are approximately 1.5-2 times better than the ones obtained from the first few moments of the convergence distribution --- provided that the latter include spatial information, either from moments of gradients, or by combining multiple smoothing scales. Including either a set of these moments or the MFs can significantly tighten constraints on cosmological parameters, compared to the conventional method of using the power spectrum alone. △ Less

Submitted 5 January, 2014; v1 submitted 17 September, 2013; originally announced September 2013.

Comments: 16 pages, 9 figures, 6 tables

Journal ref: Phys. Rev. D 88, 123002 (2013)

Asymmetric galaxy correlation functions

Authors: Camille Bonvin, Lam Hui, Enrique Gaztanaga

Abstract: We study the two-point cross-correlation function between two populations of galaxies: for instance a bright population and a faint population. We show that this cross-correlation is asymmetric under the exchange of the line-of-sight coordinate of the galaxies, i.e. that the correlation is different if the bright galaxy is in front of, or behind, the faint galaxy. We give an intuitive, quasi-Newto… ▽ More We study the two-point cross-correlation function between two populations of galaxies: for instance a bright population and a faint population. We show that this cross-correlation is asymmetric under the exchange of the line-of-sight coordinate of the galaxies, i.e. that the correlation is different if the bright galaxy is in front of, or behind, the faint galaxy. We give an intuitive, quasi-Newtonian derivation of all the effects that contribute to such an asymmetry in large-scale structure: gravitational redshift, Doppler shift, lensing, light-cone, evolution and Alcock-Paczynski effects - interestingly, the gravitational redshift term is exactly canceled by some of the others, assuming geodesic motion. Most of these effects are captured by previous calculations of general relativistic corrections to the observed galaxy density fluctuation; the asymmetry arises from terms that are suppressed by the ratio H/k - H is the Hubble constant and k is the wavenumber - which are more readily observable than the terms suppressed by (H/k)^2. Some of the contributions to the asymmetry, however, arise from terms that are generally considered 'Newtonian' - the lensing and evolution - and thus represent a contaminant in the search for general relativistic corrections. We propose methods to disentangle these different contributions. A simple method reduces the contamination to a level of < 10% for redshifts z<1. We also clarify the relation to recent work on measuring gravitational redshifts by stacking clusters. △ Less

Submitted 29 April, 2014; v1 submitted 5 September, 2013; originally announced September 2013.

Comments: v1: 22 pages, 12 figures; v2: discussion about the Alcock-Paczynski effect and about dust extinction added, new figure 13 explaining the sign of the dipole, discussion in appendix E about the comparison with clusters measurement extended. This version matches the published version

Journal ref: Phys. Rev. D 89, 083535 (2014)

Dynamics of Strongly Twisted Relativistic Magnetospheres

Authors: Kyle Parfrey, Andrei M. Beloborodov, Lam Hui

Abstract: Magnetar magnetospheres are believed to be strongly twisted, due to shearing of the stellar crust by internal magnetic stresses. We present time-dependent axisymmetric simulations showing in detail the evolution of relativistic force-free magnetospheres subjected to slow twisting through large angles. When the twist amplitude is small, the magnetosphere moves quasi-statically through a sequence of… ▽ More Magnetar magnetospheres are believed to be strongly twisted, due to shearing of the stellar crust by internal magnetic stresses. We present time-dependent axisymmetric simulations showing in detail the evolution of relativistic force-free magnetospheres subjected to slow twisting through large angles. When the twist amplitude is small, the magnetosphere moves quasi-statically through a sequence of equilibria of increasing free energy. At some twist amplitude the magnetosphere becomes tearing-mode unstable to forming a resistive current sheet, initiating large-scale magnetic reconnection in which a significant fraction of the magnetic free energy can be dissipated. This "critical" twist angle is insensitive to the resistive length scale. Rapid shearing temporarily stabilizes the magnetosphere beyond the critical angle, allowing the magnetosphere of a rapidly differentially rotating star to store, and dissipate, more free energy. In addition to these effects, shearing the surface of a rotating star increases the spindown torque applied to the star. If shearing is much slower than rotation, the resulting spikes in spindown rate can occur on timescales anywhere from the long twisting timescale to the stellar spin period or shorter, depending both on the stellar shear distribution and the existing distribution of magnetospheric twists. A model in which energy is stored in the magnetosphere, and released by a magnetospheric instability, therefore predicts large changes in the measured spindown rate before SGR giant flares. △ Less

Submitted 11 December, 2013; v1 submitted 18 June, 2013; originally announced June 2013.

Comments: 28 pages, 23 figures, published in ApJ

An Infinite Set of Ward Identities for Adiabatic Modes in Cosmology

Authors: Kurt Hinterbichler, Lam Hui, Justin Khoury

Abstract: We show that the correlation functions of any single-field cosmological model with constant growing-modes are constrained by an infinite number of novel consistency relations, which relate (N+1)-point correlation functions with a soft-momentum scalar or tensor mode to a symmetry transformation on N-point correlation functions of hard-momentum modes. We derive these consistency relations from Ward… ▽ More We show that the correlation functions of any single-field cosmological model with constant growing-modes are constrained by an infinite number of novel consistency relations, which relate (N+1)-point correlation functions with a soft-momentum scalar or tensor mode to a symmetry transformation on N-point correlation functions of hard-momentum modes. We derive these consistency relations from Ward identities for an infinite tower of non-linearly realized global symmetries governing scalar and tensor perturbations. These symmetries can be labeled by an integer n. At each order n, the consistency relations constrain - completely for n=0,1, and partially for n>= 2 - the q^n behavior of the soft limits. The identities at n=0 recover Maldacena's original consistency relations for a soft scalar and tensor mode, n=1 gives the recently-discovered conformal consistency relations, and the identities for n>= 2 are new. As a check, we verify directly that the n=2 identity is satisfied by known correlation functions in slow-roll inflation. △ Less

Submitted 27 January, 2014; v1 submitted 19 April, 2013; originally announced April 2013.

Comments: 47 pages. v3 typos corrected, matches published version

Journal ref: JCAP 01 (2014) 039

On Cosmological Perturbations of Quasidilaton

Authors: Guido D'Amico, Gregory Gabadadze, Lam Hui, David Pirtskhalava

Abstract: A theory of the quasidilaton is an extension of massive gravity by a scalar field, nonlinearly realizing a certain new global symmetry of the Lagrangian. It has been shown that unlike pure massive gravity, this theory does admit homogeneous and isotropic spatially flat solutions. Among the latter, selfaccelerated solutions attract a special attention. Previous studies of perturbations, performed i… ▽ More A theory of the quasidilaton is an extension of massive gravity by a scalar field, nonlinearly realizing a certain new global symmetry of the Lagrangian. It has been shown that unlike pure massive gravity, this theory does admit homogeneous and isotropic spatially flat solutions. Among the latter, selfaccelerated solutions attract a special attention. Previous studies of perturbations, performed in the decoupling limit, revealed one healthy scalar mode, while the second relevant scalar was not captured in that limit. Here we study full cosmological perturbations above the simplest selfaccelerated background. We show that the fluctuations of a mixed state of the quasidilaton and the helicity-0 graviton necessarily have a negative kinetic term at short distances, making this background unphysical. In addition, these cosmologies exhibit an order one sensitivity to higher dimensional terms suppressed by an energy scale that is parametrically higher than the strong coupling scale of the quasidilaton effective theory: such terms include Galileons, Goldstone-like selfinteractions and derivatives of the quasidilaton coupled to curvature, none of which introduce extra Ostrogradsky states. As one consequence, cosmology at the Hubble distances for this particular class of solutions depends on an unknown extension of the quasidilaton below its strong coupling distance scale. We note that non-FRW solutions that are similar to those of pure massive gravity should not necessarily suffer from these pathologies. △ Less

Submitted 2 April, 2013; originally announced April 2013.

Comments: 15 pages

Report number: UCSD/PTH 13-03

One-particle-irreducible consistency relations for cosmological perturbations

Authors: Walter D. Goldberger, Lam Hui, Alberto Nicolis

Abstract: We derive consistency relations for correlators of scalar cosmological perturbations which hold in the "squeezed limit" in which one or more of the external momenta become soft. Our results are formulated as relations between suitably defined one-particle irreducible N-point and (N-1)-point functions that follow from residual spatial conformal diffeomorphisms of the unitary gauge Lagrangian. As su… ▽ More We derive consistency relations for correlators of scalar cosmological perturbations which hold in the "squeezed limit" in which one or more of the external momenta become soft. Our results are formulated as relations between suitably defined one-particle irreducible N-point and (N-1)-point functions that follow from residual spatial conformal diffeomorphisms of the unitary gauge Lagrangian. As such, some of these relations are exact to all orders in perturbation theory, and do not rely on approximate deSitter invariance or other dynamical assumptions (e.g., properties of the operator product expansion or the behavior of modes at horizon crossing). The consistency relations apply model-independently to cosmological scenarios where the time evolution is driven by a single scalar field. Besides reproducing the known results for single-field inflation in the slow roll limit, we verify that our consistency relations hold more generally, for instance in ghost condensate models in flat space. We comment on possible extensions of our results to multi-field models. △ Less

Submitted 5 March, 2013; originally announced March 2013.

Comments: 12 pages

Binary Systems as Resonance Detectors for Gravitational Waves

Authors: Lam Hui, Sean T. McWilliams, I-Sheng Yang

Abstract: Gravitational waves at suitable frequencies can resonantly interact with a binary system, inducing changes to its orbit. A stochastic gravitational-wave background causes the orbital elements of the binary to execute a classic random walk, with the variance of orbital elements growing with time. The lack of such a random walk in binaries that have been monitored with high precision over long time-… ▽ More Gravitational waves at suitable frequencies can resonantly interact with a binary system, inducing changes to its orbit. A stochastic gravitational-wave background causes the orbital elements of the binary to execute a classic random walk, with the variance of orbital elements growing with time. The lack of such a random walk in binaries that have been monitored with high precision over long time-scales can thus be used to place an upper bound on the gravitational-wave background. Using periastron time data from the Hulse-Taylor binary pulsar spanning ~30 years, we obtain a bound of h_c < 7.9*10^(-14) at ~10^(-4) Hz, where h_c is the strain amplitude per logarithmic frequency interval. Our constraint complements those from pulsar timing arrays, which probe much lower frequencies, and ground-based gravitational-wave observations, which probe much higher frequencies. Interesting sources in our frequency band, which overlaps the lower sensitive frequencies of proposed space-based observatories, include white-dwarf/supermassive black-hole binaries in the early/late stages of inspiral, and TeV scale preheating or phase transitions. The bound improves as (time span)^(-2) and (sampling rate)^(-1/2). The Hulse-Taylor constraint can be improved to ~3.8*10^(-15) with a suitable observational campaign over the next decade. Our approach can also be applied to other binaries, including (with suitable care) the Earth-Moon system, to obtain constraints at different frequencies. The observation of additional binary pulsars with the SKA could reach a sensitivity of h_c ~ 3*10^(-17). △ Less

Submitted 1 March, 2013; v1 submitted 11 December, 2012; originally announced December 2012.

Comments: 8 pages, 1 figure. version 2, clarified several points thanks to private communications with Marc Kamionkowski