Links to: arXiv, form interface, find, astro-ph, recent, 1710, contact, help (Access key information)
With the advent of wide-field cosmological surveys, we are approaching samples of hundreds of thousands of galaxy clusters. While such large numbers will help reduce statistical uncertainties, the control of systematics in cluster masses becomes ever more crucial. Here we examine the effects of an important source of systematic uncertainty in galaxy-based cluster mass estimation techniques: the presence of significant dynamical substructure. Dynamical substructure manifests as dynamically distinct subgroups in phase-space, indicating an 'unrelaxed' state. This issue affects around a quarter of clusters in a generally selected sample. We employ a set of mock clusters whose masses have been measured homogeneously with commonly-used galaxy-based mass estimation techniques (kinematic, richness, caustic, radial methods). We use these to study how the relation between observationally estimated and true cluster mass depends on the presence of substructure, as identified by various popular diagnostics. We find that the scatter for an ensemble of clusters does not increase dramatically for clusters with dynamical substructure. However, we find a systematic bias for all methods, such that clusters with significant substructure have higher measured masses than their relaxed counterparts. This bias depends on cluster mass: the most massive clusters are largely unaffected by the presence of significant substructure, but masses are significantly overestimated for lower mass clusters, by $\sim10\%$ at $10^{14}$ and $\geq20\%$ for $\leq10^{13.5}$. The use of cluster samples with different levels of substructure can, therefore, bias certain cosmological parameters up to a level comparable to the typical uncertainties in current cosmological studies.
Measurements of the cosmic microwave background (CMB) temperature anisotropies have revealed a dipolar asymmetry in power at the largest scales, in apparent contradiction with the statistical isotropy of standard cosmological models. The significance of the effect is not very high, and is dependent on a posteriori choices. Nevertheless, a number of models have been proposed that produce a scale-dependent asymmetry. We confront several such models for a physical, position-space modulation with CMB temperature observations. We find that, while some models that maintain the standard isotropic power spectrum are allowed, others, such as those with modulated tensor or uncorrelated isocurvature modes, can be ruled out on the basis of the overproduction of isotropic power. This remains the case even when an extra isocurvature mode fully anti-correlated with the adiabatic perturbations is added to suppress power on large scales.
We present an algorithm enabling computation of the anisotropic redshift-space galaxy 3-point correlation function (3PCF) scaling as $N^2$, with $N$ the number of galaxies. Our previous work showed how to compute the isotropic 3PCF with this scaling by expanding the radially-binned density field around each galaxy in the survey into spherical harmonics and combining these coefficients to form multipole moments. The $N^2$ scaling occurred because this approach never explicitly required the relative angle between a galaxy pair about the primary galaxy. Here we generalize this work, demonstrating that in the presence of azimuthally-symmetric anisotropy produced by redshift-space distortions (RSD) the 3PCF can be described by two triangle side lengths, two independent total angular momenta, and a spin. This basis for the anisotropic 3PCF allows its computation with negligible additional work over the isotropic 3PCF. We also present the covariance matrix of the anisotropic 3PCF measured in this basis. Our algorithm tracks the full 5-D redshift-space 3PCF, uses an accurate line of sight to each triplet, is exact in angle, and easily handles edge correction. It will enable use of the anisotropic large-scale 3PCF as a probe of RSD in current and upcoming large-scale redshift surveys.
Due to its accuracy and generality, Monte Carlo radiative transfer (MCRT) has emerged as the prevalent method for Ly$\alpha$ radiative transfer in arbitrary geometries. The standard MCRT encounters a significant efficiency barrier in the high optical depth, diffusion regime. Multiple acceleration schemes have been developed to improve the efficiency of MCRT but the noise from photon packet discretization remains a challenge. The discrete diffusion Monte Carlo (DDMC) scheme has been successfully applied in state-of-the-art radiation hydrodynamics (RHD) simulations. Still, the established framework is not optimal for resonant line transfer. Inspired by the DDMC paradigm, we present a novel extension to resonant DDMC in which diffusion in space and frequency are treated on equal footing. We explore the robustness of our new method and demonstrate a level of performance that justifies incorporating the method into existing Ly$\alpha$ codes. We present computational speedups of $\sim 10^2$-$10^6$ relative to contemporary MCRT implementations with aggressive core-skipping. This is because the resonant DDMC runtime scales with the spatial and frequency resolution rather than the number of scatterings - the latter is typically $\propto \tau_0$ for static media, or $\propto (a \tau_0)^{2/3}$ with core-skipping. We anticipate new frontiers in which on-the-fly Ly$\alpha$ radiative transfer calculations are feasible in 3D RHD. More generally, resonant DDMC is transferable to any computationally demanding problem amenable to a Fokker-Planck approximation of frequency redistribution.
Giga-parsec scale alignments of the quasar optical polarization vectors have been proven to be robust against a scenario of contamination by the Galactic interstellar medium (ISM). This claim has been established by means of optical polarization measurements of the star light surrounding the lines of sight of the 355 quasars for which reliable optical polarization measurements are available. Taking advantage of the full-sky and high quality polarization data released by the $Planck$ satellite, we provide an independent, complementary and up-to-date estimation of the contamination level of the quasar optical polarization data by the diffuse Galactic thermal dust emission. Our analysis reveals contamination signatures of the quasar optical polarization data by Galactic dust for about twenty per cent of the sample. We provide new arguments showing that Galactic thermal dust cannot account for the reported quasar optical polarization alignments. We provide new and independent quality criteria to apply to the polarization sample for future inquiries in the framework of the cosmological-scale correlations of quasar polarization vector orientations that could compete with the isotropic principle of the cosmological paradigm.
Observations of galaxies and galaxy clusters in the local universe can account for only $10\%$ of the baryon content inferred from measurements of the cosmic microwave background and from nuclear reactions in the early Universe. Locating the remaining $90\%$ of baryons has been one of the major challenges in modern cosmology. Cosmological simulations predict that the 'missing baryons' are spread throughout filamentary structures in the cosmic web, forming a low density gas with temperatures of $10^5-10^7$ K. Previous attempts to observe this warm-hot filamentary gas via X-ray emission or absorption in quasar spectra have proven difficult due to its diffuse and low-temperature nature. Here we report a $5.1 \sigma$ detection of warm-hot baryons in stacked filaments through the thermal Sunyaev-Zel'dovich (SZ) effect, which arises from the distortion in the cosmic microwave background spectrum due to ionised gas. The estimated gas density in these 15 Megaparsec-long filaments is approximately 6 times the mean universal baryon density, and overall this can account for $\sim 30\%$ of the total baryon content of the Universe. This result establishes the presence of ionised gas in large-scale filaments, and suggests that the missing baryons problem may be resolved via observations of the cosmic web.
The efficient simulation of isotropic Gaussian random fields on the unit sphere is a task encountered frequently in numerical applications. A fast algorithm based on Markov properties and Fast Fourier Transforms in 1d is presented that generates samples on an n x n grid in O(n^2 log n). Furthermore, an efficient method to set up the necessary conditional covariance matrices is derived and simulations demonstrate the performance of the algorithm.
We discuss a cosmological phase transition within the Standard Model which incorporates spontaneously broken scale invariance as a low-energy theory. In addition to the Standard Model fields, the minimal model involves a light dilaton, which acquires a large vacuum expectation value (VEV) through the mechanism of dimensional transmutation. Under the assumption of the cancellation of the vacuum energy, the dilaton develops a very small mass at 2-loop order. As a result, a flat direction is present in the classical dilaton-Higgs potential at zero temperature while the quantum potential admits two (almost) degenerate local minima with unbroken and broken eletroweak symmetry. We found that the cosmological electroweak phase transition in this model can only be triggered by a QCD chiral symmetry breaking phase transition at low temperatures, $T\lesssim 132$ MeV. Furthermore, unlike the standard case, the universe settles into the chiral symmetry breaking vacuum via a first-order phase transition which gives rise to a stochastic gravitational background with a peak frequency $\sim 10^{-8}$ Hz as well as triggers the production of approximately solar mass primordial black holes. The observation of these signatures of cosmological phase transitions together with the detection of a light dilaton would provide a strong hint of the fundamental role of scale invariance in particle physics.
The Fermi Collaboration has recently updated their analysis of gamma rays from the center of the Galaxy. They reconfirm the presence of an unexplained emission feature which is most prominent in the region of $1-10$ GeV, known as the Galactic Center GeV excess (GCE). Although the GCE is now firmly detected, an interpretation of this emission as a signal of self-annihilating dark matter (DM) particles is not unambiguously possible due to systematic effects in the gamma-ray modeling estimated in the Galactic Plane. In this paper we build a covariance matrix, collecting different systematic uncertainties investigated in the Fermi Collaboration's paper that affect the GCE spectrum. We show that models where part of the GCE is due to annihilating DM can still be consistent with the new data. We also re-evaluate the parameter space regions of the minimal supersymmetric Standard Model (MSSM) that can contribute dominantly to the GCE via neutralino DM annihilation. All recent constraints from DM direct detection experiments such as PICO, LUX, PandaX and Xenon1T, limits on the annihilation cross section from dwarf spheroidal galaxies and Large Hadron Collider limits are considered in this analysis. Due to a slight shift in the energy spectrum of the GC excess with respect to the previous Fermi analysis, and the recent limits from direct detection experiments, we find a slightly shifted parameter region of the MSSM compared to our previous analysis that is consistent with the GCE. Neutralinos with a mass between $85-220$ GeV can describe the excess via annihilation into a pair of $W$-bosons or top quarks. Remarkably, there are low fine-tuning models among the regions that we have found. The complete set of solutions will be probed by upcoming direct detection experiments and with dedicated searches in the upcoming data of the Large Hadron Collider.
In this work we examine what are the cosmological implications of allowing the geometrical curvature density to behave independently from the energy density contents. Using the full data extracted by Planck mission from CMB, combined with BAO and SNIa measurements, we derive, in the light of this approach, new constraints on the cosmological parameters. In particular we determine the behavior of the curvature dark energy degeneracy when allowing a varying equation of state for the latter. We also examine whether this approach could bridge the gap recently found between the Hubble parameter value determined from CMB and that from the local universe measurements
Links to: arXiv, form interface, find, astro-ph, recent, 1710, contact, help (Access key information)
We present a Bayesian hierarchical inference formalism to study the relation between the properties of dark matter halos and those of their central galaxies using weak gravitational lensing. Unlike traditional methods, this technique does not resort to stacking the weak lensing signal in bins, and thus allows for a more efficient use of the information content in the data. Our method is particularly useful for constraining scaling relations between two or more galaxy properties and dark matter halo mass, and can also be used to constrain the intrinsic scatter in these scaling relations. We show that, if observational scatter is not properly accounted for, the traditional stacking method can produce biased results when exploring correlations between multiple galaxy properties and halo mass. For example, this bias can affect studies of the joint correlation between galaxy mass, halo mass, and galaxy size, or galaxy color. In contrast, our method easily and efficiently handles the intrinsic and observational scatter in multiple galaxy properties and halo mass. We test our method on mocks with varying degrees of complexity. We find that we can recover the mean halo mass and concentration, each with a $0.1$ dex accuracy, and the intrinsic scatter in halo mass with a $0.05$ dex accuracy. In its current version, our method will be most useful for studying the weak lensing signal around central galaxies in groups and clusters, as well as massive galaxies samples with $\log{M_*} > 11$, which have low satellite fractions.
The recent discovery of fast transient events near critical curves of massive galaxy clusters, which are interpreted as highly magnified individual stars in giant arcs due to caustic crossing, opens up the possibility of using such microlensing events to constrain a range of dark matter models such as primordial black holes and scalar field dark matter. Based on a simple analytic model, we study lensing properties of a point mass lens embedded in a high magnification region, and derive the dependence of the peak brightness, microlensing time scales, and event rates on the mass of the point mass lens as well as the radius of a source star that is magnified, assuming the transverse velocity of $500~{\rm km\,s^{-1}}$. We find that the lens mass and source radius of the first event MACS J1149 Lensed Star 1 (LS1) are constrained, with the lens mass range of $0.1~M_\odot \lesssim M \lesssim 4\times 10^3M_\odot$ and the source radius range of $40~R_\odot \lesssim R \lesssim 260~R_\odot$. In the most plausible case with $M\approx 0.3~M_\odot$ and $R\approx 180~R_\odot$, the source star should have been magnified by a factor of $\approx 4300$ at the peak. The derived lens properties are fully consistent with the interpretation that MACS J1149 LS1 is a microlensing event produced by a star that contributes to the intra-cluster light. We argue that compact dark matter models with high fractional mass densities for the mass range $10^{-5}M_\odot \lesssim M\lesssim 10^2M_\odot$ are inconsistent with the observation of MACS J1149 LS1 because such models predict too low magnifications. Our work demonstrates a potential use of caustic crossing events in giant arcs to constrain compact dark matter.
We propose a multi-scale edge-detection algorithm to search for the Gott-Kaiser-Stebbins imprints of a cosmic string (CS) network on the Cosmic Microwave Background (CMB) anisotropies. Curvelet decomposition and extended Canny algorithm are used to enhance the string detectability. Various statistical tools are then applied to quantify the deviation of CMB maps having a cosmic string contribution with respect to pure Gaussian anisotropies of inflationary origin. These statistical measures include the one-point probability density function, the weighted two-point correlation function (TPCF) of the anisotropies, the unweighted TPCF of the peaks and of the up-crossing map, as well as their cross-correlation. We use this algorithm on a hundred of simulated Nambu-Goto CMB flat sky maps, covering approximately $10\%$ of the sky, and for different string tensions $G\mu$. On noiseless sky maps with an angular resolution of $0.9'$, we show that our pipeline detects CSs with $G\mu$ as low as $G\mu\gtrsim 4.3\times 10^{-10}$. At the same resolution, but with a noise level typical to a CMB-S4 phase II experiment, the detection threshold would be to $G\mu\gtrsim 1.2 \times 10^{-7}$.
We initiate a study on various cosmological imprints of string axions whose scalar potentials have plateau regions. In such cases, we show that a delayed onset of oscillation generically leads to a parametric resonance instability. In particular, for ultralight axions, the parametric resonance can enhance the power spectrum slightly below the Jeans scale, alleviating the tension with the Lyman $\alpha$ forest observations. We also argue that a long-lasting resonance can lead to an emission of gravitational waves at the frequency bands which are detectable by gravitational wave interferometers and pulsar timing arrays.
We report results from 21-cm intensity maps acquired from the Parkes radio telescope and cross-correlated with galaxy maps from the 2dF galaxy survey. The data span the redshift range $0.057<z<0.098$ and cover approximately 1,300 square degrees over two long fields. Cross correlation is detected at a significance of $5.18\sigma$. The amplitude of the cross-power spectrum is low relative to the expected dark matter power spectrum, assuming a neutral hydrogen (HI) bias and mass density equal to measurements from the ALFALFA survey. The decrement is pronounced and statistically significant at small scales. At $k\sim1.5$ $ h \mathrm{Mpc^{-1}}$, the cross power spectrum is more than a factor of 6 lower than expected, with a significance of $14.8\,\sigma$. This decrement indicates either a lack of clustering of neutral hydrogen (HI), a small correlation coefficient between optical galaxies and HI, or some combination of the two. Separating 2dF into red and blue galaxies, we find that red galaxies are much more weakly correlated with HI on $k\sim1.5$ $h \mathrm{Mpc^{-1}}$ scales, suggesting that HI is more associated with blue star-forming galaxies and tends to avoid red galaxies.
A persistent signal of power asymmetry on opposite hemispheres of CMB sky was seen in full-sky temperature measurements made so far. This asymmetry was seen in microwave sky from WMAP as well as PLANCK satellites, and calls for attention the larger question of \emph{statistical isotropy}, one of the foundational principles of modern cosmology. In this work we present an analysis of polarized CMB maps from PLANCK 2015 full mission data. We apply the local variance estimator on low resolution $E-$mode maps from PLANCK 2015 polarization \texttt{Commander} solution. We find a significant hemispherical power asymmetry in polarization data on large angular scales, at the level of $\sim 2.6-3.9\%$ depending on the galactic mask, and the circular disc radius used for computing local variance maps. However the direction is found to be pointing broadly towards CMB kinetic dipole direction. Precise measurements of CMB polarization in future will shed light on this apparent discrepancy in the anisotropy axis seen in temperature and polarized CMB sky, and likely influence of systematics on our findings.
We decompose the Lyman-{\alpha} (Ly{\alpha}) forest of an extensive sample of 74 high signal-to-noise ratio and high-resolution quasar spectra into a collection of Voigt profiles. Absorbers located near caustics in the peculiar velocity field have the smallest Doppler parameters, resulting in a low-$b$ cutoff in the $b$-$N_{\text{HI}}$ set by the thermal state of intergalactic medium (IGM). We fit this cutoff as a function of redshift over the range $2.0\leq z \leq 3.4$, which allows us to measure the evolution of the IGM temperature-density ($T= T_0 (\rho/ \rho_0)^{\gamma-1}$) relation parameters $T_0$ and $\gamma$. We calibrate our measurements against Ly$\alpha$ forest simulations, using 21 different thermal models of the IGM at each redshift, also allowing for different values of the IGM pressure smoothing scale. We adopt a forward-modeling approach and self-consistently apply the same algorithms to both data and simulations, propagating both statistical and modeling uncertainties via Monte Carlo. The redshift evolution of $T_0$ shows a suggestive peak at $z=2.8$, while our evolution of $\gamma$ is consistent with $\gamma\simeq 1.4$ and disfavors inverted temperature-density relations. Our measured evolution of $T_0$ and $\gamma$ are generally in good agreement with previous determinations in the literature. Both the peak in the evolution of $T_0$ at $z = 2.8$, as well as the high temperatures $T_0\simeq 15000-20000\,$K that we observe at $2.4 < z < 3.4$, strongly suggest that a significant episode of heating occurred after the end of HI reionization, which was most likely the cosmic reionization of HeII.
We present proof of principle for a two way interplay between physics at very early Universe and late time observations. We find a relation between primordial mechanisms responsible for large scale power suppression in the primordial power spectrum and the value of reionization optical depth $\tau$. Such mechanisms can affect the estimation of $\tau$. We show that using future measurements of $\tau$, one can obtain constraints on the pre-inflationary dynamics, providing a new window on the physics of the very early Universe. Furthermore, the new, re-estimated $\tau$ can potentially resolve moderate discrepancy between high and low-$\ell$ measurements of $\tau$, hence providing empirical support for the power suppression anomaly and its primordial origin.
We study the morphologies and sizes of galaxies at z>5 using high-resolution cosmological zoom-in simulations from the Feedback In Realistic Environments project. The galaxies show a variety of morphologies, from compact to clumpy to irregular. The simulated galaxies have more extended morphologies and larger sizes when measured using rest-frame optical B-band light than rest-frame UV light; sizes measured from stellar mass surface density are even larger. The UV morphologies are usually dominated by several small, bright young stellar clumps that are not always associated with significant stellar mass. The B-band light traces stellar mass better than the UV, but it can also be biased by the bright clumps. At all redshifts, galaxy size correlates with stellar mass/luminosity with large scatter. The half-light radii range from 0.01 to 0.2 arcsec (0.05-1 kpc physical) at fixed magnitude. At z>5, the size of galaxies at fixed stellar mass/luminosity evolves as (1+z)^{-m}, with m~1-2. For galaxies less massive than M_star~10^8 M_sun, the ratio of the half-mass radius to the halo virial radius is ~10% and does not evolve significantly at z=5-10; this ratio is typically 1-5% for more massive galaxies. A galaxy's "observed" size decreases dramatically at shallower surface brightness limits. This effect may account for the extremely small sizes of z>5 galaxies measured in the Hubble Frontier Fields. We provide predictions for the cumulative light distribution as a function of surface brightness for typical galaxies at z=6.
We study a vector dark matter (VDM) model in which the dark sector couples to the Standard Model sector via a Higgs portal. If the portal coupling is small enough the VDM can be produced via the freeze-in mechanism. It turns out that the electroweak phase transition have a substantial impact on the prediction of the VDM relic density. We further assume that the dark Higgs boson which gives the VDM mass is so light that it can induce strong VDM self-interactions and solve the small-scale structure problems of the Universe. As illustrated by the latest LUX data, the extreme smallness of the Higgs portal coupling required by the freeze-in mechanism implies that the dark matter direct detection bounds are easily satisfied. However, the model is well constrained by the indirect detections of VDM from BBN, CMB, AMS-02, and diffuse $\gamma$/X-rays. Consequently, only when the dark Higgs boson mass is at most of ${\cal O}({\rm keV})$ does there exist a parameter region which leads to a right amount of VDM relic abundance and an appropriate VDM self-scattering while satisfying all other constraints simultaneously.
We present a Bayesian model for multi-resolution CMB component separation based on Wiener filtering and/or computation of constrained realizations, extending a previously developed framework. We also develop an efficient solver for the corresponding linear system for the associated signal amplitudes. The core of this new solver is an efficient preconditioner based on the pseudo-inverse of the coefficient matrix of the linear system. In the full sky coverage case, the method gives a speed-up of 2--3x in compute time compared to a simple diagonal preconditioner, and it is easier to implement in terms of practical computer code. In the case where a mask is applied and prior-driven constrained realization is sought within the mask, this is the first time full convergence has been achieved at the full resolution of the Planck dataset. Prototype benchmark code is available at https://github.com/dagss/cmbcr.
We investigate in detail the implications of the constant-roll condition on the inflationary era of a scalar field coupled to a teleparallel $f(T)$ gravity. The resulting cosmological equations constitute a reconstruction technique which enables us to find either the $f(T)$ gravity which corresponds to a given cosmological evolution, or the Hubble rate of the cosmological evolution generated by a fixed $f(T)$ gravity. We also analyze in some detail the phase space of the constant-roll teleparallel gravity and we discuss the physical significance of the resulting fixed points and trajectories. Also we calculate the observational indices of a theory with given $f(T)$ gravity, and we discuss all the implications of the constant-roll condition on these. As we demonstrate, the resulting theory can be compatible with the current observational data, for a wide range of values of the free parameters of the theory.
I discuss the role of the Higgs boson as a probe of dark matter and inflation.
Links to: arXiv, form interface, find, astro-ph, recent, 1710, contact, help (Access key information)
We present an analysis of 15 Type Ia supernovae (SNe Ia) at redshift z > 1 (9 at 1.5 < z < 2.3) recently discovered in the CANDELS and CLASH Multi-Cycle Treasury programs using WFC3 on the Hubble Space Telescope. We combine these SNe Ia with a new compilation of 1050 SNe Ia, jointly calibrated and corrected for simulated survey biases to produce accurate distance measurements. We present unbiased constraints on the expansion rate at six redshifts in the range 0.07 < z < 1.5 based only on this combined SN Ia sample. The added leverage of our new sample at z > 1.5 leads to a factor of ~3 improvement in the determination of the expansion rate at z = 1.5, reducing its uncertainty to ~20%, a measurement of H(z=1.5)/H0=2.67 (+0.83,-0.52). We then demonstrate that these six measurements alone provide a nearly identical characterization of dark energy as the full SN sample, making them an efficient compression of the SN Ia data. The new sample of SNe Ia at z > 1 usefully distinguishes between alternative cosmological models and unmodeled evolution of the SN Ia distance indicators, placing empirical limits on the latter. Finally, employing a realistic simulation of a potential WFIRST SN survey observing strategy, we forecast optimistic future constraints on the expansion rate from SNe Ia.
We present optical light curves, redshifts, and classifications for 361 spectroscopically confirmed Type Ia supernovae (SNeIa) discovered by the Pan-STARRS1 (PS1) Medium Deep Survey. We detail improvements to the PS1 SN photometry, astrometry and calibration that reduce the systematic uncertainties in the PS1 SN Ia distances. We combine the subset of 276 PS1 SN Ia ($0.03 < z < 0.65$) with distance estimates of SN Ia from SDSS, SNLS, various low-z and HST samples to form the largest combined sample of SN Ia consisting of a total of 1049 SN Ia ranging from $0.01 < z < 2.3$, which we call the 'Pantheon Sample'. Photometric calibration uncertainties have long dominated the systematic error budget of every major analysis of cosmological parameters with SNIa. We have reduced these calibration systematics to the point where they are similar in size to the other major sources of known systematic uncertainties: the nature of the intrinsic scatter of SNIa and modeling of selection effects. When combining Planck 2015 CMB measurements with the Pantheon SN sample, we find $\Omega_m = 0.306 \pm 0.012$ and $w = -1.031 \pm 0.040$ for the wCDMmodel. When the SN and CMB constraints are combined with constraints from BAO and local H0 measurements, the analysis yields the most precise measurement of dark energy to date: $w_0 = -1.011 \pm 0.087$ and $w_a = -0.215 \pm 0.402$ for the w0waCDM model. Tension with a cosmological constant previously seen in an analysis of PS1 and low-z SNe has diminished after an increase of $2\times$ in the statistics of the PS1 sample, improved calibration and photometry, and stricter light-curve quality cuts. We find the systematic uncertainties in our measurements of dark energy are almost as large as the statistical uncertainties, primarily due to limitations of modeling the low-redshift sample. This must be addressed for future progress in using SN Ia to measure dark energy.
We use nearly 1,200 supernovae (SNe) from Pan-STARRS and $\sim$200 low-$z$ ($z < 0.1$) SNe Ia to measure cosmological parameters. Though most of these SNe lack spectroscopic classifications, in a previous paper (I) we demonstrated that photometrically classified SNe can still be used to infer unbiased cosmological parameters by using a Bayesian methodology that marginalizes over core-collapse (CC) SN contamination. Our sample contains nearly twice as many SNe as the largest previous compilation of SNe Ia. Combining SNe with Cosmic Microwave Background (CMB) constraints from the Planck satellite, we measure the dark energy equation of state parameter $w$ to be -0.986$\pm$0.058 (stat$+$sys). If we allow $w$ to evolve with redshift as $w(a) = w_0 + w_a(1-a)$, we find $w_0 = -0.923 \pm 0.148$ and $w_a =$ -0.404$\pm$0.797. These results are consistent with measurements of cosmological parameters from the JLA and from a new analysis of 1049 spectroscopically confirmed SNe Ia (Scolnic et al. 2017). We try four different photometric classification priors for Pan-STARRS SNe and two alternate ways of modeling the CC SN contamination, finding that none of these variants gives a $w$ that differs by more than 1% from the baseline measurement. The systematic uncertainty on $w$ due to marginalizing over the CC SN contamination, $\sigma_w^{\textrm{CC}} = 0.019$, is approximately equal to the photometric calibration uncertainty and is lower than the systematic uncertainty in the SN\,Ia dispersion model ($\sigma_w^{\textrm{disp}} = 0.024$). Our data provide one of the best current constraints on $w$, demonstrating that samples with $\sim$5% CC SN contamination can give competitive cosmological constraints when the contaminating distribution is marginalized over in a Bayesian framework.
We present results from a set of simulations designed to constrain the weak lensing shear calibration for the Hyper Suprime-Cam (HSC) survey. These simulations include HSC observing conditions and galaxy images from the Hubble Space Telescope (HST), with fully realistic galaxy morphologies and the impact of nearby galaxies included. We find that the inclusion of nearby galaxies in the images is critical to reproducing the observed distributions of galaxy sizes and magnitudes, due to the non-negligible fraction of unrecognized blends in ground-based data, even with the excellent typical seeing of the HSC survey (0.58" in the $i$-band). Using these simulations, we detect and remove the impact of selection biases due to the correlation of weights and the quantities used to define the sample (S/N and apparent size) with the lensing shear. We quantify and remove galaxy property-dependent multiplicative and additive shear biases that are intrinsic to our shear estimation method, including a $\sim 10$ per cent-level multiplicative bias due to the impact of nearby galaxies and unrecognized blends. Finally, we check the sensitivity of our shear calibration estimates to other cuts made on the simulated samples, and find that the changes in shear calibration are well within the requirements for HSC weak lensing analysis. Overall, the simulations suggest that the weak lensing multiplicative biases in the first-year HSC shear catalog are controlled at the 1 per cent level.
We consider the prospects of probing features in the primordial power spectrum with future Cosmic Microwave Background (CMB) polarization measurements. In the scope of the inflationary scenario, such features in the spectrum can be produced by local non-smooth pieces in an inflaton potential (smooth and quasi-flat in general) which in turn may originate from fast phase transitions during inflation in other quantum fields interacting with the inflaton. They can fit some outliers in the CMB temperature power spectrum which are unaddressed within the standard inflationary ${\mathrm{\Lambda}}$CDM model. We consider Wiggly Whipped Inflation (WWI) as a theoretical framework leading to improvements in the fit to the Planck 2015 temperature and polarization data in comparison with the standard inflationary models, although not at a statistically significant level. We show that some type of features in the potential within the WWI models, leading to oscillations in the primordial power spectrum that extend to intermediate and small scales can be constrained with high confidence (at 3$\sigma$ or higher confidence level) by an instrument as the Cosmic ORigins Explorer (CORE). In order to investigate the possible confusion between inflationary features and footprints from the reionization era, we consider an extended reionization history with monotonic increase of free electrons with decrease in redshift. We discuss the present constraints on this model of extended reionization and future predictions with CORE. We also project, to what extent, this extended reionization can create confusion in identifying inflationary features in the data.
Planckian Interacting Dark Matter (PIDM) is a minimal scenario of dark matter assuming only gravitational interactions with the standard model and with only one free parameter, the PIDM mass. PIDM can be successfully produced by gravitational scattering in the thermal plasma of the Standard Model sector after inflation in the PIDM mass range from TeV up to the GUT scale, if the reheating temperature is sufficiently high. The minimal assumption of a GUT scale PIDM mass can be tested in the future by measurements of the primordial tensor-to-scalar ratio. While large primordial tensor modes would be in tension with the QCD axion as dark matter in a large mass range, it would favour the PIDM as a minimal alternative to WIMPs. Here we generalise the previously studied scalar PIDM scenario to the case of fermion, vector and tensor PIDM scenarios, and show that the phenomenology is nearly identical, independent of the spin of the PIDM. We also consider the specific realisation of the PIDM as the Kaluza Klein excitation of the graviton in orbifold compactifications of string theory, as well as in models of monodromy inflation and in Higgs inflation. Finally we discuss the possibility of indirect detection of PIDM through non-perturbative decay.
Cold ($T\sim 10^{4} \ \mathrm{K}$) gas is very commonly found in both galactic and cluster halos. There is no clear consensus on its origin. Such gas could be uplifted from the central galaxy by galactic or AGN winds. Alternatively, it could form in situ by thermal instability. Fragmentation into a multi-phase medium has previously been shown in hydrodynamic simulations to take place once $t_\mathrm{cool}/t_\mathrm{ff}$, the ratio of the cooling time to the free-fall time, falls below a threshold value. Here, we use 3D plane-parallel MHD simulations to investigate the influence of magnetic fields. We find that because magnetic tension suppresses buoyant oscillations of condensing gas, it destabilizes all scales below $l_\mathrm{A}^\mathrm{cool} \sim v_\mathrm{A} t_\mathrm{cool}$, enhancing thermal instability. This effect is surprisingly independent of magnetic field orientation or cooling curve shape, and sets in even at very low magnetic field strengths. Magnetic fields critically modify both the amplitude and morphology of thermal instability, with $\delta \rho/\rho \propto \beta^{-1/2}$, where $\beta$ is the ratio of thermal to magnetic pressure. In galactic halos, magnetic fields can render gas throughout the entire halo thermally unstable, and may be an attractive explanation for the ubiquity of cold gas, even in the halos of passive, quenched galaxies.
We present the discovery of an active galactic nucleus (AGN) that is turning off and then on again in the z=0.06 galaxy SDSS J1354+1327. This episodic nuclear activity is the result of discrete accretion events, which could have been triggered by a past interaction with the companion galaxy that is currently located 12.5 kpc away. We originally targeted SDSS J1354+1327 because its Sloan Digital Sky Survey spectrum has narrow AGN emission lines that exhibit a velocity offset of 69 km s$^{-1}$ relative to systemic. To determine the nature of the galaxy and its velocity-offset emission lines, we observed SDSS J1354+1327 with Chandra/ACIS, Hubble Space Telescope/Wide Field Camera 3, Apache Point Observatory optical longslit spectroscopy, and Keck/OSIRIS integral-field spectroscopy. We find a ~10 kpc cone of photoionized gas south of the galaxy center and a ~1 kpc semi-spherical front of shocked gas, which is responsible for the velocity offset in the emission lines, north of the galaxy center. We interpret these two outflows as the result of two separate AGN accretion events; the first AGN outburst created the southern outflow, and then $<10^5$ yrs later the second AGN outburst launched the northern shock front. SDSS J1354+1327 is the galaxy with the strongest evidence for an AGN that has turned off and then on again, and it fits into the broader context of AGN flickering that includes observations of AGN light echoes.
Knowledge of the number density of H$\alpha$ emitting galaxies is vital for assessing the scientific impact of the Euclid and WFIRST missions. In this work we present predictions from a galaxy formation model, Galacticus, for the cumulative number counts of H$\alpha$-emitting galaxies. We couple Galacticus to three different dust attenuation methods and examine the counts using each method. A $\chi^2$ minimisation approach is used to compare the model predictions to observed galaxy counts and calibrate the dust parameters. We find that weak dust attenuation is required for the Galacticus counts to be broadly consistent with the observations, though the optimum dust parameters return large values for $\chi^2$, suggesting that further calibration of Galacticus is necessary. The model predictions are also consistent with observed estimates for the optical depth and the H$\alpha$ luminosity function. Finally we present forecasts for the redshift distributions and number counts for two Euclid-like and one WFIRST-like survey. For a Euclid-like survey with redshift range $0.9\leqslant z\leqslant 1.8$ and H$\alpha+{\rm [NII]}$ blended flux limit of $2\times 10^{-16}{\rm erg}\,{\rm s}^{-1}\,{\rm cm}^{-2}$ we predict a number density between 3900--4800 galaxies per square degree. For a WFIRST-like survey with redshift range $1\leqslant z\leqslant 2$ and blended flux limit of $1\times 10^{-16}{\rm erg}\,{\rm s}^{-1}\,{\rm cm}^{-2}$ we predict a number density between 10400--15200 galaxies per square degree.
We present an empirical parameterization of the [NII]/H$\alpha$ flux ratio as a function of stellar mass and redshift valid at 0 < z < 2.7 and 8.5 < log(M) < 11.0. This description can easily be applied to (i) simulations for modeling [NII] line emission, (ii) deblend [NII] and H$\alpha$ in current low-resolution grism and narrow-band observations to derive intrinsic H$\alpha$ fluxes, and (iii) to reliably forecast the number counts of H$\alpha$ emission line galaxies for future surveys such as those planned for Euclid and WFIRST. Our model combines the evolution of the locus on the BPT diagram measured in spectroscopic data out to z ~ 2.5 with the strong dependence of [NII]/H$\alpha$ on stellar mass and [OIII]/H$\beta$ observed in local galaxy samples. We find large variations in the [NII]/H$\alpha$ flux ratio at a fixed redshift due to its dependency on stellar mass, hence the assumption of a constant [NII] contamination fraction can lead to a significant under- or over-estimate of H$\alpha$ luminosities. Specifically, measurements of the intrinsic H$\alpha$ luminosity function derived from current low-resolution grism spectroscopy assuming a constant 29% contamination of [NII] are likely over-estimated by factors of 2-4 at log(L) > 43.0 and systematically under-estimated by ~50 at log(L) < 42.5 at redshifts z ~ 1.5. This has implications on the prediction of H-alpha emitters for Euclid and WFIRST. We also study the impact of blended H$\alpha$ and [NII] on the accuracy of measured spectroscopic redshifts.
We investigate the length of the period of validity of a classical description for the cosmic axion field. To this end, we first show that we can understand the oscillating axion solution as expectation value over an underlying coherent quantum state. Once we include self-interaction of the axion, the quantum state evolves so that the expectation value over it starts to deviate from the classical solution. The time-scale of this process defines the quantum break-time. For the hypothetical dark matter axion field in our Universe, we show that quantum break-time exceeds the age of the Universe by many orders of magnitude. This conclusion is independent of specific properties of the axion model. Thus, experimental searches based on the classical approximation of the oscillating cosmic axion field are fully justified. Additionally, we point out that the distinction of classical nonlinearities and true quantum effects is crucial for calculating the quantum break-time in any system. Our analysis can also be applied to other types of dark matter that are described as classical fluids in the mean field approximation.
The progenitors of astronomical transients are linked to a specific stellar population and galactic environment, and observing their host galaxies hence constrains the physical nature of the transient itself. Here, we use imaging from the Hubble Space Telescope, and spatially-resolved, medium resolution spectroscopy from the Very Large Telescope obtained with X-Shooter and MUSE to study the host of the very luminous transient ASASSN-15lh. The dominant stellar population at the transient site is old (around 1 to 2 Gyr), without signs of recent star-formation. We also detect emission from ionized gas, originating from three different, time-invariable, narrow components of collisionally-excited metal and Balmer lines. The ratios of emission lines in the Baldwin-Phillips-Terlevich diagnostic diagram indicate that the ionization source is a weak Active Galactic Nucleus with a black hole mass of $M_\bullet = 5_{-3}^{+8}\cdot10^{8} M_\odot$, derived through the $M_\bullet$-$\sigma$ relation. The narrow line components show spatial and velocity offsets on scales of 1 kpc and 500 km/s, respectively, that we explain by gas kinematics in the narrow line region. The location of the central component, which we argue is also the position of the supermassive black hole, aligns with that of the transient within an uncertainty of 170 pc. Using this positional coincidence as well as other similarities with the hosts of Tidal Disruption Events, we strengthen the argument that the transient emission observed as ASASSN-15lh is related to the disruption of a star around a supermassive black hole, most probably spinning with a Kerr parameter $a_\bullet\gtrsim0.5$.
The global 21 cm signal from Cosmic Dawn (CD) and the Epoch of Reionization (EoR), at redshifts $z \sim 6-30$, probes the nature of first sources of radiation as well as physics of the Inter-Galactic Medium (IGM). Given that the signal is predicted to be extremely weak, of wide fractional bandwidth, and lies in a frequency range that is dominated by Galactic and Extragalactic foregrounds as well as Radio Frequency Interference, detection of the signal is a daunting task. Critical to the experiment is the manner in which the sky signal is represented through the instrument. It is of utmost importance to design a system whose spectral bandpass and additive spurious can be well calibrated and any calibration residual does not mimic the signal. SARAS is an ongoing experiment that aims to detect the global 21 cm signal. Here we present the design philosophy of the SARAS 2 system and discuss its performance and limitations based on laboratory and field measurements. Laboratory tests with the antenna replaced with a variety of terminations, including a network model for the antenna impedance, show that the gain calibration and modeling of internal additives leave no residuals with Fourier amplitudes exceeding 2~mK, or residual Gaussians of 25 MHz width with amplitudes exceeding 2~mK. Thus, even accounting for reflection and radiation efficiency losses in the antenna, the SARAS~2 system is capable of detection of complex 21-cm profiles at the level predicted by currently favoured models for thermal baryon evolution.
We study a rapidly-oscillating scalar field with potential $V(\phi) = k|\phi|^n$ nonminally coupled to the Ricci scalar $R$ via a term of the form $(1- 8 \pi G_0 \xi \phi^2) R$ in the action. In the weak coupling limit, we calculate the effect of the nonminimal coupling on the time-averaged equation of state parameter $\gamma = (p + \rho)/\rho$. The change in $\langle \gamma \rangle$ is always negative for $n \ge 2$ and always positive for $n < 0.71$ (which includes the case where the oscillating scalar field could serve as dark energy), while it can be either positive or negative for intermediate values of $n$. Constraints on the time-variation of $G$ force this change to be infinitesimally small at the present time whenever the scalar field dominates the expansion, but constraints in the early universe are not as stringent. The rapid oscillation induced in $G$ also produces an additional contribution to the Friedman equation that behaves like an effective energy density with a stiff equation of state, but we show that, under reasonable assumptions, this effective energy density is always smaller than the density of the scalar field itself.
The quantization of a spherically symmetric null shells is performed and extended to the framework of phase-space noncommutative (NC) quantum mechanics. The encountered properties are investigated making use of the Israel junction conditions on the shell, considering that it is the boundary between two spherically symmetric spacetimes. Using this method, and considering two different Kantowski-Sachs spacetimes as a representation for the Schwarzschild spacetime, the relevant quantities on the shell are computed, such as its stress-energy tensor and the action for the whole spacetime. From the obtained action, the Wheeler-deWitt equation is deduced in order to provide the quantum framework for the system. Solutions for the wavefunction of the system are found on both the commutative and NC scenarios. It is shown that, on the commutative version, the wave function has a purely oscillatory behavior in the interior of the shell. In the NC setting, it is shown that the wavefunction vanishes at the singularity, as well as, at the event horizon of the black hole.
In this work, we introduce a bottom-up $F(R)$ gravity reconstruction technique, in which we fix the observational indices and we seek for the $F(R)$ gravity which may realize them. Particularly, as an exemplification of our method, we shall assume that the scalar-to-tensor ratio has a specific form, and from it we shall reconstruct the $F(R)$ gravity that may realize it, focusing on special values of the parameters in order to obtain analytical results. The observational indices we study are compatible with the latest observational data, and we discuss how the functional form of the observational indices may affect the viability of the model.
Links to: arXiv, form interface, find, astro-ph, recent, 1710, contact, help (Access key information)
We explore a new technique to measure cosmic shear using Einstein rings as standard shapes. Einstein rings are formed when a strong lensing system has circular symmetry. These rings can become elliptical due to weak lensing by foreground structures. In Birrer et al. (2017), we showed that the detailed modelling of Einstein rings can be used to measure external shear to high precision. In this letter, we explore how a collection of Einstein rings can be used as a statistical probe of cosmic shear. We present a forecast of the cosmic shear information available in Einstein rings for different strong lensing survey configurations. We find that, assuming that the number density of Einstein rings in the COSMOS survey is representative, future strong lensing surveys should have a cosmological precision comparable to the current ground based weak lensing surveys. We discuss how this technique is complementary to the standard cosmic shear analyses since it is sensitive to different systematic and can be used for cross-calibration.
Near-infrared (IR) diffuse Galactic light (DGL) consists of scattered light and thermal emission from interstellar dust grains illuminated by interstellar radiation field (ISRF). At 1.25 and 2.2um, recent observational study shows that intensity ratios of the DGL to interstellar 100um dust emission steeply decrease toward high Galactic latitudes (b). In this paper, we investigate origin(s) of the b-dependence on the basis of models of thermal emission and scattered light. Combining a thermal emission model with regional variation of the polycyclic aromatic hydrocarbon abundance observed with Planck, we show that contribution of the near-IR thermal emission component to the observed DGL is less than ~20%. We also examine the b-dependence of the scattered light, assuming a plane-parallel Galaxy with smooth distributions of the ISRF and dust density along vertical direction, and assuming a scattering phase function according to a recently developed model of interstellar dust. We normalize the scattered light intensity to the 100um intensity corrected for deviation from the cosecant-b law according to the Planck observation. As the result, the present model taking all the b-dependence of dust and ISRF properties can account for the observed b-dependence of the near-IR DGL. However, uncertainty of the correction for the 100um emission is large and other normalizing quantities may be appropriate for more robust analysis of the DGL.
Cosmic background neutrinos have a large velocity dispersion, which causes the evolution of long-wavelength density perturbations to depend on scale. This scale-dependent growth leads to the well-known suppression in the linear theory matter power spectrum that is used to probe neutrino mass. In this paper, we study the impact of long-wavelength density perturbations on small-scale structure formation. By performing separate universe simulations where the long-wavelength mode is absorbed into the local expansion, we measure the responses of the cold dark matter (CDM) power spectrum and halo mass function, which correspond to the squeezed-limit bispectrum and halo bias. We find that the scale-dependent evolution of the long-wavelength modes causes these quantities to depend on scale and provide simple expressions to model them in terms of scale and the amount of massive neutrinos. Importantly, this scale-dependent bias reduces the suppression in the linear halo power spectrum due to massive neutrinos by 13 and 26% for objects of bias $\bar{b}=2$ and $\bar{b} \gg1$, respectively. We demonstrate with high statistical significance that the scale-dependent halo bias ${\it cannot}$ be modeled by the CDM and neutrino density transfer functions at the time when the halos are identified. This reinforces the importance of the temporal nonlocality of structure formation, especially when the growth is scale dependent.
We simulate the propagation of cosmic rays at ultra-high energies, $\gtrsim 10^{18}$ eV, in models of extragalactic magnetic fields in constrained simulations of the local Universe. We use constrained initial conditions with the cosmological magnetohydrodynamics code ENZO. The resulting models of the distribution of magnetic fields in the local Universe are used in the CRPROPA code to simulate the propagation of ultra-high energy cosmic rays. We investigate the impact of six different magneto-genesis scenarios, both primordial and astrophysical, on the propagation of cosmic rays over cosmological distances. Moreover, we study the influence of different source distributions around the Milky Way. Our study shows that different scenarios of magneto-genesis do not have a large impact on the anisotropy measurements of ultra-high energy cosmic rays. However, at high energies above the GZK-limit, there is anisotropy caused by the distribution of nearby sources, independent of the magnetic field model. This provides a chance to identify cosmic ray sources with future full-sky measurements and high number statistics at the highest energies. We further compare to the dipole signal recently reported by the Pierre Auger Collaboration. None of our source models could mimic the observed signal. Our results hint on a strong dipolar component in the distribution of sources necessary to account for that signal.
Galaxy clusters, employed by Zwicky to demonstrate the existence of dark matter, pose new stringent tests. If merging clusters demonstrate that dark matter is self-interacting with cross section $\sigma/m\sim 2$ cm$^2$/gr, MACHOs, primordial black holes and light axions that build MACHOs are ruled out as cluster dark matter. Recent strong lensing and X-ray gas data of the quite relaxed and quite spherical cluster A1835 allow to test the cases of dark matter with Maxwell-Boltzmann, Bose-Einstein and Fermi-Dirac distribution, next to Navarro-Frenck-White profiles. Fits to all these profiles are formally rejected at over $5\sigma$, except in the fermionic situation. The interpretation in terms of (nearly) Dirac neutrinos with mass of $1.61^{+0.19}_{-0.30}$ eV/$c^2$ is consistent with results on the cluster A1689, with the WMAP, Planck and DES dark matter fractions and with the nondetection of neutrinoless double $\beta$-decay. The case will be tested in the 2018 KATRIN experiment.
Following extensive tests on analytical toy-catalogues (paper III), we present the results of a more realistic study over a 711 deg2 template-based cosmological simulation. Dark matter halos from the Aardvark simulation have been ascribed luminosities, temperatures and core-radii, using local scaling relations and assuming self-similar evolution. Predicted X-ray sky-maps are converted into XMM event lists using an instrumental simulator. The XXL pipeline run on the resulting sky images, produces an 'observed' cluster catalogue over which the tests have been performed. This allowed us to investigate the relative power of various combinations of count-rate, hardness-ratio, apparent-size and redshift information. Two fitting methods were used : a traditional MCMC approach and a simple minimisation procedure (Amoeba) whose mean uncertainties are a posteriori evaluated by means of synthetic catalogues. Results were analysed and compared to the predictions from a Fisher analysis (FA). For this particular catalogue realisation, assuming perfectly known scaling relations, the CR-HR combination gives $\sigma_8$ and $\Omega_m$ at the 10% level, while CR-HR-r$_c$-z improves this to $\leq$ 3%. Adding a second hardness ratio improves results from the CR-HR1-r$_c$ combination, but to a lesser extent than adding redshift information. When all coefficients of the M-T relation (including scatter) are also fitted, cosmological parameters are constrained to within 5-10%, and larger for the M-T coefficients (up to a factor of two for the scatter). Errors returned by the MCMC, by Amoeba, and by FA predictions, are in most cases in excellent agreement and always within a factor of two. We also study the impact of the scatter of the M-Rc relation on the number of detected clusters: for the cluster typical sizes usually assumed, the larger the scatter, the lower the number of detected objects.
The current available CMB data show an anomalously low value of the CMB temperature fluctuations at large angular scales (l < 40). This lack of power is not explained by the minimal LCDM model, and one of the possible mechanisms explored in the literature to address this problem is the presence of features in the primordial power spectrum (PPS) motivated by the early universe physics. In this paper, we analyse a set of cutoff inflationary PPS models using a Bayesian model comparison approach in light of the latest Cosmic Microwave Background (CMB) data from the Planck Collaboration. Our results show that the standard power-law parameterisation is preferred over all models considered in the analysis, which motivates the search for alternative explanations for the observed lack of power in the CMB anisotropy spectrum.
We studied some statistical properties of the spatial point process displayed by GRBs of known redshift. To find ring like point patterns we developed an algorithm and defined parameters to characterize the level of compactness and regularity of the rings found in this procedure. Applying this algorithm to the GRB sample we identified three more ring like point patterns. Although, they had the same regularity but much less level of compactness than the original GRB ring. Assuming a stochastic independence of the angular and radial positions of the GRBs we obtained 1502 additional samples, altogether 542222 data points, by bootstrapping the original one. None of these data points participated in rings having similar level of compactness and regularity as the original one. Using an appropriate kernel we estimated the joint probability density of the angular and radial variables of the GRBs. Performing MCMC simulations we obtained 1502 new samples, altogether 542222 data points. Among these data points only three represented ring like patterns having similar parameters as the original one. By defining a new statistical variable we tested the independence of the angular and radial variables of the GRBs. We concluded that despite the existence of local irregularities in the GRBs' spatial distribution (e.g. the GGR) one can not reject the Cosmological Principle, based on their spatial distribution as a whole. We pointed out the large scale spatial pattern of the GRB activity does not necessarily re ects the large scale distribution of the cosmic matter.
Upcoming cosmic microwave background (CMB) surveys will soon make the first detection of the polarized Sunyaev-Zel'dovich effect, the linear polarization generated by the scattering of CMB photons on the free electrons present in collapsed objects. Measurement of this polarization along with knowledge of the electron density of the objects allows a determination of the quadrupolar temperature anisotropy of the CMB as viewed from the space-time location of the objects. Maps of these remote temperature quadrupoles have several cosmological applications. Here we propose a new application: reconstruction of the cosmological reionization history. We show that with quadrupole measurements out to redshift 3, constraints on the mean optical depth can be improved by an order of magnitude beyond the CMB cosmic variance limit.
This article provides a method for quick computation of galaxy two-point correlation function (2pCF) using Python. One of the salient features of this approach is that it can be used for calculating galaxy clustering for any arbitrary geometry (or Cosmology) model. Being efficient enough to run fast on a low-spec desktop computer, this `recipe' can be used for quick validation of alternative models and for pedagogical purposes.
By using a novel interface between the modern smoothed particle hydrodynamics code GASOLINE2 and the chemistry package KROME, we follow the hydrodynamical and chemical evolution of an isolated galaxy. In order to assess the relevance of different physical parameters and prescriptions, we constructed a suite of ten simulations, in which we vary the chemical network (primordial and metal species), how metal cooling is modelled (non-equilibrium versus equilibrium), the initial gas metallicity (from ten to hundred per cent solar), and how molecular hydrogen forms on dust. This is the first work in which metal injection from supernovae, turbulent metal diffusion, and a metal network with non-equilibrium metal cooling are self-consistently included in a galaxy simulation. We find that modelling the chemical evolution of several metal species and the corresponding non-equilibrium metal cooling has important effects on the thermodynamics of the gas, the chemical abundances, and the appearance of the galaxy: the gas is typically warmer, has a larger molecular gas mass fraction, and has a smoother disc. We also conclude that, at relatively high metallicity, the choice of molecular hydrogen formation rates on dust is not crucial. Moreover, we confirm that a higher initial metallicity produces a colder gas and a larger fraction of molecular gas, with the low-metallicity simulation best matching the observed molecular Kennicutt-Schmidt relation. Finally, our simulations agree quite well with observations which link star formation rate to metal emission lines.
We perform a general-relativistic magnetohydrodynamics simulation for $\approx 30$ ms after merger of a binary neutron star to a remnant massive neutron star (RMNS) with a high spatial resolution of the finest grid resolution $12.5$ m. First, we estimate that the Kelvin-Helmholtz instability at merger could amplify the magnetic-field energy up to $\sim 1\%$ of the thermal energy. Second, we find that the magnetorotational instability in the RMNS envelope and torus with $\rho < 10^{13}~{\rm g~cm^{-3}}$ sustains magneto-turbulent state and the effective viscous parameter in these regions is likely to converge to $\approx 0.01$--$0.02$ with respect to the grid resolution. Third, the current grid resolution is not still fine enough to sustain magneto-turbulent state in the RMNS with $\rho \ge 10^{13}~{\rm g~cm^{-3}}$.
Motivated by the observed differences in the nebular emission of nearby and high-redshift galaxies, we carry out a set of direct numerical simulations of turbulent astrophysical media exposed to a UV background. The simulations assume a metallicity of $Z/Z_{\odot}$=0.5 and explicitly track ionization, recombination, charge transfer, and ion-by-ion radiative cooling for several astrophysically important elements. Each model is run to a global steady state that depends on the ionization parameter $U$, and the one-dimensional turbulent velocity dispersion, $\sigma_{\rm 1D}$, and the turbulent driving scale. We carry out a suite of models with a T=42,000K blackbody spectrum, $n_e$ = 100 cm$^{-3}$ and $\sigma_{\rm 1D}$ ranging between 0.7 to 42 km s$^{-1},$ corresponding to turbulent Mach numbers varying between 0.05 and 2.6. We report our results as several nebular diagnostic diagrams and compare them to observations of star-forming galaxies at a redshift of $z\approx$2.5, whose higher surface densities may also lead to more turbulent interstellar media. We find that subsonic, transsonic turbulence, and turbulence driven on scales of 1 parsec or greater, have little or no effect on the line ratios. Supersonic, small-scale turbulence, on the other hand, generally increases the computed line emission. In fact with a driving scale $\approx 0.1$ pc, a moderate amount of turbulence, $\sigma_{\rm 1D}$=21-28 km s$^{-1},$ can reproduce many of the differences between high and low redshift observations without resorting to harder spectral shapes.
We present photometric redshifts for 1,031 X-ray sources in the X-ATLAS field, using the machine learning technique TPZ (Carrasco Kind & Brunner 2013). X-ATLAS covers 7.1 deg2 observed with the XMM-Newton within the Science Demonstration Phase (SDP) of the H-ATLAS field, making it one of the largest contiguous areas of the sky with both XMMNewton and Herschel coverage. All of the sources have available SDSS photometry while 810 have additionally mid-IR and/or near-IR photometry. A spectroscopic sample of 5,157 sources primarily in the XMM/XXL field, but also from several X-ray surveys and the SDSS DR13 redshift catalogue, is used for the training of the algorithm. Our analysis reveals that the algorithm performs best when the sources are split, based on their optical morphology, into point-like and extended sources. Optical photometry alone is not enough for the estimation of accurate photometric redshifts, but the results greatly improve when, at least, mid-IR photometry is added in the training process. In particular, our measurements show that the estimated photometric redshifts for the X-ray sources of the training sample, have a normalized absolute median deviation, n_mad=0.06, and the percentage of outliers, eta=10-14 percent, depending on whether the sources are extended or point-like. Our final catalogue contains photometric redshifts for 933 out of the 1,031 X-ray sources with a median redshift of 0.9.
In this study, we present a suite of high-resolution numerical simulations of an isolated galaxy to test a sub-grid framework to consistently follow the formation and dissociation of H$_2$ with non-equilibrium chemistry. The latter is solved via the package KROME, coupled to the mesh-less hydrodynamic code GIZMO. We include the effect of star formation (SF), modelled with a physically motivated prescription independent of H$_2$, supernova feedback and mass losses from low-mass stars, extragalactic and local stellar radiation, and dust and H$_2$ shielding, to investigate the emergence of the observed correlation between H$_2$ and SF rate surface densities. We present two different sub-grid models and compare them with on-the-fly radiative transfer (RT) calculations, to assess the main differences and limits of the different approaches. We also discuss a sub-grid clumping factor model to enhance the H$_2$ formation, consistent with our SF prescription, which is crucial, at the achieved resolution, to reproduce the correlation with H$_2$. We find that both sub-grid models perform very well relative to the RT simulation, giving comparable results, with moderate differences, but at much lower computational cost. We also find that, while the Kennicutt-Schmidt relation for the total gas is not strongly affected by the different ingredients included in the simulations, the H$_2$-based counterpart is much more sensitive, because of the crucial role played by the dissociating radiative flux and the gas shielding.
The globular cluster systems of galaxies are well-known to extend to large galactocentric radii. Here we quantify the size of GC systems using the half number radius of 22 GC systems around early-type galaxies from the literature. We compare GC system sizes to the sizes and masses of their host galaxies. We find that GC systems typically extend to 4$\times$ that of the host galaxy size, however this factor varies with galaxy stellar mass from about 3$\times$ for M$^{\ast}$ galaxies to 5$\times$ for the most massive galaxies in the universe. The size of a GC system scales approximately linearly with the virial radius (R$_{200}$) and with the halo mass (M$_{200}$) to the 1/3 power. The GC system of the Milky Way follows the same relations as for early-type galaxies. For Ultra Diffuse Galaxies their GC system size scales with halo mass and virial radius as for more massive, larger galaxies. UDGs indicate that the linear scaling of GC system size with stellar mass for massive galaxies flattens out for low stellar mass galaxies. Our scalings are different to those reported recently by Hudson \& Robison (2017).
I studied giant discy galaxies with optical radii more than 30 kpc. The comparison of these systems with discy galaxies of moderate sizes revealed that they tend to have higher rotation velocities, B-band luminosities, HI masses and dark-to-luminous mass ratios. The giant discs follow the trend $\log(M_{\rm HI})(R_{25})$ found for normal size galaxies. It indicates the absence of the peculiarities of evolution of star formation in these galaxies. The HI mass to luminosity ratio of giant galaxies appears not to differ from that of normal size galaxies, giving evidences in favor of similar star formation efficiency. I also found that the bars and rings occur more frequently among giant discs. I performed mass-modelling of the subsample of 18 giant galaxies with available rotation curves and surface photometry data and constructed $\chi^2$ maps for the parameters of their dark matter haloes. These estimates indicate that giant discs tend to be formed in larger more massive and rarified dark haloes in comparison to moderate size galaxies. However giant galaxies do not deviate significantly from the relations between the optical sizes and dark halo parameters for moderate size galaxies. These findings can rule out the catastrophic scenario of the formation of at least most of giant discs, since they follow the same relations as normal discy galaxies. The giant sizes of the discs can be due to the high radial scale of the dark matter haloes in which they were formed.
Cosmological observations over past couple of decades favor our universe with a tiny positive cosmological constant. Presence of cosmological constant not only imposes theoretical challenges in gravitational waves physics, it has also observational relevance. Inclusion of cosmological constant in linearized theory of gravitational waves modifies the power radiated quadrupole formula. There are two types of observations which can be impacted by the modified quadrupole formula. One is the orbital decay of an inspiraling binary and other is the modification of the waveform at the detector. Modelling a compact binary system in an elliptic orbit on de Sitter background we obtain quadrupolar power radiation and we also investigate its impact on orbital decay rate.
Higher derivative corrections are ubiquitous in effective field theories, which seemingly introduces new degrees of freedom at successive order. This is actually an artefact of the implicit local derivative expansion defining effective field theories. We argue that higher derivative corrections that introduce additional degrees of freedom should be removed and their effects captured either by lower derivative corrections, or special combinations of higher derivative corrections not propagating extra derees of freedom. Three methods adapted for this task are examined and field redefinitions are found to be most appropriate. First order higher derivative corrections in a scalar tensor theory are removed by field redefinition and it is found that their effects are captured by a subset of Horndeski theories. A case is made for restricting the effective field theory expansions in principle to only terms not introducing additional degrees of freedom.
We consider decaying hydromagnetic turbulence with initial kinetic helicity in an electrically conducting fluid and show that a weak nonhelical magnetic field eventually becomes fully helical. Already before this happens, the magnetic field undergoes inverse cascading with the magnetic energy decaying approximately like $t^{-0.5}$, which is even slower than in the fully helical case. In this parameter range, the product of magnetic energy and correlation length to the power 1.3 is approximately constant. Our result has applications to a wide range of experimental dynamos and astrophysical time-dependent plasmas, including primordial turbulence in the early universe.
The NGC 5903 galaxy group is a nearby (~30 Mpc) system of ~30 members, dominated by the giant ellipticals NGC 5903 and NGC 5898. The group contains two unusual structures, a ~110 kpc long HI filament crossing NGC 5903, and a ~75 kpc wide diffuse, steep-spectrum radio source of unknown origin which overlaps NGC 5903 and appears to be partly enclosed by the HI filament. Using a combination of Chandra, XMM-Newton, GMRT and VLA observations, we detect a previously unknown ~0.65 keV intra-group medium filling the volume within 145 kpc of NGC 5903, and find a loop of enhanced X-ray emission extending ~35 kpc southwest from the galaxy, enclosing the brightest part of the radio source. The northern and eastern parts of this X-ray structure are also strongly correlated with the southern parts of the HI filament. We determine the spectral index of the bright radio emission to be $\alpha_{150}^{612}$=1.03$\pm$0.08, indicating a radiative age >360 Myr. We discuss the origin of the correlated radio, X-ray and HI structures, either through an interaction-triggered AGN outburst with enthalpy 1.8x10$^{57}$ erg, or via a high-velocity collision between a galaxy and the HI filament. While neither scenario provides a complete explanation, we find that an AGN outburst is the most likely source of the principal X-ray and radio structures. However, it is clear that galaxy interactions continue to play an important role in the development of this relatively highly evolved galaxy group. We also resolve the question of whether the group member galaxy ESO 514-3 hosts a double-lobed radio source, confirming that the source is a superposed background AGN.
Galaxy clusters are the largest virialized structures in the observable
Universe. The knowledge of their properties provides many useful astrophysical
and cosmological information.
Our aim is to derive the luminosity and stellar mass profiles of the nearby
galaxy clusters of the Omega-WINGS survey and to study the main scaling
relations valid for such systems.
We have merged the data of the WINGS and Omega-WINGS databases, sorted the
sources according to the distance from the brightest cluster galaxy (BCG) and
calculated the integrated luminosity profiles in the $B$ and $V$ bands, taking
into account extinction, photometric and spatial completeness, K-correction and
background contribution. Then, by exploiting the spectroscopic sample we
derived the stellar mass profiles of the clusters.
We got the luminosity profiles of 46 galaxy clusters, reaching $r_{200}$ in
30 cases, and the stellar mass profiles of 42 of our objects. We successfully
fitted all the integrated luminosity growth profiles with one or two embedded
S\'ersic components, deriving the main clusters parameters.
Finally, we checked the main scaling relation among the clusters parameters
in comparison with those obtained for a selected sample of early-type galaxies
(ETGs) of the same clusters.
We found that the nearby galaxy clusters are non-homologous structures like
ETGs and exhibit a color-magnitude (CM) red-sequence relation very similar to
that observed for galaxies in clusters. These properties are not expected in
the current cluster formation scenarios. In particular the existence of a CM
relation for clusters, shown here for the first time, suggests that the
baryonic structures grow and evolve in a similar way at all scales.
Links to: arXiv, form interface, find, astro-ph, recent, 1710, contact, help (Access key information)
We present a systematic comparison of several existing and new void finding algorithms, focusing on their potential power to test a particular class of modified gravity models - chameleon $f(R)$ gravity. These models deviate from standard General Relativity (GR) more strongly in low-density regions and thus voids are a promising venue to test them. We use Halo Occupation Distribution (HOD) prescriptions to populate haloes with galaxies, and tune the HOD parameters such that the galaxy two-point correlation functions are the same in both f(R) and GR models. We identify both 3D voids as well as 2D underdensities in the plane-of-the-sky to find the same void abundance and void galaxy number density profiles across all models, which suggests that they do not contain much information beyond galaxy clustering. However, the underlying void dark matter density profiles are significantly different, with f(R) voids being more underdense than GR ones, which leads to f(R) voids having a larger tangential shear signal than their GR analogues. We investigate the potential of each void finder to test f(R) models with near-future lensing surveys such as EUCLID and LSST. The 2D voids have the largest power to probe f(R) gravity, with a LSST analysis of tunnel lensing distinguishing at 65 and 10$\sigma$ (statistical error) f(R) models with $|f_{R0}|=10^{-5}$ and $10^{-6}$ from GR.
We make progress towards an analytical understanding of the regime of validity of perturbation theory for large scale structures and the nature of some non-perturbative corrections. We restrict ourselves to 1D gravitational collapse, for which exact solutions before shell crossing are known. We review the convergence of perturbation theory for the power spectrum, recently proven by McQuinn and White, and extend it to non-Gaussian initial conditions and the bispectrum. In contrast, we prove that perturbation theory diverges for the real space two-point correlation function and for the probability density function (PDF) of the density averaged in cells and all the cumulants derived from it. We attribute these divergences to the statistical averaging intrinsic to cosmological observables, which, even on very large and "perturbative" scales, gives non-vanishing weight to all extreme fluctuations. Finally, we discuss some general properties of non-perturbative effects in real space and Fourier space.
We investigate the impact of a variety of analysis assumptions that influence
cluster identification and location on the kSZ pairwise momentum signal and
covariance estimation. Photometric and spectroscopic galaxy tracers from SDSS,
WISE, and DECaLs, spanning redshifts $0.05<z<0.7$, are considered in
combination with CMB data from Planck and WMAP. With two complementary
techniques, analytic offset modeling and direct comparisons of redMaPPer
brightest and central catalog samples, we find that miscentering uncertainties
average to $0.4-0.7\sigma$ for the Planck kSZ statistical error budget obtained
with a jackknife estimator. We also find that jackknife covariance estimates
are significantly more conservative than those obtained by CMB rotation
methods. Using redMaPPer data, we concurrently compare the impact of
photometric redshift errors and miscentering. At separations $<\sim 50$ Mpc,
where the kSZ signal is largest, miscentering uncertainties can be comparable
to JK errors, while photometric redshifts are lower but still significant.
For the next generation of CMB and LSS surveys the statistical and
photometric errors will shrink markedly. Our results demonstrate that
uncertainties introduced through using galaxy proxies for cluster locations
will need to be fully incorporated, and actively mitigated, for the kSZ to
reach its full potential as a cosmological constraining tool for dark energy
and neutrino physics.
We obtain novel closed form solutions to the Friedmann equation for cosmological models containing a component whose equation of state is that of radiation $(w=1/3)$ at early times and that of cold pressureless matter $(w=0)$ at late times. The equation of state smoothly transitions from the early to late-time behavior and exactly describes the evolution of a species with a Dirac Delta function distribution in momentum magnitudes $|\vec{p}_0|$ (i.e. all particles have the same $|\vec{p}_0|$). Such a component, here termed "hot matter", is an approximate model for both neutrinos and warm dark matter. We consider it alone and in combination with cold matter and with radiation, also obtaining closed-form solutions for the growth of super-horizon perturbations in each case. The idealized model recovers $t(a)$ to better than $1.5\%$ accuracy for all $a$ relative to a Fermi-Dirac distribution (as describes neutrinos). We conclude by adding the second moment of the distribution to our exact solution and then generalizing to include all moments of an arbitrary momentum distribution in a closed form solution.
We perform three-dimensional simulations of structure formation in the early Universe, when boosting the primordial power spectrum on approximately kpc scales. We demonstrate that our simulations are capable of producing power-law profiles close to the steep $\rho\propto r^{-9/4}$ halo profiles that are commonly assumed to be a good approximation to ultracompact minihalos (UCMHs). However, we show that for more realistic initial conditions in which halos are neither perfectly symmetric nor isolated, the steep power-law profile is disrupted and we find that the Navarro-Frenk-White profile is a better fit to most halos. In the presence of background fluctuations even extreme, nearly spherical initial conditions do not remain exceptional. Nonetheless, boosting the amplitude of initial fluctuations causes all structures to form earlier and thus at larger densities. With sufficiently large amplitude of fluctuations we find that values for the concentration of typical halos in our simulations can become very large. However, despite the signal coming from dark matter annihilation inside the cores of these halos being enhanced, it is still orders-of-magnitude smaller compared to the usually assumed UCMH profile. The upper bound on the primordial power spectrum from the non-observation of UCMHs should therefore be re-evaluated.
One often hears that the strong CP problem is {\em the} one problem which cannot be solved by anthropic reasoning. We argue that this is not so. Due to nonperturbative dynamics, states with a different CP violating paramenter $\theta$ acquire different vacuum energies after the QCD phase transition. These add to the total variation of the cosmological constant in the putative landscape of Universes. An interesting possibility arises when the cosmological constant is mostly cancelled by the membrane nucleation mechanism. If the step size in the resulting discretuum of cosmological constants, $\Delta \Lambda$, is in the interval $({\rm meV})^4 < \Delta \Lambda < (100 \, {\rm MeV})^4$, the cancellation of vacuum energy can be assisted by the scanning of $\theta$. For $({\rm meV})^4 < \Delta \Lambda < ({\rm keV})^4$ this yields $\theta < 10^{-10}$, meeting the observational limits. This scenario opens up 24 orders of magnitude of acceptable parameter space for $\Delta \Lambda$ compared membrane nucleation acting alone. In such a Universe one does not need a light axion to solve the strong CP problem.
We study the effect of anisotropic radiation illumination on the alignment of polycyclic aromatic hydrocarbons (PAHs) and report that cross section mechanism of alignment earlier considered in terms of gas-grain interactions can also be efficient for the photon-grain interaction. We demonstrate this by first calculate the angle-dependence rotational damping and excitation coefficients by photon absorption followed by infrared emission and average the coefficients over internal thermal fluctuations. Then, we calculate the degree of PAH alignment for the different environments and physical parameters, including the illumination direction, ionization fraction, and magnetic field strength. For the reflection nebula (RN) conditions with unidirectional radiation field, we find that the degree of alignment tends to increase with increasing the angle $\psi$ between the illumination direction and the magnetic field, as a result of the decrease of the cross-section of photon absorption with $\psi$. We calculate the polarization of spinning PAH emission using the obtained degree of alignment for the different physical parameters. We find that the polarization of spinning PAH emission from RN can be large, between $5-20\%$ for frequency $\nu > 20$GHz, whereas the polarization is less than $3\%$ for photodissociation regions (PDRs) with much higher gas density. The polarization for the diffuse cold neutral medium (CNM) is rather low, below $1\%$ for $\nu>20$GHz, consistent with observations by WMAP and Planck. Our results demonstrate that the RNe are the favored environment to observe the polarization of spinning dust emission as well as polarized mid-IR emission from PAHs.
We extend the idea of unimodular gravity to the modified $f(R,T)$ theories. A new class of cosmological solutions, that the unimodular constraint on the metric imposes on the $f(R,T)$ theories, are studied. This extension is done in both Jordan and Einstein frames. We show that while the Lagrange multiplier (that imposes the unimodular constraint on the action) depends on the cosmic time in Jordan frame and therefore, can act as an evolving scalar field in the universe history, in the Einstein frame it acts as a cosmological constant. Then a general reconstruction method is used to realize an explicit form of the unimodular $f(R,T)$ corresponding to a given cosmological solution. By adopting a specific form of $f(R,T)$, the issue of cosmological inflation is studied in this setup. To see the observational viability of this model, a numerical analysis on the model parameter space is done in the background of Planck2015 observational data.
We present mass-loss predictions from Monte Carlo radiative transfer models for helium (He) stars as a function of stellar mass, down to 2 Msun. Our study includes both massive Wolf-Rayet (WR) stars and low-mass He stars that have lost their envelope through interaction with a companion. For these low-mass He-stars we predict mass-loss rates that are an order of magnitude smaller than by extrapolation of empirical WR mass-loss rates. Our lower mass-loss rates make it harder for these elusive stripped stars to be discovered via line emission, and we should attempt to find them through alternative methods instead. Moreover, lower mass-loss rates will make it less likely that low-mass He stars provide stripped-envelope supernovae (SNe) of type Ibc. We express our mass-loss predictions as a function of L and Z, and not as a function of the He abundance, as we do not consider this physically astute given our earlier work. The exponent of the dM/dt vs. Z dependence is found to be 0.61, which is less steep than relationships derived from recent empirical atmospheric modelling. Our shallower exponent will make it more challenging to produce "heavy" black holes of order 40 Msun, as recently discovered in the gravitational wave event GW 150914, making low metallicity for these types of events even more necessary.
Links to: arXiv, form interface, find, astro-ph, recent, 1710, contact, help (Access key information)