We have developed a parallel code called "VoIgt profile Parameter Estimation Routine (VIPER)" for automatically fitting the H I Ly-$\alpha$ forest seen in the spectra of QSOs. We obtained the H I column density distribution function (CDDF) and line width ($b$) parameter distribution for $z < 0.45$ using spectra of 82 QSOs obtained using Cosmic Origins Spectrograph and VIPER. Consistency of these with the existing measurements in the literature validate our code. By comparing this CDDF with those obtained from hydrodynamical simulation, we constrain the H I photoionization rate ($\Gamma_{\rm HI}$) at $z < 0.45$ in four redshift bins. The VIPER, together with the Code for Ionization and Temperature Evolution (CITE) we have developed for GADGET-2, allows us to explore parameter space and perform $\chi^2$ minimization to obtain $\Gamma_{\rm HI}$. We notice that the $b$ parameters from the simulations are smaller than what are derived from the observations. We show the observed $b$ parameter distribution and $b$ vs $\log {\rm N_{HI}}$ scatter can be reproduced in simulation by introducing sub-grid scale turbulence. However, it has very little influence on the derived $\Gamma_{\rm HI}$. The $\Gamma_{\rm HI}(z)$ obtained here, $(3.9 \pm 0.1) \times 10^{-14} \; (1+z)^{4.98 \pm 0.11} \;{\rm s^{-1}}$, is in good agreement with those derived by us using flux based statistics in the previous paper. These are consistent with the hydrogen ionizing ultra-violet (UV) background being dominated mainly by QSOs without needing any contribution from the non-standard sources of the UV photons.
One of the fundamental results used in observational cosmology is the distance duality relation (DDR), which relates the luminosity distance, $D_L$, with angular diameter distance, $D_A$, at a given redshift $z$. We suggest to employ the observed limits of this relation to constrain the coupling of axion like particles (ALPs) with photons. With available data we are able to provide improved mixing limit. The method can provide very stringent constraint on ALPs mixing with future improved DDR observations. Also any deviation in DDR can be conventionally explained as photons decaying to axions or vice-versa.
We study the covariance matrix of the power spectrum and bispectrum for dark matter and halos. Using a large suite of simulations from the DEUS-PUR project, we find that the non-Gaussian contributions to the covariance of the power spectrum and bispectrum are significant for both dark matter and halos already at the mildly nonlinear scales. We compute the leading disconnected non-Gaussian correction to the matter bispectrum covariance, and find that the predictions improve the agreement in the mildly nonlinear regime. The shot noise contributions to the covariance of the halo power spectrum and bispectrum are computed. When the ensemble averaged number density is used, the Poisson model covariances are in decent agreement with the measurements. When the number density is estimated and subtracted from each realization, the covariances are significantly reduced and get close to the Gaussian ones. The signal-to-noise ratio, S/N of the halo power spectrum levels off in the mildly nonlinear regime, $k \sim 0.1 - 0.2 \mathrm{Mpc}^{-1} h $. In the nonlinear regime the S/N of the matter and halo bispectrum increases but much slower than the Gaussian results suggest. The S/N for power spectrum and bispectrum are overestimated by the Gaussian covariances, but the problem being much more serious for the bispectrum. While the Gaussian covariance suggests that the S/N of the matter bispectrum should surpass that of the matter power spectrum at $k \sim 0.2 \mathrm{Mpc}^{-1} h $, the full non-Gaussian results show that this only occurs at $k \sim 1 \mathrm{Mpc}^{-1} h $. Because the bispectrum is affected strongly by nonlinearity and shot noise, inclusion of the bispectrum only adds modest amount of S/N. [abridged]
We infer gravitational lensing shear and convergence fields from galaxy ellipticity catalogs under a spatial process prior for the lensing potential. We demonstrate the performance of our algorithm with simulated Gaussian-distributed cosmological lensing shear maps and a reconstruction of the mass distribution of the merging galaxy cluster Abell 781 using galaxy ellipticities measured with the Deep Lens Survey. Given interim posterior samples of lensing shear or convergence fields on the sky, we describe an algorithm to infer cosmological parameters via lens field marginalization. In the most general formulation of our algorithm we make no assumptions about weak shear or Gaussian distributed shape noise or shears. Because we require solutions and matrix determinants of a linear system of dimension that scales with the number of galaxies, we expect our algorithm to require parallel high-performance computing resources for application to ongoing wide field lensing surveys.
21cm intensity mapping is a novel approach aimed at measuring the power spectrum of density fluctuations and deducing cosmological information, notably from the Baryonic Acoustic Oscillations (BAO). We give an update on the progress of BAO from Integrated Neutral Gas Observations (BINGO) which is a single dish intensity mapping project. First we explain the basic ideas behind intensity mapping concept before updating the instrument design for BINGO. We also outline the survey we plan to make and its projected science output including estimates of cosmological parameters.
We use weak-lensing shear measurements to determine the mean mass of optically selected galaxy clusters in Dark Energy Survey Science Verification data. In a blinded analysis, we split the sample of more than 8,000 redMaPPer clusters into 15 subsets, spanning ranges in the richness parameter $5 \leq \lambda \leq 180$ and redshift $0.2 \leq z \leq 0.8$, and fit the averaged mass density contrast profiles with a model that accounts for seven distinct sources of systematic uncertainty: shear measurement and photometric redshift errors; cluster-member contamination; miscentering; deviations from the NFW halo profile; halo triaxiality; and line-of-sight projections. We combine the inferred cluster masses to estimate the joint scaling relation between mass, richness and redshift, $\mathcal{M}(\lambda,z) \varpropto M_0 \lambda^{F} (1+z)^{G}$. We find $M_0 \equiv \langle M_{200\mathrm{m}}\,|\,\lambda=30,z=0.5\rangle=\left[ 2.35 \pm 0.22\ \rm{(stat)} \pm 0.12\ \rm{(sys)} \right] \cdot 10^{14}\ M_\odot$, with $F = 1.12\,\pm\,0.20\ \rm{(stat)}\, \pm\, 0.06\ \rm{(sys)}$ and $G = 0.18\,\pm\, 0.75\ \rm{(stat)}\, \pm\, 0.24\ \rm{(sys)}$. The amplitude of the mass-richness relation is in excellent agreement with the weak-lensing calibration of redMaPPer clusters in SDSS by Simet et al. (2016) and with the Saro et al. (2015) calibration based on abundance matching of SPT-detected clusters. Our results extend the redshift range over which the mass-richness relation of redMaPPer clusters has been calibrated with weak lensing from $z\leq 0.3$ to $z\leq0.8$. Calibration uncertainties of shear measurements and photometric redshift estimates dominate our systematic error budget and require substantial improvements for forthcoming studies.
Due to cosmic variance we cannot learn any more about large-scale inhomogeneities from the primary cosmic microwave background (CMB) alone. More information on large scales is essential for resolving large angular scale anomalies in the CMB. Here we consider cross correlating the large-scale kinetic Sunyaev Zel'dovich (kSZ) effect and probes of large-scale structure, a technique known as kSZ tomography. The statistically anisotropic component of the cross correlation encodes the CMB dipole as seen by free electrons throughout the observable Universe, providing information about long wavelength inhomogeneities. We compute the large angular scale power asymmetry, constructing the appropriate transfer functions, and estimate the cosmic variance limited signal to noise for a variety of redshift bin configurations. The signal to noise is significant over a large range of power multipoles and numbers of bins. We present a simple mode counting argument indicating that kSZ tomography can be used to estimate more modes than the primary CMB on comparable scales. This paper motivates a more systematic investigation of how close to the cosmic variance limit it will be possible to get with future observations.
The idea of the global gravitational effect as the source of cosmological redshift was considered by de Sitter (1916, 1917), Eddington (1923), Tolman (1929) and Bondi (1947). Also Hubble (1929) called the discovered distance-redshift relation as "De Sitter effect". For homogeneous matter distribution cosmological gravitational redshift is proportional to square of distance: z_grav ~ r^2. However for a fractal matter distribution having the fractal dimension D=2 the global gravitational redshift is the linear function of distance: z_grav ~ r, which gives possibility for interpretation of the Hubble law without the space expansion. Here the field gravity fractal cosmological model (FGF) is presented, which based on two initial principles. The first assumption is that the Feynman's field gravity approach describes the gravitational interaction, which delivers a natural basis for the conceptual unity of all fundamental physical interactions within the framework of the relativistic and quantum fields in Minkowski space. The second hypothesis is that the spatial distribution of gravitating matter is a fractal at all scales up to the Hubble radius. The fractal dimension of matter distribution is assumed to be D = 2, which implies that the global gravitational redshift is the explanation of the observed linear Hubble law. In the frame of the FGF all three phenomena - the cosmic background radiation, the fractal large scale structure, and the Hubble law, - could be the consequence of a unique large scale structure evolution process of the initially homogeneous ordinary matter without nonbaryonic matter and dark energy.
We consider the scale-invariant inflationary model studied in [1]. The Lagrangian includes all the scale-invariant operators that can be built with combinations of $R, R^{2}$ and one scalar field. The equations of motion show that the symmetry is spontaneously broken after an arbitrarily long inflationary period and a fundamental mass scale is generated. Upon symmetry breaking, and in the Jordan frame, both Hubble function and the scalar field undergo damped oscillations that can eventually amplify Standard Model fields and reheat the Universe. In the present work, we study in detail inflation and the reheating mechanism of this model in the Einstein frame and we compare some of the results with the latest observational data.
While the Lyman-$\alpha$ ($\mathrm{Ly}\alpha$) emission line serves as an important tool in the study of galaxies at $z\lesssim 6$, finding Ly$\alpha$ emitters (LAE) at significantly higher redshifts has been more challenging, probably because of the increasing neutrality of the intergalactic medium above $z\sim6$. Galaxies with extremely high rest-frame Ly$\alpha$ equivalent widths, EW(Ly$\alpha$) $\gtrsim 150$ \AA{}, at $z>6$ are good candidates for Ly$\alpha$ follow-up observations, and can stand out in multiband imaging surveys because of their unusual colors. We have conducted a photometric search for such objects in the Cluster Lensing And Supernova survey with Hubble (CLASH), and report here the identification of three likely gravitationally-lensed images of a single LAE candidate at $z\sim6.3$, behind the galaxy cluster Abell 2261($z = 0.225$). In the process, we also measured with Keck/MOSFIRE the first spectroscopic redshift of a multiply-imaged galaxy behind Abell 2261, at $z = 3.337$. This allows us to calibrate the lensing model, which in turn is used to study the properties of the candidate LAE. Population III galaxy spectral energy distribution (SED) model fits to the CLASH broadband photometry of the possible LAE provide a slightly better fit than Population I/II models. The best fitted model suggests intrinsic EW(Ly$\alpha$) $\approx 160$ \AA{} after absorption in the interstellar and intergalactic medium. Future spectroscopic observations will examine this prediction as well as shed more light on the morphology of this object, which indicates it may be a merger of two smaller galaxies.
We study the average Ly$\alpha$ emission associated with high-$z$ strong (log $N$(H I) $\ge$ 21) damped Ly$\alpha$ systems (DLAs). We report Ly$\alpha$ luminosities ($L_{\rm Ly\alpha}$) for the full as well as various sub-samples based on $N$(H I), $z$, $(r-i)$ colours of QSOs and rest equivalent width of Si II$\lambda$1526 line (i.e., $W_{1526}$). For the full sample, we find $L_{\rm Ly\alpha}$$< 10^{41} (3\sigma)\ \rm erg\ s^{-1}$ with a $2.8\sigma$ level detection of Ly$\alpha$ emission in the red part of the DLA trough. The $L_{\rm Ly\alpha}$ is found to be higher for systems with higher $W_{1526}$ with its peak, detected at $\geq 3\sigma$, redshifted by about 300-400 $\rm km\ s^{-1}$ with respect to the systemic absorption redshift, as seen in Lyman Break Galaxies (LBGs) and Ly$\alpha$ emitters. A clear signature of a double-hump Ly$\alpha$ profile is seen when we consider $W_{1526} \ge 0.4$ \AA\ and $(r-i) < 0.05$. Based on the known correlation between metallicity and $W_{1526}$, we interpret our results in terms of star formation rate (SFR) being higher in high metallicity (mass) galaxies with high velocity fields that facilitates easy Ly$\alpha$ escape. The measured Ly$\alpha$ surface brightness requires local ionizing radiation that is 4 to 10 times stronger than the metagalactic UV background at these redshifts. The relationship between the SFR and surface mass density of atomic gas seen in DLAs is similar to that of local dwarf and metal poor galaxies. We show that the low luminosity galaxies will contribute appreciably to the stacked spectrum if the size-luminosity relation seen for H I at low-$z$ is also present at high-$z$. Alternatively, large Ly$\alpha$ halos seen around LBGs could also explain our measurements.
The inflationary universe can be viewed as a "Cosmological Collider" with energy of Hubble scale, producing very massive particles and recording their characteristic signals in primordial non-Gaussianities. To utilize this collider to explore any new physics at very high scales, it is a prerequisite to understand the background signals from the particle physics Standard Model. In this paper we describe the Standard Model background of the Cosmological Collider.
We investigate the stability of de Sitter space as seen by a local observer in expanding space. Using the Bunch-Davies vacuum as an initial state we find for a conformal scalar field and classical vacuum energy that tracing over the unobservable states beyond the cosmological horizon leads to a thermal spectrum of particles and that such a configuration is unstable under semi-classical backreaction. It is shown that this instability results in a gradual increase in the horizon size. Comments welcome.
The knowledge of the structure of the magnetic field inside a blazar jet, as
deduced from polarization observations at radio to optical wavelengths, is
closely related to the formation and propagation of relativistic jets that
result from accretion onto supermassive black holes. However, a largely
unexplored aspect of the theoretical understanding of radiation transfer
physics in blazar jets has been the magnetic field geometry as revealed by the
polarized emission and the connection between the variability in polarization
and flux across the spectrum. Here, we explore the effects of various magnetic
geometries that can exist inside a blazar jet: parallel, oblique, toroidal, and
tangled. We investigate the effects of changing the orientation of the magnetic
field, according to the above-mentioned geometries, on the resulting
high-energy spectral energy distributions (SEDs) and spectral variability
patterns (SVPs) of a typical blazar. We use the MUlti-ZOne Radiation Feedback
(MUZORF) model of Joshi et al. (2014) to carry out this study and to relate the
geometry of the field to the observed SEDs at X-ray and gamma-ray energies.
One of the goals of the study is to understand the relationship between
synchrotron and inverse Compton peaks in blazar SEDs and the reason for the
appearance of gamma-ray "orphan flares" observed in some blazars. This can be
associated with the directionality of the magnetic field, which creates a
difference in the radiation field as seen by an observer versus that seen by
the electrons in the emission region.
An all-sky monitor with nuclear imaging spectroscopy is a promising tool for the systematic study of supernova explosions. In particular, progenitor scenarios of type-Ia supernovae, which are not yet well understood, can be resolved using light curves in the nuclear gamma-ray band. Here we report an expected result of an all-sky monitor with imaging spectroscopy using electron-tracking Compton camera, which will enable us to observe nuclear gamma-ray lines from type-Ia supernovae.
Links to: arXiv, form interface, find, astro-ph, recent, 1610, contact, help (Access key information)
Precise measurements of the temperature and polarization anisotropies of the cosmic microwave background can be used to constrain the annihilation and decay of dark matter. In this work, we demonstrate via principal component analysis that the imprint of dark matter decay on the cosmic microwave background can be approximately parameterized by a single number for any given dark matter model. We develop a simple prescription for computing this model-dependent detectability factor, and demonstrate how this approach can be used to set model-independent bounds on a large class of decaying dark matter scenarios. We repeat our analysis for decay lifetimes shorter than the age of the universe, allowing us to set constraints on metastable species other than the dark matter decaying at early times, and decays that only liberate a tiny fraction of the dark matter mass energy. We set precise bounds and validate our principal component analysis using a Markov Chain Monte Carlo approach and Planck 2015 data.
The Detection of redshifted 21 cm emission from the epoch of reionization (EoR) is a challenging task owing to strong foregrounds that dominate the signal. In this paper, we propose a general method, based on the delay spectrum approach, to extract HI power spectra that is applicable to tracking observations using an imaging radio interferometer (Delay Spectrum with Imaging Arrays (DSIA)). Our method is based on modelling the HI signal taking into account the impact of wide field effects such as the $w$-term which are then used as appropriate weights in cross-correlating the measured visibilities. Our method is applicable to any radio interferometer that tracks a phase center and could be utilized for arrays such as MWA, LOFAR, GMRT, PAPER and HERA. In the literature the delay spectrum approach has been implemented for near-redundant baselines using drift scan observations. In this paper we explore the scheme for non-redundant tracking arrays, and this is the first application of delay spectrum methodology to such data to extract the HI signal. We analyze 3 hours of MWA tracking data on the EoR1 field. We present both 2-dimensional ($k_\parallel,k_\perp$) and 1-dimensional (k) power spectra from the analysis. Our results are in agreement with the findings of other pipelines developed to analyse the MWA EoR data.
21 cm intensity mapping has emerged as a promising technique to map the large-scale structure of the Universe, at redshifts $z$ from 1 to 10. Unfortunately, many of the key cross correlations with photo-$z$ galaxies and the CMB have been thought to be impossible due to foreground contamination for radial modes with small wavenumbers. These modes are usually subtracted in the foreground subtraction process. We recover the lost 21 cm radial modes via cosmic tidal reconstruction and find more than 60\% cross-correlation signal at $\ell\lesssim100$ and even more on larger scales can be recovered from null. The tidal reconstruction method opens up a new set of possibilities to probe our Universe and is extremely valuable not only for 21 cm surveys but also for CMB and photometric redshift observations.
We study the nonlinear $E$-mode clustering in Lagrangian space by using large scale structure (LSS) $N$-body simulations and use the displacement field information in Lagrangian space to recover the primordial linear density field. We find that, compared to Eulerian nonlinear density fields, the $E$-mode displacement fields in Lagrangian space improves the cross-correlation scale $k$ with initial density field by factor of 6 $\sim$ 7, containing 2 orders of magnitude more primordial information. This illustrates ability of potential density reconstruction algorithms, to improve the baryonic acoustic oscillation (BAO) measurements from current and future large scale structure surveys.
The turn-around radii of the galaxy groups shows the imprint of a long battle between their self-gravitational forces and the accelerating space. The standard LambdaCDM cosmology based on the general relativity (GR) predicts the existence of an upper bound on the expectation value of the turn-around radius which is rarely violated by individual galaxy groups. We speculate that a deviation of the gravitational law from GR on the cosmological scale could cause an appreciable shift of the mean turn-around radius to higher values and make the occurrence of the bound violation more probable. Analyzing the data from high-resolution N-body simulations for two specific models with modified gravity (MG) and the standard GR+LambdaCDM cosmology, we determine the turn-around radii of the massive Rockstar groups from the peculiar motions of the galactic halos located in the bound zone where the fifth force generated by MG is expected to be at most partially shielded. We detect a 4 sigma signal of difference in the odds of the bound violations between a fiducial MG and the GR models, proving that the odds of the bound violations increase with the strength of the fifth force produced by the presence of MG. The advantage of using the odds of the bound violations as a complementary diagnostics to probe the nature of gravity is discussed.
We test the consistency between Planck temperature and polarization power spectra and the concordance model of $\Lambda$ Cold Dark Matter cosmology ($\Lambda$CDM) within the framework of Crossing statistics. We find that Planck TT best fit $\Lambda$CDM power spectrum is completely consistent with EE power spectrum data while EE best fit $\Lambda$CDM power spectrum is not consistent with TT data. However, this does not point to any systematic or model-data discrepancy since in the Planck EE data, uncertainties are much larger compared to the TT data. We also investigate the possibility of any deviation from $\Lambda$CDM model analyzing the Planck 2015 data. Results from both TT and EE data analysis indicate that no deviation is required beyond the flexibility of the concordance $\Lambda$CDM model. Our analysis thus rules out any strong evidence for beyond the concordance model in the Planck spectra data. We also report a mild amplitude difference comparing temperature and polarization data, where temperature data seems to have slightly lower amplitude than expected (consistently at all multiples), as we assume both temperature and polarization data are realizations of the same underlying cosmology.
Perturbation Theory to Large Scale Structure Cosmology proposes corrections to the linearly evolved density contrast and velocity in terms of a series development in which all terms are integrals of powers of the linear density contrast multiplied by kernels. We discuss the symmetry properties of these kernels and show that their full symmetrized versions can be decomposed in different classes of subkernels. We will construct classes of subkernels with improved symmetry properties, and provide recurrence relations to generate them.
We use the recently released cosmic chronometer data and the latest measured value of the local Hubble parameter, combined with the latest joint light curves of Supernovae Type Ia, and Baryon Acoustic Oscillation distance measurements, in order to impose constraints on the viable and most used $f(R)$ gravity models. We consider four $f(R)$ models, namely the Hu-Sawicki, the Starobinsky, the Tsujikawa, and the exponential one, and we parametrize them introducing a distortion parameter $b$ that quantifies the deviation from $\Lambda$CDM cosmology. Our analysis reveals that a small but non-zero deviation from $\Lambda$CDM cosmology is slightly favored, with the corresponding fittings exhibiting very efficient $AIC$ and $BIC$ Information Criteria values. Clearly, $f(R)$ gravity is consistent with observations, and it can serve as a candidate for modified gravity.
We show that the radial acceleration relation for rotationally-supported galaxies may be explained, in the absence of cold dark matter, by a non-minimally coupled scalar field, whose fifth forces are partially screened on galactic scales by the symmetron mechanism. If realised in nature, this effect could have a significant impact on the inferred density of dark matter halos.
We briefly review the motivation to search for sterile neutrinos in the keV mass scale, as dark matter candidates, and the prospects to find them in beta decay or electron capture spectra, with a global perspective. We describe the fundamentals of the neutrino flavor-mass eigenstate mismatch that opens the possibility of detecting sterile neutrinos in such ordinary nuclear processes. Results are shown and discussed for the effect of heavy neutrino emission in electron capture in Holmium 163 and in two isotopes of Lead, 202 and 205, as well as in the beta decay of Tritium. We study the de-excitation spectrum in the considered cases of electron capture and the charged lepton spectrum in the case of Tritium beta decay. For each of these cases, we define ratios of integrated transition rates over different regions of the spectrum under study, and give new results that may guide and facilitate the analysis of possible future measurements, paying particular attention to forbidden transitions in Lead isotopes.
A hidden sector with pure non-abelian gauge symmetry is an elegant and just about the simplest model of dark matter. In this model the dark matter candidate is the lightest bound state made of the confined gauge fields, the dark glueball. In spite of its simplicity, the model has been shown to have several interesting non-standard implications in cosmology. In this work, we explore the gravitational waves from binary boson stars made of self-gravitating dark glueball fields as a natural and important consequence. We derive the dark SU($N$) star mass and radius as functions of the only two fundamental parameters in the model, the glueball mass $m$ and the number of colors $N$, and identify the regions that could be probed by the LIGO and future gravitational wave observatories.
We propose a new type of Higgs-portal dark matter-production mechanism, called bosonic-seesaw portal scenario. Bosonic seesaw provides the dynamical origin of the electroweak symmetry breaking, triggered by mixing between the elementary Higgs and a composite Higgs generated by a new-color strong dynamics (hypercolor) which dynamically breaks the classical-scale invariance of the model. The composite hypercolor-baryonic matter can then be a dark matter candidate, which significantly couples to the standard-model Higgs via the bosonic seesaw, and can be produced from the thermal plasma below the decoupling temperature around the new strong coupling scale, to account for the observed relic abundance of the dark matter: the dark matter can closely be related to the mechanism of the electroweak symmetry breaking.
We conduct a 350 micron dust continuum emission survey of 17 dust-obscured galaxies (DOGs) at z = 0.05-0.08 with the Caltech Submillimeter Observatory (CSO). We detect 14 DOGs with S_350 = 114-650 mJy and S/N > 3. By including two additional DOGs with submillimeter data in the literature, we are able to study dust contents for a sample of 16 local DOGs that consists of 12 bump and 4 power-law types. We determine their physical parameters with a two-component modified blackbody function model. The derived dust temperatures are in the range 57-122 K and 22-35 K for the warm and cold dust components, respectively. The total dust mass and the mass fraction of warm dust component are 3-34$\times10^{7} M_\odot$ and 0.03-2.52%, respectively. We compare these results with those of other submillimeter-detected infrared luminous galaxies. The bump DOGs, the majority of the DOG sample, show similar distributions of dust temperatures and total dust mass to the comparison sample. The power-law DOGs show a hint of smaller dust masses than other samples, but need to be tested with a larger sample. These findings support that the reason why DOGs show heavy dust obscuration is not an overall amount of dust content, but probably the spatial distribution of dust therein.
The tree method is a widely implemented algorithm for collisionless $N$-body simulations in astrophysics well suited for GPU(s). Adopting hierarchical time stepping can accelerate $N$-body simulations; however, it is infrequently implemented and its potential remains untested in GPU implementations. We have developed a Gravitational Oct-Tree code accelerated by HIerarchical time step Controlling named \texttt{GOTHIC}, which adopts both the tree method and the hierarchical time step. The code adopts some adaptive optimizations by monitoring the execution time of each function on-the-fly and minimizes the time-to-solution by balancing the measured time of multiple functions. Results of performance measurements with realistic particle distribution performed on NVIDIA Tesla M2090, K20X, and GeForce GTX TITAN X, which are representative GPUs of the Fermi, Kepler, and Maxwell generation of GPUs, show that the hierarchical time step achieves a speedup by a factor of around 3--5 times compared to the shared time step. The measured elapsed time per step of \texttt{GOTHIC} is 0.30~s or 0.44~s on GTX TITAN X when the particle distribution represents the Andromeda galaxy or the NFW sphere, respectively, with $2^{24} =$~16,777,216 particles. The averaged performance of the code corresponds to 10--30\% of the theoretical single precision peak performance of the GPU.
The present work is based on a parametric reconstruction of the deceleration parameter $q(z)$ in a model for the spatially flat FRW universe filled with dark energy and non- relativistic matter. We have proposed a divergence-free logarithmic parametrization of $q(z)$ to probe the entire evolution history of the universe. Using the SN Ia and Hubble parameter datasets, the constraints on the arbitrary model parameters $q_{0}$ and $q_{1}$ are obtained (within $1\sigma$ and $2\sigma$ confidence limits) by $\chi^{2}$-minimization technique. We have then reconstructed the deceleration parameter, the total EoS parameter $\omega_{tot}$, the jerk parameter and have compared the reconstructed results with the spatially flat $\Lambda$CDM model. It has been found that the behavior of $q(z)$ and $\omega_{tot}$ in our model are very similar (within $1\sigma$ confidence limit) to that of the $\Lambda$CDM model if we consider Type Ia Supernova (SN Ia) dataset only, but the evolutions of $q(z)$ and $\omega_{tot}$ are different as compare to the $\Lambda$CDM model if we add Hubble parameter dataset with the SN Ia dataset. Interestingly, we have found that the present values of $q(z)$ and $\omega_{tot}$ within $1\sigma$ errors for SN Ia+Hubble dataset are in good agreement with the $\Lambda$CDM model. The best-fit model is also found to be in good agreement with the observational Hubble data and the supernova distance modulus data against redshift parameter.
We consider the existence of Noether symmetries in bigravity cosmologies in order to constrain the material content minimally coupled to the gravitational sector that we are not inhabiting. Interestingly, a Noether symmetry not only constrain the matter content of the universe we do not inhabit but also comes as a sort of bonus on the form of a very interesting dynamics of the universe we live in. In fact, by assuming that our universe is filled with standard matter and radiation, we show that the existence of a Noether symmetry implies the existence of a vacuum energy in our universe that can explain, in a natural way, the current acceleration of the universe. This vacuum energy is intrinsic to the model and can be realized either for a theory that is not properly a bigravity model or for a genuinely bimetric scenario. In fact, it would correspond to a "mono"-universe with a $\Lambda$CDM matter symmetry or to a bimetric world where our universe would have once again a $\Lambda$CDM matter symmetry while the non-inhabited universe could have a gravitational coupling with a different sign to that in ours. The main point in these pictures is that the Noether symmetry plays a major role into dynamics. The physical consequences are also briefly discussed.
Scalar cosmological perturbations in loop quantum cosmology (LQC) is revisited in a covariant manner, using self dual Ashtekar variables. For real-valued Ashtekar-Barbero variables, this `deformed algebra' approach has been shown to implement holonomy corrections from loop quantum gravity (LQG) in a consistent manner, albeit deforming the algebra of modified constraints in the process. This deformation has serious conceptual ramifications, not the least of them being an effective `signature-change' in the deep quantum regime. In this paper, we show that working with self dual variables lead to an undeformed algebra of hypersurface deformations, even after including holonomy corrections in the effective constraints. As a necessary consequence, the diffeomorphism constraint picks up non-perturbative quantum corrections thus hinting at a modification of the underlying space-time structure, a novel ingredient compared to the usual treatment of (spatial) diffeomorphisms in LQG. This work extends a similar result obtained in the context of spherically symmetric gravity coupled to a scalar field, suggesting that self dual variables could be better suited than their real counterparts to treat inhomogeneous LQG models.
Utilizing data from the ELVIS and Via Lactea-II simulations, we characterize the local dark matter subhalo population, and use this information to refine the predictions for the gamma-ray fluxes arising from annihilating dark matter in this class of objects. We find that the shapes of nearby subhalos are significantly altered by tidal effects, and are generally not well described by NFW density profiles, instead prefering power-law profiles with an exponential cutoff. From the subhalo candidates detected by the Fermi Gamma-Ray Space Telescope, we place limits on the dark matter annihilation cross section that are only modestly weaker than those based on observations of dwarf galaxies. We also calculate the fraction of observable subhalos that are predicted to be spatially extended at a level potentially discernible to Fermi.
Links to: arXiv, form interface, find, astro-ph, recent, 1610, contact, help (Access key information)
We present a method that extends the capabilities of the PINpointing Orbit-Crossing Collapsed HIerarchical Objects (PINOCCHIO) code, allowing it to generate accurate dark matter halo mock catalogues in cosmological models where the linear growth factor and the growth rate depend on scale. Such cosmologies comprise, among others, models with massive neutrinos and some classes of modified gravity theories. We validate the code by comparing the halo properties from PINOCCHIO against N-body simulations, focusing on cosmologies with massive neutrinos: $\nu\Lambda$CDM. We analyse the halo mass function, halo two-point correlation function, halo power spectrum and the moments of the halo density field, showing that PINOCCHIO reproduces the results from simulations with the same level of precision as the original code ($\sim5-10\%$). We demonstrate that the abundance of halos in cosmologies with massless and massive neutrinos from PINOCCHIO matches very well the outcome of simulations, and point out that PINOCCHIO can reproduce the $\Omega_\nu-\sigma_8$ degeneracy that affects the halo mass function. We show that the clustering properties of the halos from PINOCCHIO matches accurately those from simulations both in real and redshift-space, in the latter case up to $k=0.3~h~{\rm Mpc}^{-1}$. We finally point out that the first moments of the halo density field from simulations are precisely reproduced by PINOCCHIO. We emphasize that the computational time required by PINOCCHIO to generate mock halo catalogues is orders of magnitude lower than the one needed for N-body simulations. This makes this tool ideal for applications like covariance matrix studies within the standard $\Lambda$CDM model but also in cosmologies with massive neutrinos or some modified gravity theories.
We present \Chandra\ observations of 23 galaxy groups and low-mass galaxy clusters at $0.03<z<0.15$ with a median temperature of ~2keV. The sample is a statistically complete flux-limited subset of the 400 deg$^2$ survey. We investigated the scaling relation between X-ray luminosity ($L$) and temperature ($T$), taking selection biases fully into account. The logarithmic slope of the bolometric \LT\ relation was found to be $3.29\pm0.33$, consistent with values typically found for samples of more massive clusters. In combination with other recent studies of the \LT\ relation we show that there is no evidence for the slope, normalisation, or scatter of the \LT\ relation of galaxy groups being different than that of massive clusters. The exception to this is that in the special case of the most relaxed systems, the slope of the core-excised \LT\ relation appears to steepen from the self-similar value found for massive clusters to a steeper slope for the lower mass sample studied here. Thanks to our rigorous treatment of selection biases, these measurements provide a robust reference against which to compare predictions of models of the impact of feedback on the X-ray properties of galaxy groups.
Recent papers have shown that a small systematic redshift shift ($\Delta z\sim 10^{-5}$) in measurements of type Ia supernovae can cause a significant bias ($\sim$1\%) in the recovery of cosmological parameters. Such a redshift shift could be caused, for example, by a gravitational redshift due to the density of our local environment. The sensitivity of supernova data to redshift shifts means supernovae make excellent probes of inhomogeneities. We therefore invert the analysis, and try to diagnose the nature of our local gravitational environment by fitting for $\Delta z$ as an extra free parameter alongside the usual cosmological parameters.
Recently, the power spectrum (PS) multipoles using the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 12 (DR12) sample are analyzed \cite{160703150}. Even though the based model for the analysis is the so-called TNS quasi-linear model including the multipole up to the eighth order in the window function \cite{TNS}, the analysis provides the multipoles up to the hexadecapole. Thus, one might be able to recover the galaxy PS by using the combination of multipoles to investigate the cosmology \cite{0407214}. We provide the analytic form of this combination of multipoles of the quasi-linear PS including the Fingers of God (FoG) effect to recover the PS at the linear regime. In order to confirm the consistency of the multipole data, we compare the multipole ratios of the linear theory including the FoG effect with those of observation. The data of the ratio of quadrupole to monopole is consistent with that of the linear theory prediction even though the current observational error is too large to distinguish the linear theory prediction from the Kaiser one ({\it i.e. without the FoG effect}). However, the ratio of hexadecapole to monopole of DR12 is inconsistent with the prediction of the linear theory even at large scale. Thus, the current multipoles data is premature to be used for recovering the galaxy PS.
The existence of dark matter and dark energy in cosmology is implied by various observations, however, they are still unclear because they have not been directly detected. In this Letter, an unified model of dark energy and dark matter that can explain the evolution history of the Universe later than inflationary era, the time evolution of the growth rate function of the matter density contrast, the flat rotation curves of the spiral galaxies, and the gravitational experiments in the solar system is proposed in mimetic matter model.
Radio halos are diffuse synchrotron sources on scales of ~1 Mpc that are found in merging clusters of galaxies, and are believed to be powered by electrons re-accelerated by the merger-driven turbulence. We present measurements of extended radio emission on similarly large scales in two clusters of galaxies hosting cool cores: Abell 2390 and Abell 2261. The analysis is based on interferometric imaging with the JVLA, VLA and GMRT. We present detailed radio images of the targets, subtract the compact emission components, and measure the spectral indices for the diffuse components. The radio emission in A2390 extends beyond a known sloshing-like brightness discontinuity, and has a very steep in-band spectral slope at 1.5 GHz that is similar to some known ultra-steep spectrum radio halos. The diffuse signal in A2261 is more extended than in A2390 but has lower luminosity. X-ray morphological indicators, derived from XMM-Newton X-ray data, place these clusters in the category of relaxed or regular systems, although some asymmetric features that can indicate past minor mergers are seen in the X-ray brightness images. If these two Mpc-scale radio sources are categorized as giant radio halos, they question the common assumption of radio halos occurring exclusively in clusters undergoing violent merging activity, in addition to commonly used criteria in distinguishing between radio halos and mini-halos.
We analyse the evolution of cosmological perturbations which leads to the formation of large isolated voids in the Universe. We assume that initial perturbations are spherical and all components of the Universe (radiation, matter and dark energy) are continuous media with perfect fluid energy-momentum tensors, which interact only gravitationally. Equations of the evolution of perturbations for every component in the comoving to cosmological background reference frame are obtained from equations of energy and momentum conservation and Einstein's ones and are integrated numerically. Initial conditions are set at the early stage of evolution in the radiation-dominated epoch, when the scale of perturbation is much larger than the particle horizon. Results show how the profiles of density and velocity of matter and dark energy are formed and how they depend on parameters of dark energy and initial conditions. In particular, it is shown that final matter density and velocity amplitudes change within range $\sim$4-7\% when the value of equation of state parameter of dark energy $w$ vary in the range from -0.8 to -1.2, and change within $\sim$1\% only when the value of effective sound speed of dark energy vary over all allowable range of its values.
Even if the Higgs field does not affect the evolution of the background geometry, its massive inhomogeneities induce large-scale gauge fields whose energy density depends on the slow-roll parameters, on the effective scalar mass and, last but not least, on the dimensionless coupling to the space-time curvature. Since the non-Abelian gauge modes are screened, the non-minimal coupling to gravity predominantly affects the evolution of the hypercharge and electromagnetic fields. While in the case of minimal coupling the obtained constraints are immaterial, as soon as the coupling increases beyond one fourth the produced fields become overcritical. We chart the whole parameter space of this qualitatively new set of bounds. Whenever the limits on the curvature coupling are enforced, the magnetic field may still be partially relevant for large-scale magnetogenesis and exceed $10^{-20}$ G for the benchmark scale of the protogalactic collapse.
Future data from galaxy redshift surveys, combined with high-resolutions maps of the cosmic microwave background, will enable measurements of the pairwise kinematic Sunyaev-Zel'dovich (kSZ) signal with unprecedented statistical significance. This signal probes the matter-velocity correlation function, scaled by the average optical depth ($\tau$) of the galaxy groups and clusters in the sample, and is thus of fundamental importance for cosmology. However, in order to translate pairwise kSZ measurements into cosmological constraints, external constraints on $\tau$ are necessary. In this work, we present a new model for the intra-cluster medium, which takes into account star-formation, feedback, non-thermal pressure, and gas cooling. Our semi-analytic model is computationally efficient and can reproduce results of recent hydrodynamical simulations of galaxy cluster formation. By calibrating the model using recent X-ray measurements of gas density profiles of clusters and $M_{\mathrm{gas}}-M$ relations of groups and clusters, we show that the integrated $\tau$ is constrained to better than 6% (with a confidence level of 95%), which sets the tightest constraints on the external prior on the optical depth of groups and clusters reported in the literature to date. Our strong prior on $\tau$ should significantly improve cosmological constraints derived from future pairwise kSZ measurements and shed new insights into the nature of dark energy, modified gravity, and neutrino mass.
We show the kinematic equivalence between cosmological models driven by Dirac-Born-Infeld fields $\phi$ with constant proper velocity of the brane and exponential potential $V=V_0e^{-B\phi}$ and interactive cosmological systems with Modified Holographic Ricci type fluids as dark energy in flat Friedmann-Robertson-Walker cosmologies.
Obtaining high-sensitivity measurements of degree-scale cosmic microwave background (CMB) polarization is the most direct path to detecting primordial gravitational waves. Robustly recovering any primordial signal from the dominant foreground emission will require high-fidelity observations at multiple frequencies, with excellent control of systematics. We explore the potential for a new platform for CMB observations, the Airlander 10 hybrid air vehicle, to perform this task. We show that the Airlander 10 platform, operating at commercial airline altitudes, is well-suited to mapping frequencies above 220 GHz, which are critical for cleaning CMB maps of dust emission. Optimizing the distribution of detectors across frequencies, we forecast the ability of Airlander 10 to clean foregrounds of varying complexity as a function of altitude, demonstrating its complementarity with both existing (Planck) and ongoing (C-BASS) foreground observations. This novel platform could play a key role in defining our ultimate view of the polarized microwave sky.
Early-type galaxies provide unique tests for the predictions of the cold dark matter cosmology and the baryonic physics assumptions entering models for galaxy formation. In this work, we use the Illustris simulation to study correlations of three main properties of early-type galaxies, namely, the stellar orbital anisotropies, the central dark matter fractions and the central radial density slopes, as well as their redshift evolution since $z=1.0$. We find that lower-mass galaxies or galaxies at higher redshift tend to be bluer in rest-frame colour, have higher central gas fractions, and feature more tangentially anisotropic orbits and steeper central density slopes than their higher-mass or lower-redshift counterparts, respectively. The projected central dark matter fraction within the effective radius shows no significant mass dependence but positively correlates with galaxy effective radii due to the aperture effect. The central density slopes obtained in the simulation by combining strong lensing measurements with single aperture kinematics are found to be shallower than the true density slopes. We identify systematic biases in this measurement due to two common modelling assumptions, isotropic stellar orbital distributions and power-law density profiles. We also compare the properties of early-type galaxies in Illustris to those from the SLACS, SL2S and BOSS surveys, finding in general broad agreement but also some tension, which appears to be mostly caused by too large galaxy sizes in Illustris.
The spectrum of Weakly-Interacting-Massive-Particle (WIMP) dark matter generically possesses bound states when the WIMP mass becomes sufficiently large relative to the mass of the electroweak gauge bosons. The presence of these bound states enhances the annihilation rate via resonances in the Sommerfeld enhancement, but they can also be produced directly with the emission of a low-energy photon. In this work we compute the rate for SU(2) triplet dark matter (the wino) to bind into WIMPonium -- which is possible via single-photon emission for wino masses above 5 TeV for relative velocity v < O(10^{-2}) -- and study the subsequent decays of these bound states. We present results with applications beyond the wino case, e.g. for dark matter inhabiting a nonabelian dark sector; these include analytic capture and transition rates for general dark sectors in the limit of vanishing force carrier mass, efficient numerical routines for calculating positive and negative-energy eigenstates of a Hamiltonian containing interactions with both massive and massless force carriers, and a study of the scaling of bound state formation in the short-range Hulthen potential. In the specific case of the wino, we find that the rate for bound state formation is suppressed relative to direct annihilation, and so provides only a small correction to the overall annihilation rate. The soft photons radiated by the capture process and by bound state transitions could permit measurement of the dark matter's quantum numbers; for wino-like dark matter, such photons are rare, but might be observable by a future ground-based gamma-ray telescope combining large effective area and a low energy threshold.
Herein we analyse late-time (post-plateau; 103 < t < 1229 d) optical spectra of low-redshift (z < 0.016), hydrogen-rich Type IIP supernovae (SNe IIP). Our newly constructed sample contains 91 nebular spectra of 38 SNe IIP, which is the largest dataset of its kind ever analysed in one study, and many of the objects have complementary photometric data. We determined the peak and total luminosity, velocity of the peak, HWHM intensity, and profile shape for many emission lines. Temporal evolution of these values and various flux ratios are studied. We also investigate the correlations between these measurements and photometric observables, such as the peak and plateau absolute magnitudes and the late-time light curve decline rates in various optical bands. The strongest and most robust result we find is that the luminosities of all spectral features (except those of helium) tend to be higher in objects with steeper late-time V-band decline rates. A steep late-time V-band slope likely arises from less efficient trapping of gamma-rays and positrons, which could be caused by multidimensional effects such as clumping of the ejecta or asphericity of the explosion itself. Furthermore, if gamma-rays and positrons can escape more easily, then so can photons via the observed emission lines, leading to more luminous spectral features. It is also shown that SNe IIP with larger progenitor stars have ejecta with a more physically extended oxygen layer that is well-mixed with the hydrogen layer. In addition, we find a subset of objects with evidence for asymmetric Ni-56 ejection, likely bipolar in shape. We also compare our observations to theoretical late-time spectral models of SNe IIP from two separate groups and find moderate-to-good agreement with both sets of models. Our SNe IIP spectra are consistent with models of 12-15 M_Sun progenitor stars having relatively low metallicity (Z $\le$ 0.01).
Refractive optical elements are widely used in millimeter and sub-millimeter astronomical telescopes. High resistivity silicon is an excellent material for dielectric lenses given its low loss-tangent, high thermal conductivity and high index of refraction. The high index of refraction of silicon causes a large Fresnel reflectance at the vacuum-silicon interface (up to 30%), which can be reduced with an anti-reflection (AR) coating. In this work we report techniques for efficiently AR coating silicon at sub-millimeter wavelengths using Deep Reactive Ion Etching (DRIE) and bonding the coated silicon to another silicon optic. Silicon wafers of 100 mm diameter (1 mm thick) were coated and bonded using the Silicon Direct Bonding technique at high temperature (1100 C). No glue is used in this process. Optical tests using a Fourier Transform Spectrometer (FTS) show sub-percent reflections for a single-layer DRIE AR coating designed for use at 320 microns on a single wafer. Cryogenic (10 K) measurements of a bonded pair of AR-coated wafers also reached sub-percent reflections. A prototype two-layer DRIE AR coating to reduce reflections and increase bandwidth is presented and plans for extending this approach are discussed.
We analyze the total and baryonic acceleration profiles of a set of well-resolved galaxies identified in the EAGLE suite of hydrodynamic simulations. Our runs start from the same initial conditions but adopt different subgrid models for stellar and AGN feedback, resulting in diverse populations of galaxies by the present day. Some of them reproduce observed galaxy scaling relations, while others do not. However, regardless of the feedback implementation, all of our galaxies follow closely a simple relationship between the total and baryonic acceleration profiles, consistent with recent observations of rotationally supported galaxies. The relation has small scatter: different feedback processes -- which produce different galaxy populations -- mainly shift galaxies along the relation, rather than perpendicular to it. Furthermore, galaxies exhibit a single characteristic acceleration, $g_{\dagger}$, above which baryons dominate the mass budget, as observed. These observations have been hailed as evidence for modified Newtonian dynamics but can be accommodated within the standard cold dark matter paradigm.
Observations using the 7 mm receiver system on the Australia Telescope Compact Array have revealed large reservoirs of molecular gas in two high-redshift radio galaxies: HATLAS J090426.9+015448 (z = 2.37) and HATLAS J140930.4+003803 (z = 2.04). Optically the targets are very faint, and spectroscopy classifies them as narrow-line radio galaxies. In addition to harbouring an active galactic nucleus the targets share many characteristics of sub-mm galaxies. Far-infrared data from Herschel-ATLAS suggest high levels of dust (>10^9 M_solar) and a correspondingly large amount of obscured star formation (~1000 M_solar / yr). The molecular gas is traced via the J = 1-0 transition of 12CO, its luminosity implying total H_2 masses of (1.7 +/- 0.3) x 10^11 and (9.5 +/- 2.4) x 10^10 ({\alpha}_CO/0.8) M_solar in HATLAS J090426.9+015448 and HATLAS J140930.4+003803 respectively. Both galaxies exhibit molecular line emission over a broad (~1000 km/s) velocity range, and feature double-peaked profiles. We interpret this as evidence of either a large rotating disk or an on-going merger. Gas depletion timescales are ~100 Myr. The 1.4 GHz radio luminosities of our targets place them close to the break in the luminosity function. As such they represent `typical' z > 2 radio sources, responsible for the bulk of the energy emitted at radio wavelengths from accretion-powered sources at high redshift, and yet they rank amongst the most massive systems in terms of molecular gas and dust content. We also detect 115 GHz rest-frame continuum emission, indicating a very steep high-radio-frequency spectrum, possibly classifying the targets as compact steep spectrum objects.
The quasar Q0918+1636 (z=3.07) has an intervening high-metallicity Damped Lyman-alpha Absorber (DLA) along the line of sight, at a redshift of z=2.58. The DLA is located at a large impact parameter of 16.2 kpc, and has an almost solar metallicity. It is shown, that a novel type of cosmological galaxy formation models, invoking a new SNII feedback prescription, the Haardt & Madau (2012) UVB field and explicit treatment of UVB self-shielding, can reproduce the observed characteristics of the DLA. UV radiation from young stellar populations in the galaxy, in particular in the photon energy range 10.36-13.61 eV (relating to Sulfur II abundance), are also considered in the analysis. It is found that a) for L~L* galaxies (at z=2.58), about 10% of the sight-lines through the galaxies at impact parameter 16.2 kpc will display a Sulfur II column density N(SII)$>$ 10$^{15.82}$ cm$^{-2}$ (the observed value for the DLA), and b) considering only cases where a near-solar metallicity will be detected at 16.2 kpc impact parameter, the probability distribution of galaxy SFR peaks near the value observed for the DLA galaxy counterpart of ~27 Msun/yr. It is argued, that the bulk of the alpha-elements, like Sulfur, traced by the high metal column density, b=16.2 kpc absorption lines, were produced by evolving young stars in the inner galaxy, and later transported outward by galactic winds.
We study supersymmetric extension of the Einstein-aether gravitational model where local Lorentz invariance is broken down to the subgroup of spatial rotations by a vacuum expectation value of a timelike vector field. By restricting to the level of linear perturbations around Lorentz-violating vacuum and using the superfield formalism we construct the most general action invariant under the linearized supergravity transformations. We show that, unlike its non-supersymmetric counterpart, the model contains only a single free dimensionless parameter, besides the usual dimensionful gravitational coupling. This makes the model highly predictive. An analysis of the spectrum of physical excitations reveal superluminal velocity of gravitons. The latter property leads to the extension of the gravitational multiplet by additional fermonic and bosonic states with helicities $\pm 3/2$ and $\pm 1$. We outline the observational constraints on the model following from its low-energy phenomenology.
Einstein's theory of gravity has been extensively tested on solar system scales, and for isolated astrophysical systems, using the perturbative framework known as the parameterized post-Newtonian (PPN) formalism. This framework is designed for use in the weak-field and slow-motion limit of gravity, and can be used to constrain a large class of metric theories of gravity with data collected from the aforementioned systems. Given the potential of future surveys to probe cosmological scales to high precision, it is a topic of much contemporary interest to construct a similar framework to link Einstein's theory of gravity and its alternatives to observations on cosmological scales. Our approach to this problem is to adapt and extend the existing PPN formalism for use in cosmology. We derive a set of equations that use the same parameters to consistently model both weak fields and cosmology. This allows us to parameterize a large class of modified theories of gravity and dark energy models on cosmological scales, using just four functions of time. These four functions can be directly linked to the background expansion of the universe, first-order cosmological perturbations, and the weak-field limit of the theory. They also reduce to the standard PPN parameters on solar system scales. We illustrate how dark energy models and scalar-tensor and vector-tensor theories of gravity fit into this framework, which we refer to as "parameterized post-Newtonian cosmology" (PPNC).
Links to: arXiv, form interface, find, astro-ph, recent, 1610, contact, help (Access key information)
There currently exist many observations which are not consistent with the cosmological principle. We review these observations with a particular emphasis on those relevant for Square Kilometre Array (SKA). In particular, several different data sets indicate a preferred direction pointing approximately towards the Virgo cluster. We also observe a hemispherical anisotropy in the Cosmic Microwave Background Radiation (CMBR) temperature fluctuations. Although these inconsistencies may be attributed to systematic effects, there remains the possibility that they indicate new physics and various theories have been proposed to explain them. One possibility, which we discuss in this review, is the generation of perturbation modes during the early pre-inflationary epoch, when the Universe may not obey the cosmological principle. Better measurements will provide better constraints on these theories. In particular, we propose measurement of the dipole in number counts, sky brightness, polarized flux and polarization orientations of radio sources. We also suggest test of alignment of linear polarizations of sources as a function of their relative separation. Finally we propose measurement of hemispherical anisotropy or equivalently dipole modulation in radio sources.
Detection of individual luminous sources during the reionization epoch and
cosmic dawn through their signatures in the HI 21-cm signal is one of the
direct approaches to probe the epoch. Here, we summarize our previous works on
this and present preliminary results on the prospects of detecting such sources
using the SKA1-low experiment. We first discuss the expected HI 21-cm signal
around luminous sources at different stages of reionization and cosmic dawn. We
then introduce two visibility based estimators for detecting such signal: one
based on the matched filtering technique and the other relies on simply combing
the visibility signal from different baselines and frequency channels.
We find that that the SKA1-low should be able to detect ionized bubbles of
radius $R_b \gtrsim 10$ Mpc with $\sim 100$ hr of observations at redshift $z
\sim 8$ provided that the mean outside neutral Hydrogen fraction $ x_{\rm HI}
\gtrsim 0.5$. We also investigate the possibility of detecting HII regions
around known bright QSOs such as around ULASJ1120+0641 discovered by Mortlock
et al. 2011. We find that a $5 \sigma$ detection is possible with $600$ hr of
SKA1-low observations if the QSO age and the outside $ x_{\rm HI} $ are at
least $\sim 2 \times 10^7$ Myr and $\sim 0.2$ respectively.
Finally, we investigate the possibility of detecting the very first X-ray and
Ly-$\alpha$ sources during the cosmic dawn. We consider mini-QSOs like sources
which emits in X-ray frequency band. We find that with a total $\sim 1000$ hr
of observations, SKA1-low should be able to detect those sources individually
with a $\sim 9 \sigma$ significance at redshift $z=15$. We summarize how the
SNR changes with various parameters related to the source properties.
Studying the cosmic dawn and the epoch of reionization through the redshifted 21 cm line are among the major science goals of the SKA1. Their significance lies in the fact that they are closely related to the very first stars in the universe. Interpreting the upcoming data would require detailed modelling of the relevant physical processes. In this article, we focus on the theoretical models of reionization that have been worked out by various groups working in India with the upcoming SKA in mind. These models include purely analytical and semi-numerical calculations as well as fully numerical radiative transfer simulations. The predictions of the 21 cm signal from these models would be useful in constraining the properties of the early galaxies using the SKA data.
The line of sight direction in the redshifted 21-cm signal coming from the cosmic dawn and the epoch of reionization is quite unique in many ways compared to any other cosmological signal. Different unique effects, such as the evolution history of the signal, non-linear peculiar velocities of the matter etc will imprint their signature along the line of sight axis of the observed signal. One of the major goals of the future SKA-LOW radio interferometer is to observe the cosmic dawn and the epoch of reionization through this 21-cm signal. It is thus important to understand how these various effects affect the signal for it's actual detection and proper interpretation. For more than one and half decades, various groups in India have been actively trying to understand and quantify the different line of sight effects that are present in this signal through analytical models and simulations. In many ways the importance of this sub-field under 21-cm cosmology have been identified, highlighted and pushed forward by the Indian community. In this article we briefly describe their contribution and implication of these effects in the context of the future surveys of the cosmic dawn and the epoch of reionization that will be conducted by the SKA-LOW.
Tomographic intensity mapping of the HI using the redshifted 21 cm observations opens up a new window towards our understanding of cosmological background evolution and structure formation. This is a key science goal of several upcoming radio telescopes including the Square Kilometer Array (SKA). In this article we focus on the post-reionization signal and investigate the of cross correlating the 21 cm signal with other tracers of the large scale structure. We consider the cross-correlation of the post-reionization 21 cm signal with the Lyman-alpha forest, Lyman-break galaxies and late time anisotropies in the CMBR maps like weak lensing and the Integrated Sachs Wolfe effect. We study the feasibility of detecting the signal and explore the possibility of obtaining constraints on cosmological models using it.
The intra-cluster and inter-galactic media (ICM, IGM) that pervade the large scale structure of the Universe are known to be magnetised at sub-micro Gauss to micro Gauss levels and to contain cosmic rays (CRs). The acceleration of CRs and their evolution along with that of magnetic fields in these media is still not well understood. Diffuse radio sources of synchrotron origin associated with the ICM such as radio halos, relics and mini-halos are direct probes of the underlying mechanisms of CR acceleration. Observations with radiotelescopes such as the GMRT, the VLA and the WSRT (0.15 - 2 GHz) have revealed scaling relations between the thermal and non-thermal properties of clusters and favour the role of shocks in the formation of radio relics and of turbulent re-acceleration in the formation of radio halos and mini-halos. Due to the limitations of current radio telescopes, wide-band studies and exploration of low mass and supercluster-scale systems is difficult. The Square Kilometer Array (SKA) is a next generation radio telescope that will operate in the frequency range of 0.05 - 20 GHz with unprecedented sensitivities and resolutions. The expected detection limits of SKA will reveal a few hundred to thousand new radio halos, relics and mini-halos providing the first large and comprehensive samples for their study. The wide frequency coverage along with sensitivity to extended structures will be able to constrain the CR acceleration mechanisms. The higher frequency (> 5 GHz) observations will be able to use the Sunyaev-Zel'dovich effect to probe the ICM pressure in addition to the tracers such as lobes of head-tail radio sources. The SKA also opens prospects to detect the "off-state" radio emission from the ICM predicted by the hadronic models and the turbulent re-acceleration models. [abridged]
We discuss the prospects of using the redshifted 21~cm emission from neutral hydrogen in the post-reionization epoch to study our universe. The main aim of the article is to highlight the efforts of Indian scientists in this area with the SKA in mind. It turns out that the intensity mapping surveys from SKA can be instrumental in obtaining tighter constraints on the dark energy models. Cross-correlation of the HI intensity maps with the Ly$\alpha$ forest data can also be useful in measuring the BAO scale.
The diffuse Galactic syncrotron emission (DGSE) is the most important diffuse foreground component for future cosmological 21-cm observations. The DGSE is also an important probe of the cosmic ray electron and magnetic field distributions in the turbulent interstellar medium (ISM) of our Galaxy. In this paper we briefly review the Tapered Gridded Estimator (TGE) which can be used to quantify the angular power spectrum of the sky signal directly from the visibilities measured in radio-interferometric observations. The salient features of the TGE are (1.) it deals with the gridded data which makes it computationally very fast (2.) it avoids a positive noise bias which normally arises from the system noise inherent to the visibility data, and (3.) it allows us to taper the sky response and thereby suppresses the contribution from unsubtracted point sources in the outer parts and the sidelobes of the antenna beam pattern. We also summarize earlier work where the TGE was used to measure the C_l of the DGSE using 150 MHz GMRT data. Earlier measurements of the angular power spectrum are restricted to smaller angular multipole l less than ~ 10^3 for the DGSE, the signal at the larger l values is dominated by the residual point sources after source subtraction. The higher sensitivity of the upcoming SKA1 Low will allow the point sources to be subtracted to a fainter level than possible with existing telescopes. We predict that it will be possible to measure the angular power spectrum of the DGSE to larger values of l with SKA1 Low. Our results show that it should be possible to achieve l_{max} ~ 10^4 and ~ 10^5 with 2 minutes and 10 hours of observations respectively.
An intriguing alternative to cold dark matter (CDM) is that the dark matter is a light ( $m \sim 10^{-22}$ eV) boson having a de Broglie wavelength $\lambda \sim 1$ kpc, often called fuzzy dark matter (FDM). We describe the arguments from particle physics that motivate FDM, review previous work on its astrophysical signatures, and analyze several unexplored aspects of its behavior. In particular, (i) FDM halos smaller than about $10^7 (m/10^{-22} {\rm eV})^{-3/2} M_\odot$ do not form. (ii) FDM halos are comprised of a core that is a stationary, minimum-energy configuration called a "soliton", surrounded by an envelope that resembles a CDM halo. (iii) The transition between soliton and envelope is determined by a relaxation process analogous to two-body relaxation in gravitating systems, which proceeds as if the halo were composed of particles with mass $\sim \rho\lambda^3$ where $\rho$ is the halo density. (iv) Relaxation may have substantial effects on the stellar disk and bulge in the inner parts of disk galaxies. (v) Relaxation can produce FDM disks but an FDM disk in the solar neighborhood must have a half-thickness of at least $300 (m/10^{-22} {\rm eV})^{-2/3}$ pc. (vi) Solitonic FDM sub-halos evaporate by tunneling through the tidal radius and this limits the minimum sub-halo mass inside 30 kpc of the Milky Way to roughly $10^8 (m/10^{-22} {\rm eV})^{-3/2} M_\odot$. (vii) If the dark matter in the Fornax dwarf galaxy is composed of CDM, most of the globular clusters observed in that galaxy should have long ago spiraled to its center, and this problem is resolved if the dark matter is FDM.
I investigate the consistency of the VIMOS Public Extragalactic Redshift
Survey v7 galaxy sample with the expansion history and linear growth rate
predicted by General Relativity (GR) and a Planck (2015) cosmology. To do so, I
measure the redshift-space power spectrum, which is anisotropic due to both
redshift-space distortions (RSD) and the Alcock-Paczynski (AP) effect. In
Chapter 6, I place constraints of $f \sigma_8(0.76) = 0.44 \pm 0.04$ and $f
\sigma_8(1.05) = 0.28 \pm 0.08$, which remain consistent with GR at 95%
confidence. Marginalising over the anisotropic AP effect degrades the
constraints by a factor of three but allows $F_{AP} \equiv (1+z) D_A H/c$ to be
simultaneously constrained. The VIPERS v7 joint-posterior on $(f \sigma_8,
F_{AP})$ shows no compelling deviation from GR.
Chapter 7 investigates the inclusion of a simple density transform:
`clipping' prior to the RSD analysis. This tackles the root-cause of
non-linearity and may extend the validity of perturbation theory. Moreover,
this marked statistic would amplify signatures of shielded modified gravity
models and includes information not available to the power spectrum. I show
that a linear real-space power spectrum with a Kaiser factor and a Lorentzian
damping yields a significant bias without clipping, but that this may be
removed with a strict threshold; similar behaviour is observed for the data.
Estimates of $f \sigma_8$ for different thresholds are highly correlated, but
this may be obtained using mocks. A maximum likelihood estimate from a
combination of thresholds is shown to achieve a 16% decrease in statistical
error relative to a single-threshold estimate. The results are encouraging to
date but represent a work in progress; the final analysis will be submitted to
Astronomy & Astrophysics as Wilson et al. (2016).
The properties of large underdensities in the distribution of galaxies in the Universe, known as cosmic voids, are potentially sensitive probes of fundamental physics. We use data from the MultiDark suite of N-body simulations and multiple halo occupation distribution mocks to study the relationship between galaxy voids and the gravitational potential $\Phi$. We find that the majority of galaxy voids correspond to local density minima in larger-scale overdensities, and thus lie in potential wells. However, a subset of voids can be identified that closely trace maxima of the gravitational potential and thus stationary points of the velocity field. We identify a new void observable, $\lambda_v$, which depends on a combination of the void size and the average galaxy density contrast within the void, and show that it provides a good proxy indicator of the potential at the void location. A simple linear scaling of $\Phi$ as a function of $\lambda_v$ is found to hold, independent of the redshift and properties of the galaxies used as tracers of voids. We provide an accurate fitting formula to describe the spherically averaged potential profile $\Phi(r)$ about void centre locations. We discuss the importance of these results for the understanding of the evolution history of voids, and for their use in precision measurements of the integrated Sachs-Wolfe effect, gravitational lensing and peculiar velocity distortions in redshift space.
The black hole merging rates inferred after the gravitational-wave detection by Advanced LIGO/VIRGO and the relatively high mass of the progenitors are consistent with models of dark matter made of massive primordial black holes (PBH). PBH binaries emit gravitational waves in a broad range of frequencies that will be probed by future space interferometers (LISA) and pulsar timing arrays (PTA). The amplitude of the stochastic gravitational-wave background expected for PBH dark matter is calculated taking into account various effects such as initial eccentricity of binaries, PBH velocities, mass distribution and clustering. It allows a detection by the LISA space interferometer, and possibly by the PTA of the SKA radio-telescope. Interestingly, one can distinguish this background from the one of non-primordial massive binaries through a specific frequency dependence, resulting from the maximal impact parameter of binaries formed by PBH capture, depending on the PBH velocity distribution and their clustering properties. Moreover, we find that the gravitational wave spectrum is boosted by the width of PBH mass distribution, compared with that of the monochromatic spectrum. The current PTA constraints already rule out broad-mass PBH models covering more than three decades of masses, but evading the microlensing and CMB constraints due to clustering.
High-energy photons (above the MeV) are a powerful probe for astrophysics and for fundamental physics under extreme conditions. During the recent years, our knowledge of the high-energy gamma-ray sky has impressively progressed thanks to the advent of new detectors for cosmic gamma rays, at ground (H.E.S.S., MAGIC, VERITAS, HAWC) and in space (AGILE, Fermi). This presentation reviews the present status of the studies of fundamental physics problems with high-energy gamma rays, and discusses the expected experimental developments.
In a general effective theory description of inflation a disformal transformation can be used to set the tensor sound speed to one. After the transformation, the tensor power spectrum then automatically only depends on the Hubble parameter. We show that this disformal transformation, however, is nothing else than a change of units. It is a very useful tool for simplifying and interpreting computations, but it cannot change any physics. While the apparent parametrical dependence of the tensor power spectrum does change under a disformal transformation, the physics described is frame invariant. We further illustrate the frame invariance of the tensor power spectrum by writing it exclusively in terms of separately invariant quantities.
Links to: arXiv, form interface, find, astro-ph, recent, 1610, contact, help (Access key information)
We quantify the error in the results of mixed baryon--dark-matter hydrodynamic simulations, stemming from outdated approximations for the generation of initial conditions. The error at redshift 0 in contemporary large simulations, is of the order of few to ten percent in the power spectra of baryons and dark matter, and their combined total-matter power spectrum. After describing how to properly assign initial displacements and peculiar velocities to multiple species, we review several approximations: (1) {using the total-matter power spectrum to compute displacements and peculiar velocities of both fluids}, (2) scaling the linear redshift-zero power spectrum back to the initial power spectrum using the Newtonian growth factor ignoring homogeneous radiation, (3) using longitudinal-gauge velocities with synchronous-gauge densities, and (4) ignoring the phase-difference in the Fourier modes for the offset baryon grid, relative to the dark-matter grid. Three of these approximations do not take into account that dark matter and baryons experience a scale-dependent growth after photon decoupling, which results in directions of velocity which are not the same as their direction of displacement. We compare the outcome of hydrodynamic simulations with these four approximations to our reference simulation, all setup with the same random seed and simulated using Gadget-III.
The arrival times, positions, and fluxes of multiple images in strong lens systems can be used to infer the presence of dark subhalos in the deflector, and thus test predictions of cold dark matter models. However, gravitational lensing does not distinguish between perturbations to a smooth gravitational potential arising from baryonic and non-baryonic mass. In this work, we quantify the extent to which the stellar mass distribution of a deflector can reproduce flux ratio and astrometric anomalies typically associated with the presence of a dark matter subhalo. Using Hubble Space Telescope images of nearby galaxies, we simulate strong lens systems with real distributions of stellar mass as they would be observed at redshift $z_d=0.5$. We add a dark matter halo and external shear to account for the smooth dark matter field, omitting dark substructure, and use a Monte Carlo procedure to characterize the distributions of image positions, time delays, and flux ratios for a compact background source of diameter 5 pc. By convolving high-resolution images of real galaxies with a Gaussian PSF, we simulate the most detailed smooth potential one could construct given high quality data, and find scatter in flux ratios of $\approx 10\%$, which we interpret as a typical deviation from a smooth potential caused by large and small scale structure in the lensing galaxy. We demonstrate that the flux ratio anomalies arising from galaxy-scale baryonic structure can be minimized by selecting the most massive and round deflectors, and by simultaneously modeling flux ratio and astrometric data.
Directional detection is an important way to detect dark matter. An input to these experiments is the dark matter velocity distribution. Recent hydrodynamical simulations have shown that the dark matter velocity distribution differs substantially from the Standard Halo Model. We study the impact of some of these updated velocity distribution in dark matter directional detection experiments. We calculate the ratio of events required to confirm the forward-backward asymmetry and the existence of the ring of maximum recoil rate using different dark matter velocity distributions for $^{19}$F and Xe targets. We show that with the use of updated dark matter velocity profiles, the forward-backward asymmetry and the ring of maximum recoil rate can be confirmed using a factor of $\sim$2 -- 3 less events when compared to that using the Standard Halo Model.
Anisotropies in distortions to the frequency spectrum of the cosmic microwave background (CMB) can be created through spatially varying heating processes in the early Universe. For instance, the dissipation of small-scale acoustic modes does create distortion anisotropies, in particular for non-Gaussian primordial perturbations. In this work, we derive approximations that allow describing the associated distortion field. We provide a systematic formulation of the problem using Fourier-space window functions, clarifying and generalizing previous approximations. Our expressions highlight the fact that the amplitudes of the spectral-distortion fluctuations induced by non-Gaussianity depend also on the homogeneous value of those distortions. Absolute measurements are thus required to obtain model-independent distortion constraints on primordial non-Gaussianity. We also include a simple description for the evolution of distortions through photon diffusion, showing that these corrections can usually be neglected. Our formulation provides a systematic framework for computing higher order correlation functions of distortions with CMB temperature anisotropies and can be extended to describe correlations with polarization anisotropies.
Advanced LIGO's discovery of gravitational-wave events GW150914 and GW151226 has stimulated extensive studies on the origin of binary black holes. Supposing the gravitational-wave events could be explained by binary primordial black hole (PBH) mergers, we investigate the corresponding stochastic gravitational-wave background (SGWB) and point out the possibility to detect this SGWB spectrum, in particular from the subsolar mass PBHs, by the Advanced LIGO in the near future. We also use the non-detection of SGWB to give a new independent constraint on the abundance of PBHs in dark matter.
We explore the warm inflation scenario theoretical predictions looking at two different dissipative regimes for several representative primordial potentials. As it is well known, the warm inflation is able to decrease the tensor-to-scalar ratio value, rehabilitating several primordial potential ruled out in the cold inflation context by the recent cosmic microwave background data. Here we show that it is also able to produce a running of the running $n_s"$ positive and within the Planck data limits. This is very remarkable since the standard cold inflation model is unable to justify the current indication of a positive constraint on $n_s"$. We achieve a parameterization for the primordial power spectrum able to take into account higher order effects as the running of the spectral index and the running of the running, and we perform statistical analysis using the most up-to-date Planck data to constrain the dissipative effects. We find that the warm inflation can explain the current observables with a good statistical significance, even for those potentials ruled out in the simplest cold inflation scenario.
We study the impact of assumptions about neutrino properties on the estimation of inflationary parameters from cosmological data, with a specific focus on the allowed contours in the $n_s/r$ plane. We study the following neutrino properties: (i) the total neutrino mass $ M_\nu =\sum_i m_i$; (ii) the number of relativistic degrees of freedom $N_{eff}$; and (iii) the neutrino hierarchy: whereas previous literature assumed 3 degenerate neutrino masses or two massless neutrino species (that do not match neutrino oscillation data), we study the cases of normal and inverted hierarchy. Our basic result is that these three neutrino properties induce $< 1 \sigma$ shift of the probability contours in the $n_s/r$ plane with both current or upcoming data. We find that the choice of neutrino hierarchy has a negligible impact. However, the minimal cutoff on the total neutrino mass $M_{\nu,{min}}=0 $ that accompanies previous works using the degenerate hierarchy does introduce biases in the $n_s/r$ plane and should be replaced by $M_{\nu,min}= 0.059$ eV as required by oscillation data. Using current CMB data from Planck and Bicep/Keck, marginalizing over $ M_\nu$ and over $r$ can lead to a shift in the mean value of $n_s$ of $\sim0.3\sigma$ towards lower values. However, once BAO measurements are included, the standard contours in the $n_s/r$ plane are basically reproduced. Larger shifts of the contours in the $n_s/r$ plane (up to 0.8$\sigma$) arise for nonstandard values of $N_{eff}$. We also provide forecasts for the future CMB experiments COrE and Stage-IV and show that the incomplete knowledge of neutrino properties, taken into account by a marginalization over $M_\nu$, could induce a shift of $\sim0.4\sigma$ towards lower values in the determination of $n_s$ (or a $\sim 0.8\sigma$ shift if one marginalizes over $N_{eff}$). Comparison to specific inflationary models is shown.
This is the second in a series of papers on preheating in inflationary models comprised of multiple scalar fields coupled nonminimally to gravity. In this paper, we work in the rigid-spacetime approximation and consider field trajectories within the single-field attractor, which is a generic feature of these models. We construct the Floquet charts to find regions of parameter space in which particle production is efficient for both the adiabatic and isocurvature modes, and analyze the resonance structure using analytic and semi-analytic techniques. Particle production in the adiabatic direction is characterized by the existence of an asymptotic scaling solution at large values of the nonminimal couplings, $\xi_I \gg 1$, in which the dominant instability band arises in the long-wavelength limit, for comoving wavenumbers $k \rightarrow 0$. However, the large-$\xi_I$ regime is not reached until $\xi_I \geq {\cal O} (100)$. In the intermediate regime, with $\xi_I \sim {\cal O}(10)$, the resonance structure depends strongly on the wavenumber and couplings. The resonance structure for isocurvature perturbations is distinct and more complicated than its adiabatic counterpart, as expected. An intermediate regime, for $\xi_I \sim {\cal O} (10)$, is again evident. For large values of $\xi_I$, the Floquet chart consists of densely spaced, nearly parallel instability bands, suggesting a very efficient preheating behavior. Quantitatively, the approach to the large-$\xi_I$ asymptotic solution for isocurvature modes is slower than in the case of the adiabatic modes.
The vector-tensor (VT) theory of gravitation revisited in this article was studied in previous papers, where it was proved that VT works and deserves attention. New observational data and numerical codes have motivated further development which is presented here. New research has been planed with the essential aim of proving that current cosmological observations, including Planck data, baryon acoustic oscillations (BAO), and so on, may be explained with VT, a theory which accounts for a kind of dark energy which has the same equation of state as vacuum. New versions of the codes CAMB and COSMOMC have been designed for applications to VT, and the resulting versions have been used to get the cosmological parameters of the VT model at suitable confidence levels. The parameters to be estimated are the same as ingeneral relativity (GR), plus a new parameter $D$. For $D=0$, VT linear cosmological perturbations reduces to those of GR, but the VT background may explain dark energy. The fits between observations and VT predictions lead to non vanishing $|D|$ upper limits at the $1\sigma $ confidence level. The value $D=0$ is admissible at this level, but this value is not that of the best fit in any case. Results strongly suggest that VT may explain current observations, at least, as well as GR; with the advantage that, as it is proved in this paper, VT has an additional parameter which facilitates adjustments to current observational data.}
We propose and construct a two-parameter perturbative expansion around a Friedmann-Lema\^{i}tre-Robertson-Walker geometry that can be used to model high-order gravitational effects in the presence of non-linear structure. This framework reduces to the weak-field and slow-motion post-Newtonian treatment of gravity in the appropriate limits, but also includes the low-amplitude large-scale fluctuations that are important for cosmological modelling. We derive a set of field equations that can be applied to the late Universe, where non-linear structure exists on supercluster scales, and perform a detailed investigation of the associated gauge problem. This allows us to identify a consistent set of perturbed quantities in both the gravitational and matter sectors, and to construct a set of gauge-invariant quantities that correspond to each of them. The field equations, written in terms of these quantities, take on a relatively simple form, and allow the effects of small-scale structure on the large-scale properties of the Universe to be clearly identified. We find that inhomogeneous structures source the global expansion, that there exist new field equations at new orders, and that there is vector gravitational potential that is a hundred times larger than one might naively expect from cosmological perturbation theory. Finally, we expect our formalism to be of use for calculating relativistic effects in upcoming ultra-large-scale surveys, as the form of the gravitational coupling between small and large scales depends on the non-linearity of Einstein's equations, and occurs at what is normally thought of as first order in cosmological perturbations.
This paper concludes our semi-analytic study of preheating in inflationary models comprised of multiple scalar fields coupled nonminimally to gravity. Using the covariant framework of Ref. [1], we extend the rigid-spacetime results of Ref. [2] by considering both the expansion of the universe during preheating, as well as the effect of the coupled metric perturbations on particle production. The adiabatic and isocurvature perturbations are governed by different effective masses that scale differently with the nonminimal couplings and evolve differently in time. The effective mass for the adiabatic modes is dominated by contributions from the coupled metric perturbations immediately after inflation. The metric perturbations contribute an oscillating tachyonic term that enhances an early period of significant particle production for the adiabatic modes, which ceases on a time-scale governed by the nonminimal couplings $\xi_I$. The effective mass of the isocurvature perturbations, on the other hand, is dominated by contributions from the fields' potential and from the curvature of the field-space manifold (in the Einstein frame), the balance between which shifts on a time-scale governed by $\xi_I$. As in Refs. [1,2], we identify distinct behavior depending on whether the nonminimal couplings are small ($\xi_I \lesssim {\cal O} (1)$), intermediate ($\xi_I \sim {\cal O} (10)$), or large ($\xi_I \geq 10^2$).
Recent analyses in the literature suggest that the concordance $\Lambda$CDM model with rigid cosmological term, $\Lambda=$const., may not be the best description of the cosmic acceleration. The class of "running vacuum models", in which $\Lambda=\Lambda(H)$ evolves with the Hubble rate, has been shown to fit the string of $SNIa+BAO+H(z)+LSS+CMB$ data significantly better than the $\Lambda$CDM. Here we provide further evidence on the time-evolving nature of the dark energy (DE) by fitting the same cosmological data in terms of scalar fields. As a representative model we use the original Peebles & Ratra potential, $V\propto\Phi^{-\alpha}$. We find clear signs of dynamical DE at $\sim 4\sigma$ c.l., thus reconfirming through a nontrivial scalar field approach the strong hints formerly found with other models and parametrizations.
The accelerating expansion of the universe is one of the most profound discoveries in modern cosmology, pointing to a universe in which 70% of the mass-energy density has an unknown form spread uniformly across the universe. This result has been well established using a combination of cosmological probes (e.g., Planck Collaboration et al. 2016), resulting in a "standard model" of modern cosmology that is a combination of a cosmological constant with cold dark matter and baryons. The first compelling evidence for this acceleration came in the late 1990's, when two independent teams studying type Ia supernovae discovered that distant SNe Ia were dimmer than expected. The combined analysis of modern cosmology experiments, including SNe Ia (Betoule et al. 2014), the cosmic microwave background (Planck Collaboration et al. 2016), and baryon acoustic oscillations (Alam et al. 2016), indicate ~ 75 sigma evidence for positive Omega_Lambda. A recent study has claimed that the evidence for acceleration from SNe Ia is marginal. Here we demonstrate errors in that analysis which reduce the significance from SNe Ia, and show that constraints on the flatness or matter density of the universe greatly increase the significance of acceleration. Analyzing the Joint Light-curve Analysis supernova sample, we find 4.2 sigma evidence for acceleration with SNe Ia alone, and 11.2 sigma in a flat universe. With the correct supernova analysis and by not rejecting all other cosmological constraints, we find that acceleration is quite secure.
We study the dynamics of galaxy mergers, with emphasis on the gas feeding of nuclear regions, using a suite of hydrodynamical simulations of galaxy encounters. The high spatial and temporal resolution of the simulations allows us to not only recover the standard picture of tidal-torque induced inflows, but also to detail another, important feeding path produced by ram pressure. The induced shocks effectively decouple the dynamics of the gas from that of the stars, greatly enhancing the loss of gas angular momentum and leading to increased central inflows. The ram-pressure shocks also cause, in many cases, the entire galactic gas disc of the smaller galaxy to abruptly change its direction of rotation, causing a complete "flip" and, several $10^8$ yr later, a subsequent "counter-flip". This phenomenon results in the existence of long-lived decoupled gas-stellar and stellar-stellar discs, which could hint at a new explanation for the origin of some of the observed kinematically decoupled cores/counter-rotating discs. Lastly, we speculate, in the case of non-coplanar mergers, on the possible existence of a new class of remnant systems similar to some of the observed X-shaped radio galaxies.
We use cosmological simulations from the FIRE (Feedback In Realistic Environments) project to study the baryon cycle and galaxy mass assembly for central galaxies in the halo mass range $M_{\rm halo} \sim 10^{10} - 10^{13} M_{\odot}$. By tracing cosmic inflows, galactic outflows, gas recycling, and merger histories, we quantify the contribution of physically distinct sources of material to galaxy growth. We show that in situ star formation fueled by fresh accretion dominates the early growth of galaxies of all masses, while the re-accretion of gas previously ejected in galactic winds often dominates the gas supply for a large portion of every galaxy's evolution. Externally processed material contributes increasingly to the growth of central galaxies at lower redshifts. This includes stars formed ex situ and gas delivered by mergers, as well as smooth intergalactic transfer of gas from other galaxies, an important but previously under-appreciated growth mode. By $z=0$, wind transfer, i.e. the exchange of gas between galaxies via winds, can dominate gas accretion onto $\sim L^{*}$ galaxies over fresh accretion and standard wind recycling. Galaxies of all masses re-accrete >50% of the gas ejected in winds and recurrent recycling is common. The total mass deposited in the intergalactic medium per unit stellar mass formed increases in lower mass galaxies. Re-accretion of wind ejecta occurs over a broad range of timescales, with median recycling times ($\sim 100-350$ Myr) shorter than previously found. Wind recycling typically occurs at the scale radius of the halo, independent of halo mass and redshift, suggesting a characteristic recycling zone around galaxies that scales with the size of the inner halo and the galaxy's stellar component.
We study two models of scalar dark matter from "large" electroweak multiplets with isospin $5/2$ ($n=6$ members) and $7/2$ ($n=8$), whose scalar potentials preserve a $Z_2$ symmetry. Because of large annihilation cross sections due to electroweak interactions, these scalars can constitute all the dark matter only for masses in the multi-TeV range. For such high masses, Sommerfeld enhancement and co-annihilations play important roles in the dark matter relic abundance calculation. We determine the allowed parameter ranges including both of these effects and show that these models are as yet unconstrained by dark matter direct detection experiments, but will be probed by currently-running and proposed future experiments. We also show that a Landau pole appears in these models at energy scales below $10^9$ GeV, indicating the presence of additional new physics below that scale.
The information extracted from large galaxy surveys with the likes of DES, DESI, Euclid, LSST, SKA, and WFIRST will be greatly enhanced if the resultant galaxy catalogues can be cross-correlated with one another. Predicting the nature of the information gain, and developing the tools to realise it, depends on establishing a consistent model of how the galaxies detected by each survey trace the same underlying matter distribution. Existing analytic methods, such as halo occupation distribution (HOD) modelling, are not well-suited for this task, and can suffer from ambiguities and tuning issues when applied to multiple tracers. We construct a simple alternative that provides a common model for the connection between galaxies and dark matter halos across a wide range of wavelengths (and thus tracer populations). This is based on a chain of parametrised statistical distributions that model the connection between (a) halo mass and bulk physical properties of galaxies, such as star-formation rate; and (b) those same physical properties and a variety of emission processes. The result is a flexible parametric model that allows analytic halo model calculations to be carried out for multiple tracers, as well as providing semi-realistic galaxy properties for fast mock catalogue generation.
We extend previous results on healthy derivative self-interactions for a Proca field to the case of a set of massive vector fields. We obtain non-gauge invariant derivative self-interactions for the vector fields that maintain the appropriate number of propagating degrees of freedom. In view of the potential cosmological applications, we restrict to interactions with an internal rotational symmetry. We provide a systematical construction order by order in derivatives of the fields and making use of the antisymmetric Levi-Civita tensor. We then compare with the one single vector field case and show that the interactions can be broadly divided into two groups, namely the ones obtained from a direct extension of the generalized Proca terms and genuine multi-Proca interactions with no correspondence in the single Proca case. We also discuss the curved spacetime version of the interactions to include the necessary non-minimal couplings to gravity. Finally, we explore the cosmological applications and show that there are three different vector fields configurations giving rise to isotropic solutions. Two of them have already been considered in the literature and the third one, representing a combination of the first two, is new and offers unexplored cosmological scenarios.
Links to: arXiv, form interface, find, astro-ph, recent, 1610, contact, help (Access key information)