Redshift drift in a universe with structure III: Numerical relativity

Basic Information

• Paper ID: 2404.06242
• Title: Redshift drift in a universe with structure III: Numerical relativity
• Authors: Sofie Marie Koksbang, Asta Heinesen, Hayley J. Macpherson
• Classification: astro-ph.CO (Cosmology and Nongalactic Astrophysics)
• Publication Date: April 2024 (arXiv v2: August 9, 2024)
• Paper Link: https://arxiv.org/abs/2404.06242

Abstract

This paper presents the first calculation of cosmological redshift drift in fully relativistic cosmological simulations using the Einstein Toolkit. Redshift drift is computed across the full sky for 50 synthetic observers and compared with the Einstein-de Sitter (EdS) model representing the large-scale spatial average spacetime of the simulation. Results demonstrate that at z≈0.2, the mean redshift drift of all observers deviates from EdS predictions at the percent level, decreasing to ~0.1% at z≈1. However, fluctuations in redshift drift across the sky reach ~10-30% at z≈0.1 and several percentage points at z≈0.5, potentially exceeding the expected precision of forthcoming redshift drift measurements. Additionally, on 0.48% of light ray paths, redshift drift temporarily becomes positive at very low redshifts z≲0.02.

Research Background and Motivation

Problem Definition

Redshift drift refers to the minute temporal variation of the redshift of a source object, enabling direct detection of cosmic expansion. This represents the first observational quantity capable of directly measuring cosmic acceleration or deceleration, unlike indirect methods such as supernovae observations.

Research Significance

1. Direct Detection of Cosmic Expansion: Redshift drift provides a unique opportunity for directly measuring the cosmic expansion history
2. Next-Generation Telescope Capabilities: Facilities such as SKA and ELT will soon possess the precision necessary to measure redshift drift
3. Theoretical Preparation: Theoretical preparation is needed to interpret future observations, particularly considering the effects of realistic cosmic structure

Limitations of Existing Methods

• Previous studies primarily relied on analytical calculations within the FLRW framework or Newtonian N-body simulations combined with perturbation theory
• Lack of research considering nonlinear structure effects within a fully relativistic framework
• Insufficient understanding of how structure affects redshift drift measurements

Core Contributions

1. First Fully Relativistic Calculation: First computation of redshift drift in fully relativistic cosmological simulations
2. Systematic Statistical Analysis: Statistical analysis of 38,400 light rays from 50 observers
3. Quantification of Structure Effects: Quantification of the degree to which cosmic structure affects redshift drift measurements
4. Observational Predictions: Provision of theoretical expectations for redshift drift observations with next-generation facilities such as SKA
5. Discovery of Positive Redshift Drift: Discovery that positive redshift drift may still occur in matter-dominated universes

Methodology

Task Definition

Computation of the temporal variation rate of redshift as photons propagate along geodesics in spacetime containing realistic cosmic structure, with analysis of its statistical properties and spatial distribution.

Simulation Architecture

Einstein Toolkit Numerical Relativity Simulations

• Evolution of Einstein equations using 3+1 BSSN formalism
• Initial conditions: Linear perturbations on EdS background at z̄=1000
• Cosmological fluid: Perfect dust fluid (P≪ρ)
• Boundary conditions: Periodic boundary conditions
• Grid resolution: Primary results use N=256³, box size L=3072h⁻¹Mpc

Ray-Tracing Method

Post-processing using the mescaline code:

• Random placement of 50 observers within the simulation domain
• Emission of 768 rays per observer along HEALPix pixel directions
• Solution of geodesic and geodesic deviation equations
• Both observers and sources comoving with the fluid

Redshift Drift Calculation

Redshift drift formula employing multipole decomposition:

$\frac{dz}{d\tau}\bigg|_O = E_E \int_{\lambda_O}^{\lambda_E} d\lambda \, \Pi$

where the integral term simplifies to:

$\Pi \approx -\frac{1}{3}u_\mu u_\nu R^{\mu\nu} - e_\mu e_\nu u_\rho u_\sigma C^{\rho\mu\sigma\nu}$

comprising contributions from both the Ricci tensor and Weyl tensor.

Technical Innovations

1. Fully Nonlinear Evolution: Independent of perturbation theory, incorporating complete nonlinear effects of structure formation
2. Multi-Observer Statistics: Systematic investigation of redshift drift variations observed by different-location observers
3. Full-Sky Analysis: Computation of redshift drift distribution and power spectrum across the full sky for each observer
4. Precise Numerical Implementation: High-precision implementation of redshift drift calculations in mescaline

Experimental Setup

Simulation Parameters

• Primary Simulation: N=256³, L=3072h⁻¹Mpc
• Convergence Tests: N=128³, L=1536h⁻¹Mpc
• Cosmological Parameters: h=0.45, other parameters adopt Planck defaults
• Evolution Range: From z̄=1000 to z̄≈0

Observational Setup

• Number of Observers: 50 (primary), 10 (high-resolution sky)
• Number of Rays: 768 per observer (Nside=8), 12,288 per observer (Nside=32)
• Observation Time: δτ₀=30 years
• Redshift Range: z≈0 to z≈1

Evaluation Metrics

• Redshift drift deviation relative to EdS model: δz/δzEdS-1
• Statistical quantities: mean, 68.1% and 95.4% confidence intervals
• Angular power spectrum analysis
• Sky map analysis

Experimental Results

Primary Results

Mean Redshift Drift

• For z>0.03, mean redshift drift deviation from EdS model <0.2%
• At z≈1, 95.4% confidence interval is ±2.5%, 68.1% confidence interval is ±1%
• Maximum fluctuations exceed 10% at z≈1

Compositional Analysis of Structure Effects

• Ricci contribution tends toward negative values
• Weyl contribution tends toward positive values
• The two approximately cancel statistically, but do not completely cancel along individual light rays

Positive Redshift Drift Phenomenon

• 0.48% of light rays exhibit positive redshift drift at z≲0.02
• Distributed across 185 rays from 5 observers
• Physical mechanism: Weyl contribution dominates at low redshift

Sky Distribution Analysis

Anisotropy Characteristics

• Sky fluctuations reach tens of percentage points at z≈0.1
• Decrease to several percentage points at z≈0.5
• Dipole contribution dominates, exceeding quadrupole by an order of magnitude

Angular Power Spectrum

• At z≈0.1, mean dipole amplitude is 3.7%, quadrupole is 1.1%
• At z≈0.5, dipole amplitude is 0.8%, quadrupole is 0.1%
• Large-angle multipoles dominate

Observer-to-Observer Differences

Mean redshift drift from observers at different locations deviates from EdS predictions by ~1% level at z≈0.2, which is significant for high-precision measurements.

Related Work

Theoretical Studies

• Analytical inhomogeneous cosmological models (LTB, Szekeres, etc.)
• Perturbation theory calculations
• Anisotropic model studies

Numerical Studies

• Newtonian N-body simulations combined with perturbation theory
• This paper represents the first fully relativistic numerical study

Observational Prospects

• SKA expected precision: 1-10%
• ELT's redshift drift measurement capabilities
• Results presented here indicate structure effects may reach observational precision levels

Conclusions and Discussion

Main Conclusions

1. Small Average Effects: Structure's impact on mean redshift drift is at the sub-percent level at cosmological redshifts
2. Significant Fluctuations: Redshift drift fluctuations along individual lines of sight can exceed 10%
3. Observational Challenges: Fluctuation levels approach or exceed the expected precision of facilities such as SKA
4. Positive Drift Phenomenon: Positive redshift drift may be observed even in matter-dominated universes

Limitations

1. No Cosmological Constant: Λ=0 in simulations; structure growth would be suppressed in the real universe
2. Source Position Drift Neglected: Position drift effects of sources not included (expected <1% contribution)
3. Resolution Constraints: Small-scale structure may affect low-redshift results
4. Limited Statistical Sample: Relatively limited number of observers and rays

Future Directions

1. Inclusion of Dark Energy: Repeat calculations within ΛCDM framework
2. Constrained Simulations: Use simulations constrained by local environment
3. Enhanced Statistics: Increase number of observers and rays
4. Observational Strategy: Provide selection strategies for actual observations

In-Depth Evaluation

Strengths

1. Pioneering Work: First calculation of redshift drift within fully relativistic framework, filling an important gap
2. Rigorous Methodology: Utilizes mature Einstein Toolkit and mescaline tool chain
3. Sufficient Statistics: 50 observers and 38,400 rays provide good statistical coverage
4. Important Results: Provides critical theoretical expectations for forthcoming observations
5. Complete Verification: Includes LTB model validation and convergence tests

Limitations

1. Physical Assumptions: Neglect of cosmological constant may overestimate structure effects
2. Computational Approximations: Neglect of position drift and certain higher-order terms
3. Resolution Dependence: Low-redshift results sensitive to numerical resolution
4. Statistical Constraints: Larger error bars for certain statistical quantities

Impact

1. Theoretical Contribution: Establishes numerical foundation for redshift drift theory
2. Observational Guidance: Provides important reference for projects such as SKA and ELT
3. Methodological Exemplar: Demonstrates application of numerical relativity to cosmological observations
4. Subsequent Research: Will stimulate further related theoretical and observational studies

Applicable Scenarios

1. Observational Preparation: High-precision redshift drift measurements with SKA Phase 2, etc.
2. Theoretical Validation: Testing cosmological models and dark energy theories
3. Systematic Errors: Assessment of structure effects' impact on cosmological parameter measurements
4. Method Development: Template for numerical relativistic calculations of other observables

References

The paper cites 62 relevant references, covering:

• Early theoretical work on redshift drift 1,2
• Observational prospects for SKA and ELT 3,4
• Redshift drift calculations in various cosmological models 5-14
• Einstein Toolkit numerical relativity tools 42,43
• Related cosmological simulation work 44-47

This work provides important theoretical preparation for the forthcoming era of redshift drift observations, quantifying for the first time the effects of realistic cosmic structure on this critical observable, with significant implications for both cosmology and observational astronomy.