We propose that stellar-mass binary black holes like GW150914 will become a tool to explore the local Universe within ~100Mpc in the era of evolved Laser Interferometer Space Antenna (eLISA). High calibration accuracy and annual motion of eLISA could enable us to localize up to ~60 binaries more accurately than the error volume of ~100Mpc^3 without presumably absent electromagnetic counterparts. This accuracy will give us a fair chance to determine the host object solely by gravitational waves. By combining the luminosity distance extracted from gravitational waves with the cosmological redshift determined from the host, the local value of the Hubble parameter will be determined up to a few % without relying on the empirically-constructed distance ladder. Gravitational-wave cosmography would pave the way for resolution of the disputed Hubble tension, where the local and global measurements disagree in the value of the Hubble parameter at 3.4sigma level, which amounts to ~9%.
We consider various models realizing baryogenesis during the electroweak phase transition (EWBG). Our focus is their possible detection in future collider experiments and possible observation of gravitational waves emitted during the phase transition. We also discuss the possibility of a non-standard cosmological history which can facilitate EWBG. We show how acceptable parameter space can be extended due to such a modification and conclude that next generation precision experiments such as the ILC will be able to confirm or falsify many models realizing EWBG. We also show that, in general, collider searches are a more powerful probe than gravitational wave searches. However, observation of a deviation from the SM without any hints of gravitational waves can point to models with modified cosmological history that generically enable EWBG with weaker phase transition and thus, smaller GW signals.
We update and extend our earlier discussion of the potential of a next generation space-borne CMB experiment for studies of extragalactic sources. Our analysis has particular bearing on the definition of a future space project, CORE, that will be submitted in response to ESA's call for a Medium-size mission opportunity (M5) as the successor of the Planck satellite. Even though the effective telescope size will be similar or somewhat smaller than that of Planck, CORE will have a considerably better angular resolution at its highest frequencies, since, at variance with Planck, it will be diffraction limited. The better resolution implies a substantial decrease of the source confusion, i.e. substantially fainter detection limits. In particular, CORE will detect several thousands of strongly lensed high-z galaxies distributed over the full sky. These are the brightest (sub-)mm sources in the sky and, as such, will allow studies of high-z star forming galaxies in exquisite detail. Also, the CORE resolution matches the typical sizes of high-z galaxy proto-clusters much better than the Planck resolution, resulting in a much higher detection efficiency. These objects will be caught in an evolutionary phase beyond the reach of surveys in other wavebands. Furthermore, CORE will provide unique information on the evolution of the star formation in virialized groups and clusters of galaxies up to the highest possible redshifts. Finally, thanks to its very high sensitivity, CORE will detect the polarized emission of thousands of radio sources in a poorly explored frequency range, and, for the first time, of dusty galaxies.
Links to: arXiv, form interface, find, astro-ph, recent, 1609, contact, help (Access key information)
Despite strong theoretical arguments, the use of clusters as cosmological probes is, in practice, frequently questioned because of the many uncertainties impinging on cluster mass estimates. Our aim is to develop a fully self-consistent cosmological approach of X-ray cluster surveys, exclusively based on observable quantities, rather than masses. This procedure is justified given the possibility to directly derive the cluster properties via ab initio modelling, either analytically or by using hydrodynamical simulations. In this third paper, we evaluate the method on cluster toy-catalogues. We model the population of detected clusters in the count-rate -- hardness-ratio -- angular size -- redshift space and compare the corresponding 4-dimensional diagram with theoretical predictions. The best cosmology+physics parameter configuration is determined using a simple minimisation procedure; errors on the parameters are derived by scanning the likelihood hyper-surfaces with a wide range of starting values. The method allows a simultaneous fit of the cosmological parameters, of the cluster evolutionary physics and of the selection effects. When using information from the X-ray survey alone plus redshifts, this approach is shown to be as accurate as the mass function for the cosmological parameters and to perform better for the cluster physics, as modelled in the scaling relations. It enables the identification of degenerate combinations of parameter values. Given the considerably shorter computer times involved for running the minimisation procedure in the observed parameter space, this method appears to clearly outperform traditional mass-based approaches when X-ray survey data alone are available.
In weak gravitational lensing, weighted quadrupole moments of the brightness profile in galaxy images are a common way to estimate gravitational shear. We employ general adaptive moments (GLAM) to study causes of shear bias on a fundamental level and for a practical definition of an image ellipticity. For GLAM, the ellipticity is identical to that of isophotes of elliptical images, and this ellipticity is always an unbiased estimator of reduced shear. Our theoretical framework reiterates that moment-based techniques are similar to a model-based approach in the sense that they fit an elliptical profile to the image to obtain weighted moments. As a result, moment-based estimates of ellipticities are prone to underfitting bias. The estimation is fundamentally limited mainly by pixellation which destroys information on the original, pre-seeing image. We give an optimized estimator for the pre-seeing GLAM ellipticity and its bias for noise-free images. To deal with images where pixel noise is prominent, we consider a likelihood model of GLAM parameters in the pre-seeing frame. Similar to the noise-free case, this likelihood is biased in the presence of underfitting. The bias does not vary with the overall noise level but it depends in detail on the correlation of pixel noise as well as the noise homogeneity over the image, which could be relevant for the calibration strategies of other methodologies. We give an analytic expression for the underfitting bias and suggest means to reduce it. Moreover, within a Bayesian framework of the GLAM ellipticity, noise-dependent bias emerges after marginalisation of the likelihood over image size or centroid position, even in the absence of underfitting. Therefore, a Bayesian approach to shape measurements does not necessarily mitigate noise bias although our tests show that it can be reduced. [abriged]
We contribute another anisotropy study to this field of research using Supernovae Type Ia (SNe Ia). In this work, we utilise the power spectrum calculation method and apply it to both the current SNe Ia data and simulation. Our simulations are constructed with the characteristics of the upcoming survey of the Large Synoptic Survey Telescope (LSST), which shall bring us the largest SNe Ia collection to date. We make predictions for the amplitude of a possible dipole anisotropy or anisotropy in higher multipole moments that would be detectable by the LSST.
Third generation ground-based interferometers as well as the planned space-based interferometer LISA are expected to detect a plethora of gravitational wave signals from coalescing binaries at cosmological distance. The emitted gravitational waves propagate in the expanding universe through the inhomogeneous distribution of matter. Here we show that the acceleration of the universe and the peculiar acceleration of the binary with respect to the observer distort the gravitational chirp signal from the simplest General Relativity prediction, affecting parameter estimations for the binaries visible by LISA. We find that the effect due to peculiar acceleration can be much larger than the one due to the universe acceleration, thereby excluding the possibility of using this latter to infer the redshift of the GW source (as previously proposed). Moreover, peculiar accelerations can introduce a bias in the estimation of parameters such as the time of coalescence and the individual masses of the binary. An error in the estimation of the arrival time will have an impact in the case of sources visible first by LISA and later by ground based interferometers.
Gravitational lensing by large-scale structure significantly impacts observations of the cosmic microwave background (CMB): it smooths the acoustic peaks in temperature and $E$-mode polarization power spectra, correlating previously uncorrelated modes; and it converts $E$-mode polarization into $B$-mode polarization. The act of measuring and removing the effect of lensing from CMB maps, or delensing, has been well studied in the context of $B$ modes, but little attention has been given to the delensing of the temperature and $E$ modes. In this paper, we model the expected delensed $T$ and $E$ power spectra to all orders in the lensing potential, demonstrating the sharpening of the acoustic peaks and a significant reduction in lens-induced power spectrum covariances. We then perform cosmological forecasts, demonstrating that delensing will yield improved sensitivity to parameters with upcoming surveys. We highlight the breaking of the degeneracy between the effective number of neutrino species and primordial helium fraction as a concrete application. We also show that delensing increases cosmological information as long as the measured lensing reconstruction is included in the analysis. We conclude that with future data, delensing will be crucial not only for primordial $B$-mode science but for a range of other observables as well.
Galactic systems, and the Universe at large, exhibit large dynamical anomalies: The observed matter in them falls very short of providing enough gravity to account for their dynamics. The mainstream response to this conundrum is to invoke large quantities of `dark matter' -- which purportedly supplies the needed extra gravity -- and also of `dark energy', to account for further anomalies in cosmology, such as the observed, accelerated expansion. The MOND paradigm offers a different solution: a breakdown of standard dynamics (gravity and/or inertia) in the limit of low accelerations -- below some acceleration $a_0$. In this limit, dynamics become space-time scale invariant, and is controlled by a gravitational constant $\mathcal{A}_0\equiv Ga_0$, which replaces Newton's $G$. With the new dynamics, the various detailed manifestations of the anomalies in galaxies are predicted with no need for dark matter. The cosmological anomalies could, but need not have to do with small accelerations. For example, the need for dark matter in accounting for the expansion history of the Universe is eliminated if the relevant gravitational constant is $\approx 2\pi G$. Such a `renormalization' of $G$ could be a dimensionless parameter of a MOND theory. The constant $a_0$ turns out to carry cosmological connotations, in that $2\pi a_0\approx cH_0\approx c^2(\Lambda/3)^{1/2}$, where $H_0$ is the present expansion rate of the Universe, and $\Lambda$ the measured `cosmological constant'. There are MOND theories in which this `coincidence' is natural. I draw on enlightening historical and conceptual analogies from quantum theory to limelight aspects of MOND. I also explain how MOND may have strong connections with effects of the quantum vacuum on local dynamics.
I derive a new MOND relation for pure-disc galaxies: The `dynamical' central surface density, $\Sigma^0_D$, deduced from the measured velocities, is a universal function of only the true, `baryonic' central surface density, $\Sigma^0_B$: $\Sigma^0_D=\Sigma_M \mathcal{S}(\Sigma^0_B/\Sigma_M)$, where $\Sigma_M\equiv a_0/2\pi G$ is the MOND surface density constant. This surprising result is shown to hold in both existing, nonrelativistic MOND theories. $\mathcal{S}(y)$ is derived: $\mathcal{S}(y)=\int_0^y\nu(y')dy'$, with $\nu(y)$ the interpolating function of the theory. The relation aymptotes to $\Sigma^0_D=\Sigma^0_B$ for $\Sigma^0_B\gg\Sigma_M$, and to $\Sigma^0_D=(4\Sigma_M\Sigma^0_B)^{1/2}$ for $\Sigma^0_B\ll\Sigma_M$. This study was prompted by the recent finding of a correlation between related attributes of disc galaxies by Lelli et al. (2016). The MOND central-surface-densities relation agrees very well with these results.
We briefly review how X-ray observations of high-redshift active galactic nuclei (AGNs) at z = 4-7 have played a critical role in understanding their basic demographics as well as their physical processes; e.g., absorption by nuclear material and winds, accretion rates, and jet emission. We point out some key remaining areas of uncertainty, highlighting where further Chandra and XMM-Newton observations/analyses, combined with new multiwavelength survey data, can advance understanding over the next decade.
In the standard model extended with a seesaw mass matrix, we study the
production of sterile neutrinos from the decay of vector bosons at $T\simeq
M_{W,Z}$. We derive a general quantum kinetic equation for the production of
sterile neutrinos and their effective mixing angles valid in a wide range of
temperature, to all orders in interactions of the standard model, and to
leading order mixing angle $\theta \ll 1$.
Production rates and effective mixing angles depend sensitively on helicity.
Positive helicity states interact more weakly with the medium and their
effective mixing angle is not modified significantly whereas the mixing angle
for negative helicity is strongly suppressed by the medium.
If $M_s \lesssim 8.35\,\mathrm{MeV}$, there are fewer states with negative
helicity produced than those with positive helicity. There is an MSW resonance
in the absence of lepton asymmetry, but is screened by the damping rate,
without production enhancement. Negative helicity states freeze-out at
$T^-_f\simeq 5\,\mathrm{GeV}$ and positive helicity states freeze-out at $T^+_f
\simeq 8\,\mathrm{GeV}$, both distributions are far from thermal. Negative
helicity states feature a broader momentum distribution than that for those
with positive helicity. Sterile neutrinos produced via vector boson decay do
not satisfy the abundance, lifetime and cosmological constraints to be the sole
dark matter component in the universe. We discuss how massive sterile neutrinos
might solve the $^{7}Li$ problem, albeit at the very edge of the possible
parameter space. A heavy sterile neutrino with a mass of a few MeV could decay
into light sterile neutrinos, of a few keV in mass, that contribute to warm
dark matter. We argue that heavy sterile neutrinos with lifetime $\leq 1/H_0$
reach local thermodynamic equilibrium.
Narrow-line Seyfert 1 galaxies have been identified by the Fermi Gamma-Ray Space Telescope as a rare class of gamma-ray emitting active galactic nuclei (AGN). The lowest-redshift candidate among them is the source 1H 0323+342. Here we present quasi-simultaneous Gemini near-infrared and Keck optical spectroscopy for it, from which we derive a black hole mass based on both the broad Balmer and Paschen emission lines. We supplement these observations with a NuSTAR X-ray spectrum taken about two years earlier, from which we constrain the black hole mass based on the short timescale spectral variability. Our multiwavelength observations suggest a black hole mass of ~2x10^7 solar masses, which agrees well with previous estimates. We build the spectral energy distribution and show that it is dominated by the thermal and reprocessed emission from the accretion disc rather than the non-thermal jet component. A detailed spectral fitting with the energy-conserving accretion disc model of Done et al. constrains the Eddington ratio to L/L_Edd ~ 0.5 for a (non-rotating) Schwarzschild black hole and to L/L_Edd ~ 1 for a Kerr black hole with dimensionless spin of a*=0.8. Higher spin values and so higher Eddington ratios are excluded, since they would strongly overpredict the observed soft X-ray flux.
We present the light curves of the hydrogen-poor superluminous supernovae (SLSNe-I) PTF12dam and iPTF13dcc, discovered by the (intermediate) Palomar Transient Factory. Both show excess emission at early times and a slowly declining light curve at late times. The early bump in PTF12dam is very similar in duration (~10 days) and brightness relative to the main peak (2-3 mag fainter) compared to those observed in other SLSNe-I, such as SN2006oz, LSQ14bdq and DES14X3taz. In contrast, the long-duration (>30 days) early excess emission in iPTF13dcc, whose brightness competes with that of the main peak, appears to be of a different nature. We construct bolometric light curves for both targets, and fit a variety of light-curve models to both the early bump and main peak in an attempt to understand the nature of these explosions. Even though the slope of the late-time light-curve decline in both SLSNe is suggestively close to that expected from the radioactive decay of $^{56}$Ni and $^{56}$Co, the amount of nickel required to power the full light curves is too large considering the estimated ejecta mass. The magnetar model including an increasing escape fraction provides a reasonable description of the PTF12dam observations, even though it somewhat overpredicts the luminosity at late times. However, neither the basic nor the double-peaked magnetar model is capable of reproducing the iPTF13dcc light curve. A model combining the shock breakout in an extended envelope, responsible for the early excess emission, with late-time magnetar energy injection causing the main peak, provides a reasonable fit to the iPTF13dcc observations. Finally, we find that the light curves of both PTF12dam and iPTF13dcc can be adequately fit with the circumstellar medium (CSM) interaction model.
Links to: arXiv, form interface, find, astro-ph, recent, 1609, contact, help (Access key information)
We present a measurement of galaxy-galaxy lensing around a magnitude-limited ($i_{AB} < 22.5$) sample of galaxies selected from the Dark Energy Survey Science Verification (DES-SV) data. We split these lenses into three photometric-redshift bins from 0.2 to 0.8, and determine the product of the galaxy bias $b$ and cross-correlation coefficient between the galaxy and dark matter overdensity fields $r$ in each bin, using scales above 4 Mpc/$h$ comoving, where we find the linear bias model to be valid given our current uncertainties. We compare our galaxy bias results from galaxy-galaxy lensing with those obtained from galaxy clustering (Crocce et al. 2016) and CMB lensing (Giannantonio et al. 2016) for the same sample of galaxies, and find our measurements to be in good agreement with those in Crocce et al. (2016), while, in the lowest redshift bin ($z\sim0.3$), they show some tension with the findings in Giannantonio et al. (2016). Our results are found to be rather insensitive to a large range of systematic effects. We measure $b\cdot r$ to be $0.87\pm 0.11$, $1.12 \pm 0.16$ and $1.24\pm 0.23$, respectively for the three redshift bins of width $\Delta z = 0.2$ in the range $0.2<z <0.8$, defined with the photometric-redshift algorithm BPZ. Using a different code to split the lens sample, TPZ, leads to changes in the measured biases at the 10-20% level, but it does not alter the main conclusion of this work: when comparing with Crocce et al. (2016) we do not find strong evidence for a cross-correlation parameter significantly below one in this galaxy sample, except possibly at the lowest redshift bin ($z\sim 0.3$), where we find $r = 0.71 \pm 0.11$ when using TPZ, and $0.83 \pm 0.12$ with BPZ, assuming the difference between the results from the two probes can be solely attributed to the cross-correlation parameter.
In the weak field regime, gravitational waves can be considered as being made up of collisionless, relativistic tensor modes that travel along null geodesics of the perturbed background metric. We work in this geometric optics picture to calculate the anisotropies in gravitational wave backgrounds resulting from astrophysical and cosmological sources. Our formalism yields expressions for the angular power spectrum of the anisotropies. We show how the anisotropies are sourced by intrinsic, Doppler, Sachs-Wolfe, and Integrated Sachs-Wolfe terms in analogy with Cosmic Microwave Background photons.
We use bootstrapping to estimate the bias of concentration estimates on N-body dark matter halos as a function of particle number. We find that algorithms based on the maximum radial velocity and radial particle binning tend to overestimate the concentration by 15%-20% for halos sampled with 200 particles and by 7% - 10% for halos sampled with 500 particles. To control this bias at low particle numbers we propose a new algorithm that estimates halo concentrations based on the integrated mass profile. The method uses the full particle information without any binning, making it reliable in cases when low numerical resolution becomes a limitation for other methods. This method reduces the bias to less than 3% for halos sampled with 200-500 particles. The velocity and density methods have to use halos with at least 4000 particles in order to keep the biases down to the same low level. We also show that the mass-concentration relationship could be shallower than expected once the biases of the different concentration measurements are taken into account. These results show that bootstrapping and the concentration estimates based on the integrated mass profile are valuable tools to probe the internal structure of dark matter halos in numerical simulations.
In this paper, we consider a simple brane model with a generic dark energy component which could drive the accelerated expansion at late times of the Universe. We use the Supernovae type Ia, $H(z)$, baryon acoustic oscillations, and cosmic microwave background radiation measurements to constrain the brane tension, which is the main observable of the theory. From the study, we found an important tension between the different data sets and evidence of no gravity modifications by the existence of an extra dimension. Although this specific braneworld model is not compatible with the current cosmological observations and offers no new insights into the dark energy problem, it is not ruled out either. Our results show the need to further test of the braneworld model with appropriate correction terms.
Peculiar velocity surveys present a very promising route to measuring the growth rate of large-scale structure and its scale dependence. However, individual peculiar velocity surveys suffer from large statistical errors due to the intrinsic scatter in the relations used to infer a galaxy's true distance. In this context we use a Fisher Matrix formalism to investigate the statistical benefits of combining multiple peculiar velocity surveys. We find that for all cases we consider there is a marked improvement on constraints on the linear growth rate $f\sigma_{8}$. For example, the constraining power of only a few peculiar velocity measurements is such that the addition of the 2MASS Tully-Fisher survey (containing only $\sim2,000$ galaxies) to the full redshift and peculiar velocity samples of the 6-degree Field Galaxy Survey (containing $\sim 110,000$ redshifts and $\sim 9,000$ velocities) can improve growth rate constraints by $\sim20\%$. Furthermore, the combination of the future TAIPAN and WALLABY+WNSHS surveys has the potential to reach a $\sim3\%$ error on $f\sigma_{8}$, which will place tight limits on possible extensions to General Relativity. We then turn to look at potential systematics in growth rate measurements that can arise due to incorrect calibration of the peculiar velocity zero-point and from scale-dependent spatial and velocity bias. For next generation surveys, we find that neglecting velocity bias in particular has the potential to bias constraints on the growth rate by over $5\sigma$, but that an offset in the zero-point has negligible impact on the velocity power spectrum.
Galaxy clusters undergo mergers that can generate extended radio sources called radio relics. Radio relics are the consequence of merger-induced shocks that propagate in the intra cluster medium (ICM). In this paper we analyse the radio, optical and X-ray data from a candidate galaxy cluster that has been selected from the radio emission coming from a candidate radio relic detected in NRAO VLA Sky Survey (NVSS). Our aim is to clarify the nature of this source and prove that under certain conditions radio emission from radio relics can be used to trace relatively low-mass galaxy clusters. We have observed the candidate galaxy cluster with the Giant Meterwave Radio Telescope (GMRT) at three different frequencies. These datasets have been analysed together with archival data from ROSAT in the X-ray and with archival data from the Gamma-Ray Burst Optical/Near-Infrared Detector (GROND) telescope in four different optical bands. We confirm the presence of a 1 Mpc long radio relic located in the outskirts of a previously unknown galaxy cluster. We confirm the presence of the galaxy cluster through dedicated optical observations and using archival X-ray data. Due to its proximity and similar redshift to a known Abell cluster, we named it: Abell 3527-bis. The galaxy cluster is among the least massive cluster known to host a radio relic. We showed that radio relics can be effectively used to trace a subset of relatively low-mass galaxy clusters that might have gone undetected in X-ray or Sunyaev-Zel'dovich (SZ) surveys. This technique might be used in future deep, low-frequency surveys as those carried on by LOFAR, uGMRT and, ultimately, SKA.
We present a simple method for the identification of weak signals associated with gravitational wave events. Its application reveals a signal with the same time lag as the GW150914 event in the released LIGO strain data with a significance around $3.2\sigma$. This signal starts about 10 minutes before GW150914 and lasts for about 45 minutes. Subsequent tests suggest that this signal is likely to be due to external sources.
We use the Fisher matrix formalism to study the expansion and growth history of the Universe using galaxy clustering with 2D angular cross-correlation tomography in spectroscopic or high resolution photometric redshift surveys. The radial information is contained in the cross correlations between narrow redshift bins. We show how multiple tracers with redshift space distortions cancel sample variance and arbitrarily improve the constraints on the dark energy equation of state $\omega(z)$ and the growth parameter $\gamma$ in the noiseless limit. The improvement for multiple tracers quickly increases with the bias difference between the tracers, up to a factor $\sim4$ in $\text{FoM}_{\gamma\omega}$. We model a magnitude limited survey with realistic density and bias using a conditional luminosity function, finding a factor 1.3-9.0 improvement in $\text{FoM}_{\gamma\omega}$ -- depending on global density -- with a split in a halo mass proxy. Partly overlapping redshift bins improve the constraints in multiple tracer surveys a factor $\sim1.3$ in $\text{FoM}_{\gamma\omega}$. This findings also apply to photometric surveys, where the effect of using multiple tracers is magnified. We also show large improvement on the FoM with increasing density, which could be used as a trade-off to compensate some possible loss with radial resolution.
Understanding of the observed structure in the universe can be reached only in the theoretical framework of dark matter. N-body simulations are indispensable for the analysis of the formation and evolution of the dark matter web. Two primary fields - density and velocity fields - are used in most of studies. However dark matter provides two additional fields which are unique for collisionless media only. These are the multi- stream field in Eulerian space and flip-flop field in Lagrangian space. The flip-flop field represents the number of sign reversals of an elementary volume of each collisionless fluid element. This field can be estimated by counting the sign reversals of the Jacobian at each particle at every time step of the simulation. The Jacobian is evaluated by numerical differentiation of the Lagrangian submanifold, i.e., the three-dimensional dark matter sheet in the six-dimensional space formed by three Lagrangian and three Eulerian coordinates. We present the results of the statistical study of the evolution of the flip-flop field from z = 50 to the present time z = 0. A number of statistical characteristics show that the pattern of the flip-flop field remains remarkably stable from z = 30 to the present time. As a result the flip-flop field evaluated at z = 0 stores a wealth of information about the dynamical history of the dark matter web. In particular one of the most intriguing properties of the flip-flop is a unique capability to preserve the information about the merging history of dark matter haloes.
The motion of the solar system with respect to the cosmic rest frame modulates the monopole of the Epoch of Reionization 21-cm signal into a dipole. This dipole has a characteristic frequency dependence that is dominated by the frequency derivative of the monopole signal. We argue that although the signal is weaker by a factor of $\sim200$, there are significant benefits in measuring the dipole. Most importantly, the direction of the cosmic velocity vector is known exquisitely well from the cosmic microwave background and is not aligned with the galaxy velocity vector that modulates the foreground monopole. Moreover, an experiment designed to measure a dipole can rely on differencing patches of the sky rather than making an absolute signal measurement, which helps with some systematic effects.
We use recently published redshift space distortion measurements of the cosmological growth rate, f sigma_8(z), to examine whether the linear evolution of perturbations in the R_h=ct cosmology is consistent with the observed development of large scale structure. We find that these observations favour R_h=ct over the version of LCDM optimized with the joint analysis of Planck and linear growth rate data, particularly in the redshift range 0 < z < 1, where a significant curvature in the functional form of f sigma_8(z) predicted by the standard model---but not by R_h=ct---is absent in the data. When LCDM is optimized using solely the growth rate measurements, however, the two models fit the observations equally well though, in this case, the low-redshift measurements find a lower value for the fluctuation amplitude than is expected in Planck LCDM. Our results strongly affirm the need for more precise measurements of f sigma_8(z) at all redshifts, but especially at z < 1.
We propose a new mechanism to generate a lepton asymmetry based on the vacuum CP-violating phase transition (CPPT). This approach differs from classical thermal leptogenesis as a specific seesaw model, and its UV completion, need not be specified. The lepton asymmetry is generated via the dynamically realised coupling of the Weinberg operator during the phase transition. This mechanism provides strong connections with low-energy neutrino experiments.
We show that the effective theory describing single component continuous media with a linear and constant equation of state of the form $p=w\rho$ is invariant under a 1-parameter family of continuous disformal transformations. In the special case of $w=1/3$ (ultrarelativistic gas), such a family reduces to conformal transformations. As examples, perfect fluids, homogeneous and isotropic solids are discussed.
Sudden singularities occur in FRW spacetimes when the scale factor remains finite and different from zero while some of its derivatives diverge. After proper rescaling, the scale factor close to such a singularity at $t=0$ takes the form $a(t)=1+ c |t|^\eta$ (where $c$ and $\eta$ are parameters and $\eta\geq 0$). We investigate analytically and numerically the geodesics of free and gravitationally bound particles through such sudden singularities. We find that even though free particle geodesics go through sudden singularities for all $\eta\geq 0$, bound systems get dissociated for a wide range of the parameter $c$. For $\eta <1$ bound particles receive a diverging impulse at the singularity and get dissociated for all positive values of the parameter $c$. For $\eta > 1$ (Sudden Future Singularities (SFS)) bound systems get a finite impulse that depends on the value of $c$ and get dissociated for values of $c$ larger than a critical value $c_{cr}(\eta,\omega_0)>0$ that increases with the value of $\eta$ and the rescaled angular velocity $\omega_0$ of the bound system. We obtain an approximate equation for the analytical estimate of $c_{cr}(\eta,\omega_0)$. We also obtain its accurate form by numerical derivation of the bound system orbits through the singularities. Bound system orbits through Big Brake singularities ($c<0$, $1<\eta<2$) are also derived numerically and are found to get disrupted at the singularity. However, they remain bound for all values of the parameter $c$ considered.
Motivated by the projectable Horava-Lifshitz model/mimetic matter scenario, we consider a particular modification of standard gravity, which manifests as an imperfect low pressure fluid. While practically indistinguishable from collection of non-relativistic weakly interacting particles on cosmological scales, it leaves drastically different signatures in the Solar system. The main effect stems from gravitational focusing of the flow of {\it Imperfect Dark Matter} passing near the Sun. This entails the strong amplification of Imperfect Dark Matter energy density compared to its average value in the surrounding halo. The enhancement is many orders of magnitude larger than in the case of Cold Dark Matter, provoking deviations of the metric in the second order in the Newtonian potential. Effects of gravitational focusing are prominent enough to substantially affect the planetary dynamics. Using the existing bound on the PPN parameter $\beta_{PPN}$, we deduce the stringent constraint on the unique constant of the model.
Links to: arXiv, form interface, find, astro-ph, recent, 1609, contact, help (Access key information)
We probe the higher-order clustering of the galaxies in the final data release (DR12) of the Sloan Digital Sky Survey Baryon Oscillation Spectroscopic Survey (BOSS) using the method of germ-grain Minkowski Functionals (MFs). Our sample consists of 410,615 BOSS galaxies from the northern Galactic cap in the redshift range 0.450--0.595. We show the MFs to be sensitive to contributions up to the six-point correlation function for this data set. We ensure with a custom angular mask that the results are more independent of boundary effects than in previous analyses of this type. We extract the higher-order part of the MFs and quantify the difference to the case without higher-order correlations. The resulting $\chi^{2}$ value of over 10,000 for a modest number of degrees of freedom, O(200), indicates a 100-sigma deviation and demonstrates that we have a highly significant signal of the non-Gaussian contributions to the galaxy distribution. This statistical power can be useful in testing models with differing higher-order correlations. Comparing the galaxy data to the QPM and MultiDark-Patchy mocks, we find that the latter better describes the observed structure. From an order-by-order decomposition we expect that, for example, already a reduction of the amplitude of the MD-Patchy mock power spectrum by 5% would remove the remaining tension.
The outcome of upcoming cosmological surveys will depend on the accurate estimates of photometric redshifts. In the framework of the implementation of the photo-z algorithm for Euclid, we are exploring new avenues to improve template-fitting methods. The paper focusses on the prescription of the extinction of source light by dust in the Milky Way. Since Galactic extinction strongly correlates with wavelength and photometry is commonly obtained in broad-band filters, the amount of absorption depends on the source SED, a point often neglected as the SED is not known a-priori. A consequence of this is that the observed E(B-V) (=A_B-A_V) will be different from the E(B-V) used to normalise the absorption law k_lambda (=A_lambda/E(B-V)). Band-pass corrections are required to renormalise the law for a given SED. We assess the band-pass corrections of a range of SEDs and find they vary by up to 20%. We investigate how dust-to-reddening scaling factors depend of the sources used for their calibration. We derive scaling factors from the color excesses of z<0.4 SDSS red galaxies and show that band-pass corrections predict the observed differences. Extinction is then estimated for a range of SEDs and filters relevant to Euclid and other cosmological ground-based surveys. For high extinction line-of-sights (E(B-V)>0.1, ~8% of the Euclid survey), the variations in corrections can be ~0.1mag in the `bluer' optical filters and ~0.04mag in the NIR filters. An inaccurate correction of extinction critically affects photo-z. In particular, for high extinctions and z<0.5, the bias (mean D_z=z_phot-z_real) exceeds 0.2%(1+z), the precision required by weak-lensing analyses. Additional uncertainty on the MW extinction law further reduces the photo-z precision. We propose a new prescription of Galactic absorption for template-fitting algorithms that takes into consideration the dependence of extinction with SED.
Directionally sensitive dark matter (DM) direct detection experiments present the only way to observe the full three-dimensional velocity distribution of the Milky Way halo local to Earth. In this work we compare methods for extracting information about the local DM velocity distribution from a set of recoil directions and energies in a range of hypothetical directional and non-directional experiments. We compare a model independent empirical parameterisation of the velocity distribution based on an angular discretisation with a model dependent approach which assumes knowledge of the functional form of the distribution. The methods are tested under three distinct halo models which cover a range of possible phase space structures for the local velocity distribution: a smooth Maxwellian halo, a tidal stream and a debris flow. In each case we use simulated directional data to attempt to reconstruct the shape and parameters describing each model as well as the DM particle properties. We find that the empirical parametrisation is able to make accurate unbiased reconstructions of the DM mass and cross section as well as capture features in the underlying velocity distribution in certain directions without any assumptions about its true functional form. We also find that by extracting directionally averaged velocity parameters with this method one can discriminate between halo models with different classes of substructure.
Future large scale structure surveys will provide increasingly tight constraints on our cosmological model. These surveys will report results on the distance scale and growth rate of perturbations through measurements of Baryon Acoustic Oscillations and Redshift-Space Distortions. It is interesting to ask: what further analyses should become routine, so as to test as-yet-unknown models of cosmic acceleration? Models which aim to explain the accelerated expansion rate of the Universe by modifications to General Relativity often invoke screening mechanisms which can imprint a non-standard density dependence on their predictions. This suggests density-dependent clustering as a `generic' constraint. This paper argues that a density-marked correlation function provides a density-dependent statistic which is easy to compute and report and requires minimal additional infrastructure beyond what is routinely available to such survey analyses. We give one realization of this idea and study it using low order perturbation theory. We encourage groups developing modified gravity theories to see whether such statistics provide discriminatory power for their models.
One alternative to the CDM paradigm is the Scalar Field Dark Matter (SFDM) model, which assumes dark matter is a spin-0 ultra-light scalar field with a typical mass $m\sim10^{-22}\mathrm{eV}/c^2$ and positive self-interactions. Due to the ultra-light boson mass, the SFDM could form Bose-Einstein condensates in the very early universe, which are interpreted as the dark matter haloes. Although cosmologically the model behaves as CDM, they differ at small scales: SFDM naturally predicts fewer satellite haloes, cores in dwarf galaxies and the formation of massive galaxies at high redshifts. The ground state (or BEC) solution at zero temperature suffices to describe low-mass galaxies but fails for larger systems. A possible solution is adding finite-temperature corrections to the SF potential which allows combinations of excited states. In this work we test the finite-temperature multistate SFDM solution at galaxy cluster scales and compare our results with the NFW and BEC profiles. We achieve this by fitting the mass distribution of 13 Chandra X-ray clusters of galaxies, excluding the brightest galaxy central region. We show that the SFDM model accurately describes the clusters' DM mass distributions offering an equivalent or better agreement than the NFW profile. The complete disagreement of the BEC model with the data is also shown. We conclude that the theoretically motivated multistate SFDM profile is an interesting alternative to empirical profiles and \textit{ad hoc} fitting-functions that attempt to couple the asymptotic NFW decline with the core SFDM model.
In this paper, based on a 2.29 GHz VLBI all-sky survey of 613 milliarcsecond ultra-compact radio sources with $0.0035<z<3.787$, we describe a method of identifing the sub-sample which can serve as individual standard rulers in cosmology. If the linear size of the compact structure is assumed to depend on source luminosity and redshift as $l_m=l L^\beta (1+z)^n$, only intermediate-luminosity quasars ($10^{27}$ W/Hz$<L<$ $10^{28}$ W/Hz) show negligible dependence ($|n|\simeq 10^{-3}$, $|\beta|\simeq 10^{-4}$), and thus represent a population of such rulers with fixed characteristic length $l=11.42$ pc. With a sample of 120 such sources covering the redshift range $0.46<z<2.80$, we confirm the existence of dark energy in the Universe with high significance, and obtain stringent constraints on both the matter density $\Omega_m=0.323\pm0.195$ and the Hubble constant $H_0=66.25\pm7.75$ km sec$^{-1}$ Mpc$^{-1}$. Finally, with the angular diameter distances $D_A$ measured for quasars extending to high redshifts ($z\sim 3.0$), we reconstruct the $D_A(z)$ function using the technique of Gaussian processes. This allows us to identify the redshift corresponding to the maximum of the $D_A(z)$ function: $z_m=1.65$ and the corresponding angular diameter distance $D_A(z_m)=1732.73\pm74.66$ Mpc. Similar reconstruction of the expansion rate function $H(z)$ based on the data from cosmic chronometers and BAO gives us $H(z_m)=172.85\pm6.06$ km sec$^{-1}$ Mpc$^{-1}$. These measurements are used to estimate the speed of light at this redshift: $c=2.995(\pm0.235)\times 10^5$ km/s. This is the first measurement of the speed of light in a cosmological setting referring to the distant past.
The prospects of future galaxy surveys for non-Gaussianity measurements call for the development of robust techniques for computing the bispectrum of primordial cosmological perturbations. In this paper, we propose a novel approach to the calculation of the squeezed bispectrum in multiple-field inflation. With use of the $\delta N$ formalism, our framework sheds new light on the recently pointed out difference between the squeezed bispectrum for global observers and that for local observers, while allowing one to calculate both. For local observers in particular, the squeezed bispectrum is found to vanish in single-field inflation. Furthermore, our framework allows one to go beyond the near-equilateral ("small hierarchy") limit, and to automatically include intrinsic non-Gaussianities that do not need to be calculated separately. The explicit computational programme of our method is given and illustrated with a few examples.
We present new measurements of the power spectra of the cosmic infrared background (CIB) anisotropies using the Planck 2015 full-mission HFI data at 353, 545, and 857 GHz over 20000 square degrees. We use techniques similar to those applied for the cosmological analysis of Planck, subtracting dust emission at the power spectrum level. Our analysis gives stable solutions for the CIB power spectra with increasing sky coverage up to about 50% of the sky. These spectra agree well with Hi cleaned spectra from Planck measured on much smaller areas of sky with low Galactic dust emission. At 545 and 857 GHz our CIB spectra agree well with those measured from Herschel data. We find that the CIB spectra at l > 500 are well fitted by a power-law model for the clustered CIB, with a shallow index {\gamma}^cib = 0.53\pm0.02. This is consistent with the CIB results at 217 GHz from the cosmological parameter analysis of Planck. We show that a linear combination of the 545 and 857 GHz Planck maps is dominated by CIB fluctuations at multipoles l > 300.
In this paper, we consider two case examples of Dirac-Born-Infeld (DBI) generalizations of canonical large-field inflation models, characterized by a reduced sound speed, $c_{S} < 1$. The reduced speed of sound lowers the tensor-scalar ratio, improving the fit of the models to the data, but increases the equilateral-mode non-Gaussianity, $f^\mathrm{equil.}_\mathrm{NL}$, which the latest results from the Planck satellite constrain by a new upper bound. We examine constraints on these models in light of the most recent Planck and BICEP/Keck results, and find that they have a greatly decreased window of viability. The upper bound on $f^\mathrm{equil.}_\mathrm{NL}$ corresponds to a lower bound on the sound speed and a corresponding lower bound on the tensor-scalar ratio of $r \sim 0.01$, so that near-future Cosmic Microwave Background observations may be capable of ruling out entire classes of DBI inflation models. The result is, however, not universal: infrared-type DBI inflation models, where the speed of sound increases with time, are not subject to the bound.
I calculate the rate of WIMP capture and annihilation in the Earth in the non-relativistic effective theory of dark matter-nucleon interactions. Neglecting operator interference, I consider all Galilean invariant interaction operators that can arise from the exchange of a heavy particle of spin less than or equal to one when WIMPs have spin 0, 1/2 or 1. I compute position and shape of the expected resonances in the mass - capture rate plane and show that Iron is not the most important element in the capture process for many currently ignored interaction operators. I compare these predictions with the recent results of an Earth WIMP analysis of IceCube in the 86-string configuration and set limits on all isoscalar and isovector coupling constants of the effective theory of dark matter-nucleon interactions. For certain interaction operators and for a dark matter particle mass of about 50 GeV, I find that these limits are stronger than those I have previously derived in an analysis of the solar WIMP search performed at IceCube in the 79-string configuration.
Astrophysical techniques have pioneered the discovery of neutrino mass properties. Current cosmological observations give an upper bound on neutrino masses by attempting to disentangle the small neutrino contribution from the sum of all matter using precise theoretical models. We discover the differential neutrino condensation effect in our TianNu N-body simulation. Neutrino masses can be inferred using this effect by comparing galaxy properties in regions of the universe with different neutrino relative abundance (i.e. the local neutrino to cold dark matter density ratio). In "neutrino-rich"' regions, more neutrinos can be captured by massive halos compared to "neutrino-poor" regions. This effect differentially skews the halo mass function and opens up the path to independent neutrino mass measurements in current or future galaxy surveys.
We study non--perturbatively the effect of the deflection angle on the BAO wiggles of the matter power spectrum in real space. We show that from redshift z~2 this introduces a dispersion of roughly 1 Mpc at BAO scale, which corresponds approximately to a 1% effect. The lensing effect induced by the deflection angle, which is completely geometrical and survey independent, smears out the BAO wiggles. The effect on the power spectrum amplitude at BAO scale is about 0.1% for z~2 and 0.2% for z~4. We compare the smoothing effects induced by the lensing potential and non--linear structure formation, showing that the two effects become comparable at z~4, while the lensing effect dominates for sources at higher redshifts. We note that this effect is not accounted through BAO reconstruction techniques.
The number of nonrelativistic axions can be changed by inelastic reactions that produce photons or relativistic axions. Any odd number of axions can annihilate into two photons. Any even number of nonrelativistic axions can scatter into two relativistic axions. We calculate the rate at which axions are lost from axion stars from these inelastic reactions. In dilute systems of axions, the dominant inelastic reaction is axion decay into two photons. In sufficiently dense systems of axions, the dominant inelastic reaction is the scattering of four nonrelativistic axions into two relativistic axions. The scattering of odd numbers of axions into two photons produces monochromatic radio-frequency signals at odd-integer harmonics of the fundamental frequency set by the axion mass. This provides a unique signature for dense systems of axions, such as a dense axion star or a collapsing dilute axion star.
We report on the successful completion of a 2 trillion particle cosmological simulation to z=0 run on the Piz Daint supercomputer (CSCS, Switzerland), using 4000+ GPU nodes for a little less than 80h of wall-clock time or 350,000 node hours. Using multiple benchmarks and performance measurements on the US Oak Ridge National Laboratory Titan supercomputer, we demonstrate that our code PKDGRAV3, delivers, to our knowledge, the fastest time-to-solution for large-scale cosmological N-body simulations. This was made possible by using the Fast Multipole Method in conjunction with individual and adaptive particle time steps, both deployed efficiently (and for the first time) on supercomputers with GPU-accelerated nodes. The very low memory footprint of PKDGRAV3 allowed us to run the first ever benchmark with 8 trillion particles on Titan, and to achieve perfect scaling up to 18000 nodes and a peak performance of 10 Pflops.
We study the dust content of galaxies from z $=$ 0 to z $=$ 9 in semi-analytic models of galaxy formation that include new recipes to track the production and destruction of dust. We include condensation of dust in stellar ejecta, the growth of dust in the interstellar medium (ISM), the destruction of dust by supernovae and in the hot halo, and dusty winds and inflows. The rate of dust growth in the ISM depends on the metallicity and density of molecular clouds. Our fiducial model reproduces the relation between dust mass and stellar mass from z $=$ 0 to z $=$ 7, the dust-to-gas ratio of local galaxies as a function of stellar mass, the double power law trend between dust-to- gas ratio and gas-phase metallicity, the number density of galaxies with dust masses less than $10^{8.3} M_\odot$, and the cosmic density of dust at z $=$ 0. The dominant mode of dust formation is dust growth in the ISM, except for galaxies with $M_* < 10^7 M_\odot$, where condensation of dust in supernova ejecta dominates. The dust-to-metal ratio of galaxies evolves as a function of gas-phase metallicity, unlike what is typically assumed in cosmological simulations. Model variants including higher condensation efficiencies, a fixed timescale for dust growth in the ISM, or no growth at all reproduce some of the observed constraints, but fail to reproduce the shape of dust scaling relations and the dust mass of high-redshift galaxies simultaneously.
We investigate the effects of dense environments on galaxy evolution by examining how the properties of galaxies in the z = 1.6 protocluster Cl 0218.3-0510 depend on their location. We determine galaxy properties using spectral energy distribution fitting to 14-band photometry, including data at three wavelengths that tightly bracket the Balmer and 4000A breaks of the protocluster galaxies. We find that two-thirds of the protocluster galaxies, which lie between several compact groups, are indistinguishable from field galaxies. The other third, which reside within the groups, differ significantly from the intergroup galaxies in both colour and specific star formation rate. We find that the fraction of red galaxies within the massive protocluster groups is twice that of the intergroup region. These excess red galaxies are due to enhanced fractions of both passive galaxies (1.7 times that of the intergroup region) and dusty star-forming galaxies (3 times that of the intergroup region). We infer that some protocluster galaxies are processed in the groups before the cluster collapses. These processes act to suppress star formation and change the mode of star formation from unobscured to obscured.
SDSS J2222+2745 is a galaxy cluster at z=0.49, strongly lensing a quasar at z=2.805 into six widely separated images. In recent HST imaging of the field, we identify additional multiply lensed galaxies, and confirm the sixth quasar image that was identified by Dahle et al. (2013). We used the Gemini North telescope to measure a spectroscopic redshift of z=4.56 of one of the secondary lensed galaxies. These data are used to refine the lens model of SDSS J2222+2745, compute the time delay and magnifications of the lensed quasar images, and reconstruct the source image of the quasar host and a second lensed galaxy at z=2.3. This second galaxy also appears in absorption in our Gemini spectra of the lensed quasar, at a projected distance of 34 kpc. Our model is in agreement with the recent time delay measurements of Dahle et al. (2015), who found tAB=47.7+/-6.0 days and tAC=-722+/-24 days. We use the observed time delays to further constrain the model, and find that the model-predicted time delays of the three faint images of the quasar are tAD=502+/-68 days, tAE=611+/-75 days, and tAF=415+/-72 days. We have initiated a follow-up campaign to measure these time delays with Gemini North. Finally, we present initial results from an X-ray monitoring program with Swift, indicating the presence of hard X-ray emission from the lensed quasar, as well as extended X-ray emission from the cluster itself, which is consistent with the lensing mass measurement and the cluster velocity dispersion.
Links to: arXiv, form interface, find, astro-ph, recent, 1609, contact, help (Access key information)
We probe the higher-order clustering of the galaxies in the final data release (DR12) of the Sloan Digital Sky Survey Baryon Oscillation Spectroscopic Survey (BOSS) using the method of germ-grain Minkowski Functionals (MFs). Our sample consists of 410,615 BOSS galaxies from the northern Galactic cap in the redshift range 0.450--0.595. We show the MFs to be sensitive to contributions up to the six-point correlation function for this data set. We ensure with a custom angular mask that the results are more independent of boundary effects than in previous analyses of this type. We extract the higher-order part of the MFs and quantify the difference to the case without higher-order correlations. The resulting $\chi^{2}$ value of over 10,000 for a modest number of degrees of freedom, O(200), indicates a 100-sigma deviation and demonstrates that we have a highly significant signal of the non-Gaussian contributions to the galaxy distribution. This statistical power can be useful in testing models with differing higher-order correlations. Comparing the galaxy data to the QPM and MultiDark-Patchy mocks, we find that the latter better describes the observed structure. From an order-by-order decomposition we expect that, for example, already a reduction of the amplitude of the MD-Patchy mock power spectrum by 5% would remove the remaining tension.
The outcome of upcoming cosmological surveys will depend on the accurate estimates of photometric redshifts. In the framework of the implementation of the photo-z algorithm for Euclid, we are exploring new avenues to improve template-fitting methods. The paper focusses on the prescription of the extinction of source light by dust in the Milky Way. Since Galactic extinction strongly correlates with wavelength and photometry is commonly obtained in broad-band filters, the amount of absorption depends on the source SED, a point often neglected as the SED is not known a-priori. A consequence of this is that the observed E(B-V) (=A_B-A_V) will be different from the E(B-V) used to normalise the absorption law k_lambda (=A_lambda/E(B-V)). Band-pass corrections are required to renormalise the law for a given SED. We assess the band-pass corrections of a range of SEDs and find they vary by up to 20%. We investigate how dust-to-reddening scaling factors depend of the sources used for their calibration. We derive scaling factors from the color excesses of z<0.4 SDSS red galaxies and show that band-pass corrections predict the observed differences. Extinction is then estimated for a range of SEDs and filters relevant to Euclid and other cosmological ground-based surveys. For high extinction line-of-sights (E(B-V)>0.1, ~8% of the Euclid survey), the variations in corrections can be ~0.1mag in the `bluer' optical filters and ~0.04mag in the NIR filters. An inaccurate correction of extinction critically affects photo-z. In particular, for high extinctions and z<0.5, the bias (mean D_z=z_phot-z_real) exceeds 0.2%(1+z), the precision required by weak-lensing analyses. Additional uncertainty on the MW extinction law further reduces the photo-z precision. We propose a new prescription of Galactic absorption for template-fitting algorithms that takes into consideration the dependence of extinction with SED.
Directionally sensitive dark matter (DM) direct detection experiments present the only way to observe the full three-dimensional velocity distribution of the Milky Way halo local to Earth. In this work we compare methods for extracting information about the local DM velocity distribution from a set of recoil directions and energies in a range of hypothetical directional and non-directional experiments. We compare a model independent empirical parameterisation of the velocity distribution based on an angular discretisation with a model dependent approach which assumes knowledge of the functional form of the distribution. The methods are tested under three distinct halo models which cover a range of possible phase space structures for the local velocity distribution: a smooth Maxwellian halo, a tidal stream and a debris flow. In each case we use simulated directional data to attempt to reconstruct the shape and parameters describing each model as well as the DM particle properties. We find that the empirical parametrisation is able to make accurate unbiased reconstructions of the DM mass and cross section as well as capture features in the underlying velocity distribution in certain directions without any assumptions about its true functional form. We also find that by extracting directionally averaged velocity parameters with this method one can discriminate between halo models with different classes of substructure.
Future large scale structure surveys will provide increasingly tight constraints on our cosmological model. These surveys will report results on the distance scale and growth rate of perturbations through measurements of Baryon Acoustic Oscillations and Redshift-Space Distortions. It is interesting to ask: what further analyses should become routine, so as to test as-yet-unknown models of cosmic acceleration? Models which aim to explain the accelerated expansion rate of the Universe by modifications to General Relativity often invoke screening mechanisms which can imprint a non-standard density dependence on their predictions. This suggests density-dependent clustering as a `generic' constraint. This paper argues that a density-marked correlation function provides a density-dependent statistic which is easy to compute and report and requires minimal additional infrastructure beyond what is routinely available to such survey analyses. We give one realization of this idea and study it using low order perturbation theory. We encourage groups developing modified gravity theories to see whether such statistics provide discriminatory power for their models.
One alternative to the CDM paradigm is the Scalar Field Dark Matter (SFDM) model, which assumes dark matter is a spin-0 ultra-light scalar field with a typical mass $m\sim10^{-22}\mathrm{eV}/c^2$ and positive self-interactions. Due to the ultra-light boson mass, the SFDM could form Bose-Einstein condensates in the very early universe, which are interpreted as the dark matter haloes. Although cosmologically the model behaves as CDM, they differ at small scales: SFDM naturally predicts fewer satellite haloes, cores in dwarf galaxies and the formation of massive galaxies at high redshifts. The ground state (or BEC) solution at zero temperature suffices to describe low-mass galaxies but fails for larger systems. A possible solution is adding finite-temperature corrections to the SF potential which allows combinations of excited states. In this work we test the finite-temperature multistate SFDM solution at galaxy cluster scales and compare our results with the NFW and BEC profiles. We achieve this by fitting the mass distribution of 13 Chandra X-ray clusters of galaxies, excluding the brightest galaxy central region. We show that the SFDM model accurately describes the clusters' DM mass distributions offering an equivalent or better agreement than the NFW profile. The complete disagreement of the BEC model with the data is also shown. We conclude that the theoretically motivated multistate SFDM profile is an interesting alternative to empirical profiles and \textit{ad hoc} fitting-functions that attempt to couple the asymptotic NFW decline with the core SFDM model.
In this paper, based on a 2.29 GHz VLBI all-sky survey of 613 milliarcsecond ultra-compact radio sources with $0.0035<z<3.787$, we describe a method of identifing the sub-sample which can serve as individual standard rulers in cosmology. If the linear size of the compact structure is assumed to depend on source luminosity and redshift as $l_m=l L^\beta (1+z)^n$, only intermediate-luminosity quasars ($10^{27}$ W/Hz$<L<$ $10^{28}$ W/Hz) show negligible dependence ($|n|\simeq 10^{-3}$, $|\beta|\simeq 10^{-4}$), and thus represent a population of such rulers with fixed characteristic length $l=11.42$ pc. With a sample of 120 such sources covering the redshift range $0.46<z<2.80$, we confirm the existence of dark energy in the Universe with high significance, and obtain stringent constraints on both the matter density $\Omega_m=0.323\pm0.195$ and the Hubble constant $H_0=66.25\pm7.75$ km sec$^{-1}$ Mpc$^{-1}$. Finally, with the angular diameter distances $D_A$ measured for quasars extending to high redshifts ($z\sim 3.0$), we reconstruct the $D_A(z)$ function using the technique of Gaussian processes. This allows us to identify the redshift corresponding to the maximum of the $D_A(z)$ function: $z_m=1.65$ and the corresponding angular diameter distance $D_A(z_m)=1732.73\pm74.66$ Mpc. Similar reconstruction of the expansion rate function $H(z)$ based on the data from cosmic chronometers and BAO gives us $H(z_m)=172.85\pm6.06$ km sec$^{-1}$ Mpc$^{-1}$. These measurements are used to estimate the speed of light at this redshift: $c=2.995(\pm0.235)\times 10^5$ km/s. This is the first measurement of the speed of light in a cosmological setting referring to the distant past.
The prospects of future galaxy surveys for non-Gaussianity measurements call for the development of robust techniques for computing the bispectrum of primordial cosmological perturbations. In this paper, we propose a novel approach to the calculation of the squeezed bispectrum in multiple-field inflation. With use of the $\delta N$ formalism, our framework sheds new light on the recently pointed out difference between the squeezed bispectrum for global observers and that for local observers, while allowing one to calculate both. For local observers in particular, the squeezed bispectrum is found to vanish in single-field inflation. Furthermore, our framework allows one to go beyond the near-equilateral ("small hierarchy") limit, and to automatically include intrinsic non-Gaussianities that do not need to be calculated separately. The explicit computational programme of our method is given and illustrated with a few examples.
We present new measurements of the power spectra of the cosmic infrared background (CIB) anisotropies using the Planck 2015 full-mission HFI data at 353, 545, and 857 GHz over 20000 square degrees. We use techniques similar to those applied for the cosmological analysis of Planck, subtracting dust emission at the power spectrum level. Our analysis gives stable solutions for the CIB power spectra with increasing sky coverage up to about 50% of the sky. These spectra agree well with Hi cleaned spectra from Planck measured on much smaller areas of sky with low Galactic dust emission. At 545 and 857 GHz our CIB spectra agree well with those measured from Herschel data. We find that the CIB spectra at l > 500 are well fitted by a power-law model for the clustered CIB, with a shallow index {\gamma}^cib = 0.53\pm0.02. This is consistent with the CIB results at 217 GHz from the cosmological parameter analysis of Planck. We show that a linear combination of the 545 and 857 GHz Planck maps is dominated by CIB fluctuations at multipoles l > 300.
In this paper, we consider two case examples of Dirac-Born-Infeld (DBI) generalizations of canonical large-field inflation models, characterized by a reduced sound speed, $c_{S} < 1$. The reduced speed of sound lowers the tensor-scalar ratio, improving the fit of the models to the data, but increases the equilateral-mode non-Gaussianity, $f^\mathrm{equil.}_\mathrm{NL}$, which the latest results from the Planck satellite constrain by a new upper bound. We examine constraints on these models in light of the most recent Planck and BICEP/Keck results, and find that they have a greatly decreased window of viability. The upper bound on $f^\mathrm{equil.}_\mathrm{NL}$ corresponds to a lower bound on the sound speed and a corresponding lower bound on the tensor-scalar ratio of $r \sim 0.01$, so that near-future Cosmic Microwave Background observations may be capable of ruling out entire classes of DBI inflation models. The result is, however, not universal: infrared-type DBI inflation models, where the speed of sound increases with time, are not subject to the bound.
I calculate the rate of WIMP capture and annihilation in the Earth in the non-relativistic effective theory of dark matter-nucleon interactions. Neglecting operator interference, I consider all Galilean invariant interaction operators that can arise from the exchange of a heavy particle of spin less than or equal to one when WIMPs have spin 0, 1/2 or 1. I compute position and shape of the expected resonances in the mass - capture rate plane and show that Iron is not the most important element in the capture process for many currently ignored interaction operators. I compare these predictions with the recent results of an Earth WIMP analysis of IceCube in the 86-string configuration and set limits on all isoscalar and isovector coupling constants of the effective theory of dark matter-nucleon interactions. For certain interaction operators and for a dark matter particle mass of about 50 GeV, I find that these limits are stronger than those I have previously derived in an analysis of the solar WIMP search performed at IceCube in the 79-string configuration.
Astrophysical techniques have pioneered the discovery of neutrino mass properties. Current cosmological observations give an upper bound on neutrino masses by attempting to disentangle the small neutrino contribution from the sum of all matter using precise theoretical models. We discover the differential neutrino condensation effect in our TianNu N-body simulation. Neutrino masses can be inferred using this effect by comparing galaxy properties in regions of the universe with different neutrino relative abundance (i.e. the local neutrino to cold dark matter density ratio). In "neutrino-rich"' regions, more neutrinos can be captured by massive halos compared to "neutrino-poor" regions. This effect differentially skews the halo mass function and opens up the path to independent neutrino mass measurements in current or future galaxy surveys.
We study non--perturbatively the effect of the deflection angle on the BAO wiggles of the matter power spectrum in real space. We show that from redshift z~2 this introduces a dispersion of roughly 1 Mpc at BAO scale, which corresponds approximately to a 1% effect. The lensing effect induced by the deflection angle, which is completely geometrical and survey independent, smears out the BAO wiggles. The effect on the power spectrum amplitude at BAO scale is about 0.1% for z~2 and 0.2% for z~4. We compare the smoothing effects induced by the lensing potential and non--linear structure formation, showing that the two effects become comparable at z~4, while the lensing effect dominates for sources at higher redshifts. We note that this effect is not accounted through BAO reconstruction techniques.
The number of nonrelativistic axions can be changed by inelastic reactions that produce photons or relativistic axions. Any odd number of axions can annihilate into two photons. Any even number of nonrelativistic axions can scatter into two relativistic axions. We calculate the rate at which axions are lost from axion stars from these inelastic reactions. In dilute systems of axions, the dominant inelastic reaction is axion decay into two photons. In sufficiently dense systems of axions, the dominant inelastic reaction is the scattering of four nonrelativistic axions into two relativistic axions. The scattering of odd numbers of axions into two photons produces monochromatic radio-frequency signals at odd-integer harmonics of the fundamental frequency set by the axion mass. This provides a unique signature for dense systems of axions, such as a dense axion star or a collapsing dilute axion star.
We report on the successful completion of a 2 trillion particle cosmological simulation to z=0 run on the Piz Daint supercomputer (CSCS, Switzerland), using 4000+ GPU nodes for a little less than 80h of wall-clock time or 350,000 node hours. Using multiple benchmarks and performance measurements on the US Oak Ridge National Laboratory Titan supercomputer, we demonstrate that our code PKDGRAV3, delivers, to our knowledge, the fastest time-to-solution for large-scale cosmological N-body simulations. This was made possible by using the Fast Multipole Method in conjunction with individual and adaptive particle time steps, both deployed efficiently (and for the first time) on supercomputers with GPU-accelerated nodes. The very low memory footprint of PKDGRAV3 allowed us to run the first ever benchmark with 8 trillion particles on Titan, and to achieve perfect scaling up to 18000 nodes and a peak performance of 10 Pflops.
We study the dust content of galaxies from z $=$ 0 to z $=$ 9 in semi-analytic models of galaxy formation that include new recipes to track the production and destruction of dust. We include condensation of dust in stellar ejecta, the growth of dust in the interstellar medium (ISM), the destruction of dust by supernovae and in the hot halo, and dusty winds and inflows. The rate of dust growth in the ISM depends on the metallicity and density of molecular clouds. Our fiducial model reproduces the relation between dust mass and stellar mass from z $=$ 0 to z $=$ 7, the dust-to-gas ratio of local galaxies as a function of stellar mass, the double power law trend between dust-to- gas ratio and gas-phase metallicity, the number density of galaxies with dust masses less than $10^{8.3} M_\odot$, and the cosmic density of dust at z $=$ 0. The dominant mode of dust formation is dust growth in the ISM, except for galaxies with $M_* < 10^7 M_\odot$, where condensation of dust in supernova ejecta dominates. The dust-to-metal ratio of galaxies evolves as a function of gas-phase metallicity, unlike what is typically assumed in cosmological simulations. Model variants including higher condensation efficiencies, a fixed timescale for dust growth in the ISM, or no growth at all reproduce some of the observed constraints, but fail to reproduce the shape of dust scaling relations and the dust mass of high-redshift galaxies simultaneously.
We investigate the effects of dense environments on galaxy evolution by examining how the properties of galaxies in the z = 1.6 protocluster Cl 0218.3-0510 depend on their location. We determine galaxy properties using spectral energy distribution fitting to 14-band photometry, including data at three wavelengths that tightly bracket the Balmer and 4000A breaks of the protocluster galaxies. We find that two-thirds of the protocluster galaxies, which lie between several compact groups, are indistinguishable from field galaxies. The other third, which reside within the groups, differ significantly from the intergroup galaxies in both colour and specific star formation rate. We find that the fraction of red galaxies within the massive protocluster groups is twice that of the intergroup region. These excess red galaxies are due to enhanced fractions of both passive galaxies (1.7 times that of the intergroup region) and dusty star-forming galaxies (3 times that of the intergroup region). We infer that some protocluster galaxies are processed in the groups before the cluster collapses. These processes act to suppress star formation and change the mode of star formation from unobscured to obscured.
SDSS J2222+2745 is a galaxy cluster at z=0.49, strongly lensing a quasar at z=2.805 into six widely separated images. In recent HST imaging of the field, we identify additional multiply lensed galaxies, and confirm the sixth quasar image that was identified by Dahle et al. (2013). We used the Gemini North telescope to measure a spectroscopic redshift of z=4.56 of one of the secondary lensed galaxies. These data are used to refine the lens model of SDSS J2222+2745, compute the time delay and magnifications of the lensed quasar images, and reconstruct the source image of the quasar host and a second lensed galaxy at z=2.3. This second galaxy also appears in absorption in our Gemini spectra of the lensed quasar, at a projected distance of 34 kpc. Our model is in agreement with the recent time delay measurements of Dahle et al. (2015), who found tAB=47.7+/-6.0 days and tAC=-722+/-24 days. We use the observed time delays to further constrain the model, and find that the model-predicted time delays of the three faint images of the quasar are tAD=502+/-68 days, tAE=611+/-75 days, and tAF=415+/-72 days. We have initiated a follow-up campaign to measure these time delays with Gemini North. Finally, we present initial results from an X-ray monitoring program with Swift, indicating the presence of hard X-ray emission from the lensed quasar, as well as extended X-ray emission from the cluster itself, which is consistent with the lensing mass measurement and the cluster velocity dispersion.
Links to: arXiv, form interface, find, astro-ph, recent, 1609, contact, help (Access key information)