Links to: arXiv, form interface, find, astro-ph, recent, 1712, contact, help (Access key information)
We provide a general procedure for optimally compressing $N$ data down to $n$ summary statistics, where $n$ is equal to the number of parameters of interest. We show that compression to the score function -- the gradient of the log-likelihood with respect to the parameters -- yields $n$ compressed statistics that are optimal in the sense that they preserve the Fisher information content of the data. Our method generalizes earlier work on linear Karhunen-Lo\'{e}ve compression for Gaussian data whilst recovering both lossless linear compression and quadratic estimation as special cases when they are optimal. We give a unified treatment that also includes the general non-Gaussian case as long as mild regularity conditions are satisfied, producing optimal non-linear summary statistics when appropriate. As a worked example, we derive explicitly the $n$ optimal compressed statistics for Gaussian data in the general case where both the mean and covariance depend on the parameters.
We study reionization in two non-flat $\Lambda$CDM inflation models that best fit the Planck 2015 cosmic microwave background anisotropy observations, ignoring or in conjunction with baryon acoustic oscillation distance measurements. We implement a principal component analysis (PCA) to estimate the uncertainties in the reionization history from a joint quasar-CMB dataset. A thorough Markov Chain Monte Carlo analysis is done over the parameter space of PCA modes for both non-flat $\Lambda$CDM inflation models as well as the original Planck 2016 tilted, spatially-flat $\Lambda$CDM inflation model. Although both flat and non-flat models can closely match the low-redshift ($z\lesssim6$) observations, we notice a possible tension between high-redshift ($z\sim8$) Lyman-$\alpha$ emitter data and the non-flat models. This is solely due to the fact that the closed models have a relatively higher reionization optical depth compared to the flat one, which in turn demands more high-redshift ionizing sources and favors an extended reionization starting as early as $z\approx14$. We also indicate the promising prospects of near-future surveys to decisively resolve this issue.
We apply a friends-of-friends algorithm to the HectoMAP redshift survey and cross-identify associated X-ray emission in the ROSAT All-Sky Survey data (RASS). The resulting flux limited catalog of X-ray cluster survey is complete to a limiting flux of $\sim3 \times10^{-13}$ erg s$^{-1}$ cm$^{-2}$ and includes 15 clusters (7 newly discovered) with redshift $z \leq 0.4$. HectoMAP is a dense survey ($\sim1200$ galaxies deg$^{-2}$) that provides $\sim50$ members (median) in each X-ray cluster. We provide redshifts for the 1036 cluster members. Subaru/Hyper Suprime-Cam imaging covers three of the X-ray systems and confirms that they are impressive clusters. The HectoMAP X-ray clusters have an $L_{X} - {\sigma}_{cl}$ scaling relation similar to that of known massive X-ray clusters. The HectoMAP X-ray cluster sample predicts $\sim 12000 \pm3000$ detectable X-ray clusters in the RASS to the limiting flux, comparable with previous estimates.
We compute a general expression for the contribution of vector perturbations to the redshift-space distortion of galaxy surveys. We show that they contribute to the same multipoles of the correlation function as scalar perturbations and should thus in principle be taken into account in data analysis. We derive constraints for next-generation surveys on the amplitude of two sources of vector perturbations, namely non-linear clustering and topological defects. While topological defects leave a very small imprint on redshift-space distortions, we show that the multipoles of the correlation function are sensitive to vorticity induced by non-linear clustering. Therefore future redshift surveys such as DESI or the SKA should be capable of measuring such vector modes, especially with the hexadecapole which appears to be the most sensitive to the presence of vorticity.
We show that there is an infinite family of slow-roll parameters histories
which can produce the same spectrum of comoving curvature perturbations. After
expressing the slow-roll parameters in terms of the scale factor this
degeneracy can be shown to be related to the freedom in the choice of the
initial conditions for the second order differential equation relating the
coefficients of the curvature perturbation equations to the scale factor. This
freedom implies that in general there is no one-to-one correspondence between
the spectrum and higher order correlation functions, unless some special
conditions are satisfied by the slow-roll parameters.
We give some numerical example of expansion histories with the same spectrum
but different bispectra. We also compute what kind of perturbations of the
de-Sitter scale factor can produce the same spectrum but different slow-roll
parameters and higher correlation functions.
We study the dependence of surface mass density profiles, which can be directly measured by weak gravitational lensing, on the orientation of haloes with respect to the line-of-sight direction, using a suite of $N$-body simulations. We find that, when major axes of haloes are aligned with the line-of-sight direction, surface mass density profiles have higher amplitudes than those averaged over all halo orientations, over all scales from $0.1$ to $100\, \mathrm{Mpc}/h$ we studied. While the orientation dependence at small scales is ascribed to the halo triaxiality, our results indicate even stronger orientation dependence in the so-called two-halo regime, up to $100\,\mathrm{Mpc}/h$. The orientation dependence for the two-halo term is well approximated by a constant shift of the amplitude and therefore a shift in the halo bias parameter value. The halo bias from the two-halo term can be overestimated or underestimated by up to $\sim 30 \%$ depending on the viewing angle, which translates into the bias in estimated halo masses by up to a factor of two from halo bias measurements. The orientation dependence at large scales originates from the anisotropic halo-matter correlation function, which has an elliptical shape with the axis ratio of $\sim 0.55$ up to $100\,\mathrm{Mpc}/h$. We discuss potential impacts of halo orientation bias on other observables such as optically selected cluster samples and a clustering analysis of large-scale structure tracers such as quasars.
We use dense redshift surveys to explore the properties of galaxy clusters selected from the redMaPPer catalog of overdensities of red galaxies. Our new survey, HeCS-red (Hectospec Cluster Survey of red-sequence selected clusters), includes 10,589 new or remeasured redshifts from MMT/Hectospec observations of redMaPPer clusters at redshifts $z$=0.08-0.25 with large estimated richnesses (richness estimate $\lambda >64$). Our spectra confirm that each of these candidate clusters corresponds to an overdensity in redshift space. The redMaPPer photometric redshifts have a slight bias towards higher redshifts. We measure the scaling relation between velocity dispersion $\sigma_p$ and redMaPPer richness estimates $\lambda$. The observed relation shows intrinsic scatter of 24\% in velocity dispersion at fixed richness, and a range of a factor of two in measured $\sigma_p$ at fixed richness. We extend our analysis to HeCS-red-ext, a sample that includes several clusters selected by X-ray flux or SZ signal. The heterogeneous sample of 121 clusters in HeCS-red-ext shows similar intrinsic scatter, but the range of $\sigma_p$ at fixed richness increases to a factor of three. We evaluate the membership probability estimates $P_{mem}$ for individual galaxies provided by redMaPPer. The spectroscopic membership fraction is larger than $P_{mem}$ for $0.05\leq P_{mem}\leq 0.7$; conversely, it is smaller than $P_{mem}$ at $P_{mem}\geq 0.8$. We compare spectroscopic richness estimates to redMaPPer richness estimates and find good agreement on average, but a range of a factor of two in spectroscopic richness at fixed redMaPPer richness. Overall, within the high-richness and low-redshift cut of our sample, spectroscopically estimated parameters such as velocity dispersion correlate well with photometric richness estimates, although the relations contain substantial scatter.
With upcoming high-quality data from surveys such as the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) or the Dark Energy Spectroscopic Instrument (DESI), improving the theoretical modeling and gaining a deeper understanding of the effects of neutrinos and dark radiation on structure formation at small scales are necessary, to obtain robust constraints free from systematic biases. Using a novel suite of hydrodynamical simulations that incorporate dark matter, baryons, massive neutrinos, and dark radiation, we present a detailed study of their impact on Lyman-Alpha forest observables. In particular, we accurately measure the tomographic evolution of the shape and amplitude of the small-scale matter and flux power spectra and search for unique signatures along with preferred scales where a neutrino mass detection may be feasible. We then investigate the thermal state of the intergalactic medium (IGM) through the temperature-density relation. Our findings suggest that at k~5h/Mpc the suppression on the matter power spectrum induced by M_nu=0.1 eV neutrinos can reach ~4% at z~3 when compared to a massless neutrino cosmology, and ~10% if a massless sterile neutrino is included; surprisingly, we also find good agreement (~2%) with some analytic predictions. For the 1D flux power spectrum, the highest response to free-streaming effects is achieved at k~0.005 s/km when M_nu=0.1 eV; this k-limit falls in the Lyman-Alpha forest regime, making the small-scale 1D flux power spectrum an excellent probe for detecting neutrino and dark radiation imprints. Our results indicate that the IGM at z~3 provides the best sensitivity to active and sterile neutrinos.
We use a high-resolution N-body simulation to study the mass accretion history of infall dark matter halos before they are accreted by larger halos. We find that their formation time distribution is bimodal. Infall halos are dominated by a young population at high redshift and by an old population at low redshift. This bimodal distribution is found to be closely connected to the two phases in the accretion histories of halos. While members of the young population are still in the fast accretion phase at the accretion time, those of the old population have already entered the slow accretion phase at the time of accretion. This bimodal distribution is not found for normal halos, nor is it seen in halo merger trees generated with the extended Press-Schechter formalism. The infall halos at the time of accretion are, on average, younger than the other halos of similar mass identified at the same time. We discuss the implications of our findings in connection to the bimodal color distribution of observed galaxies and to the link between central and satellite galaxies.
Phenomenological functions $\Sigma$ and $\mu$, also known as $G_{\rm light}/G$ and $G_{\rm matter}/G$, are commonly used to parameterize modifications of the growth of large-scale structure in alternative theories of gravity. We study the values these functions can take in Horndeski theories, i.e. the class of scalar-tensor theories with second order equations of motion. We restrict our attention to models that are in a broad agreement with tests of gravity and the observed cosmic expansion history. In particular, we require the speed of gravity to be equal to the speed of light today, as required by the recent detection of gravitational waves and electromagnetic emission from a binary neutron star merger. We examine the correlations between the values of $\Sigma$ and $\mu$ analytically within the quasi-static approximation, and numerically, by sampling the space of allowed solutions. We confirm the conjecture made in [Pogosian:2016pwr] that $(\Sigma-1)(\mu -1) \ge 0$ in viable Horndeski theories. Along with that, we check the validity of the quasi-static approximation within different corners of Horndeski theory. Our results show that, even with the tight bound on the present day speed of gravitational waves, there is room within Horndeski theories for non-trivial signatures of modified gravity at the level of linear perturbations.
We show that simulations of magnetohydrodynamic (MHD) turbulence in the multiphase interstellar medium (ISM) yield an $E/B$ ratio for polarized emission from Galactic dust in broad agreement with recent $Planck$ measurements. In addition, the $B$-mode spectra display a scale dependence that is consistent with observations over the range of scales resolved in the simulations. The simulations present an opportunity to understand the physical origin of the $E/B$ ratio, and a starting point for more refined models of Galactic emission of use for both current and future CMB experiments.
The radio spectral index is a powerful probe for classifying cosmic radio sources and understanding the origin of the radio emission. Combining data at 147 MHz and 1.4 GHz from the TIFR GMRT Sky Survey (TGSS) and the NRAO VLA Sky Survey (NVSS), we produced a large-area radio spectral index map of ~80 per cent of the sky (Dec > -40 deg), as well as a radio spectral index catalogue containing 1,396,515 sources, of which 503,647 are not upper or lower limits. Almost every TGSS source has a detected counterpart, while this is true only for 36 per cent of NVSS sources. We released both the map and the catalogue to the astronomical community. The catalogue is analysed to discover systematic behaviours in the cosmic radio population. We find a differential spectral behaviour between faint and bright sources as well as between compact and extended sources. These trends are explained in terms of radio galaxy evolution. We also confirm earlier reports of an excess of steep-spectrum sources along the galactic plane. This corresponds to 86 compact and steep-spectrum source in excess compared to expectations. The properties of this excess are consistent with normal non-recycled pulsars, which may have been missed by pulsation searches due to larger than average scattering along the line of sight.
Studying the flow of baryons into and out of galaxies is an important part of understanding the evolution of galaxies over time. We present a detailed case study of the environment around an intervening Ly $\alpha$ absorption line system at $z_{\rm abs} = 0.633$, seen towards the quasar J0423$-$0130 ($z_{\rm QSO} = 0.915$). We detect with ALMA the $^{12}$CO(2--1), $^{12}$CO(3--2) and $1.2$~mm continuum emission from a galaxy at the redshift of the Ly $\alpha$ absorber at a projected distance of $135$ kpc. From the ALMA detections, we infer ISM conditions similar to those in low redshift Luminous Infrared Galaxies. DDT MUSE integral field unit observations reveal the optical counterpart of the $^{12}$CO emission line source and three additional emission line galaxies at the absorber redshift, which together form a galaxy group. The $^{12}$CO emission line detections originate from the most massive galaxy in this group. While we cannot exclude that we miss a fainter host, we reach a dust-uncorrected star-formation rate (SFR) limit of > $0.3 \text{M}_{\odot} \text{ yr}^{-1}$ within $100$ kpc from the sightline to the background quasar. We measure the dust-corrected SFR (ranging from $3$ to $50$ M$_{\odot}$ yr$^{-1}$), the morpho-kinematics and the metallicities of the four group galaxies to understand the relation between the group and the neutral gas probed in absorption. We find that the Ly $\alpha$ absorber traces either an outflow from the most massive galaxy or intra-group gas. This case study illustrates the power of combining ALMA and MUSE to obtain a census of the cool baryons in a bounded structure at intermediate redshift.
We explore the abundance, spatial distribution, and physical properties of the OVI, OVII, and OVIII ions of oxygen in circumgalactic and intergalactic media (the CGM, IGM, and WHIM). We use the TNG100 and TNG300 large volume cosmological magneto-hydrodynamical simulations. Modeling the ionization states of simulated oxygen, we find good agreement with observations of the low-redshift OVI column density distribution function (CDDF), and present its evolution for all three ions from z=0 to z=4. Producing mock quasar absorption line spectral surveys, we show that the IllustrisTNG simulations are fully consistent with constraints on the OVI content of the CGM from COS-Halos and other low redshift observations, producing columns as high as observed. We measure the total amount of mass and average column densities of each ion using hundreds of thousands of simulated galaxies spanning 10^11 < Mhalo/Msun < 10^15 corresponding to 10^9 < M*/Msun < 10^12 in stellar mass. The stacked radial profiles of OVI around halos of different masses are computed in 3D number density as well as 2D projected column, decomposing into the 1-halo and 2-halo terms, the latter of which begins to dominate for Milky Way mass halos in the WHIM just beyond the virial radius. Relating halo OVI to properties of the central galaxy, we find a correlation between the (g-r) color of a galaxy and the total amount of OVI in its CGM. In comparison to the COS-Halos finding, this leads to a dichotomy of columns around star-forming versus passive galaxies at fixed stellar (or halo) mass. We demonstrate that this correlation is a direct result of blackhole feedback associated with quenching, which also produces additional trends with other galaxy properties, and represents a causal consequence of galactic-scale baryonic feedback impacting the physical state of the circumgalactic medium.
SN 2017dio shows both spectral characteristics of a type-Ic supernova (SN) and signs of a hydrogen-rich circumstellar medium (CSM). Prominent, narrow emission lines of H and He are superposed on the continuum. Subsequent evolution revealed that the SN ejecta are interacting with the CSM. The initial SN Ic identification was confirmed by removing the CSM interaction component from the spectrum and comparing with known SNe Ic, and reversely, adding a CSM interaction component to the spectra of known SNe Ic and comparing them to SN 2017dio. Excellent agreement was obtained with both procedures, reinforcing the SN Ic classification. The light curve constrains the pre-interaction SN Ic peak absolute magnitude to be around $M_g = -17.6$ mag. No evidence of significant extinction is found, ruling out a brighter luminosity required by a SN Ia classification. These pieces of evidence support the view that SN 2017dio is a SN Ic, and therefore the first firm case of a SN Ic with signatures of hydrogen-rich CSM in the early spectrum. The CSM is unlikely to have been shaped by steady-state stellar winds. The mass loss of the progenitor star must have been intense, $\dot{M} \sim 0.02$ $(\epsilon_{H\alpha}/0.01)^{-1}$ $(v_\textrm{wind}/500$ km s$^{-1}$) $(v_\textrm{shock}/10 000$ km s$^{-1})^{-3}$ $M_\odot$ yr$^{-1}$, peaking at a few decades before the SN. Such a high mass loss rate might have been experienced by the progenitor through eruptions or binary stripping.
Easy Parameter Inference in Cosmology (EPIC) is another Markov Chain Monte Carlo (MCMC) sampler for Cosmology. It is implemented in Python (with support for legacy Python 2) and provides Bayesian parameter inference and model comparison based on the Bayesian evidence. The Parallel Tempering algorithm is included, which can help in the exploration of posterior distributions with two or more separated peaks. Adaptive routines for obtaining better efficiency with fine-tuned algorithms are being developed and will be available in future versions. In this user's guide, I give general instructions for installation and usage, including examples, and show how to modify the code in order to add new datasets and models.
Links to: arXiv, form interface, find, astro-ph, recent, 1712, contact, help (Access key information)
In a previous communication we showed that a joint analysis of Cosmic Microwave Background (CMB) data and the current measurement of the local expansion rate favours a model with a scale invariant spectrum (HZ) over the minimal $\Lambda$CDM scenario provided that the effective number of relativistic degrees of freedom, $N_{eff}$, is taken as a free parameter. Such a result is basically obtained due to the Hubble Space Telescope (HST) value of the Hubble constant, $H_0 = 73.24 \pm 1.74$ $\rm{km.s^{-1}.Mpc^{-1}}$ (68\% C.L.), as the CMB data alone discard the HZ+$N_{eff}$ model. Although such a model is not physically motivated by current scenarios of the early universe, observations pointing to a scale invariant spectrum may indicate that the origin of cosmic perturbations lies in an unknown physical process. Here, we extend the previous results performing a Bayesian analysis using joint CMB, HST, and Baryon Acoustic Oscillations (BAO) measurements. In order to take into account the well-known tension on the value of the fluctuation amplitude parameter, $\sigma_8$, we also consider Cluster Number counts (CN) and Weak Lensing (WL) data. We use two different samples of BAO data, which are obtained using two-point spatial (BAO 2PCF) and angular (BAO 2PACF) correlation functions. Our results show that a joint CMB+HST+BAO 2PCF analysis discards the HZP$+N_{eff}$ model with respect to the minimal $\Lambda$CDM scenario whereas the combination CMB+HST+BAO 2PACF favours the former model, even when an extended dataset with NC and WL is considered.
We constrain the dark energy equation of state parameter, $w$, using the power spectrum of the thermal Sunyaev-Zeldovich (tSZ) effect. We improve upon previous analyses by taking into account the trispectrum in the covariance matrix and marginalising over the foreground parameters, the correlated noise, the mass bias $B$ in the Planck universal pressure profile, and all the relevant cosmological parameters (i.e., not just $\Omega_{\mathrm{m}}$ and $\sigma_8$). We find that the amplitude of the tSZ power spectrum at $\ell\lesssim 10^3$ depends primarily on $F\equiv \sigma_{8}(\Omega_{{\mathrm{m}}}/B)^{3/8}h^{-1/5}$, where $B$ is related to more commonly used variable $b$ by $B=(1-b)^{-1}$. We measure this parameter with 2.5% precision, $F=0.476\pm 0.012$ (68% CL). By fixing the bias to $B=1.25$ and adding the local determination of the Hubble constant $H_0$ and the amplitude of the primordial power spectrum constrained by the Planck Cosmic Microwave Background (CMB) data, we find $w=-1.10\pm0.12$, $\sigma_{\mathrm{8}}=0.802\pm0.037$, and $\Omega_{{\mathrm{m}}}=0.265\pm0.022$ (68% CL). Our limit on $w$ is consistent with and is as tight as that from the distance-alone constraint from the CMB and $H_0$. Finally, by combining the tSZ power spectrum and the CMB data we find, in the $\Lambda$ Cold Dark Matter (CDM) model, the mass bias of $B=1.71\pm 0.17$, i.e., $1-b=0.58\pm 0.06$ (68% CL).
We use the dense HectoMAP redshift survey to explore the properties of 104 redMaPPer cluster candidates. The redMaPPer systems in HectoMAP cover the full range of richness and redshift (0.08 $< z <$ 0.60). Fifteen systems included in the Subaru/Hyper Suprime-Cam public data release are bona fide clusters. The median number of spectroscopic members per cluster is $\sim20$. We include redshifts of 3547 member candidates listed in the redMaPPer catalog whether they are cluster members or not. We evaluate the redMaPPer membership probability spectroscopically. The scaled richness ({\lambda}rich/S) provided by redMaPPer correlates tightly with the spectroscopic richness regardless of the cluster redshift and appears to be a better mass proxy than the original richness, {\lambda}rich. The purity (number of real systems) in redMaPPer exceeds 90% even at the lowest richness; however, there is some incompleteness. Five massive galaxy clusters (M $\gtrsim 2 \times 10$^{13}$ M$_{\odot}$) associated with X-ray emission in the HectoMAP region are missing from the catalog.
The gravitational wave (GW) as standard siren directly determines the luminosity distance from the gravitational waveform without reference to specific cosmological model, of which the redshift can be obtained separately by means of the electromagnetic counterpart like GW events from binary neutron star (BNS) and massive black hole binaries (MBHBs). To see to what extent the standard siren can reproduce the presumed dipole anisotropy written in the simulated data of standard siren events from typical configurations of GW detectors, we find that, (1) for Laser Interferometer Space Antenna (LISA) with different MBHB models during five-year observations, the cosmic isotropy can be ruled out at $3\sigma$ confidence level (CL) and the dipole direction can be constrained roughly around $20\%$ at $2\sigma$ CL, as long as the dipole amplitude is larger than $0.06$, $0.05$ and $0.03$ for MBHB models Q3d, pop III and Q3nod with increasing constraining ability, respectively; (2) for Einstein Telescope (ET) with no-less than $200$ standard siren events, the cosmic isotropy can be ruled out at $3\sigma$ CL if the dipole amplitude is larger than $0.06$, and the dipole direction can be constrained within $20\%$ at $3\sigma$ CL if the dipole amplitude is near $0.1$; (3) for Deci-Hertz Interferometer Gravitational wave Observatory (DECIGO) with no-less than $100$ standard siren events, the cosmic isotropy can be ruled out at $3\sigma$ CL for dipole amplitude larger than $0.03$ , and the dipole direction can even be constrained within $10\%$ at $3\sigma$ CL if the dipole amplitude is larger than $0.07$. Our work manifests the promising perspective of the constraint ability on the cosmic anisotropy from standard siren approach.
The Baryon Acoustic Oscillations (BAO) refer to the ripples of material density in the Universe. As the most direct density tracers in the universe, galaxies have been commonly used in studies of BAO peak detection. The spatial number density of galaxies, to a certain extent, reflects the distribution of the material density of our Universe. Using galaxies as matter tracers, we can construct more overlapping empty spheres (DT voids) than the matter tracers, via Delaunay Triangulation technique. We show that their radii excellently reflect the galaxy number density round them, and they can serve as reliable different density region tracers. Using the data from an unprecedented large-scale $N$-body simulation "TianNu", we conduct some fundamental statistical studies and clustering analysis of the DT voids. We discuss in detail the representative features of two-point correlation functions of different DT void populations. We show that the peak, the position of which corresponds to the average radius of data samples, is the most representative feature of the two-point correlation function of the DT voids. In addition, we also construct another voids, the disjoint voids, and investigate their some statistical properties and clustering properties. And we find that the occupied space of all disjoint voids accounts for about $45\%$ of the volume of the simulation box, regardless of the number density of mock galaxies. We also investigate the BAO detections based on different tracers, i.e. mock galaxies, low-density region tracers, and high-density region tracers respectively. Our results show that BAO intensities detected by low/high-density region tracers are enhanced significantly compared to the BAO detection by mock galaxies, for the mock galaxy catalogue with the number density of $7.52\times10^{-5}$ $h^3$ Mpc$^{-3}$.
This paper aims to put constraints on the transition redshift $z_t$, which determines the onset of cosmic acceleration, in cosmological-model independent frameworks. In order to do that, we consider a flat universe and suppose a parametrization for the comoving distance $D_C(z)$ up to third degree on $z$, a second degree parametrization for the Hubble parameter $H(z)$ and a linear parametrization for the deceleration parameter $q(z)$. For each case, we show that the type supernovae Ia and $H(z)$ data complement each other on the parameter spaces and tighter constrains for the transition redshift are obtained. By combing the supernovae type Ia observations and Hubble parameter measurements it is possible to constrain the values of $z_t$ as $0.806\pm 0.094$, $0.870\pm 0.063$ and $0.973\pm 0.058$ at 1$\sigma$ c.l., for each framework, respectively. Such approaches provide reasonably model-independent estimates of this cosmological parameter.
The standard model of cosmology predicts the existence of cosmic neutrino background in the present Universe. Were we plan to detect relic neutrinos in the vicinity of the Earth, as proposed in the upcoming PTOLEMY experiment, it would be necessary to evaluate the gravitational clustering effects on cosmic relic neutrinos in the Milky Way. In this work we introduce a reweighting technique in the $N$-one-body simulation method, so that a single run of simulation can yield neutrino density profiles for a variety of neutrino masses and different phase space distributions. In light of current experimental results that favor small neutrino masses ($\lesssim 0.1~\mathrm{eV}$), the number density contrast of relic neutrinos around the Earth is found to be almost inversely proportional to the square of neutrino mass. The density contrast--mass relation and the reweighting technique simplify the investigation on the gravitational clustering of relic neutrinos, and thus are useful for the future detection of cosmic neutrino background.
We explore the phenomenology of having a second epoch of dark matter annihilation into dark radiation long after the standard thermal freeze-out. Such a hidden \emph{reannihilation} process could affect visible sectors only gravitationally. As a concrete realization we consider self-interacting dark matter (SIDM) with a light force mediator coupled to dark radiation. We demonstrate how resonantly Sommerfeld enhanced cross sections emerge to induce the reannihilation epoch. The effect is a local modification of the Hubble expansion rate and we show that the Cosmic Microwave Background (CMB) measurements -- as well as other observations -- have a high sensitivity to observe this phenomenon. Special attention is given to the model region where late kinetic decoupling and strong self-interactions can alleviate several small-scale problems in the cold dark matter paradigm at the same time. Interestingly, we find that reannihilation might here also simultaneously lower the tension between CMB and low-redshift astronomical observations of $H_0$ and $\sigma_8$. Moreover, we identify reannihilation as a clear signature to discriminate between the phenomenologically otherwise almost identical vector and scalar mediator realizations of SIDM.
Effective field theory of dark matter fluid on large scales predicts the presence of viscosity of the order of $10^{-6} H_0 M_P^2$. It has been shown that this magnitude of viscosities can resolve the discordance between large scale structure observations and Planck CMB data in the $\sigma_8$-$\Omega_m^0$ and $H_0$-$\Omega_m^0$ parameters space. Massive neutrinos suppresses the matter power spectrum on the small length scales similar to the viscosities. We show that by including the effective viscosity, which arises from summing over non linear perturbations at small length scales, severely constrains the cosmological bound on neutrino masses. Under a joint analysis of Planck CMB and different large scale observation data, we find that upper bound on the sum of the neutrino masses at 2-$\sigma$ level, decreases from $\sum m_\nu \le 0.396\,$eV (normal hierarchy) and $\sum m_\nu \le 0.378 \,$eV (inverted hierarchy) to $\sum m_\nu \le 0.267\,$eV (normal hierarchy) and $\sum m_\nu \le 0.146\,$eV (inverted hierarchy) when the effective viscosities are included.
Could there be a large population of intermediate mass black holes (IMBHs) formed in the early universe? Whether primordial or formed in Population III, these are likely to be very subdominant compared to the dark matter density, but could seed early dwarf galaxy/globular cluster and supermassive black hole formation. Via survival of dark matter density spikes, we show here that a centrally concentrated relic population of IMBHs, along with ambient dark matter, could account for the Fermi gamma-ray "excess" in the Galactic center because of dark matter particle annihilations.
We present Self-Destructing Dark Matter (SDDM), a new class of dark matter models which are detectable in large neutrino detectors. In this class of models, a component of dark matter can transition from a long-lived state to a short-lived one by scattering off of a nucleus or an electron in the Earth. The short-lived state then decays to Standard Model particles, generating a dark matter signal with a visible energy of order the dark matter mass rather than just its recoil. This leads to striking signals in large detectors with high energy thresholds. We present a few examples of models which exhibit self destruction, all inspired by bound state dynamics in the Standard Model. The models under consideration exhibit a rich phenomenology, possibly featuring events with one, two, or even three lepton pairs, each with a fixed invariant mass and a fixed energy, as well as non-trivial directional distributions. This motivates dedicated searches for dark matter in large underground detectors such as Super-K, Borexino, SNO+, and DUNE.
The existence of a light or massive scalar field with a coupling to matter weaker than gravitational strength is a possible source of violation of the weak equivalence principle. We use the first results on the E\"otv\"os parameter by the MICROSCOPE experiment to set new constraints on such scalar fields. For a massive scalar field of mass smaller than $10^{-12}$~eV (i.e. range larger than a few $10^5$~m) we improve existing constraints by one order of magnitude to $|\alpha|<10^{-11}$ if the scalar field couples to the baryon number and to $|\alpha|<10^{-12}$ if the scalar field couples to the difference between the baryon and the lepton numbers. We also consider a model describing the coupling of a generic dilaton to the standard matter fields with five parameters, for a light field: we find that for masses smaller than $10^{-12}$eV, the constraints on the dilaton coupling parameters are improved by one order of magnitude compared to previous equivalence principle tests.
In arXiv:1705.10172 it was proposed that string theory replaces Schwarzschild black holes with horizonless thin shells with an AdS interior. In this paper we extend the analysis to slowly rotating black holes, solving the Israel-Lanczos-Sen junction conditions for a rotating shell composed of stringy matter to determine the metric. Outside of the shell we find a vacuum solution that differs from Kerr with a 39% larger quadrupole moment. We discuss the observational consequences and explore the possibility to distinguish between a black shell and a black hole. Promising methods include imaging of the black hole at the center of the Milky Way using the Event Horizon Telescope, precision measurements of stars in close orbits around the central black hole, and future observations of colliding super massive black holes using the space based gravitational wave observatory LISA.
The recent discovery of gravitational wave events has offered us unique testbeds of gravity in the strong and dynamical field regime. One possible modification to General Relativity is the gravitational parity violation that gives rise to the amplitude birefringence in gravitational waves, where one of the circularly-polarized mode is amplified while the other one is suppressed during their propagation. In this paper, we study how well one can measure gravitational parity violation via the amplitude birefringence of gravitational waves from stellar-mass black hole binaries. We choose Chern-Simons gravity as an example and work within an effective field theory formalism. We consider gravitational waves from both individual sources and stochastic gravitational wave backgrounds. Regarding bounds from individual sources, we estimate them using a Fisher analysis and carry out Monte Carlo simulations by randomly distributing sources over their sky location and binary orientation. We find that the bounds on the scalar field evolution in Chern-Simons gravity from GW150914 are too weak to satisfy the weak Chern-Simons approximation, while aLIGO with its design sensitivity can place meaningful bounds. Regarding bounds from stochastic gravitational wave backgrounds, we set the threshold signal-to-noise ratio for detection of the parity-violation mode as 5 and estimate projected bounds with future detectors assuming that signals are consistent with no parity violation. We find that a network of two third-generation detectors is able to place bounds that are slightly stronger than current binary pulsar bounds. Since gravitational wave observations probe either the difference in parity violation between the source and the detector or the line-of-sight integration of the scalar field, such bounds are complementary to local measurements from solar system experiments and binary pulsar observations.
Recently a new inflationary scenario was proposed in arXiv:1703.09020 which can be applicable to an inflaton having multiple vacua. In this letter, we consider a more general situation where the inflaton potential has a (UV) saddle point around the Planck scale. This class of models can be regarded as a natural generalization of the hillclimbing Higgs inflation (arXiv:1705.03696).
We investigate cosmological dynamics based on $f(R)$ gravity in the Palatini formulation. In this study we use the dynamical system methods. We show that the evolution of the Friedmann equation reduces to the form of the piece-wise smooth dynamical system. This system is is reduced to a 2D dynamical system of the Newtonian type. We demonstrate how the trajectories can be sewn to guarantee $C^0$ extendibility of the metric similarly as `Milne-like' FLRW spacetimes are $C^0$-extendible. We point out that importance of dynamical system of Newtonian type with non-smooth right-hand sides in the context of Palatini cosmology. In this framework we can investigate singularities which appear in the past and future of the cosmic evolution. We consider cosmological systems in both Einstein and Jordan frames. We show that at each frame the topological structures of phase space are different.
The DArk Matter Particle Explorer (DAMPE) has reported a measurement of the flux of high energy cosmic ray electrons plus positrons (CREs) in the energy range between $25$ GeV and $4.6$ TeV. With unprecedented high energy resolution, the DAMPE data exhibit an excess of the CREs flux at an energy of around $1.4$ TeV. In this letter, we discuss how the observed excess can be understood in a minimal framework where the Standard Model (SM) is supplemented by a stable SM singlet scalar as dark matter (DM) and type II seesaw for generating the neutrino mass matrix. In our framework, a pair of DM particles annihilates into a pair of the SM SU(2) triplet scalars ($\Delta$s) in type II seesaw, and the subsequent $\Delta$ decays create the primary source of the excessive CREs around $1.4$ TeV. The lepton flavor structure of the primary source of CREs has a direct relationship with the neutrino oscillation data. We find that the DM interpretation of the DAMPE excess determines the pattern of neutrino mass spectrum to be the inverted hierarchy type, taking into account the constraints from the Fermi-LAT observations of dwarf spheroidal galaxies.
We consider a cosmological setup with the inflaton field in the presence of a redshift dependent Lorentz-violating time-like background to address the inflationary regime and other phases of the Universe. We also show that the regime of dark energy at large distances (low redshifts) is essentially dominated by the presence of the Lorentz-violating background.
Cosmic strings are topological defects which can be formed in GUT-scale phase transitions in the early universe. They are also predicted to form in the context of string theory. The main mechanism for a network of Nambu-Goto cosmic strings to lose energy is through the production of loops and the subsequent emission of gravitational waves, thus offering an experimental signature for the existence of cosmic strings. Here we report on the analysis conducted to specifically search for gravitational-wave bursts from cosmic string loops in the data of Advanced LIGO 2015-2016 observing run (O1). No evidence of such signals was found in the data, and as a result we set upper limits on the cosmic string parameters for three recent loop distribution models. In this paper, we initially derive constraints on the string tension $G\mu$ and the intercommutation probability, using not only the burst analysis performed on the O1 data set, but also results from the previously published LIGO stochastic O1 analysis, pulsar timing arrays, cosmic microwave background and Big-Bang nucleosynthesis experiments. We show that these data sets are complementary in that they probe gravitational waves produced by cosmic string loops during very different epochs. Finally, we show that the data sets exclude large parts of the parameter space of the three loop distribution models we consider.
According to the Weak Equivalence Principle, all bodies should fall at the same rate in a gravitational field. The MICROSCOPE satellite, launched in April 2016, aims to test its validity at the $10^{-15}$ precision level, by measuring the force required to maintain two test masses (of titanium and platinum alloys) exactly in the same orbit. A non-vanishing result would correspond to a violation of the Equivalence Principle, or to the discovery of a new long-range force. Analysis of the first data gives $\delta\rm{(Ti,Pt)}= [-1 \pm 9 (\mathrm{stat}) \pm 9 (\mathrm{syst})] \times 10^{-15}$ (1$\sigma$ statistical uncertainty) for the titanium-platinum E\"otv\"os parameter characterizing the relative difference in their free-fall accelerations.
Links to: arXiv, form interface, find, astro-ph, recent, 1712, contact, help (Access key information)
We analyze the sources of free electrons that produce the large dispersion measures, DM $\approx 300-1600$ (in units cm$^{-3}$ pc), observed toward fast radio bursts (FRBs). Individual galaxies typically produce DM $\sim 25-60$ cm$^{-3}$ pc from ionized gas in their disk, disk-halo interface, and circumgalactic medium. Toward an FRB source at redshift $z$, a homogeneous IGM containing a fraction $f_{\rm IGM}$ of cosmological baryons will produce DM $= (935~{\rm cm}^{-3}~{\rm pc}) f_{\rm IGM} \, h_{70}^{-1} I(z)$, where $I(z) = (2/3 \Omega_m)[ \{ \Omega_m(1+z)^3 + \Omega_{\Lambda} \}^{1/2} - 1 ]$. A structured IGM of photoionized Ly-alpha absorbers in the cosmic web produces similar dispersion, modeled from the observed distribution, $f_b(N,z)$, of H I (Lya-forest) absorbers in column density and redshift with ionization corrections and scaling relations from cosmological simulations. An analytic formula for DM($z$) applied to observed FRB dispersions suggests that $z_{\rm FRB} \approx 0.2-1.5$ for an IGM containing a significant baryon fraction, $f_{\rm IGM} = 0.6\pm0.1$. Future surveys of the statistical distribution, DM($z)$, of FRBs identified with specific galaxies and redshifts, can be used to calibrate the IGM baryon fraction and distribution of Ly-alpha absorbers. Fluctuations in DM at the level $\pm10$ cm$^{-3}$ pc will arise from filaments and voids in the cosmic web.
The evolution of neutral hydrogen (HI) across redshifts is a powerful probe of cosmology, large scale structure in the universe and the intergalactic medium. Using a data-driven halo model to describe the distribution of HI in the post-reionization universe ($z \sim $ 5 to 0), we obtain the best-fitting parameters from a rich sample of observational data: low redshift 21-cm emission line studies, intermediate redshift intensity mapping experiments, and higher redshift Damped Lyman Alpha (DLA) observations. Our model describes the abundance and clustering of neutral hydrogen across redshifts 0 - 5, and is useful for investigating different aspects of galaxy evolution and for comparison with hydrodynamical simulations. The framework can be applied for forecasting future observations with neutral hydrogen, and extended to the case of intensity mapping with molecular and other line transitions at intermediate redshifts.
It is broadly accepted that Supermassive Black Holes (SMBHs) are located in the centers of most massive galaxies, although there is still no convincing scenario for the origin of their massive seeds. It has been suggested that primordial black holes (PBHs) of masses $\gtrsim 10^{2} M_\odot$ may provide such seeds, which would grow to become SMBHs. We suggest an observational test to constrain this hypothesis: gas accretion around PBHs during the cosmic dark ages powers the emission of high energy photons which would modify the spin temperature as measured by 21cm Intensity Mapping (IM) observations. We model and compute their contribution to the standard sky-averaged signal and power spectrum of 21cm IM, accounting for its substructure and angular dependence for the first time. If PBHs exist, the sky-averaged 21cm IM signal in absorption would be higher, while we expect an increase in the power spectrum for $\ell~\gtrsim 10^2-10^3$. We also forecast PBH detectability and measurement errors in the abundance and Eddington ratios for different fiducial parameter configurations for various future experiments, ranging from SKA to a futuristic radio array on the dark side of the Moon. While the SKA could provide a detection, only a more ambitious experiment would provide accurate measurements.
The observed tightness of the mass discrepancy-acceleration relation (MDAR) poses a fine-tuning challenge to current models of galaxy formation. We propose that this relation could arise from collisional interactions between baryons and dark matter (DM) particles, without the need for modification of gravity or ad hoc feedback processes. We assume that these interactions satisfy the following three conditions: (i) the relaxation time of DM particles is comparable to the dynamical time in disk galaxies; (ii) DM exchanges energy with baryons due to elastic collisions; (iii) the product between the baryon-DM cross section and the typical energy exchanged in a collision is inversely proportional to the DM number density. We present an example of a particle physics model that gives a DM-baryon cross section with the desired density and velocity dependence. Direct detection constraints require our DM particles to be either very light ($m << m_b$) or very heavy ($m >> m_b$), corresponding respectively to heating and cooling of DM by baryons. In both cases, our mechanism applies and an equilibrium configuration can in principle be reached. Here, we focus on the heavy DM/cooling case as it is technically simpler. Under these assumptions, we find that rotationally-supported disk galaxies could naturally settle to equilibrium configurations satisfying a MDAR at all radii without invoking finely tuned feedback processes. We also discuss issues related to the small scale clumpiness of baryons, as well as predictions for pressure-supported systems. We argue in particular that galaxy clusters do not follow the MDAR despite being DM-dominated because they have not reached their equilibrium configuration. Finally, we revisit existing phenomenological, astrophysical and cosmological constraints on baryon-DM interactions in light of the unusual density dependence of the cross section.
Utilizing gravitational-wave (GW) lensing opens a new way to understand the history and structure of the universe. In spite of coarse angular resolution and short duration of observation, we show that LIGO can detect the GW lensing induced by small structures, in particular by compact dark matter (DM) of $10 - 10^5 M_{\rm sun}$, which remains an interesting DM candidate. The lensing is detected through GW frequency chirping, creating the natural and rapid change of lensing patterns: frequency-dependent amplification and modulation of waveforms. As a highest-frequency GW detector, LIGO is a unique GW lab to probe such light compact DM. With design sensitivity of Advanced LIGO, one-year observation can detect as many as 1000 lensed GWs and constrain compact DM fraction as small as $f_{\rm DM} \gtrsim 10^{-3}$.
Aiming at exploring the nature of dark energy (DE), we use forty-three observational Hubble parameter data (OHD) in the redshift range $0 \leqslant z \leqslant 2.36$ to make a cosmological model-independent test of the two-point $Omh^2(z_{2};z_{1})$ diagnostic. In $\Lambda$CDM model, with equation of state (EoS) $w=-1$, $Omh^2 \equiv \Omega_m h^2$ holds to be tenable, where $\Omega_m$ is the matter density parameter, and $h$ is Hubble parameter at present. Since the direct exploitation of OHD to the $Omh^2(z_{2};z_{1})$ diagnostic generate an obscure images, which is bad for analysis, we utilize two methods: the weighted mean and median statistics to bin the data categorically. The binning methods turn out to be promising and considered to be robust. By applying the $Omh^2(z_{2};z_{1})$ diagnostic into binned data, we find that the values of $Omh^2$ fluctuate as the consistent redshift intervals change, i.e., not being constant, and especially significant for the weighted mean case. Therefore, we conclude that the $\Lambda$CDM model is in doubt and other dynamical DE models with an evolving EoS should be considered more likely to be the candidates that interpret the acceleration expansion of the universe.
We investigate a compelling model of quintessential inflation in the context of $\alpha$-attractors, which naturally result in a scalar potential featuring two flat regions, the inflationary plateau and the quintessential tail. The "asymptotic freedom" of $\alpha$-attractors, near the kinetic poles, suppresses radiative corrections and interactions, which would otherwise threaten to lift the flatness of the quintessential tail and cause a 5th-force problem respectively. Since this is a non-oscillatory inflation model, we reheat the Universe through instant preheating. The parameter space is constrained by both inflation and dark energy requirements. We find an excellent correlation between the inflationary observables and model predictions, in agreement with the $\alpha$-attractors set-up. We also obtain successful quintessence for natural values of the parameters. Our model predicts potentially sizeable tensor perturbations (at the level of 1%) and a slightly varying equation of state for dark energy, to be probed in the near future.
In this paper we develop a formalism for studying the nonrelativistic limit of relativistic field theories in a systematic way. By introducing a simple, nonlocal field redefinition, we may transform the description of a given relativistic theory in terms of a real-valued, self-interacting scalar field into an equivalent theory, describing a complex scalar field. Our low-energy effective theory incorporates corrections from nonvanishing momentum as well as from the backreaction of fast-oscillating terms on the behavior of the dominant, slowly varying component of the field. Possible applications of our new approach include axion dark matter, though the methods developed here should be applicable to the low-energy limits of other field theories as well.
An exact solution for the spatially flat scale-invariant Cosmology, recently proposed by Maeder (2017) is deduced. No deviation from the numerical solution was detected. The exact solution yields transparency for the dynamical equations and faster cosmological constraints may be performed.
We describe a major update to the public GIZMO code. GIZMO has been used in simulations of cosmology; galaxy and star formation and evolution; black hole accretion and feedback; proto-stellar disk dynamics and planet formation; fluid dynamics and plasma physics; dust-gas dynamics; giant impacts and solid-body interactions; collisionless gravitational dynamics; and more. This release of the public code supports: hydrodynamics (using various mesh-free finite-volume Godunov methods or SPH), ideal and non-ideal MHD, anisotropic conduction and viscosity, radiative cooling and chemistry, star and black hole formation and feedback, sink particles, dust-gas (aero)-dynamics (with or without magnetic fields), elastic/plastic dynamics, arbitrary (gas, stellar, degenerate, solid/liquid material) equations of state, passive scalar/turbulent diffusion, large-eddy and shearing boxes, self-gravity with fully-adaptive force softenings, arbitrary cosmological expansion, and on-the-fly group-finding. It is massively-parallel with hybrid MPI+OpenMP scaling verified up to >1 million threads. The code is extensively documented, with test problems and tutorials provided for these different physics modules.
The Weak Gravity Conjecture (WGC) was proposed to constrain Effective Field Theories (EFTs) with Abelian gauge symmetry coupled to gravity. In this article, I study the WGC from low energy observers' perspective, and revisit the issue of to what extent the WGC actually constrains EFTs. For this purpose, for a given EFT, I introduce associated idealized low energy observers who only have access to the energy scale below the UV cut-off scale of the EFT. In the framework of EFT, there is a clear difference between the particles lighter than the UV cut-off scale and the particles which are heavier than the UV cut-off scale, as the lighter particles can be created below the UV cut-off scale while the heavier particles are not. This difference implies that the knowledge of the low energy observers on the stable heavy particles can be limited, as the availability of the stable heavy particles is determined by the environment prepared by some UV theory unknown to the low energy observers. The limitation of the knowledge of the low energy observers regarding the stable heavy particles whose mass is above the UV cut-off scale of the EFT leads to the limitation of the WGC for constraining EFTs. To illustrate these points in an example, I analyze a model proposed by Saraswat arXiv:1608.06951 which respects the WGC at high energy, but which may appear to violate the WGC for the low energy observers. Implications of the analysis to the bottom-up model buildings using EFTs are discussed.
The tentative 1.4 TeV excess in the $e^+e^-$ spectrum measured by The DArk Matter Particle Explorer (DAMPE) motivates the possible existence of one or more local dark matter concentrated regions. In particular, Ultra-compact Micro Halos (UCMHs) seeded by large density perturbations in the early universe, allocated within ~0.3 kpc from the solar system, could provide the potential source of electrons and positrons produced from dark matter annihilation, enough to explain the DAMPE signal. Here we consider a UCMH with density profile assuming radial in-fall and explore the preferred halo parameters to explain the 1.4 TeV "DAMPE excess". We find that typical parameter space of UCMHs can easily explain the "DAMPE excess" with usual thermal-averaged annihilation cross section of WIMP. The fraction of dark matter stored in such UCMHs in the Galactic-scale halo can be reduced to as small as $O(10^{-5})$, well within the current cosmological and astrophysical constraints.
Links to: arXiv, form interface, find, astro-ph, recent, 1712, contact, help (Access key information)
We investigate the cosmological information contained in the cross-correlation between the Integrated Sachs-Wolfe (ISW) of the Cosmic Microwave Background (CMB) anisotropy pattern and galaxy clusters from future wide surveys. Future surveys will provide cluster catalogues with a number of objects comparable with galaxy catalogues currently used for the detection of the ISW signal by cross-correlation with the CMB anisotropy pattern. By computing the angular power spectra of clusters and the corresponding cross-correlation with CMB, we perform a signal-to-noise ratio (SNR) analysis for the ISW detection as expected from the eROSITA and the Euclid space missions. We discuss the dependence of the SNR of the ISW-cluster cross-correlation on the specifications of the catalogues and on the reference cosmology. We forecast that the SNRs for ISW-cluster cross-correlation are alightly smaller compared to those which can be obtained from future galaxy surveys but the signal is expected to be detected at high significance, i.e. more than $> 3\,\sigma$. We also forecast the joint constraints on parameters of model extensions of the concordance $\Lambda$CDM cosmology by combining CMB and the ISW-cluster cross-correlation.
Recent studies have presented evidence for tension between the constraints on Omega_m and sigma_8 from the cosmic microwave background (CMB) and measurements of large-scale structure (LSS). This tension can potentially be resolved by appealing to extensions of the standard model of cosmology and/or untreated systematic errors in the modelling of LSS, of which baryonic physics has been frequently suggested. We revisit this tension using, for the first time, carefully-calibrated cosmological hydrodynamical simulations, which thus capture the back reaction of the baryons on the total matter distribution. We have extended the BAHAMAS simulations to include a treatment of massive neutrinos, which currently represents the best motivated extension to the standard model. We make synthetic thermal Sunyaev-Zel'dovich effect, weak galaxy lensing, and CMB lensing maps and compare to observed auto- and cross-power spectra from a wide range of recent observational surveys. We conclude that: i) in general there is tension between the primary CMB and LSS when adopting the standard model with minimal neutrino mass; ii) after calibrating feedback processes to match the gas fractions of clusters, the remaining uncertainties in the baryonic physics modelling are insufficient to reconcile this tension; and iii) invoking a non-minimal neutrino mass, typically of 0.2-0.4 eV (depending on the priors on the other relevant cosmological parameters and the datasets being modelled), can resolve the tension. This solution is fully consistent with separate constraints on the summed neutrino mass from the primary CMB and baryon acoustic oscillations, given the internal tensions in the Planck primary CMB dataset.
The discovery of the electromagnetic counterpart to GW170817 severely constrains the tensor mode propagation speed, eliminating a large model space of Horndeski theory. We use the cosmic microwave background data from Planck and the joint analysis of the BICEP2/Keck Array and Planck, galaxy clustering data from the SDSS LRG survey, BOSS baryon acoustic oscillation data, and redshift space distortion measurements to place constraints on the remaining Horndeski parameters. We evolve the Horndeski parameters as power laws with both the amplitude and power law index free. We find a 95% CL upper bound on the present-day coefficient of the Hubble friction term in the cosmological propagation of gravitational waves is 2.38, whereas General Relativity gives 2 at all times. While an enhanced friction suppresses the amplitude of the reionization bump of the primordial B-mode power spectrum at $\ell < 10$, our result limits the suppression to be less than 0.8%. This constraint is primarily due to the scalar integrated Sachs-Wolfe effect in temperature fluctuations at low multipoles.
We cross-correlate the positions of damped Lyman-$\alpha$ systems (DLAs) and their parent quasar catalog with a convergence map derived from the Planck cosmic microwave background (CMB) temperature data. We make consistent detections of the lensing signal of both samples in both Fourier and configuration space. By interpreting the excess signal present in the DLA catalog with respect to the parent quasar catalog as caused by the large scale structure traced by DLAs, we are able to infer the bias of these objects: $b_{\rm DLA}=2.6\pm0.9$. These results are consistent with previous measurements made in cross-correlation with the Lyman-$\alpha$ forest, although the current noise in the lensing data and the low number density of DLAs limits the constraining power of this measurement. We discuss the robustness of the analysis with respect to a number different systematic effects and forecast prospects of carrying out this measurement with data from future experiments.
We investigate the nonlinear evolution of structure in variants of the standard cosmological model which display damped density fluctuations relative to cold dark matter (e.g. in which cold dark matter is replaced by warm or interacting DM). Using N-body simulations, we address the question of how much information is retained from the initial linear power spectrum following the nonlinear growth of structure. We run a suite of N-body simulations with different initial linear matter power spectra to show that, once the system undergoes nonlinear evolution, the shape of the linear power spectrum at high wavenumbers does not affect the non-linear power spectrum, while it still matters for the halo mass function. Indeed, we find that linear power spectra which differ from one another only at wavenumbers larger than their half-mode wavenumber give rise to (almost) identical nonlinear power spectra at late times, regardless of the fact that they originate from different models with damped fluctuations. On the other hand, the halo mass function is more sensitive to the form of the linear power spectrum. Exploiting this result, we propose a two parameter model of the transfer function in generic damped scenarios, and show that this parametrisation works as well as the standard three parameter models for the scales on which the linear spectrum is relevant.
We develop an approach to parametrize cosmological perturbations beyond linear order for general dark energy and modified gravity models characterized by a single scalar degree of freedom. We derive the full nonlinear action, focusing on Horndeski theories. In the quasi- static, non-relativistic limit, there are a total of six independent relevant operators, three of which start at nonlinear order. The new nonlinear couplings modify, beyond linear order, the generalized Poisson equation relating the Newtonian potential to the matter density contrast. We derive this equation up to cubic order in perturbations and, in a companion article, we apply it to compute the one-loop matter power spectrum. Within this approach, we also discuss the Vainshtein regime around spherical sources and the relation between the Vainshtein scale and the nonlinear scale for structure formation.
We develop an approach to compute observables beyond the linear regime of dark matter perturbations for general dark energy and modified gravity models. We do so by combining the Effective Field Theory of Dark Energy and Effective Field Theory of Large-Scale Structure approaches. In particular, we parametrize the linear and nonlinear effects of dark energy on dark matter clustering in terms of the Lagrangian terms introduced in a companion paper, focusing on Horndeski theories and assuming the quasi-static approximation. The Euler equation for dark matter is sourced, via the Newtonian potential, by new nonlinear vertices due to modified gravity and, as in the pure dark matter case, by the effects of short-scale physics in the form of the divergence of an effective stress tensor. The effective fluid introduces a counterterm in the solution to the matter continuity and Euler equations, which allows a controlled expansion of clustering statistics on mildly nonlinear scales. We use this setup to compute the one-loop dark-matter power spectrum.
In the framework of the minimal quartet-metric gravity/systogravity, a scalar graviton/systolon is stated as a universal dark component, with supplementary manifestations in the different contexts either as dark matter or dark energy. An ensuing extension to the standard {\Lambda}CDM model is developed. A modification of the late expansion of the Universe, with an attractor of a scalar master equation defining an effective cosmological constant, which supersedes the true one, is proposed. A new partial solution to the cosmological constant problem is discussed.
Analog condensed matter systems present an exciting opportunity to simulate early Universe models in table-top experiments. We consider a recent proposal for an analog condensed matter experiment to simulate the relativistic quantum decay of the false vacuum. In the proposed experiment, two ultra-cold condensates are coupled via a time-varying radio-frequency field. The relative phase of the two condensates in this system is approximately described by a relativistic scalar field with a potential possessing a series of false and true vacuum local minima. If the system is set up in a false vacuum, it would then decay to a true vacuum via quantum mechanical tunnelling. Should such an experiment be realized, it would be possible to answer a number of open questions regarding non-perturbative phenomena in quantum field theory and early Universe cosmology. In this paper, we illustrate a possible obstruction: the time-varying coupling that is invoked to create a false vacuum for the long-wavelength modes of the condensate leads to a destabilization of shorter wavelength modes within the system via parametric resonance. We focus on an idealized setup in which the two condensates have identical properties and identical background densities. Describing the system by the coupled Gross-Pitaevskii equations (GPE), we use the machinery of Floquet theory to perform a linear stability analysis, calculating the wavenumber associated with the first instability band for a variety of experimental parameters. However, we demonstrate that, by tuning the frequency of the time-varying coupling, it may be possible to push the first instability band outside the validity of the GPE, where dissipative effects are expected to damp any instabilities. This provides a viable range of experimental parameters to perform analog experiments of false vacuum decay.
Very long baseline interferometry (VLBI) at millimeter (mm) wavelengths is being employed to resolve event-horizon scale structure of the environment surrounding the Milky-Way black hole, at an angular resolution of a few tens of micro-arcseconds. The same approach could also resolve the orbital separation of a population of massive black hole binaries (MBHBs). Modeling the inspiral of binaries due to gravitational wave emission and gas and requiring binary orbital periods of less than 10 years, we estimate that there may exist ~100 resolvable MBHBs that are bright enough to be observed by mm-wavelength VLBI instruments over the entire sky, at redshifts z<0.5. We propose to search for these resolvable MBHBs by identifying binaries with the required orbital separations from periodic quasar light curves identified in optical and near-IR surveys. These periodic-light-curve candidates can be followed up with radio observations to determine their promise for observation with VLBI at mm wavelengths. VLBI observations over the timescale of a binary orbit can allow unprecedented precision in the measurement of the binary mass, to within 30%. In combination with an independent binary mass measurement, VLBI observation would allow a novel, of order 10%, measurement of the Hubble constant, independent from those currently proposed and employed.
In this work we clarify aspects of renormalization on curved backgrounds focussing on the potential ramifications on the amplitude of inflationary perturbations. We provide an alternate view of the often used adiabatic prescription by deriving a correspondence between the adiabatic subtraction terms and traditional renormalization. Specifically, we show how adiabatic subtraction can be expressed as a set of counter terms that are introduced by redefining the bare parameters of the action. Our novel representation of adiabatic subtraction then allows us to easily find other renormalization prescriptions differing only in the finite parts of the counter terms. As our main result, we present for quadratic inflation how one may consistently express the renormalization of the spectrum of perturbations from inflation as a redefinition of the bare cosmological constant and Planck mass such that the observable predictions coincide with the unrenormalized result.
We propose to interpret the DAMPE electron excess at 1.5 TeV through scalar (or Dirac fermion) dark matter(DM) annihilation with doubly charged scalar mediators that have lepton-specific Yukawa couplings. Hierarchy of such lepton-specific Yukawa couplings are generated through Froggatt-Nielsen mechanism, so that the dark matter annihilation products are dominantly electrons (or $e,\mu$, $e,\tau$). Stringent constraints from LEP2 on intermediate vector boson production will no longer hold in our scenarios. In the case of scalar DM, we discuss one scenario with DM annihilating directly to leptons and the other one with DM annihilating to scalar mediators followed by mediator late decay. We also discuss the Breit-Wigner resonant enhancement and Sommerfeld enhancement in case the s-wave annihilation process is small or helicity suppressed. With both types of enhancement, constraints on the parameters can be relaxed and new ways for model building will be open in explaining the DAMPE results.
The inflationary paradigm is extremely successful regarding predictions of temperature anisotropies in the CMB. However, inflation also makes predictions for a CMB B-mode polarization, which has not been detected. Moreover, the standard inflationary paradigm is unable to accommodate the evolution from the initial state, which is assumed to be symmetric, into a non-symmetric aftermath. In Phys. Rev. D 96, 101301(R), we show that the incorporation of an element capable of explaining such a transition drastically changes the prediction for the shape and size of the B-mode spectrum. In particular, employing a realistic objective collapse model in a well-defined semiclassical context, we find that, while predictions of temperature anisotropies are nor altered (with respect to standard predictions), the B-mode spectrum gets strongly suppressed---in accordance with observations. Here we present an in-depth discussion of that analysis, together with the details of the calculation.
Understanding infrared (IR) luminosity is fundamental to understanding the
cosmic star formation history and AGN evolution. Japanese infrared satellite,
AKARI, provided unique data sets to probe this both at low and high redshift;
the AKARI all sky survey in 6 bands (9-160 $\mu$m), and the AKARI NEP survey in
9 bands (2-24$\mu$m). The AKARI performed all sky survey in 6 IR bands (9, 18,
65, 90, 140, and 160 $\mu$m) with 3-10 times better sensitivity than IRAS,
covering the crucial far-IR wavelengths across the peak of the dust emission.
Combined with a better spatial resolution, we measure the total infrared
luminosity ($L_{TIR}$) of individual galaxies, and thus, the total infrared
luminosity density of the local Universe much more precisely than previous
work. In the AKARI NEP wide field, AKARI has obtained deep images in the
mid-infrared (IR), covering 5.4 deg$^2$. However, our previous work was limited
to the central area of 0.25 deg$^2$ due to the lack of deep optical coverage.
To rectify the situation, we used the newly advent Subaru telescope's Hyper
Suprime-Cam to obtain deep optical images over the entire 5.4 deg$^2$ of the
AKARI NEP wide field.
With this deep and wide optical data, we, for the first time, can use the
entire AKARI NEP wide data to construct restframe 8$\mu$m, 12$\mu$m, and total
infrared (TIR) luminosity functions (LFs) at 0.15$<z<$2.2. A continuous 9-band
filter coverage in the mid-IR wavelength (2.4, 3.2, 4.1, 7, 9, 11, 15, 18, and
24$\mu$m) by the AKARI satellite allowed us to estimate restframe 8$\mu$m and
12$\mu$m luminosities without using a large extrapolation based on a SED fit,
which was the largest uncertainty in previous work. By combining these two
results, we reveal dust-hidden cosmic star formation history and AGN evolution
from z=0 to z=2.2, all probed by the AKARI satellite.
In this proceedings application of a fuzzy Support Vector Machine (FSVM) learning algorithm, to classify mid-infrared (MIR) sources from the AKARI NEP Deep field into three classes: stars, galaxies and AGNs, is presented. FSVM is an improved version of the classical SVM algorithm, incorporating measurement errors into the classification process; this is the first successful application of this algorithm in the astronomy. We created reliable catalogues of galaxies, stars and AGNs consisting of objects with MIR measurements, some of them with no optical counterparts. Some examples of identified objects are shown, among them O-rich and C-rich AGB stars.
We study a simple TeV-scale model of baryon number violation which explains the observed proximity of the dark matter and baryon abundances. The model has constraints arising from both low and high-energy processes, and in particular, predicts a sizeable rate for the neutron-antineutron ($n-\bar{n}$) oscillation at low energy and the monojet signal at the LHC. We find an interesting complementarity among the constraints arising from the observed baryon asymmetry, ratio of dark matter and baryon abundances, $n-\bar{n}$ oscillation lifetime and the LHC monojet signal. There are regions in the parameter space where the $n-\bar{n}$ oscillation lifetime is found to be more constraining than the LHC constraints, which illustrates the importance of the next-generation $n-\bar{n}$ oscillation experiments.
We systematically investigate shear-free cosmological models realized by p-form gauge fields; a scenario in which anisotropic spatial sections expand isotropically with expansion histories equivalent to standard FLRW models. Specifically, we present a complete list of general relativistic shear-free solutions in a class of anisotropic, spatially homogeneous and orthogonal cosmological models containing a collection of $n$ independent $p$-form gauge fields, where $p\in\{0,1,2,3\}$, in addidtion to standard LCDM matter fields modelled as perfect fluids. Here a (collection of) gauge field(s) balances anisotropic spatial curvature on the right-hand side of the shear propagation equation. The result is a class of solutions dynamically equivalent to standard FLRW cosmologies, with an effective curvature constant $K_\text{eff}$ that depends both on spatial curvature and the energy density of the gauge field(s). In the case of a single gauge field ($n=1$) we show that the only spacetimes that admit such solutions are the LRS Bianchi type III, Bianchi type VI$_0$ and Kantowski-Sachs metric, which are dynamically equivalent to open ($K_\text{eff}<0$), flat ($K_\text{eff}=0$) and closed ($K_\text{eff}>0$) FLRW models, respectively. With a collection of gauge fields ($n>1$) also Bianchi type II admits a shear-free solution ($K_\text{eff}>0$). We identify the LRS Bianchi type III solution to be the unique shear-free solution with a gauge field Hamiltonian bounded from below in the entire class of models. This is a generalization and unification of the shear-free solutions discovered by Carneiro et. al. (2001) with a massless scalar field and by Koivisto et. al. (2011) with a 2-form gauge field, which we show are physically equivalent at the field strength $p+1$ level. Along the way we develop strategies and a framework that can be utilized in a broader class of models.
Links to: arXiv, form interface, find, astro-ph, recent, 1712, contact, help (Access key information)