% ========================================================================
% © 2026 José Antonio Sánchez Lázaro
%
% This work is licensed under the Creative Commons Attribution 4.0
% International License (CC BY 4.0).
%
% Original deposit:
% Zenodo DOI: https://doi.org/10.5281/zenodo.16235702
% Date: 20 April 2026 (version v1.0.3)
%
% You are free to share, adapt and use this material for any purpose,
% provided that appropriate credit is given to the original author,
% a link to the license is provided, and any changes are indicated.
%
% Full license: https://creativecommons.org/licenses/by/4.0/
% Contact: research@darcysoft.com
% ========================================================================

\section{Detailed Analysis of Perceptual t-Echoes in LIGO O5}
\label{sec:O5_t_echo_analysis}

This section provides a realistic, end-to-end forecast for the detectability of **perceptual t-echoes** (predicted amplitude \(\sim 10^{-3}\), characteristic delay \(\sim 10^{-20}\) s) during the LIGO-Virgo-KAGRA O5 observing run, using the official LVK 2025 noise curves, a realistic binary black hole (BBH) population, and both coherent and incoherent stacking strategies.

\subsection{1. Waveform Model for t-Echoes}

A perceptual t-echo arises because the 4D worldline of the merger continues its oscillation in the real temporal coordinate \(t\) after the main coalescence. When projected onto our locally biased +t slice, this produces a faint, delayed replica of the main gravitational-wave signal.

We model the echo as a scaled and phase-shifted copy of the main waveform:
\[
h_{\rm echo}(t) = A_{\rm echo} \cdot h_{\rm main}(t - \Delta t) \cdot e^{i\phi_{\rm jitter}},
\]
where:
- \(A_{\rm echo} = 10^{-3}\) (relative amplitude from elastic strain leakage),
- \(\Delta t \approx 10^{-20}\) s (t-offset set by the coherence length \(\lambda_{\rm coh} \approx 10^{10} l_P\)),
- \(\phi_{\rm jitter} \sim \mathcal{N}(0, 0.3\,{\rm rad})\) (phase jitter from mass and redshift uncertainty).

The main waveform is generated with the IMRPhenomXHM approximant (or a simplified chirp+ringdown for speed) for a typical 30+30 \(M_\odot\) BBH at luminosity distance \(D_L \approx 400\) Mpc (median O5 event).

\subsection{2. LIGO O5 Noise Curve (LVK 2025 Official Approximation)}

We use the mid-band O5 noise power spectral density (PSD) from the LVK Observing Run Plans white paper (T2400403-v2, 2025):
\[
S_n(f) \approx 6.25 \times 10^{-48} \left(\frac{f}{100\,{\rm Hz}}\right)^{-1} \left[1 + \left(\frac{f}{300\,{\rm Hz}}\right)^4\right] \quad (10\,{\rm Hz} < f < 3000\,{\rm Hz}).
\]
This corresponds to a strain sensitivity of \(\approx 2.5 \times 10^{-24}\,{\rm Hz}^{-1/2}\) at 100 Hz, consistent with the projected O5 curve.

\subsection{3. Single-Event Sensitivity}

For a typical 30+30 \(M_\odot\) BBH at 400 Mpc:
- Matched-filter SNR of the main signal: \(\rho_{\rm main} \approx 12\)--15 (median for O5 population).
- Expected SNR of the t-echo (before stacking): \(\rho_{\rm echo} \approx 0.012\)--0.015.

This is **below the single-event detection threshold** (\(\rho > 8\)), as expected.

\subsection{4. Realistic BBH Population in O5}

From the LVK 2025 white paper:
- Expected BBH detection rate: \(\approx 300\) events per year.
- O5 duration: 4 years (with 70\% duty cycle) \(\Rightarrow\) total \(\approx 840\)--1200 events.
- Fraction with high-quality parameters (for coherent stacking): \(\approx 55\)--65\%.

We adopt a conservative stackable sample of **\(N_{\rm stack} = 720\) events**.

\subsection{5. Stacking Strategies and Expected Significance}

\subsubsection{Coherent Stacking (phase-known)}
If the t-offset phase can be predicted from the main signal parameters (mass, spin, redshift) to within \(\sigma_\phi \approx 0.3\) rad, the stacking gain is
\[
G_{\rm coh} = \sqrt{N_{\rm stack}} \times e^{-\sigma_\phi^2 / 2} \approx 25.7.
\]
Stacked echo SNR:
\[
\rho_{\rm stacked} = 0.013 \times 25.7 \approx 0.33 \quad \text{(still marginal)}.
\]

\subsubsection{Incoherent Stacking (more robust)}
For unknown phase (or large jitter), we stack the excess power:
\[
\rho_{\rm inc}^2 = N_{\rm stack} \times \rho_{\rm single}^2 \quad \Rightarrow \quad \rho_{\rm inc} \approx 0.35.
\]
Still below 3\(\sigma\).

**Conclusion from baseline calculation**: With the nominal \(A_{\rm echo} = 10^{-3}\), even optimistic coherent stacking yields only \(\sim 0.3\sigma\) — **not detectable** in O5 under current assumptions.

\subsection{6. Parameter Space Where Detection Becomes Feasible}

We vary the echo amplitude \(A_{\rm echo}\) and delay \(\Delta t\) to find the region where 3\(\sigma\) or 5\(\sigma\) is reachable:

| \(A_{\rm echo}\) | Required \(N_{\rm stack}\) for 3\(\sigma\) (coherent) | Detectable in O5? |
|------------------|-------------------------------------------------------|-------------------|
| \(10^{-3}\)      | \(\gtrsim 50\,000\)                                   | No                |
| \(3 \times 10^{-3}\) | \(\approx 5\,500\)                                 | Marginal (if rate higher) |
| \(5 \times 10^{-3}\) | \(\approx 2\,000\)                                 | Yes (with 3G)     |
| \(10^{-2}\)      | \(\approx 500\)                                        | Yes (O5 possible) |

**Key insight**: The theory predicts \(A_{\rm echo} \sim 10^{-3}\). Detection in O5 is **unlikely** unless:
- The actual elastic leakage is larger (\(\sim 3\)--\(5 \times 10^{-3}\)), or
- We use third-generation detectors (Einstein Telescope / Cosmic Explorer) with \(\sim 10\times\) better sensitivity.

\subsection{7. Recommended Analysis Strategy for LVK}

1. **Targeted search** on the 720 highest-SNR events (known parameters).
2. Use a **two-stage pipeline**:
   - Stage 1: Incoherent excess-power stack (robust to phase jitter).
   - Stage 2: Coherent matched-filter stack with phase marginalization over \(\pm 0.5\) rad.
3. Background estimation via time-slides (10 000 realizations) to set false-alarm probability.
4. Cross-check with Virgo and KAGRA data (different noise curves and antenna patterns).

\subsection{8. Conclusions and Outlook}

- With the nominal theory prediction (\(A_{\rm echo} = 10^{-3}\)), **t-echoes are not expected to be detectable in LIGO O5** even with aggressive stacking.
- A factor of 3--5 improvement in effective echo amplitude (or the use of 3G detectors) would make them detectable at 3--5\(\sigma\).
- The non-detection in O5 would **not falsify** the theory, but would tighten the upper bound on the elastic leakage parameter to \(A_{\rm echo} \lesssim 4 \times 10^{-3}\).
- Positive detection in O5 or early 3G runs would constitute **smoking-gun evidence** for perceptual 4D Euclidean geometry.

This analysis is fully reproducible with the accompanying code `O5\_t\_echo\_stacking\_analysis.py` (available on Zenodo) and can be updated as soon as the final O5 noise curves and event catalogs are released.