Source author record

Yacine Ali-Haïmoud

Yacine Ali-Haïmoud 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

26works
8topics
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

26 published item(s)

preprint2022arXiv

Astrophysical and Cosmological Probes of Dark Matter

While astrophysical and cosmological probes provide a remarkably precise and consistent picture of the quantity and general properties of dark matter, its fundamental nature remains one of the most significant open questions in physics. Obtaining a more comprehensive understanding of dark matter within the next decade will require overcoming a number of theoretical challenges: the groundwork for these strides is being laid now, yet much remains to be done. Chief among the upcoming challenges is establishing the theoretical foundation needed to harness the full potential of new observables in the astrophysical and cosmological domains, spanning the early Universe to the inner portions of galaxies and the stars therein. Identifying the nature of dark matter will also entail repurposing and implementing a wide range of theoretical techniques from outside the typical toolkit of astrophysics, ranging from effective field theory to the dramatically evolving world of machine learning and artificial-intelligence-based statistical inference. Through this work, the theory frontier will be at the heart of dark matter discoveries in the upcoming decade.

preprint2021arXiv

Insights into searches for anisotropies in the nanohertz gravitational-wave background

Within the next several years pulsar timing arrays (PTAs) are positioned to detect the stochastic gravitational-wave background (GWB) likely produced by the collection of inspiralling super-massive black holes binaries, and potentially constrain some exotic physics. So far most of the pulsar timing data analysis has focused on the monopole of the GWB, assuming it is perfectly isotropic. The natural next step is to search for anisotropies in the GWB. In this paper, we use the recently developed PTA Fisher matrix to gain insights into optimal search strategies for GWB anisotropies. For concreteness, we apply our results to EPTA data, using realistic noise characteristics of its pulsars. We project the detectability of a GWB whose angular dependence is assumed to be a linear combination of predetermined maps, such as spherical harmonics or coarse pixels. We find that the GWB monopole is always statistically correlated with these maps, implying a loss of sensitivity to the monopole when searching simultaneously for anisotropies. We then derive the angular distributions of the GWB intensity to which a PTA is most sensitive, and illustrate how one may use these "principal maps" to approximately reconstruct the angular dependence of the GWB. Since the principal maps are neither perfectly anisotropic nor uncorrelated with the monopole, we also develop a frequentist criterion to specifically search for anisotropies in the GWB without any prior knowledge about their angular distribution. Lastly, we show how to recover existing EPTA results with our Fisher formalism, and clarify their meaning. The tools presented here will be valuable in guiding and optimizing the computationally demanding analyses of pulsar timing data.

preprint2021arXiv

Testing dark matter interactions with CMB spectral distortions

Possible interactions of dark matter (DM) with Standard Model (SM) particles can be tested with spectral distortions (SDs) of the cosmic microwave background (CMB). In particular, a non-relativistic DM particle that scatters elastically with photons, electrons or nuclei imprints a negative chemical potential $μ$ to the CMB spectrum. This article revisits the first study of this effect, with an accurate treatment of heat exchange between DM and SM particles. We show that the instantaneous-decoupling approximation made in the original study systematically and significantly underestimates the amplitude of SDs. As a consequence, we derive tighter upper bounds to the DM-SM elastic-scattering cross section for DM masses $m_χ\lesssim 0.1$ MeV, from the non-detection of $μ$-distortions by FIRAS. We also show that a future instrument like PIXIE, sensitive to $|μ| \sim 10^{-8}$, would be able to probe DM-SM cross sections much smaller than first forecasted, and orders of magnitude below current upper limits from CMB-anisotropy data, up to DM masses of $\sim 1$ GeV. Lastly, we study the sensitivity of SDs to the electric and magnetic dipole moments of the DM. Although SDs can place non-trivial constraints on these models, we find that even future SD experiments are unlikely to improve upon the best current bounds. This article is accompanied by the public code DMDIST, which allows one to compute CMB SDs for generic particle-DM models, specified by their cross sections for elastic scattering with and annihilation into SM particles.

preprint2020arXiv

Fisher formalism for anisotropic gravitational-wave background searches with pulsar timing arrays

Pulsar timing arrays (PTAs) are currently the only experiments directly sensitive to gravitational waves with decade-long periods. Within the next five to ten years, PTAs are expected to detect the stochastic gravitational-wave background (SGWB) collectively sourced by inspiralling supermassive black hole binaries. It is expected that this background is mostly isotropic, and current searches focus on the monopole part of the SGWB. Looking ahead, anisotropies in the SGWB may provide a trove of additional information both on known and unknown astrophysical and cosmological sources. In this paper, we build a simple yet realistic Fisher formalism for anisotropic SGWB searches with PTAs. Our formalism is able to accommodate realistic properties of PTAs, and allows simple and accurate forecasts. We illustrate our approach with an idealized PTA consisting of identical, isotropically distributed pulsars. In a companion paper, we apply our formalism to current PTAs and show that it can be a powerful tool to guide and optimize real data analysis.

preprint2020arXiv

H0 tension or T0 tension?

We study if the discrepancy between the local and cosmological measurements of the Hubble constant $H_0$ can be reformulated as a tension in the cosmic microwave background (CMB) monopole temperature $T_0$. The latter is customarily fixed to the FIRAS best-fit value in CMB data analyses. Although this value was confirmed by several independent experiments, it is interesting to see how much parameter constraints depend on this prior. We first provide a detailed pedagogical description of the $T_0$ effects on cosmological observables. We show that the recombination history and transfer functions do not depend on $T_0$, provided they are parametrized by the energy scale rather than redshift, and at fixed dark matter and baryon densities per CMB photon. Thus, $T_0$ is only a property of the observer, quantifying the amount of expansion between key cosmological events and today. As a consequence, the sole effect of $T_0$ on small-scale primary CMB anisotropies is through the angular diameter distance to the epoch of last scattering, resulting in a near-perfect degeneracy between $T_0$ and $H_0$. This geometric degeneracy is partially lifted by the late-time integrated Sachs-Wolfe effect and CMB lensing. Still, Planck data alone is consistent with a broad region in the $H_0-T_0$ plane, implying that removing the FIRAS prior on $T_0$ can make Planck and SH0ES less discrepant, without introducing new physics beyond $Λ$CDM. One may break the degeneracy by combining Planck with SH0ES, yielding an independent measurement of $T_0$, which happens to be in $3σ$ tension with FIRAS. Therefore, the Hubble tension indeed can be recast into the $T_0$ tension. The agreement with FIRAS is restored when combining Planck with the baryon acoustic oscillation data instead of SH0ES. Thus, the tension between SH0ES and cosmological measurements of $H_0$ persists even if we discard the FIRAS $T_0$ prior.

preprint2016arXiv

Did LIGO detect dark matter?

We consider the possibility that the black-hole (BH) binary detected by LIGO may be a signature of dark matter. Interestingly enough, there remains a window for masses $20\,M_\odot \lesssim M_{\rm bh} \lesssim 100\, M_\odot$ where primordial black holes (PBHs) may constitute the dark matter. If two BHs in a galactic halo pass sufficiently close, they radiate enough energy in gravitational waves to become gravitationally bound. The bound BHs will rapidly spiral inward due to emission of gravitational radiation and ultimately merge. Uncertainties in the rate for such events arise from our imprecise knowledge of the phase-space structure of galactic halos on the smallest scales. Still, reasonable estimates span a range that overlaps the $2-53$ Gpc$^{-3}$ yr$^{-1}$ rate estimated from GW150914, thus raising the possibility that LIGO has detected PBH dark matter. PBH mergers are likely to be distributed spatially more like dark matter than luminous matter and have no optical nor neutrino counterparts. They may be distinguished from mergers of BHs from more traditional astrophysical sources through the observed mass spectrum, their high ellipticities, or their stochastic gravitational wave background. Next generation experiments will be invaluable in performing these tests.

preprint2016arXiv

Flaring of tidally compressed dark-matter clumps

We explore the physics and observational consequences of tidal compression events (TCEs) of dark-matter clumps (DMCs) by supermassive black holes (SMBHs). Our analytic calculations show that a DMC approaching a SMBH much closer than the tidal radius undergoes significant compression along the axis perpendicular to the orbital plane, shortly after pericenter passage. For DMCs composed of self-annihilating dark-matter particles, we find that the boosted DMC density and velocity dispersion lead to a flaring of the annihilation rate, most pronounced for a velocity- dependent annihilation cross section. If the end products of the annihilation are photons, this results in a gamma-ray flare, detectable (and possibly already detected) by the Fermi telescope for a range of model parameters. If the end products of dark-matter annihilation are relativistic electrons and positrons and the local magnetic field is large enough, TCEs of DMCs can lead to flares of synchrotron radiation. Finally, TCEs of DMCs lead to a burst of gravitational waves, in addition to the ones radiated by the orbital motion alone, and with a different frequency spectrum. These transient phenomena provide interesting new avenues to explore the properties of dark matter.

preprint2016arXiv

On CMB B-Mode Non-Gaussianity

We study the degree to which the cosmic microwave background (CMB) can be used to constrain primordial non-Gaussianity involving one tensor and two scalar fluctuations, focusing on the correlation of one polarization $B$ mode with two temperature modes. In the simplest models of inflation, the tensor-scalar-scalar primordial bispectrum is non-vanishing and is of the same order in slow-roll parameters as the scalar-scalar-scalar bispectrum. We calculate the $\langle BTT\rangle$ correlation arising from a primordial tensor-scalar-scalar bispectrum, and show that constraints from an experiment like CMB-Stage IV using this observable are more than an order of magnitude better than those on the same primordial coupling obtained from temperature measurements alone. We argue that $B$-mode non-Gaussianity opens up an as-yet-unexplored window into the early Universe, demonstrating that significant information on primordial physics remains to be harvested from CMB anisotropies.

preprint2016arXiv

Orbital eccentricities in primordial black holes binaries

It was recently suggested that the merger of $\sim30\,M_\odot$ primordial black holes (PBHs) may provide a significant number of events in gravitational-wave observatories over the next decade, if they make up an appreciable fraction of the dark matter. Here we show that measurement of the eccentricities of the inspiralling binary black holes can be used to distinguish these binaries from those produced by more traditional astrophysical mechanisms. These PBH binaries are formed on highly eccentric orbits and can then merge on timescales that in some cases are years or less, retaining some eccentricity in the last seconds before the merger. This is to be contrasted with massive-stellar-binary, globular-cluster, or other astrophysical origins for binary black holes (BBHs) in which the orbits have very effectively circularized by the time the BBH enters the observable LIGO window. Here we discuss the features of the gravitational-wave signals that indicate this eccentricity and forecast the sensitivity of LIGO and the Einstein Telescope to such effects. We show that if PBHs make up the dark matter, then roughly one event should have a detectable eccentricity given LIGO's expected sensitivity and observing time of six years. The Einstein Telescope should see $O(10)$ such events after ten years.

preprint2015arXiv

Constraints on dark matter interactions with standard model particles from CMB spectral distortions

We propose a new method to constrain elastic scattering between dark matter (DM) and standard model particles in the early Universe. Direct or indirect thermal coupling of non-relativistic DM with photons leads to a heat sink for the latter. This results in spectral distortions of the cosmic microwave background (CMB), the amplitude of which can be as large as a few times the DM-to-photon number ratio. We compute CMB spectral distortions due to DM-proton, DM-electron and DM-photon scattering for generic energy-dependent cross sections and DM mass m_DM >~ 1 keV. Using FIRAS measurements we set constraints on the cross sections for m_DM <~ 0.1 MeV. In particular, for energy-independent scattering we obtain sigma[DM-proton] <~ 10^(-24) cm^2 (keV/m_DM)^(1/2), sigma[DM-electron] <~ 10^(-27) cm^2 (keV/m_DM)^(1/2) and sigma[DM-photon] <~ 10^(-39) cm^2 (m_DM/keV). An experiment with the characteristics of PIXIE would extend the regime of sensitivity up to masses m_DM ~ 1 GeV.

preprint2015arXiv

Heating of Baryons due to Scattering with Dark Matter During the Dark Ages

We explore the effects of elastic scattering between dark matter and baryons on the 21-cm signal during the dark ages. In particular, we consider a dark-matter---baryon interaction with a cross section of the form $σ= σ_0 v^{-4}$, in which case the effect of the drag force between the dark matter and baryon fluids grows with time. We show that, as opposed to what was previously thought, this effect heats up the baryons due to the relative velocity between dark matter and baryons. This creates an additional source of fluctuations, which can potentially make interactions easier to detect by 21-cm measurements than by using the cosmic microwave background and the Lyman-$α$ forest. Our forecasts show that the magnitude of the cross section can be probed to $σ_0\sim 3\times 10^{-42}$ cm$^2$ for $m_χ\ll 1$ GeV and $σ_0\sim 2 \times 10^{-41}\ (m_χ/10\, \rm GeV)$ cm$^2$ for $m_χ\gg 1$ GeV with next generation experiments, and improved to $σ_0\sim 4\times 10^{-44}$ cm$^2$ for $m_χ \ll 1$ GeV and $σ_0\sim 4 \times 10^{-43}\ (m_χ/10\, \rm GeV)$ cm$^2$ for $m_χ\gg 1$ GeV with futuristic experiments.

preprint2015arXiv

Perturbative interaction approach to cosmological structure formation

A new approach to cosmological perturbation theory has been recently introduced by Bartelmann et al., relying on nonequilibrium statistical theory of classical particles, and treating the gravitational interaction perturbatively. They compute analytic expressions for the nonlinear matter power spectrum, to first order in the interaction, and at one-loop order in the linear power spectrum. The resulting power spectrum is well behaved even at large wavenumbers and seems in good agreement with results from numerical simulations. In this paper, we rederive their results concisely with a different approach, starting from the implicit integral solution to particle trajectories. We derive the matter power spectrum to first order in the interaction, but to arbitrary order in the linear power spectrum, from which the one-loop result follows. We also show that standard linear perturbation theory can only be recovered at infinite order in the gravitational interaction. At finite order in the interaction, we find that the linear power spectrum is systematically and significantly underestimated. A comprehensive study of the convergence of the theory with the order of the interaction for nonlinear scales will be the subject of future work.

preprint2015arXiv

Primordial non-gaussianity from the bispectrum of 21-cm fluctuations in the dark ages

A measurement of primordial non-gaussianity will be of paramount importance to distinguish between different models of inflation. Cosmic microwave background (CMB) anisotropy observations have set unprecedented bounds on the non-gaussianity parameter f_NL but the interesting regime f_NL <~ 1 is beyond their reach. Brightness-temperature fluctuations in the 21-cm line during the dark ages (z ~ 30-100) are a promising successor to CMB studies, giving access to a much larger number of modes. They are, however, intrinsically non-linear, which results in secondary non-gaussianities orders of magnitude larger than the sought-after primordial signal. In this paper we carefully compute the primary and secondary bispectra of 21-cm fluctuations on small scales. We use the flat-sky formalism, which greatly simplifies the analysis, while still being very accurate on small angular scales. We show that the secondary bispectrum is highly degenerate with the primordial one, and argue that even percent-level uncertainties in the amplitude of the former lead to a bias of order Delta f_NL ~ 10. To tackle this problem we carry out a detailed Fisher analysis, marginalizing over the amplitudes of a few smooth redshift-dependent coefficients characterizing the secondary bispectrum. We find that the signal-to-noise ratio for a single redshift slice is reduced by a factor of ~5 in comparison to a case without secondary non-gaussianities. Setting aside foreground contamination, we forecast that a cosmic-variance-limited experiment observing 21-cm fluctuations over 30 < z < 100 with a 0.1-MHz bandwidth and 0.1-arcminute angular resolution could achieve a sensitivity of order f_NL[local] ~ 0.03, f_NL[equilateral] ~ 0.04, and f_NL[orthogonal] ~ 0.03.

preprint2014arXiv

Search for PAHs in the Perseus molecular cloud with the Green Bank Telescope

Polycyclic Aromatic Hydrocarbons (PAHs) are believed to be the small-size tail of the interstellar carbonaceous dust grain population. Their vibrational emission is the most widely accepted source of the aromatic near-infrared features, and their rotational radiation is a likely explanation for the dust-correlated anomalous microwave emission (AME). Yet, no individual interstellar PAH molecule has been identified to date. It was recently recognized that quasi-symmetric planar PAHs ought to have a well identifiable comb-like rotational spectrum, and suggested to search for them in spectroscopic data with matched-filtering techniques. We report the results of the first such search, carried out with the Green Bank Telescope, and targeting the star-forming region IC348 in the Perseus molecular cloud, a known source of AME. Our observations amounted to 16.75 hours and spanned a 3 GHz-wide band extending from 23.3 to 26.3 GHz. Using frequency switching, we achieved a sensitivity of 0.4 mJy per 0.4 MHz channel (1-sigma). The non-detection of comb-like spectra allowed us to set upper bounds on the abundance of specific quasi-symmetric PAH molecules (specified uniquely by their moments of inertia) of approximately 0.1% of the total PAH abundance. This bound generically applies to PAHs with approximately 15 to 100 carbon atoms.

preprint2013arXiv

A new light on 21 cm intensity fluctuations from the dark ages

Fluctuations of the 21 cm brightness temperature before the formation of the first stars hold the promise of becoming a high-precision cosmological probe in the future. The growth of over densities is very well described by perturbation theory at that epoch and the signal can in principle be predicted to arbitrary accuracy for given cosmological parameters. Recently, Tseliakhovich and Hirata pointed out a previously neglected and important physical effect, due to the fact that baryons and cold dark matter (CDM) have supersonic relative velocities after recombination. This relative velocity suppresses the growth of matter fluctuations on scales $k \sim 10-10^3$ Mpc$^{-1}$. In addition, the amplitude of the small-scale power spectrum is modulated on the large scales over which the relative velocity varies, corresponding to $k \sim 0.005-1$ Mpc$^{-1}$. In this paper, the effect of the relative velocity on 21 cm brightness temperature fluctuations from redshifts $z \geq 30$ is computed. We show that the 21 cm power spectrum is affected on ${\it most}$ scales. On small scales, the signal is typically suppressed several tens of percent, except for extremely small scales ($k \gtrsim 2000$ Mpc$^{-1}$) for which the fluctuations are boosted by resonant excitation of acoustic waves. On large scales, 21 cm fluctuations are enhanced due to the non-linear dependence of the brightness temperature on the underlying gas density and temperature. The enhancement of the 21 cm power spectrum is of a few percent at $k \sim 0.1$ Mpc$^{-1}$ and up to tens of percent at $k \lesssim 0.005$ Mpc$^{-1}$, for standard $Λ$CDM cosmology. In principle this effect allows to probe the small-scale matter power spectrum not only through a measurement of small angular scales but also through its effect on large angular scales.

preprint2013arXiv

Rotational spectroscopy of interstellar PAHs

Polycyclic aromatic hydrocarbons (PAHs) are believed to be ubiquitous in the interstellar medium. Yet, to date no specific PAH molecule has been identified. In this paper, a new observational avenue is suggested to detect individual PAHs, using their rotational line emission at radio frequencies. Previous PAH searches based on rotational spectroscopy have only targeted the bowl-shaped corannulene molecule, with the underlying assumption that other polar PAHs are triaxial and have a complex and diluted spectrum unusable for identification purposes. In this paper the rotational spectrum of quasi-symmetric PAHs is computed analytically. It is shown that the asymmetry of planar, nitrogen-substituted symmetric PAHs is small enough that their rotational spectrum, when observed with a resolution of about a MHz, has the appearance of a "comb" of evenly spaced stacks of lines. The simple pattern of these "comb" spectra allows for the use of matched-filtering techniques, which can result in a significantly enhanced signal-to-noise ratio. Detection forecasts are discussed for regions harbouring "anomalous microwave emission", believed to originate from the collective PAH rotational emission. A systematic search for PAH lines in various environments is advocated. If detected, PAH "combs" would allow to the conclusive and unambiguous identification of specific, free-floating interstellar PAHs.

preprint2012arXiv

An efficient implementation of massive neutrinos in non-linear structure formation simulations

Massive neutrinos make up a fraction of the dark matter, but due to their large thermal velocities, cluster significantly less than cold dark matter (CDM) on small scales. An accurate theoretical modelling of their effect during the non-linear regime of structure formation is required in order to properly analyse current and upcoming high-precision large-scale structure data, and constrain the neutrino mass. Taking advantage of the fact that massive neutrinos remain linearly clustered up to late times, this paper treats the linear growth of neutrino overdensities in a non-linear CDM background. The evolution of the CDM component is obtained via N-body computations. The smooth neutrino component is evaluated from that background by solving the Boltzmann equation linearised with respect to the neutrino overdensity. CDM and neutrinos are simultaneously evolved in time, consistently accounting for their mutual gravitational influence. This method avoids the issue of shot-noise inherent to particle-based neutrino simulations, and, in contrast with standard Fourier-space methods, properly accounts for the non-linear potential wells in which the neutrinos evolve. Inside the most massive late-time clusters, where the escape velocity is larger than the neutrino thermal velocity, neutrinos can clump non-linearly, causing the method to formally break down. It is shown that this does not affect the total matter power spectrum, which can be very accurately computed on all relevant scales up to the present time.

preprint2012arXiv

CMB photons shedding light on dark matter

The annihilation or decay of Dark Matter (DM) particles could affect the thermal history of the universe and leave an observable signature in Cosmic Microwave Background (CMB) anisotropies. We update constraints on the annihilation rate of DM particles in the smooth cosmological background, using WMAP7 and recent small-scale CMB data. With a systematic analysis based on the Press-Schechter formalism, we also show that DM annihilation in halos at small redshift may explain entirely the reionization patterns observed in the CMB, under reasonable assumptions concerning the concentration and formation redshift of halos. We find that a mixed reionization model based on DM annihilation in halos as well as star formation at a redshift z~6.5 could simultaneously account for CMB observations and satisfy constraints inferred from the Gunn-Peterson effect. However, these models tend to reheat the inter-galactic medium (IGM) well above observational bounds: by including a realistic prior on the IGM temperature at low redshift, we find stronger cosmological bounds on the annihilation cross-section than with the CMB alone.

preprint2012arXiv

Effective conductance method for the primordial recombination spectrum

As atoms formed for the first time during primordial recombination, they emitted bound-bound and free-bound radiation leading to spectral distortions to the cosmic microwave background. These distortions might become observable in the future with high-sensitivity spectrometers, and provide a new window into physical conditions in the early universe. The standard multilevel atom method habitually used to compute the recombination spectrum is computationally expensive, impeding a detailed quantitative exploration of the information contained in spectral distortions thus far. In this work it is shown that the emissivity in optically thin allowed transitions can be factored into a computationally expensive but cosmology-independent part and a computationally cheap, cosmology-dependent part. The slow part of the computation consists in pre-computing temperature-dependent effective "conductances", linearly relating line or continuum intensity to departures from Saha equilibrium of the lowest-order excited states (2s and 2p), that can be seen as "voltages". The computation of these departures from equilibrium as a function of redshift is itself very fast, thanks to the effective multilevel atom method introduced in an earlier work. With this factorization, the recurring cost of a single computation of the recombination spectrum is only a fraction of a second on a standard laptop, more than four orders of magnitude shorter than standard computations. The spectrum from helium recombination can be efficiently computed in an identical way, and a fast code computing the full primordial recombination spectrum with this method will be made publicly available soon.

preprint2012arXiv

Spinning dust radiation: a review of the theory

This article reviews the current status of theoretical modeling of electric dipole radiation from spinning dust grains. The fundamentally simple problem of dust grain rotation appeals to a rich set of concepts of classical and quantum physics, owing the the diversity of processes involved. Rotational excitation and damping rates through various mechanisms are discussed, as well as methods of computing the grain angular momentum distribution function. Assumptions on grain properties are reviewed. The robustness of theoretical predictions now seems mostly limited by the uncertainties regarding the grains themselves, namely their abundance, dipole moments, size and shape distribution.

preprint2011arXiv

HyRec: A fast and highly accurate primordial hydrogen and helium recombination code

We present a state-of-the-art primordial recombination code, HyRec, including all the physical effects that have been shown to significantly affect recombination. The computation of helium recombination includes simple analytic treatments of hydrogen continuum opacity in the He I 2 1P - 1 1S line, the He I] 2 3P - 1 1S line, and treats feedback between these lines within the on-the-spot approximation. Hydrogen recombination is computed using the effective multilevel atom method, virtually accounting for an infinite number of excited states. We account for two-photon transitions from 2s and higher levels as well as frequency diffusion in Lyman-alpha with a full radiative transfer calculation. We present a new method to evolve the radiation field simultaneously with the level populations and the free electron fraction. These computations are sped up by taking advantage of the particular sparseness pattern of the equations describing the radiative transfer. The computation time for a full recombination history is ~2 seconds. This makes our code well suited for inclusion in Monte Carlo Markov chains for cosmological parameter estimation from upcoming high-precision cosmic microwave background anisotropy measurements.

preprint2011arXiv

Metals at the surface of last scatter

Standard big-bang nucleosynthesis (BBN) predicts only a trace abundance of lithium and no heavier elements, but some alternatives predict a nonzero primordial metallicity. Here we explore whether CMB measurements may set useful constraints to the primordial metallicity and/or whether the standard CMB calculations are robust, within the tolerance of forthcoming CMB maps, to the possibility of primordial metals. Metals would affect the recombination history (and thus CMB power spectra) in three ways: (1) Lyman-alpha photons can be removed (and recombination thus accelerated) by photoionizing metals. (2) The Bowen resonance-fluorescence mechanism may degrade Lyman-beta photons and thus enhance the Lyman-beta escape probability and speed up recombination. (3) Metals could affect the low-redshift tail of the CMB visibility function by providing additional free electrons. The last two of these provide the strongest CMB signal. However, the effects are detectable in the Planck satellite only if the primordial metal abundance is at least a few hundredths of solar for (2) and a few tenths of solar for (3). We thus conclude that Planck will not be able to improve upon current constraints to primordial metallicity, at the level of a thousandth of solar, from the Lyman-alpha forest and ultra-metal-poor halo stars, and that the CMB power-spectrum predictions for Planck suffer no uncertainty arising from the possibility that there may be primordial metals.

preprint2011arXiv

Revisiting the double-binary-pulsar probe of non-dynamical Chern-Simons gravity

One of the popular modifications to the theory of general relativity is non-dynamical Chern-Simons (CS) gravity, in which the metric is coupled to an externally prescribed scalar field. Setting accurate constraints to the parameters of the theory is important owing to their implications for the scalar field and/or the underlying fundamental theory. The current best constraints rely on measurements of the periastron precession rate in the double-binary-pulsar system and place a very tight bound on the characteristic CS lengthscale k_cs^{-1} <~ 3*10^{-9} km. This paper considers several effects that were not accounted for when deriving this bound and lead to a substantial suppression of the predicted rate of periastron precession. It is shown, in particular, that the point mass approximation for extended test bodies does not apply in this case. The constraint to the characteristic CS lengthscale is revised to k_cs^{-1} <~ 0.4 km, eight orders of magnitude weaker than what was previously found.

preprint2011arXiv

Slowly-rotating stars and black holes in dynamical Chern-Simons gravity

Chern-Simons (CS) modified gravity is an extension to general relativity (GR) in which the metric is coupled to a scalar field, resulting in modified Einstein field equations. In the dynamical theory, the scalar field is itself sourced by the Pontryagin density of the space-time. In this paper, the coupled system of equations for the metric and the scalar field is solved numerically for slowly-rotating neutron stars described with realistic equations of state and for slowly-rotating black holes. An analytic solution for a constant-density nonrelativistic object is also presented. It is shown that the black hole solution cannot be used to describe the exterior spacetime of a star as was previously assumed. In addition, whereas previous analysis were limited to the small-coupling regime, this paper considers arbitrarily large coupling strengths. It is found that the CS modification leads to two effects on the gravitomagnetic sector of the metric: (i) Near the surface of a star or the horizon of a black hole, the magnitude of the gravitomagnetic potential is decreased and frame-dragging effects are reduced in comparison to GR. (ii) In the case of a star, the angular momentum J, as measured by distant observers, is enhanced in CS gravity as compared to standard GR. For a large coupling strength, the near-zone frame-dragging effects become significantly screened, whereas the far-zone enhancement saturate at a maximum value max(Delta J) ~ (M/R) J. Using measurements of frame-dragging effects around the Earth by Gravity Probe B and the LAGEOS satellites, a weak but robust constraint is set to the characteristic CS lengthscale, xi^{1/4} <~ 10^8 km.

preprint2010arXiv

Radiative transfer effects in primordial hydrogen recombination

The calculation of a highly accurate cosmological recombination history has been the object of particular attention recently, as it constitutes the major theoretical uncertainty when predicting the angular power spectrum of Cosmic Microwave Background anisotropies. Lyman transitions, in particular the Lyman-alpha line, have long been recognized as one of the bottlenecks of recombination, due to their very low escape probabilities. The Sobolev approximation does not describe radiative transfer in the vicinity of Lyman lines to a sufficient degree of accuracy, and several corrections have already been computed in other works. In this paper, the impact of some previously ignored radiative transfer effects is calculated. First, the effect of Thomson scattering in the vicinity of the Lyman-alpha line is evaluated, using a full redistribution kernel incorporated into a radiative transfer code. The effect of feedback of distortions generated by the optically thick deuterium Lyman-alpha line blueward of the hydrogen line is investigated with an analytic approximation. It is shown that both effects are negligible during cosmological hydrogen recombination. Secondly, the importance of high-lying, non overlapping Lyman transitions is assessed. It is shown that escape from lines above Ly-gamma and frequency diffusion in Ly-beta and higher lines can be neglected without loss of accuracy. Thirdly, a formalism generalizing the Sobolev approximation is developed to account for the overlap of the high-lying Lyman lines, which is shown to lead to negligible changes to the recombination history. Finally, the possibility of a cosmological hydrogen recombination maser is investigated. It is shown that there is no such maser in the purely radiative treatment presented here.

preprint2010arXiv

Ultrafast effective multi-level atom method for primordial hydrogen recombination

Cosmological hydrogen recombination has recently been the subject of renewed attention because of its importance for predicting the power spectrum of cosmic microwave background anisotropies. It has become clear that it is necessary to account for a large number n >~ 100 of energy shells of the hydrogen atom, separately following the angular momentum substates in order to obtain sufficiently accurate recombination histories. However, the multi-level atom codes that follow the populations of all these levels are computationally expensive, limiting recent analyses to only a few points in parameter space. In this paper, we present a new method for solving the multi-level atom recombination problem, which splits the problem into a computationally expensive atomic physics component that is independent of the cosmology, and an ultrafast cosmological evolution component. The atomic physics component follows the network of bound-bound and bound-free transitions among excited states and computes the resulting effective transition rates for the small set of "interface" states radiatively connected to the ground state. The cosmological evolution component only follows the populations of the interface states. By pre-tabulating the effective rates, we can reduce the recurring cost of multi-level atom calculations by more than 5 orders of magnitude. The resulting code is fast enough for inclusion in Markov Chain Monte Carlo parameter estimation algorithms. It does not yet include the radiative transfer or high-n two-photon processes considered in some recent papers. Further work on analytic treatments for these effects will be required in order to produce a recombination code usable for Planck data analysis.