Links to: arXiv, form interface, find, astro-ph, recent, 1712, contact, help (Access key information)
Ultracompact minihalos (UCMHs) have emerged as a valuable probe of the primordial power spectrum of density fluctuations at small scales. UCMHs are expected to form at early times in regions with $\delta\rho/\rho \gtrsim 10^{-3}$, and they are theorized to possess an extremely compact $\rho\propto r^{-9/4}$ radial density profile, which enhances their observable signatures. Non-observation of UCMHs can thus constrain the primordial power spectrum. Using N-body simulations to study the collapse of extreme density peaks at $z \simeq 1000$, we show that UCMHs forming under realistic conditions do not develop the $\rho\propto r^{-9/4}$ profile, and instead develop either $\rho\propto r^{-3/2}$ or $\rho\propto r^{-1}$ inner density profiles depending on the shape of the power spectrum. We also demonstrate via idealized simulations that self-similarity -- the absence of a scale length -- is necessary to produce a halo with the $\rho\propto r^{-9/4}$ profile, and we argue that this implies such halos cannot form from a Gaussian primordial density field. Prior constraints derived from UCMH non-observation must be reworked in light of this discovery. Although the shallower density profile reduces UCMH visibility, our findings reduce their signal by as little as $\mathcal O(10^{-2})$ while allowing later-forming halos to be considered, which suggests that new constraints could be significantly stronger.
By combining cosmological probes at low, intermediate and high redshifts, we investigate the observational viability of a class of models with interaction in the dark sector. We perform a Bayesian analysis using the latest data sets of type Ia supernovae, baryon acoustic oscillations, the angular acoustic scale of the cosmic microwave background, and measurements of the expansion rate. When combined with the current measurement of the local expansion rate obtained by the Hubble Space Telescope, we find that these observations provide evidence in favour of interacting models with respect to the standard cosmology.
A fraction of galaxy clusters host diffuse radio sources whose origins are investigated through multi-wavelength studies of cluster samples. We investigate the presence of diffuse radio emission in a sample of seven galaxy clusters in the largely unexplored intermediate redshift range (0.3 < z < 0.44). In search of diffuse emission, deep radio imaging of the clusters are presented from wide band (1.1-3.1 GHz), full resolution ($\sim$ 5 arcsec) observations with the Australia Telescope Compact Array (ATCA). The visibilities were also imaged at lower resolution after point source modelling and subtraction and after a taper was applied to achieve better sensitivity to low surface brightness diffuse radio emission. In case of non-detection of diffuse sources, we set upper limits for the radio power of injected diffuse radio sources in the field of our observations. Furthermore, we discuss the dynamical state of the observed clusters based on an X-ray morphological analysis with XMM-Newton. We detect a giant radio halo in PSZ2 G284.97-23.69 (z=0.39) and a possible diffuse source in the nearly relaxed cluster PSZ2 G262.73-40.92 (z=0.421). Our sample contains three highly disturbed massive clusters without clear traces of diffuse emission at the observed frequencies. We were able to inject modelled radio halos with low values of total flux density to set upper detection limits; however, with our high-frequency observations we cannot exclude the presence of RH in these systems because of the sensitivity of our observations in combination with the high z of the observed clusters.
We study the scaling relations between the baryonic content and total mass of groups of galaxies, as these systems provide a unique way to examine the role of non-gravitational processes in structure formation. Using Planck and ROSAT data, we conduct detailed comparisons of the stacked thermal Sunyaev-Zel'dovich (tSZ) effect and X-ray scaling relations of galaxy groups found in the the Galaxy And Mass Assembly (GAMA) survey and the BAHAMAS hydrodynamical simulation. We use weak gravitational lensing data from the Kilo Degree Survey (KiDS) to determine the average halo mass of the studied systems. We analyse the simulation in the same way, using realistic weak lensing, X-ray, and tSZ synthetic observations. Furthermore, to keep selection biases under control, we employ exactly the same galaxy selection and group identification procedures to the observations and simulation. Applying this careful comparison, we find that the simulations are in agreement with the observations, particularly with regards to the scaling relations of the lensing and tSZ results. This finding demonstrates that hydrodynamical simulation have reached the level of realism that is required to interpret observational survey data and study the baryon physics within dark matter haloes, where analytical modelling is challenging. Finally, using simulated data, we demonstrate that our observational processing of the X-ray and tSZ signals is free of significant biases. We find that our optical group selection procedure has, however, some room for improvement.
We jointly analyze Bolocam Sunyaev-Zeldovich (SZ) effect and Chandra X-ray data for a set of 45 clusters to derive gas density and temperature profiles without using spectroscopic information. The sample spans the mass and redshift range $3 \times 10^{14} M_{\odot} \le M_{500} \le 25 \times 10^{14} M_{\odot}$ and $0.15\le z \le 0.89$. We define cool-core (CC) and non-cool core (NCC) subsamples based on the central X-ray luminosity, and 17/45 clusters are classified as CC. In general, the profiles derived from our analysis are found to be in good agreement with previous analyses, and profile constraints beyond $r_{500}$ are obtained for 34/45 clusters. In approximately 30% of the CC clusters our analysis shows a central temperature drop with a statistical significance of $>3\sigma$; this modest detection fraction is due mainly to a combination of coarse angular resolution and modest S/N in the SZ data. Most clusters are consistent with an isothermal profile at the largest radii near $r_{500}$, although 9/45 show a significant temperature decrease with increasing radius. The sample mean density profile is in good agreement with previous studies, and shows a minimum intrinsic scatter of approximately 10% near $0.5 \times r_{500}$. The sample mean temperature profile is consistent with isothermal, and has an intrinsic scatter of approximately 50% independent of radius. This scatter is significantly higher compared to earlier X-ray-only studies, which find intrinsic scatters near 10%, likely due to a combination of unaccounted for non-idealities in the SZ noise, projection effects, and sample selection.
We study the potential of the kinematic SZ effect as a probe for cosmology, focusing on the pairwise method. The main challenge is disentangling the cosmologically interesting mean pairwise velocity from the cluster optical depth and the associated uncertainties on the baryonic physics in clusters. Furthermore, the pairwise kSZ signal might be affected by internal cluster motions or correlations between velocity and optical depth. We investigate these effects using the Magneticum cosmological hydrodynamical simulations, one of the largest simulations of this kind performed to date. We produce tSZ and kSZ maps with an area of $\simeq 1600~\mathrm{deg}^2$, and the corresponding cluster catalogues with $M_{500c} \gtrsim 3 \times 10^{13}~h^{-1}M_\odot$ and $z \lesssim 2$. From these data sets we calibrate a scaling relation between the average Compton-$y$ parameter and optical depth. We show that this relation can be used to recover an accurate estimate of the mean pairwise velocity from the kSZ effect, and that this effect can be used as an important probe of cosmology. We demonstrate that the residual systematic effects seen in our analysis are well below the remaining uncertainties on the sub-grid feedback models implemented in hydrodynamical simulations.
Neutrinos are key astronomical messengers, because they are undeflected by magnetic field and unattenuated by electromagnetic interaction. After the first detection of extraterrestrial neutrinos in the TeV-PeV region by Neutrino Telescopes we are entering a new epoch where neutrino astronomy becomes possible. In this paper I briefly review the main issues concerning cosmological neutrinos and their experimental observation.
We perform a detailed confrontation of various oscillating dark-energy parame-trizations with the latest sets of observational data. In particular, we use data from Joint Light Curve analysis (JLA) sample from Supernoave Type Ia, Baryon Acoustic Oscillations (BAO) distance measurements, Cosmic Microwave Background (CMB) observations, redshift space distortion, weak gravitational lensing, Hubble parameter measurements from cosmic chronometers, and the local Hubble constant value, and we impose constraints on four oscillating models. We find that all models are bent towards the phantom region, nevertheless in the three of them the quintessential regime is also allowed within 1$\sigma$ confidence-level. Furthermore, the deviations from $\Lambda$CDM cosmology are small, however for two of the models they could be visible at large scales, through the impact on the temperature anisotropy of the CMB spectra and on the matter power spectra.
Shocks produced by the accretion of infalling gas in the outskirt of galaxy clusters are expected in the hierarchical structure formation scenario, as found in cosmological hydrodynamical simulations. Here, we report the detection of a shock front at a large radius in the pressure profile of the galaxy cluster A2319 at a significance of $8.6\sigma$, using Planck thermal Sunyaev-Zel'dovich data. The shock is located at $(2.93 \pm 0.05) \times R_{500}$ and is not dominated by any preferential radial direction. Using a parametric model of the pressure profile, we derive a lower limit on the Mach number of the infalling gas, $\mathcal{M} > 3.25$ at 95\% confidence level. These results are consistent with expectations derived from hydrodynamical simulations. Finally, we use the shock location to constrain the accretion rate of A2319 to $\dot{M} \simeq (1.4 \pm 0.4) \times 10^{14}$ M$_\odot$ Gyr$^{-1}$, for a total mass, $M_{200} \simeq 10^{15}$ M$_\odot$.
We present a public data release of halo catalogs from a suite of 125 cosmological $N$-body simulations from the Abacus project. The simulations span 40 $w$CDM cosmologies centered on the Planck 2015 cosmology at two mass resolutions, $4\times 10^{10}\;h^{-1}M_\odot$ and $1\times 10^{10}\;h^{-1}M_\odot$, in $1.1\;h^{-1}\mathrm{Gpc}$ and $720\;h^{-1}\mathrm{Mpc}$ boxes, respectively. The boxes are phase-matched to suppress sample variance and isolate cosmology dependence. Additional volume is available via 16 boxes of fixed cosmology and varied phase; a few boxes of single-parameter excursions from Planck 2015 are also provided. Catalogs spanning $z=1.5$ to $0.1$ are available for friends-of-friends and Rockstar halo finders and include particle subsamples. All data products are available at https://lgarrison.github.io/AbacusCosmos
In preparation for deep extragalactic imaging with the James Webb Space Telescope, we explore the clustering of massive halos at $z=8$ and $10$ using a large N-body simulation. We find that halos with masses $10^9$ to $10^{11}$ $h^{-1}\;M_\odot$, which are those expected to host galaxies detectable with JWST, are highly clustered with bias factors ranging from 5 and 30 depending strongly on mass, as well as on redshift and scale. This results in correlation lengths of 5--10$h^{-1}\;{\rm Mpc}$, similar to that of today's galaxies. Our results are based on a simulation of 130 billion particles in a box of $250h^{-1}\;{\rm Mpc}$ size using our new high-accuracy ABACUS simulation code, the corrections to cosmological initial conditions of (Garrison et al. 2016, 2016MNRAS.461.4125G), and the Planck 2015 cosmology. We use variations between sub-volumes to estimate the detectability of the clustering. Because of the very strong inter-halo clustering, we find that surveys of order 25$h^{-1}\;{\rm Mpc}$ comoving transverse size may be able to detect the clustering of $z=8$--$10$ galaxies with only 500-1000 survey objects if the galaxies indeed occupy the most massive dark matter halos.
The Milky Way dark matter halo is formed from the accretion of smaller subhalos. These sub-units also harbor stars---typically old and metal-poor---that are deposited in the Galactic inner regions by disruption events. In this Letter, we show that the dark matter and metal-poor stars in the Solar neighborhood share similar kinematics due to their common origin. Using the high-resolution Eris simulation, which traces the evolution of both the dark matter and baryons in a realistic Milky-Way analog galaxy, we demonstrate that metal-poor stars are indeed effective tracers for the local, virialized dark matter velocity distribution. The local dark matter velocities can therefore be inferred from observations of the stellar halo made by the Sloan Digital Sky Survey within 4 kpc of the Sun. This empirical distribution differs from the Standard Halo Model in important ways and suggests that the bounds on the spin-independent scattering cross section may be weakened for dark matter masses below $\sim$10 GeV. Data from Gaia will allow us to further refine the expected distribution for the smooth dark matter component, and to test for the presence of local substructure.
It has been pointed out that a large non-minimal coupling $\xi$ between the Higgs and the Ricci scalar can source higher derivative operators, which may change the predictions of Higgs inflation. A variant, called critical Higgs inflation, employs the near-criticality of the top mass to introduce an inflection point in the potential and lower drastically the value of $\xi$. We here study whether critical Higgs inflation can occur even if the pre-inflationary initial conditions do not satisfy the slow-roll behaviour (retaining translation and rotation symmetries). A positive answer is found: inflation turns out to be an attractor and therefore no fine-tuning of the initial conditions is necessary. A very large initial Higgs time-derivative (as compared to the potential energy density) is compensated by a moderate increase in the initial field value. These conclusions are reached by solving the exact Higgs equation without using the slow-roll approximation. This also allows us to treat consistently the inflection point, where the standard slow-roll approximation breaks down. Here we make use of an approach that is independent of the UV completion of gravity, by taking initial conditions that always involve sub-planckian energies.
We study the incidence of group and filamentary dwarf galaxy accretion into Milky Way (MW) mass haloes using two types of hydrodynamical simulations: EAGLE, which resolves a large cosmological volume, and the AURIGA suite, which are very high resolution zoom-in simulations of individual MW-sized haloes. The present-day 11 most massive satellites are predominantly (75%) accreted in single events, 14% in pairs and 6% in triplets, with higher group multiplicities being unlikely. Group accretion becomes more common for fainter satellites, with 60% of the top 50 satellites accreted singly, 12% in pairs, and 28% in richer groups. A group similar in stellar mass to the Large Magellanic Cloud (LMC) would bring on average 15 members with stellar mass larger than $10^4{~\rm M_\odot}$. Half of the top 11 satellites are accreted along the two richest filaments. The accretion of dwarf galaxies is highly anisotropic, taking place preferentially perpendicular to the halo minor axis, and, within this plane, preferentially along the halo major axis. The satellite entry points tend to be aligned with the present-day central galaxy disc and satellite plane, but to a lesser extent than with the halo shape. Dwarfs accreted in groups or along the richest filament have entry points that show an even larger degree of alignment with the host halo than the full satellite population. We also find that having most satellites accreted as a single group or along a single filament is unlikely to explain the MW disc of satellites.
If dark energy (DE) is a dynamical field rather than a cosmological constant, an interaction between DE and the neutrino sector could exist, modifying the neutrino oscillation phenomenology and causing CP and apparent Lorentz violating effects. The terms in the Hamiltonian for flavor propagation induced by the DE-neutrino coupling do not depend on the neutrino energy, while the ordinary components decrease as $\Delta m^2/E_{\nu}$. Therefore, the DE-induced effects are absent at lower neutrino energies, but become significant at higher energies, allowing to be searched for by neutrino observatories. We explore the impact of the DE-neutrino coupling on the oscillation probability and the flavor transition in the three-flavor framework, and investigate the CP-violating and apparent Lorentz violating effects. We find that DE-induced effects become observable for $E_{\nu}m_{\text{eff}} \sim 10^{-20}~ \text{GeV}^2$, where $m_{\rm eff}$ is the effective mass parameter in the DE-induced oscillation probability, and CP is violated over a wide energy range. We also show that current and future experiments have the sensitivity to detect anomalous effects induced by a DE-neutrino coupling and probe the new mixing parameters. The DE-induced effects on neutrino oscillation can be distinguished from other new physics possibilities with similar effects, through the detection of the directional dependence of the interaction, which is specific to this interaction with DE. However, current experiments will not yet be able to measure the small changes of $\sim 0.03\%$ in the flavor composition due to this directional effect.
Within a galaxy the stellar mass-to-light ratio $\Upsilon_*$ is not constant. We show that ignoring $\Upsilon_*$ gradients can have a more dramatic effect on dynamical ($M_*^{\rm dyn}$) compared to stellar population ($M_*^{\rm SP}$) based estimates of early-type galaxy stellar masses, because $M_*^{\rm dyn}$ is usually calibrated using the velocity dispersion measured in the central regions. If $\Upsilon_*$ is greater there, then ignoring the gradient will lead to an overestimate of $M_*^{\rm dyn}$. Spatially resolved kinematics of nearby early-type galaxies suggests that these gradients are driven by gradients in the initial mass function (IMF). Accounting for recent estimates of the IMF-driven $\Upsilon_*$ gradient reduces $M_*^{\rm dyn}$ substantially ($\sim$ a factor of two), and may be accompanied by a (smaller) change in $M_*^{\rm SP}$. Our results suggest that $M_*^{\rm dyn}$ estimates in the literature should be revised downwards, rather than revising $M_*^{\rm SP}$ estimates upwards. This has three consequences. First, if gradients in $\Upsilon_*$ are present, then $M_*^{\rm dyn}$ cannot be estimated independently of stellar population synthesis models. Second, accounting for $\Upsilon_*$ gradients changes the slope of the stellar mass function $\phi(M_*^{\rm dyn})$, and reduces the associated stellar mass density, especially at high masses. Third, if gradients are stronger in more massive galaxies, then accounting for this reduces the slope of the correlation between the ratio of the dynamical and stellar population mass estimates of a galaxy with its velocity dispersion. These conclusions potentially impact estimates of the need for feedback and adiabatic contraction, so our results highlight the importance of measurements of $\Upsilon_*$ gradients in larger samples.
Links to: arXiv, form interface, find, astro-ph, recent, 1712, contact, help (Access key information)
Hubble Space Telescope observations from December 11 2015 detected the expected fifth counter image of SN Refsdal at $z = 1.49$. In this letter, we compare the time delay predictions from numerous models with the measured value derived by Kelly et al. (2016) from very early data in the light curve of the SN Refsdal, and find a best value for $H_0 = 66^{+12}_{-12}~\mathrm{km~s^{-1}~Mpc^{-1}}$ (68\% CL), in excellent agreement with predictions from CMB and recent weak lensing data + BAO + BBN (from the DES Collaboration). This is the first constraint on $H_0$ derived from time delays between multiple lensed SN images, and the first with a galaxy cluster lens, so subject to systematic effects different from other time-delay $H_0$ estimates. Additional time delay measurements from new multiply-imaged SNe will allow derivation of competitive constraints on $H_0$.
The massive stars that ionised the Universe have short lifetimes and can only be studied near the time of formation, but any low mass star that formed contemporaneously might be observable in the local Universe today. We study the abundance pattern and spatial distribution of these `siblings of reionizers' (SoRs) in the EAGLE cosmological hydrodynamical simulation. SoRs are enriched to super-solar levels in $\alpha$-elements compared to iron. In particular we find that a large fraction of SoRs are carbon-enhanced meal poor stars that are not overabundant in neutron-capture elements - CEMP-no stars - and vice versa. Most SoRs today are hosted by halos with mass $M_h\gtrapprox 10^{12}$ M$_\odot$, and 50 percent of them are in the halo of their central galaxy (distance $>10$ kpc), mainly because they were accreted onto their host rather than formed in-situ. To a good approximation, the siblings of reionizers are CEMP-no stars that reside today in the stellar halos of massive galaxies, with nearly half of them contributing to the intracluster light in groups and clusters. Wherever possible we have compared the results of the simulation with observations.
We review the Equation of State (EoS) approach to dark sector perturbations and apply it to $f(\mathcal{R})$ gravity models of dark energy. We show that the EoS approach is numerically stable and use it to set observational constraints on designer models. Within the EoS approach we build an analytical understanding of the dynamics of cosmological perturbations for the designer class of $f(\mathcal{R})$ gravity models, characterised by the parameter $B_0$ and the background equation of state of dark energy $w$. When we use the Planck Cosmic Microwave Background (CMB) temperature anisotropy, polarisation and lensing data as well as the Baryonic Acoustic Oscillation (BAO) data from SDSS and WiggleZ, we find $B_0<0.006$ (95\%CL) for the designer models with $w=-1$. Furthermore, we find $B_0<0.0045$ and $|w+1|<0.002$ (95\%CL) for the designer models with $w\neq -1$. Previous analyses found similar results for designer and Hu-Sawicki $f(\mathcal{R})$ gravity models using the Effective Field Theory (EFT) approach (Raveri et al. 2014, Hu et al. 2016), therefore this hints for the fact that generic $f(\mathcal{R})$ models with $w\neq-1$ can be tightly constrained by current cosmological data, complementary to solar system tests (Brax et al.2008, Faulkner et al. 2007).
Cosmological large scale structure $N$-body simulations are computation-light, memory-heavy problems in supercomputing. The considerable amount of memory is usually dominated by an inefficient way of storing more than sufficient phase space information of particles. We present a new parallel, information-optimized, particle-mesh based $N$-body code CUBE, in which information- and memory-efficiency is increased by nearly an order of magnitude. This is accomplished by storing particle's relative phase space coordinates instead of global values, and in the format of fixed point as light as 1-byte. The remaining information is given by complimentary density and velocity fields (negligible in memory space) and proper ordering of particles (no extra memory). Our numerical experiments show that this information-optimized $N$-body algorithm provides accurate results within the error of the particle-mesh algorithm. This significant lowering of memory-to-computation ratio breaks the bottleneck of scaling up and speeding up large cosmological $N$-body simulations on multi-core and heterogenous computing systems.
We present angular diameter distance measurements obtained by locating the BAO scale in the distribution of galaxies selected from the first year of Dark Energy Survey data. We consider a sample of over 1.3 million galaxies distributed over a footprint of 1318 deg$^2$ with $0.6 < z_{\rm photo} < 1$ and a typical redshift uncertainty of $0.03(1+z)$. This sample was selected, as fully described in a companion paper, using a color/magnitude selection that optimizes trade-offs between number density and redshift uncertainty. We investigate the BAO signal in the projected clustering using three conventions, the angular separation, the co-moving transverse separation, and spherical harmonics. Further, we compare results obtained from template based and machine learning photometric redshift determinations. We use 1800 simulations that approximate our sample in order to produce covariance matrices and allow us to validate our distance scale measurement methodology. We measure the angular diameter distance, $D_A$, at the effective redshift of our sample divided by the true physical scale of the BAO feature, $r_{\rm d}$. We obtain close to a 4 per cent distance measurement of $D_A(z_{\rm eff}=0.81)/r_{\rm d} = 10.75\pm 0.43 $. These results are consistent with the flat $\Lambda$CDM concordance cosmological model supported by numerous other recent experimental results.
We define and characterise a sample of 1.3 million galaxies extracted from the first year of Dark Energy Survey data, optimised to measure Baryon Acoustic Oscillations in the presence of significant redshift uncertainties. The sample is dominated by luminous red galaxies located at redshifts $z \gtrsim 0.6$. We define the exact selection using color and magnitude cuts that balance the need of high number densities and small photometric redshift uncertainties, using the corresponding forecasted BAO distance error as a figure-of-merit in the process. The typical photo-$z$ uncertainty varies from $2.3\%$ to $3.6\%$ (in units of 1+$z$) from $z=0.6$ to $1$, with number densities from $200$ to $130$ galaxies per deg$^2$ in tomographic bins of width $\Delta z = 0.1$. Next we summarise the validation of the photometric redshift estimation. We characterise and mitigate observational systematics including stellar contamination, and show that the clustering on large scales is robust in front of those contaminants. We show that the clustering signal in the auto-correlations and cross-correlations is generally consistent with theoretical models, which serves as an additional test of the redshift distributions.
Mock catalogues are a crucial tool in the analysis of galaxy surveys data, both for the accurate computation of covariance matrices, and for the optimisation of analysis methodology and validation of data sets. In this paper, we present a set of 1800 galaxy mock catalogues designed to match the Dark Energy Survey Year-1 BAO sample (Crocce et al. 2017) in abundance, observational volume, redshift distribution and uncertainty, and redshift dependent clustering. The simulated samples were built upon HALOGEN (Avila et al. 2015) halo catalogues, based on a $2LPT$ density field with an exponential bias. For each of them, a lightcone is constructed by the superposition of snapshots in the redshift range $0.45<z<1.4$. Uncertainties introduced by so-called photometric redshifts estimators were modelled with a \textit{double-skewed-Gaussian} curve fitted to the data. We also introduce a hybrid HOD-HAM model with two free parameters that are adjusted to achieve a galaxy bias evolution $b(z_{\rm ph})$ that matches the data at the 1-$\sigma$ level in the range $0.6<z_{\rm ph}<1.0$. We further analyse the galaxy mock catalogues and compare their clustering to the data using the angular correlation function $ w(\theta)$, the comoving transverse separation clustering $\xi_{\mu<0.8}(s_{\perp})$ and the angular power spectrum $C_\ell$.
In this thesis, we have studied the generation and evolution of the magnetic field in the early Universe. We investigated the generation of magnetic fields in the presence of chiral symmetry and gravitational anomaly. We have used modified kinetic theory by Berry curvature and derived expressions for magnetic and vortical conductivities. We have also shown that there are two length scales: kinetic and viscous, which is one of the important results. In the latter part of the thesis, we have shown the effect of second-order viscous effects on the modes of the magnetic fields.
We discuss the constraints on the Hubble constant to be expected from standard siren sources in ground-based gravitational wave detectors. We consider binary neutron star and binary black hole sources, and focus on the role of golden sirens (the loudest and best constrained sources) to constrain cosmological parameters. We consider two approaches: the counterpart case, where electromagnetic observations provide an independent measurement of the redshift to the sources, and the statistical case, incorporating an analysis over all potential host galaxies within the localization volumes. Our analysis includes realistic measurement uncertainties and selection biases. Although the specific results depend on the configuration and sensitivity of the detector networks, we find that the statistical method would constrain $H_0$ to $4\%$ with $\sim100$ detections of binary neutron star mergers, while an equivalent statistical measurement can be accomplished with $\sim30$ golden events (all fractional uncertainties are quoted as half of the width of the symmetric 68% credible interval divided by the median of the $H_0$ posterior). Alternatively, with $\sim10/60/200$ binary neutron star standard sirens with electromagnetic counterparts, $H_0$ would be constrained to $4/2/1\%$. Given current rate uncertainties, the 5% measurement may happen within the 9 month O3 LIGO/Virgo run (starting Fall 2018), or if the rate is low this may not happen until one full year of the LIGO/Virgo network at design sensitivity (starting $\sim$2021+). Similarly, the 1% measurement may happen within two years of running at design sensitivity ($\sim$2024+), or may not happen until 3+ years of operation of a full five-detector network. Although the rates, and thus precise timetable, remain uncertain, precision standard siren cosmology can be expected in the foreseeable future.
The recent simultaneous detection of gravitational waves and a gamma ray burst from a neutron star merger significantly shrank the space of viable scalar-tensor theories by demanding that the speed of gravity is equal to that of light. The survived theories belong to the class of degenerate higher order scalar-tensor theories. We study whether these theories are suitable as dark energy candidates. We find scaling solutions in the matter dominated universe that lead to de Sitter solutions at late times without the cosmological constant, realising self-acceleration. We evaluate quasi-static perturbations around self-accelerating solutions and show that the stringent constraints coming from astrophysical objects and gravitational waves can be satisfied, leaving interesting possibilities to test these theories by cosmological observations.
It has been claimed in Ref.[arXiv:1712.02240] that massive primordial black holes (PBH) cannot constitute all of the dark matter (DM), because their gravitational-lensing imprint on the Hubble diagram of type Ia supernovae (SN) would be incompatible with present observations. In this paper, we critically review those constraints and find several caveats on the analysis. First of all, the constraints on the fraction $\alpha$ of PBH in matter seem to be driven by a very restrictive choice of priors on the cosmological parameters. In particular, the degeneracy between $\Omega_{\rm M}$ and $\alpha$ is ignored and thus, by fixing $\Omega_{\rm M}$, transferred the constraining power of SN magnitudes to $\alpha$. Furthermore, by considering more realistic physical sizes for the type-Ia supernovae, we find an effect on the SN lensing magnification distribution that leads to significantly looser constraints. Moreover, considering a wide mass spectrum of PBH, such as a lognormal distribution, further softens the constraints from SN lensing. Finally, we find that the fraction of PBH that could constitute DM today is bounded by $f_{\rm PBH} < 1.09\ (1.38)$, for JLA (Union 2.1) catalogs, and thus it is perfectly compatible with an all-PBH dark matter scenario in the LIGO band.
In a recent study, using the distribution of galaxies in the north galactic pole of SDSS DR7 region enclosed in a 500$\mpch$ box, we carried out our ELUCID simulation (Wang et al. 2016, ELUCID III). Here we {\it light} the dark matter halos and subhalos in the reconstructed region in the simulation with galaxies in the SDSS observations using a novel {\it neighborhood} abundance matching method. Before we make use of thus established galaxy-subhalo connections in the ELUCID simulation to evaluate galaxy formation models, we set out to explore the reliability of such a link. For this purpose, we focus on the following a few aspects of galaxies: (1) the central-subhalo luminosity and mass relations; (2) the satellite fraction of galaxies; (3) the conditional luminosity function (CLF) and conditional stellar mass function (CSMF) of galaxies; and (4) the cross correlation functions between galaxies and the dark matter particles, most of which are measured separately for all, red and blue galaxy populations. We find that our neighborhood abundance matching method accurately reproduces the central-subhalo relations, satellite fraction, the CLFs and CSMFs and the biases of galaxies. These features ensure that thus established galaxy-subhalo connections will be very useful in constraining galaxy formation processes. And we provide some suggestions on the three levels of using the galaxy-subhalo pairs for galaxy formation constraints. The galaxy-subhalo links and the subhalo merger trees in the SDSS DR7 region extracted from our ELUCID simulation are available upon request.
While low mass, star forming galaxies are often considered as the primary driver of reionization, their actual contribution to the cosmic ultraviolet background is still uncertain, mostly because the escape fraction of ionizing photons is only poorly constrained. Theoretical studies have shown that efficient supernova feedback is a necessary condition to create paths through which ionizing radiation can escape into the intergalactic medium. We investigate the possibility that accreting supermassive black holes in early dwarf galaxies may provide additional feedback and enhance the leakage of ionizing radiation. We use a series of high resolution cosmological radiation hydrodynamics simulations where we isolate the different sources of feedback. We find that supernova feedback prevents the growth of the black hole, thus quenching its associated feedback. Even in cases where the black hole can grow, the structure of the interstellar medium is strongly dominated by supernova feedback. We conclude that, in the dwarf galaxy regime, supermassive black holes do not appear to play a significant role in enhancing the escape fraction and in contributing to the early UV background.
We present nbodykit, an open-source, massively parallel Python toolkit for analyzing large-scale structure (LSS) data. Using Python bindings of the Message Passing Interface (MPI), we provide parallel implementations of many commonly used algorithms in LSS. nbodykit is both an interactive and scalable piece of scientific software, performing well in a supercomputing environment while still taking advantage of the interactive tools provided by the Python ecosystem. Existing functionality includes estimators of the power spectrum, 2 and 3-point correlation functions, a Friends-of-Friends grouping algorithm, mock catalog creation via the halo occupation distribution technique, and approximate N-body simulations via the FastPM scheme. The package also provides a set of distributed data containers, insulated from the algorithms themselves, that enable nbodykit to provide a unified treatment of both simulation and observational data sets. nbodykit can be easily deployed in a high performance computing environment, overcoming some of the traditional difficulties of using Python on supercomputers. We provide performance benchmarks illustrating the scalability of the software. The modular, component-based approach of nbodykit allows researchers to easily build complex applications using its tools. The package is extensively documented at this http URL, which also includes an interactive set of example recipes for new users to explore. As open-source software, we hope nbodykit provides a common framework for the community to use and develop in confronting the analysis challenges of future LSS surveys.
Radio-loud active galactic nuclei (AGNs), hosting powerful relativistic jet outflows, provide an excellent laboratory for studying jet physics. Very long baseline interferometry (VLBI) enables high-resolution imaging on milli-arcsecond (mas) and sub-mas scales, making it a powerful tool to explore the inner jet structure, shedding light on the formation, acceleration and collimation of AGN jets. In this paper, we present Very Long Baseline Array (VLBA) observations of ten radio-loud AGNs at 43 and 86~GHz, which were selected from the {\it Planck} catalogue of compact sources and are among the brightest in published VLBI images at and below 15 GHz. The image noise levels in our observations are typically 0.3 mJy beam$^{-1}$ and 1.5 mJy beam$^{-1}$ at 43 and 86 GHz, respectively. Compared with the VLBI data observed at lower frequencies from the literature, our observations with higher resolution (the highest resolution up to 0.07 mas at 86 GHz and 0.18 mas at 43 GHz) and at higher frequencies detected new jet components at sub-parsec scales, offering valuable data for studies of the physical properties of innermost jets. These include compactness factor of the radio structure (the ratio of core flux density to total flux density), and core brightness temperature ($T_{\rm b}$). In all these sources, the compact core accounts for a significant fraction ($> 60\%$) of the total flux density. Their correlated flux density at the longest baselines is higher than 0.16 Jy. The compactness of these sources make them good phase calibrators of mm-wavelength ground-based and space VLBI.
The neutrino asymmetry in the early universe plasma, $n_\nu - n_{\bar \nu}$, is calculated both before and after the electroweak phase transition (EWPT). In the Standard Model before EWPT, the leptogenesis is well known to be driven by the abelian anomaly in a massless hypercharge field. The generation of the neutrino asymmetry in the Higgs phase after EWPT, in its turn, has not been considered previously because of the absence of any quantum anomaly in an external electromagnetic field for such electroneutral particles as neutrino, unlike the Adler-Bell-Jackiw anomaly for charged left and right polarized massless electrons in the same electromagnetic field. Using the neutrino Boltzmann equation, modified by the Berry curvature term in the momentum space, we establish the violation of the macroscopic neutrino current in plasma after EWPT and exactly reproduce the nonconservation of the lepton current in the symmetric phase before EWPT arising in quantum field theory due to the nonzero lepton hypercharge and corresponding triangle anomaly in an external hypercharge field. In the last case, the non-conservation of the lepton current is derived through the kinetic approach without a computation of corresponding Feynman diagrams. Then the new kinetic equation is applied for the calculation of the neutrino asymmetry accounting for the Berry curvature and the electroweak interaction with background fermions in the Higgs phase. Such an interaction generates a neutrino asymmetry through the electroweak coupling of neutrino currents with electromagnetic fields in plasma, which is $\sim G_\mathrm{F}^2$. It turns out that this effect is especially efficient for maximally helical magnetic fields.
The present day spectrum of the extragalactic background light (EBL) in UV, optical and IR wavelengths is the integral result of multiple astrophysical processes going on throughout the evolution of the Universe. The relevant processes include star formation, stellar evolution, light absorption and emission by the cosmic dust. The properties of these processes are known with uncertainties which contribute to the EBL spectrum precision. In the present paper we develop a numerical model of the EBL spectrum while maintaining the explicit dependence on the astrophysical parameters involved. We constructed a Markov Chain in the parameter space by using the likelihood function built with the up-to-date upper and lower bounds on the EBL intensity. The posterior distributions built with the Markov Chain Monte Carlo method are used to determine an allowed range of the individual parameters of the model. Consequently, the star formation rate multiplication factor is constrained in the range $0.94 < C_{sfr} < 1.41$ at $68\%$ C.L. The method also results in the bounds on the lifetime, radius, dust particle density and opacity of the molecular clouds that have large ambiguity otherwise. It is shown that there is a reasonable agreement between the model and the intensity bounds while the astrophysical parameters of the best fit model are close to their estimates from literature.
Links to: arXiv, form interface, find, astro-ph, recent, 1712, contact, help (Access key information)
The observed surface densities of dark matter halos are known to follow a simple scaling law, ranging from dwarf galaxies to galaxy clusters, with a weak dependence on their virial mass. Here we point out that this can not only be used to provide a method to determine the standard relation between halo mass and concentration, but also to use large samples of objects in order to place constraints on dark matter self-interactions that can be more robust than constraints derived from individual objects. We demonstrate our method by considering a sample of about 50 objects distributed across the whole halo mass range. Modelling the effect of self-interactions in a way similar to what has been previously used in the literature, we derive a constraint on the self-interaction cross section per unit dark matter mass of about $\sigma/m_\chi\lesssim 0.2$\,cm$^2$/g. Formally, this is already somewhat smaller than the range $0.5-5$\,cm$^2$/g that has been invoked to explain potential inconsistencies between small-scale observations and expectations in the standard collisionless cold dark matter paradigm. We expect that these constraints can be significantly improved in the future, and made more robust, by i) an improved theoretical modelling of the effect of self-interactions, ii) taking into account a larger sample of objects and iii) by reducing the currently still relatively large uncertainties that we conservatively assign to the surface densities of individual objects. The latter can be achieved in particular by using kinematic observations to directly constrain the average halo density inside a given radius, rather than fitting the data to a pre-selected profile and then reconstruct the profile.
The combination of X-ray and SZ observations can potentially improve the cluster detection efficiency when compared to using only one of these probes, since both probe the same medium: the hot ionized gas of the intra-cluster medium. We present a method based on matched multifrequency filters (MMF) for detecting galaxy clusters from SZ and X-ray surveys. This method builds on a previously proposed joint X-ray-SZ extraction method (Tarr\'io et al. 2016) and allows to blindly detect clusters, that is finding new clusters without knowing their position, size or redshift, by searching on SZ and X-ray maps simultaneously. The proposed method is tested using data from the ROSAT all-sky survey and from the Planck survey. The evaluation is done by comparison with existing cluster catalogues in the area of the sky covered by the deep SPT survey. Thanks to the addition of the X-ray information, the joint detection method is able to achieve simultaneously better purity, better detection efficiency and better position accuracy than its predecessor Planck MMF, which is based on SZ maps only. For a purity of 85%, the X-ray-SZ method detects 141 confirmed clusters in the SPT region, whereas to detect the same number of confirmed clusters with Planck MMF, we would need to decrease its purity to 70%. We provide a catalogue of 225 sources selected by the proposed method in the SPT footprint, with masses ranging between 0.7 and 14.5 $\cdot 10^{14}$ Msun and redshifts between 0.01 and 1.2.
The nature of dark matter is one of the most pressing questions in particle physics. Yet all our present knowledge of the dark sector to date comes from its gravitational interactions with astrophysical systems. Moreover, astronomical results still have immense potential to constrain the particle properties of dark matter. We introduce a simple 2D parameter space which classifies models in terms of a particle physics interaction strength and a characteristic astrophysical scale on which new physics appears, in order to facilitate communication between the fields of particle physics and astronomy. We survey the known astrophysical anomalies that are suggestive of non-trivial dark matter particle physics, and present a theoretical and observational program for future astrophysical measurements that will shed light on the nature of dark matter.
We investigate how the nature of the galaxies that reionized the Universe affects the duration of reionization. We contrast two models: one in which galaxies on the faint side of the luminosity function dominate the ionizing emissivity, and a second in which the galaxies on the bright side of the luminosity function dominate. The faint-end of the luminosity function evolves slowly, therefore the transition from mostly neutral to mostly ionized state takes a much longer time in the first model compared to the second. Existing observational constraints on the duration of this transition are relatively weak, but taken at face value prefer the model in which galaxies on the bright side play a major role. Measurements of the kinetic Sunyaev Zeldovich effect in the cosmic microwave background from the epoch of reionization also point in the same direction.
We perform simulations in a simple model that aims to mimic the hydrodynamic evolution of a relativistic fluid during a cosmological first-order phase transitions. The observable we are concerned with is hereby the spectrum of gravitational radiation produced by colliding fluid shells. We present simple parameterizations of our results as functions of the wall velocity, the duration of the phase transition and the latent heat. We also improve on previous results in the envelope approximation and compare with hydrodynamic simulations.
We use hydrodynamical/N-body simulations to interpret the newly discovered Bullet-cluster-like merging cluster, ZwCl 0008.8+5215 (ZwCl 0008 hereafter), where a dramatic collision is apparent from multi-wavelength observations. We have been able to find a self-consistent solution for the radio, X-ray, and lensing phenomena by projecting an off-axis, binary cluster encounter viewed just after first core passage. A pair radio relics traces well the leading and trailing shock fronts that our simulation predict, providing constraints on the collision parameters. We can also account for the observed distinctive comet-like X-ray morphology and the positions of the X-ray peaks relative to the two lensing mass centroids and the two shock front locations. Relative to the Bullet cluster, the total mass is about 70% lower, ($1.2\pm0.1) \times 10^{15}$ Msun, with a correspondingly lower infall velocity, $1800\pm300$ km/s, and an impact parameter of $400\pm100$ kpc. As a result, the gas component of the infalling cluster is not trailing significantly behind the associated dark matter as in the case of the Bullet cluster. The degree of agreement we find between all the observables provides strong evidence that dark matter is effectively collisionless on large scales calling into question other claims and theories that advocate modified gravity.
One of the major goals of future Cosmic Microwave Background measurements is the accurate determination of the effective number of neutrinos $N_{\rm eff}$. Reaching an experimental sensitivity of $\Delta N_{\rm eff} = 0.013$ could indeed falsify the presence of any non-standard relativistic particles at $95 \%$ c.l.. In this paper, we test how this future constraint can be affected by the removal of two common assumptions: a negligible running of the inflationary spectral index $n_{\rm run}$ and a precise determination of the neutron lifetime $\tau_n$. We first show that the constraints on $N_{\rm eff}$ could be significantly biased by the unaccounted presence of a running of the spectral index. Considering the Stage-IV experiment, a negative running of ${\rm d}n/{\rm d}\ln k= - 0.002$ could mimic a positive variation of $\Delta N_{\rm eff} = 0.03$. Moreover, given the current discrepancies between experimental measurements of the neutron lifetime $\tau_n$, we show that the assumption of a conservative error of $\Delta\tau_n \sim 10$s could bring to a systematic error of $\Delta N_{\rm eff} = 0.02$. Complementary cosmological constraints on the running of the spectral index and a solution to the neutron lifetime discrepancy are therefore needed for an accurate and reliable future CMB bound of $N_{\rm eff}$ at percent level.
Cosmological observations are a powerful probe of neutrino properties, and in particular of their mass. In this review, we first discuss the role of neutrinos in shaping the cosmological evolution at both the background and perturbation level, and describe their effects on cosmological observables such as the cosmic microwave background and the distribution of matter at large scale. We then present the state of the art concerning the constraints on neutrino masses from those observables, and also review the prospects for future experiments. We also briefly discuss the prospects for determining the neutrino hierarchy from cosmology, the complementarity with laboratory experiments, and the constraints on neutrino properties beyond their mass.
The effect of baryons on large-scale structure, if not accounted for, can significantly bias constraints on dark energy. As the detailed physics of the baryons is not yet well understood, correcting for baryon effects introduces additional parameters which must be marginalised over, increasing the uncertainties on the inferred cosmological parameters. As forthcoming Stage IV weak lensing surveys are aiming for percent-level precision on the dark energy equation of state, the extent of this problem must be thoroughly examined. In this paper, we use a halo model with analytic modifications which capture the impact of adiabatic contraction of baryons and feedback on the matter power spectrum, and generalize the Navarro-Frenk-White profile to account for a possible inner halo core. A Fisher analysis predicts a 40% degradation in the $w_0$-$w_a$ Figure of Merit for a Stage IV survey, and degradations of up to 80% in the Figure of Merit for other cosmological parameters. We find that a Euclid-like survey can constrain an inner core to a few $\mathrm{kpc}$, while if the inner core is fixed, we can constrain adiabatic concentration and feedback parameters to a few percent. We explore the scales on which baryons and dark energy contribute most to the Fisher information, finding that simply probing deeper into the non-linear regime does little to reduce the degradation. We show that including external baryon information improves forecasts, although strong priors are required to limit degradation to 1%. Adding Planck cosmic microwave background priors improves the Figure of Merit by a factor of 2.7 and almost completely recovers the individual marginalized errors on $w_0$ and $w_a$. We also quantify the calibration of baryon modelling required to reduce biases of dark energy forecasts to acceptable levels for forthcoming lensing surveys.
The decay rate of the Bunch-Davies state of a massive scalar field in the expanding flat spatial sections of de Sitter space is determined by an analysis of the particle pair creation process in real time. The Feynman definition of particle and antiparticle Fourier mode solutions of the scalar wave equation, and their adiabatic phase analytically continued to the complexified time domain, show conclusively that the Bunch-Davies state is not the vacuum state at late times. The closely analogous creation of charged particle pairs in a uniform electric field is reviewed and Schwinger's result for the vacuum decay rate is recovered by the real time analysis. The vacuum decay rate in each case is also calculated by switching the background field on adiabatically, allowing it to act for a very long time, and then adiabatically switching it off again. In both the uniform electric field and de Sitter cases the particles created while the field is switched on are verified to be real, in the sense that they persist in the final asymptotic flat zero-field region. In the de Sitter case there is an interesting residual dependence of the rate on how the de Sitter phase is ended, indicating a greater sensitivity to spatial boundary conditions. The electric current of the created particles in the E-field case and their energy density and pressure in the de Sitter case are also computed, and the magnitude of their backreaction effects on the background field estimated. Possible consequences of the Hubble scale instability of the de Sitter vacuum for cosmology, vacuum dark energy, and the cosmological `constant' problem are discussed.
We show that polar materials are ideal targets for direct detection of sub-GeV dark matter due to the presence of gapped optical phonons as well as acoustic phonons with high sound speed. We take the example of Gallium Arsenide (GaAs), which has the properties needed for experimental realization, and where many results can be estimated analytically. We find GaAs has excellent reach to dark photon absorption, can completely cover the freeze-in benchmark for scattering via an ultralight dark photon, and is competitive with other proposals to detect sub-MeV dark matter scattering off nuclei.
Following up on our recent study, we consider the regime of graviton masses and gravitational wave propagation distances at which decoherence of the wave packets plays a major role for phenomenology. This regime is of particular interest, as it can lead to very striking phenomena of echo events in the gravitational waves coming from coalescence events. The power of the experimental search in this case lies in the fact that it becomes sensitive to a large range of graviton masses, while not relying on a specific production mechanism. We are thus able to place new relevant limits on the parameter space of the graviton mixing angle.
The representation theory of de Sitter space allows for a category of partially massless particles which have no flat space analog, but could have existed during inflation. We study the couplings of these exotic particles to inflationary perturbations and determine the resulting signatures in cosmological correlators. When inflationary perturbations interact through the exchange of these fields, their correlation functions inherit scalings that cannot be mimicked by extra massive fields. We discuss in detail the squeezed limit of the tensor-scalar-scalar bispectrum, and show that certain partially massless fields can violate the tensor consistency relation of single-field inflation. We also consider the collapsed limit of the scalar trispectrum, and find that the exchange of partially massless fields enhances its magnitude, while giving no contribution to the scalar bispectrum. These characteristic signatures provide clean detection channels for partially massless fields during inflation.
If there exist higher-spin particles during inflation which are light compared to the Hubble rate, they may leave distinct statistical anisotropic imprints on the correlators involving scalar and graviton fluctuations. We characterise such signatures using the dS/CFT$_3$ correspondence and the operator product expansion techniques. In particular, we obtain generic results for the case of partially massless higher-spin states.
Up to now searches for Dark Matter (DM) detection have not been successful, either because our paradigm in how DM signals should look like are wrong or the detector sensitivity is still too low in spite of the large progress made in recent years. We discuss both possibilities starting with what we know about DM from cosmology and why Supersymmetry provides such an interesting paradigm for cosmology and particle physics in order to appreciate what it means to give up this paradigm. In addition, we compare the predicted cross sections for direct and indirect DM detection with observations with emphasis on the latest developments. Especially, we discuss the possible origins of the two hotly debated candidates for a DM annihilation signal, namely the positron excess and the Fermi GeV excess, which are unfortunately incompatible with each other and more mundane astrophysical explanations exist.
Even if the fundamental action of gravity is local, the corresponding quantum effective action, that includes the effect of quantum fluctuations, is a nonlocal object. These nonlocalities are well understood in the ultraviolet regime but much less in the infrared, where they could in principle give rise to important cosmological effects. Here we systematize and extend previous work of our group, in which it is assumed that a mass scale $\Lambda$ is dynamically generated in the infrared, giving rise to nonlocal terms in the quantum effective action of gravity. We give a detailed discussion of conceptual aspects related to nonlocal gravity and of the cosmological consequences of these models. The requirement of providing a viable cosmological evolution severely restricts the form of the nonlocal terms, and selects a model (the so-called RR model) that corresponds to a dynamical mass generation for the conformal mode. For such a model: (1) there is a FRW background evolution, where the nonlocal term acts as an effective dark energy with a phantom equation of state, providing accelerated expansion without a cosmological constant. (2) Cosmological perturbations are well behaved. (3) Implementing the model in a Boltzmann code and comparing with observations we find that the RR model fits the CMB, BAO, SNe, structure formation data and local $H_0$ measurements at a level statistically equivalent to $\Lambda$CDM. (4) Bayesian parameter estimation shows that the value of $H_0$ obtained in the RR model is higher than in $\Lambda$CDM, reducing to $2.0\sigma$ the tension with the value from local measurements. (5) The RR model provides a prediction for the sum of neutrino masses that falls within the limits set by oscillation and terrestrial experiments. (6) Gravitational waves propagate at the speed of light, complying with the limit from GW170817/GRB 170817A.
We present GAMER-2, a GPU-accelerated adaptive mesh refinement (AMR) code for astrophysics. It provides a rich set of features, including adaptive time-stepping, several hydrodynamic schemes, magnetohydrodynamics, self-gravity, particles, star formation, chemistry and radiative processes with GRACKLE, data analysis with yt, and memory pool for efficient object allocation. GAMER-2 is fully bitwise reproducible. For the performance optimization, it adopts hybrid OpenMP/MPI/GPU parallelization and utilizes the simultaneity between CPU computation, GPU computation, and CPU-GPU communication. Load balancing is achieved using a Hilbert space-filling curve on a level-by-level basis without the need to duplicate the entire AMR hierarchy on each MPI process. To provide convincing demonstrations of the accuracy and performance of GAMER-2, we directly compare with Enzo on isolated disk galaxy simulations and with FLASH on galaxy cluster merger simulations. We show that the physical results obtained by different codes are in very good agreement, and GAMER-2 outperforms Enzo and FLASH by nearly one and two orders of magnitude, respectively, on the Blue Waters supercomputers using $1-256$ nodes. More importantly, GAMER-2 exhibits similar or even better parallel scalability compared to the other two codes. We also demonstrate good weak and strong scaling using up to 4096 GPUs and 65,536 CPU cores, and achieve a uniform resolution as high as $10{,}240^3$ cells. Furthermore, GAMER-2 can be adopted as an AMR+GPUs framework and has been extensively used for the wave dark matter ($\psi$DM) simulations. GAMER-2 is open source (available at https://github.com/gamer-project/gamer) and new contributions are welcome.
Scalar fields provide attractive modifications of pre-BBN cosmology, which have interesting implications for dark matter abundances. We analyse the effect of these modifications on the number of e-foldings between horizon exit of CMB modes and the end of inflation $(N_k)$, and examine the consequences for inflationary predictions of various models. We find significant effects in the predictions in the $(n_s, r)$ plane. For a large part of the parameter space, the shift in $N_k$ is positive; this ameliorates the tension of $m^{2} \varphi^{2}$ and natural inflation with data.
Links to: arXiv, form interface, find, astro-ph, recent, 1712, contact, help (Access key information)
We present the first cosmological constraint on dark matter scattering with protons in the early Universe for the entire range of dark matter masses between 1 keV and 1 TeV. This constraint is derived from the Planck measurements of the cosmic microwave background (CMB) temperature and polarization anisotropy, and the CMB lensing anisotropy. It improves upon previous CMB constraints by many orders of magnitude, where limits are available, and closes the gap in coverage for low-mass dark matter candidates. We focus on two canonical interaction scenarios: spin-independent and spin-dependent scattering with no velocity dependence. Our results exclude (with 95% confidence) spin-independent interactions with cross sections greater than $1.3 \times 10^{-27}$ cm$^2$ for 1 keV, $7.3 \times 10^{-27}$ cm$^2$ for 1 MeV, $5.1 \times 10^{-26}$ cm$^2$ for 1 GeV, and $6.9 \times 10^{-24}$ cm$^2$ for 1 TeV dark matter mass. Finally, we discuss the implications of this study for dark matter physics and future observations.
We present calibrations of the redshift distributions of redMaGiC galaxies in the Dark Energy Survey Year 1 (DES Y1) and Sloan Digital Sky Survey (SDSS) DR8 data. These results determine the priors of the redshift distribution of redMaGiC galaxies, which were used for galaxy clustering measurements and as lenses for galaxy-galaxy lensing measurements in DES Y1 cosmological analyses. We empirically determine the bias in redMaGiC photometric redshift estimates using angular cross-correlations with Baryon Oscillation Spectroscopic Survey (BOSS) galaxies. For DES, we calibrate a single parameter redshift bias in three photometric redshift bins: $z \in[0.15,0.3]$, [0.3,0.45], and [0.45,0.6]. Our best fit results in each bin give photometric redshift biases of $|\Delta z|<0.01$. To further test the redMaGiC algorithm, we apply our calibration procedure to SDSS redMaGiC galaxies, where the statistical precision of the cross-correlation measurement is much higher due to a greater overlap with BOSS galaxies. For SDSS, we also find best fit results of $|\Delta z|<0.01$. We compare our results to other analyses of redMaGiC photometric redshifts.
We investigate the possibility of using future photometric and radio surveys to constrain the power spectrum of primordial fluctuations that is predicted by inflationary models with a violation of the slow-roll phase. We forecast constraints with a Fisher analysis on the amplitude of the parametrized features on ultra-large scales, in order to assess whether these could be distinguishable over the cosmic variance. We find that the next generation of photometric and radio surveys has the potential to test these models at a sensitivity better than current CMB experiments and that the synergy between galaxy and CMB observations is able to constrain models with many extra parameters. In particular, an SKA continuum survey with a huge sky coverage and a flux threshold of a few $\mu$Jy could confirm the presence of a new phase in the early Universe at more than 3$\sigma$.
We apply the Minkowski Tensor statistics to two dimensional slices of the three dimensional density field. The Minkowski Tensors are a set of functions that are sensitive to directionally dependent signals in the data, and furthermore can be used to quantify the mean shape of density peaks. We begin by introducing our algorithm for constructing bounding perimeters around subsets of a two dimensional field, and reviewing the definition of Minkowski Tensors. Focusing on the translational invariant statistic $W^{1,1}_{2}$ - a $2 \times 2$ matrix - we calculate its eigenvalues for both the entire excursion set ($\Lambda_{1},\Lambda_{2}$) and for individual connected regions and holes within the set ($\lambda_{1},\lambda_{2}$). The ratio of eigenvalues $\Lambda_{2}/\Lambda_{1}$ informs us of the presence of global anisotropies in the data, and $\langle \lambda_{2}/\lambda_{1} \rangle$ is a measure of the mean shape of peaks and troughs in the density field. We study these quantities for a Gaussian field, then consider how they are modified by the effect of gravitational collapse using the latest Horizon Run 4 cosmological simulation. We find $\Lambda_{1,2}$ are essentially independent of gravitational collapse, as the process maintains statistical isotropy. However, the mean shape of peaks is modified significantly - overdensities become relatively more circular compared to underdensities of the same area. When applying the statistic to a redshift space distorted density field, we find a significant signal in the eigenvalues $\Lambda_{1,2}$, suggesting that they can be used to probe the large-scale velocity field.
We present the first direct measurements of cubic bias parameters of dark matter halos from the halo-matter-matter-matter trispectrum. We measure this statistic efficiently by cross-correlating the halo field measured in N-body simulations with specific third-order nonlocal transformations of the initial density field in the same simulation. We also obtain constraints on the quadratic bias parameters. For the quadratic tidal bias ($b_{K^2}$), we find broad consistency with the prediction of the Lagrangian local-in-matter-density ansatz, with a small but significant detection of a negative Lagrangian tidal bias. Interesting trends are further found for the remaining three cubic bias parameters involving the tidal field.
We report constraints on cosmological parameters from the angular power spectrum of a cosmic microwave background (CMB) gravitational lensing potential map created using temperature data from 2500 deg$^2$ of South Pole Telescope (SPT) data supplemented with data from Planck in the same sky region, with the statistical power in the combined map primarily from the SPT data. We fit the corresponding lensing angular power spectrum to a model including cold dark matter and a cosmological constant ($\Lambda$CDM), and to models with single-parameter extensions to $\Lambda$CDM. We find constraints that are comparable to and consistent with constraints found using the full-sky Planck CMB lensing data. Specifically, we find $\sigma_8 \Omega_{\rm m}^{0.25}=0.598 \pm 0.024$ from the lensing data alone with relatively weak priors placed on the other $\Lambda$CDM parameters. In combination with primary CMB data from Planck, we explore single-parameter extensions to the $\Lambda$CDM model. We find $\Omega_k = -0.012^{+0.021}_{-0.023}$ or $M_{\nu}< 0.70$eV both at 95% confidence, all in good agreement with results that include the lensing potential as measured by Planck over the full sky. We include two independent free parameters that scale the effect of lensing on the CMB: $A_{L}$, which scales the lensing power spectrum in both the lens reconstruction power and in the smearing of the acoustic peaks, and $A^{\phi \phi}$, which scales only the amplitude of the CMB lensing reconstruction power spectrum. We find $A^{\phi \phi} \times A_{L} =1.01 \pm 0.08$ for the lensing map made from combined SPT and Planck temperature data, indicating that the amount of lensing is in excellent agreement with what is expected from the observed CMB angular power spectrum when not including the information from smearing of the acoustic peaks.
The cosmological redshift drift could lead to the next step in high-precision cosmic geometric observations, becoming a direct and irrefutable test for cosmic acceleration. In order to test the viability and possible properties of this effect, also called Sandage-Loeb (SL) test, we generate a model independent mock data set so as to compare its constraining power with that of the future mock data sets of Type Ia Supernovae (SNe) and Baryon Acoustic Oscillations (BAO). The performance of those data sets is analyzed by testing several cosmological models with the Markov chain Monte Carlo (MCMC) method, both independently and combining all data sets. Final results show that, in general, SL data sets allow for remarkable constraints on the matter density parameter today $\Omega_m$ on every tested model, showing also a great complementarity with SNe and BAO data regarding dark energy (DE) parameters.
Redshift space distortions within voids provide a unique method to test for environmental dependence of the growth rate of structures in low density regions, where effects of modified gravity theories might be important. We derive a linear theory model for the redshift space void-galaxy correlation that is valid at all pair separations, including deep within the void, and use this to obtain expressions for the monopole $\xi^s_0$ and quadrupole $\xi^s_2$ contributions. Our derivation highlights terms that have previously been neglected but are important within the void interior. As a result our model differs from previous works and predicts new physical effects, including a change in the sign of the quadrupole term within the void radius. We show how the model can be generalised to include a velocity dispersion. We compare our model predictions to measurements of the correlation function using mock void and galaxy catalogues modelled after the BOSS CMASS galaxy sample using the Big MultiDark N-body simulation, and show that the linear model with dispersion provides an excellent fit to the data at all scales, $0\leq s\leq120\;h^{-1}$Mpc. While the RSD model matches simulations, the linear bias approximation does not hold within voids, and care is needed in fitting for the growth rate. Assuming a constant velocity dispersion, fits to the quadrupole recover the growth rate $f$ to a precision of $4\%$ using the simulation volume of $(2.5\;h^{-1}\mathrm{Gpc})^3$ at $z=0.52$.
Prior to photon decoupling magnetic random fields of comoving intensity in the nano-Gauss range distort the temperature and the polarization anisotropies of the microwave background, potentially induce a peculiar B-mode power spectrum and may even generate a frequency-dependent circularly polarized V-mode. We critically analyze the theoretical foundations and the recent achievements of an interesting trialogue involving plasma physics, general relativity and astrophysics.
The $f(R)$ gravity and scalar-tensor theory are known to be equivalent at the classical level. We study if this equivalence is valid at the quantum level. There are two descriptions of the scalar-tensor theory in the Jordan and Einstein frames. It is shown that these three formulations of the theories give the same determinant or effective action on shell, and thus they are equivalent at the quantum level on shell in arbitrary dimensions. We also compute the one-loop divergence in $f(R)$ gravity.
We implement inflation in a non-supersymmetric SU(5) model based on a non-minimal coupling of the axion field to gravity. The isocurvature fluctuations are adequately suppressed, axions comprise the dark matter, proton lifetime estimates are of order 8e34-3e35 years, and the observed baryon asymmetry arises via non-thermal leptogenesis. The presence of low scale colored scalars ensures unification of the Standard Model gauge couplings and also helps in stabilizing the electroweak vacuum.
We use the semi-analytical model of galaxy formation GALFORM to characterise an indirect signature of AGN feedback in the environment of radio galaxies at high redshifts. The predicted environment of radio galaxies is denser than that of radio-quiet galaxies with the same stellar mass. This is consistent with observational results from the CARLA survey. Our model shows that the differences in environment are due to radio galaxies being hosted by dark matter haloes that are ~1.5 dex more massive than those hosting radio-quiet galaxies with the same stellar mass. By running a control-simulation in which AGN feedback is switched-off, we identify AGN feedback as the primary mechanism affecting the build-up of the stellar component of radio galaxies, thus explaining the different environment in radio galaxies and their radio-quiet counterparts. The difference in host halo mass between radio loud and radio quiet galaxies translates into different galaxies populating each environment. We predict a higher fraction of passive galaxies around radio loud galaxies compared to their radio-quiet counterparts. Furthermore, such a high fraction of passive galaxies shapes the predicted infrared luminosity function in the environment of radio galaxies in a way that is consistent with observational findings. Our results suggest that the impact of AGN feedback at high redshifts and environmental mechanisms affecting galaxies in high halo masses can be revealed by studying the environment of radio galaxies, thus providing new constraints on galaxy formation physics at high redshifts.
21 cm cosmology is a promising new probe of the evolution of visible matter in our universe, especially during the poorly-constrained Cosmic Dawn and Epoch of Reionization. However, in order to separate the 21 cm signal from bright astrophysical foregrounds, we need an exquisite understanding of our telescopes so as to avoid adding spectral structure to spectrally-smooth foregrounds. One powerful calibration method relies on repeated simultaneous measurements of the same interferometric baseline to solve for the sky signal and for instrumental parameters simultaneously. However, certain degrees of freedom are not constrained by asserting internal consistency between redundant measurements. In this paper, we review the origin of these "degeneracies" of redundant-baseline calibration, demonstrate how they can source unwanted spectral structure in our measurements, and show how to eliminate that structure. We also generalize redundant calibration to dual-polarization instruments, derive the new degeneracies, and explore the unique challenges to calibration and spectral smoothness presented by a polarized measurement.
Extending our previous work on the robustness of inflation to perturbations in the scalar field, we investigate the effects of perturbations in the transverse traceless part of the extrinsic curvature on the evolution of an inhomogeneous inflaton field. Focusing on small field models, we show that these additional metric inhomogeneities initially reduce the total number of e-folds as the amplitude increases, but that the reduction saturates and even reverses above a certain amplitude. We present an argument that this is due to the presence of a large initial Hubble friction when metric perturbations are large.
We aim to present 70 spectra of 68 new high-redshift type Ia supernovae
(SNeIa) measured at ESO's VLT during the final two years of operation
(2006-2008) of the Supernova Legacy Survey (SNLS). We use the full five year
SNLS VLT spectral set to investigate a possible spectral evolution of SNeIa
populations with redshift and study spectral properties as a function of
lightcurve fit parameters and the mass of the host-galaxy.
Reduction and extraction are based on both IRAF standard tasks and our own
reduction pipeline. Redshifts are estimated from host-galaxy lines whenever
possible or alternatively from supernova features. We used the
spectrophotometric SNIa model SALT2 combined with a set of galaxy templates
that model the host-galaxy contamination to assess the type Ia nature of the
candidates.
We identify 68 new SNeIa with redshift ranging from z=0.207 to z=0.98
(<z>=0.62). Each spectrum is presented individually along with its best-fit
SALT2 model. The five year dataset contains 209 spectra corresponding to 192
SNeIa identified at the VLT. We also publish the redshifts of other candidates
(host galaxies or other transients) whose spectra were obtained at the same
time as the spectra of live SNe Ia. Using the full VLT SNeIa sample, we build
composite spectra around maximum light with cuts in color, lightcurve shape
parameter ('stretch'), host-galaxy mass and redshift. We find that high-z SNeIa
are bluer, brighter and have weaker intermediate mass element absorption lines
than their low-z counterparts at a level consistent with what is expected from
selection effects. We also find a flux excess in the range [3000-3400] A for
SNeIa in low mass host-galaxies or with locally blue U-V colors, and suggest
that the UV flux (or local color) may be used in future cosmological studies as
a third standardization parameter in addition to stretch and color.
If the dark matter is produced in the early universe prior to Big Bang nucleosynthesis, a modified cosmological history can drastically affect the abundance of relic dark matter particles. Here, we assume that an additional species to radiation dominates at early times, causing the expansion rate at a given temperature to be larger than in the standard radiation-dominated case. We demonstrate that, if this is the case, dark matter production via freeze-in (a scenario when dark matter interacts very weakly, and is dumped in the early universe out of equilibrium by decay or scattering processes involving particles in the thermal bath) is dramatically suppressed. We illustrate and quantitatively and analytically study this phenomenon for three different paradigmatic classes of freeze-in scenarios. For the frozen-in dark matter abundance to be as large as observations, couplings between the dark matter and visible-sector particles must be enhanced by several orders of magnitude. This sheds some optimistic prospects for the otherwise dire experimental and observational outlook of detecting dark matter produced by freeze-in.
Axion-like particles (ALPs) are interesting dark matter candidates both from the theoretical as well as from the experimental perspective. Usually they are motivated as pseudo-Nambu-Goldstone bosons. In this case one of their most important features is that their coupling to other particles is suppressed by a large scale, the vacuum expectation value of the field breaking the symmetry that gives rise to them. This naturally endows them with very weak interactions but also restricts the maximal field value and therefore the regions where sufficient dark matter is produced. In this paper we investigate deviations from this simplest setup, where the potential and interactions are as expected for a pseudo-Nambu-Goldstone boson, but the kinetic term has singularities. This leads to a significantly increased area in parameter space where such particles can be dark matter and can be probed by current and near future experiments. We discuss cosmological limits and in the course of this give a simple derivation of a formula for isocurvature fluctuations in models with general anharmonic potentials. As an application of this formula we give an update of the isocurvature constraints for QCD axion dark matter models, using the most recent results for the QCD topological susceptibility and the newest Planck data.
The cosmological model is at present not tested between the redshift of the farthest observed supernovae (z ~ 1.4) and that of the Cosmic Microwave Background (z ~ 1,100). Here we introduce a new method to measure the cosmological parameters: we show that quasars can be used as "standard candles" by employing the non-linear relation between their intrinsic UV and X-ray emission as an absolute distance indicator. We built a sample of ~ 1,900 quasars with available UV and X-ray observations, and produced a Hubble Diagram up to z ~ 5. The analysis of the quasar Hubble Diagram, when used in combination with supernovae, provides robust constraints on the matter and energy content in the cosmos. The application of this method to forthcoming, larger quasar samples, will also provide tight constraints on the dark energy equation of state and its possible evolution with time.
Links to: arXiv, form interface, find, astro-ph, recent, 1712, contact, help (Access key information)
Primordial non-Gaussianities enhanced at small wavevectors can induce a power spectrum of the galaxy overdensity that differs greatly from that of the matter overdensity at large length scales. In previous work, it was shown that "squeezed" three-point and "collapsed" four-point functions of the curvature perturbation $\zeta$ can generate these non-Gaussianities and give rise to so-called scale-dependent and stochastic bias in the galaxy overdensity power spectrum. We explore a third way to generate non-Gaussianities enhanced at small wavevectors: the infrared behavior of quantum loop contributions to the four-point correlations of $\zeta$. We show that these loop effects lead to stochastic bias, which can be observable in the context of quasi-single field inflation.
The dispersion relation of graviton is a fundamental issue for fundamental physics about gravity. In this paper we investigate how the modified dispersion relation of graviton affects the cosmic microwave background (CMB) power spectra, in particular the B-mode polarization. Our results will be useful to test the dispersion relation of graviton at the energy scale around $10^{-29}$ eV.
The lack of large-angle correlations in the fluctuations of the cosmic microwave background (CMB) conflicts with predictions of slow-roll inflation. But though probabilities (< 0.24%) for the missing correlations disfavor the conventional picture at > 3 sigma, factors not associated with the model itself may be contributing to the tension. Here we aim to show that the absence of large-angle correlations is best explained with the introduction of a non-zero minimum wavenumber k_min for the fluctuation power spectrum P(k). We assume that quantum fluctuations were generated in the early Universe with a well-defined power spectrum P(k), though with a cutoff k_min not equal to 0. We then re-calculate the angular correlation function of the CMB and compare it with Planck observations. The Planck 2013 data rule out a zero k_min at a confidence level exceeding 8 sigma. Whereas purely slow-roll inflation would have stretched all fluctuations beyond the horizon, producing a P(k) with k_min=0---and therefore strong correlations at all angles---a k_min > 0 would signal the presence of a maximum wavelength at the time (t_dec) of decoupling. This argues against the basic inflationary paradigm---perhaps even suggesting non-inflationary alternatives---for the origin and growth of perturbations in the early Universe. In at least one competing cosmology, the R_h=ct universe, the inferred k_min corresponds to the gravitational radius at t_dec.
The extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 14 sample includes 80,118 Luminous Red Galaxies. By combining these galaxies with the high-redshift tail of the BOSS galaxy sample, we form a sample of LRGs at an effective redshift $z=0.72$, covering an effective volume of 0.9~Gpc$^3$. We introduce new techniques to account for spurious fluctuations caused by targeting and by redshift failures which were validated on a set of mock catalogs. This analysis is sufficient to provide a $2.6$\% measurement of spherically averaged BAO, $D_V(z=0.72) = 2353^{+63}_{-61} (r_d/r_{d,\rm{fid}}) h^{-1}$Mpc, at 2.8$\sigma$ of significance. Together with the recent quasar-based BAO measurement at $z=1.5$, and forthcoming Emission Line Galaxy-based measurements, this measurement demonstrates that eBOSS is fulfilling its remit of extending the range of redshifts covered by such measurements, laying the ground work for forthcoming surveys such as the Dark Energy Spectroscopic Survey and Euclid.
Light thermal relics of the hot big bang, often quantified by the parameter $N_\mathrm{eff}$, are one of the primary targets of cosmological measurements. At present, the energy density in such relics is constrained to be less than ten percent of the total energy density in radiation. Upcoming cosmic microwave background (CMB) experiments, however, have the potential to measure the radiation density at the one-percent level, which is close to well-motivated theoretical targets. In this paper, we explore to what degree the CMB observations can be enhanced by future large-scale structure surveys. We carefully isolate the information encoded in the shape of the galaxy power spectrum and in the spectrum of baryon acoustic oscillations (BAO). We find that measurements of the shape of the power spectrum can significantly improve on current and near-term CMB experiments. We also show that the phase shift of the BAO spectrum induced by relic neutrinos can be detected at high significance in future experiments.
In modified gravity theories, the luminosity distance for gravitational waves (GWs) is in general different from that for electromagnetic signals. This opens the possibility of model-independent tests of GR, by comparing the results from standard candles with those from standard sirens. We illustrate this fact using a nonlocal modification of gravity, that has been shown to fit remarkably well CMB, SNe, BAO and structure formation data. We find that, for rather general reasons, the relative difference between the GW luminosity distance of the modified gravity model and the luminosity distance of $\Lambda$CDM is larger by one order of magnitude than the relative difference of electromagnetic luminosity distances. This effect would be missed in phenomenological studies that parametrize dark energy using only an equation of state $w_{\rm DE}(z)$. We discuss the prospects for distinguishing nonlocal gravity from $\Lambda$CDM at third-generation GW detectors such as the Einstein Telescope and we find that, depending on the exact sensitivity, a few tens of standard sirens with measured redshift at $z\sim 0.4$, or a few hundreds at $1< z < 2$, could suffice.
In this paper, we introduce the topological persistence diagram as a statistic for Cosmic Microwave Background (CMB) temperature anisotropy maps. A central concept in `Topological Data Analysis' (TDA), the idea of persistence is to represent a data set by a family of topological spaces. One then examines how long topological features `persist' as the family of spaces is traversed. We compute persistence diagrams for simulated CMB temperature anisotropy maps featuring various levels of primordial non-Gaussianity of local type. Postponing the analysis of observational effects, we show that persistence diagrams are more sensitive to local non-Gaussianity than previous topological statistics including the genus and Betti number curves, and can constrain $\Delta f_{NL}^{\rm loc}= 35.8$ at the 68\% confidence level on the simulation set, compared to $\Delta f_{NL}^{\rm loc}= 60.6$ for the Betti number curves. Given the resolution of our simulations, we expect applying persistence diagrams to observational data will give constraints competitive with those of the Minkowski Functionals. This is the first in a series of papers where we plan to apply TDA to different shapes of non-Gaussianity in the CMB and Large Scale Structure.
We construct an effective field theory (EFT) description of the hard photon spectrum for heavy WIMP annihilation. This facilitates precision predictions relevant for line searches, and allows the incorporation of non-trivial energy resolution effects. Our framework combines techniques from non-relativistic EFTs and soft-collinear effective theory (SCET), as well as its multi-scale extensions that have been recently introduced for studying jet substructure. We find a number of interesting features, including the simultaneous presence of SCET$_{\text{I}}$ and SCET$_{\text{II}}$ modes, as well as collinear-soft modes at the electroweak scale. We derive a factorization formula that enables both the resummation of the leading large Sudakov double logarithms that appear in the perturbative spectrum, and the inclusion of Sommerfeld enhancement effects. Consistency of this factorization is demonstrated to leading logarithmic order through explicit calculation. Our final result contains both the exclusive and the inclusive limits, thereby providing a unifying description of these two previously-considered approximations. We estimate the impact on experimental sensitivity, focusing for concreteness on an SU(2)$_{W}$ triplet fermion dark matter - the pure wino - where the strongest constraints are due to a search for gamma-ray lines from the Galactic Center. We find numerically significant corrections compared to previous results, thereby highlighting the importance of accounting for the photon spectrum when interpreting data from current and future indirect detection experiments.
We present a purely gravitational infra-red-calculable production mechanism for dark matter (DM). The source of both the DM relic abundance and the hot Standard Model (SM) plasma is a primordial density of micro black holes (BHs), which evaporate via Hawking emission into both the dark and SM sectors. The mechanism has four qualitatively different regimes depending upon whether the BH evaporation is `fast' or `slow' relative to the initial Hubble rate, and whether the mass of the DM particle is `light' or `heavy' compared to the initial BH temperature. For each of these regimes we calculate the DM yield, $Y$, as a function of the initial state and DM mass and spin. In the `slow' regime $Y$ depends on only the initial BH mass over a wide range of initial conditions, including scenarios where the BHs are a small fraction of the initial energy density. The DM is produced with a highly non-thermal energy spectrum, leading in the `light' DM mass regime ($\sim260\,\mathrm{eV}$ and above depending on DM spin) to a strong constraint from free-streaming, but also possible observational signatures in structure formation in the spin 3/2 and 2 cases. The `heavy' regime ($\sim1.2\times 10^8\,\mathrm{GeV}$ to $M_{\mathrm{Pl}}$ depending on spin) is free of these constraints and provides new possibilities for DM detection. In all cases there is a dark radiation component predicted.
We present the results of a new Chandra study of the galaxy cluster A2626. The radio emission of the cluster shows a complex system of four symmetric arcs without known correlations with the X-ray emission. The mirror symmetry of the radio arcs toward the center and the presence of two optical cores in the central galaxy suggested that they may be created by pairs of precessing radio jets powered by dual AGNs inside the cD galaxy. However, previous observations failed to observe the second jetted AGN and the spectral trend due to radiative age along the radio arcs, thus challenging this interpretation. The new Chandra observation had several scientific objectives, including the search for the second AGN that would support the jet precession model. We focus here on the detailed study of the local properties of the thermal and non-thermal emission in the proximity of the radio arcs, in order to get more insights into their origin. We performed a standard data reduction of the Chandra dataset deriving the radial profiles of temperature, density, pressure and cooling time of the intra-cluster medium. We further analyzed the 2D distribution of the gas temperature, discovering that the south-western junction of the radio arcs surrounds the cool core of the cluster. We studied the X-ray SB and spectral profiles across the junction, finding a cold front spatially coincident with the radio arcs. This may suggest a connection between the sloshing of the thermal gas and the nature of the radio filaments, raising new scenarios for their origin. A possibility is that the radio arcs trace the projection of a complex surface connecting the sites where electrons are most efficiently reaccelerated by the turbulence that is generated by the gas sloshing. In this case, diffuse emission embedded by the arcs and with extremely steep spectrum should be most visible at very low radio frequencies.
The discovery of dark matter (DM) at XENONnT or LZ would place constraints on DM particle mass and coupling constants. It is interesting to ask when these constraints can be compatible with the DM thermal production mechanism. We address this question within the most general set of renormalisable models that preserve Lorentz and gauge symmetry, and that extend the Standard Model by one DM candidate of mass $m_{\rm DM}$ and one particle of mass $M_{\rm med}$ mediating DM-quark interactions. Our analysis divides into two parts. First, we postulate that XENONnT/LZ has detected $\mu_S\sim\mathcal{O}(100)$ signal events, and use this input to calculate the DM relic density, $\Omega_{\rm DM} h^2$. Then, we identify the regions in the $M_{\rm med} - \Omega_{\rm DM} h^2$ plane which are compatible with the observed signal and with current CMB data. We find that for most of the models considered here, $\mathcal{O}(100)$ signal events at XENONnT/LZ and the DM thermal production are only compatible for resonant DM annihilations, i.e. for $M_{\rm med}\simeq2 m_{\rm DM}$. In this case, XENONnT/LZ would be able to simultaneously measure $m_{\rm DM}$ and $M_{\rm med}$. We also discuss the dependence of our results on $m_{\rm DM}$, $\mu_S$ and the DM spin, and provide analytic expressions for annihilation cross-sections and mediator decay widths for all models considered in this study.
We consider the effect of the Gibbons-Hawking radiation on the inflaton in the situation where it is coupled to a large number of spectator fields. We argue that this will lead to two important effects - a thermal contribution to the potential and a gradual change in parameters in the Lagrangian which results from thermodynamic and energy conservation arguments. We present a scenario of hilltop inflation where the field starts trapped at the origin before slowly experiencing a phase transition during which the field extremely slowly moves towards its zero temperature expectation value. We show that it is possible to obtain enough e-folds of expansion as well as the correct spectrum of perturbations without hugely fine-tuned parameters in the potential (albeit with many spectator fields). We also comment on how initial conditions for inflation can arise naturally in this situation.
Supersymmetry breaking close to the scale of grand unification can explain cosmic inflation. As we demonstrate in this paper, this can be achieved in strongly coupled supersymmetric gauge theories, such that the energy scales of inflation and supersymmmetry breaking are generated dynamically. As a consequence, both scales are related to each other and exponentially suppressed compared to the Planck scale. As an example, we consider a dynamical model in which gauging a global flavor symmetry in the supersymmetry-breaking sector gives rise to a Fayet-Iliopoulos D term. This results in successful D-term hybrid inflation in agreement with all theoretical and phenomenological constraints. The gauged flavor symmetry can be identified with $U(1)_{B-L}$, where $B$ and $L$ denote baryon and lepton number, respectively. In the end, we arrive at a consistent cosmological scenario that provides a unified picture of high-scale supersymmetry breaking, viable D-term hybrid inflation, spontaneous $B$$-$$L$ breaking at the scale of grand unification, baryogenesis via leptogenesis, and standard model neutrino masses due to the type-I seesaw mechanism.
We extend the phenomenology of loop quantum cosmology (LQC) to second order in perturbations. Our motivation is twofold. On the one hand, since LQC predicts a cosmic bounce that takes place at the Planck scale, the second order contributions could be large enough to jeopardize the validity of the perturbative expansion on which previous results rest. On the other hand, the upper bounds on primordial non-Gaussianity obtained by the Planck Collaboration are expected to play a significant role on explorations of the LQC phenomenology. We find that the bounce in LQC produces an enhancement of non-Gaussianity of several orders of magnitude, on length scales that were larger than the curvature radius at the bounce. Nonetheless, we find that one can still rely on the perturbative expansion to make predictions about primordial perturbations. We discuss the consequences of our results for LQC and its predictions for the cosmic microwave background.
Links to: arXiv, form interface, find, astro-ph, recent, 1712, contact, help (Access key information)