Showing 1–50 of 66 results for author: Wojtak, R

Consistent extinction model for type Ia supernovae in Cepheid-based calibration galaxies and its impact on $H_{0}$

Authors: Radosław Wojtak, Jens Hjorth

Abstract: The most recent SH0ES measurement of the Hubble constant, based on type Ia supernovae from the Pantheon+ compilation, employs corrections of supernova peak magnitudes which effectively accounts for extinction in the supernova host galaxies. These corrections are estimated using a probabilistic model which is trained on Hubble flow (z>0.03) supernovae and extrapolated to the calibration galaxies (t… ▽ More The most recent SH0ES measurement of the Hubble constant, based on type Ia supernovae from the Pantheon+ compilation, employs corrections of supernova peak magnitudes which effectively accounts for extinction in the supernova host galaxies. These corrections are estimated using a probabilistic model which is trained on Hubble flow (z>0.03) supernovae and extrapolated to the calibration galaxies (those with observed Cepheid distances). By comparing the corrected peak magnitudes to distance moduli from Cepheids, we show that this standard approach underestimates the brightness of reddened supernovae in the high stellar-mass ($M_{\star}>10^{10}M_{\odot}$) calibration galaxies. This can be traced back to the fact that for these galaxies, a low total-to-selective extinction coefficient (R_B~3) is assumed, while for the low stellar-mass analogues a more standard R_B~4 is assumed. We propose a minimalistic modification of the Pantheon+ extinction model in order to alleviate this systematic effect. The modification is twofold and it involves: (i) the same, Milky Way-like distribution of R_B in all calibration galaxies (with mean R_B of 4.3 -- consistent with the extinction curve used for colour corrections of the Cepheids -- and scatter 0.4) and (ii) a modified shape of the E(B-V) reddening distribution while keeping the same effective slope of the supernova peak magnitude-colour relation and the same mean E(B-V) reddening as measured for supernovae in the Hubble flow. We show that this new approach yields a significantly better fit ($Δ$BIC=-11) to the calibration data and results in a lower value of the derived Hubble constant through a stronger extinction correction of supernovae in the calibration galaxies. Our result is $H_{0}=70.5\pm1$ km/s/Mpc implying a reduction of the tension with the Planck $H_{0}$ measurement assuming a flat LCDM cosmology from $5.2σ$ to $2.8σ$. △ Less

Submitted 15 March, 2024; originally announced March 2024.

Comments: 16 pages, 7 figures, 2 tables; submitted to MNRAS

Detecting strongly-lensed type Ia supernovae with LSST

Authors: Nikki Arendse, Suhail Dhawan, Ana Sagués Carracedo, Hiranya V. Peiris, Ariel Goobar, Radek Wojtak, Catarina Alves, Rahul Biswas, Simon Huber, Simon Birrer, The LSST Dark Energy Science Collaboration

Abstract: Strongly-lensed supernovae are rare and valuable probes of cosmology and astrophysics. Upcoming wide-field time-domain surveys, such as the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST), are expected to discover an order-of-magnitude more lensed supernovae than have previously been observed. In this work, we investigate the cosmological prospects of lensed type Ia supernovae (… ▽ More Strongly-lensed supernovae are rare and valuable probes of cosmology and astrophysics. Upcoming wide-field time-domain surveys, such as the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST), are expected to discover an order-of-magnitude more lensed supernovae than have previously been observed. In this work, we investigate the cosmological prospects of lensed type Ia supernovae (SNIa) in LSST by quantifying the expected annual number of detections, the impact of stellar microlensing, follow-up feasibility, and how to best separate lensed and unlensed SNIa. We simulate SNIa lensed by galaxies, using the current LSST baseline v3.0 cadence, and find an expected number of 44 lensed SNIa detections per year. Microlensing effects by stars in the lens galaxy are predicted to lower the lensed SNIa detections by $\sim 8 \%$. The lensed events can be separated from the unlensed ones by jointly considering their colours and peak magnitudes. We define a `gold sample' of $\sim 10$ lensed SNIa per year with time delay $> 10$ days, $> 5$ detections before light-curve peak, and sufficiently bright ($m_i < 22.5$ mag) for follow-up observations. In three years of LSST operations, such a sample is expected to yield a $1.5\%$ measurement of the Hubble constant. △ Less

Submitted 16 May, 2024; v1 submitted 7 December, 2023; originally announced December 2023.

Comments: Accepted for publication in MNRAS; 15 pages, 13 figures

Measuring the Hubble constant with kilonovae using the Expanding Photosphere Method

Authors: Albert Sneppen, Darach Watson, Dovi Poznanski, Oliver Just, Andreas Bauswein, Radosław Wojtak

Abstract: While gravitational wave (GW) standard sirens from neutron star (NS) mergers have been proposed to offer good measurements of the Hubble constant, we show in this paper how a variation of the expanding photosphere method (EPM) or spectral-fitting expanding atmosphere method, applied to the kilonovae (KNe) associated with the mergers, can provide an independent distance measurement to individual me… ▽ More While gravitational wave (GW) standard sirens from neutron star (NS) mergers have been proposed to offer good measurements of the Hubble constant, we show in this paper how a variation of the expanding photosphere method (EPM) or spectral-fitting expanding atmosphere method, applied to the kilonovae (KNe) associated with the mergers, can provide an independent distance measurement to individual mergers that is potentially accurate to within a few percent. There are four reasons why the KN-EPM overcomes the major uncertainties commonly associated with this method in supernovae: 1) the early continuum is very well-reproduced by a blackbody spectrum, 2) the dilution effect from electron scattering opacity is likely negligible, 3) the explosion times are exactly known due to the GW detection, and 4) the ejecta geometry is, at least in some cases, highly spherical and can be constrained from line-shape analysis. We provide an analysis of the early VLT/X-shooter spectra AT2017gfo showing how the luminosity distance can be determined, and find a luminosity distance of $D_L = 44.5\pm0.8$ Mpc in agreement with, but more precise than, previous methods. We investigate the dominant systematic uncertainties, but our simple framework, which assumes a blackbody photosphere, does not account for the full time-dependent three-dimensional radiative transfer effects, so this distance should be treated as preliminary. The luminosity distance corresponds to an estimated Hubble constant of $H_0 = 67.0\pm 3.6$ km s$^{-1}$ Mpc$^{-1}$, where the dominant uncertainty is due to the modelling of the host peculiar velocity. We also estimate the expected constraints on $H_0$ from future KN-EPM-analysis with the upcoming O4 and O5 runs of the LIGO collaboration GW-detectors, where five to ten similar KNe would yield 1\% precision cosmological constraints. △ Less

Submitted 25 August, 2023; v1 submitted 21 June, 2023; originally announced June 2023.

Comments: A&A, In press

Journal ref: A&A 678, A14 (2023)

A new derivation of the Hubble constant from $γ$-ray attenuation using improved optical depths for the Fermi and CTA era

Authors: A. Domínguez, P. Østergaard Kirkeberg, R. Wojtak, A. Saldana-Lopez, A. Desai, J. R. Primack, J. Finke, M. Ajello, P. G. Pérez-González, V. S. Paliya, D. Hartmann

Abstract: We present $γ$-ray optical-depth calculations from a recently published extragalactic background light (EBL) model built from multiwavelength galaxy data from the Hubble Space Telescope Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (HST/CANDELS). CANDELS gathers one of the deepest and most complete observations of stellar and dust emissions in galaxies. This model resulted in a ro… ▽ More We present $γ$-ray optical-depth calculations from a recently published extragalactic background light (EBL) model built from multiwavelength galaxy data from the Hubble Space Telescope Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (HST/CANDELS). CANDELS gathers one of the deepest and most complete observations of stellar and dust emissions in galaxies. This model resulted in a robust derivation of the evolving EBL spectral energy distribution up to $z\sim 6$, including the far-infrared peak. Therefore, the optical depths derived from this model will be useful for determining the attenuation of $γ$-ray photons coming from high-redshift sources, such as those detected by the Large Area Telescope on board the Fermi Gamma-ray Space Telescope, and for multi-TeV photons that will be detected from nearby sources by the future Cherenkov Telescope Array. From these newly calculated optical depths, we derive the cosmic $γ$-ray horizon and also measure the expansion rate and matter content of the Universe including an assessment of the impact of the EBL uncertainties. We find $H_{0}=61.9$ $^{+2.9}_{-2.4}$ km s$^{-1}$ Mpc$^{-1}$ when fixing $Ω_{m}=0.32$, and $H_{0}=65.6$ $^{+5.6}_{-5.0}$ km s$^{-1}$ Mpc$^{-1}$ and $Ω_{m}=0.19\pm 0.07$, when exploring these two parameters simultaneously. △ Less

Submitted 3 November, 2023; v1 submitted 16 June, 2023; originally announced June 2023.

Comments: 11 pages, 8 figures, 1 tables; Accepted by MNRAS

Report number: IPARCOS-UCM-23-040

Two-population Bayesian hierarchical model of type Ia supernovae

Authors: Radosław Wojtak, Jens Hjorth, Jacob Osman Hjortlund

Abstract: The currently used standardisation of type Ia supernovae results in Hubble residuals whose physical origin is unaccounted for. Here, we present a complete physical interpretation of the Hubble residuals based on a novel Bayesian hierarchical model of type Ia supernovae in which latent variables describing intrinsic and extrinsic (dust related) supernova properties originate from two populations. F… ▽ More The currently used standardisation of type Ia supernovae results in Hubble residuals whose physical origin is unaccounted for. Here, we present a complete physical interpretation of the Hubble residuals based on a novel Bayesian hierarchical model of type Ia supernovae in which latent variables describing intrinsic and extrinsic (dust related) supernova properties originate from two populations. Fitting the model to SALT2 light curve parameters of supernovae in the Hubble flow we find strong ($4σ$) evidence for the presence of two overlapping, but distinct, populations differentiated primarily by their mean SALT2 shape parameter (stretch) x_1. Supernovae from the population with predominantly slow decliners (higher average x_1) are found to be intrinsically bluer (mean SALT2 colour c=-0.11) and twice as reddened by dust (mean reddening E(B-V)=0.10) than those from the opposite population dominated by fast decliners (lower average x_1) with c=-0.04$ and E(B-V)=0.05. The inferred extinction coefficient R_B in both supernova populations follows a broad (scatter 0.9) distribution with a mean of 4.1, which coincides closely with the mean extinction law measured in the Milky Way. We also find that the supernova data favour a peaked (two-tailed) distribution of selective extinction E(B-V) over the commonly adopted exponential model (one-tailed with maximum fixed at 0). Our model provides a complete explanation of the distribution of supernova light curve parameters in terms of extinction properties and the above-mentioned differences between the two populations, without the need for introducing any intrinsic scatter. △ Less

Submitted 28 August, 2023; v1 submitted 2 February, 2023; originally announced February 2023.

Comments: 19 pages, 6 figures, 3 tables; accepted for publication in MNRAS

The Young Supernova Experiment Data Release 1 (YSE DR1): Light Curves and Photometric Classification of 1975 Supernovae

Authors: P. D. Aleo, K. Malanchev, S. Sharief, D. O. Jones, G. Narayan, R. J. Foley, V. A. Villar, C. R. Angus, V. F. Baldassare, M. J. Bustamante-Rosell, D. Chatterjee, C. Cold, D. A. Coulter, K. W. Davis, S. Dhawan, M. R. Drout, A. Engel, K. D. French, A. Gagliano, C. Gall, J. Hjorth, M. E. Huber, W. V. Jacobson-Galán, C. D. Kilpatrick, D. Langeroodi , et al. (58 additional authors not shown)

Abstract: We present the Young Supernova Experiment Data Release 1 (YSE DR1), comprised of processed multi-color Pan-STARRS1 (PS1) griz and Zwicky Transient Facility (ZTF) gr photometry of 1975 transients with host-galaxy associations, redshifts, spectroscopic/photometric classifications, and additional data products from 2019 November 24 to 2021 December 20. YSE DR1 spans discoveries and observations from… ▽ More We present the Young Supernova Experiment Data Release 1 (YSE DR1), comprised of processed multi-color Pan-STARRS1 (PS1) griz and Zwicky Transient Facility (ZTF) gr photometry of 1975 transients with host-galaxy associations, redshifts, spectroscopic/photometric classifications, and additional data products from 2019 November 24 to 2021 December 20. YSE DR1 spans discoveries and observations from young and fast-rising supernovae (SNe) to transients that persist for over a year, with a redshift distribution reaching z~0.5. We present relative SN rates from YSE's magnitude- and volume-limited surveys, which are consistent with previously published values within estimated uncertainties for untargeted surveys. We combine YSE and ZTF data, and create multi-survey SN simulations to train the ParSNIP and SuperRAENN photometric classification algorithms; when validating our ParSNIP classifier on 472 spectroscopically classified YSE DR1 SNe, we achieve 82% accuracy across three SN classes (SNe Ia, II, Ib/Ic) and 90% accuracy across two SN classes (SNe Ia, core-collapse SNe). Our classifier performs particularly well on SNe Ia, with high (>90%) individual completeness and purity, which will help build an anchor photometric SNe Ia sample for cosmology. We then use our photometric classifier to characterize our photometric sample of 1483 SNe, labeling 1048 (~71%) SNe Ia, 339 (~23%) SNe II, and 96 (~6%) SNe Ib/Ic. YSE DR1 provides a training ground for building discovery, anomaly detection, and classification algorithms, performing cosmological analyses, understanding the nature of red and rare transients, exploring tidal disruption events and nuclear variability, and preparing for the forthcoming Vera C. Rubin Observatory Legacy Survey of Space and Time. △ Less

Submitted 21 February, 2023; v1 submitted 14 November, 2022; originally announced November 2022.

Comments: Accepted to ApJS; 64 pages; 35 figures; 10 tables

LensWatch: I. Resolved HST Observations and Constraints on the Strongly-Lensed Type Ia Supernova 2022qmx ("SN Zwicky")

Authors: J. D. R. Pierel, N. Arendse, S. Ertl, X. Huang, L. A. Moustakas, S. Schuldt, A. J. Shajib, Y. Shu, S. Birrer, M. Bronikowski, J. Hjorth, S. H. Suyu, S. Agarwal, A. Agnello, A. S. Bolton, S. Chakrabarti, C. Cold, F. Courbin, J. M. Della Costa, S. Dhawan, M. Engesser, O. D. Fox, C. Gall, S. Gomez, A. Goobar , et al. (17 additional authors not shown)

Abstract: Supernovae (SNe) that have been multiply-imaged by gravitational lensing are rare and powerful probes for cosmology. Each detection is an opportunity to develop the critical tools and methodologies needed as the sample of lensed SNe increases by orders of magnitude with the upcoming Vera C. Rubin Observatory and Nancy Grace Roman Space Telescope. The latest such discovery is of the quadruply-image… ▽ More Supernovae (SNe) that have been multiply-imaged by gravitational lensing are rare and powerful probes for cosmology. Each detection is an opportunity to develop the critical tools and methodologies needed as the sample of lensed SNe increases by orders of magnitude with the upcoming Vera C. Rubin Observatory and Nancy Grace Roman Space Telescope. The latest such discovery is of the quadruply-imaged Type Ia SN 2022qmx (aka, "SN Zwicky"; Goobar et al. 2022) at z = 0.3544. SN Zwicky was discovered by the Zwicky Transient Facility (ZTF) in spatially unresolved data. Here we present follow-up Hubble Space Telescope observations of SN Zwicky, the first from the multi-cycle "LensWatch" program (www.lenswatch.org). We measure photometry for each of the four images of SN Zwicky, which are resolved in three WFC3/UVIS filters (F475W, F625W, F814W) but unresolved with WFC3/IR F160W, and produce an analysis of the lensing system using a variety of independent lens modeling methods. We find consistency between time delays estimated with the single epoch of HST photometry and the lens model predictions constrained through the multiple image positions, with both inferring time delays of <1 day. Our lens models converge to an Einstein radius of (0.168+0.009-0.005)", the smallest yet seen in a lensed SN. The "standard candle" nature of SN Zwicky provides magnification estimates independent of the lens modeling that are brighter by ~1.5 mag and ~0.8 mag for two of the four images, suggesting significant microlensing and/or additional substructure beyond the flexibility of our image-position mass models. △ Less

Submitted 8 November, 2022; v1 submitted 7 November, 2022; originally announced November 2022.

Relative intrinsic scatter in hierarchical Type Ia supernova siblings analyses: Application to SNe 2021hpr, 1997bq & 2008fv in NGC 3147

Authors: Sam M. Ward, Stephen Thorp, Kaisey S. Mandel, Suhail Dhawan, David O. Jones, Kirsty Taggart, Ryan J. Foley, Gautham Narayan, Kenneth C. Chambers, David A. Coulter, Kyle W. Davis, Thomas de Boer, Kaylee de Soto, Nicholas Earl, Alex Gagliano, Hua Gao, Jens Hjorth, Mark E. Huber, Luca Izzo, Danial Langeroodi, Eugene A. Magnier, Peter McGill, Armin Rest, César Rojas-Bravo, Radosław Wojtak

Abstract: We present Young Supernova Experiment $grizy$ photometry of SN 2021hpr, the third Type Ia supernova sibling to explode in the Cepheid calibrator galaxy, NGC 3147. Siblings are useful for improving SN-host distance estimates, and investigating the contributions towards the SN Ia intrinsic scatter (post-standardisation residual scatter in distance estimates). We thus develop a principled Bayesian fr… ▽ More We present Young Supernova Experiment $grizy$ photometry of SN 2021hpr, the third Type Ia supernova sibling to explode in the Cepheid calibrator galaxy, NGC 3147. Siblings are useful for improving SN-host distance estimates, and investigating the contributions towards the SN Ia intrinsic scatter (post-standardisation residual scatter in distance estimates). We thus develop a principled Bayesian framework for analyzing SN Ia siblings. At its core is the cosmology-independent relative intrinsic scatter parameter, $σ_{Rel}$: the dispersion of siblings distance estimates relative to one another within a galaxy. It quantifies the contribution towards the total intrinsic scatter, $σ_0$, from within-galaxy variations about the siblings' common properties. It also affects the combined-distance uncertainty. We present analytic formulae for computing a $σ_{Rel}$-posterior from individual siblings distances (estimated using any SN-model). Applying a newly trained BayeSN model, we fit the light curves of each sibling in NGC 3147 individually, to yield consistent distance estimates. However, the wide $σ_{Rel}$-posterior means $σ_{Rel}\approxσ_0$ is not ruled out. We thus combine the distances by marginalizing over $σ_{Rel}$ with an informative prior: $σ_{Rel}\sim U(0,σ_0)$. Simultaneously fitting the trio's light curves improves constraints on distance, and each sibling's individual dust parameters, compared to individual fits. Higher correlation also tightens dust parameter constraints. Therefore, $σ_{Rel}$-marginalization yields robust estimates of siblings distances for cosmology, and dust parameters for siblings-host correlation studies. Incorporating NGC 3147's Cepheid-distance yields $H_0=78.4\pm 6.5\,$km/s/Mpc. Our work motivates analyses of homogeneous siblings samples, to constrain $σ_{Rel}$, and its SN-model dependence. △ Less

Submitted 1 September, 2023; v1 submitted 21 September, 2022; originally announced September 2022.

Comments: Accepted for publication in ApJ. 22 pages, 11 figures, 8 tables. Single-galaxy siblings analysis

A fast rising tidal disruption event from a candidate intermediate mass black hole

Authors: C. R. Angus, V. F. Baldassare, B. Mockler, R. J. Foley, E. Ramirez-Ruiz, S. I. Raimundo, K. D. French, K. Auchettl, H. Pfister, C. Gall, J. Hjorth, M. R. Drout, K. D. Alexander, G. Dimitriadis, T. Hung, D. O. Jones, A. Rest, M. R. Siebert, K. Taggart, G. Terreran, S. Tinyanont, C. M. Carroll, L. DeMarchi, N. Earl, A. Gagliano , et al. (14 additional authors not shown)

Abstract: Massive black holes (BHs) at the centres of massive galaxies are ubiquitous. The population of BHs within dwarf galaxies, on the other hand, is evasive. Dwarf galaxies are thought to harbour BHs with proportionally small masses, including intermediate mass BHs, with masses $10^{2} < M_{BH} < 10^{6} M_{\odot}$. Identification of these systems has historically relied upon the detection of light emit… ▽ More Massive black holes (BHs) at the centres of massive galaxies are ubiquitous. The population of BHs within dwarf galaxies, on the other hand, is evasive. Dwarf galaxies are thought to harbour BHs with proportionally small masses, including intermediate mass BHs, with masses $10^{2} < M_{BH} < 10^{6} M_{\odot}$. Identification of these systems has historically relied upon the detection of light emitted from accreting gaseous discs close to the BHs. Without this light, they are difficult to detect. Tidal disruption events (TDEs), the luminous flares produced when a star strays close to a BH and is shredded, are a direct way to probe massive BHs. The rise times of these flares theoretically correlate with the BH mass. Here we present AT2020neh, a fast rising TDE candidate, hosted by a dwarf galaxy. AT2020neh can be described by the tidal disruption of a main sequence star by a 10$^{4.7} - 10^{5.9} M_{\odot}$ BH. We find the observable rate of fast rising nuclear transients like AT2020neh to be rare, at $\lesssim 2 \times 10^{-8}$ events Mpc$^{-3}$ yr$^{-1}$. Finding non-accreting BHs in dwarf galaxies is important to determine how prevalent BHs are within these galaxies, and constrain models of BH formation. AT2020neh-like events may provide a galaxy-independent method of measuring IMBH masses. △ Less

Submitted 5 September, 2022; v1 submitted 31 August, 2022; originally announced September 2022.

Comments: Accepted for publication in Nature Astronomy

Intrinsic tension in the supernova sector of the local Hubble constant measurement and its implications

Authors: Radosław Wojtak, Jens Hjorth

Abstract: We reanalyse observations of type Ia supernovae (SNe) and Cepheids used in the local determination of the Hubble constant and find strong evidence that SN standardisation in the calibration sample (galaxies with observed Cepheids) require a steeper slope of the colour correction than in the cosmological sample (galaxies in the Hubble flow). The colour correction in the calibration sample is consis… ▽ More We reanalyse observations of type Ia supernovae (SNe) and Cepheids used in the local determination of the Hubble constant and find strong evidence that SN standardisation in the calibration sample (galaxies with observed Cepheids) require a steeper slope of the colour correction than in the cosmological sample (galaxies in the Hubble flow). The colour correction in the calibration sample is consistent with being entirely due to an extinction correction due to dust with properties similar to that of the Milky Way (R_B~4.6+/-0.4) and there is no evidence for intrinsic scatter in the SN peak magnitudes. An immediate consequence of this finding is that the local measurement of the Hubble constant becomes dependent on the choice of SN reference colour, i.e., the colour of an unreddened SN. Specifically, the Hubble constant inferred from the same observations decreases gradually with the reference colour assumed in the SN standardisation. We recover the Hubble constant measured by SH0ES for the standard choice of reference colour (SALT2 colour parameter c=0) while for a reference colour which coincides with the blue end of the observed SN colour distribution (c~-0.13), the Hubble constant from Planck observations of the CMB (assuming a flat LCDM cosmological model) is recovered. These results are intriguing in that they may provide an avenue for resolving the Hubble tension. However, since there is no obvious physical basis for the differences in colour corrections in the two SN samples, the origin of these require further investigations. △ Less

Submitted 7 July, 2022; v1 submitted 16 June, 2022; originally announced June 2022.

Comments: 12 pages, 4 figures, 3 tables; accepted for publication in MNRAS

Cosmology Intertwined: A Review of the Particle Physics, Astrophysics, and Cosmology Associated with the Cosmological Tensions and Anomalies

Authors: Elcio Abdalla, Guillermo Franco Abellán, Amin Aboubrahim, Adriano Agnello, Ozgur Akarsu, Yashar Akrami, George Alestas, Daniel Aloni, Luca Amendola, Luis A. Anchordoqui, Richard I. Anderson, Nikki Arendse, Marika Asgari, Mario Ballardini, Vernon Barger, Spyros Basilakos, Ronaldo C. Batista, Elia S. Battistelli, Richard Battye, Micol Benetti, David Benisty, Asher Berlin, Paolo de Bernardis, Emanuele Berti, Bohdan Bidenko , et al. (178 additional authors not shown)

Abstract: In this paper we will list a few important goals that need to be addressed in the next decade, also taking into account the current discordances between the different cosmological probes, such as the disagreement in the value of the Hubble constant $H_0$, the $σ_8$--$S_8$ tension, and other less statistically significant anomalies. While these discordances can still be in part the result of system… ▽ More In this paper we will list a few important goals that need to be addressed in the next decade, also taking into account the current discordances between the different cosmological probes, such as the disagreement in the value of the Hubble constant $H_0$, the $σ_8$--$S_8$ tension, and other less statistically significant anomalies. While these discordances can still be in part the result of systematic errors, their persistence after several years of accurate analysis strongly hints at cracks in the standard cosmological scenario and the necessity for new physics or generalisations beyond the standard model. In this paper, we focus on the $5.0\,σ$ tension between the {\it Planck} CMB estimate of the Hubble constant $H_0$ and the SH0ES collaboration measurements. After showing the $H_0$ evaluations made from different teams using different methods and geometric calibrations, we list a few interesting new physics models that could alleviate this tension and discuss how the next decade's experiments will be crucial. Moreover, we focus on the tension of the {\it Planck} CMB data with weak lensing measurements and redshift surveys, about the value of the matter energy density $Ω_m$, and the amplitude or rate of the growth of structure ($σ_8,fσ_8$). We list a few interesting models proposed for alleviating this tension, and we discuss the importance of trying to fit a full array of data with a single model and not just one parameter at a time. Additionally, we present a wide range of other less discussed anomalies at a statistical significance level lower than the $H_0$--$S_8$ tensions which may also constitute hints towards new physics, and we discuss possible generic theoretical approaches that can collectively explain the non-standard nature of these signals.[Abridged] △ Less

Submitted 24 April, 2022; v1 submitted 11 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021. 224 pages, 27 figures. Accepted for publication in JHEAp

Journal ref: J. High En. Astrophys. 2204, 002 (2022)

AI-driven spatio-temporal engine for finding gravitationally lensed type Ia supernovae

Authors: Doogesh Kodi Ramanah, Nikki Arendse, Radosław Wojtak

Abstract: We present a spatio-temporal AI framework that concurrently exploits both the spatial and time-variable features of gravitationally lensed supernovae in optical images to ultimately aid in future discoveries of such exotic transients in wide-field surveys. Our spatio-temporal engine is designed using recurrent convolutional layers, while drawing from recent advances in variational inference to qua… ▽ More We present a spatio-temporal AI framework that concurrently exploits both the spatial and time-variable features of gravitationally lensed supernovae in optical images to ultimately aid in future discoveries of such exotic transients in wide-field surveys. Our spatio-temporal engine is designed using recurrent convolutional layers, while drawing from recent advances in variational inference to quantify approximate Bayesian uncertainties via a confidence score. Using simulated Young Supernova Experiment (YSE) images of lensed and non-lensed supernovae as a showcase, we find that the use of time-series images adds relevant information from time variability of spatial light distribution of partially blended images of lensed supernova, yielding a substantial gain of around 20 per cent in classification accuracy over single-epoch observations. Preliminary application of our network to mock observations from the Legacy Survey of Space and Time (LSST) results in detections with accuracy reaching around 99 percent. Our innovative deep learning machinery is versatile and can be employed to search for any class of sources which exhibit variability both in flux and spatial distribution of light. △ Less

Submitted 24 March, 2022; v1 submitted 26 July, 2021; originally announced July 2021.

Comments: 14 pages, 13 figures, 2 tables. Accepted for publication in MNRAS

The Young Supernova Experiment: Survey Goals, Overview, and Operations

Authors: D. O. Jones, R. J. Foley, G. Narayan, J. Hjorth, M. E. Huber, P. D. Aleo, K. D. Alexander, C. R. Angus, K. Auchettl, V. F. Baldassare, S. H. Bruun, K. C. Chambers, D. Chatterjee, D. L. Coppejans, D. A. Coulter, L. DeMarchi, G. Dimitriadis, M. R. Drout, A. Engel, K. D. French, A. Gagliano, C. Gall, T. Hung, L. Izzo, W. V. Jacobson-Galán , et al. (46 additional authors not shown)

Abstract: Time domain science has undergone a revolution over the past decade, with tens of thousands of new supernovae (SNe) discovered each year. However, several observational domains, including SNe within days or hours of explosion and faint, red transients, are just beginning to be explored. Here, we present the Young Supernova Experiment (YSE), a novel optical time-domain survey on the Pan-STARRS tele… ▽ More Time domain science has undergone a revolution over the past decade, with tens of thousands of new supernovae (SNe) discovered each year. However, several observational domains, including SNe within days or hours of explosion and faint, red transients, are just beginning to be explored. Here, we present the Young Supernova Experiment (YSE), a novel optical time-domain survey on the Pan-STARRS telescopes. Our survey is designed to obtain well-sampled $griz$ light curves for thousands of transient events up to $z \approx 0.2$. This large sample of transients with 4-band light curves will lay the foundation for the Vera C. Rubin Observatory and the Nancy Grace Roman Space Telescope, providing a critical training set in similar filters and a well-calibrated low-redshift anchor of cosmologically useful SNe Ia to benefit dark energy science. As the name suggests, YSE complements and extends other ongoing time-domain surveys by discovering fast-rising SNe within a few hours to days of explosion. YSE is the only current four-band time-domain survey and is able to discover transients as faint $\sim$21.5 mag in $gri$ and $\sim$20.5 mag in $z$, depths that allow us to probe the earliest epochs of stellar explosions. YSE is currently observing approximately 750 square degrees of sky every three days and we plan to increase the area to 1500 square degrees in the near future. When operating at full capacity, survey simulations show that YSE will find $\sim$5000 new SNe per year and at least two SNe within three days of explosion per month. To date, YSE has discovered or observed 8.3% of the transient candidates reported to the International Astronomical Union in 2020. We present an overview of YSE, including science goals, survey characteristics and a summary of our transient discoveries to date. △ Less

Submitted 5 January, 2021; v1 submitted 19 October, 2020; originally announced October 2020.

Comments: ApJ, in press; more information at https://yse.ucsc.edu/

Simulation-based inference of dynamical galaxy cluster masses with 3D convolutional neural networks

Authors: Doogesh Kodi Ramanah, Radosław Wojtak, Nikki Arendse

Abstract: We present a simulation-based inference framework using a convolutional neural network to infer dynamical masses of galaxy clusters from their observed 3D projected phase-space distribution, which consists of the projected galaxy positions in the sky and their line-of-sight velocities. By formulating the mass estimation problem within this simulation-based inference framework, we are able to quant… ▽ More We present a simulation-based inference framework using a convolutional neural network to infer dynamical masses of galaxy clusters from their observed 3D projected phase-space distribution, which consists of the projected galaxy positions in the sky and their line-of-sight velocities. By formulating the mass estimation problem within this simulation-based inference framework, we are able to quantify the uncertainties on the inferred masses in a straightforward and robust way. We generate a realistic mock catalogue emulating the Sloan Digital Sky Survey (SDSS) Legacy spectroscopic observations (the main galaxy sample) for redshifts $z \lesssim 0.09$ and explicitly illustrate the challenges posed by interloper (non-member) galaxies for cluster mass estimation from actual observations. Our approach constitutes the first optimal machine learning-based exploitation of the information content of the full 3D projected phase-space distribution, including both the virialized and infall cluster regions, for the inference of dynamical cluster masses. We also present, for the first time, the application of a simulation-based inference machinery to obtain dynamical masses of around $800$ galaxy clusters found in the SDSS Legacy Survey, and show that the resulting mass estimates are consistent with mass measurements from the literature. △ Less

Submitted 17 December, 2020; v1 submitted 7 September, 2020; originally announced September 2020.

Comments: 14 pages, 11 figures. Accepted for publication in MNRAS. Contains non-peer reviewed supplementary material on cluster mass function in appendix

A new measurement of the Hubble constant using Type Ia supernovae calibrated with surface brightness fluctuations

Authors: N. Khetan, L. Izzo, M. Branchesi, R. Wojtak, M. Cantiello, C. Murugeshan, A. Agnello, E. Cappellaro, M. Della Valle, C. Gall, J. Hjorth, S. Benetti, E. Brocato, J. Burke, D. Hiramatsu, D. Andrew Howell, L. Tomasella, S. Valenti

Abstract: We present a new calibration of the peak absolute magnitude of SNe Type Ia based on the Surface Brightness Fluctuations (SBF) method, aimed at measuring the value of the Hubble constant. We build a sample of calibrating anchors consisting of 24 SNe hosted in galaxies having SBF distance measurements. Applying a hierarchical Bayesian approach, we calibrate the SNe luminosity and extend it into the… ▽ More We present a new calibration of the peak absolute magnitude of SNe Type Ia based on the Surface Brightness Fluctuations (SBF) method, aimed at measuring the value of the Hubble constant. We build a sample of calibrating anchors consisting of 24 SNe hosted in galaxies having SBF distance measurements. Applying a hierarchical Bayesian approach, we calibrate the SNe luminosity and extend it into the Hubble flow by using a sample of 96 SNe Ia in the redshift range $0.02 < z < 0.075$, extracted from the Combined Pantheon Sample. We estimate a value of $H_0 = 70.50 \pm 2.37(stat) \pm 3.38(sys)$ $\text{km}\ \text{s}^{-1}\ \text{Mpc}^{-1}$ (i.e. $3.4\% stat, 4.8\% sys$), which is in agreement with the value obtained using the tip of the red giant branch calibration, and consistent within the errors with the value obtained from SNe Type Ia calibrated with Cepheids and the one inferred from the analysis of the cosmic microwave background. We find that the SNe Ia distance moduli calibrated with SBF are on average larger by 0.07 mag than the ones calibrated with Cepheids. Our results point to possible differences among SNe in different types of galaxies, which could originate from different local environments and/or SNe Ia progenitor properties. Sampling different host galaxy type, SBF offers a complementary approach to Cepheids which is important in addressing possible systematics. As the SBF method has the ability to reach larger distances than Cepheids, the impending entry of LSST and JWST into operation will increase the number of SNe Ia hosted in galaxies where SBF distances can be measured, making SBF measurements attractive for improving the calibration of SNe Ia, and in the estimation of $H_0$. △ Less

Submitted 3 March, 2021; v1 submitted 18 August, 2020; originally announced August 2020.

Comments: 22 pages, 14 figures, 9 tables. Comments are very welcome

Journal ref: A&A 647, A72 (2021)

Dynamical mass inference of galaxy clusters with neural flows

Authors: Doogesh Kodi Ramanah, Radosław Wojtak, Zoe Ansari, Christa Gall, Jens Hjorth

Abstract: We present an algorithm for inferring the dynamical mass of galaxy clusters directly from their respective phase-space distributions, i.e. the observed line-of-sight velocities and projected distances of galaxies from the cluster centre. Our method employs normalizing flows, a deep neural network capable of learning arbitrary high-dimensional probability distributions, and inherently accounts, to… ▽ More We present an algorithm for inferring the dynamical mass of galaxy clusters directly from their respective phase-space distributions, i.e. the observed line-of-sight velocities and projected distances of galaxies from the cluster centre. Our method employs normalizing flows, a deep neural network capable of learning arbitrary high-dimensional probability distributions, and inherently accounts, to an adequate extent, for the presence of interloper galaxies which are not bounded to a given cluster, the primary contaminant of dynamical mass measurements. We validate and showcase the performance of our neural flow approach to robustly infer the dynamical mass of clusters from a realistic mock cluster catalogue. A key aspect of our novel algorithm is that it yields the probability density function of the mass of a particular cluster, thereby providing a principled way of quantifying uncertainties, in contrast to conventional machine learning approaches. The neural network mass predictions, when applied to a contaminated catalogue with interlopers, have a mean overall logarithmic residual scatter of 0.028 dex, with a log-normal scatter of 0.126 dex, which goes down to 0.089 dex for clusters in the intermediate to high mass range. This is an improvement by nearly a factor of four relative to the classical cluster mass scaling relation with the velocity dispersion, and outperforms recently proposed machine learning approaches. We also apply our neural flow mass estimator to a compilation of galaxy observations of some well-studied clusters with robust dynamical mass estimates, further substantiating the efficacy of our algorithm. △ Less

Submitted 8 September, 2020; v1 submitted 12 March, 2020; originally announced March 2020.

Comments: 14 pages, 9 figures, 1 table. Improved approach, saliency maps added

The lowest of the low: discovery of SN 2019gsc and the nature of faint Iax supernovae

Authors: Shubham Srivastav, Stephen J. Smartt, Giorgos Leloudas, Mark E. Huber, Ken Chambers, Daniele B. Malesani, Jens Hjorth, James H. Gillanders, A. Schultz, Stuart A. Sim, Katie Auchettl, Johan P. U. Fynbo, Christa Gall, Owen R. McBrien, Armin Rest, Ken W. Smith, Radoslaw Wojtak, David R. Young

Abstract: We present the discovery and optical follow-up of the faintest supernova-like transient known. The event (SN 2019gsc) was discovered in a star-forming host at 53\,Mpc by ATLAS. A detailed multi-colour light curve was gathered with Pan-STARRS1 and follow-up spectroscopy was obtained with the NOT and Gemini-North. The spectra near maximum light show narrow features at low velocities of 3000 to 4000… ▽ More We present the discovery and optical follow-up of the faintest supernova-like transient known. The event (SN 2019gsc) was discovered in a star-forming host at 53\,Mpc by ATLAS. A detailed multi-colour light curve was gathered with Pan-STARRS1 and follow-up spectroscopy was obtained with the NOT and Gemini-North. The spectra near maximum light show narrow features at low velocities of 3000 to 4000 km s$^{-1}$, similar to the extremely low luminosity SNe 2010ae and 2008ha, and the light curve displays a similar fast decline (\dmr $0.91 \pm 0.10$ mag). SNe 2010ae and 2008ha have been classified as type Iax supernovae, and together the three either make up a distinct physical class of their own or are at the extreme low luminosity end of this diverse supernova population. The bolometric light curve is consistent with a low kinetic energy of explosion ($E_{\rm k} \sim 10^{49}$ erg s$^{-1}$), a modest ejected mass ($M_{\rm ej} \sim 0.2$ \msol) and radioactive powering by $^{56}$Ni ($M_{\rm Ni} \sim 2 \times 10^{-3}$ \msol). The spectra are quite well reproduced with radiative transfer models (TARDIS) and a composition dominated by carbon, oxygen, magnesium, silicon and sulphur. Remarkably, all three of these extreme Iax events are in similar low-metallicity star-forming environments. The combination of the observational constraints for all three may be best explained by deflagrations of near $M_{\rm Ch}$ hybrid carbon-oxygen-neon white dwarfs which have short evolutionary pathways to formation. △ Less

Submitted 18 February, 2020; v1 submitted 27 January, 2020; originally announced January 2020.

Comments: 12 pages, 5 figures, accepted to ApJL, minor changes to submitted version

Cosmic dissonance: new physics or systematics behind a short sound horizon?

Authors: Nikki Arendse, Radosław J. Wojtak, Adriano Agnello, Geoff C. -F. Chen, Christopher D. Fassnacht, Dominique Sluse, Stefan Hilbert, Martin Millon, Vivien Bonvin, Kenneth C. Wong, Frédéric Courbin, Sherry H. Suyu, Simon Birrer, Tommaso Treu, Leon V. E. Koopmans

Abstract: Persistent tension between low-redshift observations and the Cosmic Microwave Background radiation (CMB), in terms of two fundamental distance scales set by the sound horizon $r_d$ and the Hubble constant $H_0$, suggests new physics beyond the Standard Model or residual systematics. We examine recently updated distance calibrations from Cepheids, gravitational lensing time-delay observations, and… ▽ More Persistent tension between low-redshift observations and the Cosmic Microwave Background radiation (CMB), in terms of two fundamental distance scales set by the sound horizon $r_d$ and the Hubble constant $H_0$, suggests new physics beyond the Standard Model or residual systematics. We examine recently updated distance calibrations from Cepheids, gravitational lensing time-delay observations, and the Tip of the Red Giant Branch. Calibrating the Baryon Acoustic Oscillations (BAO) and Type Ia supernovae with combinations of the distance indicators, we obtain a joint and self-consistent measurement of $H_0$ and $r_d$ at low redshift, independent of cosmological models and CMB inference. In an attempt to alleviate the tension between late-time and CMB-based measurements, we consider four extensions of the standard $Λ$CDM model. The sound horizon from our different measurements is $r_d=(137\pm3^{stat.}\pm2^{syst.})$~Mpc. Depending on the adopted distance indicators, the $combined$ tension in $H_0$ and $r_d$ ranges between 2.3 and 5.1 $σ$. We find that modifications of $Λ$CDM that change the physics after recombination fail to solve the problem, for the reason that they only resolve the tension in $H_0$, while the tension in $r_d$ remains unchanged. Pre-recombination extensions (with early dark energy or the effective number of neutrinos $\rm{N}_{\rm{eff}}=3.24 \pm 0.16$) are allowed by the data, unless the calibration from Cepheids is included. Results from time-delay lenses are consistent with those from distance-ladder calibrations and point to a discrepancy between absolute distance scales measured from the CMB (assuming the standard cosmological model) and late-time observations. New proposals to resolve this tension should be examined with respect to reconciling not only the Hubble constant but also the sound horizon derived from the CMB and other cosmological probes. △ Less

Submitted 13 May, 2020; v1 submitted 17 September, 2019; originally announced September 2019.

Comments: Accepted for publication in A&A; 14 pages, 8 figures, 4 tables

Journal ref: A&A 639, A57 (2020)

Einstein-Vlasov Calculations of Structure Formation

Authors: William E. East, Radosław Wojtak, Frans Pretorius

Abstract: We study the dynamics of small inhomogeneities in an expanding universe collapsing to form bound structures using full solutions of the Einstein-Vlasov (N-body) equations. We compare these to standard Newtonian N-body solutions using quantities defined with respect to fiducial observers in order to bound relativistic effects. We focus on simplified initial conditions containing a limited range of… ▽ More We study the dynamics of small inhomogeneities in an expanding universe collapsing to form bound structures using full solutions of the Einstein-Vlasov (N-body) equations. We compare these to standard Newtonian N-body solutions using quantities defined with respect to fiducial observers in order to bound relativistic effects. We focus on simplified initial conditions containing a limited range of length scales, but vary the inhomogeneities from small magnitude, where the Newtonian and general-relativistic calculations agree quite well, to large magnitude, where the background metric receives an order one correction. For large inhomogeneities, we find that the collapse of overdensities tends to happen faster in Newtonian calculations relative to fully general-relativistic ones. Even in this extreme regime, the differences in the spacetime evolution outside the regions of large gravitational potential and velocity are small. For standard cosmological values, we corroborate the robustness of Newtonian N-body simulations to model large scale perturbations and the related cosmic variance in the local expansion rate. △ Less

Submitted 22 November, 2019; v1 submitted 15 August, 2019; originally announced August 2019.

Comments: 11 pages, 8 figures; revised to match PRD version

Journal ref: Phys. Rev. D 100, 103533 (2019)

The Hubble-Lemaître constant and sound horizon from low-redshift probes

Authors: Radosław Wojtak, Adriano Agnello

Abstract: We revisit the claimed tension, or lack thereof, of measured values of the Hubble-Lemaître parameter $H_0$ from Cosmic Microwave Background (CMB) data and low-redshift indicators. Baryon Acoustic Oscillations (BAO) rely on the scale of the sound horizon at recombination $r_s$ to convert angular measurements into angular-diameter distances, so fixing $r_s$ from CMB measurements already constrains… ▽ More We revisit the claimed tension, or lack thereof, of measured values of the Hubble-Lemaître parameter $H_0$ from Cosmic Microwave Background (CMB) data and low-redshift indicators. Baryon Acoustic Oscillations (BAO) rely on the scale of the sound horizon at recombination $r_s$ to convert angular measurements into angular-diameter distances, so fixing $r_s$ from CMB measurements already constrains $H_0.$ If departures from concordance cosmology are to be constrained, truly independent measurements of $H_0$ are needed. We use the angular-diameter distances to three time-delay lenses from the H0LiCOW collaboration to calibrate the distance ladder, combine them with relative distances from Supernovae Ia and BAO, leaving $r_s$ completely free, and provide the inferred coefficients ($q_{0},j_{0},s_{0}$) in the polynomial expansion of H(z). We obtain $H_{0}r_{s}=(9895\pm161)$km/s and $H_0=(72\pm7)$km/s/Mpc. Combined with $H_0$ from H0LiCOW, then $r_s=(137\pm4.5)$Mpc is consistent with previous work and systematically lower than the CMB-inferred value. Our results are independent of the adopted cosmology, and removing Supernovae with z<0.1 has a negligible effect. △ Less

Submitted 6 August, 2019; originally announced August 2019.

Comments: 6 pages, 2 figures, 2 tables; published in MNRAS

Journal ref: MNRAS, 486, 5046 (2019)

Low-redshift measurement of the sound horizon through gravitational time-delays

Authors: Nikki Arendse, Adriano Agnello, Radosław Wojtak

Abstract: The matter sound horizon can be inferred from the cosmic microwave background within the Standard Model. Independent direct measurements of the sound horizon are then a probe of possible deviations from the Standard Model. We aim at measuring the sound horizon $r_s$ from low-redshift indicators, which are completely independent of CMB inference. We used the measured product $H(z)r_s$ from baryon a… ▽ More The matter sound horizon can be inferred from the cosmic microwave background within the Standard Model. Independent direct measurements of the sound horizon are then a probe of possible deviations from the Standard Model. We aim at measuring the sound horizon $r_s$ from low-redshift indicators, which are completely independent of CMB inference. We used the measured product $H(z)r_s$ from baryon acoustic oscillations (BAO) together with supernovae~\textsc{I}a to constrain $H(z)/H_{0}$ and time-delay lenses analysed by the H0LiCOW collaboration to anchor cosmological distances ($\propto H_{0}^{-1}$). {Additionally, we investigated the influence of adding a sample of quasars with higher redshift with standardisable UV-Xray luminosity distances. We adopted polynomial expansions in $H(z)$ or in comoving distances} so that our inference was completely independent of any cosmological model on which the expansion history might be based. Our measurements are independent of Cepheids and systematics from peculiar motions {to within percent-level accuracy.} The inferred sound horizon $r_s$ varies between $(133 \pm 8)$~Mpc and $(138 \pm 5)$~Mpc across different models. The discrepancy with CMB measurements is robust against model choice. Statistical uncertainties are comparable to systematics. The combination of time-delay lenses, supernovae, and BAO yields a distance ladder that is independent of cosmology (and of Cepheid calibration) and a measurement of $r_s $ that is independent of the CMB. These cosmographic measurements are then a competitive test of the Standard Model, regardless of the hypotheses on which the cosmology is based. △ Less

Submitted 19 February, 2020; v1 submitted 28 May, 2019; originally announced May 2019.

Comments: 2019A&A...632A..91A (subm. 28/05/2019), 6 pages, 3 figures

Journal ref: Astronomy & Astrophysics, 2019, Volume 632, id.A91, 6 pp. Astronomy & Astrophysics, Volume 632, id.A91, 6 pp

A new measurement of the Hubble constant and matter content of the Universe using extragalactic background light $γ$-ray attenuation

Authors: Alberto Domínguez, Radoslaw Wojtak, Justin Finke, Marco Ajello, Kari Helgason, Francisco Prada, Abhishek Desai, Vaidehi Paliya, Lea Marcotulli, Dieter Hartmann

Abstract: The Hubble constant $H_{0}$ and matter density $Ω_{m}$ of the Universe are measured using the latest $γ$-ray attenuation results from Fermi-LAT and Cherenkov telescopes. This methodology is based upon the fact that the extragalactic background light supplies opacity for very high energy photons via photon-photon interaction. The amount of $γ$-ray attenuation along the line of sight depends on the… ▽ More The Hubble constant $H_{0}$ and matter density $Ω_{m}$ of the Universe are measured using the latest $γ$-ray attenuation results from Fermi-LAT and Cherenkov telescopes. This methodology is based upon the fact that the extragalactic background light supplies opacity for very high energy photons via photon-photon interaction. The amount of $γ$-ray attenuation along the line of sight depends on the expansion rate and matter content of the Universe. This novel strategy results in a value of $H_{0}=67.4_{-6.2}^{+6.0}$~km~s$^{-1}$~Mpc$^{-1}$ and $Ω_{m}=0.14_{-0.07}^{+0.06}$. These estimates are independent and complementary to those based on the distance ladder, cosmic microwave background (CMB), clustering with weak lensing, and strong lensing data. We also produce a joint likelihood analysis of our results from $γ$ rays and these from more mature methodologies, excluding the CMB, yielding a combined value of $H_{0}=66.6\pm 1.6$~km~s$^{-1}$~Mpc$^{-1}$ and $Ω_{m}=0.29\pm 0.02$. △ Less

Submitted 8 November, 2019; v1 submitted 28 March, 2019; originally announced March 2019.

Comments: 9 pages, 6 figures, 1 table. Accepted by ApJ

Magnified or multiply imaged? - Search strategies for gravitationally lensed supernovae in wide-field surveys

Authors: Radosław Wojtak, Jens Hjorth, Christa Gall

Abstract: Strongly lensed supernovae can be detected as multiply imaged or highly magnified transients. In order to compare the performances of these two observational strategies, we calculate expected discovery rates as a function of survey depth in five grizy filters and for different classes of supernovae (Ia, IIP, IIL, Ibc and IIn). We find that detections via magnification is the only effective strateg… ▽ More Strongly lensed supernovae can be detected as multiply imaged or highly magnified transients. In order to compare the performances of these two observational strategies, we calculate expected discovery rates as a function of survey depth in five grizy filters and for different classes of supernovae (Ia, IIP, IIL, Ibc and IIn). We find that detections via magnification is the only effective strategy for relatively shallow pre-LSST surveys. For survey depths about the LSST capacity, both strategies yield comparable numbers of lensed supernovae. Supernova samples from the two methods are to a large extent independent and combining them increases detection rates by about 50 per cent. While the number of lensed supernovae detectable via magnification saturates at the limiting magnitudes of LSST, detection rates of multiply imaged supernova still go up drastically at increasing survey depth. Comparing potential discovery spaces, we find that lensed supernovae found via image multiplicity exhibit longer time delays and larger image separations making them more suitable for cosmological constraints than their counterparts found via magnification. We find that the ZTF will find about 2 type Ia and 4 core-collapse lensed supernovae per year at a limiting magnitude of 20.6 in the r band. Applying a hybrid method which combines searching for highly magnified or multiply imaged transients, we find that LSST will detect 89 type Ia and 254 core-collapse lensed supernovae per year. In all cases, lensed core-collapsed supernovae will be dominated by type IIn supernovae contributing to 80 per cent of the total counts, although this prediction relies quite strongly on the adopted spectral templates for this class of supernovae. Revisiting the case of the lensed supernova iPTF16geu, we find that it is consistent within the 2σcontours of predicted redshifts and magnifications for the iPTF survey. △ Less

Submitted 3 June, 2019; v1 submitted 18 March, 2019; originally announced March 2019.

Comments: 16 pages, 11 figures, 2 tables; accepted for publication in MNRAS

Journal ref: MNRAS, 487, 3342 (2019)

Galaxy Cluster Mass Reconstruction Project - IV. Understanding the effects of imperfect membership on cluster mass estimation

Authors: R. Wojtak, L. Old, G. A. Mamon, F. R. Pearce, R. de Carvalho, C. Sifón, M. E. Gray, R. A. Skibba, D. Croton, S. Bamford, D. Gifford, A. von der Linden, J. C. Muñoz-Cuartas, V. Müller, R. J. Pearson, E. Rozo, E. Rykoff, A. Saro, T. Sepp, E. Tempel

Abstract: The primary difficulty in measuring dynamical masses of galaxy clusters from galaxy data lies in the separation between true cluster members from interloping galaxies along the line of sight. We study the impact of membership contamination and incompleteness on cluster mass estimates obtained with 25 commonly used techniques applied to nearly 1000 mock clusters. We show that all methods overestima… ▽ More The primary difficulty in measuring dynamical masses of galaxy clusters from galaxy data lies in the separation between true cluster members from interloping galaxies along the line of sight. We study the impact of membership contamination and incompleteness on cluster mass estimates obtained with 25 commonly used techniques applied to nearly 1000 mock clusters. We show that all methods overestimate or underestimate cluster masses when applied to contaminated or incomplete galaxy samples respectively. This appears to be the main source of the intrinsic scatter in the mass scaling relation. Applying corrections based on a prior knowledge of contamination and incompleteness can reduce the scatter to the level of shot noise expected for poorly sampled clusters. We establish an empirical model quantifying the effect of imperfect membership on cluster mass estimation and discuss its universal and method-dependent features. We find that both imperfect membership and the response of the mass estimators depend on cluster mass, effectively causing a flattening of the estimated - true mass relation. Imperfect membership thus alters cluster counts determined from spectroscopic surveys, hence the cosmological parameters that depend on such counts. △ Less

Submitted 16 August, 2018; v1 submitted 8 June, 2018; originally announced June 2018.

Comments: 18 pages, 11 figures, 3tables; accepted for publication in MNRAS

Journal ref: MNRAS, 481, 324 (2018)

Comparing Fully General Relativistic and Newtonian Calculations of Structure Formation

Authors: William E. East, Radosław Wojtak, Tom Abel

Abstract: In the standard approach to studying cosmological structure formation, the overall expansion of the Universe is assumed to be homogeneous, with the gravitational effect of inhomogeneities encoded entirely in a Newtonian potential. A topic of ongoing debate is to what degree this fully captures the dynamics dictated by general relativity, especially in the era of precision cosmology. To quantitativ… ▽ More In the standard approach to studying cosmological structure formation, the overall expansion of the Universe is assumed to be homogeneous, with the gravitational effect of inhomogeneities encoded entirely in a Newtonian potential. A topic of ongoing debate is to what degree this fully captures the dynamics dictated by general relativity, especially in the era of precision cosmology. To quantitatively assess this, we directly compare standard N-body Newtonian calculations to full numerical solutions of the Einstein equations, for cold matter with various magnitude initial inhomogeneities on scales comparable to the Hubble horizon. We analyze the differences in the evolution of density, luminosity distance, and other quantities defined with respect to fiducial observers. This is carried out by reconstructing the effective spacetime and matter fields dictated by the Newtonian quantities, and by taking care to distinguish effects of numerical resolution. We find that the fully general relativistic and Newtonian calculations show excellent agreement, even well into the nonlinear regime. They only notably differ in regions where the weak gravity assumption breaks down, which arise when considering extreme cases with perturbations exceeding standard values. △ Less

Submitted 9 February, 2018; v1 submitted 17 November, 2017; originally announced November 2017.

Comments: 17 pages, 14 figures; revised to match PRD version

Journal ref: Phys. Rev. D 97, 043509 (2018)

The Distance to NGC 4993: The Host Galaxy of the Gravitational-wave Event GW170817

Authors: Jens Hjorth, Andrew J. Levan, Nial R. Tanvir, Joe D. Lyman, Radosław Wojtak, Sophie L. Schrøder, Ilya Mandel, Christa Gall, Sofie H. Bruun

Abstract: The historic detection of gravitational waves from a binary neutron star merger (GW170817) and its electromagnetic counterpart led to the first accurate (sub-arcsecond) localization of a gravitational-wave event. The transient was found to be $\sim$10" from the nucleus of the S0 galaxy NGC 4993. We report here the luminosity distance to this galaxy using two independent methods. (1) Based on our M… ▽ More The historic detection of gravitational waves from a binary neutron star merger (GW170817) and its electromagnetic counterpart led to the first accurate (sub-arcsecond) localization of a gravitational-wave event. The transient was found to be $\sim$10" from the nucleus of the S0 galaxy NGC 4993. We report here the luminosity distance to this galaxy using two independent methods. (1) Based on our MUSE/VLT measurement of the heliocentric redshift ($z_{\rm helio}=0.009783\pm0.000023$) we infer the systemic recession velocity of the NGC 4993 group of galaxies in the cosmic microwave background (CMB) frame to be $v_{\rm CMB}=3231 \pm 53$ km s$^{-1}$. Using constrained cosmological simulations we estimate the line-of-sight peculiar velocity to be $v_{\rm pec}=307 \pm 230$ km s$^{-1}$, resulting in a cosmic velocity of $v_{\rm cosmic}=2924 \pm 236$ km s$^{-1}$ ($z_{\rm cosmic}=0.00980\pm 0.00079$) and a distance of $D_z=40.4\pm 3.4$ Mpc assuming a local Hubble constant of $H_0=73.24\pm 1.74$ km s$^{-1}$ Mpc$^{-1}$. (2) Using Hubble Space Telescope measurements of the effective radius (15.5" $\pm$ 1.5") and contained intensity and MUSE/VLT measurements of the velocity dispersion, we place NGC 4993 on the Fundamental Plane (FP) of E and S0 galaxies. Comparing to a frame of 10 clusters containing 226 galaxies, this yields a distance estimate of $D_{\rm FP}=44.0\pm 7.5$ Mpc. The combined redshift and FP distance is $D_{\rm NGC 4993}= 41.0\pm 3.1$ Mpc. This 'electromagnetic' distance estimate is consistent with the independent measurement of the distance to GW170817 as obtained from the gravitational-wave signal ($D_{\rm GW}= 43.8^{+2.9}_{-6.9}$ Mpc) and confirms that GW170817 occurred in NGC 4993. △ Less

Submitted 16 October, 2017; originally announced October 2017.

Comments: 9 pages, 5 figures

Journal ref: The Astrophysical Journal Letters 848:L31 (2017)

Galaxy Cluster Mass Reconstruction Project: III. The impact of dynamical substructure on cluster mass estimates

Authors: L. Old, R. Wojtak, F. R. Pearce, M. E. Gray, G. A. Mamon, C. Sifón, E. Tempel, A. Biviano, H. K. C. Yee, R. de Carvalho, V. Müller, T. Sepp, R. A. Skibba, D. Croton, S. P. Bamford C. Power, A. von der Linden, A. Saro

Abstract: 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 pr… ▽ More 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. △ Less

Submitted 28 September, 2017; originally announced September 2017.

Comments: 14 pages, 8 figures, 3 tables. Submitted to MNRAS

Redshift remapping and cosmic acceleration in dark-matter-dominated cosmological models

Authors: Radosław Wojtak, Francisco Prada

Abstract: The standard relation between the cosmological redshift and cosmic scale factor underlies cosmological inference from virtually all kinds of cosmological observations, leading to the emergence of the LambdaCDM cosmological model. This relation is not a fundamental theory and thus observational determination of this function (redshift remapping) should be regarded as an insightful alternative to ho… ▽ More The standard relation between the cosmological redshift and cosmic scale factor underlies cosmological inference from virtually all kinds of cosmological observations, leading to the emergence of the LambdaCDM cosmological model. This relation is not a fundamental theory and thus observational determination of this function (redshift remapping) should be regarded as an insightful alternative to holding its standard form in analyses of cosmological data. Here we present non-parametric reconstructions of redshift remapping in dark-matter-dominated models and constraints on cosmological parameters from a joint analysis of all primary cosmological probes including the local measurement of the Hubble constant, Type Ia supernovae, baryonic acoustic oscillations (BAO), Planck observations of the cosmic microwave background (CMB) radiation (temperature power spectrum) and cosmic chronometers. The reconstructed redshift remapping points to an additional boost of redshift operating in late epoch of cosmic evolution, but affecting both low-redshift observations and the CMB. The model predicts a significant difference between the actual Hubble constant, h=0.48+/-0.02, and its local determination, h_obs=0.73+/-0.02. The ratio of these two values coincides closely with the maximum expansion rate inside voids formed in the corresponding open cosmological model with Omega_m=0.87+/-0.03, whereas the actual value of the Hubble constant implies the age of the Universe that is compatible with the Planck LambdaCDM cosmology. The new dark-matter-dominated model with redshift remapping provides excellent fits to all data and eliminates recently reported tensions between the Planck LambdaCDM cosmology, the local determination of the Hubble constant and the BAO measurements from the Ly-alpha forest of high-redshift quasars. △ Less

Submitted 26 June, 2017; v1 submitted 12 October, 2016; originally announced October 2016.

Comments: 21 pages, 11 figures, 4 tables; accepted for publication in MNRAS

Journal ref: MNRAS, 470, 4493 (2017)

Voids in cosmological simulations over cosmic time

Authors: Radosław Wojtak, Devon Powell, Tom Abel

Abstract: We study evolution of voids in cosmological simulations using a new method for tracing voids over cosmic time. The method is based on tracking watershed basins (contiguous regions around density minima) of well developed voids at low redshift, on a regular grid of density field. It enables us to construct a robust and continuous mapping between voids at different redshifts, from initial conditions… ▽ More We study evolution of voids in cosmological simulations using a new method for tracing voids over cosmic time. The method is based on tracking watershed basins (contiguous regions around density minima) of well developed voids at low redshift, on a regular grid of density field. It enables us to construct a robust and continuous mapping between voids at different redshifts, from initial conditions to the present time. We discuss how the new approach eliminates strong spurious effects of numerical origin when voids evolution is traced by matching voids between successive snapshots (by analogy to halo merger trees). We apply the new method to a cosmological simulation of a standard LambdaCDM cosmological model and study evolution of basic properties of typical voids (with effective radii between 6Mpc/h and 20Mpc/h at redshift z=0) such as volumes, shapes, matter density distributions and relative alignments. The final voids at low redshifts appear to retain a significant part of the configuration acquired in initial conditions. Shapes of voids evolve in a collective way which barely modifies the overall distribution of the axial ratios. The evolution appears to have a weak impact on mutual alignments of voids implying that the present state is in large part set up by the primordial density field. We present evolution of dark matter density profiles computed on iso-density surfaces which comply with the actual shapes of voids. Unlike spherical density profiles, this approach enables us to demonstrate development of theoretically predicted bucket-like shape of the final density profiles indicating a wide flat core and a sharp transition to high-density void walls. △ Less

Submitted 30 March, 2016; v1 submitted 26 February, 2016; originally announced February 2016.

Comments: 13 pages, 13 figures; accepted for publication in MNRAS

Journal ref: MNRAS, 458, 4431 (2016)

Testing the mapping between redshift and cosmic scale factor

Authors: Radosław Wojtak, Francisco Prada

Abstract: The canonical redshift-scale factor relation, 1/a=1+z, is a key element in the standard LambdaCDM model of the big bang cosmology. Despite its fundamental role, this relation has not yet undergone any observational tests since Lemaitre and Hubble established the expansion of the Universe. It is strictly based on the assumption of the Friedmann-Lemaitre-Robertson-Walker metric describing a locally… ▽ More The canonical redshift-scale factor relation, 1/a=1+z, is a key element in the standard LambdaCDM model of the big bang cosmology. Despite its fundamental role, this relation has not yet undergone any observational tests since Lemaitre and Hubble established the expansion of the Universe. It is strictly based on the assumption of the Friedmann-Lemaitre-Robertson-Walker metric describing a locally homogeneous and isotropic universe and that photons move on null geodesics of the metric. Thus any violation of this assumption, within general relativity or modified gravity, can yield a different mapping between the model redshift z=1/a-1 and the actually observed redshift z_obs, i.e. z_obs neq z. Here we perform a simple test of consistency for the standard redshift-scale factor relation by determining simultaneous observational constraints on the concordance LambdaCDM cosmological parameters and a generalized redshift mapping z=f(z_obs). Using current baryon acoustic oscillations (BAO) and Type Ia supernova (SN) data we demonstrate that the generalized redshift mapping is strongly degenerated with dark energy. Marginalization over a class of monotonic functions f(z_obs) changes substantially degeneracy between matter and dark energy density: the density parameters become anti correlated with nearly vertical axis of degeneracy. Furthermore, we show that current SN and BAO data, analysed in a framework with the generalized redshift mapping, do not constrain dark energy unless the BAO data include the measurements from the Ly-alpha forest of high-redshift quasars. △ Less

Submitted 17 March, 2016; v1 submitted 6 February, 2016; originally announced February 2016.

Comments: 11 pages, 5 figures, 3 tables; accepted for publication in MNRAS

Journal ref: MNRAS, 458, 3331 (2016)

Non-universality of dark-matter halos: cusps, cores, and the central potential

Authors: Jens Hjorth, Liliya L. R. Williams, Radosław Wojtak, Michael McLaughlin

Abstract: Dark-matter halos grown in cosmological simulations appear to have central NFW-like density cusps with mean values of $d\logρ/d\log r \approx -1$, and some dispersion, which is generally parametrized by the varying index $α$ in the Einasto density profile fitting function. Non-universality in profile shapes is also seen in observed galaxy clusters and possibly dwarf galaxies. Here we show that non… ▽ More Dark-matter halos grown in cosmological simulations appear to have central NFW-like density cusps with mean values of $d\logρ/d\log r \approx -1$, and some dispersion, which is generally parametrized by the varying index $α$ in the Einasto density profile fitting function. Non-universality in profile shapes is also seen in observed galaxy clusters and possibly dwarf galaxies. Here we show that non-universality, at any given mass scale, is an intrinsic property of DARKexp, a theoretically derived model for collisionless self-gravitating systems. We demonstrate that DARKexp - which has only one shape parameter, $φ_0$ - fits the dispersion in profile shapes of massive simulated halos as well as observed clusters very well. DARKexp also allows for cored dark-matter profiles, such as those found for dwarf spheroidal galaxies. We provide approximate analytical relations between DARKexp $φ_0$, Einasto $α$, or the central logarithmic slope in the Dehnen-Tremaine analytical $γ$-models. The range in halo parameters reflects a substantial variation in the binding energies per unit mass of dark-matter halos. △ Less

Submitted 10 August, 2015; originally announced August 2015.

Comments: ApJ, in press, 10 pages, 7 figures

Journal ref: ApJ 811, 2 (2015)

Local gravitational redshifts can bias cosmological measurements

Authors: Radoslaw Wojtak, Tamara M. Davis, Jophiel Wiis

Abstract: Measurements of cosmological parameters via the distance-redshift relation usually rely on models that assume a homogenous universe. It is commonly presumed that the large-scale structure evident in our Universe has a negligible impact on the measurement if distances probed in observations are sufficiently large (compared to the scale of inhomogeneities) and are averaged over different directions… ▽ More Measurements of cosmological parameters via the distance-redshift relation usually rely on models that assume a homogenous universe. It is commonly presumed that the large-scale structure evident in our Universe has a negligible impact on the measurement if distances probed in observations are sufficiently large (compared to the scale of inhomogeneities) and are averaged over different directions on the sky. This presumption does not hold when considering the effect of the gravitational redshift caused by our local gravitational potential, which alters light coming from all distances and directions in the same way. Despite its small magnitude, this local gravitational redshift gives rise to noticeable effects in cosmological inference using SN Ia data. Assuming conservative prior knowledge of the local potential given by sampling a range of gravitational potentials at locations of Milky-Way-like galaxies identified in cosmological simulations, we show that ignoring the gravitational redshift effect in a standard data analysis leads to an additional systematic error of ~1 per cent in the determination of density parameters and the dark energy equation of state. We conclude that our local gravitational field affects our cosmological inference at a level that is important in future observations aiming to achieve percent-level accuracy. △ Less

Submitted 23 June, 2015; v1 submitted 2 April, 2015; originally announced April 2015.

Comments: 18 pages, 8 figures, accepted for publication in JCAP

Journal ref: JCAP, 07, 025 (2015)

Galaxy Cluster Mass Reconstruction Project: II. Quantifying scatter and bias using contrasting mock catalogues

Authors: L. Old, R. Wojtak, G. A. Mamon, R. A. Skibba, F. R. Pearce, D. Croton, S. Bamford, P. Behroozi, R. de Carvalho, J. C. Muñoz-Cuartas, D. Gifford, M. E. Gray, A. von der Linden, M. R. Merrifield, S. I. Muldrew, V. Müller, R. J. Pearson, T. J. Ponman, E. Rozo, E. Rykoff, A. Saro, T. Sepp, C. Sifón, E. Tempel

Abstract: This article is the second in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our aim is to quantify the scatter, systematic bias and completeness of cluster masses derived from a diverse set of 25 galaxy-based methods using two contrasting mock galaxy catalogues based on… ▽ More This article is the second in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our aim is to quantify the scatter, systematic bias and completeness of cluster masses derived from a diverse set of 25 galaxy-based methods using two contrasting mock galaxy catalogues based on a sophisticated halo occupation model and a semi-analytic model. Analysing 968 clusters, we find a wide range in the RMS errors in log M200c delivered by the different methods (0.18 to 1.08 dex, i.e., a factor of ~1.5 to 12), with abundance matching and richness methods providing the best results, irrespective of the input model assumptions. In addition, certain methods produce a significant number of catastrophic cases where the mass is under- or over-estimated by a factor greater than 10. Given the steeply falling high-mass end of the cluster mass function, we recommend that richness or abundance matching-based methods are used in conjunction with these methods as a sanity check for studies selecting high mass clusters. We see a stronger correlation of the recovered to input number of galaxies for both catalogues in comparison with the group/cluster mass, however, this does not guarantee that the correct member galaxies are being selected. We do not observe significantly higher scatter for either mock galaxy catalogues. Our results have implications for cosmological analyses that utilise the masses, richnesses, or abundances of clusters, which have different uncertainties when different methods are used. △ Less

Submitted 25 February, 2015; originally announced February 2015.

Comments: 25 pages, 19 figures, 7 tables. Accepted for publication in MNRAS

Dark matter inner slope and concentration in galaxies: from the Fornax dwarf to M87

Authors: G. A. Mamon, J. Chevalier, A. J. Romanowsky, R. Wojtak

Abstract: We apply two new state-of-the-art methods that model the distribution of observed tracers in projected phase space to lift the mass / velocity anisotropy (VA) degeneracy and deduce constraints on the mass profiles of galaxies, as well as their VA. We first show how a distribution function based method applied to the satellite kinematics of otherwise isolated SDSS galaxies shows convincing observat… ▽ More We apply two new state-of-the-art methods that model the distribution of observed tracers in projected phase space to lift the mass / velocity anisotropy (VA) degeneracy and deduce constraints on the mass profiles of galaxies, as well as their VA. We first show how a distribution function based method applied to the satellite kinematics of otherwise isolated SDSS galaxies shows convincing observational evidence of age matching: red galaxies have more concentrated dark matter (DM) halos than blue galaxies of the same stellar or halo mass. Then, applying the MAMPOSSt technique to M87 (traced by its red and blue globular clusters) we find that very cuspy DM is favored, unless we release priors on DM concentration or stellar mass (leading to unconstrained slope). For the Fornax dwarf spheroidal (traced by its metal-rich and metal-poor stars), the inner DM slope is unconstrained, with weak evidence for a core if the stellar mass is fixed. This highlights how priors are crucial for DM modeling. Finally, we find that blue GCs around M87 and metal-rich stars in Fornax have tangential outer VA. △ Less

Submitted 27 October, 2014; originally announced October 2014.

Comments: 4 pages, 4 figures, to appear in proceeding of IAU 311 meeting on Galaxy Masses as Constraints for Formation Models

Journal ref: IAU Symp 311, Galaxy Masses as Constraints of Formation Models (2015), ed. M. Cappellari & S. Courteau, p. 16-19

Cosmological Parameters From Supernovae Associated With Gamma-ray Bursts

Authors: Xue Li, Jens Hjorth, Radosław Wojtak

Abstract: We report estimates of the cosmological parameters $Ω_m$ and $Ω_Λ$ obtained using supernovae (SNe) associated with gamma-ray bursts (GRBs) at redshifts up to 0.606. Eight high-fidelity GRB-SNe with well-sampled light curves across the peak are used. We correct their peak magnitudes for a luminosity-decline rate relation to turn them into accurate standard candles with dispersion $σ= 0.18$ mag. We… ▽ More We report estimates of the cosmological parameters $Ω_m$ and $Ω_Λ$ obtained using supernovae (SNe) associated with gamma-ray bursts (GRBs) at redshifts up to 0.606. Eight high-fidelity GRB-SNe with well-sampled light curves across the peak are used. We correct their peak magnitudes for a luminosity-decline rate relation to turn them into accurate standard candles with dispersion $σ= 0.18$ mag. We also estimate the peculiar velocity of the low-redshift host galaxy of SN 1998bw, using constrained cosmological simulations. In a flat universe, the resulting Hubble diagram leads to best-fit cosmological parameters of $(Ω_m, Ω_Λ) = (0.58^{+0.22}_{-0.25},0.42 ^{+0.25}_{-0.22})$. This exploratory study suggests that GRB-SNe can potentially be used as standardizable candles to high redshifts to measure distances in the universe and constrain cosmological parameters. △ Less

Submitted 15 October, 2014; v1 submitted 11 September, 2014; originally announced September 2014.

Comments: 7 pages, 3 figures, Accepted for publication in ApJL

Radial orbit instability in dwarf dark matter haloes

Authors: Grzegorz Gajda, Ewa L. Lokas, Radoslaw Wojtak

Abstract: Using N-body simulations we study the phenomenon of radial orbit instability occurring in dark matter haloes of the size of a dwarf galaxy. We carried out simulations of seven spherical models, with the same standard NFW density profile but different anisotropy profiles of particle orbits. Four of them underwent instability: two with a constant positive anisotropy, one with an anisotropic core and… ▽ More Using N-body simulations we study the phenomenon of radial orbit instability occurring in dark matter haloes of the size of a dwarf galaxy. We carried out simulations of seven spherical models, with the same standard NFW density profile but different anisotropy profiles of particle orbits. Four of them underwent instability: two with a constant positive anisotropy, one with an anisotropic core and an isotropic envelope and one with a very small isotropic core and an anisotropic envelope. Haloes affected by the instability become approximately axisymmetric and prolate, with the profile of the shortest-to-longest axis ratio increasing with radius. The lower limit for the central value of this axis ratio is 0.3 for an NFW halo. The density profiles of the haloes did not change significantly, whereas the velocity distributions became axisymmetric. The total angular momentum rose due to large-amplitude oscillations of its components perpendicular to the symmetry axis of the halo. We also studied orbits of individual particles assigning them to classical orbit families in triaxial potentials. We find that the membership of a given particle in a family depends on its initial total angular momentum and its components along the principal axes of matter distribution. △ Less

Submitted 20 November, 2014; v1 submitted 24 July, 2014; originally announced July 2014.

Comments: 14 pages, 10 figures, accepted by MNRAS

Journal ref: MNRAS, 447, 97 (2015)

Effect of asphericity in caustic mass estimates of galaxy clusters

Authors: Jacob Svensmark, Radoslaw Wojtak, Steen H. Hansen

Abstract: The caustic technique for measuring mass profiles of galaxy clusters relies on the assumption of spherical symmetry. When applied to aspherical galaxy clusters, the method yields mass estimates affected by the cluster orientation. Here we employ mock redshift catalogues generated from cosmological simulations to study the effect of clusters intrinsic shape and surrounding filamentary structures on… ▽ More The caustic technique for measuring mass profiles of galaxy clusters relies on the assumption of spherical symmetry. When applied to aspherical galaxy clusters, the method yields mass estimates affected by the cluster orientation. Here we employ mock redshift catalogues generated from cosmological simulations to study the effect of clusters intrinsic shape and surrounding filamentary structures on the caustic mass estimates. To this end, we develop a new method for removing perturbations from large-scale structures, modelled as the two-halo term, in a caustic analysis of stacked cluster data. We find that the cluster masses inferred from kinematical data of ~10^14 Msun clusters observed along the major axis are larger than masses from those observed along the minor axis by a factor of 1.7 within the virial radius, increasing to 1.8 within three virial radii. This discrepancy increases by 20% for the most massive clusters. In addition a smaller but still significant mass discrepancy arises when filamentary structures are present near a galaxy cluster. We find that the mean cluster mass from random sightlines is unbiased at all radii and their scatter ranges from 0.14 to 0.17 within one and three virial radii, with a 40% increase for the most massive clusters. We provide tables which estimate the caustic mass bias given observational constraints on the cluster orientation. △ Less

Submitted 24 December, 2014; v1 submitted 1 May, 2014; originally announced May 2014.

Comments: 19 pages, 9 figures, 6 tables, accepted for publication in MNRAS

Galaxy Cluster Mass Reconstruction Project: I. Methods and first results on galaxy-based techniques

Authors: L. Old, R. A. Skibba, F. R. Pearce, D. Croton, S. I. Muldrew, J. C. Muñoz-Cuartas, D. Gifford, M. E. Gray, A. von der Linden, G. A. Mamon, M. R. Merrifield, V. Müller, R. J. Pearson, T. J. Ponman, A. Saro, T. Sepp, C. Sifón, E. Tempel, E. Tundo, Y. O. Wang, R. Wojtak

Abstract: This paper is the first in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our primary aim is to test the performance of these cluster mass estimation techniques on a diverse set of models that will increase in complexity. We begin by providing participating methods with… ▽ More This paper is the first in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our primary aim is to test the performance of these cluster mass estimation techniques on a diverse set of models that will increase in complexity. We begin by providing participating methods with data from a simple model that delivers idealised clusters, enabling us to quantify the underlying scatter intrinsic to these mass estimation techniques. The mock catalogue is based on a Halo Occupation Distribution (HOD) model that assumes spherical Navarro, Frenk and White (NFW) haloes truncated at R_200, with no substructure nor colour segregation, and with isotropic, isothermal Maxwellian velocities. We find that, above 10^14 M_solar, recovered cluster masses are correlated with the true underlying cluster mass with an intrinsic scatter of typically a factor of two. Below 10^14 M_solar, the scatter rises as the number of member galaxies drops and rapidly approaches an order of magnitude. We find that richness-based methods deliver the lowest scatter, but it is not clear whether such accuracy may simply be the result of using an over-simplistic model to populate the galaxies in their haloes. Even when given the true cluster membership, large scatter is observed for the majority non-richness-based approaches, suggesting that mass reconstruction with a low number of dynamical tracers is inherently problematic. △ Less

Submitted 18 March, 2014; originally announced March 2014.

Comments: 25 pages, 15 figures, 5 tables. Accepted for publication in MNRAS

Statistical mechanics of collisionless orbits. IV. Distribution of angular momentum

Authors: Liliya L. R. Williams, Jens Hjorth, Radoslaw Wojtak

Abstract: It has been shown in previous work that DARKexp, which is a theoretically derived, maximum entropy, one shape parameter model for isotropic collisionless systems, provides very good fits to simulated and observed dark-matter halos. Specifically, it fits the energy distribution, N(E), and the density profiles, including the central cusp. Here, we extend DARKexp N(E) to include the distribution in a… ▽ More It has been shown in previous work that DARKexp, which is a theoretically derived, maximum entropy, one shape parameter model for isotropic collisionless systems, provides very good fits to simulated and observed dark-matter halos. Specifically, it fits the energy distribution, N(E), and the density profiles, including the central cusp. Here, we extend DARKexp N(E) to include the distribution in angular momentum, L^2, for spherically symmetric systems. First, we argue, based on theoretical, semi-analytical, and simulation results, that while dark-matter halos are relaxed in energy, they are not nearly as relaxed in angular momentum, which precludes using maximum entropy to uniquely derive N(E,L^2). Instead, we require that when integrating N(E,L^2) over squared angular momenta one retrieves the DARKexp N(E). Starting with a general expression for N(E,L^2) we show how the distribution of particles in L^2 is related to the shape of the velocity distribution function, VDF, and velocity anisotropy profile, β(r). We then demonstrate that astrophysically realistic halos, as judged by the VDF shape and β(r), must have linear or convex distributions in L^2, for each separate energy bin. The distribution in energy of the most bound particles must be nearly flat, and become more tilted in favor of radial orbits for less bound particles. These results are consistent with numerical simulations and represent an important step towards deriving the full distribution function for spherically symmetric dark-matter halos. △ Less

Submitted 20 January, 2014; originally announced January 2014.

Comments: 17 pages, 10 figures, accepted to ApJ. The first three in this series are: Paper I: arXiv:1010.0265, Paper II: arXiv:1010.0266, Paper III: arXiv:1010.0267

Journal ref: Astrophysical Journal 783 (2014) 13

Cosmic variance of the local Hubble flow in large-scale cosmological simulations

Authors: Radoslaw Wojtak, Alexander Knebe, William A. Watson, Ilian T. Iliev, Steffen Hess, David Rapetti, Gustavo Yepes, Stefan Gottloeber

Abstract: The increasing precision in the determination of the Hubble parameter has reached a per cent level at which large-scale cosmic flows induced by inhomogeneities of the matter distribution become non-negligible. Here we use large-scale cosmological N-body simulations to study statistical properties of the local Hubble parameter as measured by local observers. We show that the distribution of the loc… ▽ More The increasing precision in the determination of the Hubble parameter has reached a per cent level at which large-scale cosmic flows induced by inhomogeneities of the matter distribution become non-negligible. Here we use large-scale cosmological N-body simulations to study statistical properties of the local Hubble parameter as measured by local observers. We show that the distribution of the local Hubble parameter depends not only on the scale of inhomogeneities, but also on how one defines the positions of observers in the cosmic web and what reference frame is used. Observers located in random dark matter haloes measure on average lower expansion rates than those at random positions in space or in the centres of cosmic voids, and this effect is stronger from the halo rest frames compared to the CMB rest frame. We compare the predictions for the local Hubble parameter with observational constraints based on type Ia supernovae (SNIa) and CMB observations. Due to cosmic variance, for observers located in random haloes we show that the Hubble constant determined from nearby SNIa may differ from that measured from the CMB by 0.8 per cent at 1sigma statistical significance. This scatter is too small to significantly alleviate a recently claimed discrepancy between current measurements assuming a flat LCDM model. However, for observers located in the centres of the largest voids permitted by the standard LCDM model, we find that Hubble constant measurements from SNIa would be biased high by 5 per cent, rendering this tension inexistent in this extreme case. △ Less

Submitted 23 December, 2013; v1 submitted 1 December, 2013; originally announced December 2013.

Comments: 9 pages, 6 figures; accepted for publication in MNRAS; updated references

Journal ref: MNRAS, 438, 1805 (2014)

Anisotropic q-Gaussian velocity distributions in LambdaCDM halos

Authors: Leandro J. Beraldo e Silva, Gary A. Mamon, Manuel Duarte, Radoslaw Wojtak, Sébastien Peirani, Gwenaël Boué

Abstract: The velocity distribution function (VDF) of dark matter (DM) halos in $Λ$CDM dissipationless cosmological simulations, which must be non-separable in its radial and tangential components, is still poorly known. We present the first single-parameter, non-separable, anisotropic model for the VDF in $Λ$CDM halos, built from an isotropic $q$-Gaussian (Tsallis) VDF of the isotropic set of dimensionless… ▽ More The velocity distribution function (VDF) of dark matter (DM) halos in $Λ$CDM dissipationless cosmological simulations, which must be non-separable in its radial and tangential components, is still poorly known. We present the first single-parameter, non-separable, anisotropic model for the VDF in $Λ$CDM halos, built from an isotropic $q$-Gaussian (Tsallis) VDF of the isotropic set of dimensionless spherical velocity components (after subtraction of streaming motions), normalized by the respective velocity dispersions. We test our VDF on 90 cluster-mass halos of a dissipationless cosmological simulation. Beyond the virial radius, $r_{\rm vir}$, our model VDF adequately reproduces that measured in the simulated halos, but no $q$-Gaussian model can adequately represent the VDF within $r_{\rm vir}$, as the speed distribution function is then flatter-topped than any $q$-Gaussian can allow. Nevertheless, our VDF fits significantly better the simulations than the commonly used Maxwellian (Gaussian) distribution, at virtually all radii within $5\,r_{\rm vir}$. Within 0.4 (1) $r_{\rm vir}$, the non-Gaussianity index $q$ is (roughly) linearly related to the slope of the density profile and also to the velocity anisotropy profile. We provide a parametrization of the modulation of $q$ with radius for both the median fits and the fit of the stacked halo. At radii of a few percent of $r_{\rm vir}$, corresponding to the Solar position in the Milky Way, our best-fit VDF, although fitting better the simulations than the Gaussian one, overproduces significantly the fraction of high velocity objects, indicating that one should not blindly use these $q$-Gaussian fits to make predictions on the direct detection rate of DM particles. △ Less

Submitted 20 December, 2016; v1 submitted 24 October, 2013; originally announced October 2013.

Comments: This version consolidates the published version and the Erratum (changes in red)

Journal ref: MNRAS 452, 944 (2015)

Phase-space shapes of clusters and rich groups of galaxies

Authors: Radoslaw Wojtak

Abstract: Clusters and groups of galaxies are highly aspherical, with shapes approximated by nearly prolate ellipsoids of revolution. An equally fundamental property is the shape of these objects in velocity space which is the anisotropy of the global velocity dispersion tensor. Here we make use of kinematical data comprising around 600 nearby clusters and rich groups of galaxies from the SDSS to place cons… ▽ More Clusters and groups of galaxies are highly aspherical, with shapes approximated by nearly prolate ellipsoids of revolution. An equally fundamental property is the shape of these objects in velocity space which is the anisotropy of the global velocity dispersion tensor. Here we make use of kinematical data comprising around 600 nearby clusters and rich groups of galaxies from the SDSS to place constraints on the phase-space shapes of these objects, i.e. their shapes in both position and velocity space. We show that the line of sight velocity dispersion normalised by a mass dependent velocity scale correlates with the apparent elongation, with circular (elongated) clusters exhibiting an excessive (decremental) normalised velocity dispersion. This correlation holds for dynamically young or old clusters and, therefore, it originates from projecting their intrinsic phase-space shapes rather than from dynamical evolution. It signifies that clusters are preferentially prolate not only in position space, but also in velocity space. The distribution of the axial ratios in position space is found to be well approximated by a Gaussian with a mean 0.66+/-0.01 and a dispersion 0.07+/-0.008. The velocity ellipsoids representing the shapes in velocity space are more spherical, with a mean axial ratio of 0.78+/-0.03. This finding has important implications for mass measurements based on the line of sight velocity dispersion profiles in individual clusters. For typical axial ratios of the velocity ellipsoids in the analysed cluster sample, systematic errors on the mass estimates inferred from the line of sight velocity dispersions become comparable to statistical uncertainties for galaxy clusters with as few as 40 spectroscopic redshifts. △ Less

Submitted 29 January, 2014; v1 submitted 14 October, 2013; originally announced October 2013.

Comments: 9 pages, 7 figures; published in A&A; typo in eq. 5 corrected

Journal ref: Astronomy & Astrophysics, Volume 559, id.A89 (8 pages), 2013

A new method to measure the mass of galaxy clusters

Authors: Martina Falco, Steen H. Hansen, Radoslaw Wojtak, Thejs Brinckmann, Mikkel Lindholmer, Stefania Pandolfi

Abstract: The mass measurement of galaxy clusters is an important tool for the determination of cosmological parameters describing the matter and energy content of the Universe. However, the standard methods rely on various assumptions about the shape or the level of equilibrium of the cluster. We present a novel method of measuring cluster masses. It is complementary to most of the other methods, since it… ▽ More The mass measurement of galaxy clusters is an important tool for the determination of cosmological parameters describing the matter and energy content of the Universe. However, the standard methods rely on various assumptions about the shape or the level of equilibrium of the cluster. We present a novel method of measuring cluster masses. It is complementary to most of the other methods, since it only uses kinematical information from outside the virialized cluster. Our method identifies objects, as galaxy sheets or filaments, in the cluster outer region, and infers the cluster mass by modeling how the massive cluster perturbs the motion of the structures from the Hubble flow. At the same time, this technique allows to constrain the three-dimensional orientation of the detected structures with a good accuracy. We use a cosmological numerical simulation to test the method. We then apply the method to the Coma cluster, where we find two galaxy sheets, and measure the mass of Coma to be Mvir=(9.2\pm2.4)10^{14} Msol, in good agreement with previous measurements obtained with the standard methods. △ Less

Submitted 11 September, 2013; originally announced September 2013.

Comments: 10 pages, 12 figures, submitted to MNRAS

Dynamical signatures of infall around galaxy clusters: a generalized Jeans equation

Authors: Martina Falco, Gary A. Mamon, Radoslaw Wojtak, Steen H. Hansen, Stefan Gottlöber

Abstract: We study the internal kinematics of galaxy clusters in the region beyond the sphere of virialization. Galaxies around a virialized cluster are infalling towards the cluster centre with a non-zero mean radial velocity. We develop a new formalism for describing the dynamical state of clusters, by generalizing the standard Jeans formalism with the inclusion of the peculiar infall motions of galaxies… ▽ More We study the internal kinematics of galaxy clusters in the region beyond the sphere of virialization. Galaxies around a virialized cluster are infalling towards the cluster centre with a non-zero mean radial velocity. We develop a new formalism for describing the dynamical state of clusters, by generalizing the standard Jeans formalism with the inclusion of the peculiar infall motions of galaxies and the Hubble expansion as well as the contributions due to background cosmology. Using empirical fits to the radial profiles of density, mean radial velocity and velocity anisotropy of both a stacked cluster-mass halo and two isolated halos of a cosmological dark matter only simulation, we verify that our generalized Jeans equation correctly predicts the radial velocity dispersion out to 4 virial radii. We find that the radial velocity dispersion inferred from the standard Jeans equation is accurate up to 2 virial radii, but overestimated by \approx 20% for the stacked halo and by \approx 40% for the isolated halos, in the range \approx 2-3 virial radii. Our model depends on the logarithmic growth rate of the virial radius (function of halo mass or concentration), which we estimate in 7 different ways, and on the departure from self-similarity of the evolution of the peculiar velocity profile in virial units. △ Less

Submitted 22 February, 2016; v1 submitted 27 June, 2013; originally announced June 2013.

Comments: 13 pages, 11 figures, published in MNRAS

Ursa Major II - Reproducing the observed properties through tidal disruption

Authors: R. Smith, M. Fellhauer, G. N. Candlish, R. Wojtak, J. P. Farias, M. Blaña

Abstract: Recent deep photometry of the dwarf spheroidal Ursa Major II's morphology, and spectroscopy of individual stars, have provided a number of new constraints on its properties. With a velocity dispersion $\sim$6 km s$^{-1}$, and under the assumption that the galaxy is virialised, the mass-to-light ratio is found to be approaching $\sim$2000 - apparently heavily dark matter dominated. Using N-Body sim… ▽ More Recent deep photometry of the dwarf spheroidal Ursa Major II's morphology, and spectroscopy of individual stars, have provided a number of new constraints on its properties. With a velocity dispersion $\sim$6 km s$^{-1}$, and under the assumption that the galaxy is virialised, the mass-to-light ratio is found to be approaching $\sim$2000 - apparently heavily dark matter dominated. Using N-Body simulations, we demonstrate that the observed luminosity, ellipticity, irregular morphology, velocity gradient, and the velocity dispersion can be well reproduced through processes associated with tidal mass loss, and in the absence of dark matter. These results highlight the considerable uncertainty that exists in measurements of the dark matter content of Ursa Major II. The dynamics of the inner tidal tails, and tidal stream, causes the observed velocity dispersion of stars to be boosted to values of $>$5 km s$^{-1}$ ($>$20 km s$^{-1}$ at times). This effect is responsible for raising the velocity dispersion of our model to the observed values in UMaII. We test an iterative rejection technique for removing unbound stars from samples of UMaII stars whose positions on the sky, and line-of-sight velocities, are provided. We find this technique is very effective at providing an accurate bound mass from this information, and only fails when the galaxy has a bound mass less than 10$%$ of its initial mass. However when $<2%$ mass remains bound, mass overestimation by $>$3 orders of magnitude are seen. Additionally we find that mass measurements are sensitive to measurement uncertainty in line-of-sight velocities. Measurement uncertainties of 1-4 km s$^{-1}$ result in mass overestimates by a factor of $\sim$1.3-5.7. △ Less

Submitted 23 May, 2013; originally announced May 2013.

Comments: 17 pages, 12 figures, accepted to MNRAS: 23rd, May, 2013

Orbital anisotropy in cosmological haloes revisited

Authors: Radoslaw Wojtak, Stefan Gottloeber, Anatoly Klypin

Abstract: The velocity anisotropy of particles inside dark matter (DM) haloes is an important physical quantity, which is required for the accurate modelling of mass profiles of galaxies and clusters of galaxies. It is typically measured using the ratio of the radial-to-tangential velocity dispersions at a given distance from the halo centre. However, this measure is insufficient to describe the dynamics of… ▽ More The velocity anisotropy of particles inside dark matter (DM) haloes is an important physical quantity, which is required for the accurate modelling of mass profiles of galaxies and clusters of galaxies. It is typically measured using the ratio of the radial-to-tangential velocity dispersions at a given distance from the halo centre. However, this measure is insufficient to describe the dynamics of realistic haloes, which are typically quite elongated. Studying the velocity distribution in massive DM haloes in cosmological simulations, we find that in the inner parts of the haloes the local velocity ellipsoids are strongly aligned with the major axis of the halo, the alignment being stronger for more relaxed haloes. In the outer regions of the haloes, the alignment becomes gradually weaker and the orientation is more random. These two distinct regions of different degree of the alignment coincide with two characteristic regimes of the DM density profile: shallower and steeper than ρr^{-2}. This alignment of the local velocity ellipsoids requires reinterpretation of features found in measurements based on the spherically averaged ratio of the radial-to-tangential velocity dispersions. In particular, we show that the velocity distribution in the central halo regions is highly anisotropic. For cluster-size haloes with mass 10^{14}-10^{15} h^-1 Msun, the velocity anisotropy along the major axis is nearly independent of radius and is equal to β=1-σ^2_{perp}/σ^2_{radial}=0.4, which is significantly larger than the previously estimated spherically averaged velocity anisotropy. The alignment of density and velocity anisotropies, and the radial trends may also have some implications for the mass modelling based on kinematical data of such objects as galaxy clusters or dwarf spheroidals, where the orbital anisotropy is a key element in an unbiased mass inference. △ Less

Submitted 21 June, 2013; v1 submitted 8 March, 2013; originally announced March 2013.

Comments: 11 pages, 8 figures; accepted for publication in MNRAS

Why does the Jeans Swindle work?

Authors: M. Falco, S. H. Hansen, R. Wojtak, G. A. Mamon

Abstract: When measuring the mass profile of any given cosmological structure through internal kinematics, the distant background density is always ignored. This trick is often refereed to as the "Jeans Swindle". Without this trick a divergent term from the background density renders the mass profile undefined, however, this trick has no formal justification. We show that when one includes the expansion of… ▽ More When measuring the mass profile of any given cosmological structure through internal kinematics, the distant background density is always ignored. This trick is often refereed to as the "Jeans Swindle". Without this trick a divergent term from the background density renders the mass profile undefined, however, this trick has no formal justification. We show that when one includes the expansion of the Universe in the Jeans equation, a term appears which exactly cancels the divergent term from the background. We thereby establish a formal justification for using the Jeans Swindle. △ Less

Submitted 9 January, 2013; v1 submitted 11 October, 2012; originally announced October 2012.

Comments: 5 pages, 2 figures, Accepted for publication in MNRAS Letters

Spatial Anisotropy of Galaxy Kinematics in Sloan Digital Sky Survey Galaxy Clusters

Authors: Andreas Skielboe, Radoslaw Wojtak, Kristian Pedersen, Eduardo Rozo, Eli S. Rykoff

Abstract: Measurements of galaxy cluster kinematics are important in understanding the dynamical state and evolution of clusters of galaxies, as well as constraining cosmological models. While it is well established that clusters exhibit non-spherical geometries, evident in the distribution of galaxies on the sky, azimuthal variations of galaxy kinematics within clusters have yet to be observed. Here we mea… ▽ More Measurements of galaxy cluster kinematics are important in understanding the dynamical state and evolution of clusters of galaxies, as well as constraining cosmological models. While it is well established that clusters exhibit non-spherical geometries, evident in the distribution of galaxies on the sky, azimuthal variations of galaxy kinematics within clusters have yet to be observed. Here we measure the azimuthal dependence of the line-of-sight velocity dispersion profile in a stacked sample of 1743 galaxy clusters from the Sloan Digital Sky Survey (SDSS). The clusters are drawn from the SDSS DR8 redMaPPer catalog. We find that the line-of-sight velocity dispersion of galaxies lying along the major axis of the central galaxy is larger than those that lie along the minor axis. This is the first observational detection of anisotropic kinematics of galaxies in clusters. We show that the result is consistent with predictions from numerical simulations. Furthermore we find that the degree of projected anisotropy is strongly dependent on the line-of-sight orientation of the galaxy cluster, opening new possibilities for assessing systematics in optical cluster finding. △ Less

Submitted 24 September, 2012; v1 submitted 22 August, 2012; originally announced August 2012.

Comments: 6 pages, 4 figures; final version as published in ApJ Letters

Journal ref: 2012 ApJ 758 L16

Physical properties underlying observed kinematics of satellite galaxies

Authors: Radoslaw Wojtak, Gary A. Mamon

Abstract: We study the kinematics of satellites around isolated galaxies selected from the Sloan Digital Sky Survey (SDSS) spectroscopic catalog. Using a model of the phase-space density previously measured for the halos of LCDM dark matter cosmological simulations, we determine the properties of the halo mass distribution and the orbital anisotropy of the satellites as a function of the colour-based morpho… ▽ More We study the kinematics of satellites around isolated galaxies selected from the Sloan Digital Sky Survey (SDSS) spectroscopic catalog. Using a model of the phase-space density previously measured for the halos of LCDM dark matter cosmological simulations, we determine the properties of the halo mass distribution and the orbital anisotropy of the satellites as a function of the colour-based morphological type and the stellar mass of the central host galaxy. We place constraints on the halo mass and the concentration parameter of dark matter and the satellite number density profiles. We obtain a concentration-mass relation for galactic dark matter haloes that is consistent with predictions of a standard LCDM cosmological model. At given halo or stellar mass, red galaxies have more concentrated halos than their blue counterparts. The fraction of dark matter within a few effective radii is minimal for 11.25<log M_star<11.5. The number density profile of the satellites appears to be shallower than of dark matter, with the scale radius typically 60 per cent larger than of dark matter. The orbital anisotropy around red hosts exhibits a mild excess of radial motions, in agreement with the typical anisotropy profiles found in cosmological simulations, whereas blue galaxies are found to be consistent with an isotropic velocity distribution. Our new constraints on the halo masses of galaxies are used to provide analytic approximations of the halo-to-stellar mass relation for red and blue galaxies. △ Less

Submitted 4 January, 2013; v1 submitted 6 July, 2012; originally announced July 2012.

Comments: 12 pages, 11 figures; accepted for publication in MNRAS (updated references)

Journal ref: Mon.Not.Roy.Astron.Soc. 428 (2013) 2407-2417

Gravitational redshift of galaxies in clusters as predicted by general relativity

Authors: Radoslaw Wojtak, Steen H. Hansen, Jens Hjorth

Abstract: The theoretical framework of cosmology is mainly defined by gravity, of which general relativity is the current model. Recent tests of general relativity within the ΛCold Dark Matter (CDM) model have found a concordance between predictions and the observations of the growth rate and clustering of the cosmic web. General relativity has not hitherto been tested on cosmological scales independent of… ▽ More The theoretical framework of cosmology is mainly defined by gravity, of which general relativity is the current model. Recent tests of general relativity within the ΛCold Dark Matter (CDM) model have found a concordance between predictions and the observations of the growth rate and clustering of the cosmic web. General relativity has not hitherto been tested on cosmological scales independent of the assumptions of the ΛCDM model. Here we report observation of the gravitational redshift of light coming from galaxies in clusters at the 99 per cent confidence level, based upon archival data. The measurement agrees with the predictions of general relativity and its modification created to explain cosmic acceleration without the need for dark energy (f(R) theory), but is inconsistent with alternative models designed to avoid the presence of dark matter. △ Less

Submitted 29 September, 2011; originally announced September 2011.

Comments: Published in Nature issued on 29 September 2011. This version includes the Letter published there as well as the Supplementary Information. 23 pages, 7 figures

Journal ref: Nature 477:567-569, 2011