Critical overdensity $\delta_c$ is a key concept in estimating the number count of halos for different redshift and halo-mass bins, and therefore, it is a powerful tool to compare cosmological models to observations. There are currently two different prescriptions in the literature for its calculation, namely, the differential-radius and the constant-infinity methods. In this work we show that the latter yields precise results {\it only} if we are careful in the definition of the so-called numerical infinities. Although the subtleties we point out are crucial ingredients for an accurate determination of $\delta_c$ both in general relativity and in any other gravity theory, we focus on $f(R)$ modified-gravity models in the metric approach; in particular, we use the so-called large ($F=1/3$) and small-field ($F=0$) limits. For both of them, we calculate the relative errors (between our method and the others) in the critical density $\delta_c$, in the comoving number density of halos per logarithmic mass interval $n_{\ln M}$ and in the number of clusters at a given redshift in a given mass bin $N_{\rm bin}$, as functions of the redshift. We have also derived an analytical expression for the density contrast in the linear regime as a function of the collapse redshift $z_c$ and $\Omega_{m0}$ for any $F$.
Discordance in the $\Lambda$CDM cosmological model can be seen by comparing parameters constrained by CMB measurements to those inferred by probes of large scale structure. Recent improvements in observations, including final data releases from both Planck and SDSS-III BOSS, as well as improved astrophysical uncertainty analysis of CFHTLenS, allows for an update in the quantification of any tension between large and small scales. This paper is intended, primarily, as a discussion on the quantifications of discordance when comparing the parameter constraints of a model when given two different data sets. We consider KL-divergence, comparison of Bayesian evidences and other statistics which are sensitive to the mean, variance and shape of the distributions. However, as a by-product, we present an update to the similar analysis in (Battye, Charnock and Moss; 2015) where we find that, considering new data and treatment of priors, the constraints from the CMB and from a combination of LSS probes are in greater agreement and any tension only persists to a minor degree. In particular, we find the parameter constraints from the combination of LSS probes which are most discrepant with the Planck2015+Pol+BAO parameter distributions can be quantified at a 2.55$\sigma$ tension using the method introduced in (Battye, Charnock and Moss; 2015). If instead we use the distributions constrained by the combination of LSS probes which are in greatest agreement with those from Planck2015+Pol+BAO this tension is only 0.76$\sigma$.
We introduce the Lomonosov suite of high-resolution N-body cosmological simulations covering a full box of size 32 $h^{-1}$ Mpc with low-mass resolution particles ($2\times10^7$ $h^{-1}$ M$_\odot$) and three zoom-in simulations of overdense, underdense and mean density regions at much higher particle resolution ($4\times10^4$ $h^{-1}$ M$_\odot$). The main purpose of this simulation suite is to extend the concentration-mass relation of dark matter halos down to masses below those typically available in large cosmological simulations. The three different density regions available at higher resolution provide a better understanding of the effect of the local environment on halo concentration, known to be critical for small simulation boxes. Indeed, we show that the smaller the halo mass in our simulations the more important becomes to account for environmental effects if we aim for an accurate determination of halo concentrations. A precise characterization of this effect allows us to develop a robust technique to reliably extrapolate the concentration values found in zoom simulations to larger volumes. All together, Lomonosov provides a measure of the concentration-mass relation in the halo mass range $10^7-10^{10}$ $h^{-1}$ M$_\odot$ with superb halo statistics. This work represents a first important step to measure halo concentrations at intermediate, yet vastly unexplored halo mass scales, down to the smallest ones. All Lomonosov data and files are public for community's use.
We addressed the so far unexplored issue of outflows induced by exponentially
growing power sources, focusing on early supermassive black holes (BHs). We
assumed that these objects grow to $10^9\;M_{\odot}$ by z=6 by
Eddington-limited accretion and convert 5% of their bolometric output into a
wind. We first considered the case of energy-driven and momentum-driven
outflows expanding in a region where the gas and total mass densities are
uniform and equal to the average values in the Universe at $z>6$. We derived
analytic solutions for the evolution of the outflow, finding that, for an
exponentially growing power with e-folding time $t_{Sal}$, the late time
expansion of the outflow radius is also exponential, with e-folding time of
$5t_{Sal}$ and $4t_{Sal}$ in the energy-driven and momentum-driven limit,
respectively.
We then considered energy-driven outflows produced by QSOs at the center of
early dark matter halos of different masses and powered by BHs growing from
different seeds. We followed the evolution of the source power and of the gas
and dark matter density profiles in the halos from the beginning of the
accretion until $z=6$. The final bubble radius and velocity do not depend on
the seed BH mass but are instead smaller for larger halo masses. At z=6, bubble
radii in the range 50-180 kpc and velocities in the range 400-1000 km s$^{-1}$
are expected for QSOs hosted by halos in the mass range
$3\times10^{11}-10^{13}\;M_{\odot}$.
By the time the QSO is observed, we found that the total thermal energy
injected within the bubble in the case of an energy-driven outflow is
$E_{th}\sim5 \times 10^{60}$ erg. This is in excellent agreement with the value
of $E_{th}=(6.2\pm 1.7)\times 10^{60}$ erg measured through the detection of
the thermal Sunyaev-Zeldovich effect around a large population of luminous QSOs
at lower redshift. [abridged]
We investigate the behaviour of bouncing Bianchi type IX `Mixmaster' universes in general relativity. This generalises all previous studies of the cyclic behaviour of closed spatially homogeneous universes with and without entropy increase. We determine the behaviour of models containing radiation by analytic and numerical integration and show that increase of radiation entropy leads to increasing cycle size and duration. We introduce a null energy condition violating ghost field to create a smooth, non-singular bounce of finite size at the end of each cycle and compute the evolution through many cycles with and without entropy increase injected at the start of each cycle. In the presence of increasing entropy we find that the cycles grow larger and longer and the dynamics approach flatness, as in the isotropic case. However, successive cycles become increasingly anisotropic at the expansion maxima which is dominated by the general-relativistic effects of anisotropic 3-curvature. However, it becomes positive after expansion drives the dynamics close enough to isotropy for the curvature to become positive and for gravitational collapse to ensue. In the presence of a positive cosmological constant, radiation and a ghost field we show that, for a very wide range of cosmological constant values, the growing oscillations always cease and the dynamics subsequently approach those of the isotropic de Sitter universe at late times. This model is not included in the scope of earlier cosmic no-hair theorems because the 3-curvature can be positive. In the case of negative cosmological constant, radiation and an ultra-stiff field (to create non-singular bounces) we show that a sequence of chaotic oscillations also occurs, with sensitive dependence on initial conditions. In all cases, we follow the oscillatory evolution of the scale factors, the shear, and the 3-curvature from cycle to cycle.
We study generalized $\alpha$-attractor models whose rescaled scalar manifold is the triply-punctured Riemann sphere $Y(2)$ endowed with its complete hyperbolic metric. Using an explicit embedding into the end compactification, we compute solutions of the cosmological evolution equations for a few globally well-behaved scalar potentials, displaying particular trajectories with inflationary behavior. In such models, the orientation-preserving isometry group of the scalar manifold is isomorphic with the permutation group on three elements, acting on $Y(2)$ as the group of anharmonic transformations. When the scalar potential is preserved by this action, $\alpha$-attractor models of this type provide a geometric description of two-field "modular invariant $j$-models" in terms of gravity coupled to a non-linear sigma model with topologically non-trivial target and with a finite (as opposed to discrete but infinite) group of symmetries. The relation between the two perspectives is provided by the elliptic modular function $\lambda$, viewed as a field redefinition which eliminates most of the infinite unphysical ambiguity present in the Poincare half-plane description of such models.
Links to: arXiv, form interface, find, astro-ph, recent, 1703, contact, help (Access key information)
We show how LIGO is expected to detect coalescing binary black holes at $z>1$, that are lensed by the intervening galaxy population. Gravitational magnification, $\mu$, strengthens gravitational wave signals by $\sqrt{\mu}$, without altering their frequencies, which if unrecognised leads to an underestimate of the event redshift and hence an overestimate of the binary mass. High magnifications can be reached for coalescing binaries because the region of intense gravitational wave emission during coalescence is so small ($\sim$100km), permitting very close projections between lensing caustics and gravitational-wave events. Our simulations incorporate accurate waveforms convolved with the LIGO power spectral density. Importantly, we include the detection dependence on sky position and orbital orientation, which for the LIGO configuration translates into a wide spread in observed redshifts and chirp masses. Currently we estimate a detectable rate of lensed events \rateEarly{}, that rises to \rateDesign{}, at LIGO's design sensitivity limit, depending on the high redshift rate of black hole coalescence.
The colliding cluster, CIZA J2242.8+5301, displays a spectacular, almost 2 Mpc long shock front with a radio based Mach number M ~ 5, that is puzzlingly large compared with the X-ray estimate of M ~ 2.5. The extent to which the X-ray temperature jump is diluted by cooler unshocked gas projected through the cluster currently lacks quantification. Thus, here we apply our self-consistent N-body/hydro-dynamical code (based on FLASH) to model this binary cluster encounter. We can account for the location of the shock front and also the elongated X-ray emission by tidal stretching of the gas and dark matter between the two cluster centers. The required total mass is $8.9 \times 10^{14}$ Msun with a 1.3:1 mass ratio favoring the southern cluster component. The relative velocity we derive is $\simeq 2500$ km/s initially between the two main cluster components, with an impact parameter of 120 kpc. This solution implies that the shock temperature jump derived from the low angular resolution X-ray satellite SUZAKU is underestimated by a factor of two, due to cool gas in projection, bringing the observed X-ray and radio estimates into agreement. We propose that the complex southern relics in CIZA J2242.8+5301, have been broken up as the southerly moving "back" shocked gas impacts the gas still falling in along the collision axis. Finally, we use our model to generate Compton-y maps to estimate the reduction in radio flux caused by the thermal Sunyaev-Zel'dovich (SZ) effect. At 30 GHz, this amounts to $\Delta S_n = -0.072$ mJy/arcmin$^2$ and $\Delta S_s = -0.075$ mJy/arcmin$^2$ at the locations of the northern and southern shock fronts respectively. Our model estimate agrees with previous empirical estimates that have inferred the measured radio spectra can be significantly affected by the SZ effect, with implications for charged particle acceleration models of the radio relics.
IC 1613 is an isolated dwarf galaxy within the Local Group. Low foreground and internal extinction, low metallicity, and low crowding make it an invaluable testbed for the calibration of the local distance ladder. We present here new, high-fidelity distance estimates to IC 1613 via its Tip of the Red Giant Branch (TRGB) and its RR Lyrae (RRL) variables as part of the Carnegie-Chicago Hubble Program, which seeks an alternate local route to $H_0$ using Population II stars. We have measured an extinction-corrected $I_{ACS}^{\mathrm{TRGB}} = 20.32 \pm 0.02_{stat} \pm 0.01_{sys}$ mag using wide-field observations obtained from the IMACS camera on the Magellan-Baade telescope. We have further constructed $VIH$ RRL period-luminosity relations using archival and new observations from the ACS/WFC and WFC3/IR instruments aboard the $Hubble~Space~Telescope~(HST)$. In advance of future $Gaia$ data releases, we set provisional values for the TRGB luminosity via the Large Magellanic Cloud and Galactic RRL zero-points via $HST$ parallaxes. We find corresponding true distance moduli $\mu_{I}^{\mathrm{TRGB}} = 24.27 \pm 0.05$ mag and $\langle\mu^{\mathrm{RRL}}\rangle = 24.30 \pm 0.07$ mag. We compare our results to a body of recent publications on IC 1613 and find no statistically significant difference between the distances derived from stars of Population I and II.
In this paper, we apply reconstruction techniques to recover the potential parameters for a particular class of single-field models, the $\alpha$-attractor (supergravity) models of inflation. This also allows to derive the inflaton vacuum expectation value at horizon crossing. We show how to use this value as one of the input variables to constrain the postaccelerated inflationary phase. We assume that the tensor-to-scalar ratio $r$ is of the order of $10^{-3}$ , a level reachable by the expected sensitivity of the next-generation CMB experiments.
Long wavelength spectral distortions in the Cosmic Microwave Background arising from the 21-cm transition in neutral Hydrogen are a key probe of Cosmic Dawn and the Epoch of Reionization. These features may reveal the nature of the first stars and ultra-faint galaxies that transformed the spin temperature and ionization state of the primordial gas. SARAS 2 is a spectral radiometer purposely designed for precision measurement of these monopole or all-sky global 21-cm spectral distortions. We use a 4-hr observation of the radio background in the frequency band 110-200 MHz with the radiometer deployed at the Timbaktu Collective in Southern India to derive likelihoods for plausible redshifted 21-cm signals predicted by theoretical models. First light with SARAS 2 disfavors models that feature weak X-ray heating along with rapid reionization.
The primordial deuterium abundance probes fundamental physics during the Big Bang Nucleosynthesis and can be used to infer cosmological parameters. Observationally, the abundance can be measured using absorbing clouds along the lines of sight to distant quasars. Observations of the quasar PKS1937--101 contain two absorbers for which the deuterium abundance has previously been determined. Here we focus on the higher redshift one at $z_{abs} = 3.572$. We present new observations with significantly increased signal-to-noise ratio which enable a far more precise and robust measurement of the deuterium to hydrogen column density ratio, resulting in D/H = $2.62\pm0.05\times10^{-5}$. This particular measurement is of interest because it is among the most precise assessments to date and it has been derived from the second lowest column-density absorber [N(HI) $=17.9\mathrm{cm}^{-2}$] that has so-far been utilised for deuterium abundance measurements. The majority of existing high-precision measurements were obtained from considerably higher column density systems [i.e. N(HI) $>19.4\mathrm{cm}^{-2}$]. This bodes well for future observations as low column density systems are more common.
We investigate possible signatures of halo assembly bias for spectroscopically selected galaxy groups from the GAMA survey using weak lensing measurements from the spatially overlapping regions of the deeper, high-imaging-quality photometric KiDS survey. We use GAMA groups with an apparent richness larger than 4 to identify samples with comparable mean host halo masses but with a different radial distribution of satellite galaxies, which is a proxy for the formation time of the haloes. We measure the weak lensing signal for groups with a steeper than average and with a shallower than average satellite distribution and find no sign of halo assembly bias, with the bias ratio of $0.85^{+0.37}_{-0.25}$, which is consistent with the $\Lambda$CDM prediction. Our galaxy groups have typical masses of $10^{13} M_{\odot}/h$, naturally complementing previous studies of halo assembly bias on galaxy cluster scales.
We investigate cosmic string networks in the Abelian Higgs model using data from a campaign of large-scale numerical simulations on lattices of up to $4096^3$ grid points. We observe scaling or self-similarity of the networks over a wide range of scales, and estimate the asymptotic values of the mean string separation in horizon length units $\dot{\xi}$ and of the mean square string velocity $\bar v^2$ in the continuum and large time limits. The scaling occurs because the strings lose energy into classical radiation of the scalar and gauge fields of the Abelian Higgs model. We quantify the energy loss with a dimensionless radiative efficiency parameter, and show that it does not vary significantly with lattice spacing or string separation. This implies that the radiative energy loss underlying the scaling behaviour is not a lattice artefact, and justifies the extrapolation of measured network properties to large times for computations of cosmological perturbations. We also show that the core growth method, which increases the defect core width with time to extend the dynamic range of simulations, does not introduce significant systematic error. We compare $\dot{\xi}$ and $\bar v^2$ to values measured in simulations using the Nambu-Goto approximation, finding that the latter underestimate the mean string separation by about 25%, and overestimate $\bar v^2$ by about 10%. The scaling of the string separation implies that string loops decay by the emission of massive radiation within a Hubble time in field theory simulations, in contrast to the Nambu-Goto scenario which neglects this energy loss mechanism. String loops surviving for only one Hubble time emit much less gravitational radiation than in the Nambu-Goto scenario, and are consequently subject to much weaker gravitational wave constraints on their tension.
We report on a spectral study at radio frequencies of the giant radio halo in A2142 (z=0.0909), which we performed to explore its nature and origin. A2142 is not a major merger and the presence of a giant radio halo is somewhat surprising. We performed deep radio observations with the GMRT at 608 MHz, 322 MHz, and 234 MHz and with the VLA in the 1-2 GHz band. We obtained high-quality images at all frequencies in a wide range of resolutions. The radio halo is well detected at all frequencies and extends out to the most distant cold front in A2142. We studied the spectral index in two regions: the central part of the halo and a second region in the direction of the most distant south-eastern cold front, selected to follow the bright part of the halo and X-ray emission. We complemented our observations with a preliminary LOFAR image at 118 MHz and with the re-analysis of archival VLA data at 1.4 GHz. The two components of the radio halo show different observational properties. The central brightest part has higher surface brightess and a spectrum whose steepness is similar to those of the known radio halos, i.e. $\alpha^{\rm 1.78~GHz}_{\rm 118~MHz}=1.33\pm 0.08$. The ridge, which fades into the larger scale emission, is broader in size and has considerably lower surface brightess and a moderately steeper spectrum, i.e. $\alpha^{\rm 1.78~GHz}_{\rm 118~MHz}\sim 1.5$. We propose that the brightest part of the radio halo is powered by the central sloshing in A2142, similar to what has been suggested for mini-halos, or by secondary electrons generated by hadronic collisions in the ICM. On the other hand, the steeper ridge may probe particle re-acceleration by turbulence generated either by stirring the gas and magnetic fields on a larger scale or by less energetic mechanisms, such as continuous infall of galaxy groups or an off-axis merger.
Observations of the Milky Way (MW), M31, and their vicinity, known as the Local Group (LG), can provide clues about the sources of reionization. We present a suite of radiative transfer simulations based on initial conditions provided by the Constrained Local UniversE Simulations (CLUES) project that are designed to recreate the Local Universe, including a realistic MW-M31 pair and a nearby Virgo. Our box size (91 Mpc) is large enough to incorporate the relevant sources of ionizing photons for the LG. We employ a range of source models, mimicking the potential effects of radiative feedback for dark matter haloes between $10^{8}-10^{9}$ M$_{\odot}$. Although the LG mostly reionizes in an inside-out fashion, the final 40 per cent of its ionization shows some outside influence. For the LG satellites, we find no evidence that their redshift of reionization is related to the present-day mass of the satellite or the distance from the central galaxy. We find that less than 20 per cent of present-day satellites for MW and M31 have undergone any star formation prior to the end of global reionization. Approximately five per cent of these satellites could be classified as fossils, meaning the majority of star formation occurred at these early times. The more massive satellites have more cumulative star formation prior to the end of global reionization, but the scatter is significant, especially at the low-mass end. Present-day mass and distance from the central galaxy are poor predictors for the presence of ancient stellar populations in satellite galaxies.
A unique feature of gravity is its ability to control the information accessible to any specific observer. We quantify the notion of cosmic information ('CosmIn') for an eternal observer in the universe. Demanding the finiteness of CosmIn requires the universe to have a late-time accelerated expansion. Combining the introduction of CosmIn with generic features of the quantum structure of spacetime (e.g., the holographic principle), we present a holistic model for cosmology. We show that (i) the numerical value of the cosmological constant, as well as (ii) the amplitude of the primordial, scale invariant, perturbation spectrum can be determined in terms of a single free parameter, which specifies the energy scale at which the universe makes a transition from a pre-geometric phase to the classical phase. For a specific value of the parameter, we obtain the correct results for both (i) and (ii). This formalism also shows that the quantum gravitational information content of spacetime can be tested using precision cosmology.
We carry out a systematic investigation of the total mass density profile of massive (Mstar>2e11 Msun) early-type galaxies and its dependence on galactic properties and host halo mass with the aid of a variety of lensing/dynamical data and large mock galaxy catalogs. The latter are produced via semi-empirical models that, by design, are based on just a few basic input assumptions. Galaxies, with measured stellar masses, effective radii and S\'{e}rsic indices, are assigned, via abundance matching relations, host dark matter halos characterized by a typical LCDM profile. Our main results are as follows: (i) In line with observational evidence, our semi-empirical models naturally predict that the total, mass-weighted density slope at the effective radius gamma' is not universal, steepening for more compact and/or massive galaxies, but flattening with increasing host halo mass. (ii) Models characterized by a Salpeter or variable initial mass function and uncontracted dark matter profiles are in good agreement with the data, while a Chabrier initial mass function and/or adiabatic contractions/expansions of the dark matter halos are highly disfavored. (iii) Currently available data on the mass density profiles of very massive galaxies (Mstar>1e12 Msun), with Mhalo>3e14 Msun, favor instead models with a stellar profile flatter than a S\'{e}rsic one in the very inner regions (r<3-5 kpc), and a cored NFW or Einasto dark matter profile with median halo concentration a factor of ~2 or <1.3, respectively, higher than those typically predicted by N-body numerical simulations.
We extend the idea of mimetic gravity to a Randall-Sundrum II braneworld model. As for the 4-dimensional mimetic gravity, we isolate the conformal degree of freedom of 5-dimensional gravity in a covariant manner. We assume the bulk metric to be made up of a non-dynamical scalar field $\Phi$ and an auxiliary metric $\tilde{{\cal{G}}}_{AB}$ so that ${\cal{G}}_{AB}= \tilde{{\cal{G}}}^{CD}\,\Phi_{,C}\,\Phi_{,D}\,\tilde{{\cal{G}}}_{AB}$ where $A, B, ...$ are the bulk spacetime indices. Then we show that the induced conformal degree of freedom on the brane as an induced scalar field, plays the role of a mimetic field on the brane. In fact, we suppose that the scalar degree of freedom which mimics the dark sectors on the brane has its origin on the bulk scalar field, $\Phi$. By adopting some suitable mimetic potentials on the brane, we show that this brane mimetic field explains the late time cosmic expansion in the favor of observational data: the equation of state parameter of this field crosses the cosmological constant line in near past from quintessence to phantom phase in a redshift well in the range of observation. We show also that this induced mimetic scalar field has the capability to explain initial time cosmological inflation. We study parameter space of the models numerically in order to constraint the models with Planck2015 data set.
We proposed the generalized holographic dark energy model where infrared cutoff is identified with the combination of the FRW universe parameters: the Hubble rate, particle and future horizons, cosmological constant, the universe life-time (if finite) and their derivatives. It is demonstrated that with the corresponding choice of the cutoff one can map such holographic dark energy to modified gravity or gravity with general fluid. Explicitly, F(R) gravity and general perfect fluid are worked out in detail and corresponding infrared cutoff is found. Using this correspondence, we get realistic inflation or viable dark energy or unified inflationary-dark energy universe in terms of covariant holographic dark energy.
We study an inverse seesaw model of neutrino mass within the framework of $S_4$ flavour symmetry from the requirement of generating non-zero reactor mixing angle $\theta_{13}$ along with correct dark matter relic abundance. The leading order $S_4$ model gives rise to tri-bimaximal type leptonic mixing resulting in $\theta_{13}=0$. Non-zero $\theta_{13}$ is generated at one loop level by extending the model with additional scalar and fermion fields which take part in the loop correction. The particles going inside the loop are odd under an in-built $Z^{\text{Dark}}_2$ symmetry such that the lightest $Z^{\text{Dark}}_2$ odd particle can be a dark matter candidate. Correct neutrino and dark matter phenomenology can be achieved for such one loop corrections either to the light neutrino mass matrix or to the charged lepton mass matrix although the latter case is found to be more predictive. The predictions for neutrinoless double beta decay is also discussed and inverted hierarchy in the charged lepton correction case is found to be disfavoured by the latest KamLAND-Zen data.
Links to: arXiv, form interface, find, astro-ph, recent, 1703, contact, help (Access key information)
The unbiased selection process in surveys of the Sunyaev Zel'dovich effect can unveil new populations of galaxy clusters. We performed a weak lensing analysis of the PSZ2LenS sample, i.e. the PSZ2 galaxy clusters detected by the Planck mission in the sky portion covered by the lensing surveys CFHTLenS and RCSLenS. PSZ2LenS is a statistically complete and homogeneous subsample of the PSZ2 catalogue. The Planck selected clusters appear to be unbiased tracers of the massive end of the cosmological haloes. The mass concentration relation of the sample is in excellent agreement with predictions from the Lambda cold dark matter model. The stacked lensing signal is detected at $14\sigma$ significance over the radial range 0.1<R<3.2 Mpc/h, and is well described by the cuspy dark halo models predicted by numerical simulations. We confirmed that Planck estimated masses are biased low by ~27+-11 per cent with respect to weak lensing masses.
In this paper, extending past works of Del Popolo, we show how a high precision mass function (MF) can be obtained using the excursion set approach and an improved barrier taking implicitly into account a non-zero cosmological constant, the angular momentum acquired by tidal interaction of proto-structures and dynamical friction. In the case of the $\Lambda$CDM paradigm, we find that our MF is in agreement at the 3\% level to Klypin's Bolshoi simulation, in the mass range $M_{\rm vir} = 5 \times 10^9 h^{-1} M_{\odot} -- 5 \times 10^{14} h^{-1} M_{\odot}$ and redshift range $0 \lesssim z \lesssim 10$. For $z=0$ we also compared our MF to several fitting formulae, and found in particular agreement with Bhattacharya's within 3\% in the mass range $10^{12}-10^{16} h^{-1} M_{\odot}$. Moreover, we discuss our MF validity for different cosmologies.
We compare the statistics of parity even and odd multipoles of the cosmic microwave background (CMB) sky from PLANCK full mission temperature measurements. An excess power in odd multipoles compared to even multipoles has previously been found on large angular scales. Motivated by this apparent parity asymmetry, we evaluate directional statistics associated with even compared to odd multipoles, along with their significances. Primary tools are the \emph{Power Tensor} and \emph{Alignment Tensor} statistics. We limit our analysis to the first sixty multipoles i.e., $l=[2,61]$. We find no evidence for statistically unusual alignments of even parity multipoles. More than one independent statistic finds evidence for alignments of anisotropy axes of odd multipoles, with a significance equivalent to $\sim 2 \sigma$ or more. The robustness of alignment axes is tested by making galactic cuts and varying the multipole range. Very interestingly, the region spanned by the (a)symmetry axes is found to broadly contain other parity (a)symmetry axes previously observed in the literature.
Stellar nucleosynthesis proceeds via the deuteron (D), but only a small change in the fundamental constants of nature is required to unbind it. Here, we investigate the effect of altering the binding energy of the deuteron on proton burning in stars. We find that the most definitive boundary in parameter space that divides probably life-permitting universes from probably life-prohibiting ones is between a bound and unbound deuteron. Due to neutrino losses, a ball of gas will undergo rapid cooling or stabilization by electron degeneracy pressure before it can form a stable, nuclear reaction-sustaining star. We also consider a less-bound deuteron, which changes the energetics of the $pp$ and $pep$ reactions. The transition to endothermic pp and pep reactions, and the resulting beta-decay instability of the deuteron, do not seem to represent catastrophic problems for life.
We construct a toy a model which demonstrates that large field single scalar inflation can produce an arbitrarily small scalar to tensor ratio in the window of e-foldings recoverable from CMB experiments. This is done by generalizing the $\alpha$-attractor models to allow the potential to approach a constant as rapidly as we desire for super-planckian field values. This implies that a non-detection of r alone can never rule out entirely the theory of large field inflation.
We perform a forecast analysis of the ability of the LISA space-based interferometer to reconstruct the dark sector interaction using gravitational wave (GW) standard sirens at high redshift. We employ Gaussian process methods to reconstruct the distance-redshift relation in a model independent way. We adopt simulated catalogues of standard sirens given by merging massive black hole binaries (MBHBs) visible by LISA, with an electromagnetic (EM) counterpart detectable by future telescopes. The catalogues are constructed considering three different astrophysical scenarios for the evolution of MBHB mergers based on the semi-analytic model of E. Barausse (2012). We first use these standard siren datasets to assess the potential of LISA in reconstructing a possible interaction between vacuum dark energy and dark matter. Then we combine the LISA cosmological data with supernovae data simulated for the Dark Energy Survey (DES). We consider two scenarios distinguished by the time duration of the LISA mission: 5 and 10 years. Using only LISA standard siren data, the dark sector interaction can be well reconstructed from redshift $z\sim1$ to $z\sim3$ (5 yr) and $z\sim1$ up to $z\sim5$ (10 yr), though the reconstruction is inefficient at lower redshift. When combined with the DES datasets, the interaction is well reconstructed in the whole redshift region from $z\sim0$ to $z\sim3$ (5 yr) and $z\sim0$ to $z\sim5$ (10 yr). MBHB standard sirens can thus be used to constrain the dark sector interaction at redshift ranges not reachable by usual supernovae datasets which probe only the $z\lesssim 1.5$ range. GW standard sirens will not only constitute a complementary and alternative way, with respect to familiar EM observations, to probe the cosmic expansion, but will also provide new tests to constrain possible deviations from the standard $\Lambda$CDM dynamics, especially at high redshift.
A framework is developed enabling the global analysis of the stability of cosmological models using the local geometric characteristics of the infinite-dimensional superspace, i.e. using the generalised Jacobi equation reformulated for pseudo-Riemannian manifolds. We give a direct formalism for dynamical analysis in the superspace, the requisite equation pertinent for stability analysis of the universe by means of generalized covariant and Fermi derivative is derived. Then, the relevant definitions and formulae are retrieved for cosmological models with a scalar field.
Based on thermodynamics, we study the galactic clustering of an expanding Universe by considering the logarithmic and volume (quantum) corrections to Newton's law along with the repulsive effect of a harmonic force induced by the cosmological constant ($\Lambda$) in the formation of the large scale structure of the Universe. We derive the $N$-body partition function for extended-mass galaxies (galaxies with halos) analytically. For this partition function, we compute the exact equations of states, which exhibit the logarithmic, volume and cosmological constant corrections. In this setting, a modified correlation (clustering) parameter (due to these corrections) emerges naturally from the exact equations of state. We compute a corrected grand canonical distribution function for this system. Furthermore, we obtain a deviation in differential forms of the two-point correlation functions for both the point-mass and extended-mass cases. The consequences of these deviations on the correlation function's power law are also discussed.
Supermassive black holes are found at the centre of massive galaxies. During the growth of these black holes they light up to become visible as active galactic nuclei (AGN) and release extraordinary amounts of energy across the electromagnetic spectrum. This energy is widely believed to regulate the rate of star formation in the black holes' host galaxies via so-called "AGN feedback". However, the details of how and when this occurs remains uncertain from both an observational and theoretical perspective. I review some of the observational results and discuss possible observational signatures of the impact of super-massive black hole growth on star formation.
We present an updated halo-dependent and halo-independent analysis of viable light WIMP dark matter candidates which could account for the excess observed in CDMS-II-Si. We include recent constraints from LUX, PandaX-II, and PICO-60, as well as projected sensitivities for XENON1T, SuperCDMS SNOLAB, LZ, DARWIN, DarkSide-20k, and PICO-250, on candidates with spin-independent isospin conserving and isospin-violating interactions, and either elastic or exothermic scattering. We show that there exist dark matter candidates which can explain the CDMS-II-Si data and remain very marginally consistent with the null results of all current experiments, however such models are highly tuned, making a dark matter interpretation of CDMS-II-Si very unlikely. We find that these models can only be ruled out in the future by an experiment comparable to LZ or PICO-250.
A semi-classical analysis of backreaction in an expanding Universe with a conformally coupled scalar field and vacuum energy is presented. It is shown that a local observer perceives de Sitter space to contain a constant thermal energy density despite the dilution from expansion due to a continuous flux of energy radiated from the horizon. The self-consistent solution for the Hubble rate is found to be gradually evolving and at late times deviates significantly from de Sitter. Our results hence imply de Sitter space to be unstable. The solution also suggests that vacuum energy is not strictly constant: the energy flux from the horizon is balanced by the generation of negative vacuum energy, which slowly decreases the overall contribution. The concept of dynamical vacuum energy furthermore leads to a complete thermodynamic interpretation of de Sitter space which provides a simple alternate derivation of the results. We conclude that the fate of a Universe with a cosmological constant is not eternal expansion, but flat spacetime.
The two messenger results of the Galactic center gamma-ray excess and the recent analysis of the antiproton flux from the AMS-02 observations suggest the signals may be owing to the same origin of the dark matter (DM) annihilating into $b \bar b$, while these results seem in tension with the dwarf spheroidal galaxy observations. To give a compatible explanation about it, we consider the pseudoscalar DM particles $S_d^+ S_d^-$ annihilating via $S_d^+ S_d^- \rightarrow S_d^0 S_d^0$, with the process mediated by a new scalar $\phi$ and $S_d^0$ quickly decaying into $b \bar{b}$. The particles $S_d^+$, $S_d^-$ and $S_d^0$ are in triplet in hidden sector with degenerate masses, and the annihilation cross section of DM today is linear dependent on the relative velocity $v_r$. Thus, constraints from the dwarf spheroidal galaxies are relaxed. Considering the constraints of the DM relic density, the thermal equilibrium in the early Universe and the search at collider, limits on the corresponding parameters are derived. Though traces from the new sector seem challenging to be disclosed at collider and in DM direct detections, the indirect search of the gamma ray line from $S_d^0$'s decay has the potential to shed light on DM annihilations, with the energy of the gamma-ray line $\sim m_{S_d^0} /2$, i.e. about 50$-$75 GeV.
We extend the idea of conformal attractors in inflation to non-canonical sectors by developing a non-canonical conformally invariant theory from two different approaches. In the first approach, namely, ${\cal N}=1$ supergravity, the construction is more or less phenomenological, where the non-canonical kinetic sector is derived from a particular form of the K$\ddot{a}$hler potential respecting shift symmetry. In the second approach i.e., superconformal theory, we derive the form of the potential from a superconformal action and it turns out to be exactly of the same form as in the first approach. Conformal breaking of these theories results in a new class of non-canonical models which can govern inflation with modulated shape of the T-models. We further employ this framework to explore inflationary phenomenology with a representative example and show how the form of the K$\ddot{a}$hler potential can possibly be constrained in non-canonical models using the latest confidence contour in the $n_s-r$ plane given by Planck.
Large-scale extragalactic magnetic fields may induce conversions between very-high-energy photons and axion-like particles (ALPs), thereby shielding the photons from absorption on the extragalactic background light. However, in simplified "cell" models, used so far to represent extragalactic magnetic fields, this mechanism would be strongly suppressed by current astrophysical bounds. Here we consider realistic models of extragalactic magnetic fields obtained from large-scale cosmological simulations. Such simulated magnetic fields would have large enhancement in the filaments of matter. As a result, photon-ALP conversions would produce a significant spectral hardening for cosmic TeV photons. This effect would be probed with the upcoming Cherenkov Telescope Array detector. This possible detection would give a unique chance to perform a tomography of the magnetised cosmic web with ALPs.
Links to: arXiv, form interface, find, astro-ph, recent, 1703, contact, help (Access key information)
We study the properties of dark matter haloes in a wide range of modified gravity models, namely, $f(R)$, DGP, and interacting dark energy models. We study the effects of modified gravity and dark energy on the internal properties of haloes, such as the spin and the structural parameters. We find that $f(R)$ gravity enhance the median value of the Bullock spin parameter, but could not detect such effects for DGP and coupled dark energy. $f(R)$ also yields a lower median sphericity and oblateness, while coupled dark energy has the opposite effect. However, these effects are very small. We then study the interaction rate of haloes in different gravity, and find that only strongly coupled dark energy models enhance the interaction rate. We then quantify the enhancement of the alignment of the spins of interacting halo pairs by modified gravity. Finally, we study the alignment of the major axes of haloes with the large-scale structures. The alignment of the spins of interacting pairs of haloes in DGP and coupled dark energy models show no discrepancy with GR, while $f(R)$ shows a weaker alignment. Strongly coupled dark energy shows a stronger alignment of the halo shape with the large-scale structures.
A striking signal of dark matter beyond the standard model is the existence of cores in the centre of galaxy clusters. Recent simulations predict that a Brightest Cluster Galaxy (BCG) inside a cored galaxy cluster will exhibit residual wobbling due to previous major mergers, long after the relaxation of the overall cluster. This phenomena is absent with standard cold dark matter where a cuspy density profile keeps a BCG tightly bound at the centre. We test this hypothesis using cosmological simulations and deep observations of 10 galaxy clusters acting as strong gravitational lenses. Modelling the BCG wobble as a simple harmonic oscillator, we measure the wobble amplitude, A_w, in the BAHAMAS suite of cosmological hydrodynamical simulations, finding an upper limit for the CDM paradigm of $A_w < 2$ kpc at the 95% confidence limit. We carry out the same test on the data finding a non-zero amplitude of $A_w = 11.82^{+7.3}_{-3.0}$~kpc with the observations dis-favouring $A_w = 0$ at the $3\sigma$ confidence level. This detection of BCG wobbling is evidence for a dark matter core at the heart of galaxy clusters. It also shows that strong lensing models of clusters cannot assume that the BCG is exactly coincident with the large scale halo. While our small sample of galaxy clusters already indicates a non-zero Aw, with larger surveys, e.g. Euclid, we will be able to not only to confirm the effect but also to use it to determine whether or not the wobbling finds its origin in new fundamental physics or astrophysical process.
I study the process of dark matter capture by the Sun, under the assumption of a Weakly Interacting Massive Particle (WIMP), in the framework of non-relativistic effective field theory. Hypothetically, WIMPs from the galactic halo can scatter against atomic nuclei in the solar interior, settle to thermal equilibrium with the solar core and annihilate to produce an observable flux of neutrinos. In particular, I examine the thermalization process using Monte-Carlo integration of WIMP trajectories. I consider WIMPs in a mass range of 10--1000 GeV and WIMP-nucleon interaction operators with different dependence on spin and transferred momentum. I find that the density profiles of captured WIMPs are in accordance with a thermal profile described by the Sun's gravitational potential and core temperature. Depending on the operator that governs the interaction, the majority of the thermalization time is spent in either the solar interior or exterior. If normalizing the WIMP-nuclei interaction strength to a specific capture rate, I find that the thermalization time differs at most by 3 orders of magnitude between operators. In most cases of interest, the thermalization time is many orders of magnitude shorter than the age of the solar system.
Axion-like particles (ALPs) and photons can quantum mechanically interconvert when propagating through magnetic fields, and ALP-photon conversion may induce oscillatory features in the spectra of astrophysical sources. We use deep (370 ks), short frame time Chandra observations of the bright nucleus at the centre of the radio galaxy M87 in the Virgo cluster to search for signatures of light ALPs. The absence of substantial irregularities in the X-ray power-law spectrum leads to a new upper limit on the photon-ALP coupling, $g_{a\gamma}$: using a conservative model of the cluster magnetic field consistent with Faraday rotation measurements from M87 and M84, we find $g_{a \gamma} < 1.5\times10^{-12}$ GeV$^{-1}$ at 95% confidence level for ALP masses $m_a \leq 10^{-13}$ eV. This constraint is a factor of $\gtrsim 3$ stronger than the bound inferred from the absence of a gamma-ray burst from SN1987A, and it rules out a substantial fraction of the parameter space accessible to future experiments such as ALPS-II and IAXO.
This is the introductory Chapter in the monograph Loop Quantum Gravity: The First 30 Years, edited by the authors, that was just published in the series "100 Years of General Relativity. The 8 invited Chapters that follow provide fresh perspectives on the current status of the field from some of the younger and most active leaders who are currently shaping its development. The purpose of this Chapter is to provide a global overview by bridging the material covered in subsequent Chapters. The goal and scope of the monograph is described in the Preface which can be read by following the Front Matter link at the website listed below.
Recently, the author proposed an alternative vector theory of gravity. To the best of our knowledge, vector gravity also passes available tests of gravity, and, in addition, predicts the correct value of the cosmological constant without free parameters. It is important to find a new feasible test which can distinguish between vector gravity and general relativity and determine whether gravity has a vector or a tensor origin. Here we propose such an experiment based on measurement of propagation direction of gravitational waves relative to the perpendicular arms of a laser interferometer. We show that transverse gravitational wave in vector gravity produces no signal when it propagates in the direction perpendicular to the interferometer plane or along one of the arms. In contrast, general relativistic wave yields no signal when it propagates parallel to the interferometer plane at $45^{\circ }$ angle relative to an interferometer arm. The test can be performed in the nearest years in a joint run of the two LIGO and one Virgo interferometers.
Supermassive primordial stars are now suspected to be the progenitors of the most massive quasars at z~6. Previous studies of such stars were either unable to resolve hydrodynamical timescales or considered stars in isolation, not in the extreme accretion flows in which they actually form. Therefore, they could not self-consistently predict their final masses at collapse, or those of the resulting supermassive black hole seeds, but rather invoked comparison to simple polytropic models. Here, we systematically examine the birth, evolution and collapse of accreting supermassive stars under accretion rates of 0.01-10 solar masses per year using the stellar evolution code KEPLER. KEPLER includes post-Newtonian corrections to the stellar structure and an adaptive nuclear network, and is capable of transitioning to following the hydrodynamic evolution of supermassive stars after they encounter the general relativistic instability. We find that this instability triggers the collapse of the star at 150,000-330,000 solar masses for accretion rates of 0.1-10 solar masses per year, and that the final mass of the star scales roughly logarithmically with the rate. Collapse is sensitive to the mass of the convective core during accretion, so any departures from the treatment of convection, the heat content of the outer accreted envelope, or the boundary conditions in our study may lead to deviations in the final mass of the star that worsen with accretion rate. Since these stars collapse directly to black holes, our models place an upper limit of ~300,000 solar masses on the masses of the first quasars at birth.
Links to: arXiv, form interface, find, astro-ph, recent, 1703, contact, help (Access key information)
We study the properties of dark matter haloes in a wide range of modified gravity models, namely, $f(R)$, DGP, and interacting dark energy models. We study the effects of modified gravity and dark energy on the internal properties of haloes, such as the spin and the structural parameters. We find that $f(R)$ gravity enhance the median value of the Bullock spin parameter, but could not detect such effects for DGP and coupled dark energy. $f(R)$ also yields a lower median sphericity and oblateness, while coupled dark energy has the opposite effect. However, these effects are very small. We then study the interaction rate of haloes in different gravity, and find that only strongly coupled dark energy models enhance the interaction rate. We then quantify the enhancement of the alignment of the spins of interacting halo pairs by modified gravity. Finally, we study the alignment of the major axes of haloes with the large-scale structures. The alignment of the spins of interacting pairs of haloes in DGP and coupled dark energy models show no discrepancy with GR, while $f(R)$ shows a weaker alignment. Strongly coupled dark energy shows a stronger alignment of the halo shape with the large-scale structures.
A striking signal of dark matter beyond the standard model is the existence of cores in the centre of galaxy clusters. Recent simulations predict that a Brightest Cluster Galaxy (BCG) inside a cored galaxy cluster will exhibit residual wobbling due to previous major mergers, long after the relaxation of the overall cluster. This phenomena is absent with standard cold dark matter where a cuspy density profile keeps a BCG tightly bound at the centre. We test this hypothesis using cosmological simulations and deep observations of 10 galaxy clusters acting as strong gravitational lenses. Modelling the BCG wobble as a simple harmonic oscillator, we measure the wobble amplitude, A_w, in the BAHAMAS suite of cosmological hydrodynamical simulations, finding an upper limit for the CDM paradigm of $A_w < 2$ kpc at the 95% confidence limit. We carry out the same test on the data finding a non-zero amplitude of $A_w = 11.82^{+7.3}_{-3.0}$~kpc with the observations dis-favouring $A_w = 0$ at the $3\sigma$ confidence level. This detection of BCG wobbling is evidence for a dark matter core at the heart of galaxy clusters. It also shows that strong lensing models of clusters cannot assume that the BCG is exactly coincident with the large scale halo. While our small sample of galaxy clusters already indicates a non-zero Aw, with larger surveys, e.g. Euclid, we will be able to not only to confirm the effect but also to use it to determine whether or not the wobbling finds its origin in new fundamental physics or astrophysical process.
I study the process of dark matter capture by the Sun, under the assumption of a Weakly Interacting Massive Particle (WIMP), in the framework of non-relativistic effective field theory. Hypothetically, WIMPs from the galactic halo can scatter against atomic nuclei in the solar interior, settle to thermal equilibrium with the solar core and annihilate to produce an observable flux of neutrinos. In particular, I examine the thermalization process using Monte-Carlo integration of WIMP trajectories. I consider WIMPs in a mass range of 10--1000 GeV and WIMP-nucleon interaction operators with different dependence on spin and transferred momentum. I find that the density profiles of captured WIMPs are in accordance with a thermal profile described by the Sun's gravitational potential and core temperature. Depending on the operator that governs the interaction, the majority of the thermalization time is spent in either the solar interior or exterior. If normalizing the WIMP-nuclei interaction strength to a specific capture rate, I find that the thermalization time differs at most by 3 orders of magnitude between operators. In most cases of interest, the thermalization time is many orders of magnitude shorter than the age of the solar system.
Axion-like particles (ALPs) and photons can quantum mechanically interconvert when propagating through magnetic fields, and ALP-photon conversion may induce oscillatory features in the spectra of astrophysical sources. We use deep (370 ks), short frame time Chandra observations of the bright nucleus at the centre of the radio galaxy M87 in the Virgo cluster to search for signatures of light ALPs. The absence of substantial irregularities in the X-ray power-law spectrum leads to a new upper limit on the photon-ALP coupling, $g_{a\gamma}$: using a conservative model of the cluster magnetic field consistent with Faraday rotation measurements from M87 and M84, we find $g_{a \gamma} < 1.5\times10^{-12}$ GeV$^{-1}$ at 95% confidence level for ALP masses $m_a \leq 10^{-13}$ eV. This constraint is a factor of $\gtrsim 3$ stronger than the bound inferred from the absence of a gamma-ray burst from SN1987A, and it rules out a substantial fraction of the parameter space accessible to future experiments such as ALPS-II and IAXO.
This is the introductory Chapter in the monograph Loop Quantum Gravity: The First 30 Years, edited by the authors, that was just published in the series "100 Years of General Relativity. The 8 invited Chapters that follow provide fresh perspectives on the current status of the field from some of the younger and most active leaders who are currently shaping its development. The purpose of this Chapter is to provide a global overview by bridging the material covered in subsequent Chapters. The goal and scope of the monograph is described in the Preface which can be read by following the Front Matter link at the website listed below.
Recently, the author proposed an alternative vector theory of gravity. To the best of our knowledge, vector gravity also passes available tests of gravity, and, in addition, predicts the correct value of the cosmological constant without free parameters. It is important to find a new feasible test which can distinguish between vector gravity and general relativity and determine whether gravity has a vector or a tensor origin. Here we propose such an experiment based on measurement of propagation direction of gravitational waves relative to the perpendicular arms of a laser interferometer. We show that transverse gravitational wave in vector gravity produces no signal when it propagates in the direction perpendicular to the interferometer plane or along one of the arms. In contrast, general relativistic wave yields no signal when it propagates parallel to the interferometer plane at $45^{\circ }$ angle relative to an interferometer arm. The test can be performed in the nearest years in a joint run of the two LIGO and one Virgo interferometers.
Supermassive primordial stars are now suspected to be the progenitors of the most massive quasars at z~6. Previous studies of such stars were either unable to resolve hydrodynamical timescales or considered stars in isolation, not in the extreme accretion flows in which they actually form. Therefore, they could not self-consistently predict their final masses at collapse, or those of the resulting supermassive black hole seeds, but rather invoked comparison to simple polytropic models. Here, we systematically examine the birth, evolution and collapse of accreting supermassive stars under accretion rates of 0.01-10 solar masses per year using the stellar evolution code KEPLER. KEPLER includes post-Newtonian corrections to the stellar structure and an adaptive nuclear network, and is capable of transitioning to following the hydrodynamic evolution of supermassive stars after they encounter the general relativistic instability. We find that this instability triggers the collapse of the star at 150,000-330,000 solar masses for accretion rates of 0.1-10 solar masses per year, and that the final mass of the star scales roughly logarithmically with the rate. Collapse is sensitive to the mass of the convective core during accretion, so any departures from the treatment of convection, the heat content of the outer accreted envelope, or the boundary conditions in our study may lead to deviations in the final mass of the star that worsen with accretion rate. Since these stars collapse directly to black holes, our models place an upper limit of ~300,000 solar masses on the masses of the first quasars at birth.
Links to: arXiv, form interface, find, astro-ph, recent, 1703, contact, help (Access key information)