Source author record

Lam Hui

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

51works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

51 published item(s)

preprint2026arXiv

Ladder Symmetries of Higher Dimensional Black Holes

We compute the ladder operators for static tidal perturbations to higher-dimensional black holes. These operators map between solutions of the relevant equation of motion at different multipole orders. We focus on spin 0, 1 and 2 perturbations to the Schwarzschild-Tangherlini black hole and on spin 0 perturbations to the 5D Myers-Perry black hole. The ladder structure, used in conjunction with the existence of special ground state solutions, explains why the Love numbers of these higher-dimensional black holes vanish for specific combinations of the multipole moment and number of spacetime dimensions. This generalizes previous work on a ladder symmetry explanation for the vanishing of 4D black hole static Love numbers to higher dimensions.

preprint2023arXiv

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

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.

preprint2023arXiv

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

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.

preprint2022arXiv

Effective Field Theory for the Perturbations of a Slowly Rotating Black Hole

We develop the effective theory for perturbations around black holes with scalar hair, in two directions. First, we show that the scalar-Gauss--Bonnet theory, often used as an example exhibiting scalar black hole hair, can be deformed by galileon operators leading to order unity changes to its predictions. The effective theory for perturbations thus provides an efficient framework for describing and constraining broad classes of scalar-tensor theories, of which the addition of galileon operators is an example. Second, we extend the effective theory to perturbations around an axisymmetric, slowly rotating black hole, at linear order in the black hole spin. We also discuss the inclusion of parity-breaking operators in the effective theory.

preprint2020arXiv

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

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.

preprint2020arXiv

Vortices and waves in light dark matter

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.

preprint2016arXiv

Eötvös Experiments with Supermassive Black Holes

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.

preprint2015arXiv

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

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}$.

preprint2015arXiv

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

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.

preprint2015arXiv

Lagrangian space consistency relation for large scale structure

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.

preprint2014arXiv

An Infinite Set of Ward Identities for Adiabatic Modes in Cosmology

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.

preprint2014arXiv

Asymmetric galaxy correlation functions

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.

preprint2014arXiv

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

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.

preprint2014arXiv

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

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.

preprint2014arXiv

Soft-Pion Theorems for Large Scale Structure

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.

preprint2013arXiv

Binary Systems as Resonance Detectors for Gravitational Waves

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).

preprint2013arXiv

Complementarity + Back-reaction is enough

We investigate a recent development of the black hole information problem, in which a practical paradox has been formulated to show that complementarity is insufficient. A crucial ingredient in this practical paradox is to distill information from the early Hawking radiation within the past lightcone of the black hole. By causality this action can back-react on the black hole. Taking this back-reaction into account, the paradox could be resolved without invoking any new physics beyond complementarity. This resolution requires a certain constraint on the S-matrix to be satisfied. Further insights into the S-matrix could potentially be obtained by effective-field-theory computations of the back-reaction on the nice slice.

preprint2013arXiv

Dynamics of Strongly Twisted Relativistic Magnetospheres

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.

preprint2013arXiv

Novel Probes of Gravity and Dark Energy

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.

preprint2013arXiv

On Cosmological Perturbations of Quasidilaton

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.

preprint2013arXiv

One-particle-irreducible consistency relations for cosmological perturbations

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.

preprint2013arXiv

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

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.

preprint2012arXiv

An observational test of the Vainshtein mechanism

Modified gravity theories capable of genuine self-acceleration typically invoke a galileon scalar which mediates a long range force, but is screened by the Vainshtein mechanism on small scales. In such theories, non-relativistic stars carry the full scalar charge (proportional to their mass), while black holes carry none. Thus, for a galaxy free-falling in some external gravitational field, its central massive black hole is expected to lag behind the stars. To look for this effect, and to distinguish it from other astrophysical effects, one can correlate the gravitational pull from the surrounding structure with the offset between the stellar center and the black hole. The expected offset depends on the central density of the galaxy, and ranges up to ~0.1 kpc for small galaxies. The observed offset in M87 cannot be explained by this effect unless the scalar force is significantly stronger than gravity. We also discuss the systematic offset of compact objects from the galactic plane as another possible signature.

preprint2012arXiv

CMB Power Asymmetry from Non-Gaussian Modulation

Non-Gaussianity in the inflationary perturbations can couple observable scales to modes of much longer wavelength (even superhorizon), leaving as a signature a large-angle modulation of the observed cosmic microwave background (CMB) power spectrum. This provides an alternative origin for a power asymmetry which is otherwise often ascribed to a breaking of statistical isotropy. The non-Gaussian modulation effect can be significant even for typical ~10^{-5} perturbations, while respecting current constraints on non-Gaussianity, if the squeezed limit of the bispectrum is sufficiently infrared divergent. Just such a strongly infrared-divergent bispectrum has been claimed for inflation models with a non-Bunch-Davies initial state, for instance. Upper limits on the observed CMB power asymmetry place stringent constraints on the duration of inflation in such models.

preprint2012arXiv

Conformal Symmetries of Adiabatic Modes in Cosmology

We remark on the existence of non-linearly realized conformal symmetries for scalar adiabatic perturbations in cosmology. These conformal symmetries are present for any cosmological background, beyond any slow-roll or quasi-de Sitter approximation. The dilatation transformation shifts the curvature perturbation by a constant, and corresponds to the well-known symmetry under spatial rescaling. We argue that the scalar sector is also invariant under special conformal transformations, which shift the curvature perturbation by a term linear in the spatial coordinates. We discuss whether these conformal symmetries can be extended to include tensor perturbations. Tensor modes introduce their own set of non-linearly realized symmetries. We identify an infinite set of large gauge transformations which maintain the transverse, traceless gauge condition, while shifting the tensor mode non-trivially.

preprint2012arXiv

Introducing PHAEDRA: a new spectral code for simulations of relativistic magnetospheres

We describe a new scheme for evolving the equations of force-free electrodynamics, the vanishing-inertia limit of magnetohydrodynamics. This pseudospectral code uses global orthogonal basis function expansions to take accurate spatial derivatives, allowing the use of an unstaggered mesh and the complete force-free current density. The method has low numerical dissipation and diffusion outside of singular current sheets. We present a range of one- and two-dimensional tests, and demonstrate convergence to both smooth and discontinuous analytic solutions. As a first application, we revisit the aligned rotator problem, obtaining a steady solution with resistivity localised in the equatorial current sheet outside the light cylinder.

preprint2012arXiv

Large-scale Bias and Efficient Generation of Initial Conditions for Non-Local Primordial Non-Gaussianity

We study the scale-dependence of halo bias in generic (non-local) primordial non-Gaussian (PNG) initial conditions of the type motivated by inflation, parametrized by an arbitrary quadratic kernel. We first show how to generate non-local PNG initial conditions with minimal overhead compared to local PNG models for a general class of primordial bispectra that can be written as linear combinations of separable templates. We run cosmological simulations for the local, and non-local equilateral and orthogonal models and present results on the scale-dependence of halo bias. We also derive a general formula for the Fourier-space bias using the peak-background split (PBS) in the context of the excursion set approach to halos and discuss the difference and similarities with the known corresponding result from local bias models. Our PBS bias formula generalizes previous results in the literature to include non-Markovian effects and non-universality of the mass function and are in better agreement with measurements in numerical simulations than previous results for a variety of halo masses, redshifts and halo definitions. We also derive for the first time quadratic bias results for arbitrary non-local PNG, and show that non-linear bias loops give small corrections at large-scales. The resulting well-behaved perturbation theory paves the way to constrain non-local PNG from measurements of the power spectrum and bispectrum in galaxy redshift surveys.

preprint2012arXiv

No-Go Theorems for Generalized Chameleon Field Theories

The chameleon, or generalizations thereof, is a light scalar that couple to matter with gravitational strength, but whose manifestation depends on the ambient matter density. A key feature is that the screening mechanism suppressing its effects in high-density environments is determined by the local scalar field value. Under very general conditions, we prove two theorems limiting its cosmological impact: i) the Compton wavelength of such a scalar can be at most Mpc at present cosmic density, which restricts its impact to non-linear scales; ii) the conformal factor relating Einstein- and Jordan-frame scale factors is essentially constant over the last Hubble time, which precludes the possibility of self-acceleration. These results imply that chameleon-like scalar fields have a negligible effect on the linear-scale growth history; theories that invoke a chameleon-like scalar to explain cosmic acceleration rely on a form of dark energy rather than a genuine modified gravity effect. Our analysis applies to a broad class of chameleon, symmetron and dilaton theories.

preprint2012arXiv

Twisting, reconnecting magnetospheres and magnetar spindown

We present the first simulations of evolving, strongly twisted magnetar magnetospheres. Slow shearing of the magnetar crust is seen to lead to a series of magnetospheric expansion and reconnection events, corresponding to X-ray flares and bursts. The axisymmetric simulations include rotation of the neutron star and the magnetic wind through the light cylinder. We study how the increasing twist affects the spindown rate of the star, finding that a dramatic increase in spindown occurs. Particularly spectacular are explosive events caused by the sudden opening of large amounts of overtwisted magnetic flux, which may be associated with the observed giant flares. These events are accompanied by a short period of ultra-strong spindown, resulting in an abrupt increase in spin period, such as was observed in the giant flare of SGR 1900+14.

preprint2011arXiv

Anisotropic Extinction Distortion of the Galaxy Correlation Function

Similar to the magnification of the galaxies' fluxes by gravitational lensing, the extinction of the fluxes by comic dust, whose existence is recently detected by Menard et al (2009), also modify the distribution of a flux-selected galaxy sample. We study the anisotropic distortion by dust extinction to the 3D galaxy correlation function, including magnification bias and redshift distortion at the same time. We find the extinction distortion is most significant along the line of sight and at large separations, similar to that by magnification bias. The correction from dust extinction is negative except at sufficiently large transverse separations, which is almost always opposite to that from magnification bias (we consider a number count slope s > 0.4). Hence, the distortions from these two effects tend to reduce each other. At low z (~<1), the distortion by extinction is stronger than that by magnification bias, but at high z, the reverse holds. We also study how dust extinction affects probes in real space of the baryon acoustic oscillations (BAO) and the linear redshift distortion parameter beta. We find its effect on BAO is negligible. However, it introduces a positive scale-dependent correction to beta that can be as large as a few percent. At the same time, we also find a negative scale-dependent correction from magnification bias, which is up to percent level at low z, but to ~40% at high z. These corrections are non-negligible for precision cosmology, and should be considered when testing General Relativity through the scale-dependence of beta.

preprint2011arXiv

Effective field theory for hydrodynamics: thermodynamics, and the derivative expansion

We consider the low-energy effective field theory describing the infrared dynamics of non-dissipative fluids. We extend previous work to accommodate conserved charges, and we clarify the matching between field theory variables and thermodynamical ones. We discuss the systematics of the derivative expansion, for which field theory offers a conceptually clear and technically neat scheme. As an example, we compute the correction to the sound-wave dispersion relation coming from a sample second-order term. This formalism forms the basis for a study of anomalies in hydrodynamics via effective field theory, which is initiated in a companion paper.

preprint2011arXiv

Effective field theory for hydrodynamics: Wess-Zumino term and anomalies in two spacetime dimensions

We develop the formalism that incorporates quantum anomalies in the effective field theory of non-dissipative fluids. We consider the effect of adding a Wess-Zumino-like term to the low-energy effective action to account for anomalies. In this paper we restrict to two spacetime dimensions. We find modifications to the constitutive relations for the current and the stress-energy tensor, and, more interestingly, half a new propagating mode (one-and-a-halfth sound): a left- or right-moving wave with propagation speed that goes to zero with the anomaly coefficient. Unlike for the chiral magnetic wave in four dimensions, this mode propagates even in the absence of external fields. We check our results against a more standard, purely hydrodynamical derivation. Unitarity of the effective field theory suggests an upper bound on the anomaly coefficient in hydrodynamics.

preprint2011arXiv

Scale-dependent halo bias from scale-dependent growth

We derive a general expression for the large-scale halo bias, in theories with a scale-dependent linear growth, using the excursion set formalism. Such theories include modified gravity models, and models in which the dark energy clustering is non-negligible. A scale dependence is imprinted in both the formation and evolved biases by the scale-dependent growth. Mergers are accounted for in our derivation, which thus extends earlier work which focused on passive evolution. There is a simple analytic form for the bias for those theories in which the nonlinear collapse of perturbations is approximately the same as in general relativity. As an illustration, we apply our results to a simple Yukawa modification of gravity, and use SDSS measurements of the clustering of luminous red galaxies to constrain the theory's parameters.

preprint2011arXiv

The effect of peculiar velocities on supernova cosmology

We present an analysis of peculiar velocities and their effect on supernova cosmology. In particular, we study (a) the corrections due to our own motion, (b) the effects of correlations in peculiar velocities induced by large-scale structure, and (c) uncertainties arising from a possible local under- or over-density. For all of these effects we present a case study of their impact on the cosmology derived by the Sloan Digital Sky Survey-II Supernova Survey (SDSS-II SN Survey). Correcting supernova redshifts for the CMB dipole slightly over-corrects nearby supernovae that share some of our local motion. We show that while neglecting the CMB dipole would cause a shift in the derived equation of state of Delta w ~ 0.04 (at fixed matter density) the additional local-motion correction is currently negligible (Delta w<0.01). We use a covariance-matrix approach to statistically account for correlated peculiar velocities. This down-weights nearby supernovae and effectively acts as a graduated version of the usual sharp low-redshift cut. Neglecting coherent velocities in the current sample causes a systematic shift of ~2% in the preferred value of w and will therefore have to be considered carefully when future surveys aim for percent-level accuracy. Finally, we perform n-body simulations to estimate the likely magnitude of any local density fluctuation (monopole) and estimate the impact as a function of the low-redshift cutoff. We see that for this aspect the low-z cutoff of z=0.02 is well-justified theoretically, but that living in a putative local density fluctuation leaves an indelible imprint on the magnitude-redshift relation.

preprint2010arXiv

An equivalence principle for scalar forces

The equivalence of inertial and gravitational masses is a defining feature of general relativity. Here, we clarify the status of the equivalence principle for interactions mediated by a universally coupled scalar, motivated partly by recent attempts to modify gravity at cosmological distances. Although a universal scalar-matter coupling is not mandatory, once postulated, it is stable against classical and quantum renormalizations in the matter sector. The coupling strength itself is subject to renormalization of course. The scalar equivalence principle is violated only for objects for which either the graviton self-interaction or the scalar self-interaction is important---the first applies to black holes, while the second type of violation is avoided if the scalar is Galilean-symmetric.

preprint2010arXiv

Derrick's theorem beyond a potential

Scalar field theories with derivative interactions are known to possess solitonic excitations, but such solitons are generally unsatisfactory because the effective theory fails precisely where nonlinearities responsible for the solitons are important. A new class of theories possessing (internal) galilean invariance can in principle bypass this difficulty. Here, we show that these galileon theories do not possess stable solitonic solutions. As a by-product, we show that no stable solitons exist for a different class of derivatively coupled theories, describing for instance the infrared dynamics of superfluids, fluids, solids and some k-essence models.

preprint2010arXiv

Gravitational Lensing as Signal and Noise in Lyman-alpha Forest Measurements

In Lyman-alpha forest measurements it is generally assumed that quasars are mere background light sources which are uncorrelated with the forest. Gravitational lensing of the quasars violates this assumption. This effect leads to a measurement bias, but more interestingly it provides a valuable signal. The lensing signal can be extracted by correlating quasar magnitudes with the flux power spectrum and with the flux decrement. These correlations will be challenging to measure but their detection provides a direct measure of how features in the Lyman-alpha forest trace the underlying mass density field. Observing them will test the fundamental hypothesis that fluctuations in the forest are predominantly driven by fluctuations in mass, rather than in the ionizing background, helium reionization or winds. We discuss ways to disentangle the lensing signal from other sources of such correlations, including dust, continuum and background residuals. The lensing-induced measurement bias arises from sample selection: one preferentially collects spectra of magnified quasars which are behind overdense regions. This measurement bias is ~0.1-1% for the flux power spectrum, optical depth and the flux probability distribution. Since the effect is systematic, quantities such as the amplitude of the flux power spectrum averaged across scales should be interpreted with care.

preprint2010arXiv

How to Run Through Walls: Dynamics of Bubble and Soliton Collisions

It has recently been shown in high resolution numerical simulations that relativistic collisions of bubbles in the context of a multi-vacua potential may lead to the creation of bubbles in a new vacuum. In this paper, we show that scalar fields with only potential interactions behave like free fields during high-speed collisions; the kick received by them in a collision can be deduced simply by a linear superposition of the bubble wall profiles. This process is equivalent to the scattering of solitons in 1+1 dimensions. We deduce an expression for the field excursion (shortly after a collision), which is related simply to the field difference between the parent and bubble vacua, i.e. contrary to expectations, the excursion cannot be made arbitrarily large by raising the collision energy. There is however a minimum energy threshold for this excursion to be realized. We verify these predictions using a number of 3+1 and 1+1 numerical simulations. A rich phenomenology follows from these collision induced excursions - they provide a new mechanism for scanning the landscape, they might end/begin inflation, and they might constitute our very own big bang, leaving behind a potentially observable anisotropy.

preprint2010arXiv

Stellar Structure and Tests of Modified Gravity

Theories that attempt to explain cosmic acceleration by modifying gravity typically introduces a long-range scalar force that needs to be screened on small scales. One common screening mechanism is the chameleon, where the scalar force is screened in environments with a sufficiently deep gravitational potential, but acts unimpeded in regions with a shallow gravitational potential. This leads to a variation in the overall gravitational G with environment. We show such a variation can occur within a star itself, significantly affecting its evolution and structure, provided that the host galaxy is unscreened. The effect is most pronounced for red giants, which would be smaller by a factor of tens of percent and thus hotter by 100's of K, depending on the parameters of the underlying scalar-tensor theory. Careful measurements of these stars in suitable environments (nearby dwarf galaxies not associated with groups or clusters) would provide constraints on the chameleon mechanism that are four orders of magnitude better than current large scale structure limits, and two orders of magnitude better than present solar system tests.

preprint2009arXiv

Clustering of Luminous Red Galaxies IV: Baryon Acoustic Peak in the Line-of-Sight Direction and a Direct Measurement of H(z)

We study the clustering of LRG galaxies in the latest spectroscopic SDSS data releases, DR6 and DR7, which sample over 1 Gpc^3/h^3 to z=0.47. The 2-point correlation function $\xisp$ is estimated as a function of perpendicular $σ$ and line-of-sight $π$ (radial) directions. We find a significant detection of a peak at $r\simeq 110$Mpc/h, which shows as a circular ring in the $σ-π$ plane. There is also significant evidence for a peak along the radial direction whose shape is consistent with its originating from the recombination-epoch baryon acoustic oscillations (BAO). A $\xisp$ model with no radial BAO peak is disfavored at $3.2σ$, whereas a model with no magnification bias is disfavored at $2σ$. The radial data enable, for the first time, a direct measurement of the Hubble parameter $H(z)$ as a function of redshift. This is independent from earlier BAO measurements which used the spherically averaged (monopole) correlation to constrain an integral of $H(z)$. Using the BAO peak position as a standard ruler in the radial direction, we find: $H(z=0.24)= 79.69 \pm 2.32 (\pm 1.29)$ km/s/Mpc for z=0.15-0.30 and $H(z=0.43)= 86.45 \pm 3.27 (\pm 1.69)$ km/s/Mpc for $z=0.40-0.47$. The first error is a model independent statistical estimation and the second accounts for systematics both in the measurements and in the model. For the full sample, $z=0.15-0.47$, we find $H(z=0.34)= 83.80 \pm 2.96 (\pm 1.59)$ km/s/Mpc.

preprint2008arXiv

Challenges to the DGP Model from Horizon-Scale Growth and Geometry

We conduct a Markov Chain Monte Carlo study of the Dvali-Gabadadze-Porrati (DGP) self-accelerating braneworld scenario given the cosmic microwave background (CMB) anisotropy, supernovae and Hubble constant data by implementing an effective dark energy prescription for modified gravity into a standard Einstein-Boltzmann code. We find no way to alleviate the tension between distance measures and horizon scale growth in this model. Growth alterations due to perturbations propagating into the bulk appear as excess CMB anisotropy at the lowest multipoles. In a flat cosmology, the maximum likelihood DGP model is nominally a 5.3 sigma poorer fit than Lambda CDM. Curvature can reduce the tension between distance measures but only at the expense of exacerbating the problem with growth leading to a 4.8 sigma result that is dominated by the low multipole CMB temperature spectrum. While changing the initial conditions to reduce large scale power can flatten the temperature spectrum, this also suppresses the large angle polarization spectrum in violation of recent results from WMAP5. The failure of this model highlights the power of combining growth and distance measures in cosmology as a test of gravity on the largest scales.

preprint2003arXiv

Line versus Flux Statistics -- Considerations for the Low Redshift Lyman-alpha Forest

The flux/transmission power spectrum has become a popular statistical tool in studies of the high redshift ($z > 2$) Lyman-alpha forest. At low redshifts, where the forest has thinned out into a series of well-isolated absorption lines, the motivation for flux statistics is less obvious. Here, we study the relative merits of flux versus line correlations, and derive a simple condition under which one is favored over the other on purely statistical grounds. Systematic errors probably play an important role in this discussion, and they are outlined as well.

preprint2000arXiv

Constraining the Lifetime of Quasars from their Spatial Clustering

The lifetime t_Q of the luminous phase of quasars is constrained by current observations to be between 10^6 and 10^8 years, but is otherwise unkown. We model the quasar luminosity function in detail in the optical and X-ray bands using the Press-Schechter formalism, and show that the expected clustering of quasars depends strongly on their assumed lifetime. We quantify this dependence, and find that existing measurements of the correlation length of quasars are consistent with the range 10^6 < t_Q < 10^8 years. We then show that future measurements of the power spectrum of quasars out to z=3, from the 2dF or Sloan Digital Sky Survey, can significantly improve this constraint, and in principle allow a precise determination of t_Q. We estimate the systematic errors introduced by uncertainties in the modeling of the quasar-halo relationship, as well as by the possible existence of obscured quasars.

preprint1999arXiv

Multiple-streaming and the Probability Distribution of Density in Redshift Space

We examine several aspects of redshift distortions by expressing the redshift-space density in terms of the eigenvalues and orientation of the local Lagrangian deformation tensor. We explore the importance of multiple-streaming using the Zel'dovich approximation (ZA), and compute the average number of streams in real and redshift-space. It is found that multiple-streaming can be significant in redshift-space but negligible in real-space, even at moderate values of the linear fluctuation amplitude ($σ< 1$). Moreover, unlike their real-space counter-parts, redshift-space multiple-streams can flow past each other with minimal interactions. Such nonlinear redshift-space effects, which operate even when the real-space density field is quite linear, could suppress the classic compression of redshift-structures predicted by linear theory (Kaiser 1987). We also compute using the ZA the probability distribution function (PDF) of density, as well as $S_3$, in real and redshift-space, and compare it with the PDF measured from N-body simulations. The role of caustics in defining the character of the high density tail is examined. It is found that (non-Lagrangian) smoothing, due to both finite resolution or discreteness and small-scale velocity dispersions, is very effective in erasing caustic structures, unless the initial power spectrum is sufficiently truncated.

preprint1999arXiv

Superheavy Dark Matter from Thermal Inflation

It is quite plausible that the mass of the dark matter particle increases significantly after its freeze-out, due to a scalar field rolling to large values. We describe a realization of this scenario in the context of thermal inflation which naturally gives a cold dark matter particle with the correct cosmological abundance and a mass around $10^{10} {\rm GeV}$, evading the conventional upper bound of $10^5 {\rm GeV}$. We also discuss another realization which could produce a cosmologically interesting abundance of near Planck mass, possibly electromagnetically charged, particles. The detection and observational consequences of superheavy cold dark matter or WIMPZILLAs are briefly examined.

preprint1999arXiv

Weighing the Cosmological Energy Contents with Weak Gravitational Lensing

Bernardeau et al. (1997), using perturbation theory, showed that the skewness of the large-scale lensing-convergence, or projected mass density, could be used to constrain $Ω_m$, the matter content of the universe. On the other hand, deep weak-lensing field surveys in the near future will likely measure the convergence on small angular scales (< 10 arcmin.), where the signal will be dominated by highly nonlinear fluctuations. We develop a method to compute the small-scale convergence skewness, using a prescription for the highly nonlinear three-point function developed by Scoccimarro and Frieman (1998). This method gives predictions that agree well with existing results from ray-tracing N-body simulations, but is significantly faster, allowing the exploration of a large number of models. We demonstrate that the small-scale convergence skewness is insensitive to the shape and normalization of the primordial (CDM-type) power spectrum, making it dependent almost entirely on the cosmological energy contents, through their influence on the global geometrical distances and fluctuation growth rate. Moreover, nonlinear clustering appears to enhance the differences between predictions of the convergence skewness for a range of models. Hence, in addition to constraining $Ω_m$, the small-scale convergence skewness from future deep several- degree-wide surveys can be used to differentiate between curvature dominated and cosmological constant ($Λ$) dominated models, as well as to constrain the equation of state of a quintessence component, thereby distinguishing $Λ$ from quintessence as well. Finally, our method can be easily generalized to other measures such as aperture mass statistics.

preprint1996arXiv

Equation of State of the Photoionized Intergalactic Medium

We develop an efficient method to study the effects of reionization history on the temperature-density relation of the intergalactic medium in the low density limit (overdensity less than 5). It is applied to the study of photo-reionization models in which the amplitude, spectrum and onset epoch of the ionizing flux, as well as the cosmology, are systematically varied. We find that the mean temperature-density relation at z=2-4 is well approximated by a power-law equation of state for uniform reionization models. We derive analytical expressions for its evolution and exhibit its asymptotic behavior: it is found that for sufficiently early reionization, imprints of reionization history prior to z=10 on the temperature-density relation are washed out. In this limit the temperature at cosmic mean density is proportional to (Ω_b h/\sqrtΩ_0)^{1/1.7}. While the amplitude of the radiation flux at the ionizing frequency of HI is found to have a negligible effect on the temperature-density relation as long as the universe reionizes before z=5, the spectrum can change the overall temperature by about 20%, through variations in the abundances of helium species. However the slope of the mean equation of state is found to lie within a narrow range for all reionization models we study, where reionization takes place before z=5. We discuss the implications of these findings for the observational properties of the Lyman-alpha forest. In particular, uncertainties in the temperature of the intergalactic medium, due to the uncertain reionization history of our universe, introduces a 30% scaling in the amplitude of the column density distribution while the the slope of the distribution is only affected by about 5%. Finally, we discuss how a fluctuating ionizing field affects the above results. We argue that under

preprint1995arXiv

Local Approximations to the Gravitational Collapse of Cold Matter

We investigate three different local approximations for nonlinear gravitational instability in the framework of cosmological Lagrangian fluid dynamics of cold dust. They include the Zel'dovich approximation (ZA), the ``non-magnetic'' approximation of Bertschinger \& Jain (1994, NMA), and a new ``local tidal'' approximation (LTA). The LTA is exact for any perturbations whose gravitational and velocity equipotentials have the same constant shape with time, including spherical, cylindrical, and plane-parallel perturbations. We tested all three local approximations with the collapse of a homogeneous triaxial ellipsoid, for which an exact solution exists for an ellipsoid embedded in empty space and an excellent approximation is known in the cosmological context. We find that the LTA is significantly more accurate in general than the ZA and the NMA. Like the ZA, but unlike the NMA, the LTA generically leads to pancake collapse. For a randomly chosen mass element in an Einstein-de Sitter universe, assuming a Gaussian random field of initial density fluctuations, the LTA predicts that at least 78\% of initially underdense regions collapse owing to nonlinear effects of shear and tides.

preprint1995arXiv

Standard Cosmology and the BATSE Number vs. Peak Flux Distribution

The observed 2B BATSE distribution is consistent with the faintest GRBs in our sample originating from a redshift of Zmax ~ 0.8-3.0 (90\%), with the most likely values in the range of 1.0-2.2, and is largely insensitive to Omega for models with no evolution. To constrain the model parameter Omega to the range 0.1-1.0 using only Log N -- Log P distributions, more than 4000 GRBs, with a most likely value of ~ 9,000 GRBs to BATSE sensitivity. This requires a live integration time of >6 years with BATSE. Detectors sensitive to much lower limits (~ 70-400 in sensitivity) require ~ 200 GRBs, with <0.03 year 4pi ster coverage. We place limits on the amount of frequency density and, separately, peak luminosity evolution in the sample of GRBs. We find that frequency density evolution models can place the faintest GRBs at Zmax ~ 10-200, without conflicting with the observations of relative time dilation of ~ 2.