We provide an analytical description of the line broadening of HI absorbers in the Lyman-alpha forest resulting from Doppler broadening and Jeans smoothing. We demonstrate that our relation captures the dependence of the line-width on column density for narrow lines in z~3 mock spectra remarkably well. Broad lines at a given column density arise when the underlying density structure is more complex, and such clustering is not captured by our model. Our understanding of the line broadening opens the way to a new method to characterise the thermal state of the intergalactic medium and to determine the sizes of the absorbing structures.
The observational status of inflation after the Planck 2013 and 2015 results and the BICEP2/Keck Array and Planck joint analysis is discussed. These pedagogical lecture notes are intended to serve as a technical guide filling the gap between the theoretical articles on inflation and the experimental works on astrophysical and cosmological data. After a short discussion of the central tenets at the basis of inflation (negative self-gravitating pressure) and its experimental verifications, it reviews how the most recent Cosmic Microwave Background (CMB) anisotropy measurements constrain cosmic inflation. The fact that vanilla inflationary models are, so far, preferred by the observations is discussed and the reason why plateau-like potential versions of inflation are favored within this subclass of scenarios is explained. Finally, how well the future measurements, in particular of $B$-Mode CMB polarization or primordial gravity waves, will help to improve our knowledge about inflation is also investigated.
We compare the efficiency with which 2D and 3D weak lensing mass mapping techniques are able to detect clusters of galaxies using two state-of-the-art mass reconstruction techniques: MRLens in 2D and GLIMPSE in 3D. We simulate otherwise-empty cluster fields for 96 different virial mass-redshift combinations spanning the ranges $3\times10^{13}h^{-1}M_\odot \le M_{vir}\le 10^{15}h^{-1}M_\odot$ and $0.05 \le z_{\rm cl} \le 0.75$, and for each generate 1000 realisations of noisy shear data in 2D and 3D. For each field, we then compute the cluster (false) detection rate as the mean number of cluster (false) detections per reconstruction over the sample of 1000 reconstructions. We show that both MRLens and GLIMPSE are effective tools for the detection of clusters from weak lensing measurements, and provide comparable quality reconstructions at low redshift. At high redshift, GLIMPSE reconstructions offer increased sensitivity in the detection of clusters, yielding cluster detection rates up to a factor of $\sim 10\times$ that seen in 2D reconstructions using MRLens. We conclude that 3D mass mapping techniques are more efficient for the detection of clusters of galaxies in weak lensing surveys than 2D methods, particularly since 3D reconstructions yield unbiased estimators of both the mass and redshift of the detected clusters directly.
During this work an interacting chameleon-like scalar field scenario, by considering SNeIa, CMB, BAO and OHD data sets is investigated. Some cosmological parameters includes of Hubble, deceleration and coincidence parameters in such mechanism are analysed. It is realized for estimation the free parameters of a theoretical model, it is better all mentioned observational data sets be considered. In fact if one considers SNeIa, CMB and BAO but ignores OHD it maybe leads to an incorrect result. Also it will find out, when we margin the free parameters, the $\chi _{\rm{T}}^2$ function should be re-weighted, this fact arises from the abundance of SNeIa and OHD sources in comparison to CMB and BAO data sets. We margin the likelihood $\mathcal{L} (\Omega_{\rm{m0}} ,\omega_1 , \beta)$ with respect to $\omega_1$, $\beta$ and $\Omega_{\rm{m0}}$ respectively and by means of two dimensional confidence levels $68.3\%$, $90\%$ and $95.4\%$, the relative probability functions are plotted. Also the quantities which maximize the marginalized likelihoods using mentioned confidence levels are obtained. In addition, $\chi _{\rm{T}}^2 = \chi _{{\rm{SNe}}}^2 + \chi _{{\rm{OHD}}}^2 + 3 \, \chi _{{\rm{CMB}}}^2 + 3 \, \chi _{{\rm{BAO}}}^2$, we margin the relative likelihood functions in one dimension, and based on these calculations the best fitted free parameters of the model will be obtained.
We present foreground-reduced CMB maps derived from the full Planck data set in both temperature and polarization. Compared to the corresponding Planck 2013 temperature sky maps, the total data volume is larger by a factor of 3.2 for frequencies between 30 and 70 GHz, and by 1.9 for frequencies between 100 and 857 GHz. In addition, systematic errors in the forms of temperature-to-polarization leakage, analogue-to-digital conversion uncertainties, and very long time constant errors have been dramatically reduced, to the extent that the cosmological polarization signal may now be robustly recovered on angular scales $\ell\gtrsim40$. On the very largest scales, instrumental systematic residuals are still non-negligible compared to the expected cosmological signal, and modes with $\ell < 20$ are accordingly suppressed in the current polarization maps by high-pass filtering. As in 2013, four different CMB component separation algorithms are applied to these observations, providing a measure of stability with respect to algorithmic and modelling choices. The resulting polarization maps have rms instrumental noise ranging between 0.21 and 0.27$\,\mu\textrm{K}$ averaged over 55 arcmin pixels, and between 4.5 and 6.1$\,\mu\textrm{K}$ averaged over 3.4 arcmin pixels. The cosmological parameters derived from the analysis of temperature power spectra are in agreement at the $1\sigma$ level with the Planck 2015 likelihood. Unresolved mismatches between the noise properties of the data and simulations prevent a satisfactory description of the higher-order statistical properties of the polarization maps. Thus, the primary applications of these polarization maps are those that do not require massive simulations for accurate estimation of uncertainties, for instance estimation of cross-spectra and cross-correlations, or stacking analyses.
We study the phenomenology of the CMSSM/mSUGRA with non-thermal neutralino dark matter. Besides the standard parameters of the CMSSM we include the reheating temperature as an extra parameter. Imposing radiative electroweak symmetry breaking with a Higgs mass around 125 GeV and no dark matter overproduction, we contrast the scenario with different experimental bounds from colliders (LEP, LHC), cosmic microwave background (Planck), direct (LUX, XENON100, CDMS, IceCube) and indirect (Fermi) dark matter searches. The allowed parameter space is characterised by a Higgsino-like LSP with a mass around 300 GeV. The observed dark matter abundance can be saturated for reheating temperatures around 2 GeV while larger temperatures require extra non-neutralino dark matter candidates and extend the allowed parameter space. Sfermion and gluino masses are in the few TeV region. These scenarios can be achieved in string models of sequestered supersymmetry breaking which avoid cosmological moduli problems and are compatible with gauge coupling unification. Astrophysics and particle physics experiments will fully investigate this non-thermal scenario in the near future.
Massive black hole binaries are the primary source of gravitational waves (GW) for the future eLISA observatory. The detection and parameter estimation of these sources to high redshift would provide invaluable information on the formation mechanisms of seed black holes, and on the evolution of massive black holes and their host galaxies through cosmic time. The Fisher information matrix has been the standard tool for GW parameter estimation in the last two decades. However, recent studies have questioned the validity of using the Fisher matrix approach. For example, the Fisher matrix approach sometimes predicts errors of $\geq100\%$ in the estimation of parameters such as the luminosity distance and sky position. With advances in computing power, Bayesian inference is beginning to replace the Fisher matrix approximation in parameter estimation studies. In this work, we conduct a Bayesian inference analysis for 120 sources situated at redshifts of between $0.1\leq z\leq 13.2$, and compare the results with those from a Fisher matrix analysis. The Fisher matrix results suggest that for this particular selection of sources, eLISA would be unable to localize sources at redshifts of $z\lesssim6$. In contrast, Bayesian inference provides finite error estimations for all sources in the study, and shows that we can establish minimum closest distances for all sources. The study further predicts that we should be capable with eLISA, out to a redshift of at least $z\leq13$, of predicting a maximum error in the chirp mass of $\lesssim 1\%$, the reduced mass of $\lesssim20\%$, the time to coalescence of 2 hours, and to a redshift of $z\sim5$, the inclination of the source with a maximum error of $\sim60$ degrees.
Being able to reliably track perturbations across bounces and turnarounds in cyclic and bouncing cosmology lies at the heart of being able to compare the predictions of these models with the Cosmic Microwave Background observations. This has been a challenging task due to the unknown nature of the physics involved during the bounce as well as the technical challenge of matching perturbations precisely between the expansion and contraction phases. In this paper, we will present general techniques (analytical and numerical) that can be applied to understand the physics of the fluctuations, especially those with "long" wavelengths, and test its validity in some simple bouncing/cyclic toy models where the physics is well understood. We will then apply our techniques to more interesting cosmological models such as the bounce inflation and cyclic inflation.
We consider Horndeski cosmological models, with a minisuperspace Lagrangian linear in the field derivative, that are able to screen any vacuum energy and material content leading to a spatially flat de Sitter vacuum fixed by the theory itself. Furthermore, we investigate particular models with a cosmic evolution independent of the material content and use them to understand the general characteristics of this framework. We also consider more realistic models, which we denote the "term-by-term" and "tripod" models, focusing attention on cases in which the critical point is indeed an attractor solution and the cosmological history is of particular interest.
We point out that the (pseudo-)conformal Universe scenario may be realized as decay of conformally invariant, metastable vacuum, which proceeds via spontaneous nucleation and subsequent growth of a bubble of a putative new phase. We study perturbations about the bubble and show that their leading late-time properties coincide with those inherent in the original models with homogeneously rolling backgrounds. In particular, the perturbations of a spectator dimension-zero field have flat power spectrum.
Dark energy is often assumed to be composed by a single scalar field. The background cosmic expansion is not sufficient to determine whether this is true or not. We study multi-field scalar-tensor models with a general dark matter source and write the observable modified gravity parameters (effective gravitational constant and anisotropic stress) in the form of a ratio of polynomials in the Fourier wavenumber k of order 2N, where N is the number of scalar fields. By comparing these observables to real data it is in principle possible to determine the number of dark energy scalar fields coupled to gravity. We also show that there are no realistic non-trivial cases in which the order of the polynomials is reduced.
Links to: arXiv, form interface, find, astro-ph, recent, 1502, contact, help (Access key information)
In this paper, we report new limits on 21cm emission from cosmic reionization based on a 135-day observing campaign with a 64-element deployment of the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER) in South Africa. This work extends the work presented in Parsons et al. (2014) with more collecting area, a longer observing period, improved redundancy-based calibration, optimal fringe-rate filtering, and improved power-spectral analysis using optimal quadratic estimators. The result is a new $2\sigma$ upper limit on $\Delta^{2}(k)$ of (22.4 mK)$^2$ in the range $0.15 < k < 0.5h\ {\rm Mpc}^{-1}$ at $z = 8.4$. This represents a three-fold improvement over the previous best upper limit. As we discuss in more depth in a forthcoming paper (Pober et al. 2015, in prep), this upper limit supports and extends previous evidence against extremely cold reionization scenarios. We conclude with a discussion of implications for future 21cm reionization experiments, including the newly funded Hydrogen Epoch of Reionization Array (HERA).
This is the first in a series of papers studying the astrophysics and cosmology of massive, dynamically relaxed galaxy clusters. Here we present a new, automated method for identifying relaxed clusters based on their morphologies in X-ray imaging data. While broadly similar to others in the literature, the morphological quantities that we measure are specifically designed to provide a fair basis for comparison across a range of data quality and cluster redshifts, to be robust against missing data due to point-source masks and gaps between detectors, and to avoid strong assumptions about the cosmological background and cluster masses. Based on three morphological indicators - Symmetry, Peakiness and Alignment - we develop the SPA criterion for relaxation. This analysis was applied to a large sample of cluster observations from the Chandra and ROSAT archives. Of the 361 clusters which received the SPA treatment, 57 (16 per cent) were subsequently found to be relaxed according to our criterion. We compare our measurements to similar estimators in the literature, as well as projected ellipticity and other image measures, and comment on trends in the relaxed cluster fraction with redshift, temperature, and survey selection method. Code implementing our morphological analysis will be made available on the web.
We have observed damped longitudinal cosmological-scale oscillations in a unique model-independent plot of scale factor against lookback time for Type Ia supernovae data. We found several first-derivative relative maxima/minima spanning the range of reported transition-redshifts. These extrema comprise 2 full cycles with a period of approximately 0.15 Hubble times (H0=68 km/s/Mpc). This period corresponds to a fundamental frequency of approximately 7 cycles over the Hubble time. Transition-z values quoted in the literature generally fall near these minima and may explain the reported wide spread up to the predicted LCDM value of approximately z = 0.77. We also observe second and third harmonics of the fundamental. The scale factor data is analyzed several different ways including smoothing, Fourier transform and autocorrelation. We propose a cosmological scalar field harmonic oscillator model for the observation. On this time scale, for a quantum scalar field, the scalar field mass is extraordinarily small at 3x10-32 eV. Our scalar field density parameter precisely replaces the LCDM dark matter density parameter in the Friedmann equations, resulting in essentially identical data fits, and its present value matches the Planck value. Thus the wave is fundamentally a dark matter wave. We therefore posit that this scalar field manifests itself as the dark matter.
We show that the universal $\alpha$-attractor models of inflation can be realized by including an auxiliary vector field $A_{\mu}$ for the Starobinsky model with the Lagrangian $f(R)=R+R^2/(6M^2)$. If the same procedure is applied to general modified $f(R)$ theories in which the Ricci scalar $R$ is replaced by $R+A_{\mu} A^{\mu}+\beta \nabla_{\mu}A^{\mu}$ with constant $\beta$, we obtain the Brans-Dicke theory with a scalar potential and the Brans-Dicke parameter $\omega_{\rm BD}=\beta^2/4$. We also place observational constraints on inflationary models based on auxiliary vector modified $f(R)$ theories from the latest Planck measurements of the Cosmic Microwave Background (CMB) anisotropies in both temperature and polarization. In the modified Starobinsky model, we find that the parameter $\beta$ is constrained to be $\beta<25$ (68 % confidence level) from the bounds of the scalar spectral index and the tensor-to-scalar ratio.
The observed galaxy power spectrum acquires relativistic corrections from lightcone effects, and these corrections grow on very large scales. Future galaxy surveys in optical, infrared and radio bands will probe increasingly large wavelength modes and reach higher redshifts. In order to exploit the new data on large scales, an accurate analysis requires inclusion of the relativistic effects. This is especially the case for primordial non-Gaussianity and for extending tests of dark energy models to horizon scales. Here we investigate the latter, focusing on models where the dark energy interacts non-gravitationally with dark matter. Interaction in the dark sector can also lead to large-scale deviations in the power spectrum. If the relativistic effects are ignored, the imprint of interacting dark energy will be incorrectly identified and thus lead to a bias in constraints on interacting dark energy on very large scales.
We correlate the positions of radio galaxies in the FIRST survey with the CMB lensing convergence estimated from the Atacama Cosmology Telescope over 470 square degrees to determine the bias of these galaxies. We remove optically cross-matched sources below redshift $z=0.2$ to preferentially select Active Galactic Nuclei (AGN). We measure the angular cross-power spectrum $C_l^{\kappa g}$ at $4.4\sigma$ significance in the multipole range $100<l<3000$, corresponding to physical scales between $\approx$ 2--60 Mpc at an effective redshift $z_{\rm eff}= 1.5$. Modelling the AGN population with a redshift-dependent bias, the cross-spectrum is well fit by the Planck best-fit $\Lambda$CDM cosmological model. Fixing the cosmology we fit for the overall bias model normalization, finding $b(z_{\rm eff}) = 3.5 \pm 0.8$ for the full galaxy sample, and $b(z_{\rm eff})=4.0\pm1.1 (3.0\pm1.1)$ for sources brighter (fainter) than 2.5 mJy. This measurement characterizes the typical halo mass of radio-loud AGN: we find $\log(M_{\rm halo} / M_\odot) = 13.6^{+0.3}_{-0.4}$.
We study the onset of the reheating epoch at the end of axion-driven inflation where the axion is coupled to an Abelian, $U(1)$, gauge field via a Chern-Simons interaction term. We focus primarily on $m^2\phi^2$ inflation and explore the possibility that preheating can occur for a range of coupling values consistent with recent observations and bounds on the overproduction of primordial black holes. We find that for a wide range of parameters preheating is efficient. In certain cases the inflaton is seen to transfer all its energy to the gauge fields within a few oscillations. We find that the gauge fields on sub-horizon scales end in an unpolarized state, due to the existence of strong rescattering between the inflaton and gauge modes. We also present a preliminary study of an axion monodromy model coupled to $U(1)$ gauge fields, seeing a similarly efficient preheating behavior as well as indications that the coupling strength has an effect on the creation of oscillons.
Planck-2015 data seem to favour a large value of the lensing amplitude parameter, $A_{\rm L}=1.22\pm0.10$, in CMB spectra. This result is in $2\sigma$ tension with the lensing reconstruction result, $A_{\rm L}^{\phi\phi}=0.95\pm0.04$. In this paper, we simulate several CMB anisotropy and CMB lensing spectra based on Planck-2015 best-fit cosmological parameter values and Planck blue book beam and noise specifications. We analyse several modified gravity models within the effective field theory framework against these simulations and find that models whose effective Newton constant is enhanced can modulate the CMB anisotropy spectra in a way similar to that of the $A_{\rm L}$ parameter. However, in order to lens the CMB anisotropies sufficiently, like in the Planck-2015 results, the growth of matter perturbations is substantially enhanced and gives a high $\sigma_8$ value. This in turn proves to be problematic when combining these data to other probes, like weak lensing from CFHTLenS, that favour a smaller amplitude of matter fluctuations.
We perform a comprehensive study of models of dark matter (DM) in a Universe with a non-thermal cosmological history, i.e with a phase of pressure-less matter domination before the onset of big-bang nucleosynethesis (BBN). Such cosmological histories are generically predicted by UV completions that contain gravitationally coupled scalar fields (moduli). We classify the different production mechanisms for DM in this framework, generalizing previous works by considering a wide range of DM masses/couplings and allowing for DM to be in equilibrium with a "dark" sector. We identify four distinct parametric regimes for the production of relic DM, and derive accurate semi-analytic approximations for the DM relic abundance. Our results are particularly relevant for supersymmetric theories, in which the standard non-thermally produced DM candidates are disfavored by indirect detection constraints. We also comment on experimental signals in this framework, focusing on novel effects involving the power spectrum of DM density perturbations. In particular, we identify a class of models where the spectrum of DM density perturbations is sensitive to the pressure-less matter dominated era before BBN, giving rise to interesting astrophysical signatures to be looked for in the future. A worthwhile future direction would be to study well-motivated theoretical models within this framework and carry out detailed studies of the pattern of expected experimental signals.
We present deep MUSE integral-field unit (IFU) spectroscopic observations of the giant (~150 x 80 kpc) Ly-alpha halo around the z=4.1 radio galaxy TNJ J1338-1942. This 9-hr observation maps the two-dimensional kinematics of the Ly-alpha emission across the halo. We identify two HI absorbers which are seen against the Ly-alpha emission, both of which cover the full 150 x 80 kpc extent of the halo and so have covering fractions ~1. The stronger and more blue-shifted absorber (dv~1200 km/s) has dynamics that mirror that of the underlying halo emission and we suggest that this high column material (n(HI) ~ 10^19.4 /cm^2), which is also seen in CIV absorption, represents an out-flowing shell that has been driven by the AGN (or star formation) within the galaxy. The weaker (n(HI)~10^14 /cm^2) and less blue shifted (dv~500 km/s) absorber most likely represents material in the cavity between the out-flowing shell and the Ly-alpha halo. We estimate that the mass in the shell must be of order 10^10 Msol -- a significant fraction of the ISM from a galaxy at z=4. The large scales of these coherent structures illustrate the potentially powerful influence of AGN feedback on the distribution and energetics of material in their surroundings. Indeed, the discovery of high-velocity (~1000 km/s), group-halo-scale (i.e. >150 kpc) and mass-loaded winds in the vicinity of the central radio source are broadly in agreement with the requirements of models that invoke AGN-driven outflows to regulate star formation and black-hole growth in massive galaxies at early times.
This paper describes outstanding issues in astrophysics and cosmology that can be solved by astronomical observations in a broad spectral range from far infrared to millimeter wavelengths. The discussed problems related to the formation of stars and planets, galaxies and the interstellar medium, studies of black holes and the development of the cosmological model can be addressed by the planned space observatory Millimetron (the "Spectr-M" project) equipped with a cooled 10-m mirror. Millimetron can operate both as a single-dish telescope and as a part of a space-ground interferometer with very long baseline.
There has been much debate about the origin of the diffuse gamma--ray background in the energy range from a few hundred keV to 10 MeV. At lower energies, AGNs and Seyfert galaxies can explain the background, but their contribution cuts off above $\simeq$ 0.3 MeV. In the MeV range, the spectrum drops sharply for increasing energies. It flattens beyond $\sim$ 10 MeV, and blazars appear to account for the fluxes observed there. That leaves an unexplained window for which different candidate sources have been proposed, including annihilations of weakly interacting massive particles (WIMPS). One candidate are Type Ia supernovae (SNe Ia). Previous estimates of SNe Ia to the cosmic gamma--ray background were based on a restricted number of SN Ia explosion models and, on very limited measurements of the SN Ia rates as a function of redshift $z$. In the present work, we use a wide variety of explosion models and the most recent measurements of the SN Ia rates, which now cover a wide redshift interval. If we adopt the central values of the current measurements of the SN Ia rates, the SN still fall short of the observed background, which, however, has the same spectral shape as the predicted one. Only for a fit running between the 1\,$\sigma$ and 2\,$\sigma$ upper limits of the highest observed SN Ia rates, with their large error bars, it is possible an agreement between model prediction and observations of the gamma--ray background. Uncertainties are discussed and the conclusion is that SNe Ia do, in any case, make a non--negligible contribution to the cosmic gamma--ray background in the MeV range.
We provide an $F(R)$ gravity description of a $\Lambda$CDM bouncing model, without the need for matter fluids or for cosmological constant. As we explicitly demonstrate, the two cosmological eras that constitute the $\Lambda$CDM bouncing model, can be generated by $F(R)$ gravity which can lead to accelerating cosmologies. The resulting F(R) gravity has Einstein frame inflationary properties that have concordance to the latest Planck observational data. Both the $F(R)$ gravity stability properties are thoroughly investigated and also, the gravitational particle production, a feature necessary for the viability of the $\Lambda$CDM bounce scenario, is also addressed. As we will show, the $\Lambda$CDM bounce model can be successfully described by pure $F(R)$ gravity, with appealing phenomenological attributes, which we extensively discuss.
We take a first look at preheating after axion monodromy inflation, assuming a standard coupling between the inflaton field and a scalar matter field. We find that in spite of the fact that the oscillation of the inflaton about the field value which minimizes the potential is anharmonic, there is nevertheless a parametric resonance instability, and we determine the Floquet exponent which describes this instability as a function of the parameters of the axion monodromy potential.
We investigate the chaotic inflationary model using the two-loop effective potential of a self-interacting scalar field theory in curved spacetime. We use the potential which contains a non-minimal scalar curvature coupling and a quartic scalar self-interaction. We analyze the Lyapunov stability of de Sitter solution and show the stability bound. Calculating the inflationary parameters, we systematically explore the spectral index $n_s$ and the tensor-to-scalar ratio $r$, with varying the four parameters, the scalar-curvature coupling $\xi_0$, the scalar quartic coupling $\lambda_0$, the renormalization scale $\mu$ and the e-folding number $N$. It is found that the two-loop correction on $n_s$ is much larger than the leading-log correction, which has previously been studied. We show that the model is consistent with the observation by Planck with WMAP and a recent joint analysis of BICEP2.
A quasar catalogue is presented with a total of 510764 objects including 424748 type 1 QSOs and 26623 type 1 AGN complete from the literature to 25 January 2015. Also included are 25015 high-confidence SDSS-based photometric quasars with radio/X-ray associations, 1595 BL Lac objects, and 32783 type 2 objects. Each object is displayed with arcsecond-accurate astrometry, red and blue photometry, redshift, citations, and radio and X-ray associations where present. Also, 114 new spectroscopically confirmed quasars are presented.
We present a new suite of hydrodynamical simulations and use it to study, in detail, black hole and galaxy properties. The high time, spatial and mass resolution, and realistic orbits and mass ratios, down to 1:6 and 1:10, enable us to meaningfully compare star formation rate (SFR) and BH accretion rate (BHAR) timescales, temporal behaviour and relative magnitude. We find that (i) BHAR and galaxy-wide SFR are typically temporally uncorrelated, and have different variability timescales, except during the merger proper, lasting ~0.2-0.3 Gyr. BHAR and nuclear (<100 pc) SFR are better correlated, and their variability are similar. Averaging over time, the merger phase leads typically to an increase by a factor of a few in the BHAR/SFR ratio. (ii) BHAR and nuclear SFR are intrinsically proportional, but the correlation lessens if the long-term SFR is measured. (iii) Galaxies in the remnant phase are the ones most likely to be selected as systems dominated by an active galactic nucleus (AGN), because of the long time spent in this phase. (iv) The timescale over which a given diagnostic probes the SFR has a profound impact on the recovered correlations with BHAR, and on the interpretation of observational data.
Links to: arXiv, form interface, find, astro-ph, recent, 1502, contact, help (Access key information)
Current constraints on dark matter density profiles from weak lensing are typically limited to radial scales greater than 50-100 kpc. In this paper, we explore the possibility of probing the very inner regions of galaxy/halo density profiles by measuring stacked weak lensing on scales of only a few tens of kpc. Our forecasts focus on scales smaller than the equality radius (Req) where the stellar component and the dark matter component contribute equally to the lensing signal. We compute the evolution of Req as a function of lens stellar mass and redshift and show that Req=7-34 kpc for galaxies with the stellar mass of 10^{9.5}-10^{11.5} solar masses. Unbiased shear measurements will be challenging on these scales. We introduce a simple metric to quantify how many source galaxies overlap with their neighbours and for which shear measurements will be challenging. Rejecting source galaxies with close-by companions results in about a 20 per cent decrease in the overall source density. Despite this decrease, we show that Euclid and WFIRST will be able to constrain galaxy/halo density profiles at Req with signal-to-noise ratio >20 for the stellar mass of >10^{10} solar masses. Weak lensing measurements at Req, in combination with stellar kinematics on smaller scales, will be a powerful means by which to constrain both the inner slope of the dark matter density profile as well as the mass and redshift dependence of the stellar initial mass function.
We constrain the deviation of adiabatic evolution of the Universe using the data on the Cosmic Microwave Background (CMB) temperature anisotropies measured by the {\it Planck} satellite and a sample of 481 X-ray selected clusters with spectroscopically measured redshifts. To avoid antenna beam effects, we bring all the maps to the same resolution. We use a CMB template to subtract the cosmological signal while preserving the Thermal Sunyaev-Zeldovich (TSZ) anisotropies; next, we remove galactic foreground emissions around each cluster and we mask out all known point sources. If the CMB black-body temperature scales with redshift as $T(z)=T_0(1+z)^{1-\alpha}$, we constrain deviations of adiabatic evolution to be $\alpha=-0.007\pm 0.013$, consistent with the temperature-redshift relation of the standard cosmological model. This result could suffer from a potential bias associated with the CMB template, that we quantify it to be less than $-0.02$, but is free from those biases associated with using TSZ selected clusters; it represents the best constraint to date of the temperature-redshift relation of the Big-Bang model, confirming previous results.
The dense shell method for the determination of distances to type IIn supernovae has been briefly reviewed. Applying our method to SN 2006gy, SN 2009ip, and SN 2010jl supernovae, we have obtained distances in excellent agreement with the previously known distances to the parent galaxies. The dense shell method is based on the radiation hydrodynamic model of a supernova. The method of the blackbody model, as well as the correctness of its application for simple estimates of distances from observation data, has been justified.
Ghost-free theories beyond the Horndeski class exhibit a partial breaking of the Vainshtein mechanism inside non-relativistic sources of finite extent. We exploit this breaking to identify new and novel astrophysical probes of these theories. Non-relativistic objects feel a gravitational force that is weaker than that predicted by general relativity. The new equation of hydrostatic equilibrium equation is derived and solved to predict the modified behaviour of stars. It is found that main-sequence stars are dimmer and cooler than their general relativity counterparts but the red giant phase is largely indistinguishable. The rotation curves and lensing potential of Milky Way-like galaxies are calculated. The circular velocities are smaller than predicted by general relativity at fixed radius and the lensing mass is larger than the dynamical mass. We discuss potential astrophysical probes of these theories and identify strong lensing as a particularly promising candidate.
We perform a Bayesian model comparison for scenarios within the quadratic curvaton model, determining the degree to which both are disfavoured with respect to the $\Lambda$CDM concordance model and single-field quadratic inflation, using the recent Planck data release. Despite having three additional model parameters, the simplest curvaton scenario is not disfavoured relative to single-field quadratic inflation, and it becomes favoured against this single-field model when we include the joint BICEP/Keck/Planck analysis. In all cases we assume an instantaneous inflaton decay and no surviving isocurvature perturbations. We show that the current constraints on local non-Gaussianity are insufficiently precise to have any significant impact on the evidence ratios so far. We also determine the precision $\sigma(f_{\mathrm{NL}})$ required by future measurements assuming a fiducial value of $f_{\mathrm{NL}}=-5/4$ or $10.8$ to no longer disfavour the curvaton against the $\Lambda$CDM parametrisation, and we discuss the effect that the predicted increase in precision from future measurements on $f_{\mathrm{NL}}$ may have. We show that our results are not very sensitive to our choice of priors.
We study the low energy effective theory describing gravity with broken spatial diffeomorphism invariance. In the unitary gauge, the Goldstone bosons associated with broken diffeomorphisms are eaten and the graviton becomes a massive spin-2 particle with 5 well-behaved degrees of freedom. In this gauge, the most general theory is built with the lowest dimension operators invariant under only temporal diffeomorphisms. Imposing the additional shift and SO(3) internal symmetries, we analyze the perturbations on a FRW background. At linear perturbation level, the observables of this theory are characterized by six parameters, including the usual cosmological parameters and two additional coupling constants for the symmetry-breaking scalars. We discuss several examples relevant to theories of massive gravity.
We present results from the largest contiguous narrow-band survey in the near-infrared. We have used WIRCam/CFHT and the lowOH2 filter (1.187$\pm$0.005 \mu m) to survey ~10 deg$^2$ of contiguous extragalactic sky in the SA22 field. A total of ~6000 candidate emission-line galaxies are found. We use deep CFHTLS $ugriz$ and UKIDSS DXS $J$ and $K$ data to obtain robust photometric redshifts. We combine our data with HiZELS (COSMOS+UDS) and explore VVDS, VIPERS, KMOS and obtain our own spectroscopic follow-up with FMOS and MOSFIRE to derive large samples of high-redshift emission-line selected galaxies: 3471 H\alpha\ emitters at z=0.8, 1343 [OIII]+H\beta\ emitters at z=1.4 and 572 [OII] emitters at z=2.2. We probe co-moving volumes of >10$^6$ Mpc$^3$ and find significant over-densities, including an 8.5\sigma\ (spectroscopically confirmed) over-density of H\alpha\ emitters at z=0.81. We derive H\alpha, [OIII]+H\beta\ and [OII] luminosity functions at z=0.8, 1.4 and 2.2, respectively, and present implications for future surveys such as EUCLID. Our uniquely large volumes/areas allow us to sub-divide the samples in thousands of randomised combinations of areas and provide a robust empirical measurement of sample/cosmic variance. We show that surveys for star-forming/emission-line galaxies at a depth similar to ours can only overcome cosmic-variance (errors <10%) if they are based on volumes >5x10$^5$ Mpc$^3$; errors on L$^*$ and \Phi$^*$ due to sample (cosmic) variance on surveys probing ~10${^4}$ Mpc$^3$ and ~10$^5$ Mpc$^3$ are typically very high: ~300% and ~40-60%, respectively.
We present a general formulation of the theory for a non-minimally coupled perfect fluid in which both conformal and disformal couplings are present. We discuss how such non-minimal coupling is compatible with the assumptions of a perfect fluid and derive both the Einstein and the fluid equations for such model. We found that, while the Euler equation is significantly modified with the introduction of an extra force related to the local gradients of the curvature, the continuity equation is unaltered, thus allowing for the definition of conserved quantities along the fluid flow. As an application to cosmology and astrophysics we compute the effects of the non-minimal coupling on a Friedmann--Lema\^itre--Robertson--Walker background metric and on the Newtonian limit of our theory.
We demonstrate how the properties of a galaxy depend on the mass of its host dark matter subhalo, using two independent models of galaxy formation. For the cases of stellar mass and black hole mass, the median property value displays a monotonic dependence on subhalo mass. The slope of the relation changes for subhalo masses for which heating by active galactic nuclei becomes important. The median property values are predicted to be remarkably similar for central and satellite galaxies. The two models predict considerable scatter around the median property value, though the size of the scatter is model dependent. There is only modest evolution with redshift in the median galaxy property at a fixed subhalo mass. Properties such as cold gas mass and star formation rate, however, are predicted to have a complex dependence on subhalo mass. In these cases subhalo mass is not a good indicator of the value of the galaxy property. We illustrate how the predictions in the galaxy property - subhalo mass plane differ from the assumptions made in empirical models of galaxy clustering by reconstructing the model output using a subhalo abundance matching scheme. In its simplest form, abundance matching generally does not reproduce the clustering predicted by the models, typically resulting in an overprediction of the clustering signal. We show how the basic abundance matching scheme can be extended to reproduce the model predictions more faithfully, which has implications for the analysis of galaxy clustering, particularly for low abundance samples.
We have undertaken the largest systematic study of the high-mass stellar initial mass function (IMF) to date using the optical color-magnitude diagrams (CMDs) of 85 resolved, young (4 Myr < t < 25 Myr), intermediate mass star clusters (10^3-10^4 Msun), observed as part of the Panchromatic Hubble Andromeda Treasury (PHAT) program. We fit each cluster's CMD to measure its mass function (MF) slope for stars >2 Msun. For the ensemble of clusters, the distribution of stellar MF slopes is best described by $\Gamma=+1.45^{+0.03}_{-0.06}$ with a very small intrinsic scatter. The data also imply no significant dependencies of the MF slope on cluster age, mass, and size, providing direct observational evidence that the measured MF represents the IMF. This analysis implies that the high-mass IMF slope in M31 clusters is universal with a slope ($\Gamma=+1.45^{+0.03}_{-0.06}$) that is steeper than the canonical Kroupa (+1.30) and Salpeter (+1.35) values. Using our inference model on select Milky Way (MW) and LMC high-mass IMF studies from the literature, we find $\Gamma_{\rm MW} \sim+1.15\pm0.1$ and $\Gamma_{\rm LMC} \sim+1.3\pm0.1$, both with intrinsic scatter of ~0.3-0.4 dex. Thus, while the high-mass IMF in the Local Group may be universal, systematics in literature IMF studies preclude any definitive conclusions; homogenous investigations of the high-mass IMF in the local universe are needed to overcome this limitation. Consequently, the present study represents the most robust measurement of the high-mass IMF slope to date. We have grafted the M31 high-mass IMF slope onto widely used sub-solar mass Kroupa and Chabrier IMFs and show that commonly used UV- and Halpha-based star formation rates should be increased by a factor of ~1.3-1.5 and the number of stars with masses >8 Msun are ~25% fewer than expected for a Salpeter/Kroupa IMF. [abridged]
We report new deep ALMA observations aimed at investigating the [CII]158um line and continuum emission in three spectroscopically confirmed Lyman Break Galaxies at 6.8<z<7.1, i.e. well within the re-ionization epoch. With Star Formation Rates of SFR ~ 10 Msun/yr these systems are more representative of the high-z galaxy population than the extreme ones targeted in the past by millimeter observations. At the location of the optical emission, tracing both the Lyalpha line and the far-UV continuum, no [CII] emission is detected. However, for the galaxy with the deepest observation we detect [CII] emission at redshift z=7.107, fully consistent with the Lyalpha redshift, but spatially offset by 0.7" (4 kpc) from the optical emission. These results support expectations from recent models that molecular clouds in the central parts of primordial galaxies are rapidly disrupted by stellar feedback. As a result, [CII] emission mostly arises from more external accreting/satellite clumps of neutral gas. Thermal far-infrared continuum is not detected in any of the three galaxies. However, the upper limits on the infrared-to-UV emission ratio do not exceed those derived in metal- and dust-poor galaxies.
We present the results of a three years monitoring campaigns of the $z = 0.024$ type-1 active galactic nucleus (AGN) PGC50427. Through the use of Photometric Reverberation Mapping with broad and narrow band filters, we determine the size of the broad-line emitting region by measuring the time delay between the variability of the continuum and the H$\alpha$ emission line. The H$\alpha$ emission line responds to blue continuum variations with an average rest frame lag of $19.0 \pm 1.23$ days. Using single epoch spectroscopy we determined a broad-line H$\alpha$ velocity width of 1020 km s$^{-1}$ and in combination with the rest frame lag and adoption a geometric scaling factor $f = 5.5$, we calculate a black hole mass of $M_{BH} \sim 17 \times 10^{6} M_{\odot}$. Using the flux variation gradient method, we separate the host galaxy contribution from that of the AGN to calculate the rest frame 5100\AA~ luminosity at the time of our monitoring campaign. The rest frame lag and the host-subtracted luminosity permit us to derive the position of PGC50427 in the BLR size -- AGN luminosity diagram, which is remarkably close to the theoretically expected relation of $R \propto L^{0.5}$. The simultaneous optical and NIR ($J$ and $K_{s}$) observations allow us to determine the size of the dust torus through the use of dust reverberation mapping method. We find that the hot dust emission ($\sim 1800K$) lags the optical variations with an average rest frame lag of $46.2 \pm 2.60$ days. The dust reverberation radius and the nuclear NIR luminosity permit us to derive the position of PGC50427 on the known $\tau - M{V}$ diagram. The simultaneus observations for the broad-line region and dust thermal emission demonstrate that the innermost dust torus is located outside the BLR in PGC50427, supporting the unified scheme for AGNs. (Abstract shortened, see the manuscript.)
Scalar cosmological perturbations are investigated in the framework of a model with interacting dark energy and dark matter. In addition to these constituents, the inhomogeneous Universe is supposed to be filled with the standard noninteracting constituents corresponding to the conventional $\Lambda$CDM model. The interaction term is chosen in the form of a linear combination of dark sector energy densities with evolving coefficients. The methods of discrete cosmology are applied, and strong theoretical constraints on the parameters of the model are derived. A brief comparison with observational data is performed.
Links to: arXiv, form interface, find, astro-ph, recent, 1502, contact, help (Access key information)
We investigate the impact that warm dark matter (WDM) has in terms of 21cm intensity mapping in the post-reionization Universe at z = 3 - 5. We perform hydrodynamic simulations for 5 different models: cold dark matter and WDM with 1,2,3,4 keV (thermal relic) mass and assign the neutral hydrogen a-posteriori using two different methods that both reproduce observations in terms of column density distribution function of neutral hydrogen systems. Contrary to naive expectations, the suppression of power present in the linear and non-linear matter power spectra, results in an increase of power in terms of neutral hydrogen and 21cm power spectra. This is due to the fact that there is a lack of small mass halos in WDM models with respect to cold dark matter: in order to distribute a total amount of neutral hydrogen within the two cosmological models, a larger quantity has to be placed in the most massive halos, that are more biased compared to the cold dark matter cosmology. We quantify this effect and address significance for the telescope SKA1-LOW, including a realistic noise modeling. The results indicate that we will be able to rule out a 4 keV WDM model with 5000 hours of observations at z > 3, with a statistical significance of > 3 sigma, while a smaller mass of 3 keV, comparable to present day constraints, can be ruled out at more than 2 sigma confidence level with 1000 hours of observations at z > 5.
We will expose in this paper our advances towards a proof of the equivalence between FRW background expansion, during some period of time that contains primordial inflation, and the statistical isotropy of the primordial curvature perturbation $\zeta$ at the end of this period of time. Our motivation rests on the growing interest in the existence of a preferred direction in the Universe hinted by the continuous presence of anomalies in the CMB data.
We present the most general parametrisation of models of dark energy in the form of a scalar field which is explicitly coupled to dark matter. We follow and extend the Parameterized Post-Friedmannian approach, previously applied to modified gravity theories, in order to include interacting dark energy. We demonstrate its use through a number of worked examples and show how the initially large parameter space of free functions can be significantly reduced and constrained to include only a few non-zero coefficients. This paves the way for a model-independent approach to classify and test interacting dark energy theories.
We hypothesize that cosmic string loops are the seeds about which globular clusters accrete. Fixing the cosmic string tension by demanding that the peak in the distribution of masses of objects accreting onto string loops agrees with the peak in the observed mass distribution of globular clusters in our Milky Way galaxy, we then compute the expected number density and mass function of globular clusters, and compare with observations. Our hypothesis naturally explains why globular clusters are the oldest and most dense objects in a galaxy, and why they are found in the halo of the galaxy.
We introduce an alternative parametrisation for the scale dependence of the non-linearity parameter $f_{\rm NL}$ in quasi-local models of non-Gaussianity. Our parametrisation remains valid when $f_{\rm NL}$ changes sign, unlike the commonly adopted power law ansatz $f_{\rm NL}(k) \propto k^{ n_{f_{\rm NL}} }$. We motivate our alternative parametrisation by appealing to the self-interacting curvaton scenario, and as an application, we apply it to the CMB power asymmetry. Explaining the power asymmetry requires a strongly scale dependent non-Gaussianity. We show that regimes of model parameter space where $f_{\rm NL}$ is strongly scale dependent are typically associated with a large $g_{\rm NL}$ and quadrupolar power asymmetry, which can be ruled out by existing observational constraints.
Consistency relations, which relate the squeezed limit of an (N+1)-point correlation function to an N-point function, are non-perturbative symmetry statements that hold even if the associated high momentum modes are deep in the nonlinear regime and astrophysically complex. Recently, Kehagias & Riotto and Peloso & Pietroni discovered a consistency relation applicable to large scale structure. We show that this can be recast into a simple physical statement in Lagrangian space: that the squeezed correlation function (suitably normalized) vanishes. This holds regardless of whether the correlation observables are at the same time or not, and regardless of whether multiple-streaming is present. The simplicity of this statement suggests that an analytic understanding of large scale structure in the nonlinear regime may be particularly promising in Lagrangian space.
We calculate the early universe evolution of perturbations in the dark matter energy density in the context of simple dark sector models containing a GeV scale light mediator. We consider the case that the mediator is long lived, with lifetime up to a second, and before decaying it temporarily dominates the energy density of the universe. We show that for primordial perturbations that enter the horizon around this period, the interplay between linear growth during matter domination and collisional damping can generically lead to a sharp peak in the spectrum of dark matter density perturbation. As a result, the population of the smallest DM halos gets enhanced. Possible implications of this scenario are discussed.
We provide a detailed quantitative description of singular inflation. Its close analogy with finite-time future singularity which is associated to dark energy era is described. Calling and classifying the singularities of such inflation as finite-time cosmological singularities we investigate their occurrence, with special emphasis on the Type IV singularity. The study is performed in the context of a general non-canonical scalar-tensor theory. In addition, the impact of finite time singularities on the slow-roll parameters is also investigated. Particularly, we study three cases, in which the singularity occurs during the inflationary era, at the end, and also we study the case that the singularity occurs much more later than inflation ends. Using the obtained slow-roll parameters, for each case, we calculate explicitly the spectral index of primordial curvature perturbations $n_s$, the associated running of the spectral index $a_s$ and of the scalar-to-tensor ratio $r$ and compare the resulting values to the Planck and BICEP2 data. As we demonstrate, in some cases corresponding to the Type IV singularity, there might be the possibility of agreement with the observational data, when the singularity occurs at the end, or after inflation. However, absolute concordance of all observational indices is not achieved. On the contrary, if the singularity occurs during the inflationary era, this is catastrophic for the theory, since the observational indices become divergent. We also show how a Type IV singularity may be consistently accommodated in the Universe's late time evolution. Finally, we investigate which $F(R)$ gravity can generate the Type IV singularity, with special emphasis on the behavior near the finite time singularity.
It has been widely remarked that compact, massive, elliptical-like galaxies are abundant at high redshifts but exceedingly rare in the Universe today, implying significant evolution such that their sizes at z ~ 2+/-0.6 have increased by factors of 3 to 6 to become today's massive elliptical galaxies. These claims have been based on studies which measured the half-light radii of galaxies as though they are all single component systems. Here we identify 21 spheroidal stellar systems within 90 Mpc that have half-light, major-axis radii R_e < ~2 kpc, stellar masses 0.7x10^{11} < M_*/M_Sun < 1.4x10^{11}, and Sersic indices typically around a value of n=2 to 3. This abundance of compact, massive spheroids in our own backyard - with a number density of 6.9x10^{-6} / Mpc^3 (or 3.5x10^{-5} / Mpc^3 per unit dex in stellar mass) - and with the same physical properties as the high-redshift galaxies, had been over-looked because they are encased in stellar disks which usually result in `galaxy' sizes notably larger than 2 kpc. Moreover, this number density is a lower limit because it has not come from a volume-limited sample. The actual density may be closer to 10^{-4} / Mpc^3, although further work is required to confirm this. We therefore conclude that not all massive `spheroids' have undergone dramatic structural and size evolution since z ~ 2+/-0.6. Given that the bulges of local early-type disk galaxies are known to consist of predominantly old stars which existed at z ~ 2, it seems likely that some of the observed high redshift spheroids did not increase in size by building (3D) triaxial envelopes as commonly advocated, and that the growth of (2D) disks has also been important over the past 9-11 billion years.
We present an analysis of the data produced by the MaNGA prototype run (P-MaNGA), aiming to test how the radial gradients in recent star formation histories, as indicated by the 4000AA-break (D4000), Hdelta absorption (EW(Hd_A)) and Halpha emission (EW(Ha)) indices, can be useful for understanding disk growth and star formation cessation in local galaxies. We classify 12 galaxies observed on two P-MaNGA plates as either centrally quiescent (CQ) or centrally star-forming (CSF), according to whether D4000 measured in the central spaxel of each datacube exceeds 1.6. For each galaxy we generate both 2D maps and radial profiles of D4000, EW(Hd_A) and EW(Ha). We find that CSF galaxies generally show very weak or no radial variation in these diagnostics. In contrast, CQ galaxies present significant radial gradients, in the sense that D4000 decreases, while both EW(Hd_A) and EW(Ha) increase from the galactic center outward. The outer regions of the galaxies show greater scatter on diagrams relating the three parameters than their central parts. In particular, the clear separation between centrally-measured quiescent and star-forming galaxies in these diagnostic planes is largely filled in by the outer parts of galaxies whose global colors place them in the green valley, supporting the idea that the green valley represents a transition between blue-cloud and red-sequence phases, at least in our small sample. These results are consistent with a picture in which the cessation of star formation propagates from the center of a galaxy outwards as it moves to the red sequence.
A variety of supersymmetric models give rise to a split mass spectrum characterized by very heavy scalars but sub-TeV gauginos, usually with a wino-like LSP. Such models predict a thermally-produced underabundance of wino-like WIMP dark matter so that non-thermal DM production mechanisms are necessary. We examine the case where theories with a wino-like LSP are augmented by a Peccei-Quinn sector including an axion-axino-saxion supermultiplet in either the SUSY KSVZ or SUSY DFSZ models and with/without saxion decays to axions/axinos. We show allowed ranges of PQ breaking scale f_a for various cases which are generated by solving the necessary coupled Boltzmann equations. We also present results for a model with radiatively-driven naturalness but with a wino-like LSP.
"High-velocity features" (HVFs) are spectral features in Type Ia supernovae (SNe Ia) that have minima indicating significantly higher (by greater than about 6000 km/s) velocities than typical "photospheric-velocity features" (PVFs). The PVFs are absorption features with minima indicating typical photospheric (i.e., bulk ejecta) velocities (usually ~9000-15,000 km/s near B-band maximum brightness). In this work we undertake the most in-depth study of HVFs ever performed. The dataset used herein consists of 445 low-resolution optical and near-infrared (NIR) spectra (at epochs up to 5 d past maximum brightness) of 210 low-redshift SNe Ia that follow the "Phillips relation." A series of Gaussian functions is fit to the data in order to characterise possible HVFs of Ca II H&K, Si II {\lambda}6355, and the Ca II NIR triplet. The temporal evolution of the velocities and strengths of the PVFs and HVFs of these three spectral features is investigated, as are possible correlations with other SN Ia observables. We find that while HVFs of Ca II are ubiquitous (except in underluminous SNe Ia, where they are never found), HVFs of Si II {\lambda}6355 are significantly rarer, and they tend to exist at the earliest epochs and mostly in objects with large photospheric velocities. It is also shown that stronger HVFs of Si II {\lambda}6355 are found in objects that lack C II absorption at early times and that have red ultraviolet (UV)/optical colours near maximum brightness. These results lead to a self-consistent connection between the presence and strength of HVFs of Si II {\lambda}6355 and many other mutually correlated SN Ia observables, including photospheric velocity.
Links to: arXiv, form interface, find, astro-ph, recent, 1502, contact, help (Access key information)
This article is the second in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our aim is to quantify the scatter, systematic bias and completeness of cluster masses derived from a diverse set of 25 galaxy-based methods using two contrasting mock galaxy catalogues based on a sophisticated halo occupation model and a semi-analytic model. Analysing 968 clusters, we find a wide range in the RMS errors in log M200c delivered by the different methods (0.18 to 1.08 dex, i.e., a factor of ~1.5 to 12), with abundance matching and richness methods providing the best results, irrespective of the input model assumptions. In addition, certain methods produce a significant number of catastrophic cases where the mass is under- or over-estimated by a factor greater than 10. Given the steeply falling high-mass end of the cluster mass function, we recommend that richness or abundance matching-based methods are used in conjunction with these methods as a sanity check for studies selecting high mass clusters. We see a stronger correlation of the recovered to input number of galaxies for both catalogues in comparison with the group/cluster mass, however, this does not guarantee that the correct member galaxies are being selected. We do not observe significantly higher scatter for either mock galaxy catalogues. Our results have implications for cosmological analyses that utilise the masses, richnesses, or abundances of clusters, which have different uncertainties when different methods are used.
We explore the viability of a boson dark matter candidate with an asymmetry between the number densities of particles and antiparticles. A simple thermal field theory analysis confirms that, under certain general conditions, this component would develop a Bose-Einstein condensate in the early universe that, for appropriate model parameters, could survive the ensuing cosmological evolution until now. The condensation of a dark matter component in equilibrium with the thermal plasma is a relativistic process, hence the amount of matter dictated by the charge asymmetry is complemented by a hot relic density frozen out at the time of decoupling. Contrary to the case of ordinary WIMPs, dark matter particles in a condensate can be very light, $10^{-22}\,{\rm eV} \lesssim m \lesssim 10^2\,{\rm eV}$; the lower limit arises from constraints on small-scale structure formation, while the upper bound ensures that the density from thermal relics is not too large. Big-Bang nucleosynthesis constrains the temperature of decoupling to the scale of the QCD phase transition or above. This requires large dark matter-to-photon ratios and very weak interactions with standard model particles. Finally, we argue that a given boson particle that was in thermal equilibrium in the early universe may be in a condensate, or in the form of thermal relics, but we cannot have a combination of both contributing significantly to the mass density today.
We calibrate the halo mass function accounting for halo baryons and present fitting formulae for spherical overdensity masses $M_{500\textrm c}$, $M_{200\textrm c}$, and $M_{200\textrm m}$. We use the hydrodynamical Magneticum simulations, which are well suited because of their high resolution and large cosmological volumes of up to $\sim2$ Gpc$^3$. Baryonic effects globally decrease the masses of galaxy clusters, which, at given mass, results in a decrease of their number density. This effect vanishes at high redshift $z\sim2$ and for high masses $\gtrsim 5\times10^{14}M_\odot$. We perform cosmological analyses of three idealized approximations to the cluster surveys by the South Pole Telescope (SPT), Planck, and eROSITA. For the SPT-like and the Planck-like samples, we find that the impact of baryons on the cosmological results is negligible. In the eROSITA-like case, we find that neglecting the baryonic impact leads to an underestimate of $\Omega_\textrm m$ by about 0.01, which is comparable to the expected uncertainty from eROSITA. We compare our mass function fits with the literature. In particular, in the analysis of our Planck-like sample, results obtained using our mass function are shifted by $\Delta(\sigma_8)\simeq0.05$ with respect to results obtained using the Tinker et al. (2008) fit. This shift represents a large fraction of the observed difference between the latest results from Planck clusters and CMB anisotropies, and the tension is essentially removed. We discuss biases that can be introduced through inadequate mass function parametrizations that introduce false cosmological sensitivity. Additional work to calibrate the halo mass function is therefore crucial for progress in cluster cosmology.
We present an implementation of smoothed particle hydrodynamics (SPH) with improved accuracy for simulations of galaxies and the large-scale structure. In particular, we combine, implement, modify and test a vast majority of SPH improvement techniques in the latest instalment of the GADGET code. We use the Wendland kernel functions, a particle wake-up time-step limiting mechanism and a time-dependent scheme for artificial viscosity, which includes a high-order gradient computation and shear flow limiter. Additionally, we include a novel prescription for time-dependent artificial conduction, which corrects for gravitationally induced pressure gradients and largely improves the SPH performance in capturing the development of gas-dynamical instabilities. We extensively test our new implementation in a wide range of hydrodynamical standard tests including weak and strong shocks as well as shear flows, turbulent spectra, gas mixing, hydrostatic equilibria and self-gravitating gas clouds. We jointly employ all modifications; however, when necessary we study the performance of individual code modules. We approximate hydrodynamical states more accurately and with significantly less noise than standard SPH. Furthermore, the new implementation promotes the mixing of entropy between different fluid phases, also within cosmological simulations. Finally, we study the performance of the hydrodynamical solver in the context of radiative galaxy formation and non-radiative galaxy cluster formation. We find galactic disks to be colder, thinner and more extended and our results on galaxy clusters show entropy cores instead of steadily declining entropy profiles. In summary, we demonstrate that our improved SPH implementation overcomes most of the undesirable limitations of standard SPH, thus becoming the core of an efficient code for large cosmological simulations.
Many recent studies have highlighted certain failures of the standard Eulerian-space cosmological perturbation theory (SPT). Its problems include (1) not capturing large-scale bulk flows [leading to an O(1) error in the 1-loop SPT prediction for the baryon acoustic peak in the correlation function], (2) assuming that the Universe behaves as a pressureless, inviscid fluid, and (3) treating fluctuations on scales that are non-perturbative as if they were. Recent studies have highlighted the successes of perturbation theory in Lagrangian space or theories that solve equations for the effective dynamics of smoothed fields. Both approaches mitigate some or all of the aforementioned issues with SPT. We discuss these physical developments by specializing to the simplified 1D case of gravitationally interacting sheets, which allows us to substantially reduces the analytic overhead and still (as we show) maintain many of the same behaviors as in 3D. In 1D, linear-order Lagrangian perturbation theory ("the Zeldovich approximation") is exact up to shell crossing, and we prove that n^{th}-order Eulerian perturbation theory converges to the Zeldovich approximation as n goes to infinity. In no 1D cosmology that we consider (including a CDM-like case and power-law models) do these theories describe accurately the matter power spectrum on any mildly nonlinear scale. We find that theories based on effective equations are much more successful at describing the dynamics, and we test some of the assumptions that underlie prior applications of these theories. Finally, we discuss many topics that have recently appeared in the perturbation theory literature such as beat coupling, the shift and smearing of the baryon acoustic oscillation feature, and the advantages of Fourier versus configuration space. Our simplified 1D case serves as an intuitive review of these perturbation theory results.
We show that the observed rotation curves of spiral galaxies constrain the sound speed of the dark matter to be $c_s < 10^{-4} c$, where $c$ is the speed of light in vacuum. Using the Modified Chaplygin Gas as a representative example of a class of unified dark energy models incorporating an effective dark matter component with a non-zero sound speed, we determine the most stringent constraint to date on the value of the constant contribution to the equation of state parameter in this class of models. Finally, we explain the reason why previous constraints using the Cosmic Microwave Background and Baryonic Acoustic Oscillations were not as competitive as the one presented in this paper and discuss the limitations of the recently proposed Extended Chaplygin Gas.
Cosmic voids have been shown to be an effective probe of cosmology, complementary to galaxy clusters. In this work, we present a simple theoretical framework for predicting of the size of the largest voids expected within a given redshift and volume. Our model is based on the exact extreme-value statistics which has previously been successfully applied to massive galaxy clusters. We implement our formalism using the void-abundance model of Sheth and Van de Weygaert and find results that are consistent with dark-matter simulation. We discuss several ways in which our framework can be adapted to other void models.
We describe a simple class of cosmological models called $\alpha$ attractors, which provide an excellent fit to the latest Planck data. These theories are most naturally formulated in the context of supergravity with logarithmic Kahler potentials. We develop generalized versions of these models which can describe not only inflation but also dark energy and supersymmetry breaking.
We discuss metric perturbations of the relativistic diffusion equation around the homogeneous Juttner equilibrium of massless particles in a homogeneous expanding universe. The metric perturbation describes matter distribution and the gravitational wave background in an inhomogeneous universe. We show that the lowest order perturbation can be treated as a variation of temperature. We derive a formula expressing temperature fluctuations in terms of the diffusion and tensor power spectrum. We discuss the multipole expansion of the fluctuations in the presence of diffusion.
An excess in gamma-rays emanating from the galactic centre has recently been observed in the Fermi-LAT data. This signal can be interpreted as resulting from WIMP annihilation, with the spectrum well-fit by dark matter annihilating dominantly into either bottom-quark or Higgs pairs. Supersymmetric models provide a well-motivated framework to study the implications of this signal in these channels. With a neutralino dark matter candidate, the gamma-ray excess cannot be easily accommodated in the minimal supersymmetric model, which in any case requires tuning below the percent level to explain the observed Higgs mass. Instead we are naturally led to consider the next-to-minimal model with a singlet superfield. This not only allows for the annihilation channel into bottom-quark pairs to be implemented, but also provides new possibilities for annihilation into Higgs-pseudoscalar pairs. We show that the fit to the gamma-ray excess for the Higgs-pseudoscalar channel can be just as good as for annihilation into bottom-quark pairs. Moreover, in the parameter range of interest, the next-to-minimal supersymmetric model solves the mu-problem and can explain the 125 GeV Higgs mass with improved naturalness. We also consider an extension by adding a right-handed neutrino superfield with the right-handed sneutrino acting as a dark matter candidate. Interestingly, this allows for the annihilation into pseudoscalar pairs which also provide a good fit to the gamma-ray excess. Furthermore, in the case of a neutralino LSP, the late decay of a sneutrino NLSP can non-thermally produce the observed relic abundance. Finally, the WIMP annihilation into scalar pairs allows for the possibility of detecting the Higgs or pseudoscalar decay into two photons, providing a smoking-gun signal of the model.
The analyticity of response functions and scattering amplitudes implies powerful relations between low-energy observables and the underlying short-distance dynamics. These 'IR/UV' relations are rooted in basic physical principles, such as causality and unitarity. In this paper, we seek similar connections in inflation, relating cosmological observations to the physics responsible for the accelerated expansion. We assume that the inflationary theory is Lorentz invariant at short distances, but allow for non-relativistic interactions and a non-trivial speed of propagation at low energies. Focusing on forward scattering, we derive a 'sum rule' which equates a combination of low-energy parameters to an integral which is sensitive to the high-energy behavior of the theory. While for relativistic amplitudes unitarity is sufficient to prove positivity of the sum rule, this is not guaranteed in the non-relativistic case. We discuss the conditions under which positivity still applies, and show that they are satisfied by all known UV completions of single-field inflation. In that case, we obtain a consistency condition for primordial non-Gaussianity, which constrains the size and the sign of the equilateral four-point function in terms of the amplitude of the three-point function. The resulting bound rules out about half of the parameter space that is still allowed by current observations. Finding a violation of our consistency condition would point towards less conventional theories of inflation, or violations of basic physical principles.
Here we introduce GAMESH, a novel pipeline which implements self-consistent radiative and chemical feedback in a computational model of galaxy formation. By combining the cosmological chemical-evolution model GAMETE with the radiative transfer code CRASH, GAMESH can post process realistic outputs of a N-body simulation describing the redshift evolution of the forming galaxy. After introducing the GAMESH implementation and its features, we apply the code to a low-resolution N-body simulation of the Milky Way formation and we investigate the combined effects of self-consistent radiative and chemical feedback. Many physical properties, which can be directly compared with observations in the Galaxy and its surrounding satellites, are predicted by the code along the merger-tree assembly. The resulting redshift evolution of the Local Group star formation rates, reionisation and metal enrichment along with the predicted Metallicity Distribution Function of halo stars are critically compared with observations. We discuss the merits and limitations of the first release of GAMESH, also opening new directions to a full implementation of feedback processes in galaxy formation models by combining semi-analytic and numerical methods.
Using wide field narrow-band surveys, we provide a new measurement of the $z=6.6$ Lyman-$\alpha$ Emitter (LAE) luminosity function (LF), which constraints the bright end for the first time. We use a combination of archival narrow-band NB921 data in UDS and new NB921 measurements in SA22 and COSMOS/UltraVISTA, all observed with the Subaru telescope, with a total area of $\sim 5$ deg$^2$. We exclude lower redshift interlopers by using broad-band optical and near-infrared photometry and also exclude three supernovae with data split over multiple epochs. We spectroscopically confirm the two most luminous Ly$\alpha$ emitters ever found at $z=6.604$ and $6.541$ in the COSMOS field using Keck/DEIMOS and VLT/FORS2. Combining the UDS and COSMOS samples we find no evolution of the bright end of the Ly$\alpha$ LF between $z=5.7$ and $6.6$, which is supported by spectroscopic follow-up, and conclude that \emph{Himiko}-like sources are not as rare as previously thought, with number densities of $\sim 1.5\times10^{-5}$ Mpc$^{-3}$. Combined with our wide-field SA22 measurements, our results indicate a non-Schechter-like bright end of the LF at $z=6.6$ and a different evolution of \emph{observed} faint and bright LAEs. This differential evolution was not addressed in previous studies, or discarded as cosmic variance, but we argue instead that it may be an effect of re-ionisation. Using a toy-model, we show that such differential evolution of the LF is expected, since brighter sources are able to ionise their surroundings earlier, such that Ly$\alpha$ photons are able to escape. Our targets are excellent candidates for detailed follow-up studies and provide the possibility to give a unique view on the earliest stages in the formation of galaxies and re-ionisation process.
Galaxies and the dark matter halos that host them are not spherically
symmetric, yet spherical symmetry is a helpful simplifying approximation for
idealised calculations and analysis of observational data. The assumption leads
to an exact conservation of angular momentum for every particle, making the
dynamics unrealistic. But how much does that inaccuracy matter in practice for
analyses of stellar distribution functions, collisionless relaxation, or dark
matter core-creation?
We provide a general answer to this question for a wide class of aspherical
systems; specifically, we consider distribution functions that are "maximally
stable", i.e. that do not evolve at first order when external potentials (which
arise from baryons, large scale tidal fields or infalling substructure) are
applied. We show that a spherically-symmetric analysis of such systems gives
rise to the false conclusion that the density of particles in phase space is
ergodic (a function of energy alone).
Using this idea we are able to demonstrate that: (a) observational analyses
that falsely assume spherical symmetry are made more accurate by imposing a
strong prior preference for near-isotropic velocity dispersions in the centre
of spheroids; (b) numerical simulations that use an idealised
spherically-symmetric setup can yield misleading results and should be avoided
where possible; and (c) triaxial dark matter halos (formed in collisionless
cosmological simulations) nearly attain our maximally-stable limit, but their
evolution freezes out before reaching it.
We present some "back-of-the-envelope" calculations to try to understand cold dark matter, its searches, and extensions of the Standard Model. Some of the insights obtained from this exercise may be useful.
We discuss an extension of the massively parallel cosmological simulation code GADGET-2, which enables investigation of distinct gravitational force laws between particle species. In addition to simplifying investigations of a universally modified force law, the ngravs extension allows state-of-the-art collisionless cosmological simulations of quite exotic gravitational scenarios. We briefly review the algorithms used by GADGET-2, and present our extension to multiple gravities, highlighting additional features that facilitate consideration of exotic force laws. We discuss the accuracy and performance of the ngravs extension, both internally and with an unaltered GADGET-2, under all relevant operational modes. The ngravs extension is publicly released to the research community.
Detection of 21 cm emission of HI from the epoch of reionization, at redshifts z>6, is limited primarily by foreground emission. We investigate the signatures of wide-field measurements and an all-sky foreground model using the delay spectrum technique that maps the measurements to foreground object locations through signal delays between antenna pairs. We demonstrate interferometric measurements are inherently sensitive to all scales, including the largest angular scales, owing to the nature of wide-field measurements. These wide-field effects are generic to all observations but antenna shapes impact their amplitudes substantially. A dish-shaped antenna yields the most desirable features from a foreground contamination viewpoint, relative to a dipole or a phased array. Comparing data from recent Murchison Widefield Array observations, we demonstrate that the foreground signatures that have the largest impact on the HI signal arise from power received far away from the primary field of view. We identify diffuse emission near the horizon as a significant contributing factor, even on wide antenna spacings that usually represent structures on small scales. For signals entering through the primary field of view, compact emission dominates the foreground contamination. These two mechanisms imprint a characteristic pitchfork signature on the "foreground wedge" in Fourier delay space. Based on these results, we propose that selective down-weighting of data based on antenna spacing and time can mitigate foreground contamination substantially by a factor ~100 with negligible loss of sensitivity.
After inflation, a period of preheating may have produced a stochastic background of high frequency gravitational waves (GWs) that would persist until today. The nature of the inflaton's coupling to Standard Model or other fields is unknown, so it is useful to ask what features such fields may typically have, and how these affect predictions for the GW's produced. Here we consider the inflaton to be coupled to a light scalar field, and show that even a very small quartic self-interaction term will reduce the amplitude of the GW spectrum. For self-coupling $\lambda_{\chi} \gtrsim g^2$, where $g^2$ is the inflaton-scalar coupling, the peak energy density goes as $\Omega_{\rm gw}^{(\lambda_{\chi})} / \Omega_{\rm gw}^{(\lambda_{\chi}=0)} \sim (g^2/\lambda_{\chi})^{2}$. A consequence is that if the universe reheats through an inflaton-Higgs coupling then the spectrum would be suppressed but the dynamics would be sensitive to the Higgs potential near the energy scale of inflation.
In order to estimate the effects of local structure on the Hubble parameter we calculate the low-redshift expansion for $H(z)$ and $\frac{\delta H}{H}$ for an observer at the center of a spherically symmetric matter distribution in presence of a cosmological constant. We then test the accuracy of the formulae comparing them with fully relativistic non pertubative numerical calculations for different cases for the density profile. Our analytical approach does not involve perturbation theory, and is based on the use of exact solutions of Einstein's field equations. The formulae can be used to take into account the non linear effects on the Hubble expansion parameter due to the monopole component of the local structure.
We have examined a solution to the FRW model of the Einstein and de Sitter Universe, often termed the standard model of cosmology, using wide values for the normalized cosmological constant Omega_L and spacetime curvature Omega_k with proposed values of normalized matter density. These solutions were evaluated using a combination of the third type of elliptical equations and were found to display critical points for redshift z, between 1 and 3, when Omega_L is positive. These critical points occur at values for normalized cosmological constant higher than those currently thought important, though we find this solution interesting because the Omega_L term may increase in dominance as the Universe evolves bringing this discontinuity into importance. We also find positive Omega_L tends towards attractive at values of z which are commonly observed for distant galaxies.
The source-subtracted, 1.1 and 1.6 {\mu}m NICMOS images used in earlier analyses of the near-infrared Hubble Ultra Deep Field contained residual flux in extended wings of identified sources that contributed an unknown amount to fluctuation power. When compared to the original results, a reanalysis after subtracting this residual flux shows that mean-square and rms fluctuations decrease a maximum of 52 and 31 per cent at 1.6 {\mu}m and 50 and 30 per cent at 1.1 {\mu}m. However, total mean-square fluctuations above 0.5 arcsec only decrease 6.5 and 1.4 per cent at 1.6 and 1.1 {\mu}m, respectively. These changes would not affect any published conclusions based on the prior analyses. These results exclude previous suggestions that extended wings of detected galaxies may be a major contributor to the source-subtracted near-infrared background and confirm that most fluctuation power in these images must be explained by other means.
Links to: arXiv, form interface, find, astro-ph, recent, 1502, contact, help (Access key information)