Showing 1–50 of 74 results for author: Mohayaee, R

Monge-Ampère gravity, optimal transport theory and their link to the Galileons

Authors: Albert Bonnefous, Yann Brenier, Roya Mohayaee

Abstract: Mathematicians have been proposing for sometimes that Monge-Ampère equation, a nonlinear generalization of the Poisson equation, where trace of the Hessian is replaced by its determinant, provides an alternative non-relativistic description of gravity. Monge-Ampère equation is affine invariant, has rich geometric properties, connects to optimal transport theory, and remains bounded at short distan… ▽ More Mathematicians have been proposing for sometimes that Monge-Ampère equation, a nonlinear generalization of the Poisson equation, where trace of the Hessian is replaced by its determinant, provides an alternative non-relativistic description of gravity. Monge-Ampère equation is affine invariant, has rich geometric properties, connects to optimal transport theory, and remains bounded at short distances. Monge-Ampère gravity, that uses a slightly different form of the Monge-Ampère equation, naturally emerges through the application of large-deviation principle to a Brownian system of indistinguishable and independent particles. In this work we provide a physical formulation of this mathematical model, study its theoretical viability and confront it with observations. We show that Monge-Ampère gravity cannot replace the Newtonian gravity as it does not withstand the solar-system test. We then show that Monge-Ampère gravity can describe a scalar field, often evoked in modified theories of gravity such as Galileons. We show that Monge-Ampère gravity, as a nonlinear model of a new scalar field, is screened at short distances, and behaves differently from Newtonian gravity above galactic scales but approaches it asymptotically. Finally, we write a relativistic Lagrangian for Monge-Ampère gravity in flat space time, which is the field equation of a sum of the Lagrangians of all Galileons. We also show how the Monge-Ampère equation can be obtained from the fully covariant Lagrangian of quartic Galileon in the static limit. The connection between optimal transport theory and modified theories of gravity with second-order field equations, unravelled here, remains a promising domain to further explore. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Comments: 16 pages, 3 figures

Monge Ampère gravity: from the large deviation principle to cosmological simulations through optimal transport

Authors: Bruno Lévy, Yann Brenier, Roya Mohayaee

Abstract: We study Monge-Ampère gravity (MAG) as an effective theory of cosmological structure formation through optimal transport theory. MAG is based on the Monge-Ampère equation, a nonlinear version of the Poisson equation, that relates the Hessian determinant of the potential to the density field. We explain how MAG emerges from a conditioned system of independent and indistinguishable Brownian particle… ▽ More We study Monge-Ampère gravity (MAG) as an effective theory of cosmological structure formation through optimal transport theory. MAG is based on the Monge-Ampère equation, a nonlinear version of the Poisson equation, that relates the Hessian determinant of the potential to the density field. We explain how MAG emerges from a conditioned system of independent and indistinguishable Brownian particles, through the large deviation principle, in the continuum limit. To numerically explore this highly non-linear theory, we develop a novel N-body simulation method based on semi-discrete optimal transport. Our results obtained from the very first N-body simulation of Monge-Ampère gravity with over 100 millions particles show that on large scales, Monge-Ampère gravity is similar to the Newtonian gravity but favours the formation of anisotropic structures such as filaments. At small scales, MAG has a weaker clustering and is screened in high-density regions. Although here we study the Monge-Ampère gravity as an effective rather than a fundamental theory, our novel highly-performant optimal transport algorithm can be used to run high-resolution simulations of a large class of modified theories of gravity, such as Galileons, in which the equations of motion are second-order and of Monge-Amp��re type. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: 15 pages, 9 figures

Displacement Field Analysis via Optimal Transport: Multi-Tracer Approach to Cosmological Reconstruction

Authors: Farnik Nikakhtar, Ravi K. Sheth, Nikhil Padmanabhan, Bruno Lévy, Roya Mohayaee

Abstract: We demonstrate the effectiveness of one of the many multi-tracer analyses enabled by Optimal Transport (OT) reconstruction. Leveraging a semi-discrete OT algorithm, we determine the displacements between initial and observed positions of biased tracers and the remaining matter field. With only redshift-space distorted final positions of biased tracers and a simple premise for the remaining mass di… ▽ More We demonstrate the effectiveness of one of the many multi-tracer analyses enabled by Optimal Transport (OT) reconstruction. Leveraging a semi-discrete OT algorithm, we determine the displacements between initial and observed positions of biased tracers and the remaining matter field. With only redshift-space distorted final positions of biased tracers and a simple premise for the remaining mass distribution as input, OT solves the displacement field. This extracted field, assuming asymptotically uniform density and a gradient flow displacement, enables reconstruction of the initial overdensity fluctuation field. We show that the divergence of the OT displacement field is a good proxy of the linear density field, even though the method never assumes the linear theory growth. Additionally, this divergence field can be combined with the reconstructed protohalos to provide a higher signal-to-noise measurement of the BAO standard ruler than was possible with either measurement individually. △ Less

Submitted 18 March, 2024; originally announced March 2024.

Comments: 12 pages, 10 figures

Optimal Transport Reconstruction of Biased Tracers in Redshift Space

Authors: Farnik Nikakhtar, Nikhil Padmanabhan, Bruno Lévy, Ravi K. Sheth, Roya Mohayaee

Abstract: Recent research has emphasized the benefits of accurately reconstructing the initial Lagrangian positions of biased tracers from their positions at a later time, to gain cosmological information. A weighted semi-discrete optimal transport algorithm can achieve the required accuracy, provided the late-time positions are known, with minimal information about the background cosmology. The algorithm's… ▽ More Recent research has emphasized the benefits of accurately reconstructing the initial Lagrangian positions of biased tracers from their positions at a later time, to gain cosmological information. A weighted semi-discrete optimal transport algorithm can achieve the required accuracy, provided the late-time positions are known, with minimal information about the background cosmology. The algorithm's performance relies on knowing the masses of the biased tracers, and depends on how one models the distribution of the remaining mass that is not associated with these tracers. We demonstrate that simple models of the remaining mass result in accurate retrieval of the initial Lagrangian positions, which we quantify using pair statistics and the void probability function. This is true even if the input positions are affected by redshift-space distortions. The most sophisticated models assume that the masses of the tracers, and the amount and clustering of the missing mass are known; we show that the method is robust to realistic errors in the masses of the tracers and remains so as the model for the missing mass becomes increasingly crude. △ Less

Submitted 7 July, 2023; originally announced July 2023.

Comments: 15 pages, 18 figures

Distortions of the Hubble diagram: Line-of-sight signatures of local galaxy clusters

Authors: Jenny G. Sorce, Roya Mohayaee, Nabila Aghanim, Klaus Dolag, Nicola Malavasi

Abstract: The Universe expansion rate is modulated around local inhomogeneities due to their gravitational potential. Velocity waves are then observed around galaxy clusters in the Hubble diagram. This paper studies them in a ~738 Mpc wide, with 2048^3 particles, cosmological simulation of our cosmic environment (a.k.a. CLONE: Constrained LOcal \& Nesting Environment simulation). For the first time, the sim… ▽ More The Universe expansion rate is modulated around local inhomogeneities due to their gravitational potential. Velocity waves are then observed around galaxy clusters in the Hubble diagram. This paper studies them in a ~738 Mpc wide, with 2048^3 particles, cosmological simulation of our cosmic environment (a.k.a. CLONE: Constrained LOcal \& Nesting Environment simulation). For the first time, the simulation shows that velocity waves that arise in the lines-of-sight of the most massive dark matter halos agree with those observed in local galaxy velocity catalogs in the lines-of-sight of Coma and several other local (Abell) clusters. For the best-constrained clusters such as Virgo and Centaurus, i.e. those closest to us, secondary waves caused by galaxy groups, further into the non-linear regime, also stand out. This match is not utterly expected given that before being evolved into a fully non-linear z=0 state, assuming ΛCDM, CLONE initial conditions are constrained with solely linear theory, power spectrum and highly uncertain and sparse local peculiar velocities. Additionally, Gaussian fits to velocity wave envelopes show that wave properties are tightly tangled with cluster masses. This link is complex though and involves the environment and formation history of the clusters. A proposed metric, measuring the distance between the observed and several re-centered simulated lines-of-sight, waves included, is shown to be capable of providing a tight mass range estimate for massive local clusters. Using machine learning techniques to grasp more thoroughly the complex wave-mass relation, velocity waves could in the near future be used to provide additional and independent mass estimates from galaxy dynamics within large cluster radii. △ Less

Submitted 29 April, 2024; v1 submitted 3 January, 2023; originally announced January 2023.

Comments: Accepted for publications in A&A, 17 pages, 12 figures and 5 tables

Galaxy number-count dipole and superhorizon fluctuations

Authors: Guillem Domènech, Roya Mohayaee, Subodh P. Patil, Subir Sarkar

Abstract: In view of the growing tension between the dipole anisotropy of number counts of cosmologically distant sources and of the cosmic microwave background (CMB), we investigate the number count dipole induced by primordial perturbations with wavelength comparable to or exceeding the Hubble radius today. First, we find that neither adiabatic nor isocurvature superhorizon modes can generate an intrinsic… ▽ More In view of the growing tension between the dipole anisotropy of number counts of cosmologically distant sources and of the cosmic microwave background (CMB), we investigate the number count dipole induced by primordial perturbations with wavelength comparable to or exceeding the Hubble radius today. First, we find that neither adiabatic nor isocurvature superhorizon modes can generate an intrinsic number count dipole. However a superhorizon isocurvature mode does induce a relative velocity between the CMB and the (dark) matter rest frames and thereby affects the CMB dipole. We revisit the possibility that it has an intrinsic component due to such a mode, thus enabling consistency with the galaxy number count dipole if the latter is actually kinematic in origin. Although this scenario is not particularly natural, there are possible links with other anomalies and it predicts a concommitant galaxy number count quadrupole which may be measurable in future surveys. We also investigate the number count dipole induced by modes smaller than the Hubble radius, finding that subject to CMB constraints this is too small to reconcile the dipole tension. △ Less

Submitted 4 July, 2022; originally announced July 2022.

Comments: 25 pages + 15 figures + 1 table + 8 pages of appendices

Journal ref: JCAP 10 (2022) 019

A Challenge to the Standard Cosmological Model

Authors: Nathan Secrest, Sebastian von Hausegger, Mohamed Rameez, Roya Mohayaee, Subir Sarkar

Abstract: We present the first joint analysis of catalogs of radio galaxies and quasars to determine if their sky distribution is consistent with the standard $Λ$CDM model of cosmology. This model is based on the cosmological principle, which asserts that the universe is statistically isotropic and homogeneous on large scales, so the observed dipole anisotropy in the cosmic microwave background (CMB) must b… ▽ More We present the first joint analysis of catalogs of radio galaxies and quasars to determine if their sky distribution is consistent with the standard $Λ$CDM model of cosmology. This model is based on the cosmological principle, which asserts that the universe is statistically isotropic and homogeneous on large scales, so the observed dipole anisotropy in the cosmic microwave background (CMB) must be attributed to our local peculiar motion. We test the null hypothesis that there is a dipole anisotropy in the sky distribution of radio galaxies and quasars consistent with the motion inferred from the CMB, as is expected for cosmologically distant sources. Our two samples, constructed respectively from the NRAO VLA Sky Survey and the Wide-field Infrared Survey Explorer, are systematically independent and have no shared objects. Using a completely general statistic that accounts for correlation between the found dipole amplitude and its directional offset from the CMB dipole, the null hypothesis is independently rejected by the radio galaxy and quasar samples with $p$-value of $8.9\times10^{-3}$ and $1.2\times10^{-5}$, respectively, corresponding to $2.6σ$ and $4.4σ$ significance. The joint significance, using sample size-weighted $Z$-scores, is $5.1σ$. We show that the radio galaxy and quasar dipoles are consistent with each other and find no evidence for any frequency dependence of the amplitude. The consistency of the two dipoles improves if we boost to the CMB frame assuming its dipole to be fully kinematic, suggesting that cosmologically distant radio galaxies and quasars may have an intrinsic anisotropy in this frame. △ Less

Submitted 12 August, 2022; v1 submitted 11 June, 2022; originally announced June 2022.

Comments: 11 pages, 4 figures, accepted for publication in ApJ Letters. Code and data available at https://zenodo.org/record/6784602

Journal ref: Astrophys. J. Lett. 937 (2022) L31

Optimal Transport Reconstruction of Baryon Acoustic Oscillations

Authors: Farnik Nikakhtar, Ravi K. Sheth, Bruno Lévy, Roya Mohayaee

Abstract: A weighted, semi-discrete, fast optimal transport (OT) algorithm for reconstructing the Lagrangian positions of proto-halos from their evolved Eulerian positions is presented. The algorithm makes use of a mass estimate of the biased tracers and of the distribution of the remaining mass (the `dust'), but is robust to errors in the mass estimates. Tests with state-of-art cosmological simulations sho… ▽ More A weighted, semi-discrete, fast optimal transport (OT) algorithm for reconstructing the Lagrangian positions of proto-halos from their evolved Eulerian positions is presented. The algorithm makes use of a mass estimate of the biased tracers and of the distribution of the remaining mass (the `dust'), but is robust to errors in the mass estimates. Tests with state-of-art cosmological simulations show that if the dust is assumed to have a uniform spatial distribution, then the shape of the OT-reconstructed pair correlation function of the tracers is very close to linear theory, enabling sub-percent precision in the BAO distance scale that depends weakly, if at all, on a cosmological model. With a more sophisticated model for the dust, OT returns an estimate of the displacement field which yields superb reconstruction of the proto-halo positions, and hence of the shape and amplitude of the initial pair correlation function of the tracers. This enables direct and independent determinations of the bias factor $b$ and the smearing scale $Σ$, potentially providing new methods for breaking the degeneracy between $b$ and $σ_8$. △ Less

Submitted 9 January, 2023; v1 submitted 3 March, 2022; originally announced March 2022.

Comments: matches published version (Physical Review Letters)

Accurate Baryon Acoustic Oscillations reconstruction via semi-discrete optimal transport

Authors: Sebastian von Hausegger, Bruno Lévy, Roya Mohayaee

Abstract: Optimal transport theory has recently reemerged as a vastly resourceful field of mathematics with elegant applications across physics and computer science. Harnessing methods from geometry processing, we report on the efficient implementation for a specific problem in cosmology -- the reconstruction of the linear density field from low redshifts, in particular the recovery of the Baryonic Acoustic… ▽ More Optimal transport theory has recently reemerged as a vastly resourceful field of mathematics with elegant applications across physics and computer science. Harnessing methods from geometry processing, we report on the efficient implementation for a specific problem in cosmology -- the reconstruction of the linear density field from low redshifts, in particular the recovery of the Baryonic Acoustic Oscillation (BAO) scale. We demonstrate our algorithm's accuracy by retrieving the BAO scale in noise-less cosmological simulations that are dedicated to cancel cosmic variance; we find uncertainties to be reduced by a factor of 4.3 compared with performing no reconstruction, and a factor of 3.1 compared with standard reconstruction. △ Less

Submitted 6 April, 2022; v1 submitted 17 October, 2021; originally announced October 2021.

Comments: Accepted for publication in PRL. 7 pages incl. references; 2 figures; 1 table. v1 initial submission. v2 fixed affiliation. v3 matches published version, but larger figures

Do supernovae indicate an accelerating universe?

Authors: Roya Mohayaee, Mohamed Rameez, Subir Sarkar

Abstract: In the late 1990's, observations of 93 Type Ia supernovae were analysed in the framework of the FLRW cosmology assuming these to be `standard(isable) candles'. It was thus inferred that the Hubble expansion rate is accelerating as if driven by a positive Cosmological Constant $Λ$. This is still the only direct evidence for the `dark energy' that is the dominant component of the standard $Λ$CDM cos… ▽ More In the late 1990's, observations of 93 Type Ia supernovae were analysed in the framework of the FLRW cosmology assuming these to be `standard(isable) candles'. It was thus inferred that the Hubble expansion rate is accelerating as if driven by a positive Cosmological Constant $Λ$. This is still the only direct evidence for the `dark energy' that is the dominant component of the standard $Λ$CDM cosmological model. Other data such as BAO, CMB anisotropies, stellar ages, the rate of structure growth, etc are all `concordant' with this model but do not provide independent evidence for accelerated expansion. Analysis of a larger sample of 740 SNe Ia shows that these are not quite standard candles, and highlights the "corrections" applied to analyse the data in the FLRW framework. The latter holds in the reference frame in which the CMB is isotropic, whereas observations are made in our heliocentric frame in which the CMB has a large dipole anisotropy. This is assumed to be of kinematic origin i.e. due to our non-Hubble motion driven by local inhomogeneity in the matter distribution. The $Λ$CDM model predicts how this peculiar velocity should fall off as the averaging scale is raised and the universe becomes sensibly homogeneous. However observations of the local `bulk flow' are inconsistent with this expectation and convergence to the CMB frame is not seen. Moreover the kinematic interpretation implies a corresponding dipole in the sky distribution of high redshift quasars, which is rejected by observations at 4.9$σ$. The acceleration of the Hubble expansion rate is also anisotropic at 3.9$σ$ and aligned with the bulk flow. Thus dark energy may be an artefact of analysing data assuming that we are idealised observers in an FLRW universe, when in fact the real universe is inhomogeneous and anisotropic out to distances large enough to impact on cosmological analyses. △ Less

Submitted 6 June, 2021; originally announced June 2021.

Comments: 14 pages, 4 figures, 4 tables

Journal ref: Eur. Phys. J. Spec. Top. 230(9), 2067-2076 (2021)

Hints of FLRW Breakdown from Supernovae

Authors: Chethan Krishnan, Roya Mohayaee, Eoin Ó Colgáin, M. M. Sheikh-Jabbari, Lu Yin

Abstract: A 10\% difference in the scale for the Hubble parameter constitutes a clear problem for cosmology. Here, considering angular distribution of Type Ia supernovae (SN) within the Pantheon compilation and working within flat $Λ$CDM cosmology, we observe a correlation between higher $H_0$ and the CMB dipole direction, confirming our previous results for strongly-lensed quasars \cite{Krishnan:2021dyb}.… ▽ More A 10\% difference in the scale for the Hubble parameter constitutes a clear problem for cosmology. Here, considering angular distribution of Type Ia supernovae (SN) within the Pantheon compilation and working within flat $Λ$CDM cosmology, we observe a correlation between higher $H_0$ and the CMB dipole direction, confirming our previous results for strongly-lensed quasars \cite{Krishnan:2021dyb}. Concretely, we record a $\sim 1$ km/s/Mpc variation in $H_0$ at antipodal points on the sky within the Pantheon sample, which is evident in the Low $z$ subsample ($z \lesssim 0.075$) and gets enhanced by higher redshift SN. Our work raises the possibility that we may be at the precision required to probe anisotropic Hubble expansions, while providing a concrete prediction for future inferences of $H_0$. △ Less

Submitted 23 February, 2022; v1 submitted 4 June, 2021; originally announced June 2021.

Comments: 5 pages, 5 figures, comments are most welcome; v2 appendix added on general trend; v3 presentation improved and streamlined; v4 references and appendix added, to appear in PRD

Does Hubble Tension Signal a Breakdown in FLRW Cosmology?

Authors: Chethan Krishnan, Roya Mohayaee, Eoin Ó Colgáin, M. M. Sheikh-Jabbari, Lu Yin

Abstract: The tension between early and late Universe probes of the Hubble constant has motivated various new FLRW cosmologies. Here, we reanalyse the Hubble tension with a recent age of the Universe constraint. This allows us to restrict attention to matter and a dark energy sector that we treat without assuming a specific model. Assuming analyticity of the Hubble parameter $H(z)$, and a generic low redshi… ▽ More The tension between early and late Universe probes of the Hubble constant has motivated various new FLRW cosmologies. Here, we reanalyse the Hubble tension with a recent age of the Universe constraint. This allows us to restrict attention to matter and a dark energy sector that we treat without assuming a specific model. Assuming analyticity of the Hubble parameter $H(z)$, and a generic low redshift modification to flat $Λ$CDM, we find that low redshift data ($z \lesssim 2.5$) and well-motivated priors only permit a dark energy sector close to the cosmological constant $Λ$. This restriction rules out late Universe modifications within FLRW. We show that early Universe physics that alters the sound horizon can yield an upper limit of $H_0 \sim 71 \pm 1$ km/s/Mpc. Since various local determinations may be converging to $H_0 \sim 73$ km/s/Mpc, a breakdown of the FLRW framework is a plausible resolution. We outline how future data, in particular strongly lensed quasar data, could also provide further confirmations of such a resolution. △ Less

Submitted 23 August, 2021; v1 submitted 20 May, 2021; originally announced May 2021.

Comments: v1 10 pages, comments welcome, especially on appendix D & Figure 4; v2, typos corrected, Figure 4 improved, observations strengthened; v3 references added, matches version published as invited contribution to CQG "Focus Issue on the Hubble Constant Tension"

Journal ref: 2021 Class. Quantum Grav. 38 184001

A fast semi-discrete optimal transport algorithm for a unique reconstruction of the early Universe

Authors: Bruno Lévy, Roya Mohayaee, Sebastian von Hausegger

Abstract: We leverage powerful mathematical tools stemming from optimal transport theory and transform them into an efficient algorithm to reconstruct the fluctuations of the primordial density field, built on solving the Monge-Ampère-Kantorovich equation. Our algorithm computes the optimal transport between an initial uniform continuous density field, partitioned into Laguerre cells, and a final input set… ▽ More We leverage powerful mathematical tools stemming from optimal transport theory and transform them into an efficient algorithm to reconstruct the fluctuations of the primordial density field, built on solving the Monge-Ampère-Kantorovich equation. Our algorithm computes the optimal transport between an initial uniform continuous density field, partitioned into Laguerre cells, and a final input set of discrete point masses, linking the early to the late Universe. While existing early universe reconstruction algorithms based on fully discrete combinatorial methods are limited to a few hundred thousand points, our algorithm scales up well beyond this limit, since it takes the form of a well-posed smooth convex optimization problem, solved using a Newton method. We run our algorithm on cosmological $N$-body simulations, from the AbacusCosmos suite, and reconstruct the initial positions of $\mathcal{O}(10^7)$ particles within a few hours with an off-the-shelf personal computer. We show that our method allows a unique, fast and precise recovery of subtle features of the initial power spectrum, such as the baryonic acoustic oscillations. △ Less

Submitted 16 December, 2020; originally announced December 2020.

Comments: 22 pages

A Test of the Cosmological Principle with Quasars

Authors: Nathan Secrest, Sebastian von Hausegger, Mohamed Rameez, Roya Mohayaee, Subir Sarkar, Jacques Colin

Abstract: We study the large-scale anisotropy of the Universe by measuring the dipole in the angular distribution of a flux-limited, all-sky sample of 1.36 million quasars observed by the Wide-field Infrared Survey Explorer (WISE). This sample is derived from the new CatWISE2020 catalog, which contains deep photometric measurements at 3.4 and 4.6 $μ$m from the cryogenic, post-cryogenic, and reactivation pha… ▽ More We study the large-scale anisotropy of the Universe by measuring the dipole in the angular distribution of a flux-limited, all-sky sample of 1.36 million quasars observed by the Wide-field Infrared Survey Explorer (WISE). This sample is derived from the new CatWISE2020 catalog, which contains deep photometric measurements at 3.4 and 4.6 $μ$m from the cryogenic, post-cryogenic, and reactivation phases of the WISE mission. While the direction of the dipole in the quasar sky is similar to that of the cosmic microwave background (CMB), its amplitude is over twice as large as expected, rejecting the canonical, exclusively kinematic interpretation of the CMB dipole with a p-value of $5\times10^{-7}$ ($4.9σ$ for a normal distribution, one-sided), the highest significance achieved to date in such studies. Our results are in conflict with the cosmological principle, a foundational assumption of the concordance $Λ$CDM model. △ Less

Submitted 18 January, 2021; v1 submitted 30 September, 2020; originally announced September 2020.

Comments: 8 pages, 4 figures, accepted for publication in ApJ Letters. Code and data available at https://doi.org/10.5281/zenodo.4431089

Journal ref: Astrophys. J. Lett. 908 (2021) L51

Cosmological inference from within the peculiar local universe

Authors: Roya Mohayaee, Mohamed Rameez, Subir Sarkar

Abstract: The existence of 'peculiar' velocities due to the formation of cosmic structure marks a point of discord between the real Universe and the usually assumed Friedmann-Lemaítre-Robertson-Walker metric which accomodates only the smooth Hubble expansion on large scales. In the standard $Λ$CDM model framework, Type Ia supernovae data are routinely "corrected" for the peculiar velocities of both the obse… ▽ More The existence of 'peculiar' velocities due to the formation of cosmic structure marks a point of discord between the real Universe and the usually assumed Friedmann-Lemaítre-Robertson-Walker metric which accomodates only the smooth Hubble expansion on large scales. In the standard $Λ$CDM model framework, Type Ia supernovae data are routinely "corrected" for the peculiar velocities of both the observer and the supernova host galaxies relative to the cosmic rest frame, in order to infer evidence for acceleration of the expansion rate from their Hubble diagram. However observations indicate a strong, coherent local bulk flow that continues outward without decaying out to a redshift $z \gtrsim 0.1$, contrary to the $Λ$CDM expectation. By querying the halo catalogue of the Dark Sky Hubble-volume N-body simulation, we find that an observer placed in an unusual environment like our local Universe should see correlations between supernovae in the JLA catalogue that are 2-8 times stronger than seen by a typical or Copernican observer. This accounts for our finding that peculiar velocity corrections have a large impact on the value of the Cosmological Constant inferred from supernova data. We also demonstrate that local Universe-like observers will infer a downward biased value of the clustering parameter $S_8$ from comparing the density and velocity fields. More realistic modelling of the peculiar local Universe is thus essential for correctly interpreting cosmological data. △ Less

Submitted 8 May, 2024; v1 submitted 23 March, 2020; originally announced March 2020.

Comments: 15 pages, 8 figures, 1 table

Journal ref: Universe 10 (2024) 209

Subhalo sinking and off-center massive black holes in dwarf galaxies

Authors: Pierre Boldrini, Roya Mohayaee, Joseph Silk

Abstract: Using fully GPU $N$-body simulations, we demonstrate for the first time that subhalos sink and transfer energy via dynamical friction into the centres of dwarf galaxies. This dynamical heating kicks any central massive black hole (MBH) out to tens of parsecs, especially at early epochs ($z$=1.5-3). This mechanism helps explain the observed off-center BHs in dwarf galaxies and also predicts that of… ▽ More Using fully GPU $N$-body simulations, we demonstrate for the first time that subhalos sink and transfer energy via dynamical friction into the centres of dwarf galaxies. This dynamical heating kicks any central massive black hole (MBH) out to tens of parsecs, especially at early epochs ($z$=1.5-3). This mechanism helps explain the observed off-center BHs in dwarf galaxies and also predicts that off-center BHs are more common in higher mass dwarf galaxies since dynamical friction becomes significantly weaker and BHs take more time to sink back towards the centres of their host galaxies. One consequence of off-center BHs during early epochs of dwarf galaxies is to quench any BH feedback. △ Less

Submitted 9 March, 2020; v1 submitted 5 March, 2020; originally announced March 2020.

Comments: 5 pages, 5 figures. Version to appear in MNRAS Letters

Flattening of dark matter cusps during mergers: model of M31

Authors: Pierre Boldrini, Roya Mohayaee, Joseph Silk

Abstract: We run high resolution fully GPU N-body simulations to demonstrate that the dark matter distribution in M31 is well-fitted by a core-like profile. Rich observational data especially on the giant stellar stream provides stringent constraints on the initial conditions of our simulations. We demonstrate that accretion of a satellite on a highly eccentric orbit heats up the central parts of M31, cause… ▽ More We run high resolution fully GPU N-body simulations to demonstrate that the dark matter distribution in M31 is well-fitted by a core-like profile. Rich observational data especially on the giant stellar stream provides stringent constraints on the initial conditions of our simulations. We demonstrate that accretion of a satellite on a highly eccentric orbit heats up the central parts of M31, causes an outward migration of dark matter particles, flattens the central cusp over more than a decade in scale and generates a new model-independent dark matter profile that is well-fitted by a core. Our results imply that core-like central profiles could be a common feature of massive galaxies that have been initially cuspy but have accreted satellites on nearly radial orbits. △ Less

Submitted 13 July, 2021; v1 submitted 27 February, 2020; originally announced February 2020.

Comments: 14 pages, 10 figures. Version to appear in ApJ

A response to Rubin & Heitlauf: "Is the expansion of the universe accelerating? All signs still point to yes"

Authors: Jacques Colin, Roya Mohayaee, Mohamed Rameez, Subir Sarkar

Abstract: We have shown (Colin et al., 2019) that the acceleration of the Hubble expansion rate inferred from Type Ia supernovae (SNe Ia) is, at $3.9σ$ significance, a dipole approximately aligned with the CMB dipole, while its monopole component, which can be interpreted as due to a Cosmological Constant or dark energy, is consistent with zero at $1.4σ$. This has been challenged by Rubin & Heitlauf (2019)… ▽ More We have shown (Colin et al., 2019) that the acceleration of the Hubble expansion rate inferred from Type Ia supernovae (SNe Ia) is, at $3.9σ$ significance, a dipole approximately aligned with the CMB dipole, while its monopole component, which can be interpreted as due to a Cosmological Constant or dark energy, is consistent with zero at $1.4σ$. This has been challenged by Rubin & Heitlauf (2019) who assert that the dipole arises because we made an incorrect assumption about the SNe Ia light-curve parameters (viz. took them to be sample- and redshift independent), and did not allow for the motion of the Solar system (w.r.t. the 'CMB frame' in which the CMB dipole supposedly vanishes). In fact what has an even larger impact on our finding is that we reversed the inconsistent "corrections" made for the peculiar velocities of the SNe Ia host galaxies w.r.t the CMB frame, which in fact serve to bias the lever arm of the Hubble diagram towards higher inferred values of the monopole. We demonstrate that even if all such corrections are made consistently and both sample- and redshift-dependence is allowed for in the standardisation of supernova light curves, the evidence for isotropic acceleration rises to just $2.8\,σ$. Thus the criticism of Rubin & Heitlauf serves only to highlight that "corrections" must be made to the SNe Ia data assuming the standard $Λ$CDM model, in order to recover it from the data. △ Less

Submitted 23 March, 2020; v1 submitted 9 December, 2019; originally announced December 2019.

Comments: 6 pages, 2 figures, 3 tables; Expanded to respond to v.2 of arXiv:1912.02191

Embedding globular clusters in dark matter minihalos solves the cusp-core and timing problems in the Fornax dwarf galaxy

Authors: Pierre Boldrini, Roya Mohayaee, Joseph Silk

Abstract: We use a fully GPU $N$-body code to demonstrate that dark matter minihalos, as a new component of globular clusters, resolve both the timing and cusp-core problems in Fornax if the five (or six) globular clusters were recently accreted ($\leq$ 3 Gyr ago) by Fornax. Under these assumptions, infall of these globular clusters does not occur and no star clusters form in the centre of Fornax in accorda… ▽ More We use a fully GPU $N$-body code to demonstrate that dark matter minihalos, as a new component of globular clusters, resolve both the timing and cusp-core problems in Fornax if the five (or six) globular clusters were recently accreted ($\leq$ 3 Gyr ago) by Fornax. Under these assumptions, infall of these globular clusters does not occur and no star clusters form in the centre of Fornax in accordance with observations. We find that crossings of globular clusters that have DM minihalos near the Fornax centre induce a cusp-to-core transition of the dark matter halo and hence resolve the cusp-core problem in this dwarf galaxy. The dark matter core size depends on the frequency of globular cluster crossings. Our simulations clearly demonstrate also that between the passages, the dark matter halo can regenerate its cusp. Moreover, our models are in good agreement with constraints on the dark matter masses of globular clusters as our clusters lose a large fraction of their initial dark matter minihalos. These results provide circumstantial evidence for the universal existence of dark matter halos in globular clusters. △ Less

Submitted 4 January, 2020; v1 submitted 16 September, 2019; originally announced September 2019.

Comments: 11 pages, 15 figures, Accepted by MNRAS

Cusp-to-core transition in low-mass dwarf galaxies induced by dynamical heating of cold dark matter by primordial black holes

Authors: Pierre Boldrini, Yohei Miki, Alexander Y. Wagner, Roya Mohayaee, Joseph Silk, Alexandre Arbey

Abstract: We performed a series of high-resolution $N$-body simulations to examine whether dark matter candidates in the form of primordial black holes (PBHs) can solve the cusp-core problem in low-mass dwarf galaxies. If some fraction of the dark matter in low-mass dwarf galaxies consists of PBHs and the rest is cold dark matter, dynamical heating of the cold dark matter by the PBHs induces a cusp-to-core… ▽ More We performed a series of high-resolution $N$-body simulations to examine whether dark matter candidates in the form of primordial black holes (PBHs) can solve the cusp-core problem in low-mass dwarf galaxies. If some fraction of the dark matter in low-mass dwarf galaxies consists of PBHs and the rest is cold dark matter, dynamical heating of the cold dark matter by the PBHs induces a cusp-to-core transition in the total dark matter profile. The mechanism works for PBHs in the 25-100 M$_{\sun}$ mass window, consistent with the LIGO detections, but requires a lower limit on the PBH mass fraction of 1$\%$ of the total dwarf galaxy dark matter content. The cusp-to-core transition time-scale is between 1 and 8 Gyr. This time-scale is also a constant multiple of the relaxation time between cold dark matter particles and PBHs, which depends on the mass, the mass fraction and the scale radius of the initial density profile of PBHs. We conclude that dark matter cores occur naturally in halos comprised of cold dark matter and PBHs, without the need to invoke baryonic processes. △ Less

Submitted 21 January, 2020; v1 submitted 16 September, 2019; originally announced September 2019.

Comments: 8 pages, 7 figures, Accepted by MNRAS

Fornax globular cluster distributions: implications for the cusp-core problem

Authors: P. Boldrini, R. Mohayaee, J. Silk

Abstract: We re-investigate the Fornax cusp-core problem using observational results on the spatial and mass distributions of globular clusters (GCs) in order to put constraints on the dark matter profile. We model Fornax using high-resolution N-body simulations with entirely live systems, i.e. self-gravitating systems composed of stars and dark matter, which account correctly for dynamical friction and tid… ▽ More We re-investigate the Fornax cusp-core problem using observational results on the spatial and mass distributions of globular clusters (GCs) in order to put constraints on the dark matter profile. We model Fornax using high-resolution N-body simulations with entirely live systems, i.e. self-gravitating systems composed of stars and dark matter, which account correctly for dynamical friction and tidal effects between Fornax and the globular clusters. We test two alternative hypotheses, which are a cored and a cuspy halo for Fornax by exploring a reasonable range of initial conditions on globular clusters. For Fornax cored dark matter halo, we derive a lower limit on the core size of $r_{c}\gtrsim$ 0.5 kpc. Contrary to many previous works, we show also that for different initial conditions, a cuspy halo is not ruled out in our simulations based on observations of Fornax globular clusters. △ Less

Submitted 9 March, 2019; v1 submitted 28 February, 2019; originally announced March 2019.

Comments: Accepted by MNRAS. Correction on Figure 1. arXiv admin note: substantial text overlap with arXiv:1806.09591

Evidence for anisotropy of cosmic acceleration

Authors: Jacques Colin, Roya Mohayaee, Mohamed Rameez, Subir Sarkar

Abstract: Observations reveal a `bulk flow' in the local Universe which is faster and extends to much larger scales than is expected around a typical observer in the standard $Λ$CDM cosmology. This is expected to result in a scale-dependent dipolar modulation of the acceleration of the expansion rate inferred from observations of objects within the bulk flow. From a maximum-likelihood analysis of the Joint… ▽ More Observations reveal a `bulk flow' in the local Universe which is faster and extends to much larger scales than is expected around a typical observer in the standard $Λ$CDM cosmology. This is expected to result in a scale-dependent dipolar modulation of the acceleration of the expansion rate inferred from observations of objects within the bulk flow. From a maximum-likelihood analysis of the Joint Lightcurve Analysis (JLA) catalogue of Type Ia supernovae we find that the deceleration parameter, in addition to a small monopole, indeed has a much bigger dipole component aligned with the CMB dipole which falls exponentially with redshift $z$: $q_0 = q_\mathrm{m} + \vec{q}_\mathrm{d}.\hat{n}\exp(-z/S)$. The best fit to data yields $q_\mathrm{d} = -8.03$ and $S = 0.0262~(\Rightarrow d \sim 100~\mathrm{Mpc})$, rejecting isotropy ($q_\mathrm{d} = 0$) with $3.9σ$ statistical significance, while $q_\mathrm{m} = -0.157$ and consistent with no acceleration ($q_\mathrm{m} = 0$) at $1.4σ$. Thus the cosmic acceleration deduced from supernovae may be an artefact of our being non-Copernican observers, rather than evidence for a dominant component of `dark energy' in the Universe. △ Less

Submitted 18 October, 2019; v1 submitted 14 August, 2018; originally announced August 2018.

Comments: 7 pages, 4 figures, 4 tables; Title changed, Bayesian information criterion values added, conclusions strengthened; Results can be reproduced using code at: https://github.com/rameez3333/Dipole_JLA ; Accepted for publication in Astronomy & Astrophysics Letters

Journal ref: A&A 631, L13 (2019)

Does Fornax have a cored halo? Implications for the nature of dark matter

Authors: P. Boldrini, R. Mohayaee, J. Silk

Abstract: Fornax is the most massive of the Milky Way dwarf spheroidal galaxies and has five globular clusters orbiting in a dense background of dark matter. Observational analyses suggest that globular clusters were initially much more massive and lost most of their stars to the Fornax field. We re-investigate the Fornax cusp-core problem, to clarify tensions between simulations and observations concerning… ▽ More Fornax is the most massive of the Milky Way dwarf spheroidal galaxies and has five globular clusters orbiting in a dense background of dark matter. Observational analyses suggest that globular clusters were initially much more massive and lost most of their stars to the Fornax field. We re-investigate the Fornax cusp-core problem, to clarify tensions between simulations and observations concerning the dark matter halo density profile. N-body simulations predict a centrally steep power-law density profile, while observations of the globular clusters seem to prefer that the dark matter halo density is constant at the center. For the first time, we ran pure N-body simulations with entirely live systems, i.e. self-gravitating systems composed of particles (including stars and dark matter). Only this numerical approach accounts correctly for dynamical friction and tidal effects between Fornax and the globular clusters. We show that a weak cusp ($r_{s}=$ 1.5 kpc) or a large core ($r_{c}=$ 848 pc) are not compatible with the current observed positions and masses of Fornax clusters. In contrast, a small dark matter core ($r_{c}=$ 282 pc) for Fornax naturally reproduces the cluster spatial and mass distributions over a wide range of initial globular cluster masses. We derive an upper limit of $r_{c}\lesssim$ 282 pc. This core size range favors only warm dark datter (WDM). It is also possible to obtain a compatible core size range from cold dark matter (CDM) theories, if the initial halo's central cusp is heated by gas to form a small core. △ Less

Submitted 5 March, 2019; v1 submitted 25 June, 2018; originally announced June 2018.

Comments: Substituted by a new paper using much higher resolution simulation with different initial conditions on dark matter halos entitled: Fornax globular cluster distributions: implications for the cusp-core problem and same authors (arXiv:1903.00354)

The dipole anisotropy of AllWISE galaxies

Authors: M. Rameez, R. Mohayaee, S. Sarkar, J. Colin

Abstract: We determine the dipole in the WISE galaxy catalogue. After reducing star contamination to <0.1% by rejecting sources with high apparent motion and those close to the Galactic plane, we eliminate low redshift sources to suppress the non-kinematic, clustering dipole. We remove sources within {\pm}5° of the super-galactic plane, as well as those within 1'' of 2MRS sources at redshift z < 0.03. We en… ▽ More We determine the dipole in the WISE galaxy catalogue. After reducing star contamination to <0.1% by rejecting sources with high apparent motion and those close to the Galactic plane, we eliminate low redshift sources to suppress the non-kinematic, clustering dipole. We remove sources within {\pm}5° of the super-galactic plane, as well as those within 1'' of 2MRS sources at redshift z < 0.03. We enforce cuts on the source angular extent to preferentially select distant ones. As we progress along these steps, the dipole converges in direction to within 5° of the CMB dipole and its magnitude also progressively reduces but stabilises at {\sim}0.012, corresponding to a velocity >1000 km/s if it is solely of kinematic origin. However, previous studies have shown that only {\sim}70% of the velocity of the Local Group as inferred from the CMB dipole is due to sources at z < 0.03. We examine the Dark Sky simulations to quantify the prevalence of such environments and find that <2.1% of Milky Way-like observers in a ΛCDM universe should observe the bulk flow (> 240 km/s extending to z > 0.03) that we do. We construct mock catalogues in the neighbourhood of such peculiar observers in order to mimic our final galaxy selection and quantify the residual clustering dipole. After subtracting this the remaining dipole is 0.0048 {\pm} 0.0022, corresponding to a velocity of 420 {\pm} 213 km/s which is consistent with the CMB. However the sources (at z > 0.03) of such a large clustering dipole remain to be identified. △ Less

Submitted 19 March, 2018; v1 submitted 9 December, 2017; originally announced December 2017.

Comments: Accepted for publication in MNRAS; Minor changes and clarifications, References added; 11 pages, 9 figures, 3 tables

Journal ref: MNRAS 477 (2018) 1772-1781

High redshift radio galaxies and divergence from the CMB dipole

Authors: J. Colin, R. Mohayaee, M. Rameez, S. Sarkar

Abstract: Previous studies have found our velocity in the rest frame of radio galaxies at high redshift to be substantially larger than that inferred from the CMB temperature dipole anisotropy. We construct a full sky catalogue NVSUMSS, by merging the NVSS and SUMSS catalogues and removing local sources by various means including cross-correlating with the 2MRS catalogue. We take into account both aberratio… ▽ More Previous studies have found our velocity in the rest frame of radio galaxies at high redshift to be substantially larger than that inferred from the CMB temperature dipole anisotropy. We construct a full sky catalogue NVSUMSS, by merging the NVSS and SUMSS catalogues and removing local sources by various means including cross-correlating with the 2MRS catalogue. We take into account both aberration and Doppler boost to deduce our velocity from the hemispherical number count asymmetry, as well as via a 3-dimensional linear estimator. Both the magnitude and direction depend on cuts made to the catalogue, e.g. on the lowest source flux, however these effects are small. With the hemispheric number count asymmetry method we obtain a velocity of 1729 $\pm$ 187 km/s i.e. about 4 times larger than that obtained from the CMB dipole, but close in direction, towards RA=149 $\pm$ 2 degree, DEC = -17 $\pm$ 12 degree. With the 3-dimensional estimator, the derived velocity is 1355 $\pm$ 174 km/s towards RA=141 $\pm$ 11 degree, DEC=-9 $\pm$ 10 degree. We assess the statistical significance of these results by constructing catalogues of random distributions and show that they are at best significant at the $2.81 σ$ (99.95% confidence) level. △ Less

Submitted 9 August, 2017; v1 submitted 27 March, 2017; originally announced March 2017.

Comments: 13 pages, 12 figures, 10 tables; matches published version

Journal ref: Monthly Notices of the Royal Astronomical Society 471 (2017) 1045

Critical points of the cosmic velocity field and the uncertainties in the value of the Hubble constant

Authors: Hao Liu, Roya Mohayaee, Pavel Naselsky

Abstract: The existence of critical points for the peculiar velocity field is a natural feature of the correlated vector field. These points appear at the junctions of velocity domains with different orientations of their averaged velocity vectors. Since peculiar velocities are the important cause of the scatter in the Hubble expansion rate, we propose that a more precise determination of the Hubble constan… ▽ More The existence of critical points for the peculiar velocity field is a natural feature of the correlated vector field. These points appear at the junctions of velocity domains with different orientations of their averaged velocity vectors. Since peculiar velocities are the important cause of the scatter in the Hubble expansion rate, we propose that a more precise determination of the Hubble constant can be made by restricting analysis to a subsample of observational data containing only the zones around the critical points of the peculiar velocity field, associated with voids and saddle points. On large-scales the critical points, where the first derivative of the gravitational potential vanishes, can easily be identified using the density field and classified by the behavior of the Hessian of the gravitational potential. We use high-resolution N-body simulations to show that these regions are stable in time and hence are excellent tracers of the initial conditions. Furthermore, we show that the variance of the Hubble flow can be substantially minimized by restricting observations to the subsample of such regions of vanishing velocity instead of aiming at increasing the statistics by averaging indiscriminately using the full data sets, as is the common approach. △ Less

Submitted 20 May, 2016; v1 submitted 17 May, 2016; originally announced May 2016.

Comments: 17 pages, 6 figures and 2 tables, minor modifications after receiving several useful comments

Journal ref: JCAP06(2016)009

Unraveling the nature of Gravity through our clumpy Universe

Authors: Shant Baghram, Saeed Tavasoli, Farhang Habibi, Roya Mohayaee, Joseph Silk

Abstract: We propose a new probe to test the nature of gravity at various redshifts through large-scale cosmological observations. We use our void catalog, extracted from the Sloan Digital Sky Survey (SDSS, DR10), to trace the distribution of matter along the lines of sight to SNe Ia that are selected from the Union 2 catalog. We study the relation between SNe Ia luminosities and convergence and also the pe… ▽ More We propose a new probe to test the nature of gravity at various redshifts through large-scale cosmological observations. We use our void catalog, extracted from the Sloan Digital Sky Survey (SDSS, DR10), to trace the distribution of matter along the lines of sight to SNe Ia that are selected from the Union 2 catalog. We study the relation between SNe Ia luminosities and convergence and also the peculiar velocities of the sources. We show that the effects, on SNe Ia luminosities, of convergence and of peculiar velocities predicted by the theory of general relativity and theories of modified gravities are different and hence provide a new probe of gravity at various redshifts. We show that the present sparse large-scale data does not allow us to determine any statistically- significant deviation from the theory of general relativity but future more comprehensive surveys should provide us with means for such an exploration. △ Less

Submitted 25 November, 2014; originally announced November 2014.

Comments: Essay received honorable mention in the Gravity Research Foundation 2014 essay competition. To appear in IJMPD. 10 pages, 2 figures

Journal ref: International Journal of Modern Physics D Vol. 23, No. 12, 1442025 (2014)

A single-merger scenario for the formation of the giant stream and the warp of M31

Authors: Raphael Sadoun, Roya Mohayaee, Jacques Colin

Abstract: We propose that the accretion of a dwarf spheroidal galaxy provides a common origin for the giant southern stream and the warp of M31. We run about 40 full N-body simulations with live M31, infalling galaxies with varying masses and density profiles, and cosmologically-plausible initial orbital parameters. Excellent agreement with a full range of observational data is obtained for a model in which… ▽ More We propose that the accretion of a dwarf spheroidal galaxy provides a common origin for the giant southern stream and the warp of M31. We run about 40 full N-body simulations with live M31, infalling galaxies with varying masses and density profiles, and cosmologically-plausible initial orbital parameters. Excellent agreement with a full range of observational data is obtained for a model in which a dark-matter-rich dwarf spheroidal, whose trajectory lies on the thin plane of corotating satellites of M31, is accreted from its turnaround radius of about 200 kpcs into M31 at approximately 3 Gyrs ago. The satellite is disrupted as it orbits in the potential well of the galaxy and forms the giant stream and in return heats and warps the disk of M31. We show that our cosmologically-motivated model is favoured by the kinematic data over the phenomenological models in which the satellite starts its infall from a close distance of M31. Our model predicts that the remnant of the disrupted satellite resides in the region of the North-Eastern shelf of M31. The results here suggest that the surviving satellites of M31 that orbit on the same thin plane, as the disrupted satellite once did, could have all been accreted from an intergalactic filament. △ Less

Submitted 18 July, 2013; originally announced July 2013.

Comments: 18 pages, 22 figures, 3 tables

The challenge of large and empty voids in the SDSS DR7 redshift survey

Authors: Saeed Tavasoli, Kaveh Vasei, Roya Mohayaee

Abstract: We present catalogues of voids for the SDSS DR7 redshift survey and for Millennium I simulation mock data. We aim to compare the observations with simulations based on a $Λ$CDM model and a semi-analytic galaxy formation model. We use the void statistics as a test for these models. We assembled a mock catalogue that closely resembles the SDSS DR7 catalogue and carried out a parallel statistical ana… ▽ More We present catalogues of voids for the SDSS DR7 redshift survey and for Millennium I simulation mock data. We aim to compare the observations with simulations based on a $Λ$CDM model and a semi-analytic galaxy formation model. We use the void statistics as a test for these models. We assembled a mock catalogue that closely resembles the SDSS DR7 catalogue and carried out a parallel statistical analysis of the observed and simulated catalogue. We find that in the observation and the simulation, voids tend to be equally spherical. The total volume occupied by the voids and their total number are slightly larger in the simulation than in the observation. We find that large voids are less abundant in the simulation and the total luminosity of the galaxies contained in a void with a given radius is higher on average than observed by SDSS DR7 survey. We expect these discrepancies to be even more significant in reality than found here since the present value of $σ_8$ given by WMAP7 is lower than the value of 0.9 used in the Millennium I simulation. The reason why the simulation fails to produce enough large and dark voids might be the failure of certain semi-analytic galaxy formation models to reduce the small-scale power of $Λ$CDM and to produce sufficient power on large scales. △ Less

Submitted 18 February, 2013; v1 submitted 8 October, 2012; originally announced October 2012.

Comments: Kolmogrov-Smirnov test added. version to appear in A&A

Refined Second Law of Thermodynamics for fast random processes

Authors: Erik Aurell, Krzysztof Gawȩdzki, Carlos Mejía-Monasterio, Roya Mohayaee, Paolo Muratore-Ginanneschi

Abstract: We establish a refined version of the Second Law of Thermodynamics for Langevin stochastic processes describing mesoscopic systems driven by conservative or non-conservative forces and interacting with thermal noise. The refinement is based on the Monge-Kantorovich optimal mass transport. General discussion is illustrated by numerical analysis of a model for micron-size particle manipulated by opt… ▽ More We establish a refined version of the Second Law of Thermodynamics for Langevin stochastic processes describing mesoscopic systems driven by conservative or non-conservative forces and interacting with thermal noise. The refinement is based on the Monge-Kantorovich optimal mass transport. General discussion is illustrated by numerical analysis of a model for micron-size particle manipulated by optical tweezers. △ Less

Submitted 16 January, 2012; originally announced January 2012.

Comments: 17 pages

Journal ref: Journal of Statistical Physics 147, 487--505 (2012)

Enhancements to velocity-dependent dark matter interactions from tidal streams and shells in the Andromeda galaxy

Authors: Robyn E. Sanderson, Roya Mohayaee, Joe Silk

Abstract: Dark matter substructure around nearby galaxies provides an interesting opportunity for confusion-free indirect detection of dark matter. We calculate the boost over a smooth background distribution of dark matter for gamma-ray emission from dark matter self-annihilations in tidal structure in M31, assuming a cross-section inversely proportional to the relative velocities of the dark matter partic… ▽ More Dark matter substructure around nearby galaxies provides an interesting opportunity for confusion-free indirect detection of dark matter. We calculate the boost over a smooth background distribution of dark matter for gamma-ray emission from dark matter self-annihilations in tidal structure in M31, assuming a cross-section inversely proportional to the relative velocities of the dark matter particles as proposed by the Sommerfeld effect. The low velocity of the material in the structure results in a significant increase in gamma-ray emission compared to both the background halo and the predicted emission for a velocity-independent cross section. We also calculate the expected signal for Fermi, for reasonable choices of the dark matter parameters. We find that for a cross section proportional to the inverse-square of the relative velocity, the enhancement to the annihilation rate is sufficient to test the velocity dependence of the cross section by spatial correlation with the stellar component of the stream, given sufficient detector sensitivity. △ Less

Submitted 8 November, 2011; v1 submitted 21 June, 2011; originally announced June 2011.

Comments: new version revised per referee's comments

The low frequency of dual AGNs versus the high merger rate of galaxies: A phenomenological model

Authors: Qingjuan Yu, Youjun Lu, Roya Mohayaee, Jacques Colin

Abstract: Dual AGNs are natural byproducts of hierarchical mergers of galaxies in the LambdaCDM cosmogony. Recent observations have shown that only a small fraction (~ 0.1%-1%) of AGNs at redshift z<~ 0.3 are dual with kpc-scale separations, which is rather low compared to the high merger rate of galaxies. Here we construct a phenomenological model to estimate the number density of dual AGNs and its evoluti… ▽ More Dual AGNs are natural byproducts of hierarchical mergers of galaxies in the LambdaCDM cosmogony. Recent observations have shown that only a small fraction (~ 0.1%-1%) of AGNs at redshift z<~ 0.3 are dual with kpc-scale separations, which is rather low compared to the high merger rate of galaxies. Here we construct a phenomenological model to estimate the number density of dual AGNs and its evolution according to the observationally-estimated major merger rates of galaxies and various scaling relations on the properties of galaxies and their central massive black holes. We show that our model reproduces the observed frequency and separation distribution of dual AGNs provided that significant nuclear activities are triggered only in gas-rich progenitor galaxies with central massive black holes and only when the nuclei of these galaxies are roughly within the half-light radii of their companion galaxies. Under these constraints, the observed low dual AGN frequency is consistent with the relatively high merger rate of galaxies and supports the hypothesis that major mergers lead to AGN/QSO activities. We also predict that the number of kpc-scale dual AGNs decreases with increasing redshift and only about 0.02%--0.06% of AGNs are dual AGNs with double-peaked narrow line features at redshifts of z 0.5-1.2. Future observations of high-redshift dual AGNs would provide a solid test for this prediction. △ Less

Submitted 10 May, 2011; originally announced May 2011.

Comments: 8 pages, 3 figures

Probing the anisotropic local universe and beyond with SNe Ia data

Authors: Jacques Colin, Roya Mohayaee, Subir Sarkar, Arman Shafieloo

Abstract: The question of the transition to global isotropy from our anisotropic local Universe is studied using the Union 2 catalogue of Type Ia supernovae (SNe Ia). We construct a "residual" statistic sensitive to systematic shifts in their brightness in different directions and use this to search in different redshift bins for a preferred direction on the sky in which the SNe Ia are brighter or fainter r… ▽ More The question of the transition to global isotropy from our anisotropic local Universe is studied using the Union 2 catalogue of Type Ia supernovae (SNe Ia). We construct a "residual" statistic sensitive to systematic shifts in their brightness in different directions and use this to search in different redshift bins for a preferred direction on the sky in which the SNe Ia are brighter or fainter relative to the 'standard' LCDM cosmology. At low redshift (z<0.05) we find that an isotropic model such as LCDM is barely consistent with the SNe Ia data at 2-3 sigma. A complementary maximum likelihood analysis of peculiar velocities confirms this finding -- there is a bulk flow of around 260 km/sec at z \sim 0.06, which disagrees with LCDM at 1-2 sigma. Since the Shapley concentration is believed to be largely responsible for this bulk flow, we make a detailed study of the infall region: the SNe Ia falling away from the Local Group towards Shapley are indeed significantly dimmer than those falling towards us and on to Shapley. Convergence to the CMB rest frame must occur well beyond Shapley (z>0.06) so the low redshift bulk flow can systematically bias any reconstruction of the expansion history of the Universe. At high redshifts z>0.15 the agreement between the SNe Ia data and the isotropic LCDM model does improve, however, the sparseness and low quality of the data means that LCDM cannot be singled out as the preferred cosmological model. △ Less

Submitted 11 April, 2011; v1 submitted 29 November, 2010; originally announced November 2010.

Comments: 9 pages, 9 figures; A new sentence is added to the caption of Figure 8 for clarification

Journal ref: Mon. Not. Roy. Astron. Soc. 414 (2011) 264-271

Nonspherical similarity solutions for dark halo formation

Authors: Mark Vogelsberger, Roya Mohayaee, Simon D. M. White

Abstract: We carry out fully 3-dimensional simulations of evolution from self-similar, spherically symmetric linear perturbations of a Cold Dark Matter dominated Einstein-de Sitter universe. As a result of the radial orbit instability, the haloes which grow from such initial conditions are triaxial with major-to-minor axis ratios of order 3:1. They nevertheless grow approximately self-similarly in time. In… ▽ More We carry out fully 3-dimensional simulations of evolution from self-similar, spherically symmetric linear perturbations of a Cold Dark Matter dominated Einstein-de Sitter universe. As a result of the radial orbit instability, the haloes which grow from such initial conditions are triaxial with major-to-minor axis ratios of order 3:1. They nevertheless grow approximately self-similarly in time. In all cases they have power-law density profiles and near-constant velocity anisotropy in their inner regions. Both the power-law index and the value of the velocity anisotropy depend on the similarity index of the initial conditions, the former as expected from simple scaling arguments. Halo structure is thus not "universal" but remembers the initial conditions. On larger scales the density and anisotropy profiles show two characteristic scales, corresponding to particles at first pericentre and at first apocentre after infall. They are well approximated by the NFW model only for one value of the similarity index. In contrast, at all radii within the outer caustic the pseudo phase-space density can be fit by a single power law with an index which depends only very weakly on the similarity index of the initial conditions. This behaviour is very similar to that found for haloes formed from LCDM initial conditions and so can be considered approximately universal. △ Less

Submitted 23 July, 2010; originally announced July 2010.

Comments: 8 pages, 7 figures, submitted to MNRAS

Hierarchy in the Phase Space and Dark Matter Astronomy

Authors: Niayesh Afshordi, Roya Mohayaee, Edmund Bertschinger

Abstract: We develop a theoretical framework for describing the hierarchical structure of the phase space of cold dark matter haloes, due to gravitationally bound substructures. Because it includes the full hierarchy of the cold dark matter initial conditions and is hence complementary to the halo model, the stable clustering hypothesis is applied for the first time here to the small-scale phase space str… ▽ More We develop a theoretical framework for describing the hierarchical structure of the phase space of cold dark matter haloes, due to gravitationally bound substructures. Because it includes the full hierarchy of the cold dark matter initial conditions and is hence complementary to the halo model, the stable clustering hypothesis is applied for the first time here to the small-scale phase space structure. As an application, we show that the particle dark matter annihilation signal could be up to two orders of magnitude larger than that of the smooth halo within the Galactic virial radius. The local boost is inversely proportional to the smooth halo density, and thus is O(1) within the solar radius, which could translate into interesting signatures for dark matter direct detection experiments: The temporal correlation of dark matter detection can change by a factor of 2 in the span of 10 years, while there will be significant correlations in the velocity space of dark matter particles. This can introduce O(1) uncertainty in the direction of local dark matter wind, which was believed to be a benchmark of directional dark matter searches or the annual modulation signal. △ Less

Submitted 2 November, 2009; originally announced November 2009.

Comments: 5 pages, 4 figures

Journal ref: Phys.Rev.D81:101301,2010

Caustics in growing Cold Dark Matter Haloes

Authors: Mark Vogelsberger, Simon D. M. White, Roya Mohayaee, Volker Springel

Abstract: We simulate the growth of isolated dark matter haloes from self-similar and spherically symmetric initial conditions. Our N-body code integrates the geodesic deviation equation in order to track the streams and caustics associated with individual simulation particles. The radial orbit instability causes our haloes to develop major-to-minor axis ratios approaching 10 to 1 in their inner regions.… ▽ More We simulate the growth of isolated dark matter haloes from self-similar and spherically symmetric initial conditions. Our N-body code integrates the geodesic deviation equation in order to track the streams and caustics associated with individual simulation particles. The radial orbit instability causes our haloes to develop major-to-minor axis ratios approaching 10 to 1 in their inner regions. They grow similarly in time and have similar density profiles to the spherical similarity solution, but their detailed structure is very different. The higher dimensionality of the orbits causes their stream and caustic densities to drop much more rapidly than in the similarity solution. This results in a corresponding increase in the number of streams at each point. At 1% of the turnaround radius (corresponding roughly to the Sun's position in the Milky Way) we find of order 10^6 streams in our simulations, as compared to 10^2 in the similarity solution. The number of caustics in the inner halo increases by a factor of several, because a typical orbit has six turning points rather than one, but caustic densities drop by a much larger factor. This reduces the caustic contribution to the annihilation radiation. For the region between 1% and 50% of the turnaround radius, this is 4% of the total in our simulated haloes, as compared to 6.5% in the similarity solution. Caustics contribute much less at smaller radii. These numbers assume a 100 GeV c^-2 neutralino with present-day velocity dispersion 0.03 cm s^-1, but reducing the dispersion by ten orders of magnitude only doubles the caustic luminosity. We conclude that caustics will be unobservable in the inner parts of haloes. Only the outermost caustic might potentially be detectable. △ Less

Submitted 4 September, 2009; v1 submitted 24 June, 2009; originally announced June 2009.

Comments: 12 pages, 17 figures, minor corrections, accepted by MNRAS

Hierarchical Phase Space Structure of Dark Matter Haloes: Tidal debris, Caustics, and Dark Matter annihilation

Authors: Niayesh Afshordi, Roya Mohayaee, Edmund Bertschinger

Abstract: Most of the mass content of dark matter haloes is expected to be in the form of tidal debris. The density of debris is not constant, but rather can grow due to formation of caustics at the apocenters and pericenters of the orbit, or decay as a result of phase mixing. In the phase space, the debris assemble in a hierarchy which is truncated by the primordial temperature of dark matter. Understand… ▽ More Most of the mass content of dark matter haloes is expected to be in the form of tidal debris. The density of debris is not constant, but rather can grow due to formation of caustics at the apocenters and pericenters of the orbit, or decay as a result of phase mixing. In the phase space, the debris assemble in a hierarchy which is truncated by the primordial temperature of dark matter. Understanding this phase structure can be of significant importance for the interpretation of many astrophysical observations and in particular dark matter detection experiments. With this purpose in mind, we develop a general theoretical framework to describe the hierarchical structure of the phase space of cold dark matter haloes. We do not make any assumption of spherical symmetry and/or smooth and continuous accretion. Instead, working with correlation functions in the action-angle space, we can fully account for the hierarchical structure (predicting a two-point correlation function ~ (ΔJ)^{-1.6} in the action space), as well as the primordial discreteness of the phase space. As an application, we estimate the boost to the dark matter annihilation signal due to the structure of the phase space within virial radius: the boost due to the hierarchical tidal debris is of order unity, whereas the primordial discreteness of the phase structure can boost the total annihilation signal by up to an order of magnitude. The latter is dominated by the regions beyond 20% of the virial radius, and is largest for the recently formed haloes with the least degree of phase mixing. △ Less

Submitted 10 November, 2008; originally announced November 2008.

Comments: 12 pages, 7 figures, comments are welcome

Journal ref: Phys.Rev.D79:083526,2009

Cosmic flow from 2MASS redshift survey: The origin of CMB dipole and implications for LCDM cosmology

Authors: G. Lavaux, R. Brent Tully, R. Mohayaee, S. Colombi

Abstract: We generate the peculiar velocity field for the 2MASS Redshift Survey (2MRS) catalog using an orbit-reconstruction algorithm. The reconstructed velocities of individual objects in 2MRS are well-correlated with the peculiar velocities obtained from high-precision observed distances within 3,000 km/s. We estimate the mean matter density to be 0.31 +/- 0.05 by comparing observed to reconstructed velo… ▽ More We generate the peculiar velocity field for the 2MASS Redshift Survey (2MRS) catalog using an orbit-reconstruction algorithm. The reconstructed velocities of individual objects in 2MRS are well-correlated with the peculiar velocities obtained from high-precision observed distances within 3,000 km/s. We estimate the mean matter density to be 0.31 +/- 0.05 by comparing observed to reconstructed velocities in this volume. The reconstructed motion of the Local Group in the rest frame established by distances within 3,000 km/s agrees with the observed motion and is generated by fluctuations within this volume, in agreement with observations. Then, we reconstruct the velocity field of 2MRS in successively larger radii, to study the problem of convergence towards the CMB dipole. We find that less than half of the amplitude of the CMB dipole is generated within a volume enclosing the Hydra-Centaurus-Norma supercluster at around 40 Mpc/h. Although most of the amplitude of the CMB dipole seems to be recovered by 120 Mpc/h, the direction does not agree and hence we observe no convergence up to this scale. We develop a statistical model which allows us to estimate cosmological para meters from the reconstructed growth of convergence of the velocity of the Local Group towards the CMB dipole motion. For scales up to 60 Mpc/h, assuming a Local Group velocity of 627 km/s, we estimate Omega_m h^2 = 0.11 +/- 0.06 and sigma_8=0.9 +/- 0.4, in agreement with WMAP5 measurements at the 1-sigma level. However, for scales up to 100 Mpc/h, we obtain Omega_m h^2 = 0.08 +/- 0.03 and sigma_8=1.0 +/- 0.4, which agrees at the 1 to 2-sigma level with WMAP5 results. (abridged) △ Less

Submitted 7 February, 2012; v1 submitted 20 October, 2008; originally announced October 2008.

Comments: 19 pages, 9 figures, 4 tables, accepted by ApJ. No modifications since last version. The modeled velocity field is available at first author's webpage: http://guilhem.lavaux.free.fr/wordpress/?page_id=281

Journal ref: Astrophys.J.709:483-498,2010

Dark matter accretion wakes of high-redshift black holes

Authors: Roya Mohayaee, Jacques Colin

Abstract: Anisotropic emission of gravitational waves during the merger or formation of black holes can lead to the ejection of these black holes from their host galaxies. A recoiled black hole which moves on an almost radial bound orbit outside the virial radius of its central galaxy, in the cold dark matter background, reaches its apapsis in a finite time. The low value of dark matter velocity dispersio… ▽ More Anisotropic emission of gravitational waves during the merger or formation of black holes can lead to the ejection of these black holes from their host galaxies. A recoiled black hole which moves on an almost radial bound orbit outside the virial radius of its central galaxy, in the cold dark matter background, reaches its apapsis in a finite time. The low value of dark matter velocity dispersion at high redshifts and also the black hole velocity near the apapsis passage yield a high-density wake around these black holes. Gamma-ray emission can result from the enhancement of dark matter annihilation in these wakes. The diffuse high-energy gamma-ray background from the ensemble of such black holes in the Hubble volume is also evaluated. △ Less

Submitted 18 February, 2008; originally announced February 2008.

Comments: Talk presented at "Jean-Pierre Lasota, X-ray binaries, accretion disks and compact stars" (October 2007); Abramowicz, M. Ed., New Astronomy Review, in press

Journal ref: New Astron.Rev.51:898-905,2008

The Cosmic Ray Signature of Dark Matter Caustics

Authors: Roya Mohayaee, Pierre Salati

Abstract: Gravitational collapse of dark matter, merger of dark matter haloes and tidal disruption of satellites are among processes which lead to the formation of fine and dense dark matter shells, also known as dark matter caustics. The putative weakly interacting species which may form the dark matter are expected to strongly annihilate in these dense regions of the Milky Way halo and generate in parti… ▽ More Gravitational collapse of dark matter, merger of dark matter haloes and tidal disruption of satellites are among processes which lead to the formation of fine and dense dark matter shells, also known as dark matter caustics. The putative weakly interacting species which may form the dark matter are expected to strongly annihilate in these dense regions of the Milky Way halo and generate in particular antiprotons and positrons. We derive the flux of these rare antimatter particles at the Earth and show that it depends significantly on the cut-off radius of the dark matter distribution at the galactic centre. Boost factors of ~30 are found with respect to a smooth NFW profile for high-energy antiprotons and low-energy positrons if this cut-off radius is taken to be 300 pc -- a somewhat extreme value though. This yields a detectable antiproton signal around hundreds of Gev in models where the annihilation cross section today is enhanced by non--perturbative effects as in the generic case of a heavy Wino. However, dark matter caustics cannot provide a better explanation for the HEAT excess reported above ~10 GeV than a smooth NFW or isothermal cored distribution. △ Less

Submitted 23 June, 2008; v1 submitted 21 January, 2008; originally announced January 2008.

Comments: 16 pages, 11 figures, version to appear in MNRAS

Report number: LAPTH-1235/08

The Monge-Ampère-Kantorovich approach to reconstruction in cosmology

Authors: Roya Mohayaee, Andrei Sobolevskii

Abstract: Motion of a continuous fluid can be decomposed into an "incompressible" rearrangement, which preserves the volume of each infinitesimal fluid element, and a gradient map that transfers fluid elements in a way unaffected by any pressure or elasticity (the polar decomposition of Y. Brenier). The Euler equation describes a system whose kinematics is dominated by the incompressible rearrangement. Th… ▽ More Motion of a continuous fluid can be decomposed into an "incompressible" rearrangement, which preserves the volume of each infinitesimal fluid element, and a gradient map that transfers fluid elements in a way unaffected by any pressure or elasticity (the polar decomposition of Y. Brenier). The Euler equation describes a system whose kinematics is dominated by the incompressible rearrangement. The opposite limit, in which the incompressible component is negligible, corresponds to the Zel'dovich approximation, a model of motion of self-gravitating fluid in cosmology. We present a method of approximate reconstruction of the large-scale proper motions of matter in the Universe from the present-day mass density field. The method is based on recovering the corresponding gradient transfer map. We discuss its algorithmics, tests of the method against mock cosmological catalogues, and its application to observational data, which result in tight constraints on the mean mass density Omega_m and age of the Universe. △ Less

Submitted 16 December, 2007; originally announced December 2007.

Comments: 6 pages, 2 figures; based on an invited lecture at the conference "Euler's Equations: 250 Years On" (see http://www.obs-nice.fr/etc7/EE250/); to be published in a special issue of Physica D containing the proceedings of that conference

Journal ref: PhysicaD237:2145-2150,2008

Gamma-ray emission from dark matter wakes of recoiled black holes

Authors: Roya Mohayaee, Jacques Colin, Joe Silk

Abstract: A new scenario for the emission of high-energy gamma-rays from dark matter annihilation around massive black holes is presented. A black hole can leave its parent halo, by means of gravitational radiation recoil, in a merger event or in the asymmetric collapse of its progenitor star. A recoiled black hole which moves on an almost-radial orbit outside the virial radius of its central halo, in the… ▽ More A new scenario for the emission of high-energy gamma-rays from dark matter annihilation around massive black holes is presented. A black hole can leave its parent halo, by means of gravitational radiation recoil, in a merger event or in the asymmetric collapse of its progenitor star. A recoiled black hole which moves on an almost-radial orbit outside the virial radius of its central halo, in the cold dark matter background, reaches its apapsis in a finite time. Near or at the apapsis passage, a high-density wake extending over a large radius of influence, forms around the black hole. It is shown that significant gamma-ray emission can result from the enhancement of neutralino annihilation in these wakes. At its apapsis passage, a black hole is shown to produce a flash of high-energy gamma-rays whose duration is determined by the mass of the black hole and the redshift at which it is ejected. The ensemble of such black holes in the Hubble volume is shown to produce a diffuse high-energy gamma-ray background whose magnitude is compared to the diffuse emission from dark matter haloes alone. △ Less

Submitted 5 February, 2008; v1 submitted 21 September, 2007; originally announced September 2007.

Comments: version to appear in Astrophysical Journal letters (labels on Fig. 3 corrected)

Using Local Volume data to constrain Dark Matter dynamics

Authors: G. Lavaux, R. Mohayaee, S. Colombi, R. B. Tully

Abstract: The peculiar velocity reconstruction methods allow one to have a deeper insight into the distribution of dark matter: both to measure mean matter density and to obtain the primordial density fluctuations. We present here the Monge-Ampere-Kantorovitch method applied to mock catalogues mimicking in both redshift and distance catalogues. After having discussed the results obtained for a class of bi… ▽ More The peculiar velocity reconstruction methods allow one to have a deeper insight into the distribution of dark matter: both to measure mean matter density and to obtain the primordial density fluctuations. We present here the Monge-Ampere-Kantorovitch method applied to mock catalogues mimicking in both redshift and distance catalogues. After having discussed the results obtained for a class of biases that may be corrected for, we focus on the systematics coming from the unknown distribution of unobserved mass and from the statistical relationship between mass and luminosity. We then show how to use these systematics to put constraints on the dark matter distribution. Finally a preliminary application to an extended version (c z < 3000 km/s) of the Neighbour Galaxy Catalogue is presented. We recover the peculiar velocities in our neighbourhood and present a preliminary measurement of the local Omega_M. △ Less

Submitted 20 August, 2007; originally announced August 2007.

Comments: 4 pages, 2 figures. To be published in the proceedings of ``Galaxies in the Local Volume'', Sydney 8 to 13 July 2007

Observational biases in Lagrangian reconstructions of cosmic velocity fields

Authors: G. Lavaux, R. Mohayaee, S. Colombi, R. B. Tully, F. Bernardeau, J. Silk

Abstract: Lagrangian reconstruction of large-scale peculiar velocity fields can be strongly affected by observational biases. We develop a thorough analysis of these systematic effects by relying on specially selected mock catalogues. For the purpose of this paper, we use the MAK reconstruction method, although any other Lagrangian reconstruction method should be sensitive to the same problems. We extensi… ▽ More Lagrangian reconstruction of large-scale peculiar velocity fields can be strongly affected by observational biases. We develop a thorough analysis of these systematic effects by relying on specially selected mock catalogues. For the purpose of this paper, we use the MAK reconstruction method, although any other Lagrangian reconstruction method should be sensitive to the same problems. We extensively study the uncertainty in the mass-to-light assignment due to luminosity incompleteness, and the poorly-determined relation between mass and luminosity. The impact of redshift distortion corrections is analyzed in the context of MAK and we check the importance of edge and finite-volume effects on the reconstructed velocities. Using three mock catalogues with different average densities, we also study the effect of cosmic variance. In particular, one of them presents the same global features as found in observational catalogues that extend to 80 Mpc/h scales. We give recipes, checked using the aforementioned mock catalogues, to handle these particular observational effects, after having introduced them into the mock catalogues so as to quantitatively mimic the most densely sampled currently available galaxy catalogue of the nearby universe. Once biases have been taken care of, the typical resulting error in reconstructed velocities is typically about a quarter of the overall velocity dispersion, and without significant bias. We finally model our reconstruction errors to propose an improved Bayesian approach to measure Omega_m in an unbiased way by comparing the reconstructed velocities to the measured ones in distance space, even though they may be plagued by large errors. We show that, in the context of observational data, a nearly unbiased estimator of Omega_m may be built using MAK reconstruction. △ Less

Submitted 12 October, 2007; v1 submitted 24 July, 2007; originally announced July 2007.

Comments: 29 pages, 21 figures, 6 tables, Accepted by MNRAS on 2007 October 2. Received 2007 September 30; in original form 2007 July 24

Dark matter caustics and the enhancement of self-annihilation flux

Authors: Roya Mohayaee, Sergei Shandarin, Joseph Silk

Abstract: Cold dark matter haloes are populated by caustics, which are yet to be resolved in N-body simulations or observed in the Universe. Secondary infall model provides a paradigm for the study of caustics in "typical" haloes assuming that they have had no major mergers and have grown only by smooth accretion. This is a particular characteristic of the smallest dark matter haloes of about 10^{-5} Mo,… ▽ More Cold dark matter haloes are populated by caustics, which are yet to be resolved in N-body simulations or observed in the Universe. Secondary infall model provides a paradigm for the study of caustics in "typical" haloes assuming that they have had no major mergers and have grown only by smooth accretion. This is a particular characteristic of the smallest dark matter haloes of about 10^{-5} Mo, which although "atypical" contain no substructures and could have survived until now with no major mergers. Thus using this model as the first guidline, we evaluate the neutralino self-annihilation flux for these haloes. Our results show that caustics could leave a distinct sawteeth signature on the differential and cumulative fluxes coming from the outer regions of these haloes. The total annihilation signal from the regions away from the centre can be boosted by about forty percents. △ Less

Submitted 16 April, 2007; originally announced April 2007.

Comments: To appear in JCAP, 4 pages, 3 figures

Journal ref: JCAP 0705:015,2007

Caustics in Dark Matter Haloes

Authors: Roya Mohayaee, Stephane Colombi, Bernard Fort, Raphael Gavazzi, Sergei Shandarin, Jihad Touma

Abstract: Caustics are formally singular structures, with infinite density, that form in collisionless media. The non-negligible velocity dispersion of dark matter particles renders their density finite. We evaluate the maximum density of the caustics within the framework of secondary infall model of formation of dark matter haloes. The result is then used to demonstrate that caustics can be probed by pro… ▽ More Caustics are formally singular structures, with infinite density, that form in collisionless media. The non-negligible velocity dispersion of dark matter particles renders their density finite. We evaluate the maximum density of the caustics within the framework of secondary infall model of formation of dark matter haloes. The result is then used to demonstrate that caustics can be probed by properly stacking the weak-lensing signal of about 600 haloes. CFHTLS accompanied by X-ray observations and the space-based experiments like SNAP or DUNE can provide us with the required statistics and hence a way of distinguishing between the viable dark matter particle candidates. The extension of our results to more realistic models including the effects of mergers of haloes is briefly outlined. △ Less

Submitted 24 October, 2005; v1 submitted 19 October, 2005; originally announced October 2005.

Comments: Minor changes, two references added, 6 pages, 3 figures, to appear in Proc. 21st IAP Colloquium "Mass Profiles and Shapes of Cosmological Structures", Paris 4-9 July 2005, [EAS Publications Series, eds: G. Mamon, F. Combes, C. Deffayet, B. Fort]

The Cosmological Mean Density and its Local Variations Probed by Peculiar Velocities

Authors: Roya Mohayaee, R. Brent Tully

Abstract: Peculiar velocities thoughout the region of the local supercluster are reconstructed by two different orbit-retracing methods. The requirement of the optimal correlation between the radial components of reconstructed velocities and the observed peculiar velocities derived from our extensive new catalog of distances puts stringent constraints on the values of the cosmological parameters. Our cons… ▽ More Peculiar velocities thoughout the region of the local supercluster are reconstructed by two different orbit-retracing methods. The requirement of the optimal correlation between the radial components of reconstructed velocities and the observed peculiar velocities derived from our extensive new catalog of distances puts stringent constraints on the values of the cosmological parameters. Our constraints intersect those from studies of microwave background fluctuations and statistical properties of galaxy clustering: the ensemble of constraints are consistent with Omega_m=0.22\pm 0.02. While motions throughout the Local Supercluster provide a measure of the mean ratio of mass to light, there can be large local fluctuations. Our reconstruction of the infall velocities in the immediate vicinity of the Virgo Cluster shows that there is a mass-to-light anomaly of a factor of 3 to 6 between groups in the general field environment and the heavily populated Virgo Cluster. △ Less

Submitted 7 November, 2005; v1 submitted 13 September, 2005; originally announced September 2005.

Comments: 4 pages, 2 figures, version to appear in Astrophysical Journal Letters

Journal ref: Astrophys.J.635:L113-L116,2005

Probing dark matter caustics with weak lensing

Authors: Raphael Gavazzi, Roya Mohayaee, Bernard Fort

Abstract: Caustics are high-density structures that form frequently in collisionless media. Under self-gravity, cold dark matter flows focus onto caustics which are yet to be resolved in numerical simulations and or observed in the real world. If detected, caustics would provide a strong evidence for dark matter and would rule out alternative models such as those with modified dynamics. Here, we demonstra… ▽ More Caustics are high-density structures that form frequently in collisionless media. Under self-gravity, cold dark matter flows focus onto caustics which are yet to be resolved in numerical simulations and or observed in the real world. If detected, caustics would provide a strong evidence for dark matter and would rule out alternative models such as those with modified dynamics. Here, we demonstrate how they might be observed in the weak lensing data. We evaluate the shear distortion and show that its radial profile is marked by a characteristic sawtooth pattern due to the caustics in dark matter haloes that form by selfsimilar accretion. We discuss the observational complications, mainly due to the poor knowledge of the virial radii of the haloes and demonstrate that a superposition of about 300 cluster-size haloes would give a signal-to-noise ratio which is sufficiently large for the detection of caustics with ground-based observations. This number is reduced to 100 for space-based observations. These bounds can be easily achieved by the ongoing wide field optical surveys such as CFHTLS and the future space-based projects SNAP and DUNE which have to be accompanied by a X-ray follow-up of the selected clusters for a precise determination of their virial radii. △ Less

Submitted 16 March, 2006; v1 submitted 2 June, 2005; originally announced June 2005.

Journal ref: Astron.Astrophys.445:43-49,2006; Erratum-ibid.454:715,2006

Gravitational cooling and density profile near caustics in collisionless dark matter haloes

Authors: Roya Mohayaee, Sergei Shandarin

Abstract: Cold dark matter haloes are populated by high-density structures with sharply-peaked profiles known as caustics which have not yet been resolved by 3-dimensional numerical simulations. Here, we derive semi-analytic expressions for the density profiles near caustics in haloes which form by self-similar accretions of dark matter with infinitesimal velocity dispersion. A simple rescaling shows that… ▽ More Cold dark matter haloes are populated by high-density structures with sharply-peaked profiles known as caustics which have not yet been resolved by 3-dimensional numerical simulations. Here, we derive semi-analytic expressions for the density profiles near caustics in haloes which form by self-similar accretions of dark matter with infinitesimal velocity dispersion. A simple rescaling shows that similarly to the case of absolutely cold medium these profiles are universal: they are valid for all caustics and irrespective of physical parameters of the halo. We derive the maximum density of the caustics and show that it depends on the velocity dispersion and the caustic location. We show that both the absolute and relative thickness of the caustic monotonically decrease toward the center of the halo while the maximum density grows. This indicates that the radial component of the thermal velocities decreases in the inner streams, i.e. the collisionless medium cools down in the radial direction descending to the center of the halo. Finally, we demonstrate that there can be a significant contribution to the emission measure from dark matter particle annihilation in the caustics. △ Less

Submitted 18 October, 2005; v1 submitted 8 March, 2005; originally announced March 2005.

Comments: to appear in MNRAS, minor changes

Reconstruction of primordial density fields

Authors: Roya Mohayaee, Hugues Mathis, Stephane Colombi, Joseph Silk

Abstract: The Monge-Ampere-Kantorovich (MAK) reconstruction is tested against cosmological N-body simulations. Using only the present mass distribution sampled with particles, and the assumption of homogeneity of the primordial distribution, MAK recovers for each particle the non-linear displacement field between its present position and its Lagrangian position on a primordial uniform grid. To test the me… ▽ More The Monge-Ampere-Kantorovich (MAK) reconstruction is tested against cosmological N-body simulations. Using only the present mass distribution sampled with particles, and the assumption of homogeneity of the primordial distribution, MAK recovers for each particle the non-linear displacement field between its present position and its Lagrangian position on a primordial uniform grid. To test the method, we examine a standard LCDM N-body simulation with Gaussian initial conditions and 6 models with non-Gaussian initial conditions: a chi-squared model, a model with primordial voids and four weakly non-Gaussian models. Our extensive analyses of the Gaussian simulation show that the level of accuracy of the reconstruction of the nonlinear displacement field achieved by MAK is unprecedented, at scales as small as about 3 Mpc. In particular, it captures in a nontrivial way the nonlinear contribution from gravitational instability, well beyond the Zel'dovich approximation. This is also confirmed by our analyses of the non-Gaussian samples. Applying the spherical collapse model to the probability distribution function of the divergence of the displacement field, we also show that from a well-reconstructed displacement field, such as that given by MAK, it is possible to accurately disentangle dynamical contributions induced by gravitational clustering from possible initial non-Gaussianities, allowing one to efficiently test the non-Gaussian nature of the primordial fluctuations. In addition, a simple application of MAK using the Zel'dovich approximation allows one to also recover accurately the present-day peculiar velocity field on scales of about 8 Mpc. △ Less

Submitted 14 October, 2005; v1 submitted 12 January, 2005; originally announced January 2005.

Comments: Version to appear in MNRAS, 24 pages, 21 figures appearing (uses 35 figure files), 1 table

Journal ref: Mon.Not.Roy.Astron.Soc. 365 (2006) 939-959