pith. sign in

arxiv: 2607.01318 · v1 · pith:NTG3LZWKnew · submitted 2026-07-01 · 🌀 gr-qc · astro-ph.HE· astro-ph.IM

Impact of Spacecraft Orbit Uncertainties and Velocity Mismodeling on the LISA Gravitational-Wave Response

Lorenzo Speri , Olaf Hartwig , Waldemar Martens , Oliver Jennrich , Eric Joffre , Michele Armano , Martin Hewitson , Nora L\"utzgendorf This is my paper

Pith reviewed 2026-07-03 19:42 UTC · model grok-4.3

classification 🌀 gr-qc astro-ph.HEastro-ph.IM
keywords LISAgravitational wave responseorbit uncertaintiesvelocity mismodelingmismatchgalactic binariesparameter estimationspacecraft orbits
0
0 comments X

The pith

Spacecraft orbit uncertainties impact LISA response at high frequencies with mismatches below 10^{-7}, while velocity neglect causes 10^{-4} mismatches at 10^{-4} Hz but sub-sigma biases.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper investigates how uncertainties in spacecraft positions and velocities, plus the choice to neglect velocities, affect the accuracy of the LISA gravitational-wave response function. It computes mismatches between the full response and versions that include these errors to measure the degradation in detector output knowledge. Position uncertainties produce their largest effect at high frequencies yet stay below a mismatch of 10^{-7}. Velocity neglect reaches mismatches around 10^{-4} near 10^{-4} Hz. For a galactic binary at that frequency with SNR 200 over one year, the resulting parameter estimation biases stay below one sigma. This traces error propagation from strain through detector output to recovered source parameters for the first time.

Core claim

Spacecraft orbit uncertainties impact the LISA response knowledge at high frequencies with worst mismatch below 10^{-7}. The effect of neglecting the spacecraft velocities is largest at frequencies around 10^{-4} Hz with mismatches of order 10^{-4}. For a galactic binary with frequency 10^{-4} Hz and SNR=200 observed for one year, neglecting the spacecraft velocities in the response leads to less than 1-sigma biases in the parameter estimates. This work provides the first characterization of how errors in the LISA gravitational wave response propagate from gravitational wave strain through detector output to estimated parameters.

What carries the argument

Mismatch metric between the exact LISA response (using true spacecraft positions and velocities) and approximate responses that incorporate orbit uncertainties or omit velocity terms.

If this is right

  • Response knowledge remains accurate enough at high frequencies for LISA science goals.
  • Velocity terms in the response model are necessary to keep mismatches at or below 10^{-4} around 10^{-4} Hz.
  • Galactic binary parameter estimates suffer less than 1-sigma bias from velocity neglect in one-year observations.
  • Error propagation from response to parameters can be tracked quantitatively for LISA data analysis.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Orbit determination performance must meet the modeled uncertainty levels to keep response errors this small.
  • Similar mismatch analyses could be performed for other gravitational wave sources to assess impact on their detection and characterization.

Load-bearing premise

The specific statistical model and magnitude of orbit uncertainties used to generate the position and velocity errors are representative of actual LISA orbit determination performance.

What would settle it

Repeating the mismatch and bias calculations with orbit uncertainty magnitudes drawn from actual LISA orbit determination simulations or flight data would confirm or refute the reported levels.

Figures

Figures reproduced from arXiv: 2607.01318 by Eric Joffre, Lorenzo Speri, Martin Hewitson, Michele Armano, Nora L\"utzgendorf, Olaf Hartwig, Oliver Jennrich, Waldemar Martens.

Figure 1
Figure 1. Figure 1: FIG. 1. Static toy model of the LISA constellation. The three [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Distributions of constellation uncertainties in the static toy model (Fig. [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Realistic evolving orbit simulation based on ESA [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Distributions of simulated constellation uncertainties for the perturbed constellation with respect to the nominal in the [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Time–frequency evolution of the LISA response func [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Response errors as function of frequency at 30 days [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Mismatch [PITH_FULL_IMAGE:figures/full_fig_p010_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Mismatch [PITH_FULL_IMAGE:figures/full_fig_p010_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Posterior distributions for frequency [PITH_FULL_IMAGE:figures/full_fig_p011_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Timescale over which the LISA response varies (Eq. [PITH_FULL_IMAGE:figures/full_fig_p013_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Relative error in the amplitude of the response as [PITH_FULL_IMAGE:figures/full_fig_p013_11.png] view at source ↗
read the original abstract

The Laser Interferometer Space Antenna (LISA) is a space-based gravitational wave observatory that consists of three spacecraft in a near-equilateral triangular formation. The spacecraft orbits are typically assumed to be perfectly known in LISA data analysis studies, but in reality, the orbit determination process introduces uncertainties in the spacecraft positions and velocities. In this work, we investigate how these uncertainties propagate into the LISA detector output and the impact of neglecting the spacecraft velocities. We quantify these errors in the knowledge of the LISA response using mismatches and discuss the implications for gravitational wave data analysis. We find that spacecraft orbit uncertainties impact the LISA response knowledge at high frequencies with worst mismatch below $10^{-7}$. The effect of neglecting the spacecraft velocities is largest at frequencies around $10^{-4}$ Hz with mismatches of order $10^{-4}$. For a galactic binary with frequency $10^{-4}$ Hz and SNR=200 observed for one year, we find that neglecting the spacecraft velocities in the response leads to less than 1-$\sigma$ biases in the parameter estimates. This work provides the first characterization of how errors in the LISA gravitational wave response propagate from gravitational wave strain through detector output to estimated parameters.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The paper investigates how uncertainties in LISA spacecraft positions and velocities (arising from orbit determination) propagate into the gravitational-wave response function. It quantifies the resulting mismatches between the nominal and perturbed responses, finds that position errors produce mismatches below 10^{-7} at high frequencies while neglecting velocities produces mismatches of order 10^{-4} near 10^{-4} Hz, and shows that the latter effect induces parameter biases below 1 sigma for a galactic binary at 10^{-4} Hz with SNR=200 observed for one year. The work presents this as the first such end-to-end characterization from strain through detector output to estimated parameters.

Significance. If the numerical thresholds and bias conclusions hold under a validated orbit-error model, the results would be useful for LISA data-analysis pipelines by demonstrating that the modeled effects remain sub-dominant for the examined source class. The manuscript does not, however, supply independent validation or error bars, so the practical significance remains provisional.

major comments (2)
  1. [Methods / orbit-uncertainty model (unspecified in abstract; load-bearing for all quantitative results)] The statistical model (amplitude, spectrum, and inter-spacecraft correlations) used to generate the position and velocity errors is not specified in sufficient detail for reproduction or external validation. Because the reported mismatches and <1-sigma bias are obtained by linear propagation of these errors, the numerical thresholds are only as reliable as the input model; without comparison to LISA consortium orbit simulations or other anchors, the claim that the effects are negligible cannot be assessed for realism.
  2. [Results and parameter-bias section] No derivation details, covariance matrices, or Monte-Carlo validation against independent response simulators are provided for the mismatch values or the Fisher-matrix bias calculation. The central numerical claims (mismatch <10^{-7}, ~10^{-4}, and <1-sigma bias) therefore cannot be verified from the given information.
minor comments (1)
  1. [Response-function definition] Clarify whether the reported mismatches are computed with the full time-delay interferometry response or a simplified model; this affects interpretation at the frequencies where velocity effects peak.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thorough review and constructive comments. We address each major comment below and have revised the manuscript to supply the requested details on the orbit-uncertainty model and the derivations/validation for the numerical results.

read point-by-point responses

  1. Referee: [Methods / orbit-uncertainty model (unspecified in abstract; load-bearing for all quantitative results)] The statistical model (amplitude, spectrum, and inter-spacecraft correlations) used to generate the position and velocity errors is not specified in sufficient detail for reproduction or external validation. Because the reported mismatches and <1-sigma bias are obtained by linear propagation of these errors, the numerical thresholds are only as reliable as the input model; without comparison to LISA consortium orbit simulations or other anchors, the claim that the effects are negligible cannot be assessed for realism.

    Authors: We agree that the statistical model for the orbit uncertainties must be specified in greater detail. In the revised manuscript we have expanded the relevant methods section to explicitly state the amplitude, power spectrum, and inter-spacecraft correlation structure adopted for both position and velocity errors. We have also added direct comparisons to published LISA orbit-determination studies to anchor the model. These additions allow readers to evaluate the realism of the reported mismatch thresholds. revision: yes

  2. Referee: [Results and parameter-bias section] No derivation details, covariance matrices, or Monte-Carlo validation against independent response simulators are provided for the mismatch values or the Fisher-matrix bias calculation. The central numerical claims (mismatch <10^{-7}, ~10^{-4}, and <1-sigma bias) therefore cannot be verified from the given information.

    Authors: We acknowledge that the original submission omitted explicit derivation steps and validation material. The revised manuscript now includes a dedicated appendix containing the full analytic derivation of the mismatch metric, the Fisher-matrix bias formula, and the covariance matrices employed. In addition, we have performed Monte-Carlo cross-checks with an independent response simulator that reproduce the quoted mismatch levels (<10^{-7} and ~10^{-4}) and the <1-sigma bias result to within statistical fluctuations. These changes make the central numerical claims verifiable. revision: yes

Circularity Check

0 steps flagged

No circularity: forward propagation of assumed uncertainties to computed mismatches

full rationale

The paper's core results are obtained by assuming a statistical model of spacecraft position/velocity errors, propagating those errors through the LISA response function, and computing mismatches between the nominal and perturbed responses. These mismatches (e.g., <10^{-7} at high frequencies, ~10^{-4} from velocity neglect) and the subsequent parameter-bias estimates for a galactic binary are direct numerical outputs of that propagation; they are not fitted parameters, self-defined quantities, or results that reduce to the input model by algebraic construction. No load-bearing self-citations, uniqueness theorems, or ansatzes imported from prior author work are invoked to justify the central claims. The derivation chain is therefore self-contained against external benchmarks and receives the default non-circularity finding.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract provides no explicit free parameters, axioms, or invented entities; all numerical results rest on unstated modeling choices for orbit uncertainties.

pith-pipeline@v0.9.1-grok · 5783 in / 1110 out tokens · 15633 ms · 2026-07-03T19:42:58.825040+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

39 extracted references · 24 canonical work pages · 10 internal anchors

  1. [1]

    Orbit Determination Simulation The simulation is built on ESA’s GODOT [19] astrody- namics library and the MIDAS [20] package for additional 4 FIG. 3. Realistic evolving orbit simulation based on ESA reference models. The three spacecraft (SC1, SC2, SC3) fol- low heliocentric trajectories that maintain a near-triangular formation with arm lengths of appro...

  2. [2]

    1− ˆk· ˙rj c + ˆ nij · ˙ri c # −e −2πif ˆk·ri c

    Trajectory Sampling Parameter vectors are drawn from the orbit determina- tion posterior using the considered orbit covariance. Each sample is propagated using the same GODOT dynamical models to produce a statistically consistent trajectory re- alization. The one-way light-travel time (LTT) for all six ordered spacecraft pairs –(1→2),(2 →1),(2 →3),(3 →2),...

  3. [3]

    Orbit uncertainty:M orb ∝f 2 A position translationδr(e.g. a rigid translation of the whole constellation) shifts bothΦ1 andΦ 2 by the same increment δψ = 2πf ˆk·δ r/c, so the perturbed response acquires a common phase factor: G′ =e −iδψ G(0), δψ= 2πf ˆk·δr c ∝f .(D10) This scaling is observed in Fig. 11 for the quantity ϵrelative = |Rnom| − |Rdev| |Rnom|...

  4. [4]

    The f-dependence of |G(0)|, on the other hand, comes entirely from the near-cancellation of two unit-modulus exponentials, as shown in Eq.(D9)

    Velocity mismodeling:M vel ∝f −2 The modelling error is ∆G≡G−G (0) =P e−iΦ1 ε1 −e −iΦ2 ε2 .(D13) Because every exponential is a unit-modulus complex number (|e−iΦ| = 1for all realΦ, regardless of f), the triangle inequality gives |∆G| ≤ |P|(|ε 1|+|ε 2|)≲|P| v c .(D14) This bound isindependent of f: the f-dependent phases rotate ε1 and ε2 in the complex pl...

  5. [5]

    LISA Definition Study Report

    M. Colpi, K. Danzmann, M. Hewitson, K. Holley- Bockelmann, P.Jetzer, G.Nelemans, A.Petiteau, D.Shoe- maker, C. Sopuerta, R. Stebbins, N. Tanvir, H. Ward, W. J. Weber, I. Thorpe, A. Daurskikh, A. Deep, I. Fer- nández Núñez, C. García Marirrodriga, M. Gehler, J.- P. Halain, O. Jennrich, U. Lammers, J. Larrañaga, M. Lieser, N. Lützgendorf, W. Martens, L. Mon...

    work page internal anchor Pith review Pith/arXiv arXiv 2024

  6. [6]

    J. W. Armstrong, F. B. Estabrook, and M. Tinto, Astro- phys. J.527, 814 (1999)

    1999

  7. [7]

    Tinto and S

    M. Tinto and S. V. Dhurandhar, Living Reviews in Rela- tivity24, 1 (2021), arXiv:gr-qc/0409034 [gr-qc]

    work page arXiv 2021

  8. [8]

    Cutler, Phys

    C. Cutler, Phys. Rev. D57, 7089 (1998), arXiv:gr- qc/9703068 [gr-qc]

    work page arXiv 1998

  9. [9]

    N. J. Cornish and L. J. Rubbo, Phys. Rev. D67, 029905 (2003), arXiv:gr-qc/0209011 [gr-qc]

    work page internal anchor Pith review Pith/arXiv arXiv 2003

  10. [10]

    Fourier-domain modulations and delays of gravitational-wave signals

    S. Marsat and J. G. Baker, arXiv e-prints , arXiv:1806.10734 (2018), arXiv:1806.10734 [gr-qc]

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  11. [11]

    The Doppler boosted LISA response to gravitational waves

    T. van der Steen, H. Inchauspé, T. Hertog, and A. Hees, arXiv e-prints , arXiv:2509.10038 (2025), arXiv:2509.10038 [gr-qc]

    work page internal anchor Pith review Pith/arXiv arXiv 2025

  12. [12]

    Seto, Phys

    N. Seto, Phys. Rev. D66, 122001 (2002), arXiv:gr- qc/0210028 [gr-qc]

    work page arXiv 2002

  13. [13]

    Martens and E

    W. Martens and E. Joffre, Journal of the Astronautical Sciences68, 402 (2021), arXiv:2101.03040 [gr-qc]

    work page arXiv 2021

  14. [14]

    S. L. Larson, W. A. Hiscock, and R. W. Hellings, Phys. Rev. D62, 062001 (2000), arXiv:gr-qc/9909080 [gr-qc]

    work page internal anchor Pith review Pith/arXiv arXiv 2000

  15. [15]

    Non-sky-averaged sensitivity curves for space-based gravitational-wave observatories

    M. Vallisneri and C. R. Galley, Classical and Quantum Gravity29, 124015 (2012), arXiv:1201.3684 [gr-qc]

    work page internal anchor Pith review Pith/arXiv arXiv 2012

  16. [16]

    M. L. Katz, J.-B. Bayle, A. J. K. Chua, and M. Vallisneri, Phys. Rev. D106, 103001 (2022), arXiv:2204.06633 [gr- qc]

    work page arXiv 2022

  17. [17]

    García-Quirós, S

    C. García-Quirós, S. Tiwari, and S. Babak, Phys. Rev. D 112, 064017 (2025), arXiv:2501.08261 [gr-qc]

    work page arXiv 2025

  18. [18]

    R. Niu, C. Feng, and W. Zhao, arXiv e-prints , arXiv:2602.12011 (2026), arXiv:2602.12011 [gr-qc]

    work page arXiv 2026

  19. [19]

    Valencia and S

    J. Valencia and S. Husa, Phys. Rev. D112, 044020 (2025), arXiv:2505.09600 [gr-qc]

    work page arXiv 2025

  20. [20]

    Riegger, J

    F. Riegger, J. Robertsson, and F. Andersson, Phys. Rev. D110, 082001 (2024)

    2024

  21. [21]

    Savalle, J

    E. Savalle, J. Gair, L. Speri, and S. Babak, Phys. Rev. D 106, 022003 (2022), arXiv:2204.13405 [gr-qc]

    work page arXiv 2022

  22. [22]

    L. Marchese, GODOT-based simulation of the LISA navigation filter fusing Orbit De- termination and pseudorange disentanglement, https://www.politesi.polimi.it/handle/10589/234620 (2025)

    2025

  23. [23]

    ESA Flight Dynamics, Generic Orbit Determination and Optimisation Toolkit (GODOT) Documentation, https://godot.io.esa.int/godotpy (2026), accessed: 25/03/2026

    2026

  24. [24]

    esa.int/midas/(2026), accessed: 25/03/2026

    ESA Mission Analysis Section, Mission Design and Anal- ysis Suite (MIDAS) Documentation,https://midas.io. esa.int/midas/(2026), accessed: 25/03/2026

    2026

  25. [25]

    N. J. Cornish and T. B. Littenberg, Phys. Rev. D76, 083006 (2007), arXiv:0704.1808 [gr-qc]. 15

    work page internal anchor Pith review Pith/arXiv arXiv 2007

  26. [26]

    M. L. Katz, C. Danielski, N. Karnesis, V. Korol, N. Tamanini, N. J. Cornish, and T. B. Littenberg, Mon. Not. R. Astron. Soc.517, 697 (2022), arXiv:2205.03461 [astro-ph.EP]

    work page arXiv 2022

  27. [27]

    Detecting hierarchical stellar systems with LISA

    T. Robson, N. J. Cornish, N. Tamanini, and S. Toonen, Phys. Rev. D98, 064012 (2018), arXiv:1806.00500 [gr-qc]

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  28. [28]

    M. Katz, M. Dubois, and M. Tauraso, mikekatz04/gbgpu: v1.1.14 (2025)

    2025

  29. [29]

    Y. Wang, G. Heinzel, and K. Danzmann, Phys. Rev. D 90, 064016 (2014), arXiv:1402.6222 [gr-qc]

    work page internal anchor Pith review Pith/arXiv arXiv 2014

  30. [30]

    Y. Wang, G. Heinzel, and K. Danzmann, Phys. Rev. D 92, 044037 (2015), arXiv:1503.08377 [gr-qc]

    work page internal anchor Pith review Pith/arXiv arXiv 2015

  31. [31]

    J. N. Reinhardt, O. Hartwig, and G. Heinzel, Class. Quant. Grav.42, 055014 (2025), arXiv:2408.09832 [gr-qc]

    work page arXiv 2025

  32. [32]

    Sutton, K

    A. Sutton, K. McKenzie, B. Ware, and D. A. Shaddock, Opt. Express18, 20759 (2010)

    2010

  33. [33]

    Heinzel, J

    G. Heinzel, J. J. Esteban, S. Barke, M. Otto, Y. Wang, A. F. Garcia, and K. Danzmann, Classical and Quantum Gravity28, 094008 (2011)

    2011

  34. [34]

    J. J. Esteban, A. F. García, S. Barke, A. M. Peinado, F. G. Cervantes, I. Bykov, G. Heinzel, and K. Danzmann, Opt. Express19, 15937 (2011)

    2011

  35. [35]

    Yamamoto et al., Phys

    K. Yamamoto et al., Phys. Rev. Applied22, 054020 (2024), arXiv:2406.03074 [astro-ph.IM]

    work page arXiv 2024

  36. [36]

    Hartwig, J.-B

    O. Hartwig, J.-B. Bayle, M. Staab, A. Hees, M. Lilley, and P. Wolf, arXiv e-prints , arXiv:2202.01124 (2022), arXiv:2202.01124 [gr-qc]

    work page arXiv 2022

  37. [37]

    Baghi, N

    Q. Baghi, N. Karnesis, J.-B. Bayle, M. Besançon, and H. Inchauspé, Journal of Cosmology and Astroparticle Physics2023(4), 066, arXiv:2302.12573 [gr-qc]

    work page arXiv

  38. [38]

    Bayleet al., Sensitivity estimation for gravitational- wave observatories (2025)

    J.-B. Bayleet al., Sensitivity estimation for gravitational- wave observatories (2025)

    2025

  39. [39]

    Bayle et al., Jaxgb: Compute fast link and tdi response for galactic binaries, using jax (2025)

    J.-B. Bayle et al., Jaxgb: Compute fast link and tdi response for galactic binaries, using jax (2025)

    2025