We develop a method for constraining the scaling relation of optical richness ($\lambda$) with halo masses ($M$) for a sample of galaxy clusters based on a forward modeling approach for the Planck cosmology, where we model the probability distribution of optical richness for a given mass, $P(\ln \lambda| M)$. To model the abundance and the stacked lensing profiles, we employ the halo emulator that outputs the halo mass function and the stacked lensing profile for an arbitrary set of halo mass and redshift, calibrated based on a suite of high-resolution $N$-body simulations. By applying our method to 8,312 SDSS redMaPPer clusters with $20\le \lambda \le 100$ and $0.10\le z_{\lambda}\le0.33$, we show that the log-normal distribution model for $P(\ln \lambda|M)$, with four free parameters, well reproduces the measured abundances and lensing profiles simultaneously. The constraints are characterized by the mean relation, $\left\langle \ln{\lambda}\right\rangle(M)=A+B\ln(M/M_{\rm pivot})$, with $A=3.211^{+0.043}_{-0.045}$ and $B=0.999^{+0.037}_{-0.051}$ (68% CL), for the pivot mass scale $M_{\rm pivot}=3\times 10^{14} h^{-1}M_\odot$, and the scatter $\sigma_{\mathrm{\ln\lambda}|M}=\sigma_0+q\ln(M/M_{\rm pivot})$ with $\sigma_0=0.451^{+0.045}_{-0.037}$ and $q=-0.173^{+0.032}_{-0.024}$. However, we find that a large scatter for the low richness bins, especially $20\le \lambda \lesssim 30$, is required in order for the model to reproduce the measurements by the contributions from low-mass halos for the Planck cosmology. Without such a large scatter, the model prediction for the lensing profiles tends to overestimate the measured amplitudes. This might imply a possible contamination of low-richness clusters due to the projection effects. Such a low-mass halo contribution is significantly reduced when applying our method to the sample of $30\le \lambda \le 100$.
By using N-body hydrodynamical cosmological simulations in which the chemistry of major metals and molecules is consistently solved for, we study the interaction of metallic fine-structure lines with the CMB. Our analysis shows that the collisional induced emissions in the OI 145 $\mu$m and CII 158 $\mu$m lines during reionization introduce a distortion of the CMB spectrum at low frequencies ($\nu < 300$ GHz) with amplitudes up to $\Delta I_{\nu}/B_{\nu}(T_{\rm CMB})\sim 10^{-8}$-$10^{-7}$, i.e., at the $\sim 0.1$ percent level of FIRAS upper limits. Shorter wavelength fine-structure transitions (OI 63 $\mu$m, FeII 26 $\mu$m, and SiII 35 $\mu$m) typically sample the reionization epoch at higher observing frequencies ($\nu > 400$ GHz). This corresponds to the Wien tail of the CMB spectrum and the distortion level induced by those lines may be as high as $\Delta I_{\nu}/B_{\nu}(T_{\rm CMB})\sim 10^{-4}$. The angular anisotropy produced by these lines should be more relevant at higher frequencies: while practically negligible at $\nu=145 $GHz, signatures from CII 158 $\mu$m and OI 145 $\mu$m should amount to 1%-5% of the anisotropy power measured at $l \sim 5000$ and $\nu=220 $GHz by the ACT and SPT collaborations (after assuming $\Delta \nu_{\rm obs}/\nu_{\rm obs}\simeq 0.005$ for the line observations). Our simulations show that anisotropy maps from different lines (e.g., OI 145 $\mu$m and CII 158 $\mu$m) at the same redshift show a very high degree ($>0.8$) of spatial correlation, allowing for the use of observations at different frequencies to unveil the same snapshot of the reionization epoch. Finally, our simulations demonstrate that line-emission anisotropies extracted in narrow frequency/redshift shells are practically uncorrelated in frequency space, thus enabling standard methods for removal of foregrounds that vary smoothly in frequency, just as in HI 21 cm studies.
In this paper we estimate diffuse foreground minimized Cosmic Microwave Background (CMB) Stokes Q and U polarization maps based upon the fundamental concept of Gaussian nature of CMB and strong non-Gaussian nature of astrophysical polarized foregrounds using WMAP nine year published polarization maps. We excise regions of the sky that define position of the known point sources, regions that are strongly contaminated by either the detector noise or by the diffuse foregrounds or both, and then perform foreground minimizations over the surviving sky regions that constitute approximately $50\%$ of the full sky area. We critically evaluate performance of foreground minimizations in several ways and show that our foreground minimization method removes significant foregrounds from input maps. The cleaned Stokes \{Q, U\} polarization maps have less EE and BB power from relevant sky region compared to WMAP foreground-reduced Stokes \{Q, U\} polarization maps at different multipole ranges. We validate our methodology by performing detailed Monte Carlo simulations. The main driving machinery of our method is an internal-linear-combination (ILC) approach, however, unlike simple variance minimization performed in usual ILC method, the core of the method being dependent on the theoretically well motivated concept of Gaussianity of CMB polarization a direct connection is established between observations and models of inflation. Additionally, the method, like the usual ILC method, is independent on modeling uncertainties of polarized foregrounds. It will be useful to apply our method in future generation low-noise CMB polarization experiments.
The 21-cm signal from the Epoch of Reionization (EoR) is expected to be detected in the next few years, either with existing instruments or by the upcoming SKA and HERA projects. In this context there is a pressing need for publicly available high-quality templates covering a wide range of possible signals. These are needed both for end-to-end simulations of the up-coming instruments, as well as to develop signal analysis methods. In this work we present such a set of templates, publicly available, for download at https://21ssd.obspm.fr/. The database contains 21-cm brightness temperature lightcones at high and low resolution, and several derived statistical quantities for 45 models spanning our choice of 3D parameter space. These data are the result of fully coupled radiative hydrodynamic high resolution ($1024^3$) simulations performed with the LICORICE code. Both X-ray and Lyman line transfer is performed to account for heating and Wouthuysen-Field coupling fluctuations. We also present a first exploitation of the data using the power spectrum and the Pixel Distribution Function (PDF) as functions of redshifts, computed from lightcone data. We analyse how these two quantities behave when varying the model parameters while taking into account the thermal noise expected of a typical SKA survey. Finally, we show that the power spectrum and the PDF have different -- and to some extent complementary -- abilities to distinguish between different models. This opens the door to formulating an optimal sampling of the parameter space, dependant on the chosen diagnostics.
Successive releases of Planck data have demonstrated the strength of the Sunyaev--Zeldovich (SZ) effect in detecting hot baryons out to the galaxy cluster peripheries. To infer the hot gas pressure structure from nearby galaxy clusters to more distant objects, we developed a parametric method that models the spectral energy distribution and spatial anisotropies of both the Galactic thermal dust and the Cosmic Microwave Background, that are mixed-up with the cluster SZ and dust signals. Taking advantage of the best angular resolution of the High Frequency Instrument channels (5 arcmin) and using X-ray priors in the innermost cluster regions that are not resolved with Planck, this modelling allowed us to analyze a sample of 61 nearby members of the Planck catalog of SZ sources ($0 < z < 0.5$, $\tilde{z} = 0.15$) using the full mission data, as well as to examine a distant sample of 23 clusters ($0.5 < z < 1$, $\tilde{z} = 0.56$) that have been recently followed-up with XMM-Newton and Chandra observations. We find that (i) the average shape of the mass-scaled pressure profiles agrees with results obtained by the Planck collaboration in the nearby cluster sample, and that (ii) no sign of evolution is discernible between averaged pressure profiles of the low- and high-redshift cluster samples. In line with theoretical predictions for these halo masses and redshift ranges, the dispersion of individual profiles relative to a self-similar shape stays well below 10 % inside $r_{500}$ but increases in the cluster outskirts.
Lensing of the CMB is now a well-developed probe of large-scale clustering over a broad range of redshifts. By exploiting the non-Gaussian imprints of lensing in the polarization of the CMB, the CORE mission can produce a clean map of the lensing deflections over nearly the full-sky. The number of high-S/N modes in this map will exceed current CMB lensing maps by a factor of 40, and the measurement will be sample-variance limited on all scales where linear theory is valid. Here, we summarise this mission product and discuss the science that it will enable. For example, the summed mass of neutrinos will be determined to an accuracy of 17 meV combining CORE lensing and CMB two-point information with contemporaneous BAO measurements, three times smaller than the minimum total mass allowed by neutrino oscillations. In the search for B-mode polarization from primordial gravitational waves with CORE, lens-induced B-modes will dominate over instrument noise, limiting constraints on the gravitational wave power spectrum amplitude. With lensing reconstructed by CORE, one can "delens" the observed polarization internally, reducing the lensing B-mode power by 60%. This improves to 70% by combining lensing and CIB measurements from CORE, reducing the error on the gravitational wave amplitude by 2.5 compared to no delensing (in the null hypothesis). Lensing measurements from CORE will allow calibration of the halo masses of the 40000 galaxy clusters that it will find, with constraints dominated by the clean polarization-based estimators. CORE can accurately remove Galactic emission from CMB maps with its 19 frequency channels. We present initial findings that show that residual Galactic foreground contamination will not be a significant source of bias for lensing power spectrum measurements with CORE. [abridged]
Cosmological models with Galileon gravity are an alternative to the standard $\Lambda {\rm CDM}$ paradigm with testable predictions at the level of its self-accelerating solutions for the expansion history, as well as large-scale structure formation. Here, we place constraints on the full parameter space of these models using data from the cosmic microwave background (CMB) (including lensing), baryonic acoustic oscillations (BAO) and the Integrated Sachs-Wolfe (ISW) effect. We pay special attention to the ISW effect for which we use the cross-spectra, $C_\ell^{\rm T g}$, of CMB temperature maps and foreground galaxies from the WISE survey. The sign of $C_\ell^{\rm T g}$ is set by the time evolution of the lensing potential in the redshift range of the galaxy sample: it is positive if the potential decays (like in $\Lambda {\rm CDM}$), negative if it deepens. We constrain three subsets of Galileon gravity separately known as the Cubic, Quartic and Quintic Galileons. The cubic Galileon model predicts a negative $C_\ell^{\rm T g}$ and exhibits a $7.8\sigma$ tension with the data, which effectively rules it out. For the quartic and quintic models the ISW data also rule out a significant region of the parameter space but permit regions where the goodness-of-fit is comparable to $\Lambda {\rm CDM}$. The data prefers a non zero sum of the neutrino masses ($\sum m_\nu\approx 0.5$eV) with $ \sim 5\sigma$ significance in these models. The best-fitting models have values of $H_0$ consistent with local determinations, thereby avoiding the tension that exists in $\Lambda {\rm CDM}$. We also identify and discuss a $\sim 2\sigma$ tensions that Galileon gravity exhibits with recent BAO measurements. Our analysis shows overall that Galileon cosmologies cannot be ruled out by current data but future lensing, BAO and ISW data hold strong potential to do so.
We present a new Bayesian algorithm making use of Markov Chain Monte Carlo sampling that allows us to simultaneously estimate the unknown continuum level of each quasar in an ensemble of high-resolution spectra, as well as their common probability distribution function (PDF) for the transmitted Ly$\alpha$ forest flux. This fully automated PDF regulated continuum fitting method models the unknown quasar continuum with a linear Principal Component Analysis (PCA) basis, with the PCA coefficients treated as nuisance parameters. The method allows one to estimate parameters governing the thermal state of the intergalactic medium (IGM), such as the slope of the temperature-density relation $\gamma-1$, while marginalizing out continuum uncertainties in a fully Bayesian way. Using realistic mock quasar spectra created from a simplified semi-numerical model of the IGM, we show that this method recovers the underlying quasar continua to a precision of $\simeq7\%$ and $\simeq10\%$ at $z=3$ and $z=5$, respectively. Given the number of principal component spectra, this is comparable to the underlying accuracy of the PCA model itself. Most importantly, we show that we can achieve a nearly unbiased estimate of the slope $\gamma-1$ of the IGM temperature-density relation with a precision of $\pm8.6\%$ at $z=3$, $\pm6.1\%$ at $z=5$, for an ensemble of ten mock high-resolution quasar spectra. Applying this method to real quasar spectra and comparing to a more realistic IGM model from hydrodynamical simulations would enable precise measurements of the thermal and cosmological parameters governing the IGM, albeit with somewhat larger uncertainties given the increased flexibility of the model.
We consider FRW cosmology in $f(R)= R+ \gamma R^2+\delta R^3$ modified framework. The Palatini approach reduces its dynamics to the simple generalization of Friedmann equation. Thus we study the dynamics in two-dimensional phase space with some details. After reformulation of the model in the Einstein frame, it reduces to the FRW cosmological model with a homogeneous scalar field and vanishing kinetic energy term. This potential determines the running cosmological constant term as a function of the Ricci scalar. As a result we obtain the emergent dark energy parametrization from the covariant theory. We study also singularities of the model and demonstrate that in the Einstein frame some undesirable singularities disappear.
We study the dependence of galaxy clustering on atomic gas mass using a sample of $\sim$16,000 galaxies with redshift in the range of $0.0025<z<0.05$ and HI mass of $M_{\rm HI}>10^8M_{\odot}$, drawn from the 70% complete sample of the Arecibo Legacy Fast ALFA survey. We construct subsamples of galaxies with $M_{\rm HI}$ above different thresholds, and make volume-limited clustering measurements in terms of three statistics: the projected two-point correlation function, the projected cross-correlation function with respect to a reference sample selected from the Sloan Digital Sky Survey, and the redshift-space monopole moment. In contrast to previous studies, which found no/weak HI-mass dependence, we find both the clustering amplitude on scales above a few Mpc and the bias factors to increase significantly with increasing HI mass for subsamples with HI mass thresholds above $10^9M_{\odot}$. For HI mass thresholds below $10^9M_{\odot}$, while the measurements have large uncertainties caused by the limited survey volume and sample size, the inferred galaxy bias factors are systematically lower than the minimum halo bias factor from mass-selected halo samples. The simple halo model, in which galaxy content is only determined by halo mass, has difficulties in interpreting the clustering measurements of the HI-selected samples. We extend the simple model by including the halo formation time as an additional parameter. A model that puts HI-rich galaxies into halos that formed late can reproduce the clustering measurements reasonably well. We present the implications of our best-fitting model on the correlation of HI mass with halo mass and formation time, as well as the halo occupation distributions and HI mass functions for central and satellite galaxies. These results are compared with the predictions from semi-analytic galaxy formation models and hydrodynamic galaxy formation simulations.
We present a novel cosmological solution in the framework of extended quasidilaton theory which underwent scrutiny recently. We only consider terms that do not generate the Boulware-Deser degree of freedom, hence the "ghost-free" quasidilaton theory, and show three new branches of cosmological evolution therein. One of the solutions passes the perturbative stability tests. This new solution exhibits a late time self-acceleration and all graviton polarizations acquire masses that converge to a constant in the asymptotic future. Moreover, all modes propagate at the speed of light. We propose that this solution can be used as a benchmark model for future phenomenological studies.
Links to: arXiv, form interface, find, astro-ph, recent, 1707, contact, help (Access key information)
We measure the Cosmic Microwave Background (CMB) skewness power spectrum in $\textit{Planck}$, using frequency maps of the HFI instrument and the Sunyaev-Zel'dovich (SZ) component map. The two-to-one skewness power spectrum measures the cross-correlation between CMB lensing and the thermal SZ effect. We also directly measure the same cross-correlation using $\textit{Planck}$ CMB lensing map and the SZ map and compare it to the cross-correlation derived from the skewness power spectrum. We model fit the SZ power spectrum and CMB lensing-SZ cross power spectrum via the skewness power spectrum to constrain the gas pressure profile of dark matter halos. The gas pressure profile is compared to existing measurements in the literature including a direct estimate based on the stacking of SZ clusters in $\textit{Planck}$.
We study the clustering of galaxies (from ~ $0.1$ to a few tens Mpc scales) at high redshift ($z>7$) using the BLUETIDES simulation and compare it to current constraints from Hubble legacy and Hyper Suprime Cam (HSC) fields. BLUETIDES is the largest high resolution cosmological hydrodynamic simulation with a box length of $400$ $Mpc/h$ and $0.7$ trillion particles (~ $3$ million star forming galaxies). We select galaxies with magnitudes consistent with the observed fields (~ $10^5$ star forming galaxies at $z$ ~ $8$) and measure their clustering. The 1-halo term dominates on scales $r < 0.1$ $Mpc/h$ ($\Theta < 3"$) with non-linear effect at transition scales between the 1-halo and 2-halo term affecting scales $0.1 < r < 20$ $Mpc/h$ ($3" < \Theta < 90"$). The predicted angular clustering amplitude of BLUETIDES galaxies is in good agreement with current measurements. The measured linear bias at $z = 8$ is $7.7 \pm 0.5$ (evolving close to linearly with redshift). This is consistent with the inferred bias from observations. The typical halo masses are $M_H$ ~ $2 \times 10^{10}$ $M_{\odot}/h$, which is only slightly lower than that inferred from HOD models used to constrain observed clustering. Using the simulations we show that the current clustering measurements probe the transition between 1-halo to 2-halo regime so non-linear effects are important at all scales. We measure a strong scale dependence in the bias at scales below a few Mpc/h ($\Theta < 30"$). The occupation numbers for satellites in BLUETIDES is somewhat higher than standard HODs adopted in these analyses. That should imply a higher number of galaxies detected by JWST observing these fields. At $z$ ~ $8$ BLUETIDES predicts enhanced clustering (by factors ~ $2$) within the 1-halo and 1-halo to 2-halo transition scales as compared to linear HOD models; these differences increase (up to factor ~ $10$) at higher redshifts.
Baryonic effects are amongst the most severe systematics to the tomographic analysis of weak lensing data which is the principal probe in many future generations of cosmological surveys like LSST, Euclid etc.. Modeling or parameterizing these effects is essential in order to extract valuable constraints on cosmological parameters. In a recent paper, Eifler et al. (2015) suggested a reduction technique for baryonic effects by conducting a principal component analysis (PCA) and removing the largest baryonic eigenmodes from the data. In this article, we conducted the investigation further and addressed two critical aspects. Firstly, we performed the analysis by separating the simulations into training and test sets, computing a minimal set of principle components from the training set and examining the fits on the test set. We found that using only four parameters, corresponding to the four largest eigenmodes of the training set, the test sets can be fitted thoroughly with an RMS $\sim 0.0011$. Secondly, we explored the significance of outliers, the most exotic/extreme baryonic scenarios, in this method. We found that excluding the outliers from the training set results in a relatively bad fit and degraded the RMS by nearly a factor of 3. Therefore, for a direct employment of this method to the tomographic analysis of the weak lensing data, the principle components should be derived from a training set that comprises adequately exotic but reasonable models such that the reality is included inside the parameter domain sampled by the training set. The baryonic effects can be parameterized as the coefficients of these principle components and should be marginalized over the cosmological parameter space.
We explore the possibility that Fast Radio Bursts are due to the annihilation of cusps on cosmic string loops. We compute the energy released in the annihilation events in the radio region, the expected event rate, and the time scale of the bursts. We find that the energy and event rates are sufficiently high and the time scale is sufficiently small to explain the current data. We predict how the event rate will change as the resolution of telescopes improves. Since the burst rate depends on the string tension, future data will allow the determination of the tension.
Primordial Black Holes (PBH) could be the cold dark matter of the universe. They could have arisen from large (order one) curvature fluctuations produced during inflation that reentered the horizon in the radiation era. At reentry, these fluctuations source gravitational waves (GW) via second order anisotropic stresses. These GW, together with those (possibly) sourced during inflation by the same mechanism responsible for the large curvature fluctuations, constitute a primordial stochastic GW background (SGWB) that unavoidably accompanies the PBH formation. We study how the amplitude and the range of frequencies of this signal depend on the statistics (Gaussian versus $\chi^2$) of the primordial curvature fluctuations, and on the evolution of the PBH mass function due to accretion and merging. We then compare this signal with the sensitivity of present and future detectors, at PTA and LISA scales. We find that this SGWB will help to probe, or strongly constrain, the early universe mechanism of PBH production. The comparison between the peak mass of the PBH distribution and the peak frequency of this SGWB will provide important information on the merging and accretion evolution of the PBH mass distribution from their formation to the present era. Different assumptions on the statistics and on the PBH evolution also result in different amounts of CMB $\mu$-distortions. Therefore the above results can be complemented by the detection (or the absence) of $\mu$-distortions with an experiment such as PIXIE.
We study the accelerating present universe in terms of the time evolution of the equation of state $w(z)$ due to the thawing and freezing potentials in the quintessence model. From the observations of the derivatives $dw/da$ and $d^2w/da^2$ ($a=1/(1+z)$) at the scale factor of $a = 1$, the characteristic parameters of each potential are fixed and then the theoretical time development of the scalar field $Q$ in the potential has been numerically calculated from the present ($z=0$) to the past. The possible time evolution of $w(z)$ as a function of the redshift $z$ has been shown for each scalar potential and compared with available observations. So the allowed regions of $dw/da-d^2w/da^2$ space are estimated from the observational constraint of $w(z) \ (z \le 2)$ for the investigated potentials. There are some cases which show the thawing feature in the freezing type potentials. It is also discussed the dependence of the region on the parameter $\Delta \ (w=-1+\Delta)$ at present and the other observational constraint $w(z) \ (z \le 5)$. If the detailed observation of $dw/da-d^2w/da^2$ is obtained, it will be determined the potential type in the quintessence scenario.
We study models in which neutrino masses are generated dynamically at cosmologically late times. Our study is purely phenomenological and parameterized in terms of three effective parameters characterizing the redshift of mass generation, the width of the transition region, and the present day neutrino mass. We also study the possibility that neutrinos become strongly self-interacting at the time where the mass is generated. We find that in a number of cases, models with large present day neutrino masses are allowed by current CMB, BAO and supernova data. The increase in the allowed mass range makes it possible that a non-zero neutrino mass could be measured in direct detection experiments such as KATRIN. Intriguingly we also find that there are allowed models in which neutrinos become strongly self-interacting around the epoch of recombination.
We present a new exact solution to cosmology as a function of the Hubble constant $H_0$ and matter content $\omega_m$ with accelerated expansion by a dark energy $\Lambda=\omega_0^2$, where $\omega_0=\sqrt{1-q}H$ is a fundamental mode of the cosmological horizon and $q$ is the deceleration parameter. $H_0$ is herein close to the minimum of $H(z)$ consistent with values obtained from surveys of the Local Universe, associated with $Q_0>2.5$, $Q_0=Q(0)$, $Q(z)=dq(z)/dz$, distinct from $H^\prime(z)>0$ and $Q_0\lesssim1$ in $\Lambda$CDM. A dynamical dark energy hereby gives a tension-free estimate of $H_0$ consistent with $H_0$ obtained from surveys of the Local Universe with a combined result $H_0\simeq 73.75\pm 1.44\,\mbox{km\,s}^{-1}\mbox{Mpc}^{-1}.$
We update the constraints on the cosmological parameters by adopting the Planck data released in 2015 and Baryon Acoustic Oscillation (BAO) measurements including the new DR14 quasar sample measurement at redshift $z=1.52$, and we conclude that the based six-parameter $\Lambda$CDM model is preferred. Exploring some extensions to the $\Lambda$CDM models, we find that the equation of state of dark energy reads $w=-1.036\pm 0.056$ in the $w$CDM model, the effective relativistic degrees of freedom in the Universe is $N_\text{eff}=3.09_{-0.20}^{+0.18}$ in the $N_\text{eff}+\Lambda$CDM model and the spatial curvature parameter is $\Omega_k=(1.8\pm1.9)\times 10^{-3}$ in the $\Omega_k+\Lambda$CDM model at $68\%$ confidence level (C.L.), and the $95\%$ C.L. upper bounds on the sum of three active neutrinos masses are $\sum m_\nu<0.16$ eV for the normal hierarchy (NH) and $\sum m_\nu<0.19$ eV for the inverted hierarchy (IH) with $\Delta\chi^2\equiv \chi^2_\text{NH}-\chi^2_\text{IH}=-1.25$.
The lines-of-sight to highly reddened SNe Ia show peculiar continuum polarization curves, growing toward blue wavelengths and peaking at $\lambda_{max} \lesssim 0.4 \mu m$, like no other sight line to any normal Galactic star. We examined continuum polarization measurements of a sample of asymptotic giant branch (AGB) and post-AGB stars from the literature, finding that some PPNe have polarization curves similar to those observed along SNe Ia sight lines. Those polarization curves are produced by scattering on circumstellar dust. We discuss the similarity and the possibility that at least some SNe Ia might explode during the post-AGB phase of their binary companion. Furthermore, we speculate that the peculiar SNe Ia polarization curves might provide observational support to the core-degenerate progenitor model.
At TeV energies and above gamma rays can induce electromagnetic cascades, whose charged component is sensitive to intervening intergalactic magnetic fields (IGMFs). When interpreting gamma-ray measurements in the energy range between a few GeV and hundreds of TeV, one has to carefully account for effects due to IGMFs, which depend on their strength and power spectrum. Therefore, gamma-ray-induced electromagnetic cascades can be used as probes of cosmic magnetism, since their arrival distribution as well as spectral and temporal properties can provide unique information about IGMFs, whose origin and properties are currently poorly understood. In this contribution we present an efficient three-dimensional Monte Carlo code for simulations of gamma-ray propagation. We focus on the effects of different configurations of IGMFs, in particular magnetic helicity and the power spectrum of stochastic fields, on the morphology of the arrival directions of gamma rays, and discuss the prospects for detecting pair haloes around distant blazars.
The sizes of entire systems of globular clusters (GCs) depend on the formation and destruction histories of the GCs themselves, but also on the assembly, merger and accretion history of the dark matter (DM) haloes that they inhabit. Recent work has shown a linear relation between total mass of globular clusters in the globular cluster system and the mass of its host dark matter halo, calibrated from weak lensing. Here we extend this to GC system sizes, by studying the radial density profiles of GCs around galaxies in nearby galaxy groups. We find that radial density profiles of the GC systems are well fit with a de Vaucouleurs profile. Combining our results with those from the literature, we find tight relationship ($\sim 0.2$ dex scatter) between the effective radius of the GC system and the virial radius (or mass) of its host DM halo. The steep non-linear dependence of this relationship ($R_{e, GCS} \propto R_{200}^{2.5 - 3}$) is currently not well understood, but is an important clue regarding the assembly history of DM haloes and of the GC systems that they host.
In this paper we propose that artificial neural network, the basis of machine learning, is useful to generate the inflationary landscape from a cosmological point of view. Traditional numerical simulations of a global cosmic landscape typically need an exponential complexity when the number of fields is large. However, a basic application of artificial neural network could solve the problem based on the universal approximation theorem of the multilayer perceptron. A toy model in inflation with multiple light fields is investigated numerically as an example of such an application.
We have carried a detailed analysis on the impact of cosmological redshift in the non-parametric approach to automated galaxy morphology classification. We artificially redshifted each galaxy from the EFIGI 4458 sample (re-centered at $z\sim 0$) simulating SDSS, DES, LSST, and HST instruments setups over the range $0 < z < 1.5$. We then traced how the morphometry is degraded in each $z$ using \textsc{Morfometryka}. In the process we re-sampled all catalog to several resolutions and to a diverse SNR range, allowing us to understand the impact of image sampling and noise on our measurements separately. We summarize by exploring the impact of these effects on our capacity to perform automated galaxy supervised morphological classification by investigating the degradation of our classifier's metrics as a function of redshift for each instrument. The overall conclusion is that we can make reliable classification with \textsc{Morfometryka} for $z< 0.2$ with SDSS, for $z<0.5$ with DES, for $z<0.8$ with LSST and for at least $z < 1.5$ with HST.
Links to: arXiv, form interface, find, astro-ph, recent, 1707, contact, help (Access key information)
We investigate the performance of a simple Bayesian fitting approach to correct the cosmic microwave background (CMB) B-mode polarization for gravitational lensing effects in the recovered probability distribution of the tensor-to-scalar ratio. We perform a two-dimensional power spectrum fit of the amplitude of the primordial B-modes (tensor-to-scalar ratio, $r$) and the amplitude of the lensing B-modes (parameter $A_{lens}$), jointly with the estimation of the astrophysical foregrounds including both synchrotron and thermal dust emissions. Using this Bayesian framework, we forecast the ability of the proposed CMB space mission LiteBIRD to constrain $r$ in the presence of realistic lensing and foreground contributions. We compute the joint posterior distribution of $r$ and $A_{lens}$, which we improve by adopting a prior on $A_{lens}$ taken from the South Pole Telescope (SPT) measurement. As it applies to the power spectrum, this approach cannot mitigate the uncertainty on $r$ that is due to E-mode cosmic variance transferred to B-modes by lensing, unlike standard delensing techniques that are performed on maps. However, the method allows to correct for the bias on $r$ induced by lensing, at the expense of a larger uncertainty due to the increased volume of the parameter space. We quantify, for different values of the tensor-to-scalar ratio, the trade-off between bias correction and increase of uncertainty on $r$. For LiteBIRD simulations, which include foregrounds and lensing contamination, we find that correcting the foreground-cleaned CMB B-mode power spectrum for the lensing bias, not the lensing cosmic variance, still guarantees a $3\sigma$ detection of $r=5\times 10^{-3}$. The significance of the detection is increased to $6\sigma$ when the current SPT prior on $A_{lens}$ is adopted.
We investigate in details the gravitational wave (GW) from the first-order phase transition (PT) in the extended standard model of particle physics with a dimension-six operator, which is capable of exhibiting the recently discovered slow first-order PT in addition to the usually studied fast first-order PT. To simplify the discussion, it is sufficient to work with an example of a toy model with the sextic term, and we propose an unified description for both slow and fast first-order PTs. We next study the full one-loop effective potential of the model with fixed/running renormalization-group (RG) scales. Compared to the prediction of GW energy density spectrum from the fixed RG scale, we find that the presence of running RG scale could amplify the peak amplitude by amount of one order of magnitude while shift the peak frequency to the lower frequency regime, and the promising regime of detection within the sensitivity ranges of various space-based GW detectors shrinks down to a lower cut-off value of the sextic term rather than the previous expectation.
We show that an inflation model in which a spectator axion field is coupled to an SU(2) gauge field produces a large three-point function (bispectrum) of primordial gravitational waves, $B_{h}$, on the scales relevant to the cosmic microwave background experiments. The amplitude of the bispectrum at the equilateral configuration is characterized by $B_{h}/P_h^2=\mathcal{O}(10)\times \Omega_A^{-1}$, where $\Omega_A$ is a fraction of the energy density in the gauge field and $P_h$ is the power spectrum of gravitational waves produced by the gauge field.
Here we report the discovery of an extremely massive and large supercluster (called Saraswati) found in the Stripe 82 region of SDSS. This supercluster is a major concentration of galaxies and galaxy clusters, forming a wall-like structure spanning at least 200 Mpc across at the redshift $z \approx 0.3$. This enormous structure is surrounded by a network of galaxy filaments, clusters, and large, $\sim40 - 170$ Mpc diameter, voids. The mean density contrast $\delta$ (relative to the background matter density of the universe) of Saraswati is $\gtrsim 1.62$ and the main body of the supercluster comprises at least 43 massive galaxy clusters (mean $z = 0.28$) with a total mass of $\sim 2 \times 10^{16} M_{\odot}$. The spherical collapse model suggests that the central region of radius $\sim20$ Mpc and mass at least $ 4 \times 10^{15} M_{\odot}$ may be collapsing. This places it among the few largest and most massive superclusters known, comparable to the most massive `Shapley Concentration' ($z \approx 0.046$) in the nearby universe. The Saraswati supercluster and its environs reveal that some extreme large-scale, prominent matter density enhancements had formed $\sim4$ Gy in the past when dark energy had just started to dominate structure formation. This galactic concentration sheds light on the role of dark energy and cosmological initial conditions in supercluster formation, and tests the competing cosmological models.
We study structure formation using relativistic cosmological linear perturbation theory in the presence of intrinsic and relative (with respect to matter) non-adiabatic dark energy perturbations. For different dark energy models we assess the impact of non-adiabaticity on the matter growth promoting a comparison with growth rate data. The dark energy models studied lead to peculiar signatures of the (non)adiabatic nature of dark energy perturbations in the evolution of the $f \sigma_{8}(z)$ observable. We show that non-adiabatic DE models become close to be degenerated with respect to the $\Lambda$CDM model at first order in linear perturbations. This would avoid the identification of the non-adiabatic nature of dark energy using current available data. Therefore, such evidence indicates that new probes are necessary to reveal the non-adiabatic features in the dark energy sector.
A detection of B-mode polarization of the Cosmic Microwave Background (CMB) anisotropies would confirm the presence of a primordial gravitational wave background (GWB). In the inflation paradigm this would be an unprecedented probe of the energy scale of inflation as it is directly proportional to the power spectrum of the GWB. However, similar tensor perturbations can be produced by the matter fields present during inflation, breaking this simple relationship. It is therefore important to be able to distinguish between different generation mechanisms of the GWB. In this paper, we analyse the detectability of a new axion-SU(2) gauge field model using its chiral, scale-dependent tensor spectrum. We forecast the detectability of the resulting CMB TB and EB cross-correlations by the LiteBIRD satellite, considering the effects of residual foregrounds, gravitational lensing, and for the first time assess the ability of such an experiment to jointly detect primordial TB and EB spectra and self-calibrate its polarimeter. We find that LiteBIRD will be able to detect the chiral signal for $r_*>0.03$ with $r_*$ denoting the tensor-to-scalar ratio at the peak scale, and that the maximum signal-to-noise for $r_*<0.07$ is $\sim 2$. We go on to consider an advanced stage of a LISA-like mission, and find that such experiments would complement CMB observations by providing sensitivity to GWB chirality on scales inaccessible to the CMB. We conclude that in order to use the CMB to distinguish this model from a conventional vacuum fluctuation model two-point statistics provide some power, but to achieve high statistical significance we would require higher order statistics which take advantage of the model's non-Gaussianity. On the other hand, in the case of a spectrum peaked at very small scales, inaccessible to the CMB, a highly significant detection could be made using space-based laser interferometers.
If the symmetry breaking responsible for axion dark matter production occurs during the radiation-dominated epoch in the early Universe, then this produces large amplitude perturbations that collapse into dense objects known as axion miniclusters. The characteristic minicluster mass, $M_0$, is set by the mass inside the horizon when axion oscillations begin. For the QCD axion $M_0\sim 10^{-10}M_\odot$, however for an axion-like particle $M_0$ can approach $M_\odot$ or higher. Using the Press-Schechter formalism we compute the mass function of halos formed by hierarchical structure formation from these seeds. We compute the concentrations and collapse times of these halos and show that they can grow to be as massive as $10^6M_0$. Within the halos, miniclusters likely remain tightly bound, and we compute their gravitational microlensing signal taking the fraction of axion dark matter collapsed into miniclusters, $f_{\rm MC}$, as a free parameter. A large value of $f_{\rm MC}$ severely weakens constraints on axion scenarios from direct detection experiments. We take into account the non-Gaussian distribution of sizes of miniclusters and determine how this effects the number of microlensing events. We develop the tools to consider microlensing by an extended mass function of non-point-like objects, and use microlensing data to place the first observational constraints on $f_{\rm MC}$. This opens a new window for the potential discovery of the axion.
We forecast the prospective constraints on the ionized gas model $f_{\rm gas}(z)$ at different evolutionary epochs via the tomographic cross-correlation between kinetic Sunyaev-Zeldovich (kSZ) effect and the reconstructed momentum field at different redshifts. The experiments we consider are the Planck and CMB Stage-4 survey for CMB and the SDSS-III for the galaxy spectroscopic survey. We calculate the tomographic cross-correlation power spectrum, and use the Fisher matrix to forecast the detectability of different $f_{\rm gas}(z)$ models. We find that for constant $f_{\rm gas}$ model, Planck can constrain the error of $f_{\rm gas}$ ($\sigma_{f_{\rm gas}}$) at each redshift bin to $\sim 0.2$, whereas four cases of CMB-S4 can achieve $\sigma_{f_{\rm gas}} \sim 10^{-3}$. For $f_{\rm gas}(z)=f_{\rm gas,0}/(1+z)$ model the error budget will be slightly broadened. We also investigate the model $f_{\rm gas}(z)=f_{\rm gas,0}/(1+z)^{\alpha}$. Planck is unable to constrain the index of redshift evolution, but the CMB-S4 experiments can constrain the index $\alpha$ to the level of $\sigma_{\alpha} \sim 0.01$--$0.1$. The tomographic cross-correlation method will provide an accurate measurement of the ionized gas evolution at different epochs of the Universe.
The baryon-acoustic oscillation (BAO) feature in the Lyman-$\alpha$ forest is one of the key probes of the cosmic expansion rate at redshifts z~2.5, well before dark energy is believed to have become dynamically significant. A key advantage of the BAO as a standard ruler is that it is a sharp feature and hence is more robust against broadband systematic effects than other cosmological probes. However, if the Lyman-$\alpha$ forest transmission is sensitive to the initial streaming velocity of the baryons relative to the dark matter, then the BAO peak position can be shifted. Here we investigate this sensitivity using a suite of hydrodynamic simulations of small regions of the intergalactic medium with a range of box sizes and physics assumptions; each simulation starts from initial conditions at the kinematic decoupling era (z~1059), undergoes a discrete change from neutral gas to ionized gas thermal evolution at reionization (z~8), and is finally processed into a Lyman-$\alpha$ forest transmitted flux cube. Streaming velocities suppress small-scale structure, leading to less violent relaxation after reionization. The changes in the gas distribution and temperature-density relation at low redshift are more subtle, due to the convergent temperature evolution in the ionized phase. The change in the BAO scale is estimated to be of the order of 0.12% at z=2.5; some of the major uncertainties and avenues for future improvement are discussed. The predicted streaming velocity shift would be a subdominant but not negligible effect (of order $0.26\sigma$) for the upcoming DESI Lyman-$\alpha$ forest survey, and exceeds the cosmic variance floor. It is hoped that this study will motivate additional theoretical work on the magnitude of the BAO shift, both in the Lyman-$\alpha$ forest and in other tracers of large-scale structure.
The presence of asymmetry between fermions of opposite handedness in plasmas of relativistic particles can lead to exponential growth of a helical magnetic field via a small-scale chiral dynamo instability known as the chiral magnetic effect. Here we show, using dimensional arguments and numerical simulations, that this process produces through the Lorentz force chiral-magnetically driven turbulence. A k^{-2} magnetic energy spectrum emerges via inverse transfer over a certain range of wavenumbers k. The total chirality (magnetic helicity plus normalized chiral chemical potential) is conserved in this system. Therefore, as the helical magnetic field grows, most of the total chirality gets transferred into magnetic helicity until the chiral magnetic effect terminates. Quantitative results for height, slope, and extent of the spectrum are obtained. Consequences of this effect for cosmic magnetic fields are discussed.
The cosmological coincidences between the matter and radiation energy densities at recombination as well as between the densities of matter and the cosmological constant at present time are well known. We point out that moreover the third intersection between the energy densities of radiation and the cosmological constant coincides with the reionization epoch. To quantify the statistical relevance of this concurrence, we compute the Bayes factor between the concordance cosmology with free Thomson scattering optical depth and a model for which this parameter is inferred from imposing a match between the time of density equality and the epoch of reionization. This is to characterize the potential explanatory gain if one were to find a parameter-free physical connection. We find a very strong preference for such a concurrence on the Jeffreys scale from current cosmological observations. We furthermore discuss the effect of choice of priors, changes in reionization history, and free sum of neutrino masses. We also estimate the impact of adding intermediate polarization data from the Planck High Frequency Instrument and prospects for future 21 cm surveys. In the first case, preference for the correlation remains substantial, whereas future data may give results more decisive in pro or substantial in contra. Finally, we provide a discussion on different interpretations of these findings. In particular, we show how a connection between the star-formation history and the cosmological background dynamics can give rise to this concurrence.
We study gravitational-wave production from bubble dynamics (bubble collisions and sound waves) during a cosmic first-order phase transition with an analytic approach. We model the system with the thin-wall approximation but without the envelope approximation often adopted in the literature. We first write down analytic expressions for the gravitational-wave spectrum, and then evaluate them with numerical methods. It is found that, in the long-lasting limit of the collided walls, the spectrum grows from $\propto f^3$ to $\propto f^1$ for low frequencies, showing a significant enhancement compared to the one with the envelope approximation. It is also found that the spectrum saturates in the same limit, indicating a decrease in the correlation of the energy-momentum tensor at late times. The bubble walls in our setup are considered as modeling the scalar field configuration and/or the bulk motion of the fluid, and therefore our results have implications to gravitational-wave production from both bubble collisions (scalar dynamics) and sound waves (fluid dynamics).
We conduct a comprehensive study of the effects of incorporating galaxy morphology information in photometric redshift estimation. Using machine learning methods, we assess the changes in the scatter and catastrophic outlier fraction of photometric redshifts when galaxy size, ellipticity, S\'{e}rsic index and surface brightness are included in training on galaxy samples from the SDSS and the CFHT Stripe-82 Survey (CS82). We show that by adding galaxy morphological parameters to full $ugriz$ photometry, only mild improvements are obtained, while the gains are substantial in cases where fewer passbands are available. For instance, the combination of $grz$ photometry and morphological parameters almost fully recovers the metrics of $5$-band photometric redshifts. We demonstrate that with morphology it is possible to determine useful redshift distribution $N(z)$ of galaxy samples without any colour information. We also find that the inclusion of quasar redshifts and associated object sizes in training improves the quality of photometric redshift catalogues, compensating for the lack of a good star-galaxy separator. We further show that morphological information can mitigate biases and scatter due to bad photometry. As an application, we derive both point estimates and posterior distributions of redshifts for the official CS82 catalogue, training on morphology and SDSS Stripe-82 $ugriz$ bands when available. Our redshifts yield a 68th percentile error of $0.058(1+z)$, and a catastrophic outlier fraction of $5.2$ per cent. We further include a deep extension trained on morphology and single $i$-band CS82 photometry.
Links to: arXiv, form interface, find, astro-ph, recent, 1707, contact, help (Access key information)
We introduce the IllustrisTNG project, a new suite of cosmological magneto-hydrodynamical simulations performed with the moving-mesh code AREPO employing an updated Illustris galaxy formation model. Here we focus on the general properties of magnetic fields and the diffuse radio emission in galaxy clusters. Magnetic fields are prevalent in galaxies, and their build-up is closely linked to structure formation. We find that structure formation amplifies the initial seed fields ($10^{-14}$ comoving Gauss) to the values observed in low-redshift galaxies ($1-10\,\mu G$). The magnetic field topology is closely connected to galaxy morphology such that irregular fields are hosted by early-type galaxies, while large-scale, ordered fields are present in disc galaxies. Using a simple model for the energy distribution of relativistic electrons we predict the diffuse radio emission of $280$ clusters with a baryonic mass resolution of $1.1\times 10^{7}\,{\rm M_{\odot}}$, and generate mock observations for VLA, LOFAR, ASKAP and SKA. Our simulated clusters show extended radio emission, whose detectability correlates with their virial mass. We reproduce the observed scaling relations between total radio power and X-ray emission, $M_{500}$, and the Sunyaev-Zel'dovich $Y_{\rm 500}$ parameter. The radio emission surface brightness profiles of our most massive clusters are in reasonable agreement with VLA measurements of Coma and Perseus. Finally, we discuss the fraction of detected extended radio haloes as a function of virial mass and source count functions for different instruments. Overall our results agree encouragingly well with observations, but a refined analysis requires a more sophisticated treatment of relativistic particles in large-scale galaxy formation simulations.
We examine the statistics of the low-redshift Ly-alpha forest in an adaptive mesh refinement hydrodynamic cosmological simulation of sufficient volume to include distinct large-scale environments. We compare our HI column density distribution of absorbers both with recent work and between two highly-refined regions of our simulation: a large-scale overdensity and a large-scale underdensity (on scales of approximately 20 Mpc). We recover the average results presented in Kollmeier et al. (2014) using different simulation methods. We further break down these results as a function of environment to examine the detailed dependence of absorber statistics on large-scale density. We find that the slope of the HI column density distribution in the 10$^{12.5}$ $\le$ N$_{HI}$/cm$^{-2}$ $\le$ 10$^{14.5}$ range depends on environment such that the slope becomes steeper for higher environmental density, and this difference reflects distinct physical conditions of the intergalactic medium on these scales. We track this difference to the different temperature structures of filaments in varying environments. Specifically, filaments in the overdensity are hotter and, correspondingly, are composed of gas with lower HI fractions than those in underdense environments. Our results highlight that in order to understand the physics driving the HI CDD, we need not only improved accounting of the sources of ionizing UV photons, but also of the physical conditions of the IGM and how this may vary as a function of large-scale environment.
Visinelli and Gondolo (2015, hereafter VG15) derived analytic expressions for the evolution of the dark matter temperature in a generic cosmological model. They then calculated the dark matter kinetic decoupling temperature $T_{\mathrm{kd}}$ and compared their results to the Gelmini and Gondolo (2008, hereafter GG08) calculation of $T_{\mathrm{kd}}$ in an early matter-dominated era (EMDE), which occurs when the Universe is dominated by either a decaying oscillating scalar field or a semistable massive particle before Big Bang nucleosynthesis. VG15 found that dark matter decouples at a lower temperature in an EMDE than it would in a radiation-dominated era, while GG08 found that dark matter decouples at a higher temperature in an EMDE than it would in a radiation-dominated era. VG15 attributed this discrepancy to the presence of a matching constant that ensures that the dark matter temperature is continuous during the transition from the EMDE to the subsequent radiation-dominated era and concluded that the GG08 result is incorrect. We show that the disparity is due to the fact that VG15 compared $T_\mathrm{kd}$ in an EMDE to the decoupling temperature in a radiation-dominated universe that would result in the same dark matter temperature at late times. Since decoupling during an EMDE leaves the dark matter colder than it would be if it decoupled during radiation domination, this temperature is much higher than $T_\mathrm{kd}$ in a standard thermal history, which is indeed lower than $T_{\mathrm{kd}}$ in an EMDE, as stated by GG08.
We show that a subdominant component of dissipative dark matter resembling the Standard Model can form many intermediate-mass black hole seeds during the first structure formation epoch, such that a few percent of all the dwarf galaxies host one. We also observe that, in the presence of this matter sector, the black holes will grow at a much faster rate with respect to the ordinary case. These facts can explain the observed abundance of supermassive black holes feeding high-redshift quasars. The scenario will have interesting observational consequences, such as dark substructures and gravitational wave production.
To derive a power spectrum for energy density inhomogeneities in a closed universe, we study a spatially-closed inflation model. The constants of integration in the inflation epoch solutions are determined from Hawking's prescription that the quantum state of the universe only include field configurations that are regular on the Euclidean section. The power spectrum we find is not a power law, and depends on spatial wavenumber in the way expected for a generalization to the closed model of the standard flat-space scale-invariant power spectrum. (Abridged.)
We study Planck 2015 cosmic microwave background (CMB) anisotropy data using the energy density inhomogeneity power spectrum generated by quantum fluctuations during an early epoch of inflation in the non-flat $\Lambda$CDM model. Unlike earlier analyses of non-flat models, which assumed an inconsistent power-law power spectrum of energy density inhomogeneities, we find that the Planck 2015 data alone, and also in conjunction with baryon acoustic oscillation measurements, are reasonably well fit by a closed $\Lambda$CDM model in which spatial curvature contributes a few percent of the current cosmological energy density budget. In this model, the measured Hubble constant and non-relativistic matter density parameter are in good agreement with values determined using most other data. Depending on parameter values, the closed $\Lambda$CDM model has reduced power, relative to the tilted, spatially-flat $\Lambda$CDM case, and can partially alleviate the low multipole CMB temperature anisotropy deficit and can help partially reconcile the CMB anisotropy and weak lensing $\sigma_8$ constraints, at the expense of slightly worsening the fit to higher multipole CMB temperature anisotropy data. Our results are interesting but tentative; a more thorough analysis is needed to properly gauge their significance.
Results of recent experiment reinstate feasibility to the hypothesis that biomolecular homochirality originates from beta decay. Coupled with hints that this process occurred extraterrestrially suggests aluminum-26 as the most likely source. If true, then its appropriateness is highly dependent on the half-life and energy of this decay. Demanding that this mechanism hold places new constraints on the anthropically allowed range for multiple parameters, including the electron mass, difference between up and down quark masses, the fine structure constant, and the electroweak scale. These new constraints on particle masses are tighter than those previously found. However, one edge of the allowed region is nearly degenerate with an existing bound, which, using what is termed here as `the principle of noncoincident peril', is argued to be a strong indicator that the fine structure constant must be an environmental parameter in the multiverse.
We propose a new method to measure the tensor-to-scalar ratio $r$ using the circular polarization of the 21 cm radiation from the pre-reionization epoch. Our method relies on the splitting of the $F = 1$ hyperfine level of neutral hydrogen due to the quadrupole moment of the CMB. We show that unlike the Zeeman effect, where $M_{F}=\pm 1$ have opposite energy shifts, the CMB quadrupole shifts $M_{F}=\pm 1$ together relative to $M_{F}= 0$. This splitting leads to a small circular polarization of the emitted 21 cm radiation. In this paper (Paper I in a series on this effect), we present calculations on the microphysics behind this effect, accounting for all processes that affect the hyperfine transition. We conclude with an analytic formula for the circular polarization from the Dark Ages as a function of pre-reionization parameters and the value of the remote quadrupole of the CMB. We also calculate the splitting of the $F = 1$ hyperfine level due to other anisotropic radiation sources and show that they are not dominant. In a companion paper (Paper II) we make forecasts for measuring the tensor-to-scalar ratio $r$ using future radio arrays.
In the first paper of this series, we showed that the CMB quadrupole at high redshifts results in a small circular polarization of the emitted 21 cm radiation. In this paper we forecast the sensitivity of future radio experiments to measure the CMB quadrupole during the era of first cosmic light ($z\sim 20$). The tomographic measurement of 21 cm circular polarization allows us to construct a 3D remote quadrupole field. Measuring the $B$-mode component of this remote quadrupole field can be used to put bounds on the tensor-to-scalar ratio $r$. We make Fisher forecasts for a future Fast Fourier Transform Telescope (FFTT), consisting of an array of dipole antennas in a compact grid configuration, as a function of array size and observation time. We find that a FFTT with a side length of 100 km can achieve $\sigma(r)\sim 4\times 10^{-3}$ after ten years of observation and with a sky coverage $f_{\mathrm{sky}}\sim 0.7$. The forecasts are dependent on the evolution of the Lyman-$\alpha$ flux in the pre-reionization era, that remains observationally unconstrained. Finally, we calculate the typical order of magnitudes for circular polarization foregrounds and comment on their mitigation strategies. We conclude that detection of primordial gravitational waves with 21 cm observations is in principle possible, so long as the primordial magnetic field amplitude is small, but would require a very futuristic experiment with corresponding advances in calibration and foreground suppression techniques.
We detect and characterise extended, diffuse radio emission from galaxy clusters at 168 MHz within the Epoch of Reionization 0-hour field; a $45^\circ \times 45^\circ$ region of the southern sky centred on R. A. $= 0^\circ$, decl. $=-27^\circ$. We detect 31 diffuse radio sources; 3 of which are newly detected haloes in Abell 0141, Abell 2811, and Abell S1121; 2 newly detected relics in Abell 0033 and Abell 2751; 5 new halo candidates and a further 5 new relic candidates. Further, we detect a new phoenix candidate in Abell 2556 as well as 2 candidate dead radio galaxies at the centres of Abell 0122 and Abell S1136 likely associated with the brightest cluster galaxies. Beyond this we find 2 clusters with unclassifiable, diffuse steep-spectrum emission as well as a candidate double relic system associated with RXC J2351.0-1934. We present measured source properties such as their integrated flux densities, spectral indices, and sizes where possible. We find several of the diffuse sources to be ultra-steep including the halo in Abell 0141 which has $\alpha \leq -2.1 \pm 0.1$ making it one of the steepest halos known. Finally, we compare our sample of haloes with previously detected haloes and revisit established scaling relations of the radio halo power ($P_{1.4}$) with the cluster X-ray luminosity ($L_{\mathrm{X}}$) and mass ($M_{500}$). We find consistent fitting parameters for assumed power law relationships, and find that the $P_{1.4}-L_{\mathrm{X}}$ has less raw scatter than the corresponding $P_{1.4}-M_{500}$ despite inhomogeneous $L_{\mathrm{X}}$ measurements. These scaling relation properties are consistent with a sample of only non-cool core clusters hosting radio haloes.
An observational evidence for the intrinsic galaxy alignments in isolated spiral pairs is presented. From the catalog of the galaxy groups identified by Tempel et al. in the flux limited galaxy sample of the Sloan Digital Sky Survey Data Release 10, we select those groups consisting only of two spiral galaxies as isolated spiral pairs and investigate if and how strongly the spin axes of their two spiral members are aligned with each other. We detect a 4 sigma signal of intrinsic spin alignment in the isolated spiral pairs, which leads to the rejection of the null hypothesis at the 99.999 percent confidence level via the Kolmogorov Smirnov test. It is also found that those isolated pairs comprising two early-type spiral galaxies exhibit the strongest signal of intrinsic spin alignment while the weakest signal is found from the isolated pairs with one early type and one late type spiral galaxies. We also show that the strength of the alignment signal has a weak dependence on the angular separation distance as well as on the luminosity ratio of the member galaxies. Using the dark matter halos consisting of only two subhalos resolved in the EAGLE hydrodynamic simulations, we repeat the same analysis but fail to find any alignment tendency between the spin angular momentum vectors of the stellar components of the subhalos, which is in tension with the observational result. A couple of possible sources of this newly discovered local anomaly is discussed.
The correlation of weak lensing and Cosmic Microwave Anisotropy (CMB) data traces the pressure distribution of the hot, ionized gas and the underlying matter density field. The measured correlation is dominated by baryons residing in halos. Detecting the contribution from unbound gas by measuring the residual cross-correlation after masking all known halos requires a theoretical understanding of this correlation and its dependence with model parameters. Our model assumes that the gas in filaments is well described by a log-normal probability distribution function, with temperatures $10^{5-7}$K and overdensities $\xi\le 100$. The lensing-comptonization cross-correlation is dominated by gas with overdensities in the range $\xi\approx[3-33]$; the signal is generated at redshifts $z\le 1$. If only 10\% of the measured cross-correlation is due to unbound gas, then the most recent measurements set an upper limit of $\bar{T}_e\lesssim 10^6$K on the mean temperature of Inter Galactic Medium. The amplitude is proportional to the baryon fraction stored in filaments. The lensing-comptonization power spectrum peaks at a different scale than the gas in halos making it possible to distinguish both contributions. To trace the distribution of the low density and low temperature plasma on cosmological scales, the effect of halos will have to be subtracted from the data, requiring observations with larger signal-to-noise ratio than currently available.
Self-interacting dark matter has been proposed as a solution to small scale problems in cosmological structure formation, and hints of dark matter self scattering have been observed in mergers of galaxy clusters. One of the simplest models for self-interacting dark matter is a particle that is charged under dark electromagnetism, a new gauge interaction analogous to the usual electromagnetic force, but operating on the dark matter particle instead of the visible particles. In this case, the collisional behaviour of dark matter is primarily due to the formation of collisionless shocks, that should affect the distribution of DM in merging galaxy clusters. We evaluate the time and length scales of shock formation in cluster mergers, and discuss the implications for modelling charged dark matter in cosmological simulations.
We investigate quasar outflows at $z \geq 6$ by performing zoom-in cosmological hydrodynamical simulations. By employing the SPH code GADGET-3, we zoom in the $2 R_{200}$ region around a $2 \times 10^{12} M_{\odot}$ halo at $z = 6$, inside a $(500 ~ {\rm Mpc})^3$ comoving volume. We compare the results of our AGN runs with a control simulation in which only stellar/SN feedback is considered. Seeding $10^5 M_{\odot}$ BHs at the centers of $10^{9} M_{\odot}$ halos, we find the following results. BHs accrete gas at the Eddington rate over $z = 9 - 6$. At $z = 6$, our most-massive BH has grown to $M_{\rm BH} = 4 \times 10^9 M_{\odot}$. Fast ($v_{r} > 1000$ km/s), powerful ($\dot{M}_{\rm out} \sim 2000 M_{\odot}$/yr) outflows of shock-heated low-density gas form at $z \sim 7$, and propagate up to hundreds kpc. Star-formation is quenched over $z = 8 - 6$, and the total SFR (SFR surface density near the galaxy center) is reduced by a factor of $5$ ($1000$). We analyse the relative contribution of multiple physical process: (i) disrupting cosmic filamentary cold gas inflows, (ii) reducing central gas density, (iii) ejecting gas outside the galaxy; and find that AGN feedback has the following effects at $z = 6$. The inflowing gas mass fraction is reduced by $\sim 12 \%$, the high-density gas fraction is lowered by $\sim 13 \%$, and $\sim 20 \%$ of the gas outflows at a speed larger than the escape velocity ($500$ km/s). We conclude that quasar-host galaxies at $z \geq 6$ are accreting non-negligible amount of cosmic gas, nevertheless AGN feedback quenches their star formation dominantly by powerful outflows ejecting gas out of the host galaxy halo.
We introduce the first two simulations of the IllustrisTNG project, a next generation of cosmological magnetohydrodynamical simulations, focusing on the optical colors of galaxies. We explore TNG100, a rerun of the original Illustris box (~100 Mpc) and TNG300, which includes 2x2500^3 resolution elements in a volume twenty times larger (~300 Mpc); both are run using the new TNG model for galaxy formation. Here we present first results on the galaxy color bimodality at low redshift. Accounting for the attenuation of stellar light by dust, we compare the simulated (g-r) colors of 10^9 < M*/Msun < 10^12.5 galaxies to the observed distribution from the Sloan Digital Sky Survey (SDSS). We find a striking improvement with respect to the original Illustris simulation, as well as excellent quantitative agreement in comparison to the observations, with a sharp transition in median color from blue to red at a characteristic M* ~ 10^10.5 Msun. Investigating the build-up of the color-mass plane and the formation of the red sequence, we demonstrate that the primary driver of galaxy color transition in the TNG model is supermassive blackhole feedback in its low-accretion state. Across the entire population we measure a median color transition timescale dt_green of ~1.6 Gyr, a value which drops for increasingly massive galaxies. We find signatures of the physical process of quenching: at fixed stellar mass, the color of a galaxy correlates with its SFR, age, metallicity, and gas fraction, as well as the magnetic properties of both its interstellar medium and extended gaseous halo. Finally, we measure the amount of stellar mass growth on the red sequence. Galaxies with M* > 10^11 Msun which redden at z<1 accumulate on average ~25% of their final z=0 mass post-reddening; at the same time, ~18% of such massive galaxies acquire half or more of their final stellar mass while on the red sequence.
Hydrodynamical simulations of galaxy formation have now reached sufficient volume to make precision predictions for clustering on cosmologically relevant scales. Here we use our new IllustrisTNG simulations to study the non-linear correlation functions and power spectra of baryons, dark matter, galaxies and haloes over an exceptionally large range of scales. We find that baryonic effects increase the clustering of dark matter on small scales and damp the total matter power spectrum on scales up to k ~ 10 h/Mpc by 20%. The non-linear two-point correlation function of the stellar mass is close to a power-law over a wide range of scales and approximately invariant in time from very high redshift to the present. The two-point correlation function of the simulated galaxies agrees well with SDSS at its mean redshift z ~ 0.1, both as a function of stellar mass and when split according to galaxy colour, apart from a mild excess in the clustering of red galaxies in the stellar mass range 10^9-10^10 Msun/h^2. Given this agreement, the TNG simulations can make valuable theoretical predictions for the clustering bias of different galaxy samples. We find that the clustering length of the galaxy auto-correlation function depends strongly on stellar mass and redshift. Its power-law slope gamma is nearly invariant with stellar mass, but declines from gamma ~ 1.8 at redshift z=0 to gamma ~ 1.6 at redshift z ~ 1, beyond which the slope steepens again. We detect significant scale-dependencies in the bias of different observational tracers of large-scale structure, extending well into the range of the baryonic acoustic oscillations and causing nominal (yet fortunately correctable) shifts of the acoustic peaks of around ~5%.
We advocate the idea that the surprising emission of extreme ultra violet (EUV) radiation and soft x-rays from the Sun are powered externally by incident dark matter (DM) particles. The energy and the spectral shape of this otherwise unexpected solar irradiation is estimated within the quark nugget dark matter model. This model was originally invented as a natural explanation of the observed ratio $\Omega_{\rm dark} \sim \Omega_{\rm visible}$ when the DM and visible matter densities assume the same order of magnitude values. This generic consequence of the model is a result of the common origin of both types of matter which are formed during the same QCD transition and both proportional to the same fundamental dimensional parameter $\Lambda_{\rm QCD}$. We also present arguments suggesting that the transient brightening-like "nano-flares" in the Sun may be related to the annihilation events which inevitably occur in the solar atmosphere within this dark matter scenario.
The IllustrisTNG project is a new suite of cosmological magneto-hydrodynamical simulations of galaxy formation performed with the moving-mesh code Arepo and updated models for feedback physics. Here we introduce the first two simulations of the series, TNG100 and TNG300, and quantify the stellar mass content of about 4000 massive galaxy groups and clusters ($10^{13} \leq M_{\rm 200c}/M_{\rm sun} \leq 10^{15}$) at recent times ($z \leq 1$). We find that, while Milky Way-sized haloes contain about 95% of their diffuse stars within the inner tenth of their virial radii, this fraction drops to just 50% for the richest galaxy clusters. These clusters have half of their total stellar mass bound to satellite galaxies, with the other half being associated with the central galaxy and the diffuse intra-cluster light (ICL). The exact ICL fraction depends sensitively on the definition of a central galaxy's mass and varies in our most massive clusters between 20 to 40% of the total stellar mass. Haloes of $5\times 10^{14}M_{\rm sun}$ and above have more diffuse stellar mass outside 100 kpc than within 100 kpc, with power-law slopes of the radial density distribution as shallow as the dark matter's ( $-3.5 < \alpha_{\rm 3D} < -3$). We confirm that total halo mass is a very good predictor of stellar mass, and vice versa: the stellar mass measured within 30 kpc scales as $\propto (M_{\rm 500c})^{0.49}$ with a $\sim 0.12$ dex scatter. The 3D stellar half mass radii of massive galaxies fall too on a tight ($\sim$0.16 dex scatter) power-law relation with halo mass, $r^{\rm stars}_{\rm 0.5} \propto (M_{\rm 500c})^{0.53}$. Even more fundamentally, we find that halo mass alone is a good predictor for the whole stellar mass profiles of massive galaxies beyond the inner few kpc, and we show how on average these can be precisely recovered given a single mass measurement of the galaxy or its halo.
In the case of a metastable electroweak vacuum the quantum corrected effective potential plays a crucial role in the potential instability of the Standard Model. In the Early Universe, in particular during inflation and reheating, this instability can be triggered leading to catastrophic vacuum decay. We discuss how the large spacetime curvature of the Early Universe can be incorporated in the calculation and in many cases significantly modify the flat space prediction. The two key new elements are the unavoidable generation of the non-minimal coupling between the Higgs field and the scalar curvature of gravity and a curvature induced contribution to the running of the constants. For the minimal set up of the Standard Model and a decoupled inflation sector we show how a metastable vacuum can lead to very tight bounds for the non-minimal coupling. We also discuss a novel and very much related dark matter generation mechanism.
This contribution reviews scalar-tensor theories whose Lagrangian contains second-order derivatives of a scalar field but nevertheless propagate only one scalar mode (in addition to the usual two tensor modes), and are thus not plagued with the Ostrodradsky instability. These theories, which encompass the so-called Horndeski and Beyond Horndeski theories, have recently been fully classified up to cubic order in second-order derivatives. After introducing these theories, I present a few phenomenological aspects. In cosmology, these theories can be included in the unified effective description of dark energy and modified gravity. Finally, neutron star solutions in some specific models are discussed.
The large hierarchy between the Planck scale and the weak scale can be explained by the dynamical breaking of supersymmetry in strongly coupled gauge theories. Similarly, the hierarchy between the Planck scale and the energy scale of inflation may also originate from strong dynamics, which dynamically generate the inflaton potential. We present a model of the hidden sector which unifies these two ideas, i.e., in which the scales of inflation and supersymmetry breaking are provided by the dynamics of the same gauge group. The resultant inflation model is chaotic inflation with a fractional power-law potential in accord with the upper bound on the tensor-to-scalar ratio. The supersymmetry breaking scale can be much smaller than the inflation scale, so that the solution to the large hierarchy problem of the weak scale remains intact. As an intrinsic feature of our model, we find that the sgoldstino, which might disturb the inflationary dynamics, is automatically stabilized during inflation by dynamically generated corrections in the strongly coupled sector. This renders our model a field-theoretical realization of what is sometimes referred to as sgoldstino-less inflation.
Assuming both that our Universe is evolving into a de Sitter space and a vanishing cosmological constant, leaves only the option that the observed acceleration is provided by a "kinetic" energy of a scalar field. From an effective field theory point of view, the absence of Ostrogradsky instabilities restricts the choice to shift-symmetric Horndeski theories. Within these theories, we find the conditions for the existence of a de Sitter critical point in a universe filled by matter, radiation and a Horndeski scalar. Moreover, we show that this point is a universal attractor and we provide the tracking trajectory. Therefore, if a de Sitter fixed point exists within these models, our Universe will eventually evolve into a de Sitter space. As an example, we have discussed the case of the combined Galileon-Slotheon system, in which the Galileon is kinetically non-minimal coupled to the Einstein tensor. Interestingly, we have also found that the tracker trajectory of this system does not follow previous literature assumptions.
The offsets between the radial velocities of the rotational transitions of carbon monoxide and the fine structure transitions of neutral and singly ionized carbon are used to test the hypothetical variation of the fine structure constant, alpha. From the analysis of the [CI] and [CII] fine structure lines and low J rotational lines of 12CO and 13CO, emitted by the dark cloud L1599B in the Milky Way disk, we find no evidence for fractional changes in alpha at the level of |$\Delta \alpha/\alpha$| < 3*10^-7. For the neighbour galaxy M33 a stringent limit on Delta alpha/alpha is set from observations of three HII zones in [CII] and CO emission lines: |$\Delta \alpha/\alpha$| < 4*10^-7. Five systems over the redshift interval z = 5.7-6.4, showing CO J=6-5, J=7-6 and [CII] emission, yield a limit on |$\Delta \alpha/\alpha$| < 1.3*10^-5. Thus, a combination of the [CI], [CII], and CO emission lines turns out to be a powerful tool for probing the stability of the fundamental physical constants over a wide range of redshifts not accessible to optical spectral measurements.
Links to: arXiv, form interface, find, astro-ph, recent, 1707, contact, help (Access key information)
We introduce the IllustrisTNG project, a new suite of cosmological magneto-hydrodynamical simulations performed with the moving-mesh code AREPO employing an updated Illustris galaxy formation model. Here we focus on the general properties of magnetic fields and the diffuse radio emission in galaxy clusters. Magnetic fields are prevalent in galaxies, and their build-up is closely linked to structure formation. We find that structure formation amplifies the initial seed fields ($10^{-14}$ comoving Gauss) to the values observed in low-redshift galaxies ($1-10\,\mu G$). The magnetic field topology is closely connected to galaxy morphology such that irregular fields are hosted by early-type galaxies, while large-scale, ordered fields are present in disc galaxies. Using a simple model for the energy distribution of relativistic electrons we predict the diffuse radio emission of $280$ clusters with a baryonic mass resolution of $1.1\times 10^{7}\,{\rm M_{\odot}}$, and generate mock observations for VLA, LOFAR, ASKAP and SKA. Our simulated clusters show extended radio emission, whose detectability correlates with their virial mass. We reproduce the observed scaling relations between total radio power and X-ray emission, $M_{500}$, and the Sunyaev-Zel'dovich $Y_{\rm 500}$ parameter. The radio emission surface brightness profiles of our most massive clusters are in reasonable agreement with VLA measurements of Coma and Perseus. Finally, we discuss the fraction of detected extended radio haloes as a function of virial mass and source count functions for different instruments. Overall our results agree encouragingly well with observations, but a refined analysis requires a more sophisticated treatment of relativistic particles in large-scale galaxy formation simulations.
We examine the statistics of the low-redshift Ly-alpha forest in an adaptive mesh refinement hydrodynamic cosmological simulation of sufficient volume to include distinct large-scale environments. We compare our HI column density distribution of absorbers both with recent work and between two highly-refined regions of our simulation: a large-scale overdensity and a large-scale underdensity (on scales of approximately 20 Mpc). We recover the average results presented in Kollmeier et al. (2014) using different simulation methods. We further break down these results as a function of environment to examine the detailed dependence of absorber statistics on large-scale density. We find that the slope of the HI column density distribution in the 10$^{12.5}$ $\le$ N$_{HI}$/cm$^{-2}$ $\le$ 10$^{14.5}$ range depends on environment such that the slope becomes steeper for higher environmental density, and this difference reflects distinct physical conditions of the intergalactic medium on these scales. We track this difference to the different temperature structures of filaments in varying environments. Specifically, filaments in the overdensity are hotter and, correspondingly, are composed of gas with lower HI fractions than those in underdense environments. Our results highlight that in order to understand the physics driving the HI CDD, we need not only improved accounting of the sources of ionizing UV photons, but also of the physical conditions of the IGM and how this may vary as a function of large-scale environment.
Visinelli and Gondolo (2015, hereafter VG15) derived analytic expressions for the evolution of the dark matter temperature in a generic cosmological model. They then calculated the dark matter kinetic decoupling temperature $T_{\mathrm{kd}}$ and compared their results to the Gelmini and Gondolo (2008, hereafter GG08) calculation of $T_{\mathrm{kd}}$ in an early matter-dominated era (EMDE), which occurs when the Universe is dominated by either a decaying oscillating scalar field or a semistable massive particle before Big Bang nucleosynthesis. VG15 found that dark matter decouples at a lower temperature in an EMDE than it would in a radiation-dominated era, while GG08 found that dark matter decouples at a higher temperature in an EMDE than it would in a radiation-dominated era. VG15 attributed this discrepancy to the presence of a matching constant that ensures that the dark matter temperature is continuous during the transition from the EMDE to the subsequent radiation-dominated era and concluded that the GG08 result is incorrect. We show that the disparity is due to the fact that VG15 compared $T_\mathrm{kd}$ in an EMDE to the decoupling temperature in a radiation-dominated universe that would result in the same dark matter temperature at late times. Since decoupling during an EMDE leaves the dark matter colder than it would be if it decoupled during radiation domination, this temperature is much higher than $T_\mathrm{kd}$ in a standard thermal history, which is indeed lower than $T_{\mathrm{kd}}$ in an EMDE, as stated by GG08.
We show that a subdominant component of dissipative dark matter resembling the Standard Model can form many intermediate-mass black hole seeds during the first structure formation epoch, such that a few percent of all the dwarf galaxies host one. We also observe that, in the presence of this matter sector, the black holes will grow at a much faster rate with respect to the ordinary case. These facts can explain the observed abundance of supermassive black holes feeding high-redshift quasars. The scenario will have interesting observational consequences, such as dark substructures and gravitational wave production.
To derive a power spectrum for energy density inhomogeneities in a closed universe, we study a spatially-closed inflation model. The constants of integration in the inflation epoch solutions are determined from Hawking's prescription that the quantum state of the universe only include field configurations that are regular on the Euclidean section. The power spectrum we find is not a power law, and depends on spatial wavenumber in the way expected for a generalization to the closed model of the standard flat-space scale-invariant power spectrum. (Abridged.)
We study Planck 2015 cosmic microwave background (CMB) anisotropy data using the energy density inhomogeneity power spectrum generated by quantum fluctuations during an early epoch of inflation in the non-flat $\Lambda$CDM model. Unlike earlier analyses of non-flat models, which assumed an inconsistent power-law power spectrum of energy density inhomogeneities, we find that the Planck 2015 data alone, and also in conjunction with baryon acoustic oscillation measurements, are reasonably well fit by a closed $\Lambda$CDM model in which spatial curvature contributes a few percent of the current cosmological energy density budget. In this model, the measured Hubble constant and non-relativistic matter density parameter are in good agreement with values determined using most other data. Depending on parameter values, the closed $\Lambda$CDM model has reduced power, relative to the tilted, spatially-flat $\Lambda$CDM case, and can partially alleviate the low multipole CMB temperature anisotropy deficit and can help partially reconcile the CMB anisotropy and weak lensing $\sigma_8$ constraints, at the expense of slightly worsening the fit to higher multipole CMB temperature anisotropy data. Our results are interesting but tentative; a more thorough analysis is needed to properly gauge their significance.
Results of recent experiment reinstate feasibility to the hypothesis that biomolecular homochirality originates from beta decay. Coupled with hints that this process occurred extraterrestrially suggests aluminum-26 as the most likely source. If true, then its appropriateness is highly dependent on the half-life and energy of this decay. Demanding that this mechanism hold places new constraints on the anthropically allowed range for multiple parameters, including the electron mass, difference between up and down quark masses, the fine structure constant, and the electroweak scale. These new constraints on particle masses are tighter than those previously found. However, one edge of the allowed region is nearly degenerate with an existing bound, which, using what is termed here as `the principle of noncoincident peril', is argued to be a strong indicator that the fine structure constant must be an environmental parameter in the multiverse.
We propose a new method to measure the tensor-to-scalar ratio $r$ using the circular polarization of the 21 cm radiation from the pre-reionization epoch. Our method relies on the splitting of the $F = 1$ hyperfine level of neutral hydrogen due to the quadrupole moment of the CMB. We show that unlike the Zeeman effect, where $M_{F}=\pm 1$ have opposite energy shifts, the CMB quadrupole shifts $M_{F}=\pm 1$ together relative to $M_{F}= 0$. This splitting leads to a small circular polarization of the emitted 21 cm radiation. In this paper (Paper I in a series on this effect), we present calculations on the microphysics behind this effect, accounting for all processes that affect the hyperfine transition. We conclude with an analytic formula for the circular polarization from the Dark Ages as a function of pre-reionization parameters and the value of the remote quadrupole of the CMB. We also calculate the splitting of the $F = 1$ hyperfine level due to other anisotropic radiation sources and show that they are not dominant. In a companion paper (Paper II) we make forecasts for measuring the tensor-to-scalar ratio $r$ using future radio arrays.
In the first paper of this series, we showed that the CMB quadrupole at high redshifts results in a small circular polarization of the emitted 21 cm radiation. In this paper we forecast the sensitivity of future radio experiments to measure the CMB quadrupole during the era of first cosmic light ($z\sim 20$). The tomographic measurement of 21 cm circular polarization allows us to construct a 3D remote quadrupole field. Measuring the $B$-mode component of this remote quadrupole field can be used to put bounds on the tensor-to-scalar ratio $r$. We make Fisher forecasts for a future Fast Fourier Transform Telescope (FFTT), consisting of an array of dipole antennas in a compact grid configuration, as a function of array size and observation time. We find that a FFTT with a side length of 100 km can achieve $\sigma(r)\sim 4\times 10^{-3}$ after ten years of observation and with a sky coverage $f_{\mathrm{sky}}\sim 0.7$. The forecasts are dependent on the evolution of the Lyman-$\alpha$ flux in the pre-reionization era, that remains observationally unconstrained. Finally, we calculate the typical order of magnitudes for circular polarization foregrounds and comment on their mitigation strategies. We conclude that detection of primordial gravitational waves with 21 cm observations is in principle possible, so long as the primordial magnetic field amplitude is small, but would require a very futuristic experiment with corresponding advances in calibration and foreground suppression techniques.
We detect and characterise extended, diffuse radio emission from galaxy clusters at 168 MHz within the Epoch of Reionization 0-hour field; a $45^\circ \times 45^\circ$ region of the southern sky centred on R. A. $= 0^\circ$, decl. $=-27^\circ$. We detect 31 diffuse radio sources; 3 of which are newly detected haloes in Abell 0141, Abell 2811, and Abell S1121; 2 newly detected relics in Abell 0033 and Abell 2751; 5 new halo candidates and a further 5 new relic candidates. Further, we detect a new phoenix candidate in Abell 2556 as well as 2 candidate dead radio galaxies at the centres of Abell 0122 and Abell S1136 likely associated with the brightest cluster galaxies. Beyond this we find 2 clusters with unclassifiable, diffuse steep-spectrum emission as well as a candidate double relic system associated with RXC J2351.0-1934. We present measured source properties such as their integrated flux densities, spectral indices, and sizes where possible. We find several of the diffuse sources to be ultra-steep including the halo in Abell 0141 which has $\alpha \leq -2.1 \pm 0.1$ making it one of the steepest halos known. Finally, we compare our sample of haloes with previously detected haloes and revisit established scaling relations of the radio halo power ($P_{1.4}$) with the cluster X-ray luminosity ($L_{\mathrm{X}}$) and mass ($M_{500}$). We find consistent fitting parameters for assumed power law relationships, and find that the $P_{1.4}-L_{\mathrm{X}}$ has less raw scatter than the corresponding $P_{1.4}-M_{500}$ despite inhomogeneous $L_{\mathrm{X}}$ measurements. These scaling relation properties are consistent with a sample of only non-cool core clusters hosting radio haloes.
An observational evidence for the intrinsic galaxy alignments in isolated spiral pairs is presented. From the catalog of the galaxy groups identified by Tempel et al. in the flux limited galaxy sample of the Sloan Digital Sky Survey Data Release 10, we select those groups consisting only of two spiral galaxies as isolated spiral pairs and investigate if and how strongly the spin axes of their two spiral members are aligned with each other. We detect a 4 sigma signal of intrinsic spin alignment in the isolated spiral pairs, which leads to the rejection of the null hypothesis at the 99.999 percent confidence level via the Kolmogorov Smirnov test. It is also found that those isolated pairs comprising two early-type spiral galaxies exhibit the strongest signal of intrinsic spin alignment while the weakest signal is found from the isolated pairs with one early type and one late type spiral galaxies. We also show that the strength of the alignment signal has a weak dependence on the angular separation distance as well as on the luminosity ratio of the member galaxies. Using the dark matter halos consisting of only two subhalos resolved in the EAGLE hydrodynamic simulations, we repeat the same analysis but fail to find any alignment tendency between the spin angular momentum vectors of the stellar components of the subhalos, which is in tension with the observational result. A couple of possible sources of this newly discovered local anomaly is discussed.
The correlation of weak lensing and Cosmic Microwave Anisotropy (CMB) data traces the pressure distribution of the hot, ionized gas and the underlying matter density field. The measured correlation is dominated by baryons residing in halos. Detecting the contribution from unbound gas by measuring the residual cross-correlation after masking all known halos requires a theoretical understanding of this correlation and its dependence with model parameters. Our model assumes that the gas in filaments is well described by a log-normal probability distribution function, with temperatures $10^{5-7}$K and overdensities $\xi\le 100$. The lensing-comptonization cross-correlation is dominated by gas with overdensities in the range $\xi\approx[3-33]$; the signal is generated at redshifts $z\le 1$. If only 10\% of the measured cross-correlation is due to unbound gas, then the most recent measurements set an upper limit of $\bar{T}_e\lesssim 10^6$K on the mean temperature of Inter Galactic Medium. The amplitude is proportional to the baryon fraction stored in filaments. The lensing-comptonization power spectrum peaks at a different scale than the gas in halos making it possible to distinguish both contributions. To trace the distribution of the low density and low temperature plasma on cosmological scales, the effect of halos will have to be subtracted from the data, requiring observations with larger signal-to-noise ratio than currently available.
Self-interacting dark matter has been proposed as a solution to small scale problems in cosmological structure formation, and hints of dark matter self scattering have been observed in mergers of galaxy clusters. One of the simplest models for self-interacting dark matter is a particle that is charged under dark electromagnetism, a new gauge interaction analogous to the usual electromagnetic force, but operating on the dark matter particle instead of the visible particles. In this case, the collisional behaviour of dark matter is primarily due to the formation of collisionless shocks, that should affect the distribution of DM in merging galaxy clusters. We evaluate the time and length scales of shock formation in cluster mergers, and discuss the implications for modelling charged dark matter in cosmological simulations.
We investigate quasar outflows at $z \geq 6$ by performing zoom-in cosmological hydrodynamical simulations. By employing the SPH code GADGET-3, we zoom in the $2 R_{200}$ region around a $2 \times 10^{12} M_{\odot}$ halo at $z = 6$, inside a $(500 ~ {\rm Mpc})^3$ comoving volume. We compare the results of our AGN runs with a control simulation in which only stellar/SN feedback is considered. Seeding $10^5 M_{\odot}$ BHs at the centers of $10^{9} M_{\odot}$ halos, we find the following results. BHs accrete gas at the Eddington rate over $z = 9 - 6$. At $z = 6$, our most-massive BH has grown to $M_{\rm BH} = 4 \times 10^9 M_{\odot}$. Fast ($v_{r} > 1000$ km/s), powerful ($\dot{M}_{\rm out} \sim 2000 M_{\odot}$/yr) outflows of shock-heated low-density gas form at $z \sim 7$, and propagate up to hundreds kpc. Star-formation is quenched over $z = 8 - 6$, and the total SFR (SFR surface density near the galaxy center) is reduced by a factor of $5$ ($1000$). We analyse the relative contribution of multiple physical process: (i) disrupting cosmic filamentary cold gas inflows, (ii) reducing central gas density, (iii) ejecting gas outside the galaxy; and find that AGN feedback has the following effects at $z = 6$. The inflowing gas mass fraction is reduced by $\sim 12 \%$, the high-density gas fraction is lowered by $\sim 13 \%$, and $\sim 20 \%$ of the gas outflows at a speed larger than the escape velocity ($500$ km/s). We conclude that quasar-host galaxies at $z \geq 6$ are accreting non-negligible amount of cosmic gas, nevertheless AGN feedback quenches their star formation dominantly by powerful outflows ejecting gas out of the host galaxy halo.
We introduce the first two simulations of the IllustrisTNG project, a next generation of cosmological magnetohydrodynamical simulations, focusing on the optical colors of galaxies. We explore TNG100, a rerun of the original Illustris box (~100 Mpc) and TNG300, which includes 2x2500^3 resolution elements in a volume twenty times larger (~300 Mpc); both are run using the new TNG model for galaxy formation. Here we present first results on the galaxy color bimodality at low redshift. Accounting for the attenuation of stellar light by dust, we compare the simulated (g-r) colors of 10^9 < M*/Msun < 10^12.5 galaxies to the observed distribution from the Sloan Digital Sky Survey (SDSS). We find a striking improvement with respect to the original Illustris simulation, as well as excellent quantitative agreement in comparison to the observations, with a sharp transition in median color from blue to red at a characteristic M* ~ 10^10.5 Msun. Investigating the build-up of the color-mass plane and the formation of the red sequence, we demonstrate that the primary driver of galaxy color transition in the TNG model is supermassive blackhole feedback in its low-accretion state. Across the entire population we measure a median color transition timescale dt_green of ~1.6 Gyr, a value which drops for increasingly massive galaxies. We find signatures of the physical process of quenching: at fixed stellar mass, the color of a galaxy correlates with its SFR, age, metallicity, and gas fraction, as well as the magnetic properties of both its interstellar medium and extended gaseous halo. Finally, we measure the amount of stellar mass growth on the red sequence. Galaxies with M* > 10^11 Msun which redden at z<1 accumulate on average ~25% of their final z=0 mass post-reddening; at the same time, ~18% of such massive galaxies acquire half or more of their final stellar mass while on the red sequence.
Hydrodynamical simulations of galaxy formation have now reached sufficient volume to make precision predictions for clustering on cosmologically relevant scales. Here we use our new IllustrisTNG simulations to study the non-linear correlation functions and power spectra of baryons, dark matter, galaxies and haloes over an exceptionally large range of scales. We find that baryonic effects increase the clustering of dark matter on small scales and damp the total matter power spectrum on scales up to k ~ 10 h/Mpc by 20%. The non-linear two-point correlation function of the stellar mass is close to a power-law over a wide range of scales and approximately invariant in time from very high redshift to the present. The two-point correlation function of the simulated galaxies agrees well with SDSS at its mean redshift z ~ 0.1, both as a function of stellar mass and when split according to galaxy colour, apart from a mild excess in the clustering of red galaxies in the stellar mass range 10^9-10^10 Msun/h^2. Given this agreement, the TNG simulations can make valuable theoretical predictions for the clustering bias of different galaxy samples. We find that the clustering length of the galaxy auto-correlation function depends strongly on stellar mass and redshift. Its power-law slope gamma is nearly invariant with stellar mass, but declines from gamma ~ 1.8 at redshift z=0 to gamma ~ 1.6 at redshift z ~ 1, beyond which the slope steepens again. We detect significant scale-dependencies in the bias of different observational tracers of large-scale structure, extending well into the range of the baryonic acoustic oscillations and causing nominal (yet fortunately correctable) shifts of the acoustic peaks of around ~5%.
We advocate the idea that the surprising emission of extreme ultra violet (EUV) radiation and soft x-rays from the Sun are powered externally by incident dark matter (DM) particles. The energy and the spectral shape of this otherwise unexpected solar irradiation is estimated within the quark nugget dark matter model. This model was originally invented as a natural explanation of the observed ratio $\Omega_{\rm dark} \sim \Omega_{\rm visible}$ when the DM and visible matter densities assume the same order of magnitude values. This generic consequence of the model is a result of the common origin of both types of matter which are formed during the same QCD transition and both proportional to the same fundamental dimensional parameter $\Lambda_{\rm QCD}$. We also present arguments suggesting that the transient brightening-like "nano-flares" in the Sun may be related to the annihilation events which inevitably occur in the solar atmosphere within this dark matter scenario.
The IllustrisTNG project is a new suite of cosmological magneto-hydrodynamical simulations of galaxy formation performed with the moving-mesh code Arepo and updated models for feedback physics. Here we introduce the first two simulations of the series, TNG100 and TNG300, and quantify the stellar mass content of about 4000 massive galaxy groups and clusters ($10^{13} \leq M_{\rm 200c}/M_{\rm sun} \leq 10^{15}$) at recent times ($z \leq 1$). We find that, while Milky Way-sized haloes contain about 95% of their diffuse stars within the inner tenth of their virial radii, this fraction drops to just 50% for the richest galaxy clusters. These clusters have half of their total stellar mass bound to satellite galaxies, with the other half being associated with the central galaxy and the diffuse intra-cluster light (ICL). The exact ICL fraction depends sensitively on the definition of a central galaxy's mass and varies in our most massive clusters between 20 to 40% of the total stellar mass. Haloes of $5\times 10^{14}M_{\rm sun}$ and above have more diffuse stellar mass outside 100 kpc than within 100 kpc, with power-law slopes of the radial density distribution as shallow as the dark matter's ( $-3.5 < \alpha_{\rm 3D} < -3$). We confirm that total halo mass is a very good predictor of stellar mass, and vice versa: the stellar mass measured within 30 kpc scales as $\propto (M_{\rm 500c})^{0.49}$ with a $\sim 0.12$ dex scatter. The 3D stellar half mass radii of massive galaxies fall too on a tight ($\sim$0.16 dex scatter) power-law relation with halo mass, $r^{\rm stars}_{\rm 0.5} \propto (M_{\rm 500c})^{0.53}$. Even more fundamentally, we find that halo mass alone is a good predictor for the whole stellar mass profiles of massive galaxies beyond the inner few kpc, and we show how on average these can be precisely recovered given a single mass measurement of the galaxy or its halo.
In the case of a metastable electroweak vacuum the quantum corrected effective potential plays a crucial role in the potential instability of the Standard Model. In the Early Universe, in particular during inflation and reheating, this instability can be triggered leading to catastrophic vacuum decay. We discuss how the large spacetime curvature of the Early Universe can be incorporated in the calculation and in many cases significantly modify the flat space prediction. The two key new elements are the unavoidable generation of the non-minimal coupling between the Higgs field and the scalar curvature of gravity and a curvature induced contribution to the running of the constants. For the minimal set up of the Standard Model and a decoupled inflation sector we show how a metastable vacuum can lead to very tight bounds for the non-minimal coupling. We also discuss a novel and very much related dark matter generation mechanism.
This contribution reviews scalar-tensor theories whose Lagrangian contains second-order derivatives of a scalar field but nevertheless propagate only one scalar mode (in addition to the usual two tensor modes), and are thus not plagued with the Ostrodradsky instability. These theories, which encompass the so-called Horndeski and Beyond Horndeski theories, have recently been fully classified up to cubic order in second-order derivatives. After introducing these theories, I present a few phenomenological aspects. In cosmology, these theories can be included in the unified effective description of dark energy and modified gravity. Finally, neutron star solutions in some specific models are discussed.
The large hierarchy between the Planck scale and the weak scale can be explained by the dynamical breaking of supersymmetry in strongly coupled gauge theories. Similarly, the hierarchy between the Planck scale and the energy scale of inflation may also originate from strong dynamics, which dynamically generate the inflaton potential. We present a model of the hidden sector which unifies these two ideas, i.e., in which the scales of inflation and supersymmetry breaking are provided by the dynamics of the same gauge group. The resultant inflation model is chaotic inflation with a fractional power-law potential in accord with the upper bound on the tensor-to-scalar ratio. The supersymmetry breaking scale can be much smaller than the inflation scale, so that the solution to the large hierarchy problem of the weak scale remains intact. As an intrinsic feature of our model, we find that the sgoldstino, which might disturb the inflationary dynamics, is automatically stabilized during inflation by dynamically generated corrections in the strongly coupled sector. This renders our model a field-theoretical realization of what is sometimes referred to as sgoldstino-less inflation.
Assuming both that our Universe is evolving into a de Sitter space and a vanishing cosmological constant, leaves only the option that the observed acceleration is provided by a "kinetic" energy of a scalar field. From an effective field theory point of view, the absence of Ostrogradsky instabilities restricts the choice to shift-symmetric Horndeski theories. Within these theories, we find the conditions for the existence of a de Sitter critical point in a universe filled by matter, radiation and a Horndeski scalar. Moreover, we show that this point is a universal attractor and we provide the tracking trajectory. Therefore, if a de Sitter fixed point exists within these models, our Universe will eventually evolve into a de Sitter space. As an example, we have discussed the case of the combined Galileon-Slotheon system, in which the Galileon is kinetically non-minimal coupled to the Einstein tensor. Interestingly, we have also found that the tracker trajectory of this system does not follow previous literature assumptions.
The offsets between the radial velocities of the rotational transitions of carbon monoxide and the fine structure transitions of neutral and singly ionized carbon are used to test the hypothetical variation of the fine structure constant, alpha. From the analysis of the [CI] and [CII] fine structure lines and low J rotational lines of 12CO and 13CO, emitted by the dark cloud L1599B in the Milky Way disk, we find no evidence for fractional changes in alpha at the level of |$\Delta \alpha/\alpha$| < 3*10^-7. For the neighbour galaxy M33 a stringent limit on Delta alpha/alpha is set from observations of three HII zones in [CII] and CO emission lines: |$\Delta \alpha/\alpha$| < 4*10^-7. Five systems over the redshift interval z = 5.7-6.4, showing CO J=6-5, J=7-6 and [CII] emission, yield a limit on |$\Delta \alpha/\alpha$| < 1.3*10^-5. Thus, a combination of the [CI], [CII], and CO emission lines turns out to be a powerful tool for probing the stability of the fundamental physical constants over a wide range of redshifts not accessible to optical spectral measurements.
Links to: arXiv, form interface, find, astro-ph, recent, 1707, contact, help (Access key information)