Cosmological hydrodynamic simulations can accurately predict the properties of the intergalactic medium (IGM), but only under the condition of retaining high spatial resolution necessary to resolve density fluctuations in the IGM. This resolution constraint prohibits simulating large volumes, such as those probed by BOSS and future surveys, like DESI and 4MOST. To overcome this limitation, we present Iteratively Matched Statistics (IMS), a novel method to accurately model the Lyman-alpha forest with collisionless N-body simulations, where the relevant density fluctuations are unresolved. We use a small-box, high-resolution hydrodynamic simulation to obtain the probability distribution function (PDF) and the power spectrum of the real-space Lyman-alpha forest flux. These two statistics are iteratively mapped onto a pseudo-flux field of an N-body simulation, which we construct from the matter density. We demonstrate that our method can perfectly reproduce line-of-sight observables, such as the PDF and power spectrum, and accurately reproduce the 3D flux power spectrum (5-20%). We quantify the performance of the commonly used Gaussian smoothing technique and show that it has significantly lower accuracy (20-80%), especially for N-body simulations with achievable mean inter-particle separations in large-volume simulations. In addition, we show that IMS produces reasonable and smooth spectra, making it a powerful tool for modeling the IGM in large cosmological volumes and for producing realistic "mock" skies for Lyman-alpha forest surveys.
Measuring the absolute scale of the neutrino masses is one of the most exciting opportunities available with near-term cosmological datasets. Two quantities that are sensitive to neutrino mass, scale-dependent halo bias $b(k)$ and the linear growth parameter $f(k)$ inferred from redshift-space distortions, can be measured without cosmic variance. Unlike the amplitude of the matter power spectrum, which always has a finite error, the error on $b(k)$ and $f(k)$ continues to decrease as the number density of tracers increases. This paper presents forecasts for statistics of galaxy and lensing fields that are sensitive to neutrino mass via $b(k)$ and $f(k)$. The constraints on neutrino mass from the auto- and cross-power spectra of spectroscopic and photometric galaxy samples are weakened by scale-dependent bias unless a very high density of tracers is available. In the high density limit, using multiple tracers allows cosmic-variance to be beaten and the forecasted errors on neutrino mass shrink dramatically. In practice, beating the cosmic variance errors on neutrino mass with $b(k)$ will be a challenge, but this signal is nevertheless a new probe of neutrino effects on structure formation that is interesting in its own right.
It has been recently suggested that emerging tension between cosmological parameter values derived in high redshift (CMB anisotropy) and low redshift (cluster counts, Hubble constant) measurements can be reconciled in a model which contains subdominant fraction of dark matter decaying after recombination. We check the model against the CMB Planck data. We find that lensing of CMB anisotropies by the Large Scale Structure gives strong extra constraints on this model, limiting the fraction as F<8% at 2 sigma confidence level. However, investigating the combined data set of CMB and conflicting low z measurements, we obtain that the model with F=2-5 % exhibits better fit (by 1.5-3 sigma depending on the lensing priors) as compared to that of the concordance LambdaCDM cosmological model.
The dispersion measure of extragalactic radio transients, such as of recently discovered Fast Radio Burst FRB150418, can be used to measure the column density of free electrons in the intergalactic medium. The same electrons also scatter the Cosmic Microwave Background (CMB) photons, affecting precision measurements of cosmological parameters. We explore the connection between the dispersion measure of radio transients existing during the Epoch of Reionization (EoR) and the total optical depth for the CMB, $\tau_{CMB}$, showing that the existence of such transients would provide a new sensitive probe of $\tau_{CMB}$. As an example, we consider the population of FRBs. Assuming they exist during the EoR, we show that: (i) such sources can probe the reionization history by measuring $\tau_{CMB}$ to sub-percent accuracy, and (ii) they can be detected with high significance by an instrument such as the Square Kilometer Array.
We use the effective field theory of dark energy (EFT of DE) formalism to constrain dark energy models belonging to the Horndeski class with the recent Planck 2015 CMB data. The space of theories is spanned by a certain number of parameters determining the linear cosmological perturbations, while the expansion history is set to that of a standard $\Lambda$CDM model. We always demand that the theories be free of fatal instabilities. Additionally, we consider two optional conditions, namely that scalar and tensor perturbations propagate with subliminal speed. Such criteria severely restrict the allowed parameter space and are thus very effective in shaping the posteriors. As a result, we confirm that no theory performs better than $\Lambda$CDM when CMB data alone are analysed. Indeed, the healthy dark energy models considered here are not able to reproduce those phenomenological behaviours of the effective Newton constant and gravitational slip parameters that, according to previous studies, best fit the data.
The brightest cluster galaxy (BCG) in the majority of relaxed, cool core galaxy clusters is radio loud, showing non-thermal radio jets and lobes ejected by the central active galactic nucleus (AGN). Such relativistic plasma has been unambiguously shown to interact with the surrounding thermal intra-cluster medium (ICM) thanks to spectacular images where the lobe radio emission is observed to fill the cavities in the X-ray-emitting gas. This `radio-mode AGN feedback' phenomenon, which is thought to quench cooling flows, is widespread and is critical to understand the physics of the inner regions of galaxy clusters and the properties of the central BCG. At the same time, mechanically-powerful AGN are likely to drive turbulence in the central ICM which may contribute to gas heating and also play a role for the origin of non-thermal emission on cluster-scales. Diffuse non-thermal emission has been observed in a number of cool core clusters in the form of a radio mini-halo surrounding the radio-loud BCG on scales comparable to that of the cooling region. This contribution outlines the main points covered by the talk on these topics. In particular, after summarizing the cooling flow regulation by AGN heating and the non-thermal emission from cool core clusters, we present a recent study of the largest collection of known mini-halo clusters (~ 20 objects) which investigated the scenario of a common origin of radio mini-halos and gas heating. We further discuss the prospects offered by future radio surveys with the Square Kilometre Array (SKA) for building large (>> 100 objects), unbiased mini-halo samples while probing at the same time the presence of radio-AGN feedback in the host clusters.
Intrinsic galaxy shape and angular momentum alignments can arise in cosmological large-scale structure due to tidal interactions or galaxy formation processes. Cosmological hydrodynamical simulations have recently come of age as a tool to study these alignments and their contamination to weak gravitational lensing. We probe the redshift and luminosity evolution of intrinsic alignments in Horizon-AGN between $z=0$ and $z=3$ for galaxies with an $r$-band absolute magnitude of $M_r\leq-20$. Alignments transition from being radial at low redshifts and high luminosities, dominated by the contribution of ellipticals, to being tangential at high redshift and low luminosities, where discs dominate the signal. This cannot be explained by the evolution of the fraction of ellipticals and discs alone: intrinsic evolution in the amplitude of alignments is necessary. We constrain the evolution of the alignment amplitude as a function of luminosity for elliptical galaxies alone and find it to be in good agreement with current observations and the nonlinear tidal alignment model at projected separations of $\gtrsim 1$ Mpc. Alignments of discs are null in projection and consistent with current low redshift observations. The combination of the two populations yields an overall amplitude to be a factor of $\simeq 2$ lower than observed alignments of luminous red galaxies with a steeper luminosity dependence. The restriction on accurate galaxy shapes implies that the galaxy population in the simulation is complete only to $M_r\leq-20$. Higher resolution simulations will be necessary to avoid extrapolation of the intrinsic alignment predictions to the range of luminosities probed by future surveys.
We analyse solutions of the MHD equations around the electroweak transition taking into account the effects of the chiral anomaly. It is shown that a transition that is not of the first order has direct consequences on the evolution of the asymmetry between left- and right-handed leptons. Assuming an initial chiral asymmetry in the symmetric phase at temperatures higher than the transition temperature, as well as the existence of magnetic fields, it is demonstrated that the asymmetry typically grows with time, until it undergoes a fast decrease at the transition, and then eventually gets damped at lower temperatures in the broken phase. We argue that it is unlikely to have any significant magnetic field amplification as a consequence of the electroweak transition in the Standard model, even when the chiral anomaly is introduced. The presence of a chiral asymmetry between left- and right-handed charge carriers naturally leads to the creation of helical magnetic fields from non-helical fields and this can have consequences on their subsequent evolution.
Recent constrains on the sum of neutrino masses inferred by analyzing cosmological data, show that detecting a non-zero neutrino mass is within reach of forthcoming cosmological surveys, implying a direct determination of the absolute neutrino mass scale. The measurement relies on constraining the shape of the matter power spectrum below the neutrino free streaming scale: massive neutrinos erase power at these scales. Detection of a lack of small-scale power, however, could also be due to a host of other effects. It is therefore of paramount importance to validate neutrinos as the source of power suppression at small scales. We show that, independent on hierarchy, neutrinos always show a footprint on large, linear scales; the exact location and properties can be related to the measured power suppression (an astrophysical measurement) and atmospheric neutrinos mass splitting (a neutrino oscillation experiment measurement). This feature can not be easily mimicked by systematic uncertainties or modifications in the cosmological model. The measurement of such a feature, up to 1% relative change in the power spectrum, is a smoking gun for confirming the determination of the absolute neutrino mass scale from cosmological observations. It also demonstrates the synergy of astrophysics and particle physics experiments.
Recently it has been pointed out that a cosmic background of relativistic axion-like particles (ALPs) would be produced by the primordial decays of heavy fields in the post-inflation epoch, contributing to the extra-radiation content in the Universe today. Primordial magnetic fields would trigger conversions of these ALPs into sub-MeV photons during the dark ages. This photon flux would produce an early reionization of the Universe, leaving a significant imprint on the total optical depth to recombination $\tau$. Using the current measurement of $\tau$ and the limit on the extra-radiation content $\Delta N_{\rm eff} $ by the Planck experiment we put a strong bound on the ALP-photon conversions. Namely we obtain upper limits on the product of the photon-ALP coupling constant $g_{a\gamma}$ times the magnetic field strength $B$ down to $g_{a\gamma} B \gtrsim 6 \times 10^{-18} \textrm{GeV}^{-1} \textrm{nG} $ for ultralight ALPs.
Keane et al. (2016) have recently claimed to have obtained the first precise localization for a Fast Radio Burst thanks to the identification of a contemporaneous fading slow (~week-timescale) radio transient. We show that the quiescent radio luminosity of the proposed host galaxy points to the presence of an AGN, and therefore that the claimed transient may instead represent common AGN variability. We further show that the expected number of variable (rather than transient) sources in the Parkes localization region of FRB 150418 is order unity. Finally we show that the properties of the radio counterpart are incompatible with a synchrotron-emitting blastwave. Taken together, these results indicate that the claimed radio source is unlikely to be associated with FRB 150418, and hence that a precise localization and redshift determination cannot be justified.
We apply the delay in timing of FERMI GMB transient occurred in coincidence with gravitational waves event GW150914 observed by LIGO to constrain the size of the spherical brane-universe expanding in multi-dimensional space-time. A bound on spatial curvature of the brane is obtained.
We investigate how a Higgs mechanism could be responsible for the emergence of gravity in extensions of Einstein theory. In this scenario, at high energies, symmetry restoration could "turn off" gravity, with dramatic implications for cosmology and quantum gravity. The sense in which gravity is muted depends on the details of the implementation. In the most extreme case gravity's dynamical degrees of freedom would only be unleashed after the Higgs field acquires a non-trivial vacuum expectation value, with gravity reduced to a topological field theory in the symmetric phase. We might also identify the Higgs and the Brans-Dicke fields in such a way that in the unbroken phase Newton's constant vanishes, decoupling matter and gravity. We discuss the broad implications of these scenarios.
For particle physics observables at colliders such as the LHC at CERN, it has been common practice for many decades to estimate the theoretical uncertainty by studying the variations of the predicted cross sections with a priori unpredictable scales. In astroparticle physics, this has so far not been possible, since most of the observables were calculated at Born level only, so that the renormalization scheme and scale dependence could not be studied in a meaningful way. In this paper, we present the first quantitative study of the theoretical uncertainty of the neutralino dark matter relic density from scheme and scale variations. We first explain in detail how the renormalization scale enters the tree-level calculations through coupling constants, masses and mixing angles. We then demonstrate a reduction of the renormalization scale dependence through one-loop SUSY-QCD corrections in many different dark matter annihilation channels and enhanced perturbative stability of a mixed on-shell/$\bar{\rm DR}$ renormalization scheme over a pure $\bar{\rm DR}$ scheme in the top-quark sector. In the stop-stop annihilation channel, the Sommerfeld enhancement and its scale dependence are shown to be of particular importance. Finally, the impact of our higher-order SUSY-QCD corrections and their scale uncertainties are studied in three typical scenarios of the phenomenological Minimal Supersymmetric Standard Model with eleven parameters (pMSSM-11). We find that the theoretical uncertainty is reduced in many cases and can become comparable to the size of the experimental one in some scenarios.
We give, in some detail, a critical overview over recent work towards deriving a cosmological phenomenology from the fundamental quantum dynamics of group field theory (GFT), based on the picture of a macroscopic universe as a "condensate" of a large number of quanta of geometry which are given by excitations of the GFT field over a "no-space" vacuum. We emphasise conceptual foundations, relations to other research programmes in GFT and the wider context of loop quantum gravity (LQG), and connections to the quantum physics of real Bose-Einstein condensates. We show how to extract an effective dynamics for GFT condensates from the microscopic GFT physics, and how to compare it with predictions of more conventional quantum cosmology models, in particular loop quantum cosmology (LQC). No detailed familiarity with the GFT formalism is assumed.
We study the evolution of cosmological perturbations in an anti-de-Sitter (AdS) bulk through a cosmological singularity by mapping the dynamics onto the boundary conformal fields theory by means of the AdS/CFT correspondence. We consider a deformed AdS space-time obtained by considering a time-dependent dilaton which induces a curvature singularity in the bulk at a time which we call $t = 0$, and which asymptotically approaches AdS both for large positive and negative times. The boundary field theory becomes free when the bulk curvature goes to infinity. Hence, the evolution of the fluctuations is under better controle on the boundary than in the bulk. To avoid unbounded particle production across the bounce it is necessary to smooth out the curvature singularity at very high curvatures. We show how the bulk cosmological perturbations can be mapped onto boundary gauge field fluctuations. We evolve the latter and compare the spectrum of fluctuations on the infrared scales relevant for cosmological observations before and after the bounce point. We find that the index of the power spectrum of fluctuations is the same before and after the bounce.
The recent discovery of unexplained X-ray line of $3.5-3.6$ keV emitted from the Perseus cluster of galaxies and M31 and the excess X-ray line of $8.7$ keV emitted from the Milky Way center may indicate that dark matter would decay. In this article, I show that approximately 80 \% of dark matter being 7.1 keV sterile neutrinos and 20 \% of dark matter being 17.4 keV sterile neutrinos can satisfactorily explain the observed X-ray lines and account for all missing mass. No free parameter is needed in this model. This scenario is also compatible with current robust observational constraints from the matter power spectrum in large-scale structures and would alleviate the challenges faced by the existing dark matter models.
We study the inner dynamics of accreting Eddington-inspired Born-Infeld black holes using the homogeneous approximation and taking charge as a surrogate for angular momentum. We show that there is a minimum of the accretion rate below which mass inflation does not occur, and we derive an analytical expression for this threshold as a function of the fundamental scale of the theory, the accretion rate, the mass, and the charge of the black hole. Our result explicitly demonstrates that, no matter how close Eddington-inspired Born-Infeld gravity is to general relativity, there is always a minimum accretion rate below which there is no mass inflation. For larger accretion rates, mass inflation takes place inside the black hole as in general relativity until the extremely rapid density variations bring it to an abrupt end. We derive analytical scaling solutions for the value of the energy density and of the Misner-Sharp mass attained at the end of mass inflation as a function of fundamental scale of the theory, the accretion rate, the mass, and the charge of the black hole, and compare these with the corresponding numerical solutions. We find that, except for unreasonably high accretion rates, our analytical results appear to provide an accurate description of homogeneous mass inflation inside accreting Eddington-inspired Born-Infeld black holes.
We present an exact expression for the $1/f$ contribution to the noise of the CMB temperature and polarization maps for a survey in which the scan pattern is isotropic. The result for polarization applies likewise to surveys with and without a rotating half-wave plate. A representative range of survey parameters is explored and implications for the design and optimization of future surveys are discussed. These results are most directly applicable to space-based surveys, which afford considerable freedom in the choice of the scan pattern on the celestial sphere. We discuss the applicability of the methods developed here to analyzing past experiments and present some conclusions pertinent to the design of future experiments. The techniques developed here do not require that the excess low frequency noise have exactly the $1/f$ shape and readily generalize to other functional forms for the detector noise power spectrum. In the case of weakly anisotropic scanning patterns the techniques in this paper can be used to find an preconditioner for solving the map making equation efficiently using the conjugate gradient method.
We consider all degenerate scalar-tensor theories that depend quadratically on second order derivatives of a scalar field, which we have identified in a previous work. These theories, whose degeneracy in general ensures the absence of Ostrogradski instability, include the quartic Horndenski Lagrangian as well as its quartic extension beyond Horndeski, but also other families of Lagrangians. We study how all these theories transform under general conformal-disformal transformations and find that they can be separated into three main classes that are stable under these transformations. This leads to a complete classification modulo conformal-disformal transformations. Finally, we show that these higher order theories include mimetic gravity and some particular khronometric theories. They also contain theories that do not correspond, to our knowledge, to already studied theories, even up to field redefinitions.
We treat quantum creation of gravitons by small scale factor oscillations around the average of an expanding universe. Such oscillations are predicted both by modified gravity theories with a term proportional to the square of the Ricci scalar in the gravitational action, and by semiclassical gravity in which a renormalized expectation value of a quantum matter stress tensor is the source of gravity. A perturbative method due to Birrell and Davies involving an expansion in a conformal coupling parameter is used to calculate the number density and energy density of the created gravitons. Cosmological constraints on the present graviton energy density are discussed, which indicate an upper bound of the order of $4\%$ of the total mass-energy density of the universe. This constraint leads to an upper bound on the dimensionless amplitude of the oscillations. We also discuss decoherence of quantum systems produced by the spacetime geometry fluctuations due to such a graviton bath, and find a lower bound on the decoherence time resulting from this process.
High-energy neutrino (HEN) and gravitational wave (GW) can probe astrophysical sources in addition to electromagnetic observations. Multimessenger studies can reveal nature of the sources which may not be discerned from one type of signal alone. We discuss HEN emission in connection with the Advanced Laser Interferometer Gravitational-wave Observatory (ALIGO) event GW150914 which could be associated with a short gamma-ray burst (GRB) detected by the $Fermi$ Gamma-ray Burst Monitor (GBM) 0.4 s after the GW event and within localization uncertainty of the GW event. We calculate HEN flux from this short GRB, GW150914-GBM, and show that non-detection of a high-energy starting event (HESE) by the IceCube Neutrino Observatory can constrain the total isotropic-equivalent jet energy of this short burst to be less than $3\times 10^{52}$ erg.
In this article we will construct the most general torsion-free
parity-invariant covariant theory of gravity that is free from ghost-like and
tachyonic nstabilities around constant curvature space-times in four
dimensions. Specifically, this includes the Minkowski, de Sitter and anti-de
Sitter backgrounds. We will first argue in details how starting from a general
covariant action for the metric one arrives at an "equivalent" action that at
most contains terms that are quadratic in curvatures but nevertheless is
sufficient for the purpose of studying stability of the original action. We
will then briefly discuss how such a "quadratic curvature action" can be
decomposed in a covariant formalism into separate sectors involving the tensor,
vector and scalar modes of the metric tensor; most of the details of the
analysis however, will be presented in an accompanying paper. We will find that
only the transverse and trace-less spin-2 graviton with its two helicity states
and possibly a spin-0 Brans-Dicke type scalar degree of freedom are left to
propagate in 4 dimensions. This will also enable us to arrive at the
consistency conditions required to make the theory perturbatively stable around
constant curvature backgrounds.
This will be included as a chapter in the book entitled "At the Frontier of
Spacetime - Scalar-Tensor Theory, Bells Inequality, Machs Principle, Exotic
Smoothness" (Springer 2016)
Links to: arXiv, form interface, find, astro-ph, recent, 1602, contact, help (Access key information)
This work describes the implementation and application of a correlation determination method based on Self Organizing Maps and Bayesian Inference (SOMBI). SOMBI aims to automatically identify relations between different observed parameters in unstructured cosmological or astrophysical surveys by automatically identifying data clusters in high-dimensional datasets via the Self Organizing Map neural network algorithm. Parameter relations are then revealed by means of a Bayesian inference within respective identified data clusters. Specifically such relations are assumed to be parametrized as a polynomial of unknown order. The Bayesian approach results in a posterior probability distribution function for respective polynomial coefficients. To decide which polynomial order suffices to describe correlation structures in data, we include a method for model selection, the Bayesian Information Criterion, to the analysis. The performance of the SOMBI algorithm is tested with mock data. As illustration we also provide applications of our method to cosmological data. In particular, we present results of a correlation analysis between galaxy and AGN properties provided by the SDSS catalog with the cosmic large-scale-structure.
Superclusters are the largest relatively isolated systems in the cosmic web.
Using the SDSS BOSS survey we search for the largest superclusters in the
redshift range $0.43<z<0.71$.
We generate a luminosity-density field smoothed over $8 h^{-1}\mathrm{Mpc}$
to detect the large-scale over-density regions. Each individual over-density
region is defined as single supercluster in the survey. We define the
superclusters in the way that they are comparable with the superclusters found
in the SDSS main survey.
We found a system we call the BOSS Great Wall (BGW), which consists of two
walls with diameters 186 and 173 $h^{-1}$Mpc, and two other major superclusters
with diameters of 64 and 91 $h^{-1}$Mpc. As a whole, this system consists of
830 galaxies with the mean redshift 0.47. We estimate the total mass to be
approximately $2\times10^{17}h^{-1}M_\odot$. The morphology of the
superclusters in the BGW system is similar to the morphology of the
superclusters in the Sloan Great Wall region.
The BGW is one of the most extended and massive system of superclusters yet
found in the Universe.
It is common practice in cosmology to model large-scale structure observables as lognormal random fields, and this approach has been successfully applied in the past to the matter density and weak lensing convergence fields separately. We argue that this approach has fundamental limitations which prevent its use for jointly modelling these two fields since the lognormal distribution's shape can prevent certain correlations to be attainable. Given the need of ongoing and future large-scale structure surveys for fast joint simulations of clustering and weak lensing, we propose two ways of overcoming these limitations. The first approach slightly distorts the power spectra of the fields using one of two algorithms that minimises either the absolute or the fractional distortions. The second one is by obtaining more accurate convergence marginal distributions, for which we provide a fitting function, by integrating the lognormal density along the line of sight. The latter approach also provides a way to determine directly from theory the skewness of the convergence distribution and, therefore, the parameters for a lognormal fit. We present the public code Full-sky Lognormal Astro-fields Simulation Kit (FLASK) which can make tomographic realisations on the sphere of an arbitrary number of correlated lognormal or Gaussian random fields by applying either of the two proposed solutions, and show that it can create joint simulations of clustering and lensing with sub-per-cent accuracy over relevant angular scales and redshift ranges.
We study evolution of voids in cosmological simulations using a new method for tracing voids over cosmic time. The method is based on tracking watershed basins (contiguous regions around density minima) of well developed voids at low redshift, on a regular grid of density field. It enables us to construct a robust and continuous mapping between voids at different redshifts, from initial conditions to the present time. We discuss how the new approach eliminates strong spurious effects of numerical origin when voids evolution is traced by matching voids between successive snapshots (by analogy to halo merger trees). We apply the new method to a cosmological simulation of a standard LambdaCDM cosmological model and study evolution of basic properties of typical voids (with effective radii between 6Mpc/h and 20Mpc/h at redshift z=0) such as volumes, shapes, matter density distributions and relative alignments. The final voids at low redshifts appear to retain a significant part of the configuration acquired in initial conditions. Shapes of voids evolve in a collective way which barely modifies the overall distribution of the axial ratios. The evolution appears to have a weak impact on mutual alignments of voids implying that the present state is in large part set up by the primordial density field. We present evolution of dark matter density profiles computed on iso-density surfaces which comply with the actual shapes of voids. Unlike spherical density profiles, this approach enables us to demonstrate development of theoretically predicted bucket-like shape of the final density profiles indicating a wide flat core and a sharp transition to high-density walls voids.
The angular correlation function is a powerful tool for deriving the clustering properties of AGN and hence the mass of the corresponding dark matter halos in which they reside. However, studies based on the application of the angular correlation function on X-ray samples, yield results apparently inconsistent with those based on the direct estimation of the spatial correlation function. The goal of the present paper is to attempt to investigate this issue by analysing a well defined sample. To this end we use the hard-band (2-10 keV) X-ray selected sources of the Chandra AEGIS fields, chosen because of the availability of accurately derived flux sensitivity maps. In particular we use the 186 hard-band sources with spectroscopic redshifts in the range z=0.3-1.3, a range selected in order to contain the bulk of the AGN while minimizing the contribution of unknown clustering and luminosity evolution from very high redshifts. Using the projected spatial auto-correlation function, we derive a clustering comoving length of 5.4+-1.0 Mpc (for gamma=1.8), consistent with results in the literature. We further derive the angular correlation function and the corresponding spatial clustering length using the Limber's inversion equation and a novel parametrization of the clustering evolution model that also takes into account the bias evolution of the host dark matter halo. The Limber's inverted spatial comoving clustering length of 5.5+-1.2 Mpc at a median redshift of z~0.75, matches the directly measured one, from the spatial correlation function analysis, but for a significant non-linear contribution to the growing mode of perturbations, estimated independently from literature results of x_0 at different redshifts. Therefore, using this sample of X-ray AGN and our clustering evolution parametrization we have found an excellent consistency between the angular and spatial clustering analysis.
B and L violating interactions of ordinary particles with their twin particles from hypothetical mirror world can co-generate baryon asymmetries in both worlds in comparable amounts, $\Omega'_B/\Omega_B \sim 5$ or so. On the other hand, the same interactions induce the oscillation phenomena between the neutral particles of two sectors which convert e.g. mirror neutrons into our antineutrons. These oscillations are environment dependent and can have fascinating physical consequences.
The viscous inhomogeneities of a relativistic plasma determine a further class of entropic modes whose amplitude must be sufficiently small since curvature perturbations are observed to be predominantly adiabatic and Gaussian over large scales. When the viscous coefficients only depend on the energy density of the fluid the corresponding curvature fluctuations are shown to be almost adiabatic. After addressing the problem in a gauge-invariant perturbative expansion, the same analysis is repeated at a non-perturbative level by investigating the nonlinear curvature inhomogeneities induced by the spatial variation of the viscous coefficients. It is demonstrated that the quasiadiabatic modes are suppressed in comparison with a bona fide adiabatic solution. Because of its anomalously large tensor to scalar ratio the quasiadiabatic mode cannot be a substitute for the conventional adiabatic paradigm so that, ultimately, the present findings seems to exclude the possibility of a successful accelerated dynamics solely based on relativistic viscous fluids. If the dominant adiabatic mode is not affected by the viscosity of the background a sufficiently small fraction of entropic fluctuations of viscous origin cannot be a priori ruled out.
In this study, we obtain the size distribution of voids as a 3-parameter redshift independent log-normal void probability function (VPF) directly from the Cosmic Void Catalog (CVC). Although many statistical models of void distributions are based on the counts in randomly placed cells, the log-normal VPF that we here obtain is independent of the shape of the voids due to the parameter-free void finder of the CVC. We use three void populations drawn from the CVC generated by the Halo Occupation Distribution (HOD) Mocks which are tuned to three mock SDSS samples to investigate the void distribution statistically and the effects of the environments on the size distribution. As a result, it is shown that void size distributions obtained from the HOD Mock samples are satisfied by the 3-parameter log-normal distribution. In addition, we find that there may be a relation between hierarchical formation, skewness and kurtosis of the log-normal distribution for each catalog. We also show that the shape of the 3-parameter distribution from the samples is strikingly similar to the galaxy log-normal mass distribution obtained from numerical studies. This similarity of void size and galaxy mass distributions may possibly indicate evidence of nonlinear mechanisms affecting both voids and galaxies, such as large scale accretion and tidal effects. Considering in this study all voids are generated by galaxy mocks and show hierarchical structures in different levels, it may be possible that the same nonlinear mechanisms of mass distribution affect the void size distribution.
We propose the use of the kinetic Sunyaev-Zel'dovich (kSZ) effect to probe the circumgalactic medium (CGM), with the aid of a spectroscopic survey covering the same area of a SZ survey. One can design an optimal estimator of the kSZ effect of the CGM with a matched filter, and construct the cross correlation between the estimator and the peculiar velocity recovered from the galaxy survey, which can be measured by stacking a number of galaxies. We investigate two compelling profiles for the CGM, the MB profile (Maller & Bullock 2004) and the $\beta$ profile, and estimate the detectability against the synergy of a fiducial galaxy survey with number density $10^{-3}h^3\,$ Mpc$^{-3}$ and an ACT-like SZ survey. We show that the shape of the filter does not change much with redshift for the $\beta$ profile, while there are significant side lobes at $z<0.1$ for the MB profile. By stacking $\sim 10^4$ Milky Way-size halos around z $\sim 0.5$, one can get $\gtrsim$ 1 $\sigma$ signal to noise (S/N) for the both profiles. The S/N increases with decreasing redshift before it reaches a maximum ($\sim$ 7.5 at z $\simeq$ 0.15 for the MB profile, $\sim 19$ at $z\simeq 0.03$ for the $\beta$ profile). Due to the large beam size, a Planck-like CMB survey can marginally detect the kSZ signal by stacking the same number of galaxies at $z<0.1$. The search for the CGM in realistic surveys will involve dividing the galaxies into subsamples with similar redshift and mass of host halos, and scaling the results presented here to obtain the S/N.
The BOSS quasar sample is used to study cosmic homogeneity with a 3D survey over the largest volume used to date, 14 $h^{-3}$ Gpc$^3$. The selected sample includes 47,858 quasars over 5983 deg$^2$ and in the range $2.2<z<2.8$. We measure the count-in-sphere $\mathcal{N}(<\! r)$, i.e., the average number of objects around a given object, and its logarithmic derivative, the fractal correlation dimension $D_2(r)$. For a homogeneous distribution $\mathcal{N}(<\! r) \propto r^3$ and $D_2(r)=3$. With the simplest $DD/RR$ estimator, these behaviors can be checked without determining an average density, which requires homogeneity to be defined. We observe the predicted behavior over one decade in the radius $r$. In particular, the fractal correlation dimension $D_2$ is found to be compatible with three on large scales at a high accuracy: $3-\langle D_2 \rangle < 6 \times 10^{-4}$ (2 $\sigma$) over the range $250<r<1200 \, h^{-1}$Mpc. This result establishes homogeneity on large scales, therefore allowing us to define a reliable average density. We then introduce a more accurate estimator which requires a crude estimate of the average density, and find that the transition to homogeneity, described by $\mathcal{N}(<\! r)$ and $D_2(r)$, quantitatively agrees with the $\Lambda$CDM prediction and that $D_2$ is still compatible with three at an even higher accuracy: $3-\langle D_2 \rangle < 1.5 \times 10^{-4}$ (2 $\sigma$) over the range $250<r<1200 \, h^{-1}$Mpc.
Pre-recombination acoustic oscillations induce an initial relative velocity $v_{cb}$ between baryons and dark matter. We show that the leading effect on galaxy clustering at lower redshifts is induced by its divergence $\theta_{cb} = \partial_i v_{cb}^i$, and estimate the magnitude of this new effect through a spherical collapse calculation. We then derive all streaming velocity contributions to the galaxy power spectrum at 1-loop order, leading to several new terms. Including all these contributions will be essential to avoid a bias in future efforts to use the baryon acoustic oscillation feature in galaxy clustering as standard ruler.
The Starobinsky model is one of the inflation models consistent with the result of CMB observation by the Planck satellite. We consider the dynamics of the Starobinsky inflation in the presence of another scalar field with a large expectation value during inflation due to a negative Hubble-induced mass. We find that it would not be affected unless the other field has an amplitude very close to the Planck scale, and, if so, the spectral index of curvature perturbation would deviate from the Planck result.
In this letter, we consider a class of scenarios in which the dark matter is part of a heavy hidden sector that is thermally decoupled from the Standard Model in the early universe. The dark matter freezes-out by annihilating to a lighter, metastable state, whose subsequent abundance can naturally come to dominate the energy density of the universe. When this state decays, it reheats the visible sector and dilutes all relic abundances, thereby allowing the dark matter to be orders of magnitude heavier than the weak scale. For concreteness, we consider a simple realization with a Dirac fermion dark matter candidate coupled to a massive gauge boson that decays to the Standard Model through its kinetic mixing with hypercharge. We identify viable parameter space in which the dark matter can be as heavy as ~1-100 PeV without being overproduced in the early universe.
We investigate the mass content of galaxies in the core of the galaxy cluster Abell 611. We perform a strong lensing analysis of the cluster core and use velocity dispersion measurements for individual cluster members as additional constraints. Despite the small number of multiply-imaged systems and cluster members with central velocity dispersions available in the core of A611, the addition of velocity dispersion measurements leads to tighter constraints on the mass associated with the galaxy component, and as a result, on the mass associated with the dark matter halo. Without the spectroscopic velocity dispersions, we would overestimate the mass of the galaxy component by a factor of $\sim1.5$, or, equivalently, we would underestimate the mass of the cluster dark halo by $\sim5\%$. We perform an additional lensing analysis using surface brightness (SB) reconstruction of the tangential giant arc. This approach improves the constraints on the mass parameters of the 5 galaxies close to the arc by up to a factor $\sim10$. The galaxy velocity dispersions resulting from the SB analysis are consistent at the 1$\sigma$ confidence level with the spectroscopic measurements and with the prediction from the simple pointlike analysis. In contrast the truncation radii for 2-3 galaxies depart significantly from the galaxy scaling relation and suggest differences in the stripping history from galaxy to galaxy.
Gamma-ray lines from dark matter annihilation ($\chi\chi\to \gamma X$, where $X=\gamma,h,Z$) are always accompanied, at lower energies, by a continuum gamma-ray spectrum stemming both from radiative corrections ($X=\gamma$) and from the decay debris of the second particle possibly present in the final state ($X=h,Z$). This model-independent gamma-ray emission can be exploited to derive novel limits on gamma-ray lines that do not rely on the line-feature. Although such limits are not expected to be as stringent, they can be used to probe the existence of $\gamma$-ray lines for dark matter masses beyond the largest energies accessible to current telescopes. Here, we use continuous gamma-ray searches from Fermi-LAT observations of Milky Way dwarf spheroidal galaxies and from H.E.S.S. observations of the Galactic Halo to extend the limits on the annihilation cross sections into monochromatic photons to dark matter masses well beyond 500 GeV (Fermi-LAT) and 20 TeV (H.E.S.S.). In this large mass regime, our results provide the first constraints on $\gamma$-ray lines from dark matter annihilation.
Based on the Lema\^itre-Tolman-Bondi (LTB) metric we consider two flat inhomogeneous big-bang models. We aim at clarifying, as far as possible analytically, basic features of the dynamics of the simplest inhomogeneous models and to point out the potential usefulness of exact inhomogeneous solutions as generalizations of the homogeneous configurations of the cosmological standard model. We discuss explicitly partial successes but also potential pitfalls of these simplest models. Although primarily seen as toy models, the relevant free parameters are fixed by best-fit values using the Joint Light-curve Analysis (JLA)-sample data. On the basis of a likelihood analysis we find that a local hump provides a better description of the observations than a local void. Future redshift-drift measurements are discussed as a promising tool to discriminate between inhomogeneous configurations and the $\Lambda$CDM model.
We propose a novel framework for asymmetric scalar dark matter (ADM), which has interesting collider phenomenology in terms of an unstable ADM bound state (ADMonium) produced via Higgs portals. ADMonium is a natural consequence of the basic features of ADM: the (complex scalar) ADM is charged under a dark local $U(1)_d$ symmetry which is broken at a low scale and provides a light gauge boson $X$. The dark gauge coupling is strong and then ADM can annihilate away into $X$-pair effectively. Therefore, the ADM can form bound state due to its large self-interaction via $X$ mediation. To explore the collider signature of ADMonium, we propose that ADM has a two-Higgs doublet portal. The ADMonium can have a sizable mixing with the heavier Higgs boson, which admits a large cross section of ADMonium production associated with $b\bar b$. Of particular interest, our setup nicely explains the recent di-photon anomaly at 750 GeV via the events from ${\rm ADMonium}\ra 2X(\ra e^+e^-)$, where the electrons are identified as (un)converted photons. The prediction of an ADM near 375 GeV can be tested at the direct detection experiments in the near future.
We investigate the inflationary universe in a theory where two scalar fields non-minimally coupling to the scalar curvature and an extra $R^2$ term exist and the conformal invariance is broken. In particular, the slow-roll inflation is explored for the case that one scalar field is dynamical and that two scalar fields are dynamical. As a result, we show that the spectral index of the curvature perturbations and the tensor-to-scalar ratio of the density perturbations can be compatible with the Planck results. It is also demonstrated that the graceful exit from inflation can be realized.
Pure disk galaxies without any bulge component, i.e., neither classical nor
pseudo, seem to have escaped the affects of merger activity inherent to
hierarchical galaxy formation models as well as strong internal secular
evolution.
We discover that a significant fraction (15 - 18 %) of disk galaxies in the
Hubble Deep Field (0.4 < z < 1.0) as well as in the local Universe (0.02 < z <
0.05) are such pure disk systems (hereafter, PDS). The spatial distribution of
light in these PDS is well described by a single exponential function from the
outskirts to the centre and appears to have remained intact over the last 8
billion years keeping the mean central surface brightness and scale-length
nearly constant. These two disk parameters of PDS are brighter and shorter,
respectively, than of those disks which are part of disk galaxies with bulges.
Since the fraction of PDS as well as their profile defining parameters do not
change, it indicates that these galaxies have not witnessed either major
mergers or multiple minor mergers since z~1. However, there is substantial
increase in their total stellar mass and total size over the same time range.
This suggests that smooth accretion of cold gas via cosmic filaments is the
most probable mode of their evolution. We speculate that PDS are dynamically
hotter and cushioned in massive dark matter halos which may prevent them from
undergoing strong secular evolution.
In astrophysics, the two main methods traditionally in use for solving the Euler equations of ideal fluid dynamics are smoothed particle hydrodynamics and finite volume discretization on a stationary mesh. However, the goal to efficiently make use of future exascale machines with their ever higher degree of parallel concurrency motivates the search for more efficient and more accurate techniques for computing hydrodynamics. Discontinuous Galerkin (DG) methods represent a promising class of methods in this regard, as they can be straightforwardly extended to arbitrarily high order while requiring only small stencils. Especially for applications involving comparatively smooth problems, higher-order approaches promise significant gains in computational speed for reaching a desired target accuracy. Here, we introduce our new astrophysical DG code TENET designed for applications in cosmology, and discuss our first results for 3D simulations of subsonic turbulence. We show that our new DG implementation provides accurate results for subsonic turbulence, at considerably reduced computational cost compared with traditional finite volume methods. In particular, we find that DG needs about 1.8 times fewer degrees of freedom to achieve the same accuracy and at the same time is more than 1.5 times faster, confirming its substantial promise for astrophysical applications.
In this article we revisit the significance of the often debated structural similarity between the equations of electromagnetism and fluid dynamics. Although the matching of the two sets of equations has successfully been done for non-dissipative forms of the equations, little has been done for cases where the dissipative terms are non-negligible. We consider the consequence of non-negligible viscosity and diffusivity, and how the fine-tuning of these parameters could allow fluid dynamics to be used to indirectly study certain properties of magnetic fields.
We discuss our current understanding of the nature of the faint, high-frequency radio sky. The Tenth Cambridge (10C) survey at 15.7 GHz is the deepest high-frequency radio survey to date, covering 12 square degrees to a completeness limit of 0.5 mJy, making it the ideal starting point from which to study this population. In this work we have matched the 10C survey to several lower-frequency radio catalogues and a wide range of multi-wavelength data (near- and far-infrared, optical and X-ray). We find a significant increase in the proportion of flat-spectrum sources at flux densities below 1 mJy - the median radio spectral index between 15.7 GHz and 610 MHz changes from 0.75 for flux densities greater than 1.5 mJy to 0.08 for flux densities less than 0.8 mJy. The multi-wavelength analysis shows that the vast majority (> 94 percent) of the 10C sources are radio galaxies; it is therefore likely that these faint, flat spectrum sources are a result of the cores of radio galaxies becoming dominant at high frequencies. We have used new observations to extend this study to even fainter flux densities, calculating the 15.7-GHz radio source count down to 0.1 mJy, a factor of five deeper than previous studies. There is no evidence for a new population of sources, showing that the high-frequency sky continues to be dominated by radio galaxies down to at least 0.1 mJy.
The frequency-dependent time delays in fast radio bursts (FRBs) can be used to constrain the photon mass, if the FRB redshifts are known, but the similarity between the frequency dependences of dispersion due to plasma effects and a photon mass complicates the derivation of a limit on $m_\gamma$. The redshift of FRB 150418 has been measured to $\sim 2$% and its dispersion measure (DM) is known to $\sim 0.1$%, but the strength of the constraint on $m_\gamma$ is limited by uncertainties in the modelling of the host galaxy and the Milky Way, as well as possible inhomogeneities in the intergalactic medium (IGM). Allowing for these uncertainties, the recent data on FRB 150418 indicate that $m_\gamma \lesssim 1.7 \times 10^{-14}$ eV c$^{-2}$ ($4.6 \times 10^{-50}$ kg). In the future, the different redshift dependences of the plasma and photon mass contributions to DM can be used to improve the sensitivity to $m_\gamma$ if more FRB redshifts are measured. For a fixed fractional uncertainty in the extra-galactic contribution to the DM of an FRB, one with a lower redshift would provide greater sensitivity to $m_\gamma$.
Links to: arXiv, form interface, find, astro-ph, recent, 1602, contact, help (Access key information)
Modified theories of gravity provide us with a unique opportunity to generate innovative tests of gravity. In Chameleon f(R) gravity, the gravitational potential differs from the weak-field limit of general relativity (GR) in a mass dependent way. We develop a probe of gravity which compares high mass clusters, where Chameleon effects are weak, to low mass clusters, where the effects can be strong. We utilize the escape velocity edges in the radius/velocity phase space to infer the gravitational potential profiles on scales of 0.3-1 virial radii. We show that the escape edges of low mass clusters are enhanced compared to GR, where the magnitude of the difference depends on the background field value |fR0|. We validate our probe using N-body simulations and simulated light cone galaxy data. For a DESI (Dark Energy Spectroscopic Instrument) Bright Galaxy Sample, including observational systematics, projection effects, and cosmic variance, our test can differentiate between GR and Chameleon f(R) gravity models, |fR0| = 4e-6 (2e-6) at > 5{\sigma} (> 2{\sigma}), more than an order of magnitude better than current cluster-scale constraints.
The mapping of dark matter clustering from real space to redshift space introduces the anisotropic property to the measured density power spectrum in redshift space, known as the Redshift Space Distortion effect. The mapping formula is intrinsically non-linear, which is complicated by the higher order polynomials due to indefinite cross correlations between the density and velocity fields, and the Finger--of--God effect due to the randomness of the peculiar velocity field. %Furthermore, the rigorous test of this mapping formula is contaminated by the unknown non--linearity of the density and velocity fields, including their auto- and cross-correlations, for calculating which our theoretical calculation breaks down beyond some scales. Whilst the full higher order polynomials remain unknown, the other systematics can be controlled consistently within the same order truncation in the expansion of the mapping formula, as shown in this paper. The systematic due to the unknown non--linear density and velocity fields is removed by separately measuring all terms in the expansion directly using simulations. The uncertainty caused by the velocity randomness is controlled by splitting the FoG term into two pieces, 1) the non--local FoG term being independent of the separation vector between two different points, and 2) the local FoG term appearing as an indefinite polynomials which is expanded in the same order as all other perturbative polynomials. Using 100 realizations of simulations, we find that the best fitted non--local FoG function is Gaussian, with only one scale--independent free parameter, and that our new mapping formulation accurately reproduces the observed 2-dimensional density power spectrum in redshift space at the smallest scales by far, up to $k\sim 0.3h$/Mpc$^{-1}$, considering the resolution of future experiments.
We investigate possible cosmological effects of interacting scalar radiation and dark matter. After its decoupling, scalar radiation can stream freely as neutrinos or self-interact strongly as perfect fluid, highly depending on the magnitude of its self-couplings. We obtain the general and novel structure for self-scattering rate and compare it with the expansion rate of our Universe. If its trilinear/cubic coupling is non-zero, scalar radiation can be eventually treated as perfect fluid. Possible effects on CMB are also discussed. When this scalar also mediates interaction among dark matter particles, the linear matter power spectrum for large scale structure can be modified differently from other models. We propose to use Debye shielding to avoid the singularity appearing in the scattering between scalar radiation and dark matter.
In this Letter we have derived the Jeans length in the context of the Kaniadakis statistics. We have compared this result with the Jeans length obtained in the non-extensive Tsallis statistics and discussed the main differences between these two models. We have also obtained the kappa-sound velocity. Finally, we have applied the results obtained here to analyze an astrophysical system.
In this work, we present a study of the central regions of cool-core clusters hosting radio mini-halos, which are di use synchrotron sources extended on cluster-scales surrounding the radio-loud brightest galaxy. We aim to investigate the interplay between the thermal and non-thermal components in the intracluster medium in order to get more insights into these radio sources, whose nature is still unclear. It has recently been proposed that turbulence plays a role for heating the gas in cool cores. A correlation between the radio luminosity of mini-halos, $\nu P_{\nu}$, and the cooling flow power, $P_{\rm CF}$, is expected in the case that this turbulence also plays a role for the acceleration of the relativistic particles. We carried out a homogeneous re-analysis of X-ray Chandra data of the largest sample of cool-core clusters hosting radio mini-halos currently available ($\sim20$ objects), finding a quasi-linear correlation, $\nu P_{\nu} \propto P_{\rm CF}^{0.8}$. We show that the scenario of a common origin of radio mini-halos and gas heating in cool-core clusters is energetically viable, provided that mini-halos trace regions where the magnetic field strength is $B \gg 0.5 \mu G$.
Non-thermal histories for the early universe have received notable attention as they are a rich source of phenomenology, while also being well motivated by top-down approaches to beyond the Standard Model physics. The early (pre-BBN) matter phase in these models leads to enhanced growth of density perturbations on sub-Hubble scales. Here we consider whether primordial black hole formation associated with the enhanced growth is in conflict with existing observations. Such constraints depend on the tilt of the primordial power spectrum, and we find that non-thermal histories are tightly constrained in the case of a significantly blue spectrum. Alternatively, if dark matter is taken to be of non-thermal origin we can restrict the primordial power spectrum on scales inaccessible to CMB and LSS observations. We establish constraints for a wide range of scalar masses (reheat temperatures) with the most stringent bounds resulting from the formation of $10^{15}$ g black holes. These black holes would be evaporating today and are constrained by FERMI observations. We also consider whether the breakdown of the coherence of the scalar oscillations on sub-horizon scales can lead to a Jean's pressure preventing black hole formation and relaxing our constraints. Our main conclusion is that primordial black hole constraints, combined with existing constraints on non-thermal WIMPs, favor a primordial spectrum closer to scale invariance or a red tilted spectrum.
New high-resolution r band imaging of the brightest cluster galaxy (BCG) in Abell 85 (Holm 15A) was obtained using the Gemini Multi Object Spectrograph. These data were taken with the aim of deriving an accurate surface brightness profile of the BCG of Abell 85, in particular its central region. The new Gemini data show clear evidence of a previously unreported nuclear emission that is evident as a distinct light excess in the central kiloparsec of the surface brightness profile. We find that the light profile is never flat nor does it present a downward trend towards the center of the galaxy. That is, the new Gemini data show a different physical reality from the featureless, "evacuated core" recently claimed for the Abell 85 BCG. After trying different models, we find that the surface brightness profile of the BCG of Abell 85 is best fit by a double Sersic model.
Over the past few years, several occasions of large, continuous rotations of the electric vector position angle (EVPA) of linearly polarized optical emission from blazars have been reported. These events are often coincident with high energy gamma-ray flares and they have attracted considerable attention, as they could allow one to probe the magnetic field structure in the gamma-ray emitting region of the jet. The flat-spectrum radio quasar 3C279 is one of the most prominent examples showing this behaviour. Our goal is to study the observed EVPA rotations and to distinguish between a stochastic and a deterministic origin of the polarization variability. We have combined multiple data sets of R-band photometry and optical polarimetry measurements of 3C279, yielding exceptionally well-sampled flux density and polarization curves that cover a period of 2008-2012. Several large EVPA rotations are identified in the data. We introduce a quantitative measure for the EVPA curve smoothness, which is then used to test a set of simple random walk polarization variability models against the data. 3C279 shows different polarization variation characteristics during an optical low-flux state and a flaring state. The polarization variation during the flaring state, especially the smooth approx. 360 degrees rotation of the EVPA in mid-2011, is not consistent with the tested stochastic processes. We conclude that during the two different optical flux states, two different processes govern the polarization variation, possibly a stochastic process during the low-brightness state and a deterministic process during the flaring activity.
Links to: arXiv, form interface, find, astro-ph, recent, 1603, contact, help (Access key information)
We consider the possibility that the black-hole (BH) binary detected by LIGO may be a signature of dark matter. Interestingly enough, there remains a window for masses $10\,M_\odot \lesssim M_{\rm bh} \lesssim 100\, M_\odot$ where primordial black holes (PBHs) may constitute the dark matter. If two BHs in a galactic halo pass sufficiently close, they can radiate enough energy in gravitational waves to become gravitationally bound. The bound BHs will then rapidly spiral inward due to emission of gravitational radiation and ultimately merge. Uncertainties in the rate for such events arise from our imprecise knowledge of the phase-space structure of galactic halos on the smallest scales. Still, reasonable estimates span a range that overlaps the $2-53$ Gpc$^{-3}$ yr$^{-1}$ rate estimated from GW150914, thus raising the possibility that LIGO has detected PBH dark matter. PBH mergers are likely to be distributed spatially more like dark matter than luminous matter and have no optical nor neutrino counterparts. They may be distinguished from mergers of BHs from more traditional astrophysical sources through the observed mass spectrum, their high ellipticities, or their stochastic gravitational wave background. Next generation experiments will be invaluable in performing these tests.
We introduce FastPM, a highly-scalable approximated particle mesh N-body solver, which implements the particle mesh (PM) scheme enforcing correct linear evolution via a cheap low resolution broadband correction at each time step. Employing a 2-dimensional domain decomposing scheme, FastPM scales extremely well with a very large number of CPUs. In contrast to COmoving-LAgrangian (COLA) approach, we do not require to split the force or track separately the 2LPT solution, reducing the code complexity and memory requirements. We compare FastPM with different number of steps ($N_s$) and force resolution factor ($B$) against 3 benchmarks: halo mass function from Friends of Friends halo finder, halo and dark matter power spectrum, and cross correlation coefficient (or stochasticity), relative to a high resolution TreePM simulation. We show that the broadband correction scheme reduces the halo stochasticity when compared to COLA with the same number of steps and force resolution. While increasing $N_s$ and $B$ improves the transfer function and cross correlation coefficient, for many applications FastPM achieves sufficient accuracy at low $N_s$ and $B$. For example, $N_s=10$ and $B=2$ simulation provides a factor of 10 saving of computing time relative to $N_s=40$, $B=3$ simulation, yet the halo benchmarks are very similar at $z=0$. For abundance matched halos the stochasticity remains low even for $N_s=5$. FastPM compares well against less expensive schemes, being only 7 (4) times more expensive than 2LPT initial condition generator for $N_s=10$ ($N_s=5$). Some of the applications where FastPM can be useful are generating a large number of mocks, producing non-linear statistics where one varies a large number of nuisance or cosmological parameters, or serving as part of an initial conditions solver.
21 cm Epoch of Reionization observations promise to transform our understanding of galaxy formation, but these observations are impossible without unprecedented levels of instrument calibration. We present end-to-end simulations of a full EoR power spectrum analysis including all of the major components of a real data processing pipeline: models of astrophysical foregrounds and EoR signal, frequency-dependent instrument effects, sky-based antenna calibration, and the full PS analysis. This study reveals that traditional sky-based per-frequency antenna calibration can only be implemented in EoR measurement analyses if the calibration model is unrealistically accurate. For reasonable levels of catalog completeness, the calibration introduces contamination in otherwise foreground-free power spectrum modes, precluding a PS measurement. We explore the origin of this contamination and potential mitigation techniques. We show that there is a strong joint constraint on the precision of the calibration catalog and the inherent spectral smoothness of antennae, and that this has significant implications for the instrumental design of the SKA and other future EoR observatories.
In this paper we investigate a situation where relativistic particles are reaccelerated diffusing across regions of reconnection and magnetic dynamo in super-Alfvenic, incompressible large-scale turbulence. We present an exploratory study of this mechanism in the intra-cluster-medium (ICM). In view of large-scale turbulence in the ICM we adopt a reconnection scheme that is based on turbulent reconnection and MHD turbulence. In this case particles are accelerated and decelerated in a systematic way in reconnecting and magnetic-dynamo regions, respectively, and on longer time-scales undergo a stochastic process diffusing across these sites (similar to second-order Fermi). Our study extends on larger scales numerical studies that focused on the acceleration in and around turbulent reconnecting regions. We suggest that this mechanism may play a role in the reacceleration of relativistic electrons in galaxy clusters providing a new physical scenario to explain the origin of cluster-scale diffuse radio emission. Indeed differently from current turbulent reacceleration models proposed for example for radio halos this mechanism is based on the effect of large-scale incompressible and super-Alfvenic turbulence. In this new model turbulence governs the interaction between relativistic particles and magnetic field lines that diffuse, reconnect and are stretched in the turbulent ICM.
Cluster star-forming galaxies are found to have an excess of Far-Infrared emission relative to H-alpha (Ha), when compared to those in the field, which could be caused by intense AGN activity, dust and/or declining star formation histories. Here we present spectroscopic observations of Ha emitters in the Cl 0939+4713 (Abell 851) super-cluster at z=0.41, using AF2+WYFFOS on the WHT. We measure [OII], Hbeta (Hb), [OIII], Ha and [NII] for a sample of 119 Ha emitters in and around the cluster. We find that 17+-5% of the Ha emitters are AGN, irrespective of environment. For star-forming galaxies, we obtain Balmer decrements, metallicities and ionisation parameters with different methods, individually and by stacking. We find a strong mass-metallicity relation at all environments, with no significant dependence on environment. The ionisation parameter declines with increasing stellar mass for low-mass galaxies. Ha emitters residing in intermediate environments show the highest ionisation parameters (along with high [OIII]/Ha and high [OIII]/[OII] line ratios, typically twice as large as in the highest and lowest densities), which decline with increasing environmental density. Dust extinction (A$_{H\alpha}$) correlates strongly with stellar mass, but also with environmental density. Star-forming galaxies in the densest environments are found to be significantly dustier (A$_{H\alpha}$~1.5-1.6) than those residing in the lowest density environments (A$_{H\alpha}$~0.6), deviating significantly from what would be predicted given their stellar masses.
We present a new technique for wide and shallow observations using the near-infrared channel of Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST). Wide-field near-IR surveys with HST are generally inefficient, as guide star acquisitions make it impractical to observe more than one pointing per orbit. This limitation can be circumvented by guiding with gyros alone, which is possible as long as the telescope has three functional gyros. The method presented here allows us to observe mosaics of eight independent WFC3-IR pointings in a single orbit by utilizing the fact that HST drifts by only a very small amount in the 25 seconds between non-destructive reads of unguided exposures. By shifting the reads and treating them as independent exposures the full resolution of WFC3 can be restored. We use this "drift and shift" (DASH) method in the Cycle 23 COSMOS-DASH program, which will obtain 456 WFC3 $H_{160}$ pointings in 57 orbits, covering an area of 0.6 degree$^2$ in the COSMOS field down to $H_{160} = 25$. When completed, the program will more than triple the area of extra-galactic survey fields covered by near-IR imaging at HST resolution. We demonstrate the viability of the method with the first four orbits (32 pointings) of this program. We show that the resolution of the WFC3 camera is preserved, and that structural parameters of galaxies are consistent with those measured in guided observations.
We present a strong and weak gravitational lens model of the galaxy cluster MACSJ0416.1-2403, constrained using spectroscopy from the Grism Lens-Amplified Survey from Space (GLASS) and Hubble Frontier Fields (HFF) imaging data. We search for emission lines in known multiply imaged sources in the GLASS spectra, obtaining secure spectroscopic redshifts of 31 multiple images belonging to 16 distinct source galaxies. The GLASS spectra provide the first spectroscopic measurements for 6 of the source galaxies. The weak lensing signal is acquired from 884 galaxies in the F606W HFF image. By combining the weak lensing constraints with 15 multiple image systems with spectroscopic redshifts and 9 multiple image systems with photometric redshifts, we reconstruct the gravitational potential of the cluster on an adaptive grid. The resulting total mass density map is compared with a stellar mass density map obtained from the deep Spitzer Frontier Fields imaging data to study the relative distribution of stellar and total mass in the cluster. We find that the projected stellar mass to total mass ratio, $f_{\star}$, varies considerably with the stellar surface mass density. The mean projected stellar mass to total mass ratio is $\langle f_{\star} \rangle= 0.009 \pm 0.003 $ (stat.), but with a systematic error as large as $0.004-0.005$, dominated by the choice of the IMF. We find agreement with several recent measurements of $f_{\star}$ in massive cluster environments. The lensing maps of convergence, shear, and magnification are made available to the broader community in the standard HFF format.
The cold dark matter model of structure formation faces apparent problems on
galactic scales. Several threads point to excessive halo concentration,
including central densities that rise too steeply with decreasing radius. Yet,
random fluctuations in the gaseous component can 'heat' the centres of haloes,
decreasing their densities. We present a theoretical model deriving this effect
from first principles: stochastic variations in the gas density are converted
into potential fluctuations that act on the dark matter; the associated force
correlation function is calculated and the corresponding stochastic equation
solved. Assuming a power law spectrum of fluctuations with maximal and minimal
cutoff scales, we derive the velocity dispersion imparted to the halo particles
and the relevant relaxation time. We further perform numerical simulations,
with fluctuations realised as a Gaussian random field, which confirm the
formation of a core within a timescale comparable to that derived analytically.
Non-radial collective modes enhance the energy transport process that erases
the cusp, though the parametrisations of the analytical model persist.
In our model, the dominant contribution to the dynamical coupling driving the
cusp-core transformation comes from the largest scale fluctuations. Yet, the
efficiency of the transformation is independent of the value of the largest
scale and depends weakly (linearly) on the power law exponent; it effectively
depends on two parameters: the gas mass fraction and the normalisation of the
power spectrum. This suggests that cusp-core transformations observed in
hydrodynamic simulations of galaxy formation may be understood and parametrised
in simple terms, the physical and numerical complexities of the various
implementations notwithstanding.
We present a renewed look at M31's Giant Stellar Stream along with the nearby structures Stream C and Stream D, exploiting a new algorithm capable of fitting to the red giant branch (RGB) of a structure in both colour and magnitude space. Using this algorithm, we are able to generate probability distributions in distance, metallicity and RGB width for a series of subfields spanning these structures. Specifically, we confirm a distance gradient of approximately 20 kpc per degree along a 6 degree extension of the Giant Stellar Stream, with the farthest subfields from M31 lying ~ 120 kpc more distant than the inner-most subfields. Further, we find a metallicity that steadily increases from -0.7^{+0.1}_{-0.1} dex to -0.2^{+0.2}_{-0.1} dex along the inner half of the stream before steadily dropping to a value of -1.0^{+0.2}_{-0.2} dex at the farthest reaches of our coverage. The RGB width is found to increase rapidly from 0.4^{+0.1}_{-0.1} dex to 1.1^{+0.2}_{-0.1} dex in the inner portion of the stream before plateauing and decreasing marginally in the outer subfields of the stream. In addition, we estimate Stream C to lie at a distance between 794 and 862 kpc and Stream D between 758 kpc and 868 kpc. We estimate the median metallicity of Stream C to lie in the range -0.7 to -1.6 dex and a metallicity of -1.1^{+0.3}_{-0.2} dex for Stream D. RGB widths for the two structures are estimated to lie in the range 0.4 to 1.2 dex and 0.3 to 0.7 dex respectively. In total, measurements are obtained for 19 subfields along the Giant Stellar Stream, 4 along Stream C, 5 along Stream D and 3 general M31 spheroid fields for comparison. We thus provide a higher resolution coverage of the structures in these parameters than has previously been available in the literature.
Within the context of Newton's theory of gravitation, restricted to point-like test particles and central bodies, stable circular orbits in ordinary space are related to stable circular paths on a massless, unmovable, undeformable vortex-like surface, under the action of a tidal gravitational field along the symmetry axis. An interpretation is made in the light of a holographic principle, in the sense that motions in ordinary space are connected with motions on a selected surface and vice versa. Then ordinary space is conceived as a 3-hypersurface bounding a $n$-hypervolume where gravitation takes origin, within a $n$-hyperspace. The extension of the holographic principle to extra dimensions implies the existence of a minimum distance where test particles may still be considered as distinct from the central body. Below that threshold, it is inferred test particles lose theirs individuality and "glue" to the central body via unification of the four known interactions and, in addition, (i) particles can no longer be conceived as point-like but e.g., strings or membranes, and (ii) quantum effects are dominant and matter turns back to a pre-big bang state. A more detailed formulation including noncircular motions within the context of general relativity, together with further knowledge on neutron stars, quark stars and black holes, would provide further insight on the formulation of quantum gravity.
We reveal the universality of short-term anisotropic inflation. As a demonstration, we study inflation with an exponential type gauge kinetic function which is ubiquitous in models obtained by dimensional reduction from higher dimensional fundamental theory. It turns out that an anisotropic inflation universally takes place in the later stage of conventional inflation. Remarkably, we find that primordial gravitational waves with a peak amplitude around $10^{-26}$ ~ $10^{-27}$ are copiously produced in high-frequency bands 10MHz~100MHz. If we could detect such gravitational waves in future, we would be able to probe higher dimensional fundamental theory.
In this paper, we constrain four time-dependent dark energy (TDDE) models by using the Type Ia supernovae (SNe Ia), baryonic acoustic oscillations (BAO), observational Hubble parameter (OHD) data-sets as well as the single data point from the newest event GW150914. Subsequently, adopting the best fitting values of the model parameters, we apply the original statefinder, statefinder hierarchy, the growth rate of matter perturbations and $Om(z)$ diagnostics to distinguish the TDDE scenarios and the $\Lambda$CDM scenario from each other. We discover that all the TDDE models and $\Lambda$CDM model can be distinguished better at the present epoch by using the statefinder hierarchy than using the original statefinder, the growth rate of matter perturbations and $Om(z)$ diagnostics, especially, in the planes of $\{S_3^{(1)},S_4^{(1)}\}$, $\{S_3^{(2)},S_4^{(2)}\}$, $\{S_5^{(1)},S_5^{(2)}\}$ and $\{S_4^{(2)},S_5^{(2)}\}$.
The population of supermassive black holes (SMBHs) is composed by quiescent SMBHs, such as those seen in local galaxies including the Milky Way's, and active ones, resulting in quasars and active galactic nuclei (AGN). Outside our neighbourhood, all the information we have on SMBHs is derived from quasars and AGN, giving us a partial view. We study the evolution of the SMBH population, total and active, by the continuity equation, backwards in time from z=0 to z=4. Type-1 and type-2 AGN are differentiated in the model on the basis of the Eddington ratio distribution, chosen on the basis of observational estimates. The duty cycle is obtained by matching the luminosity function of quasars, and the average radiative efficiency is the only free parameter in the model. For higher radiative efficiencies (>~0.07) a large fraction of the SMBH population, most of them quiescent, must already be in place by z=4. For lower radiative efficiencies (~0.05), the duty cycle increases with the redshift and the SMBH population evolves dramatically since z=4. The mass function of active SMBHs does not depend on the choice of the radiative efficiency or of the local SMBH mass function, but it is mainly determined by the quasar luminosity function, once the Eddington ratio distribution is fixed. Only a direct measurement of the total BHMF at redshifts z>~2 could break these degeneracies giving important constraints on the average radiative efficiency. Focusing on type-1 AGN, for which observational estimates of the mass function and Eddington ratio distribution exist at various redshift, models with lower radiative efficiencies reproduce better the high-mass end of the mass function at high-z, but they tend to over-predict it at low-z, and vice-versa for models with higher radiative efficiencies.
We discuss the coupling of the electromagnetic field with a curved and torsioned Lyra manifold using the Duffin-Kemmer-Petiau theory. We will show how to obtain the equations of motion and energy-momentum and spin density tensors by means of the Schwinger Variational Principle.
Links to: arXiv, form interface, find, astro-ph, recent, 1603, contact, help (Access key information)
We present a detailed modelling of the joint covariance matrix between cluster number counts and the galaxy angular power spectrum. To this end, we use a Halo Model framework complemented by a Halo Occupation Distribution model (HOD), and we work in full-sky. We demonstrate the importance of accounting for non-Gaussianity to produce accurate covariance predictions, as the Gaussian part of the covariance can in fact become subdominant in certain configurations. We discuss in particular the case of the super-sample covariance (SSC), including the effects of galaxy shot-noise, halo second order bias and non-local bias, and demonstrating interesting mathematical properties. Using the joint covariance matrix and a Fisher matrix methodology, we examine the prospects of combining these two probes to constrain cosmological and HOD parameters. We find that the combination indeed results in noticeable better constraints, in particular because the cross-covariance introduces a synergy between the probes on small scales. We conclude that accounting for all non-Gaussian effects is important to best extract information from galaxy surveys with these two observables.
Most inflationary models predict primordial perturbations to be statistically isotropic and homogeneous. Cosmic-Microwave-Background (CMB) observations, however, indicate a possible departure from statistical isotropy in the form of a dipolar power modulation at large angular scales. Alternative models of inflation, beyond the simplest single-field slow-roll models, can generate a small power asymmetry, consistent with these observations. Observations of clustering of quasars show, however, agreement with statistical isotropy at much smaller angular scales. Here we propose to use off-diagonal components of the angular power spectrum of the 21-cm fluctuations during the dark ages to test this power asymmetry. We forecast results for the planned SKA radio array, a future radio array, and the cosmic-variance-limited case as a theoretical proof of principle. Our results show that the 21-cm-line power spectrum will enable access to information at very small scales and at different redshift slices, thus improving upon the current CMB constraints by $\sim 2$ orders of magnitude for a dipolar asymmetry, and by $\sim 1-3$ orders of magnitude for a quadrupolar asymmetry case.
Automated photometric supernova classification has become an active area of research in recent years in light of current and upcoming imaging surveys such as the Dark Energy Survey (DES) and the Large Synoptic Telescope (LSST), given that spectroscopic confirmation of type for all supernovae discovered with these surveys will be impossible. Here, we develop a multi-faceted classification pipeline, combining existing and new approaches. Our pipeline consists of two stages: extracting descriptive features from the light curves and classification using a machine learning algorithm. Our feature extraction methods vary from model-dependent techniques, namely SALT2 fits, to more independent techniques fitting parametric models to curves, to a completely model-independent wavelet approach. We cover a range of representative machine learning algorithms, including naive Bayes, k-nearest neighbors, support vector machines, artificial neural networks and boosted decision trees. We test the pipeline on simulated multi-band DES light curves from the Supernova Photometric Classification Challenge. Using the commonly-used area under the Receiver Operating Characteristic curve (AUC) as a metric, we find that the SALT2 fits and the wavelet approach, with the boosted decision trees algorithm, each achieves an AUC of 0.98, where 1 represents perfect classification. We find that a representative training set is essential for good classification, whatever the feature set or algorithm, suggesting that spectroscopic follow-up is best focused on fewer objects at a broad range of redshifts than on many bright objects. Importantly, we find that by using either one of the two best feature extraction methods (SALT2 model fits and wavelet decomposition) and a boosted decision tree algorithm, accurate classification is possible purely from light curve data, without the need for any redshift information. [Abridged]
We present measurements of the [NII]/Ha ratio as a probe of gas-phase oxygen abundance for a sample of 419 star-forming galaxies at z=0.6-2.7 from the KMOS3D near-IR multi-IFU survey. The mass-metallicity relation (MZR) is determined consistently with the same sample selection, metallicity tracer, and methodology over the wide redshift range probed by the survey. We find good agreement with long-slit surveys in the literature, except for the low-mass slope of the relation at z~2.3, where this sample is less biased than previous samples based on optical spectroscopic redshifts. In this regime we measure a steeper slope than some literature results. Excluding the AGN contribution from the MZR reduces sensitivity at the high mass end, but produces otherwise consistent results. There is no significant dependence of the [NII]/Ha ratio on SFR or environment at fixed redshift and stellar mass. The IFU data allow spatially resolved measurements of [NII]/Ha, from which we can infer abundance gradients for 180 galaxies, thus tripling the current sample in the literature. The observed gradients are on average flat, with only 15 gradients statistically offset from zero at >3sigma. We have modelled the effect of beam-smearing, assuming a smooth intrinsic radial gradient and known seeing, inclination and effective radius for each galaxy. Our seeing-limited observations can recover up to 70% of the intrinsic gradient for the largest, face-on disks, but only 30% for the smaller, more inclined galaxies. We do not find significant trends between observed or corrected gradients and any stellar population, dynamical or structural galaxy parameters, mostly in agreement with existing studies with much smaller sample sizes. In cosmological simulations, strong feedback is generally required to produce flat gradients at high redshift.
Links to: arXiv, form interface, find, astro-ph, recent, 1603, contact, help (Access key information)