Links to: arXiv, form interface, find, astro-ph, recent, 1710, contact, help (Access key information)
It has recently been proposed that massive primordial black holes (PBH) could constitute all of the dark matter, providing a novel scenario of structure formation, with early reionization and a rapid growth of the massive black holes at the center of galaxies and dark matter halos. The scenario arises from broad peaks in the primordial power spectrum that give both a spatially clustered and an extended mass distribution of PBH. The constraints from the observed microlensing events on the extended mass function have already been addressed. Here we study the impact of spatial clustering on the microlensing constraints. We find that the bounds can be relaxed significantly for relatively broad mass distributions if the number of primordial black holes within each cluster is typically above one hundred. On the other hand, even if they arise from individual black holes within the cluster, the bounds from CMB anisotropies are less stringent due to the enhanced black hole velocity in such dense clusters. This way, the window between a few and ten solar masses has opened up for PBH to comprise the totality of the dark matter.
We perform a quantitative comparison between N-body simulations and the Schr\"odinger-Poisson system in 1+1 dimensions. In particular, we study halo formation with different initial conditions. We observe the convergence of various observables in the Planck constant h-bar and also test virialization. We discuss the generation of higher order cumulants of the particle distribution function which demonstrates that the Schr\"odinger-Poisson equations should not be perceived as a generalization of the dust model with quantum pressure but rather as one way of sampling the phase space of the Vlasov-Poisson system -- just as N-body simulations. Finally, we quantitatively recover the scaling behavior of the halo density profile from N-body simulations.
The main purpose of this work is to distinguish various holographic type dark energy (DE) models, including the $\Lambda$HDE, HDE, NADE and RDE model, by using various diagnostic tools. The first diagnostic tool is the Statefinder hierarchy, in which the evolution of Statefinder hierarchy parmeter $S^{(1)}_3(z)$ and $S^{(1)}_4(z)$ are studied. The second is composite null diagnostic (CND), in which the trajectories of $\{S^{(1)}_3, \epsilon\}$ and $\{S^{(1)}_4, \epsilon\}$ are investigated, where $\epsilon$ is the fractional growth parameter. The last is $w-w'$ analysis, where $w$ is the equation of state for DE and the prime denotes derivative with respect to $ln a$. In the analysis we consider two cases: varying current fractional DE density $\Omega_{de0}$ and varying DE model parameter $C$. We find that: (1) Both the Statefinder hierarchy and the CND have qualitative impact on $\Lambda$HDE, but only have quantitative impact on HDE. (2) $S_4^{(1)}$ can lead to larger differences than $S_3^{(1)}$, while the CND pair has a stronger ability to distinguish different models than the Statefinder hierarchy. (3) For the case of varying $C$, the $\{w, w'\}$ pair has qualitative impact on $\Lambda$HDE; for the case of varying $\Omega_{de0}$, the $\{w, w'\}$ pair only has quantitative impact; these results are different from the cases of HDE, RDE and NADE, in which the $\{w, w'\}$ pair only has quantitative impact on these models. In conclusion, compared with HDE, RDE and NADE, the $\Lambda$HDE model can be easily distinguished by using these diagnostic tools.
It has been previously shown that the linear inflation appears naturally as a solution of Coleman-Weinberg inflation, provided that the inflaton has a non-minimal coupling to gravity and the Planck scale is dynamically generated. We revisit the previous study by improving the discussion of reheating and by comparing the results of the metric and the Palatini formulations of non-minimal gravity. We find that both formulations predict linear inflation but a different number of $e$-folds. If the non-minimal coupling is larger than one, future experimental sensitivity can discriminate between the two realizations.
We present a blind search for doublet intergalactic metal absorption with a method we dub `agnostic stacking'. Using forward-modelling we combine it with direct detections in the literature to measure the overall metal population. Here we apply this novel approach to the search for NeVIII in 26 high-quality COS spectra of QSOs at z>0.7. We probe an unprecedented low limit of log N>12.3 at 0.47<z<1.34 with a total pathlength ${\Delta}$z = 7.36. The method selects absorption without requiring knowledge of its source, be it observing noise, artifacts, or any line transition. Stacking this mixed population with NeVIII absorption dilutes doublet features in composite spectra in a deterministic manner. We stack potential NeVIII absorption in two regimes: absorption too weak to be statistically significant in direct line studies (12.3 < log N< 13.7), and strong absorbers (log N> 13.7). We do not detect NeVIII in either regime, and place upper limits on the population using agnostic stacking alone. Combining our measurements with direct line detections, the NeVIII population is reproduced with a single power law column density distribution of slope \b{eta} = -1.86 and normalisation log f_{13.7} = -13.99, leading to an incidence rate of strong NeVIII absorbers of dn/dz =1.38. Comparing our results with a group of 3 systems in PG1148+549, these have a 0.024% probability of arising by chance. We infer a cosmic mass density for NeVIII in the column density range 12.3 < log N < 15.0 of ${\Omega}$(NeVIII) = 2.2x10^{-8}. We translate this inferred density into an estimate of the baryon density of the NeVIII-bearing gas, and arrive at ${\Omega}b~1.8x10^{-3}$, which constitutes only 4% of the total baryonic mass. The measured NeVIII column density distribution function and cosmic density here are inconsistent with predictions of the EAGLES simulations at ${\sigma}>2.0$ significance. (abridged)
Recently there has been much interest in the spatial distribution of light scalar dark matter, especially axions, throughout the universe. When the local gravitational interactions between the scalar modes are sufficiently rapid, it can cause the field to re-organize into a BEC of gravitationally bound clumps. While these clumps are stable when only gravitation is included, the picture is complicated by the presence of the axion's attractive self-interactions, which can potentially cause the clumps to collapse. Here we perform a detailed stability analysis to determine under what conditions the clumps are stable. In this paper we focus on spherical configurations, leaving aspherical configurations for future work. We identify branches of clump solutions of the axion-gravity-self-interacting system and study their stability properties. We find that clumps that are (spatially) large are stable, while clumps that are (spatially) small are unstable and may collapse. Furthermore, there is a maximum number of particles that can be in a clump. We map out the full space of solutions, which includes quasi-stable axitons, and clarify how a recent claim in the literature of a new ultra-dense branch of stable solutions rests on an invalid use of the non-relativistic approximation. We also consider repulsive self-interactions that may arise from a generic scalar dark matter candidate, finding a single stable branch that extends to arbitrary particle number.
Though theoretically expected, the charge exchange emission from galaxy clusters has not yet been confidently detected. Accumulating hints were reported recently, including a rather marginal detection with the Hitomi data of the Perseus cluster. As suggested in Gu et al. (2015), a detection of charge exchange line emission from galaxy clusters would not only impact the interpretation of the newly-discovered 3.5 keV line, but also open up a new research topic on the interaction between hot and cold matter in clusters. We aim to perform the most systematic search for the O VIII charge exchange line in cluster spectra using the RGS on board XMM. We introduce a sample of 21 clusters observed with the RGS. The dominating thermal plasma emission is modeled and subtracted with a two-temperature CIE component, and the residuals are stacked for the line search. The systematic uncertainties in the fits are quantified by refitting the spectra with a varying continuum and line broadening. By the residual stacking, we do find a hint of a line-like feature at 14.82 A, the characteristic wavelength expected for oxygen charge exchange. This feature has a marginal significance of 2.8 sigma, and the average equivalent width is 2.5E-4 keV. We further demonstrate that the putative feature can be hardly affected by the systematic errors from continuum modelling and instrumental effects, or the atomic uncertainties of the neighbouring thermal lines. Assuming a realistic temperature and abundance pattern, the physical model implied by the possible oxygen line agrees well with the theoretical model proposed previously to explain the reported 3.5 keV line. If the charge exchange source indeed exists, we would expect that the oxygen abundance is potentially overestimated by 8-22% in previous X-ray measurements which assumed pure thermal lines.
The direct detection of gravitational waves (GW) from merging binary black holes marks the beginning of a new era in gravitational physics, and it brings forth new opportunities to test theories of gravity. To this end, it is crucial to search for anomalous deviations from general relativity in a model-independent way, irrespective of gravity theories, GW sources, and background spacetimes. In this paper, we propose a new universal framework for testing gravity with GW, based on the generalized propagation of a GW in an effective field theory that describes modification of gravity at cosmological scales. Then we perform a parameter estimation study, showing how well the future observation of GW can constrain the model parameters in the generalized models of GW propagation.
Links to: arXiv, form interface, find, astro-ph, recent, 1710, contact, help (Access key information)
We investigate how Einstein rings and magnified arcs are affected by small-mass dark-matter haloes placed along the line-of-sight to gravitational lens systems. By comparing the gravitational signature of line-of-sight haloes with that of substructures within the lensing galaxy, we derive a mass-redshift relation that allows us to rescale the detection threshold (i.e. lowest detectable mass) for substructures to a detection threshold for line-of-sight haloes at any redshift. We then quantify the line-of-sight contribution to the total number density of low-mass objects that can be detected through strong gravitational lensing. Finally, we assess the degeneracy between substructures and line-of-sight haloes of different mass and redshift to provide a statistical interpretation of current and future detections, with the aim of distinguishing between CDM and WDM. We find that line-of-sight haloes statistically dominate with respect to substructures, by an amount that strongly depends on the source and lens redshifts, and on the chosen dark matter model. Substructures represent about 30 percent of the total number of perturbers for low lens and source redshifts (as for the SLACS lenses), but less than 10 per cent for high redshift systems. We also find that for data with high enough signal-to-noise ratio and angular resolution, the non-linear effects arising from a double-lens-plane configuration are such that one is able to observationally recover the line-of-sight halo redshift with an absolute error precision of 0.15 at the 68 per cent confidence level.
The ground-breaking detections of gravitational waves from black hole mergers by LIGO have rekindled interest in primordial black holes (PBHs) and the possibility of dark matter being composed of PBHs. It has been suggested that PBHs of tens of solar masses could serve as dark matter candidates. Recent analytical studies demonstrated that compact ultra-faint dwarf galaxies can serve as a sensitive test for the PBH dark matter hypothesis, since stars in such a halo-dominated system would be heated by the more massive PBHs, their present-day distribution can provide strong constraints on PBH mass. In this study, we further explore this scenario with more detailed calculations, using a combination of dynamical simulations and Bayesian inference methods. The joint evolution of stars and PBH dark matter is followed with a Fokker-Planck code PhaseFlow. We run a large suite of such simulations for different dark matter parameters, then use a Markov Chain Monte Carlo approach to constrain the PBH properties with observations of ultra-faint galaxies. We find that two-body relaxation between the stars and PBH drives up the stellar core size, and increases the central stellar velocity dispersion. Using the observed half-light radius and velocity dispersion of stars in the compact ultra-faint dwarf galaxies as joint constraints, we infer that these dwarfs may have a cored dark matter halo with the central density in the range of 1-2 $\rm{M_{\odot}/pc^3}$, and that the PBHs may have a mass range of 2-14 $\rm{M_{\odot}}$ if they constitute all or a substantial fraction of the dark matter.
We derive cosmological constraints from the probability distribution function
(PDF) of evolved large-scale matter density fluctuations. We do this by
splitting lines of sight by density based on their count of tracer galaxies,
and by measuring both gravitational shear around and counts-in-cells in
overdense and underdense lines of sight, in Dark Energy Survey (DES) First Year
and Sloan Digital Sky Survey (SDSS) data. Our analysis uses a perturbation
theory model (see companion paper Friedrich at al.) and is validated using
N-body simulation realizations and log-normal mocks. It allows us to constrain
cosmology, bias and stochasticity of galaxies w.r.t. matter density and, in
addition, the skewness of the matter density field.
From a Bayesian model comparison, we find that the data weakly prefer a
connection of galaxies and matter that is stochastic beyond Poisson
fluctuations on <=20 arcmin angular smoothing scale. The two stochasticity
models we fit yield DES constraints on the matter density
$\Omega_m=0.26^{+0.04}_{-0.03}$ and $\Omega_m=0.28^{+0.05}_{-0.04}$ that are
consistent with each other. These values also agree with the DES analysis of
galaxy and shear two-point functions (3x2pt) that only uses second moments of
the PDF. Constraints on $\sigma_8$ are model dependent
($\sigma_8=0.97^{+0.07}_{-0.06}$ and $0.80^{+0.06}_{-0.07}$ for the two
stochasticity models), but consistent with each other and with the 3x2pt
results if stochasticity is at the low end of the posterior range.
As an additional test of gravity, counts and lensing in cells allow to
compare the skewness $S_3$ of the matter density PDF to its LCDM prediction. We
find no evidence of excess skewness in any model or data set, with better than
25 per cent relative precision in the skewness estimate from DES alone.
We present density split statistics, a framework that studies lensing and
counts-in-cells as a function of foreground galaxy density, thereby providing a
large-scale measurement of both 2-point and 3-point statistics. Our method
extends our earlier work on trough lensing and is summarized as follows: given
a foreground (low redshift) population of galaxies, we divide the sky into
subareas of equal size but distinct galaxy density. We then measure lensing
around uniformly spaced points separately in each of these subareas, as well as
counts-in-cells statistics (CiC). The lensing signals trace the matter density
contrast around regions of fixed galaxy density. Through the CiC measurements
this can be related to the density profile around regions of fixed matter
density. Together, these measurements constitute a powerful probe of cosmology,
the skewness of the density field and the connection of galaxies and matter.
In this paper we show how to model both the density split lensing signal and
CiC from basic ingredients: a non-linear power spectrum, clustering hierarchy
coefficients from perturbation theory and a parametric model for galaxy bias
and shot-noise. Using N-body simulations, we demonstrate that this model is
sufficiently accurate for a cosmological analysis on year 1 data from the Dark
Energy Survey.
We use the Fisher information matrix to investigate the angular resolution and luminosity distance uncertainty for coalescing binary neutron stars (BNSs) and neutron star-black hole binaries (NSBHs) detected by the third-generation (3G) gravitational-wave (GW) detectors. Our study focuses on an individual 3G detector and a network of up to four 3G detectors at different locations including the US, Europe, China and Australia for the proposed Einstein Telescope (ET) and Cosmic Explorer (CE) detectors and an ideal detector with a flat low-frequency sensitivity. We find that, due to the effect of the Earth's rotation, a time-dependent antenna beam-pattern function can help better localize BNS and NSBH sources, especially those edge-on ones. We use numerical simulations to study the localization for a random sample of (1.4+1.4) ${\rm M}_\odot$ BNSs and low-mass NSBHs of (1.4+10) ${\rm M}_\odot$ at various redshifts. The medium angular resolution for a network of two CE detectors in the US and Europe respectively is around 20 deg$^2$ at $z=0.2$ for our BNS and NSBH samples. A 20 deg$^2$ medium angular resolution can be achieved for a network of two ET-D detectors at a much higher redshift of $z=0.5$ than for two CEs. We discuss the implications of our results to multi-messenger astronomy and in particular to using GW sources as independent tools to constrain the cosmological parameters. We find that in general, if 10 BNSs or NSBHs at $z=0.1$ with known redshifts are detected with $\le 50\%$ distance uncertainty by 3G networks consisting of two ET-like detectors, the Hubble constant $H_0$ can be measured with an accuracy of $0.9\%$. If 1000 face-on BNSs at $z<2$ are detected with known redshifts, we are able to constrain the equation-of-state parameters of dark energy $w_0$ and $w_a$ with accuracies $\Delta w_0=0.03$ and $\Delta w_a=0.2$, respectively.(Abridged version).
Motivated by the recent detection of the gravitational wave signal emitted by a binary neutron star merger, we analyse the possible impact of dark matter on such signals. We show that dark matter cores in merging neutron stars may yield an observable supplementary peak in the gravitational wave power spectral density following the merger, which could be distinguished from the features produced by the neutron components.
In this paper we present cosmological constraints on several well known $f(R)$ models, but also on a new class of models that are variants of the Hu-Sawicki one of the form $f(R)=R-\frac{2\Lambda}{1+b\;y(R,\Lambda)}$, that interpolate between the cosmological constant model and a matter dominated universe for different values of the parameter $b$, usually expected to be small for viable models and which in practice measures the deviation from General Relativity. We use the latest growth rate, Cosmic Microwave Background, Baryon Acoustic Oscillations, Supernovae type Ia and Hubble parameter data to place stringent constraints on the models and to compare them to the cosmological constant model but also other viable $f(R)$ models such as the Starobinsky or the degenerate hypergeometric models. We find that this kind of Hu-Sawicki variant parameterizations are in general compatible with the currently available data and can provide useful toy models to explore the available functional space of $f(R)$ models, something very useful with the current and upcoming surveys that will test deviations from General Relativity.
The evolution of the linear and scale independent bias, based on the most popular dark matter bias models within the $\Lambda$CDM cosmology, is confronted to that of the Dark Energy Survey (DES) Luminous Red Galaxies (LRGs). Applying a $\chi^2$ minimization procedure between models and data we find that all the considered linear bias models reproduce well the LRG bias data. The differences among the bias models are absorbed in the predicted mass of the dark-matter halo in which LRGs live and which ranges between $\sim 8 \times 10^{12} h^{-1} M_{\odot}$ and $2.3 \times 10^{13} h^{-1} M_{\odot}$, for the different models. Similar results,as far as the DM halo mass are concerned, are found by confronting the theoretical angular clustering models, which also include the evolution of bias, and the corresponding 2SLAQ LRG clustering. This analysis provided also a value of $\Omega_{m}=0.30\pm 0.01$, which is in excellent agreement with recent joint analyses of different cosmological probes and the reanalysis of the Planck data.
Using Hamilton-Jacobi formalism, The scenario of warm inflation with viscous pressure is considered. The formalism gives a way of computing the slow-rolling parameters without extra approximation, and it is well-known as a powerful method in cold inflation. The model is studied in detail for three different cases of dissipation and bulk viscous pressure coefficients. In the first case where both coefficients are taken as a constant, it is shown that the case could not portray warm inflationary scenario compatible with observational data even it is possible to restrict the model parameters. For other cases, the results shows that the model could properly predicts the perturbation parameters in which they stay in perfect agreement with Planck data. As a further argument, $r-n_s$ and $\alpha_s-n_s$ are drown that show the required result could stand in acceptable area expressing a compatibility with observational data.
The article is dedicated to one of the most undeservedly overlooked properties of the cosmological models: the behaviour at, near and due to a jump discontinuity. It is most interesting that while the usual considerations of the cosmological dynamics deals heavily in the singularities produced by the discontinuities of the second kind (a.k.a. the essential discontinuities) of one (or more) of the physical parameters, almost no research exists to date that would turn to their natural extension/counterpart: the singularities induced by the discontinuities of the first kind (a.k.a. the jump discontinuities). It is this oversight that this article aims to amend. In fact, it demonstrates that the inclusion of such singularities allows one to produce a number of very interesting scenarios of cosmological evolution. For example, it produces the cosmological models with a finite value of the equation of state parameter $w=p/\rho$ even when both the energy density and the pressure diverge, while at the same time keeping the scale factor finite. Such a dynamics is shown to be possible only when the scale factor experiences a finite jump at some moment of time. Furthermore, if it is the first derivative of the scale factor that experiences a jump, then a whole new and different type of a sudden future singularity appears. Finally, jump discontinuities suffered by either a second or third derivatives of a scale factor lead to cosmological models experiencing a sudden dephantomization -- or avoiding the phantomization altogether. This implies that theoretically there should not be any obstacles for extending the cosmological evolution beyond the corresponding singularities; therefore, such singularities can be considered a sort of a cosmological phase transition.
We discuss the feasibility of using astrophysical observations of white dwarfs as probes of fundamental physics. We quantify the effects of varying fundamental couplings on the white dwarf mass-radius relation in a broad class of unification scenarios, both for the simple case of a polytropic stellar structure model and for more general models. Independent measurements of the mass and radius, together with direct spectroscopic measurements of the fine-structure constant in white dwarf atmospheres lead to constraints on combinations of the two phenomenological parameters describing the underlying unification scenario (one of which is related to the strong sector of the theory while the other is related to the electroweak sector). While currently available measurements do not yet provide stringent constraints, we show that forthcoming improvements, expected for example from the Gaia satellite, can break parameter degeneracies and lead to constraints that ideally complement those obtained from local laboratory tests using atomic clocks.
The detection of GW170817 in both gravitational waves and electromagnetic waves heralds the age of gravitational-wave multi-messenger astronomy. On 17 August 2017 the Advanced LIGO and Virgo detectors observed GW170817, a strong signal from the merger of a binary neutron-star system. Less than 2 seconds after the merger, a gamma-ray burst (GRB 170817A) was detected within a region of the sky consistent with the LIGO-Virgo-derived location of the gravitational-wave source. This sky region was subsequently observed by optical astronomy facilities, resulting in the identification of an optical transient signal within $\sim 10$ arcsec of the galaxy NGC 4993. These multi-messenger observations allow us to use GW170817 as a standard siren, the gravitational-wave analog of an astronomical standard candle, to measure the Hubble constant. This quantity, which represents the local expansion rate of the Universe, sets the overall scale of the Universe and is of fundamental importance to cosmology. Our measurement combines the distance to the source inferred purely from the gravitational-wave signal with the recession velocity inferred from measurements of the redshift using electromagnetic data. This approach does not require any form of cosmic "distance ladder;" the gravitational wave analysis can be used to estimate the luminosity distance out to cosmological scales directly, without the use of intermediate astronomical distance measurements. We determine the Hubble constant to be $70.0^{+12.0}_{-8.0} \, \mathrm{km} \, \mathrm{s}^{-1} \, \mathrm{Mpc}^{-1}$ (maximum a posteriori and 68% credible interval). This is consistent with existing measurements, while being completely independent of them. Additional standard-siren measurements from future gravitational-wave sources will provide precision constraints of this important cosmological parameter.
The observation of GW170817 and its electromagnatic counterpart implies that gravitational waves travel at the speed of light, with deviations smaller than a few parts in $10^{-15}$. We discuss the consequences of this experimental result for models of dark energy and modified gravity characterized by a single scalar degree of freedom. To avoid tuning, the speed of gravitational waves must be unaffected not only for our particular cosmological solution, but also for nearby solutions obtained by slightly changing the matter abundance. For this to happen the coefficients of various operators must satisfy precise relations that we discuss both in the language of the Effective Field Theory of Dark Energy and in the covariant one, for Horndeski, beyond Horndeski and degenerate higher-order theories. The simplification is dramatic: of the three functions describing quartic and quintic beyond Horndeski theories, only one remains and reduces to a standard conformal coupling to the Ricci scalar for Horndeski theories. We show that the deduced relations among operators do not introduce further tuning of the models, since they are stable under quantum corrections.
The LIGO/VIRGO collaboration has recently announced the detection of gravitational waves from a neutron star-neutron star merger (GW170817) and the simultaneous measurement of an optical counterpart (the gamma-ray burst GRB 170817A). The close arrival time of the gravitational and electromagnetic waves limits the difference in speed of photons and gravitons to be less than about one part in $10^{15}$. This has three important implications for cosmological scalar-tensor gravity theories that are often touted as dark energy candidates and alternatives to $\Lambda$CDM. First, for the most general scalar-tensor theories---beyond Horndeski models---three of the five parameters appearing in the effective theory of dark energy can now be severely constrained on astrophysical scales; we present the results of combining the new gravity wave results with galaxy cluster observations. Second, the combination with the lack of strong equivalence principle violations exhibited by the supermassive black hole in M87, constrains the quartic galileon model to be cosmologically irrelevant. Finally, we derive a new bound on the disformal coupling to photons that implies that such couplings are irrelevant for the cosmic evolution of the field.
Multi-messenger gravitational wave (GW) astronomy has commenced with the detection of the binary neutron star merger GW170817 and its associated electromagnetic counterparts. The almost coincident observation of the GW and the gamma ray burst GRB170817A constrain the speed of GWs at the level of $|c_g/c-1|\leq4.5\cdot10^{-16}$. We use this result to probe the nature of dark energy (DE), showing that scalar-tensor theories with derivative interactions with the curvature are highly disfavored. As an example we consider the case of Galileons, a well motivated gravity theory with viable cosmology, which predicts a variable GW speed at low redshift, and is hence strongly ruled out by GW170817. Our result essentially eliminates any cosmological application of these DE models and, in general, of quartic and quintic Horndeski and most beyond Horndeski theories. We identify the surviving scalar-tensor models and, in particular, present specific beyond Horndeski theories avoiding this constraint. The viable scenarios are either conformally equivalent to theories in which $c_g=c$ or rely on cancellations of the anomalous GW speed that are valid on arbitrary backgrounds. Our conclusions can be extended to any other gravity theory predicting an anomalous GW propagation speed such as Einstein-Aether, Ho\v{r}ava gravity, Generalized Proca, TeVeS and other MOND-like gravities.
Using $\sim 100$ X-ray selected clusters in the Dark Energy Survey Science Verification data, we constrain the luminosity function of cluster red sequence galaxies in the redshift range of $0.1<z<1.05$. We develop a hierarchical Bayesian method that simultaneously models redshift evolution and cluster mass dependence. The results from this method are tested by red sequence luminosity function parameters derived in cluster redshift or mass bins. We find a hint that the faint end slope of a Schechter function fit may evolve with redshift at a significance level of $\sim 1.9 \sigma$. Faint cluster red sequence galaxies possibly become more abundant at lower redshift, indicating a different formation time from the bright red sequence galaxies. Optical cluster cosmology analyses may wish to consider this effect when deriving mass proxies. We also constrain the amplitude of the luminosity function with the hierarchical Bayesian method, which strongly correlates with cluster mass and provides an improved estimate of cluster masses. This technique can be applied to a larger sample of X-ray or optically selected clusters from the Dark Energy Survey, significantly improving the sensitivity of the analysis.
We tested the performance of photometric redshifts for galaxies in the Hubble Ultra Deep field down to 30th magnitude. We compared photometric redshift estimates from three spectral fitting codes from the literature (EAZY, BPZ and BEAGLE) to high quality redshifts for 1227 galaxies from the MUSE integral field spectrograph. All these codes can return photometric redshifts with bias |Dzn|=|z-z_phot|/(1+z)<0.05 down to F775W=30 and spectroscopic incompleteness is unlikely to strongly modify this statement. We have, however, identified clear systematic biases in the determination of photometric redshifts: in the 0.4<z<1.5 range, photometric redshifts are systematically biased low by as much as Dzn=-0.04 in the median, and at z>3 they are systematically biased high by up to Dzn = 0.05, an offset that can in part be explained by adjusting the amount of intergalactic absorption applied. In agreement with previous studies we find little difference in the performance of the different codes, but in contrast to those we find that adding extensive ground-based and IRAC photometry actually can worsen photo-z performance for faint galaxies. We find an outlier fraction, defined through |Dzn|>0.15, of 8% for BPZ and 10% for EAZY and BEAGLE, and show explicitly that this is a strong function of magnitude. While this outlier fraction is high relative to numbers presented in the literature for brighter galaxies, they are very comparable to literature results when the depth of the data is taken into account. Finally, we demonstrate that while a redshift might be of high confidence, the association of a spectrum to the photometric object can be very uncertain and lead to a contamination of a few percent in spectroscopic training samples that do not show up as catastrophic outliers, a problem that must be tackled in order to have sufficiently accurate photometric redshifts for future cosmological surveys.
These lectures, presented at TASI 2016: Anticipating the Next Discoveries in Particle Physics, provide an introduction to some key methods and tools of indirect dark matter searches. Topics covered include estimation of dark matter signals, thermal freezeout and related scenarios, potential effects of dark matter annihilation on the early universe, modeling photon signals from annihilation or decay, and a brief and qualitative introduction to diffusive propagation of cosmic rays. The second half of the notes gives a status report (circa summer 2016) on selected experimental searches, the resulting constraints and some potential signal candidates. These notes are intended as an introduction to indirect dark matter searches for graduate students, focusing on back-of-the-envelope estimates and useful concepts rather than detailed quantitative computations.
I give a brief review of a few recent developments and future directions in the search for dark matter using high-energy neutrinos from the Sun. This includes the ability to recast neutrino telescope limits on nuclear scattering of dark matter to arbitrary new theories, and new calculations of the solar atmospheric background relevant to such searches. I also touch on applications to global searches for new physics, and prospects for improving searches for asymmetric dark matter in the Sun.
We construct average spectra of host galaxies of slower, faster, bluer, and redder Type Ia Supernovae (SNe Ia) from the SDSS-II supernova survey. The average spectrum of slower declining (broader light-curve width or higher stretch) SN Ia hosts shows stronger emission lines compared to the average spectrum of faster declining (narrower light-curve width or lower stretch) SN Ia hosts. Using pPXF, we find that hosts of slower declining SNe Ia have metallicities that are, on average, 0.24 dex higher than average metallicities of faster declining SN Ia hosts. Similarly, redder SN Ia hosts have slightly higher metallicities than bluer SN Ia hosts. Lick index analysis of metallic lines and Balmer lines show that faster declining SN Ia hosts have relatively higher metal content and have relatively older stellar populations compared with slower declining SN Ia hosts. We calculate average $\rm H_{\alpha}$ Star Formation Rate (SFR), stellar mass, and the specific-SFR (sSFR) of host galaxies in these subgroups of SNe Ia. We find that slower declining SN Ia hosts have significantly higher ($>5\sigma$) sSFR than faster declining SN Ia hosts. A Kolmogorov-Smirnov test shows that these two types of hosts originate from different parent distributions. Our results, when compared with the models of \cite{childress14}, indicate that slower declining SNe Ia, being hosted in actively star-forming galaxies, are young (prompt) SNe Ia, originating from similar progenitor age groups.
We present optical and ultraviolet spectra of the first electromagnetic counterpart to a gravitational wave (GW) source, the binary neutron star merger GW170817. Spectra were obtained nightly between 1.5 and 9.5 days post-merger, using the SOAR and Magellan telescopes; the UV spectrum was obtained with the \textit{Hubble Space Telescope} at 5.5 days. Our data reveal a rapidly-fading blue component ($T\approx5500$ K at 1.5 days) that quickly reddens; spectra later than $\gtrsim 4.5$ days peak beyond the optical regime. The spectra are mostly featureless, although we identify a possible weak emission line at $\sim 7900$ \AA\ at $t\lesssim 4.5$ days. The colours, rapid evolution and featureless spectrum are consistent with a "blue" kilonova from polar ejecta comprised mainly of light $r$-process nuclei with atomic mass number $A\lesssim 140$. This indicates a sight-line within $\theta_{\rm obs}\lesssim 45^{\circ}$ of the orbital axis. Comparison to models suggests $\sim0.03$ M$_\odot$ of blue ejecta, with a velocity of $\sim 0.3c$. The required lanthanide fraction is $\sim 10^{-4}$, but this drops to $<10^{-5}$ in the outermost ejecta. The large velocities point to a dynamical origin, rather than a disk wind, for this blue component, suggesting that both binary constituents are neutron stars (as opposed to a binary consisting of a neutron star and a black hole). For dynamical ejecta, the high mass favors a small neutron star radius of $\lesssim 12$ km. This mass also supports the idea that neutron star mergers are a major contributor to $r$-process nucleosynthesis.
We study the diffusion of massive particles in the space time of an Abelian Higgs string. The particles in the early universe plasma execute Brownian motion. This motion of the particles is modeled as a two dimensional random walk in the plane of the Abelian Higgs string. The particles move randomly in the space time of the string according to their geodesic equations. We observe that for certain values of their energy and angular momentum, an overdensity of particles is observed close to the string. We find that the string parameters determine the distribution of the particles. We make an estimate of the density fluctuation generated around the string as a function of the deficit angle. Though the thickness of the string is small, the length is large and the overdensity close to the string may have cosmological consequences in the early universe.
We find vacuum solutions such that massive gravitons are confined in a local spacetime region by their gravitational energy in asymptotically flat spacetimes in the context of the bigravity theory. We call such self-gravitating objects massive graviton geons. The basic equations can be reduced to the Schr\"odinger-Poisson equations with the tensor "wavefunction" in the Newtonian limit. We obtain a non-spherically symmetric solution with $j=2,\ell=0$ as well as a spherically symmetric solution with $j=0,\ell=2$ in this system where $j$ is the total angular momentum quantum number and $\ell$ is the orbital angular momentum quantum number, respectively. The energy eigenvalue of the Schr\"odinger equation in the non-spherical solution is smaller than that in the spherical solution. We then study the perturbative stability of the spherical solution and find that there is an unstable mode in the quadrupole mode perturbations which may be interpreted as the transition mode to the non-spherical solution. The results suggest that the non-spherically symmetric solution is the ground state of the massive graviton geon. The massive graviton geons may decay in time due to emissions of gravitational waves but this timescale can be quite long when the massive gravitons are non-relativistic and then the geons can be long-lived. We also argue possible prospects of the massive graviton geons: applications to the ultralight dark matter scenario, nonlinear (in)stability of the Minkowski spacetime, and a quantum transition of the spacetime.
We consider a renormalizable extension of the minimal supersymmetric standard model endowed by an R and a gauged B - L symmetry. The model incorporates chaotic inflation driven by a quartic potential associated with the Higgs field which leads to a spontaneous breaking of U(1)B-L employing quadratic Kahler potentials with a prominent shift-symmetric part proportional to c- and a tiny violation, proportional to c+, included in a logarithm with prefactor -N<0. It also offers an explanation of mu term of the MSSM provided that one related parameter in the superpotential is somewhat small. Baryogenesis occurs via non-thermal leptogenesis which is realized by the inflaton's decay to the lightest or next-to-lightest right-handed neutrino with masses lower than 1.8x10^13 GeV. Our scenario can be confronted with the current data on the inflationary observables, the baryon asymmetry of the universe, the gravitino limit on the reheating temperature and the data on the neutrino oscillation parameters, for 0.012<=c+/c-<1/N and gravitino as light as 1 TeV.
We study matter-wave interferometry in the presence of a stochastic background of gravitational waves. It is shown that if the background has a scale-invariant spectrum over a wide bandwidth (which is expected in a class of inflationary models of Big Bang cosmology), then separated-path interference cannot be observed for a lump of matter of size above a limit which is very insensitive to the strength and bandwidth of the fluctuations, unless the interferometer is servo-controlled or otherwise protected. For ordinary solid matter this limit is of order 1--10 mm. A servo-controlled or cross-correlated device would also exhibit limits to the observation of macroscopic interference, which we estimate for ordinary matter moving at speeds small compared to c.
During the second observing run of the Laser Interferometer gravitational- wave Observatory (LIGO) and Virgo Interferometer, a gravitational-wave signal consistent with a binary neutron star coalescence was detected on 2017 August 17th (GW170817), quickly followed by a coincident short gamma-ray burst trigger by the Fermi satellite. The Distance Less Than 40 (DLT40) Mpc supernova search performed pointed follow-up observations of a sample of galaxies regularly monitored by the survey which fell within the combined LIGO+Virgo localization region, and the larger Fermi gamma ray burst error box. Here we report the discovery of a new optical transient (DLT17ck, also known as SSS17a; it has also been registered as AT 2017gfo) spatially and temporally coincident with GW170817. The photometric and spectroscopic evolution of DLT17ck are unique, with an absolute peak magnitude of Mr = -15.8 \pm 0.1 and an r-band decline rate of 1.1mag/d. This fast evolution is generically consistent with kilonova models, which have been predicted as the optical counterpart to binary neutron star coalescences. Analysis of archival DLT40 data do not show any sign of transient activity at the location of DLT17ck down to r~19 mag in the time period between 8 months and 21 days prior to GW170817. This discovery represents the beginning of a new era for multi-messenger astronomy opening a new path to study and understand binary neutron star coalescences, short gamma-ray bursts and their optical counterparts.
The historic detection of gravitational waves from a binary neutron star merger (GW170817) and its electromagnetic counterpart led to the first accurate (sub-arcsecond) localization of a gravitational-wave event. The transient was found to be $\sim$10" from the nucleus of the S0 galaxy NGC 4993. We report here the luminosity distance to this galaxy using two independent methods. (1) Based on our MUSE/VLT measurement of the heliocentric redshift ($z_{\rm helio}=0.009783\pm0.000023$) we infer the systemic recession velocity of the NGC 4993 group of galaxies in the cosmic microwave background (CMB) frame to be $v_{\rm CMB}=3231 \pm 53$ km s$^{-1}$. Using constrained cosmological simulations we estimate the line-of-sight peculiar velocity to be $v_{\rm pec}=307 \pm 230$ km s$^{-1}$, resulting in a cosmic velocity of $v_{\rm cosmic}=2924 \pm 236$ km s$^{-1}$ ($z_{\rm cosmic}=0.00980\pm 0.00079$) and a distance of $D_z=40.4\pm 3.4$ Mpc assuming a local Hubble constant of $H_0=73.24\pm 1.74$ km s$^{-1}$ Mpc$^{-1}$. (2) Using Hubble Space Telescope measurements of the effective radius (15.5" $\pm$ 1.5") and contained intensity and MUSE/VLT measurements of the velocity dispersion, we place NGC 4993 on the Fundamental Plane (FP) of E and S0 galaxies. Comparing to a frame of 10 clusters containing 226 galaxies, this yields a distance estimate of $D_{\rm FP}=44.0\pm 7.5$ Mpc. The combined redshift and FP distance is $D_{\rm NGC 4993}= 41.0\pm 3.1$ Mpc. This 'electromagnetic' distance estimate is consistent with the independent measurement of the distance to GW170817 as obtained from the gravitational-wave signal ($D_{\rm GW}= 43.8^{+2.9}_{-6.9}$ Mpc) and confirms that GW170817 occurred in NGC 4993.
Recently, the optical counterpart of a gravitational wave source GW170817 has been identified in NGC 4993 galaxy. Together with evidence from observations in electromagnetic waves, the event has been suggested as a result of a merger of two neutron stars. We analyze the multi-wavelength data to characterize the host galaxy property and its distance to examine if the properties of NGC 4993 are consistent with this picture. Our analysis shows that NGC 4993 is a bulge-dominated galaxy with reff ~ 2-3 kpc and the Sersic index of n = 3-4 for the bulge component. The spectral energy distribution from 0.15 to 24 micron indicates that this galaxy has no significant ongoing star formation, the mean stellar mass of (0.3 - 1.2) times 10^11 Msun,the mean stellar age greater than ~3 Gyr, and the metallicity of about 20% to 100% of solar abundance. Optical images reveal dust lanes and extended features that suggest a past merging activity. Overall, NGC 4993 has characteristics of normal, but slightly disturbed elliptical galaxies. Furthermore, we derive the distance to NGC 4993 with the fundamental plane relation using 17 parameter sets of 7 different filters and the central stellar velocity dispersion from literature, finding an angular diameter distance of 37.7 +- 8.7 Mpc. NGC 4993 is similar to some host galaxies of short gamma-ray bursts but much different from those of long gamma-ray bursts, supporting the picture of GW170817 as a result of a merger of two NSs.
Links to: arXiv, form interface, find, astro-ph, recent, 1710, contact, help (Access key information)
Time delay lensing is a mature and competitive cosmological probe. However, it is limited in accuracy by the well-known problem of the mass-sheet degeneracy: too rigid assumptions on the density profile of the lens can potentially bias the inference on cosmological parameters. I investigate the degeneracy between the choice of the lens density profile and the inference on the Hubble constant, focusing on double image systems. By expanding lensing observables in terms of the local derivatives of the lens potential around the Einstein radius, and assuming circular symmetry, I show that three degrees of freedom in the radial direction are necessary to achieve a few percent accuracy in the time-delay distance. Additionally, while the time delay is strongly dependent on the second derivative of the potential, observables typically used to constrain lens models in time-delay studies, such as image position and radial magnification information, are mostly sensitive to the first and third derivatives, making it very challenging to accurately determine time-delay distances with lensing data alone. Tests on mock observations show that the assumption of a power-law density profile results in a 5% average bias on $H_0$, with a 6% scatter. Using a more flexible model and adding unbiased velocity dispersion constraints allows to obtain an inference with 1% accuracy. A power-law model can still provide 3% accuracy if velocity dispersion measurements are used to constrain its slope. Although this work is based on the assumption of axisymmetry, its main findings can be generalized to cases with moderate ellipticity.
The relationship between the integrated H$\beta$ line luminosity and the velocity dispersion of the ionized gas of HII galaxies and giant HII regions represents an exciting standard candle that presently can be used up to redshifts z ~ 4. Locally it is used to obtain precise measurements of the Hubble constant by combining the slope of the relation obtained from nearby ($z \leq $ 0.2) HII galaxies with the zero point determined from giant HII regions belonging to an `anchor sample' of galaxies for which accurate redshift-independent distance moduli are available. We present new data for 36 giant HII regions in 13 galaxies of the anchor sample that includes the megamaser galaxy NGC 4258. Our data is the result of the first four years of observation of our primary sample of 130 giant HII regions in 73 galaxies with Cepheid determined distances. Our best estimate of the Hubble parameter is $71.0\pm2.8(random)\pm2.1(systematic)$ km /s Mpc This result is the product of an independent approach and, although at present less precise than the latest SNIa results, it is amenable to substantial improvement.
The peculiar velocity field can be used to study the large-scale distribution of matter in the local universe and test cosmological models. However, present measurements of peculiar velocities are based on empirical distance indicators, which introduce large uncertainties. In this paper, we propose a method to derive the peculiar velocities, which bases on the distances measured from gravitational waves and the redshifts inferred from the host galaxies. Using the first gravitational wave event GW170817 from binary neutron star merger, we find that the peculiar velocity of the host galaxy NGC 4993 is $-275~\rm km~s^{-1}$, if the Hubble constant from Planck Collaboration is used. In future, with the uncertainty of the distance of GW events reducing to $0.1 \%$, the uncertainty of the peculiar velocity can be reduced to $\sim 10$ km/s at 100 Mpc. With accumulated GW events being observed, we can reconstruct the peculiar velocity field, which can be used to test the gravitational-instability paradigm and the $\Lambda$CDM model.
We study the dynamics of cosmological phase transitions in the case of small velocities of bubble walls, $v_w<0.1$. We discuss the conditions in which this scenario arises in a physical model, and we compute the development of the phase transition. We consider different kinds of approximations and refinements for relevant aspects of the dynamics, such as the dependence of the wall velocity on hydrodynamics, the distribution of the latent heat, and the variation of the nucleation rate. Although in this case the common simplifications of a constant wall velocity and an exponential nucleation rate break down due to reheating, we show that a delta-function rate and a velocity which depends linearly on the temperature give a good description of the dynamics and allow to solve the evolution analytically. We also consider a Gaussian nucleation rate, which gives a more precise result for the bubble size distribution. We discuss the implications for the computation of cosmic remnants.
The recent observations of gravitational-wave and electromagnetic emission produced by the merger of the binary neutron-star system GW170817 have opened the possibility of using standard siren to constrain the value of the Hubble constant. While the reported bound is significantly weaker than those recently derived by usual luminosity distances methods, they do not require any form of cosmic distance ladder and can be considered as complementary and, in principle, more conservative. Here we combine the new measurement with the Planck Cosmic Microwave Background observations in a 12 parameters extended LambdaCDM scenario, where the Hubble constant is weakly constrained by CMB data and bound to a low value $H_0=55^{+7}_{-20}$ km/s/Mpc at 68 % C.L. The non-Gaussian shape of the GW170817 bound makes lower values of the Hubble constant in worst agreement with observations. The inclusion of the new GW170817 Hubble constant measurement significantly reduces the allowed parameter space, improving the cosmological bounds on several parameters as the neutrino mass, curvature and the dark energy equation of state.
The detection of an electromagnetic counterpart (GRB 170817A) to the gravitational wave signal (GW170817) from the merger of two neutron stars opens a completely new arena for testing theories of gravity. We show that this measurement allows us to place stringent constraints on general scalar-tensor and vector-tensor theories, while allowing us to place an independent bound on the graviton mass in bimetric theories of gravity. These constraints severely reduce the viable range of cosmological models that have been proposed as alternatives to general relativistic cosmology.
Following the detection of the GW170817 signal and its associated electromagnetic emissions, we discuss the prospects of the local Hubble parameter measurement using double neutron stars (DNSs). The kilonova emissions of GW170817 are genuinely unique in terms of the rapid evolution of color and magnitude and we expect that, for a good fraction $\gtrsim 50\%$ of the DNS events within $\sim 200$Mpc, we could identify their host galaxies, using their kilonovae. At present, the estimated DNS merger rate $(1.5^{+3.2}_{-1.2})\times 10^{-6} {\rm Mpc^{-3} yr^{-1}}$ has a large uncertainty. But, if it is at the high end, we could measure the local Hubble parameter $H_L$ with the level of $\Delta H_L/H_L\sim 0.042$ ($1\sigma$ level), after the third observational run (O3). This accuracy is four times better than that obtained from GW170817 alone, and we will be able to examine the Hubble tension at $2.1\sigma$ level.
The luminosity distance measurement of GW170817 derived from the gravitational wave analysis in Abbott et al. 2017 (here, A17:H0) is highly correlated with the measured inclination of the NS-NS system. To improve the precision of the distance measurement, we attempt to constrain the inclination by modeling the broad-band X-ray-to-radio emission from GW170817, which is dominated by the interaction of the jet with the environment. We update our previous analysis of the GW170817 afterglow and we consider the entire radio and X-ray data set obtained at $t<40$ days since merger. We find that the broad-band X-ray to radio emission from GW170817 is consistent with an off-axis relativistic jet with kinetic energy $10^{48}\,\rm{erg}<E_{k}\le 3\times 10^{50} \,\rm{erg}$ propagating into an environment with circumbinary density $n\sim10^{-2}-10^{-4} \,\rm{cm^{-3}}$, with a preference for wider jets with opening angle $\theta_j=15$ deg. For these jets, our modeling indicates an off-axis angle $\theta_{\rm obs}\sim25-50$ deg. We combine our constraints on $\theta_{\rm obs}$ with the joint distance-inclination constraint from A17:H0. Using the same $\sim 170$ km/sec peculiar velocity uncertainty assumed in A17:H0 but with an inclination constraint from the radio and X-ray data, we get a value of $H_0=$$74.0 \pm \frac{13.7}{5.3}$ km/s/Mpc, which is higher than the value of $H_0=$$70.0 \pm \frac{12.0}{8.0}~$ km/s/Mpc found in A17:H0. Further, using a more realistic peculiar velocity uncertainty of 250 km/sec derived from previous work, we find $H_0=$$75.5 \pm \frac{14.0}{7.3}$ km/s/Mpc for H0 from this system. We note that this is in modestly better agreement with the local distance ladder than the Planck CMB, though a significant such discrimination will require $\sim 50$ such events. Future measurements at $t>100$ days of the X-ray and radio emission will lead to tighter constraints.
We study the equations of conformal gravity, as given by Mannheim, in the weak field limit, so that a linear approximation is adequate. Specializing to static fields with spherical symmetry, we obtain a second-order equation for one of the metric functions. We obtain the Green function for this equation, and represent the metric function in the form of integrals over the source. Near a compact source such as the Sun the solution no longer has Schwarzschild form. Using Flanagan's method of obtaining a conformally invariant metric tensor we attempt to get a solution of Schwarzschild type. We find, however, that the 1/r terms disappear altogether. We conclude that a solution of Mannheim type cannot exist for these field equations.
We make a theoretical prediction for the ratio of the dark energy to other components in the Universe based on the scenario of the sequestering mechanism which was recently proposed as one possible way to solve the cosmological constant problem. In order to evaluate the value of dark energy and the others, we assume a specific scale factor which describes the Big-Crunch scenario in the scalar-tensor theory. We specify the parameter region where the one can explain the observed dark energy-matter ratio.
On August 17, 2017 the LIGO interferometers detected the gravitational wave (GW) signal (GW170817) from the coalescence of binary neutron stars. This signal was also simultaneously seen throughout the electromagnetic (EM) spectrum from radio waves to gamma-rays. We point out that this simultaneous detection of GW and EM signals rules out a class of modified gravity theories, which dispense with the need for dark matter. This simultaneous observation also provides the first ever test of Einstein's Weak Equivalence Principle (WEP) between gravitons and photons. We calculate the Shapiro time delay due to the gravitational potential of the total dark matter distribution along the line of sight (complementary to the calculation in arXiv:1710.05834) to be about 1000 days. Using this estimate for the Shapiro delay and from the time difference of 1.7 seconds between the GW signal and gamma-rays, we can constrain violations of WEP using the parameterized post-Newtonian (PPN) parameter $\gamma$, and is given by $|\gamma_{\rm {GW}}-\gamma_{\rm{EM}}|<3.9 \times 10^{-8}$.
We discuss the possibility of realising a multi-component dark matter scenario with widely separated dark matter masses: one having keV scale mass and the other with GeV-TeV scale mass, within the framework of left right symmetric models. To reduce the level of fine tuning, we consider a version of left right model with universal seesaw where keV-GeV scale right handed neutrinos having tiny mixing with light neutrinos can be naturally realised. Due to gauge interactions, both the dark matter candidates are produced thermally in the early Universe but overproducing the keV mass candidate. We consider one of the right handed neutrinos to be decaying at late epochs, just before the big bang nucleosynthesis, in order to dilute the thermally overproduced keV dark matter. We constrain the parameter space from the requirement of producing sub-dominant keV-TeV dark matter, satisfying indirect detection constraints from gamma ray searches and producing the tantalising 3.55 keV monochromatic X-ray line, reported by several groups to be present in galaxy and galaxy cluster data, from the decay of a 7.1 keV dark matter on cosmological scales. We find that these requirements constrain the right sector gauge boson masses to be much heavier than the present collider sensitivity, but can have interesting signatures at experiments like the IceCube.
Lyman-$\alpha$ emitters (LAEs) are a promising probe of the large-scale structure at high redshift, $z\gtrsim 2$. In particular, the Hobby-Eberly Telescope Dark Energy Experiment aims at observing LAEs at 1.9 $<z<$ 3.5 to measure the Baryon Acoustic Oscillation (BAO) scale and the Redshift-Space Distortion (RSD). However, Zheng et al. (2011) pointed out that the complicated radiative transfer (RT) of the resonant Lyman-$\alpha$ emission line generates an anisotropic selection bias in the LAE clustering on large scales, $s\gtrsim 10$ Mpc. This effect could potentially induce a systematic error in the BAO and RSD measurements. Also, Croft et al. (2016) claims an observational evidence of the effect in the Lyman-$\alpha$ intensity map, albeit statistically insignificant. We aim at quantifying the impact of the Lyman-$\alpha$ RT on the large-scale galaxy clustering in detail. For this purpose, we study the correlations between the large-scale environment and the ratio of an apparent Lyman-$\alpha$ luminosity to an intrinsic one, which we call the `observed fraction', at $2<z<6$. We apply our Lyman-$\alpha$ RT code by post-processing the full Illustris simulations. We simply assume that the intrinsic luminosity of the Lyman-$\alpha$ emission is proportional to the star formation rate of galaxies in Illustris, yielding a sufficiently large sample of LAEs to measure the anisotropic selection bias. We find little correlations between large-scale environment and the observed fraction induced by the RT, and hence a smaller anisotropic selection bias than what was claimed by Zheng et al. (2011). We argue that the anisotropy was overestimated in the previous work due to the insufficient spatial resolution: it is important to keep the resolution such that it resolves the high density region down to the scale of the interstellar medium, $\sim1$ physical kpc. (abridged)
Links to: arXiv, form interface, find, astro-ph, recent, 1710, contact, help (Access key information)
Measuring and calibrating relations between cluster observables is critical for resource-limited studies. The mass-richness relation of clusters offers an observationally inexpensive way of estimating masses. Its calibration is essential for cluster and cosmological studies, especially for high-redshift clusters. Weak gravitational lensing magnification is a promising and complementary method to shear studies, that can be applied at higher redshifts. We employed the weak lensing magnification method to calibrate the mass-richness relation up to a redshift of 1.4. We used the Spitzer Adaptation of the Red-Sequence Cluster Survey (SpARCS) galaxy cluster candidates ($0.2<z<1.4$) and optical data from the Canada France Hawaii Telescope (CFHT) to test whether magnification can be effectively used to constrain the mass of high-redshift clusters. Lyman-Break Galaxies (LBGs) selected using the $u$-band dropout technique and their colours were used as a background sample of sources. LBG positions were cross-correlated with the centres of the sample of SpARCS clusters to estimate the magnification signal measured for cluster sub-samples, binned in both redshift and richness. We detected a weak lensing magnification signal for all bins at a detection significance of 2.6-5.5$\sigma$. In particular, the significance of the measurement for clusters with $z>1.0$ is 4.1$\sigma$; for the entire cluster sample we obtained an average M$_{200}$ of $1.28^{+0.23}_{-0.21}$ $\times 10^{14} \, \textrm{M}_{\odot}$. Our measurements demonstrated the feasibility of using weak lensing magnification as a viable tool for determining the average halo masses for samples of high redshift galaxy clusters. The results also established the success of using galaxy over-densities to select massive clusters at $z > 1$. Additional studies are necessary for further modelling of the various systematic effects we discussed.
The extragalactic background radiation produced by distant galaxies emitting in the far infrared limits the sensitivity of telescopes operating in this range due to confusion. We have constructed a model of the infrared background based on numerical simulations of the large-scale structure of the Universe and the evolution of dark matter halos. The predictions of this model agree well with the existing data on source counts. We have constructed maps of a sky field with an area of 1 deg$^2$ directly from our simulated observations and measured the confusion limit. At wavelengths $100-300$ $\mu$m the confusion limit for a 10-m telescope has been shown to be at least an order of magnitude lower than that for a 3.5-m one. A spectral analysis of the simulated infrared background maps clearly reveals the large-scale structure of the Universe. The two-dimensional power spectrum of these maps has turned out to be close to that measured by space observatories in the infrared. However, the fluctuations in the number of intensity peaks observed in the simulated field show no clear correlation with superclusters of galaxies; the large-scale structure has virtually no effect on the confusion limit.
The nonlinear interaction between one graviton and two scalars is enhanced in specific inflationary models, potentially leading to distinguishable signatures in the bispectrum of the cosmic microwave background (CMB) anisotropies. We develop the tools to examine such bispectrum signatures, and show a first application using WMAP temperature data. We consider several $\ell$-ranges, estimating the $g_{tss}$ amplitude parameter, by means of the so-called separable modal methodology. We do not find any evidence of a tensor-scalar-scalar signal at any scale. Our tightest bound on the size of the tensor-scalar-scalar correlator is derived from our measurement including all the multipoles in the range $ 2 \leq \ell \leq 500$ and it reads $g_{tss} = -48 \pm 28$ ($68\%$CL). This is the first direct observational constraint on the primordial tensor-scalar-scalar correlation, and it will be cross-checked and improved by applying the same pipeline to high-resolution temperature and polarization data from $Planck$ and forthcoming CMB experiments.
The Cosmic Neutrino Background is a prediction of the standard cosmological model, but it has been never observed directly. In the experiments with the aim of detecting relic CNB neutrinos, currently under development, the expected event rate depends on the local density of relic neutrinos. Since massive neutrinos can be attracted by the gravitational potential of our galaxy and cluster around it, a local overdensity of cosmic neutrinos should exist. Considering the minimal masses guaranteed by neutrino oscillations, we review the computation of the local density of relic neutrinos and we present realistic prospects for a PTOLEMY-like experiment.
Galaxy clusters are the most massive constituents of the large-scale structure of the Universe. While the hot thermal gas that pervades galaxy clusters is relatively well understood through observations with X-ray satellites, our understanding of the non-thermal part of the intra-cluster medium remains incomplete. With LOFAR and GMRT observations, we have identified a phenomenon that can be unveiled only at extremely low radio-frequencies and offers new insights into the non-thermal component. We propose that the interplay between radio-emitting plasma and the perturbed intra-cluster medium can gently re-energise relativistic particles initially injected by active galactic nuclei. Sources powered through this mechanism can maintain electrons at higher energies than radiative ageing would allow. If this mechanism is common for aged plasma, a population of mildly relativistic electrons can be accumulated inside galaxy clusters providing the seed population for merger-induced re-acceleration mechanisms on larger scales such as turbulence and shock waves.
Splashback refers to the process of matter that is accreting onto a dark matter halo reaching its first orbital apocenter and turning around in its orbit. The cluster-centric radius at which this process occurs, r_sp, defines a halo boundary that is connected to the dynamics of the cluster, in contrast with other common halo boundary definitions such as R_200. A rapid decline in the matter density profile of the halo is expected near r_sp. We measure the galaxy number density and weak lensing mass profiles around RedMapper galaxy clusters in the first year Dark Energy Survey (DES) data. For a cluster sample with mean mass ~2.5 x 10^14 solar masses, we find strong evidence of a splashback-like steepening of the galaxy density profile and measure r_sp=1.16 +/- 0.08 Mpc/h, consistent with earlier SDSS measurements of More et al. (2016) and Baxter et al. (2017). Moreover, our weak lensing measurement demonstrates for the first time the existence of a splashback-like steepening of the matter profile of galaxy clusters. We measure r_sp=1.28 +/- 0.18 Mpc/h from the weak lensing data, in good agreement with our galaxy density measurements. Applying our analysis to different cluster and galaxy samples, we find that consistent with LambdaCDM simulations, r_sp scales with R_200m and does not evolve with redshift over the redshift range of 0.3--0.6. We also find that potential systematic effects associated with the RedMapper algorithm may impact the location of r_sp, in particular the choice of scale used to estimate cluster richness. We discuss progress needed to understand the systematic uncertainties and fully exploit forthcoming data from DES and future surveys, emphasizing the importance of more realistic mock catalogs and independent cluster samples.
Homogeneous oscillations of the inflaton after inflation can be unstable to small spatial perturbations even without coupling to other fields. We show that for inflaton potentials $\propto |\phi|^{2n}$ near $|\phi|=0$ and flatter beyond some $|\phi|=M$, the inflaton condensate oscillations can lead to self-resonance, followed by its complete fragmentation. We find that for non-quadratic minima ($n>1$), shortly after backreaction, the equation of state parameter, $w\rightarrow1/3$. If $M\ll m_{pl}$, radiation domination is established within less than an e-fold of expansion after the end of inflation. In this case self-resonance is efficient and the condensate fragments into transient, localised spherical objects which are unstable and decay, leaving behind them a virialized field with mean kinetic and gradient energies much greater than the potential energy. This end-state yields $w=1/3$. When $M\ll m_{pl}$ we observe slow and steady, self-resonace that can last many {\it e}-folds before backreaction eventually shuts it off, followed by fragmentation and $w\rightarrow 1/3$. We provide analytical estimates for the duration to $w\rightarrow 1/3$ after inflation, which can be used as an upper bound (under certain assumptions) on the duration of the transition between the inflationary and the radiation dominated states of expansion. This upper bound can reduce uncertainties in CMB observables such as the spectral tilt $n_{\rm{s}}$, and the tensor-to-scalar ratio $r$. For quadratic minima ($n=1$), $w\rightarrow0$ regardless of the value of $M$.
We highlight the general scenario of dark matter freeze-out whilst the energy density of the universe is dominated by a decoupled non-relativistic species. Decoupling during matter domination changes the freeze-out dynamics, since the Hubble rate is parametrically different for matter and radiation domination. Furthermore, for successful Big Bang Nucleosynthesis the state dominating the early universe energy density must decay, this dilutes (or repopulates) the dark matter. As a result, the masses and couplings required to reproduce the observed dark matter relic density can differ significantly from radiation dominated freeze-out.
In this paper, we show how the proper choice of gauge is critical in analyzing the stability of non-singular cosmological bounce solutions based on Horndeski theories. We show that it is possible to construct non-singular cosmological bounce solutions with classically stable behavior for all modes with wavelengths above the Planck scale where: (a) the solution involves a stage of null-energy condition violation during which gravity is described by a modification of Einstein's general relativity; and (b) the solution reduces to Einstein gravity both before and after the null-energy condition violating stage. Similar considerations apply to galilean genesis scenarios.
To model a realistic situation for the beginning we consider massive real scalar $\phi^4$ theory in a (1+1)-dimensional asymptotically static Minkowski spacetime with an intermediate stage of expansion. To have an analytic headway we assume that scalars have a big mass. At past and future infinities of the background we have flat Minkowski regions which are joint by the inflationary expansion region. We use the tree-level Keldysh propagator in the theory in question to calculate the expectation value of the stress-energy tensor which is, thus, due to the excitations of the zero-point fluctuations. Then we show that even for large mass, if the de Sitter expansion stage is long enough, the quantum loop corrections to the expectation value of the stress-energy tensor are not negligible in comparison with the tree-level contribution. That is revealed itself via the excitation of the higher-point fluctuations of the exact modes: During the expansion stage a non-zero particle number density for the exact modes is generated. This density is not Plankian and serves as a quench which leads to a thermalization in the out Minkowski stage.
We use the energy-balance code MAGPHYS to determine stellar and dust masses, and dust corrected star-formation rates for over 200,000 GAMA galaxies, 170,000 G10-COSMOS galaxies and 200,000 3D-HST galaxies. Our values agree well with previously reported measurements and constitute a representative and homogeneous dataset spanning a broad range in stellar mass (10^8---10^12 Msol), dust mass (10^6---10^9 Msol), and star-formation rates (0.01---100 Msol per yr), and over a broad redshift range (0.0 < z < 5.0). We combine these data to measure the cosmic star-formation history (CSFH), the stellar-mass density (SMD), and the dust-mass density (DMD) over a 12 Gyr timeline. The data mostly agree with previous estimates, where they exist, and provide a quasi-homogeneous dataset using consistent mass and star-formation estimators with consistent underlying assumptions over the full time range. As a consequence our formal errors are significantly reduced when compared to the historic literature. Integrating our cosmic star-formation history we precisely reproduce the stellar-mass density with an ISM replenishment factor of 0.50 +/- 0.07, consistent with our choice of Chabrier IMF plus some modest amount of stripped stellar mass. Exploring the cosmic dust density evolution, we find a gradual increase in dust density with lookback time. We build a simple phenomenological model from the CSFH to account for the dust mass evolution, and infer two key conclusions: (1) For every unit of stellar mass which is formed 0.0065---0.004 units of dust mass is also formed; (2) Over the history of the Universe approximately 90 to 95 per cent of all dust formed has been destroyed and/or ejected.
We present a study of NGC 4993, the host galaxy of the GW170817 gravitational wave event, the GRB170817A short gamma-ray burst (sGRB) and the AT2017gfo kilonova. We use Dark Energy Camera imaging, AAT spectra and publicly available data, relating our findings to binary neutron star (BNS) formation scenarios and merger delay timescales. NGC4993 is a nearby (40 Mpc) early-type galaxy, with $i$-band S\'ersic index $n=4.0$ and low asymmetry ($A=0.04\pm 0.01$). These properties are unusual for sGRB hosts. However, NGC4993 presents shell-like structures and dust lanes indicative of a recent galaxy merger, with the optical transient located close to a shell. We constrain the star formation history (SFH) of the galaxy assuming that the galaxy merger produced a star formation burst, but find little to no on-going star formation in either spatially-resolved broadband SED or spectral fitting. We use the best-fit SFH to estimate the BNS merger rate in this type of galaxy, as $R_{NSM}^{gal}= 5.7^{+0.57}_{-3.3} \times 10^{-6} {\rm yr}^{-1}$. If star formation is the only considered BNS formation scenario, the expected number of BNS mergers from early-type galaxies detectable with LIGO during its first two observing seasons is $0.038^{+0.004}_{-0.022}$, as opposed to $\sim 0.5$ from all galaxy types. Hypothesizing that the binary system formed due to dynamical interactions during the galaxy merger, the subsequent time elapsed can constrain the delay time of the BNS coalescence. By using velocity dispersion estimates and the position of the shells, we find that the galaxy merger occurred $t_{\rm mer}\lesssim 200~{\rm Myr}$ prior to the BNS coalescence.
We study the evolution of Gravitational Waves (GWs) during and after inflation as well as the resulting observational consequences in a Lorentz-violating massive gravity theory with one scalar (inflaton) and two tensor degrees of freedom. We consider two explicit examples of the tensor mass $m_g$ that depends either on the inflaton field $\phi$ or on its time derivative $\dot{\phi}$, both of which lead to parametric excitations of GWs during reheating after inflation. The first example is Starobinsky's $R^2$ inflation model with a $\phi$-dependent $m_g$ and the second is a low-energy-scale inflation model with a $\dot{\phi}$-dependent $m_g$. We compute the energy density spectrum $\Omega_{\rm GW}(k)$ today of the GW background. In the Starobinsky's model, we show that the GWs can be amplified up to the detectable ranges of both CMB and DECIGO. In low-scale inflation with a fast transition to the reheating stage driven by the potential $V(\phi)=M^2 \phi^2/2$ around $\phi \approx M_{\rm pl}$ (where $M_{\rm pl}$ is the reduced Planck mass), we find that the peak position of $\Omega_{\rm GW}(k)$ induced by the parametric resonance can reach the sensitivity region of advanced LIGO for the Hubble parameter of order 1 GeV at the end of inflation. Thus, our massive gravity scenario offers exciting possibilities for probing the physics of primordial GWs at various different frequencies.
Links to: arXiv, form interface, find, astro-ph, recent, 1710, contact, help (Access key information)
We investigate the effect of the accelerated expansion of the Universe due to a cosmological constant, $\Lambda$, on the cosmic star formation rate. We utilise hydrodynamical simulations from the EAGLE suite, comparing a $\Lambda$CDM Universe to an Einstein-de Sitter model with $\Lambda=0$. Despite the differences in the rate of growth of structure, we find that dark energy, at its observed value, has negligible impact on star formation in the Universe. We study these effects beyond the present day by allowing the simulations to run forward into the future ($t>13.8$ Gyr). We show that the impact of $\Lambda$ becomes significant only when the Universe has already produced most of its stellar mass, only decreasing the total co-moving density of stars ever formed by ${\approx}15\%$. We develop a simple analytic model for the cosmic star formation rate that captures the suppression due to a cosmological constant. The main reason for the similarity between the models is that feedback from accreting black holes dramatically reduces the cosmic star formation at late times. Interestingly, simulations without feedback from accreting black holes predict an upturn in the cosmic star formation rate for $t>15$ Gyr due to the rejuvenation of massive ($ > 10^{11} \mathrm{M}_{\odot}$) galaxies. We briefly discuss the implication of the weak dependence of the cosmic star formation on $\Lambda$ in the context of the anthropic principle.
If primordial black holes (PBHs) form directly from inhomogeneities in the early universe, then the number in the mass range $10^6 -10^9M_{\odot}$ is severely constrained by upper limits to the $\mu$-distortion in the cosmic microwave background (CMB). This is because inhomogeneities on these scales will be dissipated by Silk damping in the redshift interval $5\times 10^4\lesssim z\lesssim2\times 10^6$. If the primordial fluctuations on a given mass-scale have a Gaussian distribution and PBHs form on the high-$\sigma$ tail, as in the simplest scenarios, then the $\mu$ constraints exclude PBHs in this mass range from playing any interesting cosmological role. Only if the fluctuations are highly non-Gaussian, or form through some mechanism unrelated to the primordial fluctuations, can this conclusion be obviated.
Type Ia supernovae (SNIe) are generally accepted to act as standardisable candles, and their use in cosmology led to the first confirmation of the as yet unexplained accelerated cosmic expansion. Many of the theoretical models to explain the cosmic acceleration assume modifications to Einsteinian General Relativity which accelerate the expansion, but the question of whether such modifications also affect the ability of SNIe to be standardisable candles has rarely been addressed. This paper is an attempt to answer this question. For this we adopt a semi-analytical model to calculate SNIe light curves in non-standard gravity. We use this model to show that the average rescaled intrinsic peak luminosity -- a quantity that is assumed to be constant with redshift in standard analyses of Type Ia supernova (SNIa) cosmology data -- depends on the strength of gravity in the supernova's local environment because the latter determines the Chandrasekhar mass -- the mass of the SNIa's white dwarf progenitor right before the explosion. This means that SNIe are no longer standardisable candles in scenarios where the strength of gravity evolves over time, and therefore the cosmology implied by the existing SNIa data will be different when analysed in the context of such models. As an example, we show that the observational SNIa cosmology data can be fitted with both a model where $(\Omega_{\rm M}, \Omega_{\Lambda})=(0.62, 0.38)$ and Newton's constant $G$ varies as $G(z)=G_0(1+z)^{-1/4}$ and the standard model where $(\Omega_{\rm M}, \Omega_{\Lambda})=(0.3, 0.7)$ and $G$ is constant, when the Universe is assumed to be flat.
We present an optimization scheme for the focal plane of a multi-frequency CMB B-modes experiment with a fixed number of detectors and apply it to the specific case of LSPE experiment. Optimal focal planes are identified on the ground of different figures of merit defined in terms of the forecasted uncertainty $\sigma_r$ on the tensor--to-scalar ratio $r$ and the expected map variance from foreground and instrumental noise residuals. We then perform a forecast analysis in order to assess the precision achievable in B-modes measurements.
Using the DAMA/LIBRA anomaly as an example, we formalise the notion of halo-independence in the context of Bayesian statistics and quantified maximum entropy. We consider an infinite set of possible profiles, weighted by an entropic prior and constrained by a likelihood describing noisy measurements of modulated moments by DAMA/LIBRA. Assuming an isotropic dark matter (DM) profile in the galactic rest frame, we find the most plausible DM profiles and predictions for unmodulated signal rates at DAMA/LIBRA. The entropic prior contains an a priori unknown regularisation factor, $\beta$, that describes the strength of our conviction that the profile is approximately Maxwellian. By varying $\beta$, we smoothly interpolate between a halo-independent and a halo-dependent analysis, thus exploring the impact of prior information about the DM profile.
Voronoi grids have been successfully used to represent density structures of gas in astronomical hydrodynamics simulations. While some codes are explicitly built around using a Voronoi grid, others, such as Smoothed Particle Hydrodynamics (SPH), use particle-based representations and can benefit from constructing a Voronoi grid for post-processing their output. So far, calculating the density of each Voronoi cell from SPH data has been done numerically, which is both slow and potentially inaccurate. This paper proposes an alternative analytic method, which is fast and accurate. We derive an expression for the integral of a cubic spline kernel over the volume of a Voronoi cell and link it to the density of the cell. Mass conservation is ensured rigorously by the procedure. The method can be applied more broadly to integrate a spherically symmetric polynomial function over the volume of a random polyhedron.
The axion causes two serious cosmological problems, domain wall and isocurvature perturbation problems. Linde pointed that the isocurvature perturbations are suppressed when the Peccei-Quinn (PQ) scalar field takes a large value $\sim M_{\rm pl}$ (Planck scale) during inflation. In this case, however, the PQ field with large amplitude starts to oscillate after inflation and large fluctuations of the PQ field are produced through parametric resonance, which leads to the formation of domain walls. We consider a supersymmetric axion model and examine whether domain walls are formed by using lattice simulation. It is found that the domain wall problem does not appear in the SUSY axion model when the initial value of the PQ field is less than $10^3\times v$ where $v$ is the PQ symmetry breaking scale.
Links to: arXiv, form interface, find, astro-ph, recent, 1710, contact, help (Access key information)