We calibrate here cosmological radiative transfer simulation with ATON/RAMSES with a range of measurements of the Lyman alpha opacity from QSO absorption spectra. We find the Lyman alpha opacity to be very sensitive to the exact timing of hydrogen reionisation. Models reproducing the measured evolution of the mean photoionisation rate and average mean free path reach overlap at z ~ 7 and predict an accelerated evolution of the Lyman alpha opacity at z > 6 consistent with the rapidly evolving luminosity function of Lyman alpha emitters in this redshift range. Similar to "optically thin" simulations our full radiative transfer simulations fail, however, to reproduce the high-opacity tail of the Lyman alpha opacity PDF at z > 5. We argue that this is due to spatial UV fluctuations in the post-overlap phase of reionisation on substantially larger scales than predicted by our source model, where the ionising emissivity is dominated by large numbers of sub-L* galaxies. We further argue that this suggests a significant contribution to the ionising UV background by much rarer bright sources at high redshift.
We present a mass map reconstructed from weak gravitational lensing shear measurements over 139 sq. deg from the Dark Energy Survey (DES) Science Verification data. The mass map probes both luminous and dark matter, thus providing a tool for studying cosmology. We find good agreement between the mass map and the distribution of massive galaxy clusters identified using a red-sequence cluster finder. Potential candidates for super-clusters and voids are identified using these maps. We measure the cross-correlation between the mass map and a magnitude-limited foreground galaxy sample and find a detection at the 5-7 sigma level on a large range of scales. These measurements are consistent with simulated galaxy catalogs based on LCDM N-body simulations, suggesting low systematics uncertainties in the map. We summarize our key findings in this letter; the detailed methodology and tests for systematics are presented in a companion paper.
Distance measurements are currently the most powerful tool to study the expansion history of the universe without specifying its matter content nor any theory of gravitation. Assuming only an isotropic, homogeneous and flat universe, in this work we introduce a model-independent method to reconstruct directly the deceleration function via a piecewise function. Including a penalty factor, we are able to vary continuously the complexity of the deceleration function from a linear case to an arbitrary $(n+1)$-knots spline interpolation. We carry out a Monte Carlo analysis to determine the best penalty factor, evaluating the bias-variance trade-off, given the uncertainties of the SDSS-II and SNLS supernova combined sample (JLA), compilations of baryon acoustic oscillation (BAO) and $H(z)$ data. We show that, evaluating a single fiducial model, the conclusions about the bias-variance ratio are misleading. We determine the reconstruction method in which the bias represents at most $10\%$ of the total uncertainty. In all statistical analyses, we fit the coefficients of the deceleration function along with four nuisance parameters of the supernova astrophysical model. For the full sample, we also fit $H_0$ and the sound horizon $r_s(z_d)$ at the drag redshift. The bias-variance trade-off analysis shows that, apart from the deceleration function, all other estimators are unbiased. Finally, we apply the Ensemble Sampler Markov Chain Monte Carlo (ESMCMC) method to explore the posterior of the deceleration function up to redshift $1.3$ (using only JLA) and $2.3$ (JLA+BAO+$H(z)$). We obtain that the standard cosmological model agrees within $3\sigma$ level with the reconstructed results in the whole studied redshift intervals. Since our method is calibrated to minimize the bias, the error bars of the reconstructed functions are a good approximation for the total uncertainty.
Guided by non-relativistic Effective Field Theory (EFT) we classify the most general spin-dependent interactions between a fermionic Weakly Interacting Massive Particle (WIMP) and nuclei, and within this class of models we discuss the viability of an interpretation of the DAMA modulation result in terms of a signal from WIMP elastic scatterings using a halo-independent approach. We find that, although several relativistic EFT's can lead to a spin-dependent cross section, in some cases with an explicit, non-negligible dependence on the WIMP incoming velocity, three main scenarios can be singled out in the non-relativistic limit which approximately encompass them all, and that only differ by their dependence on the transferred momentum. For two of them compatibility between DAMA and other constraints is possible for a WIMP mass below 30 GeV, but only for a WIMP velocity distribution in the halo of our Galaxy which departs from a Maxwellian. This is achieved by combining a suppression of the WIMP effective coupling to neutrons (to evade constraints from xenon and germanium detectors) to an explicit quadratic or quartic dependence of the cross section on the transferred momentum (that leads to a relative enhancement of the expected rate off sodium in DAMA compared to that off fluorine in droplet detectors and bubble chambers). For larger WIMP masses the same scenarios are excluded by scatterings off iodine in COUPP.
We estimate the rate of dark matter scattering in collapsed structures throughout the history of the Universe. If the scattering cross-section is velocity-independent, then the canonical picture is correct that scatterings occur mainly at late times. The scattering rate peaks slightly at redshift z~6, and remains significant today. Half the scatterings occur after z~1, in structures more massive than 10^12 M_sun. Within a factor of two, these numbers are robust to changes in the assumed astrophysics, and the scatterings would be captured in cosmological simulations. However, for particle physics models with a velocity-dependent cross-section (as for Yukawa potential interactions via a massive mediator), the scattering rate peaks before z~20, in objects with mass less than 10^4 M_sun. These precise values are sensitive to the redshift-dependent mass-concentration relation and the small-scale cutoff in the matter power spectrum. In extreme cases, the qualitative effect of early interactions may be reminiscent of warm dark matter and strongly affect the subsequent growth of structure. However, these scatterings are being missed in existing cosmological simulations with limited mass resolution.
Time is a parameter playing a central role in our most fundamental modeling of natural laws. Relativity theory shows that the comparison of times measured by different clocks depends on their relative motions and on the strength of the gravitational field in which they are embedded. In standard cosmology, the time parameter is the one measured by fundamental clocks, i.e. clocks at rest with respect to the expanding space. This proper time is assumed to flow at a constant rate throughout the whole history of the Universe. We make the alternative hypothesis that the rate at which cosmological time flows depends on the dynamical state of the Universe. In thermodynamics, the arrow of time is strongly related to the second law, which states that the entropy of an isolated system will always increase with time or, at best, stay constant. Hence, we assume that time measured by fundamental clocks is proportional to the entropy of the region of the Universe that is causally connected to them. Under that simple assumption, we build a cosmological model that explains the Type Ia Supernovae data (the best cosmological standard candles) without the need for exotic dark matter nor dark energy.
Evidence for dark matter self-interactions has recently been reported based on the observation of a spatial offset between the dark matter halo and the stars in a galaxy in the cluster Abell 3827. Interpreting the offset as due to dark matter self-interactions leads to a cross section measurement of sigma_DM/m ~ (1-1.5) cm^2/g, where m is the mass of the dark matter particle. We use this observation to constrain singlet scalar dark matter coupled to the Standard Model and to two-Higgs-doublet models. We show that the most natural scenario in this class of models is very light dark matter, below about 0.1 GeV, whose relic abundance is set by freeze-in, i.e., by slow production of dark matter in the early universe via extremely tiny interactions with the Higgs boson, never reaching thermal equilibrium. We also show that the dark matter abundance can be established through the usual thermal freeze-out mechanism in the singlet scalar extension of the Yukawa-aligned two-Higgs-doublet model, but that it requires rather severe fine tuning of the singlet scalar mass.
We study the effects on the spectrum and distribution of high-energy neutrinos due to scattering with dark matter both outside and within our galaxy, focusing on the neutrinos observed by the IceCube experiment with energies up to several PeV. If these neutrinos originate from extra-galactic astrophysical sources, then scattering in transit with dark matter particles will delay their arrival to Earth. This results in a cut-off in their spectrum at an energy set by the scattering cross section, allowing us to place an upper limit on cross sections $\sigma$ which increase with energy E at the level of $\sigma$ < 10^{-17} x (m / GeV) x (E / PeV)^2 cm^2, for dark matter particles of mass m. Once these neutrinos enter our galaxy, the large dark matter densities result in further scattering, especially towards the Galactic Centre. Intriguingly, we find that for $\sigma$ ~ 10^{-22} x (m / GeV) x (E / PeV)^2 cm^2, the distribution of the neutrinos on the sky has a small cluster of events towards the centre of the galaxy, potentially explaining the ~2 sigma excess seen by IceCube in this region without needing a galactic source.
We present a new analysis of the ionizing emissivity ($\dot{N}_{\rm{ion}}$, s$^{-1}$ Mpc$^{-3}$) for galaxies during the epoch of reionization and their potential for completing and maintaining reionization. We use extensive SED modelling -- incorporating two plausible mechanisms for the escape of Lyman continuum photon -- to explore the range and evolution of ionizing efficiencies consistent with new results on galaxy colours ($\beta$) during this epoch. We estimate $\dot{N}_{\rm{ion}}$ for the latest observations of the luminosity and star-formation rate density at $z<10$, outlining the range of emissivity histories consistent with our new model. Given the growing observational evidence for a UV colour-magnitude relation in high-redshift galaxies, we find that for any plausible evolution in galaxy properties, red (brighter) galaxies are less efficient at producing ionizing photons than their blue (fainter) counterparts. The assumption of a redshift and luminosity evolution in $\beta$ leads to two important conclusions. Firstly, the ionizing efficiency of galaxies naturally increases with redshift. Secondly, for a luminosity dependent ionizing efficiency, we find that galaxies down to a rest-frame magnitude of $M_{\rm{UV}} \approx -15$ alone can potentially produce sufficient numbers of ionizing photons to maintain reionization as early as $z\sim8$ for a clumping factor of $C_{\rm{H {\small II}}} \leq 3$.
We study collisions between nearly planar domain walls including the effects of small initial nonplanar fluctuations. These perturbations represent the small fluctuations that must exist in a quantum treatment of the problem. In a previous paper, we demonstrated that at the linear level a subset of these fluctuations experience parametric amplification as a result of their coupling to the planar symmetric background. Here we study the full three-dimensional nonlinear dynamics using lattice simulations, including both the early time regime when the fluctuations are well described by linear perturbation theory as well as the subsequent stage of fully nonlinear evolution. We find that the nonplanar fluctuations have a dramatic effect on the overall evolution of the system. Specifically, once these fluctuations begin to interact nonlinearly the split into a planar symmetric part of the field and the nonplanar fluctuations loses its utility. At this point the colliding domain walls dissolve, with the endpoint of this being the creation of a population of oscillons in the collision region. The original (nearly) planar symmetry has been completely destroyed at this point and an accurate study of the system requires the full three-dimensional simulation.
We present deep 10h VLT/XSHOOTER spectroscopy for an extraordinarily luminous and extended Lya emitter at z=6.595 referred to as Himiko and first discussed by Ouchi et al. (2009), with the purpose of constraining the mechanisms powering its strong emission. Complementary to the spectrum, we discuss NIR imaging data from the CANDELS survey. We find neither for HeII nor any metal line a significant excess, with 3 sigma upper limits of 6.8, 3.1, and 5.8x10^{-18} erg/s/cm^2 for CIV $\lambda$1549, HeII $\lambda$1640, CIII] $\lambda$1909, respectively, assuming apertures with 200 km/s widths and offset by -250 km/s w.r.t to the peak Lya redshift. These limits provide strong evidence that an AGN is not a major contribution to Himiko's Lya flux. Strong conclusions about the presence of PopIII star-formation or gravitational cooling radiation are not possible based on the obtained HeII upper limit. Our Lya spectrum confirms both spatial extent and flux (8.8+/-0.5x10^{-17} erg/s/cm^2) of previous measurements. In addition, we can unambiguously exclude any remaining chance of it being a lower redshift interloper by significantly detecting a continuum redwards of Lya, while being undetected bluewards.
In the Bona-Masso formulation, Einstein equations are written as a set of flux conservative first order hyperbolic equations that resemble fluid dynamics equations. Based on this formulation, we construct a lattice Boltzmann model for Numerical Relativity. Our model is validated with well-established tests, showing good agreement with analytical solutions. Furthermore, we show that by increasing the relaxation time, we gain stability at the cost of losing accuracy, and by decreasing the lattice spacings while keeping a constant numerical diffusivity, the accuracy and stability of our simulations improves. Finally, in order to show the potential of our approach a linear scaling law for parallelisation with respect to number of CPU cores is demonstrated. Our model represents the first step in using lattice kinetic theory to solve gravitational problems.
The N-independence observed in the evolution of massive black hole binaries (MBHBs) in recent simulation of merging stellar bulges suggests a simple interpretation beyond complex time-dependent relaxation processes. We conjecture that the MBHB hardening rate is equivalent to that of a binary immersed in a field of unbound stars with density $\rho$ and typical velocity $\sigma$, provided that $\rho$ and $\sigma$ are the stellar density and the velocity dispersion at the influence radius of the MBHB. By comparing direct N-body simulations to an hybrid model based on 3-body scattering experiments, we verify this hypothesis: when normalized to the stellar density and velocity dispersion at the binary influence radius, the N-body MBHB hardening rate approximately matches that predicted by 3-body scatterings in the investigated cases. The eccentricity evolution obtained with the two techniques is also in reasonable agreement. This result is particularly practical because it allows to estimate the lifetime of MBHBs forming in dry mergers based solely on the stellar density profile of the host galaxy. We briefly discuss some implications of our finding for the gravitational wave signal observable by pulsar timing arrays and for the expected population of MBHBs lurking in massive ellipticals.
Microvariability consists in small time scale variations of low amplitude in the photometric light curves of quasars, and represents an important tool to investigate their inner core. Detection of quasar microvariations is challenging for their non-periodicity, as well as the need for high monitoring frequency and high signal-to-noise ratio. Statistical tests developed for the analysis of quasar differential light curves usually show either low power or low reliability, or both. In this paper we compare two statistical procedures that include several stars to perform tests with enhanced power and high reliability. We perform light curve simulations of variable quasars and non-variable stars, and analyze them with statistical procedures developed from the F-test and the analysis of variance. The results show a large improvement in the power of both statistical probes, and a larger reliability, when several stars are included in the analysis. The results from the simulations agree with those obtained from observations of real quasars. The high power and high reliability of the tests discussed in this paper improve the results that can be obtained from short and long time scale variability studies. These techniques are not limited to quasar variability; on the contrary, they can be easily implemented to other sources such as variable stars. Their applications to future research and to the analysis of large field photometric monitoring archives can reveal new variable sources.
Links to: arXiv, form interface, find, astro-ph, recent, 1505, contact, help (Access key information)
We study collisions between pairs of bubbles nucleated in an ambient false vacuum. For the first time, we include the effects of small initial (quantum) fluctuations around the instanton profiles describing the most likely initial bubble profile. Past studies of this problem neglect these fluctuations and work under the assumption that the collisions posess an exact SO(2,1) symmetry. We use three-dimensional lattice simulations to demonstrate that for double-well potentials, small initial perturbations to this symmetry can be amplified as the system evolves. Initially the amplification is well-described by linear perturbation theory around the SO(2,1) background, but the onset of strong nonlinearities amongst the fluctuations quickly leads to a drastic breaking of the original SO(2,1) symmetry and the production of oscillons in the collision region. We explore several single-field models, and we find it is hard to both realize inflation inside of a bubble and produce oscillons in a collision. Finally, we extend our results to a simple two-field model. The additional freedom allowed by the second field allows us to construct viable inflationary models that allow oscillon production in collisions. The breaking of the SO(2,1) symmetry allows for a new class of observational signatures from bubble collisions that do not posess azimuthal symmetry, including the production of gravitational waves which cannot be supported by an SO(2,1) spacetime.
We study the characteristics of the galaxy cluster samples expected from the European Space Agency's Euclid satellite and forecast constraints on cosmological parameters describing a variety of cosmological models. The method used in this paper, based on the Fisher Matrix approach, is the same one used to provide the constraints presented in the Euclid Red Book (Laureijs et al.2011). We describe the analytical approach to compute the selection function of the photometric and spectroscopic cluster surveys. Based on the photometric selection function, we forecast the constraints on a number of cosmological parameter sets corresponding to different extensions of the standard LambdaCDM model. The dynamical evolution of dark energy will be constrained to Delta w_0=0.03 and Delta w_a=0.2 with free curvature Omega_k, resulting in a (w_0,w_a) Figure of Merit (FoM) of 291. Including the Planck CMB covariance matrix improves the constraints to Delta w_0=0.02, Delta w_a=0.07 and a FoM=802. The amplitude of primordial non-Gaussianity, parametrised by f_NL, will be constrained to \Delta f_NL ~ 6.6 for the local shape scenario, from Euclid clusters alone. Using only Euclid clusters, the growth factor parameter \gamma, which signals deviations from GR, will be constrained to Delta \gamma=0.02, and the neutrino density parameter to Delta Omega_\nu=0.0013 (or Delta \sum m_\nu=0.01). We emphasise that knowledge of the observable--mass scaling relation will be crucial to constrain cosmological parameters from a cluster catalogue. The Euclid mission will have a clear advantage in this respect, thanks to its imaging and spectroscopic capabilities that will enable internal mass calibration from weak lensing and the dynamics of cluster galaxies. This information will be further complemented by wide-area multi-wavelength external cluster surveys that will already be available when Euclid flies. [Abridged]
We examine the cosmological correlators induced by the simultaneous breaking of parity and of statistical isotropy, e.g., in presence of the coupling ${\cal L} = f(\phi) ( - \frac{1}{4} F^2 + \frac{\gamma}{4} F \tilde{F} )$ between the inflaton $\phi$ and a vector field with vacuum expectation value ${\bf A}$. For a suitably chosen function $f$, the energy in the vector field $\rho_{\rm A}$ does not decay during inflation. This results in nearly scale-invariant signatures of broken statistical isotropy and parity. Specifically, we find that the scalar-scalar correlator of primordial curvature perturbations includes a quadrupolar anisotropy, $P_\zeta ( {\bf k}) = P(k)[ 1 + g_* ( \hat{\bf k} \cdot \hat{\bf A})^2]$, and a (angle-averaged) scalar bispectrum that is a linear combination of the first $3$ Legendre polynomials, $B_\zeta(k_1, k_2, k_3) = \sum_L c_L P_L (\hat{\bf k}_1 \cdot \hat{\bf k}_2) P(k_1) P(k_2) + 2~{\rm perms} $, with $c_0 : c_1 : c_2 = 2 : -3 : 1$ ($c_1 \neq 0$ is a consequence of parity violation, corresponding to the constant $\gamma \neq 0$). The latter is one of the main results of this paper, which provides for the first time a clear example of an inflationary model where a non-negligible $c_1$ contribution to the bispectrum is generated. The scalar-tensor and tensor-tensor correlators induce characteristic signatures in the Cosmic Microwave Background temperature anisotropies (T) and polarization (E/B modes); namely, non-diagonal contributions to $\langle a_{\ell_1 m_1} a_{\ell_2 m_2}^* \rangle$, with $|\ell_1 - \ell_2| = 1$ in TT, TE, EE and BB, and $|\ell_1 - \ell_2| = 2$ in TB and EB. The latest CMB bounds on the scalar observables ($g_*$, $c_0$, $c_1$ and $c_2$), translate into the upper limit $\rho_{\rm A} / \rho_\phi \lesssim 10^{-9}$ at $\gamma=0$. We find that the upper limit on the vector energy density becomes much more stringent as $\gamma$ grows.
We use astrophysical and atomic clock tests of the stability of the fine-structure constant $\alpha$, together with Type Ia supernova and Hubble parameter data, to constrain the simplest class of dynamical dark energy models where the same degree of freedom is assumed to provide both the dark energy and (through a dimensionless coupling, $\zeta$, to the electromagnetic sector) the $\alpha$ variation. We show how current data tightly constrains a combination of $\zeta$ and the dark energy equation of state $w_0$. At the $95\%$ confidence level and marginalizing over $w_0$ we find $|\zeta|<5\times10^{-6}$, with the atomic clock tests dominating the constraints. The forthcoming generation of high-resolution ultra-stable spectrographs will enable significantly tighter constraints.
We study the possibility of realizing the growth rate of matter density perturbations lower than that in General Relativity. Using the approach of the effective field theory of modified gravity encompassing theories beyond Horndeski, we derive the effective gravitational coupling $G_{\rm eff}$ and the gravitational slip parameter $\eta$ for the perturbations deep inside the Hubble radius. In Horndeski theories the necessary condition for achieving weak gravity is that the tensor propagation speed $c_{\rm t}$ is smaller than the speed of light, but this is not a sufficient condition due to the presence of a scalar-matter interaction that always enhances $G_{\rm eff}$. Beyond the Horndeski domain it is possible to realize $G_{\rm eff}$ smaller than the Newton gravitational constant $G$, while the scalar and tensor perturbations satisfy no-ghost and stability conditions. We present a concrete dark energy scenario with varying $c_{\rm t}$ and numerically study the evolution of perturbations to confront the model with the observations of red-shift space distortions and weak lensing.
Usually the effects of isotropic inhomogeneities are not seriously taken into account in the determination of the cosmological parameters because of Copernican principle whose statement is that we do not live in the privileged domain in the universe. But Copernican principle has not been observationally confirmed yet in sufficient accuracy, and there is the possibility that there are non-negligible large-scale isotropic inhomogeneities in our universe. In this paper, we study the effects of the isotropic inhomogeneities on the determination of the cosmological parameters and show the probability that non-Copernican isotropic inhomogeneities mislead us into believing, for example, the phantom energy of the equation of state, $p=w\rho$ with $w<-1$, even in case that $w=-1$ is the true value.
We present a manifestly local, diffeomorphism invariant and locally Poincare invariant formulation of vacuum energy sequestering. In this theory, quantum vacuum energy generated by matter loops is cancelled by auxiliary fields. The auxiliary fields decouple from gravity almost completely. Their only residual effect is an {\it a priori} arbitrary, finite contribution to the curvature of the background geometry, which is radiatively stable. Its value is to be determined by a measurement, like the finite part of any radiatively stable UV-sensitive quantity in quantum field theory
Percolation theory is applied to the phase transition dynamics of domain pattern formation in segregating quasi-two-dimensional binary Bose--Einstein condensates. Our numerical experiments revealed that the percolation threshold is close to 0.5. A long-range open domain wall appears with a fractal dimension between two percolating domains. Such a wall can survive for a long time as a relic of the phase transition according to the dynamic finite-size-scaling hypothesis, which seems to be in contrast to the current understanding in cosmology that an infinite defect violates a scale invariance.
Wavelength-dependent point spread functions (PSFs) violate an implicit assumption in current galaxy shape measurement algorithms that deconvolve the PSF measured from stars (which have stellar spectral energy distributions (SEDs)) from images of galaxies (which have galactic SEDs). Since the absorption length of silicon depends on wavelength, CCDs are a potential source of PSF chromaticity. Here we develop two toy models to estimate the sensitivity of the cosmic shear survey from the Large Synoptic Survey Telescope to chromatic effects in CCDs. We then compare these toy models to simulated estimates of PSF chromaticity derived from the LSST photon simulator PhoSim. We find that even though sensor contributions to PSF chromaticity are subdominant to atmospheric contributions, they can still significantly bias cosmic shear results if left uncorrected, particularly in the redder filter bands and for objects that are off-axis in the field of view.
We study the proposal by Mersini et al. that the observed dark energy might be explained by the back-reaction of the set of tail modes in a theory with a dispersion relation in which the mode frequency decays exponentially in the trans-Planckian regime. The matter tail modes are frozen out, however they induce metric fluctuations. The energy-momentum tensor with which the tail modes effect the background geometry obtains contributions from both metric and matter fluctuations. We calculate the equation of state induced by the tail modes taking into account the gravitational contribution. We find that, in contrast to the case of frozen super-Hubble cosmological fluctuations, in this case the matter perturbations dominate, and they yield an equation of state which to leading order takes the form of a positive cosmological constant.
We review the status of String Gas Cosmology after the 2015 Planck data release. String gas cosmology predicts an almost scale-invariant spectrum of cosmological perturbations with a slight red tilt, like the simplest inflationary models. It also predicts a scale-invariant spectrum of gravitational waves with a slight blue tilt, unlike inflationary models which predict a red tilt of the gravitational wave spectrum. String gas cosmology yields two consistency relations which determine the tensor to scalar ratio and the slope of the gravitational wave spectrum given the amplitude and tilt of the scalar spectrum. We show that these consistency relations are in good agreement with the Planck data. We discuss future observations which will be able to differentiate between the predictions of inflation and those of string gas cosmology.
During inflation, scalar fields, including the Higgs boson, may acquire a nonzero vacuum expectation value, which must later relax to the equilibrium value during reheating. In the presence of the time-dependent condensate, the vacuum state can evolve into a state with a nonzero particle number. We show that, in the presence of lepton number violation in the neutrino sector, the particle production can explain the observed matter-antimatter asymmetry of the universe. We find that this form of leptogenesis is particularly effective when the Higgs condensate decays rapidly and at low reheat temperature. As part of the calculation, we present some exact results for the Bogoliubov transformations for Majorana fermions with a nonzero time-dependent chemical potential, in addition to a time-dependent mass.
Detection of the global redshifted 21~cm signal is an excellent means of deciphering the physical processes during the Dark Ages and subsequent Epoch of Reionization (EoR). However, the detection of this faint signal is challenging due to the high precision required in instrumental calibration and modeling of substantially brighter foregrounds and instrumental systematics. In particular, modeling and removal of receiver noise with mK accuracy remains a formidable task in experiments aiming to detect the global signal using single-element spectral radiometers. Interferometers do not respond to receiver noise; therefore, we explore here the theory of the response of interferometers to global signals. We first derive the response to uniform sky of interferometers made of element antennas, then extend the analysis to interferometers made of 1-D arrays and finally consider 2-D aperture antennas. The analysis suggests that short-spacing interferometers made of omnidirectional antennas have the best sensitivity for the detection of the global signal. These antennas may be wideband 1-D arrays of dipoles or 1-D aperture antennas; we argue that the interferometer is best configured to be EW with the elements oriented NS. However, the performance of such interferometers would be limited by crosstalk, mode-coupling of foreground continuum sources into spectral confusion, and uncertainty in its telescope filter function. We conclude that the only useful interferometer for the global EoR signal is an interferometer between two collinear wideband short dipoles, with sensitivity to the global signal realized by placing a space beam splitter between the elements to form a Zero-Spacing Interferometer.
It is known that after Affleck-Dine baryogenesis, spatial inhomogeneities of Affleck-Dine field grow into non-topological solitons called Q-balls. In gauge mediated SUSY breaking models, sufficiently large Q-balls with baryon charge are stable while Q-balls with lepton charge can always decay into leptons. For a Q-ball that carries nonzero $B$ and $L$ charges, the difference between the baryonic component and the leptonic component in decay rate may induce nonzero electric charge on the Q-ball. This implies that charged Q-ball, also called gauged Q-ball, may emerge in our universe. In this paper, we investigate two complex scalar fields, a baryonic scalar field and a leptonic one, in an Abelian gauge theory. We find stable solutions of gauged Q-balls for different baryon and lepton charges. Those solutions shows that a Coulomb potential arises and the Q-ball becomes electrically charged as expected. It is energetically favored that some amount of leptonic component decays, but there is an upper bound on its amount due to the Coulomb force. The baryonic decay also becomes possible by virtue of electrical repulsion and we find the condition to suppress it so that the charged Q-balls can survive in the universe.
In de Rham-Gabadadze-Tolley (dRGT) massive gravity and bigravity, a non-minimal matter coupling involving both metrics generically re-introduces the Boulware--Deser (BD) ghost. A non-minimal matter coupling via a simple, yet specific composite metric has been proposed, which eliminates the BD ghost below the strong coupling scale. Working explicitly in the metric formulation and for arbitrary spacetime dimensions, we show that this composite metric is the unique consistent non-minimal matter coupling below the strong coupling scale, which emerges out of two diagnostics, namely, absence of Ostrogradski ghosts in the decoupling limit and absence of the BD ghost from matter quantum loop corrections.
A large set of independent astronomical observations have provided a strong evidence for nonbaryonic dark matter in the Universe. One of the most investigated candidates is an unknown long-lived Weakly Interacting Massive Particle (WIMP) which was in thermal equilibrium with the primeval plasma. Here we investigate the WIMP abundance based on the relativistic kinetic treatment for gravitationally induced particle production recently proposed in the literature (Lima \& Baranov, Phys. Rev. D {\bf 90}, 043515, 2014). The new evolution equation is deduced and solved both numerically and also through a semi-analytical approach. The predictions of the WIMP observables are discussed and compared with the ones obtained in the standard approach.
Links to: arXiv, form interface, find, astro-ph, recent, 1505, contact, help (Access key information)
Among the possible alternatives to the standard cosmological model ($\Lambda$CDM), coupled Dark Energy models postulate that Dark Energy (DE), seen as a dynamical scalar field, may interact with Dark Matter (DM), giving rise to a "fifth-force", felt by DM particles only. In this paper, we study the impact of these cosmologies on the statistical properties of galaxy populations by combining high-resolution numerical simulations with semi-analytic models (SAM) of galaxy formation and evolution. New features have been implemented in the reference SAM in order to have it run self-consistently and calibrated on these cosmological simulations. They include an appropriate modification of the mass temperature relation and of the baryon fraction in DM haloes, due to the different virial scalings and to the gravitational bias, respectively. Our results show that the predictions of our coupled-DE SAM do not differ significantly from theoretical predictions obtained with standard SAMs applied to a reference $\Lambda$CDM simulation, implying that the statistical properties of galaxies provide only a weak probe for these alternative cosmological models. On the other hand, we show that both galaxy bias and the galaxy pairwise velocity distribution are sensitive to coupled DE models: this implies that these probes might be successfully applied to disentangle among quintessence, $f(R)$-Gravity and coupled DE models.
Looking at the region connecting two clusters is a promising way to identify and study the Warm-Hot Intergalactic Medium. Observations show that the spectrum of the bridge between A3556 and A3558 has a stronger soft X-ray emission than the nearby region. Suzaku observations could not discriminate the origin of the extra emission. In this work we analyze a dedicated Chandra observation of the same target to identify point sources and characterize the background emission in the bridge. We find that the count number of the point sources is much higher than average field population (using CDFS~4~Ms as a reference). Moreover, the shape of the cumulative distribution resembles that of galaxy distribution suggesting that the point sources are galaxies in a filament. The Suzaku extra emission is well explained by the high abundance of point sources identified by Chandra. Furthermore, we used optical/IR observations of point sources in the same field to estimate the density of the putative filament as rho~150 rho_b$, below Suzaku sensitivity.
Vector mode of cosmological perturbation theory imprints characteristic signals on the weak lensing signals such as curl- and B-modes which are never imprinted by the scalar mode. However, the vector mode is neglected in the standard first-order cosmological perturbation theory since it only has a decaying mode. This situation changes if the cosmological perturbation theory is expanded up to second order. The second-order vector and tensor modes are inevitably induced by the product of the first-order scalar modes. We study the effect of the second-order vector mode on the weak lensing curl- and B-modes. The curl-mode induced by the second-order vector mode dominates instead of the primordial gravitational waves when the tensor-to-scalar ratio is $r = 0.1$ and the second-order tensor mode at $\ell \geq 200$. Furthermore, the B-mode cosmic shear induced by the second-order vector mode dominates on almost all scales. However, we find that the observational signatures of the second-order vector and tensor modes cannot exceed the expected noise of ongoing and upcoming weak lensing measurements. We conclude that the curl- and B-modes induced by the second-order vector and tensor modes are unlikely to be detected in future experiments.
The mapping between the distributions of the observed galaxy stellar mass and the underlying dark matter halos provides the crucial link from theories of large-scale structure formation to interpreting the complex phenomena of galaxy formation and evolution. We develop a novel statistical method, based on the Halo Occupation Distribution model (HOD), to solve for this mapping by jointly fitting the galaxy clustering and the galaxy-galaxy lensing measured from the Sloan Digital Sky Survey (SDSS). The method, called the iHOD model, extracts maximum information from the survey by including ~80% more galaxies than the traditional HOD methods, and takes into account the incompleteness of the stellar mass samples in a statistically consistent manner. The derived stellar-to-halo mass relation not only explains the clustering and lensing of SDSS galaxies over almost four decades in stellar mass, but also successfully predicts the stellar mass functions observed in SDSS. Due to its capability of modelling significantly more galaxies, the iHOD is able to break the degeneracy between the logarithmic scatter in the stellar mass at fixed halo mass and the slope of the stellar-to-halo mass relation at high mass end, without the need to assume a strong prior on the scatter and/or use the stellar mass function as an input. We detect a decline of the scatter with halo mass, from 0.22 dex at below 10^{12} Msun to 0.18 dex at 10^{14} Msun. The model also enables stringent constraints on the satellite stellar mass functions at fixed halo mass, predicting a departure from the Schechter functional form in high mass halos. The iHOD model can be easily applied to existing and future spectroscopic datasets, greatly improving the statistical constraint on the stellar-to-halo mass relation compared to the traditional HOD methods within the same survey.
The Ly$\alpha$ line of high-redshift galaxies has emerged as a powerful probe of both early galaxy evolution and the epoch of reionization (EoR). Motivated by the upcoming wide-field survey with the Subaru Hyper Supreme-Cam (HSC), we study the angular correlation function (ACF) of narrow-band selected, $z\approx7$ LAEs. The clustering of LAEs is determined by both: (i) their typical host halo masses, $\bar{M}_{\rm h}$; (ii) the absorption due to a patchy EoR, characterized by an average neutral fraction of the IGM, $\bar{x}_{\rm HI}$. We bracket the allowed LAE ACF by exploring extreme scenarios for both the intrinsic Ly$\alpha$ emission and the morphology of cosmic ionized patches in physical EoR models. Current LAE ACF measurements imply that the Universe is mostly ionized at $z\approx7$, with $\bar{x}_{\rm HI}\lesssim0.5$ (1-$\sigma$) even for an extremely conservative model of intrinsic emission. The upcoming Ultra Deep campaign with the HSC could improve on these constraints by tens of percent, or $\bar{x}_{\rm HI}\lesssim0.3$ if the mean value of the ACF remains unchanged. The ACF at a fixed observed LAE number density and $\bar{x}_{\rm HI}$ is extremely insensitive to the EoR morphology; distinguishing between different EoR models would therefore require more accurate redshift determinations with spectroscopic follow-up observations. We also find that the low values of the currently-observed ACF implies that LAEs are hosted by relatively small dark matter halos, with $\bar{M}_{\rm h}\lesssim10^{10}M_\odot$. Combined with their observed number densities, this implies a duty cycle $\lesssim$ few per cent. These values are over an order of magnitude lower than the analogous ones for color-selected, Lyman break galaxies. This discrepancy could be due to the narrow-band LAEs searches preferentially selecting a population of young, star-burst galaxies, residing in less massive halos.
We present a mass estimate of the Planck-discovered cluster PLCK G100.2-30.4, derived from a weak lensing analysis of deep SUBARU griz images. We perform a careful selection of the background galaxies using the multi-band imaging data, and undertake the weak lensing analysis on the deep (1hr) r-band image. The shape measurement is based on the KSB algorithm; we adopt the PSFex software to model the Point Spread Function (PSF) across the field and correct for this in the shape measurement. The weak lensing analysis is validated through extensive image simulations. We compare the resulting weak lensing mass profile and total mass estimate to those obtained from our re-analysis of XMM-Newton observations, derived under the hypothesis of hydrostatic equilibrium. The total integrated mass profiles are in remarkably good agreement, agreeing within 1$\sigma$ across their common radial range. A mass $M_{500} \sim 7 x 10^{14} M_\odot$ is derived for the cluster from our weak lensing analysis. Comparing this value to that obtained from our reanalysis of XMM-Newton data, we obtain a bias factor of (1-b) = 0.8 $\pm$ 0.1. This is compatible within 1$\sigma$ with the value of (1-b) obtained by Planck Collaboration XXIV from their calibration of the bias factor using newly-available weak lensing reconstructed masses.
The $m$-$z$ relation for type Ia supernovae is one of the key pieces of evidence supporting the cosmological `concordance model' with $\lambda_0 \approx 0.7$ and $\Omega_0 \approx 0.3$. However, it is well known that the $m$-$z$ relation depends not only on $\lambda_0$ and $\Omega_0$ (with $H_0$ as a scale factor) but also on the density of matter along the line of sight, which is not necessarily the same as the large-scale density. I investigate to what extent the measurement of $\lambda_0$ and $\Omega_0$ depends on this density when it is characterized by the parameter $\eta$ ($0 \le \eta \le 1$), which describes the ratio of density along the line of sight to the overall density. I also discuss what constraints can be placed on $\eta$, both with and without constraints on $\lambda_0$ and $\Omega_0$ in addition to those from the $m$-$z$ relation for type~Ia supernovae.
We investigate the effects of homogeneous general dark energy on the non-linear matter perturbation in fully general relativistic context. Taking into account the next-to-leading corrections, the total power spectrum with general dark energy deviates from the LambdaCDM spectrum, which is nearly identical to that in the Einstein-de Sitter universe, as large as a few percent at scales comparable to that for the baryon acoustic oscillations and increases on smaller scales. The contribution from the curvature perturbation, which is absent in the Newtonian theory, exhibits even more drastic difference larger than 100%, while the amplitude is heavily suppressed on all scales.
We compare the shapes and intrinsic alignments of galaxies in the MassiveBlack-II cosmological hydrodynamic simulation (MBII) to those in a dark matter-only (DMO) simulation performed with the same volume (100$h^{-1}$Mpc)$^{3}$, cosmological parameters, and initial conditions. Understanding the impact of baryonic physics on galaxy shapes and alignments and their relation to the dark matter distribution should prove useful to map the intrinsic alignments of galaxies from hydrodynamic to dark matter-only simulations. We find that dark matter subhalos are typically rounder in MBII, and the shapes of stellar matter in low mass galaxies are more misaligned with the shapes of the dark matter of the corresponding subhalos in the DMO simulation. At $z=0.06$, the fractional difference in the mean misalignment angle between MBII and DMO simulations varies from $\sim 28 \% - 12 \%$ in the mass range $10^{10.8} - 6.0 \times 10^{14} h^{-1}M_{\odot}$. We study the dark matter halo shapes and alignments as a function of radius, and find that while galaxies in MBII are more aligned with the inner parts of their dark matter subhalos, there is no radial trend in their alignments with the corresponding subhalo in the DMO simulation. This result highlights the importance of baryonic physics in determining the alignment of the galaxy with respect to the inner parts of the halo. Finally, we compare the ellipticity-direction (ED) correlation for galaxies to that for dark matter halos, finding that it is suppressed on all scales by stellar-dark matter misalignment. In the projected shape-density correlation ($w_{\delta+}$), which includes ellipticity weighting, this effect is partially canceled by the higher mean ellipticities of the stellar component, but differences of order $30-40\%$ remain on scales $> 1$ Mpc over a range of subhalo masses, with scale-dependent effects below $1$ Mpc.
With the goal of deriving dissipative hydrodynamics from an action, we study classical actions for open systems, which follow from the generic structure of effective actions in the Schwinger-Keldysh Closed-Time-Path formalism with two time axes and a doubling of degrees of freedom. The central structural feature of such effective actions is the coupling between degrees of freedom on the two time axes. This reflects the fact that from an effective field theory point of view, dissipation is the loss of energy of the low-energy hydrodynamical degrees of freedom to the integrated-out, UV degrees of freedom of the environment. The dynamics of only the hydrodynamical modes may therefore not posses a conserved stress-energy tensor. After a general discussion of the CTP effective actions, we use the variational principle to derive the energy-momentum balance equation for a dissipative fluid from an effective Goldstone action of the long-range hydrodynamical modes. Despite the absence of conserved energy and momentum, we show that we can construct the first-order dissipative stress-energy tensor and derive the Navier-Stokes equations near hydrodynamical equilibrium. The shear viscosity is shown to vanish in the classical theory under consideration, while the bulk viscosity is determined by the form of the effective action. We also discuss the thermodynamics of the system and analyse the entropy production.
Cosmology is making impressive progress and it is producing stringent bounds on the sum of the neutrino masses {\Sigma}, a parameter of great importance for the current laboratory experiments. In this letter, we exploit the potential relevance of the analysis of Palanque-Delabrouille et al. [JCAP 1502, 045 (2015)] to the neutrinoless double beta decay (0\nu\beta\beta) search. In fact, this result seems to favor the normal hierarchy spectrum for the light neutrino masses. If this is confirmed, the impact on the 0\nu\beta{\beta} experiments will be tremendous, bringing the possibility of detecting a signal out of the reach of the next generation of experiments.
It has been suggested that several small-scale structure anomalies in $\Lambda$CDM cosmology can be solved by strong self-interaction between dark matter particles. It was shown by Braaten and Hammer that the presence of a near threshold S-wave resonance can make the scattering cross section at nonrelativistic speeds come close to saturating the unitarity bound. This can result in the formation of a stable bound state of two asymmetric dark matter particles (which we call darkonium). Laha and Braaten studied the nuclear recoil energy spectrum in dark matter direct detection experiments due to this incident bound state. Here we study the angular recoil spectrum, and show that it is uniquely determined up to normalization by the S-wave scattering length. Observing this angular recoil spectrum in a dark matter directional detection experiment will uniquely determine many of the low-energy properties of dark matter independent of the underlying dark matter microphysics.
The quest for a $B$-mode imprint from primordial gravity waves on the polarization of the cosmic microwave background (CMB) requires the characterization of foreground polarization from Galactic dust. We present a statistical study of the filamentary structure of the $353\,$GHz Planck Stokes maps at high Galactic latitude, relevant to the study of dust emission as a polarized foreground to the CMB. We filter the intensity and polarization maps to isolate filaments in the range of angular scales where the power asymmetry between $E$-modes and $B$-modes is observed. Using the Smoothed Hessian Major Axis Filament Finder, we identify 259 filaments at high Galactic latitude, with lengths larger or equal to $2$\deg\ (corresponding to $3.5\,$pc in length for a typical distance of $100\,$pc). These filaments show a preferred orientation parallel to the magnetic field projected onto the plane of the sky, derived from their polarization angles. We present mean maps of the filaments in Stokes $I$, $Q$, $U$, $E$, and $B$, computed by stacking individual images rotated to align the orientations of the filaments. Combining the stacked images and the histogram of relative orientations, we estimate the mean polarization fraction of the filaments to be $11\,$%. Furthermore, we show that the correlation between the filaments and the magnetic field orientations may account for the $E$ and $B$ asymmetry and the $C_{\ell}^{TE}/C_{\ell}^{EE}$ ratio, reported in the power spectra analysis of the Planck $353\,$GHz polarization maps. Future models of the dust foreground for CMB polarization studies will need to take into account the observed correlation between the dust polarization and the structure of interstellar matter.
Recently, we developed a pair of meshless finite-volume Lagrangian methods for hydrodynamics: the 'meshless finite mass' (MFM) and 'meshless finite volume' (MFV) methods. These capture advantages of both smoothed-particle hydrodynamics (SPH) and adaptive mesh-refinement (AMR) schemes. Here, we extend these to include ideal magneto-hydrodynamics (MHD). The MHD equations are second-order consistent and conservative. We augment these with a divergence-cleaning scheme, which maintains div*B~0 to high accuracy. We implement these in the code GIZMO, together with a state-of-the-art implementation of SPH MHD. In every one of a large suite of test problems, the new methods are competitive with moving-mesh and AMR schemes using constrained transport (CT) to ensure div*B=0. They are able to correctly capture the growth and structure of the magneto-rotational instability (MRI), MHD turbulence, and the launching of magnetic jets, in some cases converging more rapidly than AMR codes. Compared to SPH, the MFM/MFV methods exhibit proper convergence at fixed neighbor number, sharper shock capturing, and dramatically reduced noise, div*B errors, and diffusion. Still, 'modern' SPH is able to handle most of our tests, at the cost of much larger kernels and 'by hand' adjustment of artificial diffusion parameters. Compared to AMR, the new meshless methods exhibit enhanced 'grid noise' but reduced advection errors and numerical diffusion, velocity-independent errors, and superior angular momentum conservation and coupling to N-body gravity solvers. As a result they converge more slowly on some problems (involving smooth, slowly-moving flows) but more rapidly on others (involving advection or rotation). In all cases, divergence-control beyond the popular Powell 8-wave approach is necessary, or else all methods we consider will systematically converge to unphysical solutions.
Sterile neutrinos with mass ~1 eV and order 10% mixing with active neutrinos have been proposed as a solution to anomalies in neutrino oscillation data, but are tightly constrained by cosmological limits. It was recently shown that these constraints are avoided if sterile neutrinos couple to a new MeV-scale gauge boson A'. However, even this scenario is restricted by structure formation constraints when A'-mediated collisional processes lead to efficient active-to-sterile neutrino conversion after neutrinos have decoupled. In view of this, we reevaluate in this paper the viability of sterile neutrinos with such "secret" interactions. We carefully dissect their evolution in the early Universe, including the various production channels and the expected modifications to large scale structure formation. We argue that there are two regions in parameter space - one at very small A' coupling, one at relatively large A' coupling - where all constraints from big bang nucleosynthesis (BBN), cosmic microwave background (CMB), and large scale structure (LSS) data are satisfied. Interestingly, the large A' coupling region is precisely the region that was previously shown to have potentially important consequences for the small scale structure of dark matter halos if the A' boson couples also to the dark matter in the Universe.
LCDM is remarkably successful in predicting the cosmic microwave background and large-scale structure, and LCDM parameters have been determined with only mild tensions between different types of observations. Hydrodynamical simulations starting from cosmological initial conditions are increasingly able to capture the complex interactions between dark matter and baryonic matter in galaxy formation. Simulations with relatively low resolution now succeed in describing the overall galaxy population. For example, the EAGLE simulation in volumes up to 100 cubic Mpc reproduces the observed local galaxy mass function nearly as well as semi-analytic models. It once seemed that galaxies are pretty smooth, that they generally grow in size as they evolve, and that they are a combination of disks and spheroids. But recent HST observations combined with high-resolution hydrodynamic simulations are showing that most star-forming galaxies are very clumpy; that galaxies often undergo compaction which reduces their radius and increases their central density; and that most lower-mass star-forming galaxies are not spheroids or disks but are instead elongated when their centers are dominated by dark matter. We also review LCDM challenges on smaller scales: cusp-core, "too big to fail," and substructure issues. Although starbursts can rapidly drive gas out of galaxy centers and thereby reduce the dark matter density, it remains to be seen whether this or other baryonic physics can explain the observed rotation curves of the entire population of dwarf and low surface brightness galaxies. If not, perhaps more complicated physics such as self-interacting dark matter may be needed. But standard LCDM appears to be successful in predicting the dark matter halo substructure that is now observed via gravitational lensing and breaks in cold stellar streams, and any alternative theory must do at least as well.
In this paper, we investigate the influences of the radiation pressures of the central engines on the black hole virial masses $M_{\rm{RM}}$ for 40 active galactic nuclei (AGNs) with high accretion rates. The radiation pressures on the clouds in broad-line regions are equivalent to decrease the gravitational forces of the central black holes, i.e., the clouds undergo the decreased gravitational forces of the black holes due to the central radiation. For the AGNs with the high accretion rates, a part of the gravitational forces of the black holes is counteracted by the central radiation, and the counteracted fraction depends on the percent of ionized hydrogen in clouds (or the ionized depth ratio of clouds). The black hole masses counteracted by the radiation pressures $M_{\rm{RP}}$ are not negligible compared to, or are comparable to $M_{\rm{RM}}$ at least for a part of the AGNs. The black hole masses $M_{\rm{\bullet}}$ are underestimated at least by a factor of 30--50 percent for the AGNs with the close- and super-Eddington limit, regardless of redshifts of sources. It is more appropriate to ignore the radiation pressures of the central sources for the AGNs with lower accretion rates, but the influences of the radiation pressures shall be taken into account to estimate the black hole masses of the extremely high accretion rate AGNs based on the reverberation mapping method.
Modified Newtonian dynamics (MOND) provides a paradigm alternative to dark matter that has been successful in fitting and predicting the rich phenomenology of rotating disc galaxies. There have also been attempts to test MOND in dispersion-supported early-type galaxies, but it remains unclear whether MOND can fit the various empirical properties of early-type galaxies. As a way of rigorously testing MOND in elliptical galaxies we calculate the MOND-predicted velocity dispersion profiles (VDPs) in the inner regions of $\sim 2000$ nearly round SDSS elliptical galaxies under a variety of assumptions on VD anisotropy, and then compare the predicted distribution of VDP slopes with the observed distribution in 11 ATLAS3d galaxies selected with essentially the same criteria. We find that the MOND model parameterised with an interpolating function that works well for rotating galaxies can also reproduce the observed distribution of VDP slopes based only on the observed stellar mass distribution without DM or any other galaxy-to-galaxy varying factor. This is remarkable in view that Newtonian dynamics with DM requires a specific amount and/or profile of DM for each galaxy in order to reproduce the observed distribution of VDP slopes. When we analyse non-round galaxy samples using the MOND-based spherical Jeans equation, we do not find any systematic difference in the mean property of the VDP slope distribution compared with the nearly round sample. However, in line with previous studies of MOND through individual analyses of elliptical galaxies, varying MOND interpolating function or VD anisotropy can lead to systematic change in the VDP slope distribution, indicating that a statistical analysis of VDPs can be used to constrain specific MOND models with an accurate measurement of VDP slopes or a prior constraint on VD anisotropy.
In the scenario of single-field inflation, this field is done in terms of Jacobian elliptic functions. This approach provides, when constrained to particular cases, analytic solutions already known in the past, generalizing them to a bigger family of analytical solutions. The emergent cosmology is analysed using the Hamilton-Jacobi approach and then, the main results are contrasted with the recent measurements obtained from the Planck 2015 data.
The present paper is devoted to find a new generalization of the generalized Chaplygin gas. Therefore, starting from the Hubble parameter associated to the Chaplygin scalar field and using some elliptic identities, the elliptic generalization is straightforward. Thus, all relevant quantities that drive inflation are calculated exactly. Finally, using the measurement on inflation from the Planck 2015 results, observational constraints on the parameters are given.
Beginning with a set of simplified models for spin-0, spin-$\half$, and spin-1 dark matter candidates using completely general Lorentz invariant and renormalizable Lagrangians, we derive the full set of non-relativistic operators and nuclear matrix elements relevant for direct detection of dark matter, and use these to calculate rates and recoil spectra for scattering on various target nuclei. This allows us to explore what high energy physics constraints might be obtainable from direct detection experiments, what degeneracies exist, which operators are ubiquitous and which are unlikely or sub-dominant. We find that there are operators which are common to all spins as well operators which are unique to spin-$\half$ and spin-1 and elucidate two new operators which have not been previously considered. In addition we demonstrate how recoil energy spectra can distinguish fundamental microphysics if multiple target nuclei are used. Our work provides a complete roadmap for taking generic fundamental dark matter theories and calculating rates in direct detection experiments. This provides a useful guide for experimentalists designing experiments and theorists developing new dark matter models.
Links to: arXiv, form interface, find, astro-ph, recent, 1505, contact, help (Access key information)
We explore a chameleon type of interacting dark matter-dark energy scenario in which a scalar field adiabatically traces the minimum of an effective potential sourced by the dark matter density. We discuss extensively the effect of this coupling on cosmological observables, especially the parameter degeneracies expected to arise between the model parameters and other cosmological parameters, and then test the model against observations of the cosmic microwave background (CMB) anisotropies and other cosmological probes. We find that the chameleon parameters $\alpha$ and $\beta$, which determine respectively the slope of the scalar field potential and the dark matter-dark energy coupling strength, can be constrained to $\alpha < 0.17$ and $\beta < 0.19$ using CMB data alone. The latter parameter in particular is constrained only by the late Integrated Sachs-Wolfe effect. Adding measurements of the local Hubble expansion rate $H_0$ tightens the bound on $\alpha$ by a factor of two, although this apparent improvement is arguably an artefact of the tension between the local measurement and the $H_0$ value inferred from Planck data in the minimal $\Lambda$CDM model. The same argument also precludes chameleon models from mimicking a dark radiation component, despite a passing similarity between the two scenarios in that they both delay the epoch of matter-radiation equality. Based on the derived parameter constraints, we discuss possible signatures of the model for ongoing and future large-scale structure surveys.
The quest for the particle nature of dark matter is one of the big open questions of modern physics. The CRESST-II experiment, located at the Gran Sasso laboratory in Italy, is optimised for the detection of the elastic scattering of dark matter particles with ordinary matter. We present the result obtained with an improved detector setup with increased radiopurity and enhanced background rejection. The limit obtained in the so-called low mass region between one and three GeV/c2 is at the present among the best limits obtained for direct dark matter experiments. In addition we give an outlook of the future potential for direct dark matter detection using further improved CRESST CaWO4 cryogenic detectors.
We determine the concentration-mass relation of 19 X-ray selected galaxy clusters from the CLASH survey in theories of gravity that directly modify the lensing potential. We model the clusters as NFW haloes and fit their lensing signal, in the Cubic Galileon and Nonlocal gravity models, to the lensing convergence profiles of the clusters. We discuss a number of important issues that need to be taken into account, associated with the use of nonparametric and parametric lensing methods, as well as assumptions about the background cosmology. Our results show that the concentration and mass estimates in the modified gravity models are, within the errorbars, the same as in $\Lambda$CDM. This result demonstrates that, for the Nonlocal model, the modifications to gravity are too weak at the cluster redshifts, and for the Galileon model, the screening mechanism is very efficient inside the cluster radius. However, at distances $\sim \left[2-20\right] {\rm Mpc}/h$ from the cluster center, we find that the surrounding force profiles are enhanced by $\sim20-40\%$ in the Cubic Galileon model. This has an impact on dynamical mass estimates, which means that tests of gravity based on comparisons between lensing and dynamical masses can also be applied to the Cubic Galileon model.
A model equation for the Reynolds number dependence of the dimensionless dissipation rate in freely decaying homogeneous magnetohydrodynamic turbulence in the absence of a mean magnetic field is derived from the real-space energy balance equation, leading to $C_{\varepsilon}=C_{\varepsilon, \infty}+C/R_- +O(1/R_-^2))$, where $R_-$ is a generalized Reynolds number. The constant $C_{\varepsilon, \infty}$ describes the total energy transfer flux. This flux depends on magnetic and cross helicities, because these affect the nonlinear transfer of energy, suggesting that the value of $C_{\varepsilon,\infty}$ is not universal. Direct numerical simulations were conducted on up to $2048^3$ grid points, showing good agreement between data and the model. The model suggests that the magnitude of cosmological-scale magnetic fields is controlled by the values of the vector field correlations. The ideas introduced here can be used to derive similar model equations for other turbulent systems.
The statistical tests - done by the authors - are surveyed, which verify the null-hypothesis of the intrinsic randomness in the angular distribution of gamma-ray bursts collected at BATSE Catalog. The tests use the counts-in-cells method, an analysis of spherical harmonics, a test based on the two-point correlation function and a method based on multiscale methods. The tests suggest that the intermediate subclass of gamma-ray bursts are distributed anisotropically.
Many modern cosmological scenarios feature large volumes of spacetime in a de Sitter vacuum phase. Such models are said to be faced with a "Boltzmann Brain problem" - the overwhelming majority of observers with fixed local conditions are random fluctuations in the de Sitter vacuum, rather than arising via thermodynamically sensible evolution from a low-entropy past. We argue that this worry can be straightforwardly avoided in the Many-Worlds (Everett) approach to quantum mechanics, as long as the underlying Hilbert space is infinite-dimensional. In that case, de Sitter settles into a truly stationary quantum vacuum state. While there would be a nonzero probability for observing Boltzmann-Brain-like fluctuations in such a state, "observation" refers to a specific kind of dynamical process that does not occur in the vacuum (which is, after all, time-independent). Observers are necessarily out-of-equilibrium physical systems, which are absent in the vacuum. Hence, the fact that projection operators corresponding to states with observers in them do not annihilate the vacuum does not imply that such observers actually come into existence. The Boltzmann Brain problem is therefore much less generic than has been supposed.
In this paper we present a dynamical analysis for a coupled tachyonic dark energy with dark matter. The tachyonic field $\phi$ is considered in the presence of barothropic fluids (matter and radiation) and the autonomous system due to the evolution equations is studied. The three cosmological eras (radiation, matter and dark energy) are described through the critical points, for a generic potential $V(\phi)$.
$f(R)$ gravity is one of the simplest theories of modified gravity to explain the accelerated cosmic expansion. Although it is usually assumed that the quasi-Newtonian approach for cosmic perturbations is good enough to describe the evolution of large scale structure in $f(R)$ models, some studies have suggested that this method is not valid for all $f(R)$ models. Here, we show that in the matter-dominated era, the pressure and shear equations alone, which can be recast into four first-order equations to solve for cosmological perturbations exactly, are sufficient to solve for the Newtonian potential, $\Psi$, and the curvature potential, $\Phi$. Based on these two equations, we are able to clarify how the exact linear perturbations fit into different limits. We find that in the subhorizon limit, the so called quasi-static assumption plays no role in reducing the exact linear perturbations in any viable $f(R)$ gravity. Our findings also disagree with previous studies where we find little difference between our exact solutions and that of the quasi-Newtonian approach even up to $k=10 c^{-1} H_0$.
During the violent relaxation of a self-gravitating system a significant fraction of its mass may be ejected. If the time varying gravitational field also breaks spherical symmetry this mass can potentially carry angular momentum. Thus starting initial configurations with zero angular momentum can in principle lead to a bound virialized system with non-zero angular momentum. We explore here, using numerical simulations, how much angular momentum can be generated in a virialized structure in this way, starting from configurations of cold particles which are very close to spherically symmetric. For initial configurations in which spherical symmetry is broken only by the Poissonian fluctuations associated with the finite particle number $N$, with $N$ in range $10^3$ to $10^5$, we find that the relaxed structures have standard "spin" parameters $\lambda \sim 10^{-3}$, and decreasing slowly with $N$. For slightly ellipsoidal initial conditions, in which the finite-$N$ fluctuations break the residual reflection symmetries, we observe values $\lambda \sim 10^{-2}$ i.e. of the same order of magnitude as those reported for elliptical galaxies and dark matter halos in cosmological simulations. The net angular momentum vector is typically aligned close to normal to the major semi-axis of the triaxial relaxed structure, and also with that of the ejected mass. This simple mechanism may provide an alternative, or complement, to "tidal torque theory" for understanding the origin of angular momentum in astrophysical structures.
In a series of papers Kallosh, Linde, and collaborators have provided a unified description of single-field inflation with several types of potentials, ranging from power law to supergravity, in terms of just one parameter $\alpha$. These so-called $\alpha$-attractors predict a spectral index $n_{s}$ and a tensor-to-scalar ratio $r$, which are fully compatible with the latest Planck data. The only common feature of all $\alpha$-attractors is the analyticity of the scalar potential in the non-canonical Einstein frame. In this paper we explore the case of non-analytic potentials and we find that they lead to a class of attractors characterized by quasi-scale invariance in the Jordan frame. In the canonical Einstein frame they all converge to a model with a linear potential and a universal relation between $r$ and $n_{s}$ that can fit the observational data. We show that the breaking of exact, classical, scale invariance in the Jordan frame can be attributed to one-loop corrections, in line with previous results.
We present high-resolution observations, using both the European VLBI Network (EVN) at 1.7-GHz, and the Very Long Baseline Array (VLBA) at 5 and 8.4-GHz to image radio structures of 14 compact sources classified as broad absorption line (BAL) quasars based on the absorption index (AI). All source but one were resolved, with the majority showing core-jet morphology typical for radio-loud quasars. We discuss in details the most interesting cases. The high radio luminosities and small linear sizes of the observed objects indicate they are strong young AGNs. Nevertheless, the distribution of the radio-loudness parameter, log(Ri), of a larger sample of AI quasars shows that the objects observed by us constitute the most luminous, small subgroup of AI population. Additionally we report that for the radio-loudness parameter, the distribution of AI quasars and those selected by using the traditional balnicity index (BI), BI quasars differ significantly. Strong absorption is connected with the lower log(Ri), and thus probably with larger viewing angles. Since, the AI quasars have on average larger log(Ri), the orientation can cause that we see them less absorbed. However, we suggest that the orientation is not the only parameter that affects the detected absorption. The fact that the strong absorption is associated with the weak radio emission is equally important and worth exploring.
We present deep Hubble Space Telescope imaging at the locations of four, potentially hostless, long-faded Type Ia supernovae (SNe Ia) in low-redshift, rich galaxy clusters that were identified in the Multi-Epoch Nearby Cluster Survey. Assuming a steep faint-end slope for the galaxy cluster luminosity function ($\alpha_d=-1.5$), our data includes all but $\lesssim0.2\%$ percent of the stellar mass in cluster galaxies ($\lesssim0.005\%$ with $\alpha_d=-1.0$), a factor of 10 better than our ground-based imaging. Two of the four SNe Ia still have no possible host galaxy associated with them ($M_R>-9.2$), confirming that their progenitors belong to the intracluster stellar population. The third SNe Ia appears near a faint disk galaxy ($M_V=-12.2$) which has a relatively high probability of being a chance alignment. A faint, red, point source coincident with the fourth SN Ia's explosion position ($M_V=-8.4$) may be either a globular cluster (GC) or faint dwarf galaxy. We estimate the local surface densities of GCs and dwarfs to show that a GC is more likely, due to the proximity of an elliptical galaxy, but neither can be ruled out. This faint host implies that the SN Ia rate in dwarfs or GCs may be enhanced, but remains within previous observational constraints. We demonstrate that our results do not preclude the use of SNe Ia as bright tracers of intracluster light at higher redshifts, but that it will be necessary to first refine the constraints on their rate in dwarfs and GCs with deep imaging for a larger sample of low-redshift, apparently hostless SNe Ia.
Thusfar, there does not appear to be an agreed definition of homogeneous dark energy (DE). In this work, we argue that a correct definition of homogeneous DE is one whose density perturbation in comoving gauge vanishes. Using different DE models, we then investigate the consequence of this approach in the power spectrum -- with all the power spectra being normalized to match each other on small scales, at z = 0. We find that on super-Hubble scales, relativistic corrections in the observed galaxy power spectrum are able to distinguish a homogeneous DE from the concordance model and from a clustering DE, at low z and for high magnification bias. However, the matter power spectrum: is incapable of distinguishing a homogeneous DE from the concordance model (on all scales), at z = 0; but is able to differentiate it from a clustering DE, particularly at low z. Moreover, we found that relativistic effects become enhanced with decreasing magnification bias, and with increasing z.
Links to: arXiv, form interface, find, astro-ph, recent, 1505, contact, help (Access key information)
The next generation of large sky photometric surveys will finally be able to use arc statistics as a cosmological probe. Here we present the first of a series of papers on this topic. In particular, we study how arc counts are sensitive to the variation of two cosmological parameters: the (total) matter density parameter, $\Omega_m$, and the normalisation of the primordial power spectrum, expressed in terms of $\sigma_8$. Both these parameters influence the abundances of collapsed structures and their internal structure. We compute the expected number of gravitational arcs with various length-to-width ratios in mock light cones, by varying these cosmological parameters in the ranges $0.1\leq\Omega_m\leq0.5$ and $0.6\leq\sigma_8\leq 1$. We find that the arc counts dependence on $\Omega_m$ and $\sigma_8$ is similar, but not identical, to that of the halo counts. We investigate how the precision of the constraints on the cosmological parameters based on arc counts depends on the survey area. We find that the constraining power of arc statistics degrades critically only for surveys covering an area smaller than $10\%$ of the whole sky. Finally, we consider the case in which the search for arcs is done only in frames where galaxy clusters have been previously identified. Adopting the selection function for galaxy clusters expected to be detected from photometric data in future wide surveys, we find that less than $10\%$ of the arcs will be missed, with only a small degradation of the corresponding cosmological constraints.
We describe a novel and fully Bayesian approach to detect halos of different masses in cosmological observations and to quantify corresponding uncertainties. To demonstrate the capability of this approach we perform a Bayesian analysis of a realistic spectroscopic galaxy mock catalogue with the previously developed HADES inference algorithm. This procedure provides us with accurately inferred three-dimensional density fields and corresponding quantification of uncertainties inherent to any cosmological observation. Based upon these results we develop a novel Bayesian methodology to detect halos of different masses in cosmological observations subject to noise and systematic uncertainties. Specifically, we use a Bayesian chain rule to connect properties of halos found in simulations with actual observations. In an entirely Bayesian manner, this approach returns detection probabilities for halos above specified mass thresholds throughout the entire observed domain. We present maps of such detection probabilities and demonstrate the validity of this approach within our mock scenario. The proposed methodology can easily be extended to account for more complex scientific questions and is a promising novel tool to analyse the cosmic large-scale structure in observations.
We present and make publicly available the first data release (DR1) of the Keck Observatory Database of Ionized Absorption toward Quasars (KODIAQ) survey. The KODIAQ survey is aimed at studying galactic and circumgalactic gas in absorption at high-redshift, with a focus on highly-ionized gas traced by OVI, using the HIRES spectrograph on the Keck-I telescope. KODIAQ DR1 consists of a fully-reduced sample of 170 quasars at 0.29 < z_em < 5.29 observed with HIRES at high resolution (36,000 <= R <= 103,000) between 2004 and 2012. DR1 contains 247 spectra available in continuum normalized form, representing a sum total exposure time of ~1.6 megaseconds. These co-added spectra arise from a total of 567 individual exposures of quasars taken from the Keck Observatory Archive (KOA) in raw form and uniformly processed using a HIRES data reduction package made available through the XIDL distribution. DR1 is publicly available to the community, housed as a higher level science product at the KOA. We will provide future data releases that make further QSOs, including those with pre-2004 observations taken with the previous-generation HIRES detectors.
We propose a new and robust method to test the consistency of the cosmic evolution given by a cosmological model. It is realized by comparing the combined quantity r_d^CMB/D_V^SN, which is derived from the comoving sound horizon r_d from cosmic microwave background (CMB) measurements and the effective distance D_V derived from low-redshift Type-Ia supernovae (SNe Ia) data, with direct and independent r_d/D_V obtained by baryon acoustic oscillation (BAO) measurements at median redshifts. We apply this test method for the LCDM and wCDM models, and investigate the consistency of the derived value of r_d/D_V from Planck 2015 and the SN Ia data sets of Union2.1 and JLA (z<1.5), and the r_d/D_V directly given by BAO data from six-degree-field galaxy survey (6dFGS), Sloan Digital Sky Survey Data Release 7 Main Galaxy Survey (SDSS-DR7 MGS), DR11 of SDSS-III, WiggleZ and Ly-alpha forecast surveys from Baryon Oscillation Spectroscopic Data (BOSS) DR-11 over 0.1<z<2.36. We find that the non-flat LCDM model is well consistent with the BAO measurements within 1-sigma CL across the whole redshift range (0<z<2.4), while the flat wCDM model is in tension with BAO data at 1.5 sigma CL at z=2.34 and z=2.36. Future surveys will further tight up the constraints significantly, and provide stronger test on the consistency.
We show that emerging tension between the direct astronomical measurements at low redshifts and cosmological parameters deduced from the Planck measurements of the CMB anisotropies can be alleviated if the dark matter consists of two fractions, stable part being dominant and a smaller unstable fraction constituting of about 5 - 10 per cent at the recombination epoch. The latter decays after cosmological recombination but earlier than the present epoch.
Using N-body simulations, we measure the power spectrum of the effective dark matter density field, which is defined through the modified Poisson equation in $f(R)$ cosmologies. We find that when compared to the conventional dark matter power spectrum, the effective power spectrum deviates more significantly from the $\Lambda$CDM model. For models with $f_{R0}=-10^{-4}$, the deviation can exceed 150\% while the deviation of the conventional matter power spectrum is less than 50\%. Even for models with $f_{R0}=-10^{-6}$, for which the conventional matter power spectrum is very close to the $\Lambda$CDM prediction, the effective power spectrum shows sizeable deviations. Our results indicate that traditional analyses based on the dark matter density field may seriously underestimate the impact of $f(R)$ gravity on galaxy clustering. We therefore suggest the use of the effective density field in such studies.
We study the cross-correlation of distribution of galaxies, the Sunyaev-Zel'dovich (SZ) and X-ray power spectra of galaxies from current and upcoming surveys and show these to be excellent probes of the nature, i.e. extent, evolution and energetics, of the circumgalactic medium (CGM). For a flat pressure profile, the SZ cross power spectrum shows oscillations at $l$-values corresponding to the length scales smaller than $\sim \frac{2}{3}$ times the virial radius of the galaxy. These oscillations are sensitive to the steepness of the pressure profile of the CGM and vanish for a sufficiently steep profile. Similar oscillations are also present in the X-ray cross power spectrum which is, however, more sensitive to the density profile. We forecast the detectability of the cross-correlated galaxy distribution, SZ and X-ray signals by combining SPT-DES and eROSITA-DES/eROSITA-LSST surveys, respectively. We find that, for the SPT-DES survey, the signal-to-noise ratio (SNR) peaks at high mass and redshift with SNR $\sim 9$ around $M_h\sim 10^{13} h^{-1} M_{\odot}$ and $z\sim 1.5\hbox{--} 2$ for flat density and temperature profiles. The SNR peaks at $\sim 6 (12 )$ for the eROSITA-DES (eROSITA-LSST) surveys. We also perform a Fisher matrix analysis and find that the gas fraction in the CGM can be constrained to a precision of $\sim 34\% (23 \%)$ by the SPT-DES and $\sim 23\% (14 \%)$ by the eROSITA-DES surveys in the presence (absence) of an unknown redshift evolution of the gas fraction. Finally, we demonstrate that the cross-correlated SZ-galaxy and X-ray-galaxy power spectrum can be used as powerful probes of the CGM energetics and potentially discriminate between different feedback models recently proposed in the literature; for example, one can distinguish a `no AGN feedback' scenario from a CGM energized by `fixed-velocity hot winds' at greater than $3\sigma$.
We use three domain wall simulations from the radiation era to the late time dark energy domination era based on the PRS algorithm to calculate the energy-momentum tensor components of domain wall networks in an expanding universe. Unequal time correlators in the radiation, matter and cosmological constant epochs are calculated using the scaling regime of each of the simulations. The CMB power spectrum of a network of domain walls is determined. The first ever quantitative constraint for the domain wall surface tension is obtained using a Markov chain Monte Carlo method; an energy scale of domain walls of 0.93 MeV, which is close but below the Zel'dovich bound, is determined.
We estimate the degree of circular polarization for the gravitational waves generated during the electroweak and QCD phase transitions from the kinetic and magnetic helicity generated by bubble collisions during those cosmological phase transitions.
We obtain a constraint on the parameters of a generalized cubic Galileon gravity model exhibiting the Vainshtein mechanism by using multi-wavelength observations of the Coma Cluster. The generalized cubic Galileon model is characterized by three parameters of the turning scale associated with the Vainshtein mechanism, and the amplitude of modifying a gravitational potential and a lensing potential. X-ray and Sunyaev-Zel'dovich (SZ) observations of the intra-cluster medium are sensitive to the gravitational potential, while the weak-lensing (WL) measurement is specified by the lensing potential. A joint fit of a complementary multi-wavelength dataset of X-ray, SZ and WL measurements enables us to simultaneously constrain these three parameters of the generalized cubic Galileon model for the first time.
We revisit the kink-like parametrization of the deceleration parameter ($q(z)$) \cite{ishida08}, which considers a transition, at redshift $z_t$, from cosmic deceleration to acceleration. In this parametrization the initial ($z \gg z_t$) value of the q-parameter is $q_i$, its final ($z=-1$) value is $q_f$ and the duration of the transition is parametrized by $\tau$. We obtain constraints on the free parameters of the model using recent data from type Ia supernovae (SN Ia), baryon acoustic oscillations (BAO), cosmic microwave background (CMB) and the Hubble parameter (H(z)). The use of H(z) data introduces an explicit dependence of the combined likelihood on the present value of the Hubble parameter ($H_0$), allowing us to explore the influence of different priors when marginalizing over this parameter. We also study the importance of the CMB information in the results by considering data from WMAP7, WMAP9 (Wilkinson Microwave Anisotropy Probe - 7 and 9 years) and the Planck satellite. Assuming a flat space geometry, $q_i=1/2$ and expressing the present value of the deceleration parameter ($q_0$) as a function of the others three free parameters, we obtain $z_t=0.68^{+0.10}_{-0.09}$, $\tau=0.24^{+0.14}_{-0.10}$ and $q_0=-0.48^{+0.11}_{-0.14}$, at 68\% of confidence level, with flat prior on $H_0$. If in addition we fix $q_f=-1$, we get $z_t=0.67^{+0.03}_{-0.04}$, $\tau=0.32^{+0.04}_{-0.03}$ and $q_0=-0.55^{+0.05}_{-0.06}$.
The scale of Baryon Acoustic Oscillations (BAO) imprinted in the matter power spectrum provides an almost-perfect standard ruler: it only suffers sub-percent deviations from fixed comoving length due to non-linear effects. We study the BAO shift in the large Horndeski class of gravitational theories and compute its magnitude in momentum space using second-order perturbation theory and a peak-background split. The standard prediction is affected by the modified linear growth, as well as by non-linear gravitational effects that alter the mode-coupling kernel. For covariant Galileon models, we find a $14-45\%$ enhancement of the BAO shift with respect to standard gravity and a distinct time evolution depending on the parameters. Despite the larger values, the shift remains well below the forecasted precision of next-generation galaxy surveys. Models that produce significant BAO shift would cause large redshift-space distortions or affect the bispectrum considerably. Our computation therefore validates the use of the BAO scale as a comoving standard ruler for tests of general dark energy models.
The data analysis of the gravitational wave signals emitted by coalescing neutron star binaries requires the availability of an accurate analytical representation of the dynamics and waveforms of these systems. We propose an effective-one-body (EOB) model that describes the general relativistic dynamics of neutron star binaries from the early inspiral up to merger. Our EOB model incorporates an enhanced attractive tidal potential motivated by recent analytical advances in the post-Newtonian and gravitational self-force description of relativistic tidal interactions. No fitting parameters are introduced for the description of tidal interaction in the late, strong-field dynamics. We compare the model energetics and the gravitational wave phasing with new high-resolution multi-orbit numerical relativity simulations of equal-mass configurations with different equations of state. We find agreement within the uncertainty of the numerical data for all configurations. Our model is the first semi-analytical model which captures the tidal amplification effects close to merger. It thereby provides the most accurate analytical representation of binary neutron star dynamics and waveforms currently available.
Self-accelerating solutions in massive gravity provide explicit, calculable examples that exhibit the general interplay between superluminality, the well-posedness of the Cauchy problem, and strong coupling. For three particular classes of vacuum solutions, one of which is new to this work, we construct the conformal diagram for the characteristic surfaces on which isotropic stress-energy perturbations propagate. With one exception, all solutions necessarily possess spacelike characteristics, indicating perturbative superluminality. Foliating the spacetime with these surfaces gives a pathological frame where kinetic terms of the perturbations vanish, confusing the Hamiltonian counting of degrees of freedom. This frame dependence distinguishes the vanishing of kinetic terms from strong coupling of perturbations or an ill-posed Cauchy problem. We give examples where spacelike characteristics do and do not originate from a point where perturbation theory breaks down and where spacelike surfaces do or do not intersect all characteristics in the past light cone of a given observer. The global structure of spacetime also reveals issues that are unique to theories with two metrics: in all three classes of solutions, the Minkowski fiducial space fails to cover the entire de Sitter spacetime allowing worldlines of observers to end in finite proper time at determinant singularities. Characteristics run tangent to these surfaces requiring {\it ad hoc} rules to establish continuity across singularities.
We present statistics of 89 faint 1.2-mm continuum sources with a flux density of ~0.01-1 mJy detected by about 100 deep ALMA pointing data that include the complete deep datasets archived by 2015 March. These faint sources are identified in 50 blank fields and behind one cluster, Abell 1689, that magnifies the background sources by gravitational lensing. Evaluating various important effects including the false detection, detection completeness, and flux boosting as well as the lensing magnification by modeling and simulations, we derive number counts of 1.2 mm continuum sources. We find that the number counts are well represented by the Schechter function down to ~0.01 mJy, and that the total integrated 1.2 mm flux of the securely identified sources is 22.8^(+6.1)_(-6.4) Jy deg^(-2) that corresponds to 104^(+27)_(-30)% of the extragalactic background light (EBL) measured by COBE observations. These results suggest that the major 1.2 mm EBL contributors are sources with >~0.01 mJy, and that very faint 1.2 mm sources with <~ 0.01 mJy contribute negligibly to the EBL with the possible flattening and/or truncation of number counts in this very faint flux regime. To understand the physical origin of our faint ALMA sources, we measure the galaxy bias bg by the counts-in-cells technique under the assumption that the sources reside at z~2.5, and place a stringent upper limit of bg<4.1 that is not similar to bg values of massive DRGs and SMGs but comparable to those of UV-bright sBzKs and LBGs. Moreover, in optical and near-infrared (NIR) deep fields, we identify optical-NIR counterparts for 54% of our faint ALMA sources, majority of which have luminosities and colors same as sBzKs and LBGs. We thus conclude that about a half of our faint ALMA sources are dust-poor high-z galaxies as known as sBzKs and LBGs in optical studies.
B and V time-series photometry of the M31 dwarf spheroidal satellite Andromeda XXI (And XXI) was obtained with the Large Binocular Cameras at the Large Binocular Telescope. We have identified 50 variables in And XXI, of which 41 are RR Lyrae stars (37 fundamental-mode RRab, and 4 first-overtone RRc, pulsators) and 9 are Anomalous Cepheids (ACs). The average period of the RRab stars (<Pab> = 0.64 days) and the period-amplitude diagram place And~XXI in the class of Oosterhoff II - Oosterhoff-Intermediate objects. From the average luminosity of the RR Lyrae stars we derived the galaxy distance modulus of (m-M)$_0$=$24.40\pm0.17$ mag, which is smaller than previous literature estimates, although still consistent with them within 1 $\sigma$. The galaxy color-magnitude diagram shows evidence for the presence of three different stellar generations in And~XXI: 1) an old ($\sim$ 12 Gyr) and metal poor ([Fe/H]=$-$1.7 dex) component traced by the RR Lyrae stars; 2) a slightly younger (10-6 Gyr) and more metal rich ([Fe/H]=$-$1.5 dex) component populating the red horizontal branch, and 3) a young age ($\sim$ 1 Gyr) component with same metallicity, that produced the ACs. Finally, we provide hints that And~XXI could be the result of a minor merging event between two dwarf galaxies.
Theories of modified gravity, in both the linear and fully non-linear regime, are often studied under the assumption that the evolution of the new (often scalar) degree of freedom present in the theory is quasi-static. This approximation significantly simplifies the study of the theory, and one often has good reason to believe that it should hold. Nevertheless it is a crucial assumption that should be explicitly checked whenever possible. In this paper we do so for the Vainshtein mechanism. By solving for the full spatial and time evolution of the Dvali-Gabadadze-Porrati and the Cubic Galileon model, in a spherical symmetric spacetime, we are able to demonstrate that the Vainshtein solution is a stable attractor and forms no matter what initial conditions we take for the scalar field. Furthermore,the quasi-static approximation is also found to be a very good approximation whenever it exists. For the best-fit Cubic Galileon model, however, we find that for deep voids at late times, the numerical solution blows up at the same time as the quasi-static solution ceases to exist. We argue that this phenomenon is a true instability of the model.
In this work we consider an inflationary potential in fractional form coupled non-minimally to gravity. This potential has a dominant constant energy density at early times which can realize successful inflation. We show that this potential predicts small tensor-to-scalar ratio of the order of $r\approx 0.01$ which is fully consistent with Planck constraints. Using the lower and upper bounds on reheating temperature, we provide additional constraints on the non-minimally coupling parameter $\xi$ of the model. We also study the preheating stage predicted by this kind of potentials using numerical calculations.
We compared the total mass density profiles of three different types of galaxies using weak gravitational lensing: (i) 29 galaxies that host quasars at z~0.32 that are in a post-starburst (PSQ) phase with high star formation indicating recent merger activity, (ii) 22 large elliptical galaxies from the SLACS sample that do not host a quasar at z~0.23, and (iii) 17 galaxies that host moderately luminous quasars at z~0.36 powered by disk instabilities, but with no intense star formation. On an initial test we found no evidence for a connection between the merger state of a galaxy and the profile of the halo, with the PSQ profile comparable to that of the other two samples and consistent with the Leauthaud et al. (2014) study of moderately luminous quasars in COSMOS. Given the compatibility of the two quasar samples, we combined these and found no evidence for any connection between black hole activity and the dark matter halo. All three mass profiles remained compatible with isothermality given the present data.
Dark energy is investigated from the perspective of quantum cosmology. By treating the existence of a classical universe as a constraint, it is found that the normal ordering ambiguity factor q in Wheeler-DeWitt equation tends to take its value on domain (-1, 3). Furthermore, to ensure the existence of a classical universe, there must be dark energy in the universe. It is in this sense we propose that dark energy is the reason for the existence of a classical universe.
Links to: arXiv, form interface, find, astro-ph, recent, 1505, contact, help (Access key information)