Source author record

Roya Mohayaee

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

ResearcherUnclaimed source record

astro-ph.CO astro-ph.GA hep-ph gr-qc astro-ph astro-ph.HE astro-ph.IM cond-mat hep-ex hep-th math.AP math.NA math.OC nlin.CD Numerical Analysis physics.comp-ph

Catalog footprint

What is connected

21works

16topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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

Building this map preview

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

preprint2023arXiv

Optimal Transport Reconstruction of Baryon Acoustic Oscillations

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

preprint2022arXiv

Accurate Baryon Acoustic Oscillations reconstruction via semi-discrete optimal transport

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.

preprint2022arXiv

Hints of FLRW Breakdown from Supernovae

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

preprint2021arXiv

A Test of the Cosmological Principle with Quasars

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2016arXiv

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

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.

preprint2014arXiv

Unraveling the nature of Gravity through our clumpy Universe

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.

preprint2013arXiv

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

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.

preprint2013arXiv

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

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.

preprint2011arXiv

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

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.

preprint2011arXiv

Probing the anisotropic local universe and beyond with SNe Ia data

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.

preprint2011arXiv

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

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.

preprint2010arXiv

Nonspherical similarity solutions for dark halo formation

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.

preprint2009arXiv

Caustics in growing Cold Dark Matter Haloes

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.

preprint2009arXiv

Hierarchy in the Phase Space and Dark Matter Astronomy

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.

preprint2008arXiv

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

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.

preprint2002arXiv

A reconstruction of the initial conditions of the Universe by optimal mass transportation

Reconstructing the density fluctuations in the early Universe that evolved into the distribution of galaxies we see today is a challenge of modern cosmology [ref.]. An accurate reconstruction would allow us to test cosmological models by simulating the evolution starting from the reconstructed state and comparing it to the observations. Several reconstruction techniques have been proposed [8 refs.], but they all suffer from lack of uniqueness because the velocities of galaxies are usually not known. Here we show that reconstruction can be reduced to a well-determined problem of optimisation, and present a specific algorithm that provides excellent agreement when tested against data from N-body simulations. By applying our algorithm to the new redshift surveys now under way [ref.], we will be able to recover reliably the properties of the primeval fluctuation field of the local Universe and to determine accurately the peculiar velocities (deviations from the Hubble expansion) and the true positions of many more galaxies than is feasible by any other method. A version of the paper with higher-quality figures is available at http://www.obs-nice.fr/etc7/nature.pdf

preprint2001arXiv

Scale Invariance in a Perturbed Einstein-de Sitter Cosmology

This paper seeks to check the validity of the "apparent fractal conjecture" (Ribeiro 2001ab: gr-qc/9909093, astro-ph/0104181), which states that the observed power-law behaviour for the average density of large-scale distribution of galaxies arises when some observational quantities, selected by their relevance in average density profile determination, are calculated along the past light cone. Implementing these conditions in the proposed set of observational relations profoundly changes the behaviour of many observables in the standard cosmological models. In particular, the average density becomes observationally inhomogeneous, even in the spatially homogeneous spacetime of standard cosmology, change which was already analysed by Ribeiro (1992b, 1993, 1994, 1995: astro-ph/9910145) for a non-perturbed model. Here we derive observational relations in a perturbed Einstein-de Sitter cosmology by means of the perturbation scheme proposed by Abdalla and Mohayaee (1999: astro-ph/9810146), where the scale factor is expanded in power series to yield perturbative terms. The differential equations derived in this perturbative context, and other observables necessary in our analysis, are solved numerically. The results show that our perturbed Einstein-de Sitter cosmology can be approximately described by a decaying power-law like average density profile, meaning that the dust distribution of this cosmology has a scaling behaviour compatible with the power-law profile of the density-distance correlation observed in the galaxy catalogues. These results show that, in the context of this work, the apparent fractal conjecture is correct.

Roya Mohayaee

What is connected

Connect this record

See the researcher in context

Building this map preview

21 published item(s)

Optimal Transport Reconstruction of Baryon Acoustic Oscillations

Accurate Baryon Acoustic Oscillations reconstruction via semi-discrete optimal transport

Hints of FLRW Breakdown from Supernovae

A Test of the Cosmological Principle with Quasars

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

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

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

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

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

Unraveling the nature of Gravity through our clumpy Universe

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

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

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

Probing the anisotropic local universe and beyond with SNe Ia data

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

Nonspherical similarity solutions for dark halo formation

Caustics in growing Cold Dark Matter Haloes

Hierarchy in the Phase Space and Dark Matter Astronomy

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

A reconstruction of the initial conditions of the Universe by optimal mass transportation

Scale Invariance in a Perturbed Einstein-de Sitter Cosmology