Mitigating Systematic Errors in Parameter Estimation of Binary Black Hole Mergers in O1-O3 LIGO-Virgo Data
Pith reviewed 2026-05-19 16:57 UTC · model grok-4.3
The pith
Parametric uncertainty models for waveforms reduce systematic errors in LIGO-Virgo binary black hole parameter estimates and align results across models and data versions.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that equipping gravitational-wave parameter estimation with parametric phase and amplitude uncertainty models, together with sufficiently broad priors on those uncertainty parameters, absorbs the dominant systematic errors present in selected O1-O3 events. This absorption renders astrophysical inferences, such as spin precession and spin alignment, consistent across waveform models (IMRPhenomXPHM, IMRPhenomXO4a, NRSur7dq4) and across raw versus deglitched strain data, even when the original analyses disagreed because of glitches near the signal or model differences.
What carries the argument
Parametric models that introduce free parameters for deviations in the phase and amplitude of the gravitational waveform, used with broad priors to marginalize over systematic mismatches.
If this is right
- Inconsistent precession measurements for GW200129_065458 become consistent non-zero values of approximately 0.58 across three waveform models.
- Anti-aligned spin inferences for GW191109_010717 remain but agree between raw and deglitched frame files and across the same three waveform models.
- Systematic errors from glitches occurring near a signal or from the deglitching process are reduced without manual removal of data segments.
- Parameter estimation results become less sensitive to the choice of waveform approximant when the uncertainty parameters are included.
Where Pith is reading between the lines
- Population analyses that combine many events could adopt these uncertainty models as a standard step to reduce scatter in inferred spin distributions.
- Events that still show large uncertainty parameters after the fit could be flagged for deeper follow-up with more expensive numerical-relativity waveforms.
- The approach may generalize to other classes of signals, such as neutron-star mergers, where waveform modeling uncertainties are also large.
Load-bearing premise
The parametric uncertainty models can absorb the dominant systematic errors without introducing new biases into the astrophysical parameters.
What would settle it
Reanalysis of a simulated binary black hole signal with a known injected glitch and known true parameters; if the recovered astrophysical parameters remain biased away from the injected values even after the uncertainty parameters are marginalized, the method fails to mitigate the systematics.
Figures
read the original abstract
Systematic errors in the parameter estimation (PE) of gravitational wave (GW) mergers can arise from various sources, including waveform systematics, noise mischaracterization, data analysis artifacts, and other unknown factors. In this study, we analyze selected events from the first three observing runs of the LIGO-Virgo-KAGRA (LVK) collaboration. We choose events that have been flagged in various studies as potentially affected by systematic errors. Here, we reanalyze these events using a couple of parametric models developed in previous work that incorporate uncertainties in both the phase and amplitude of the GW waveform. In this data-driven approach, we apply sufficiently broad priors on the uncertainty parameters to account for potential systematic errors. Our findings show that the proposed method effectively reduces systematic errors, even those arising from data artifacts, such as glitches occurring near a signal and the deglitching process in GW frame files. Similarly, inconsistent results from different waveform models become much more consistent in our framework. One noteworthy event we examine is GW191109\_010717, which is particularly interesting due to its anti-aligned spin properties. We report that, within our framework, the event still exhibits anti-aligned spin characteristics, but the inference results become consistent across raw and deglitched frame files, as well as across the waveform models used for this event (IMRPhenomXPHM, IMRPhenomXO4a, and NRSur7dq4). A similar trend is observed for the event GW200129\_065458, which previously yielded a high, but inconsistent precession parameter among different waveform models. In contrast, we observe a non-zero and consistent value of $\chi_{p}=0.60^{+0.31}_{-0.33}, 0.58^{+0.30}_{-0.29}$ and $0.56^{+0.31}_{-0.28}$ for the IMRPhenomXPHM, IMRPhenomXO4a, and NRSur7dq4 waveform models, respectively.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes a data-driven method to mitigate systematic errors in parameter estimation of binary black hole mergers from LIGO-Virgo O1-O3 data. It reanalyzes selected events flagged for potential systematics (including GW191109_010717 and GW200129_065458) by augmenting standard waveform models with parametric phase and amplitude uncertainty models equipped with broad priors. The central claim is that this approach reduces discrepancies arising from waveform model differences and data artifacts such as glitches and deglitching, yielding consistent astrophysical inferences (e.g., non-zero precession parameter χ_p ≈ 0.58 across IMRPhenomXPHM, IMRPhenomXO4a, and NRSur7dq4 for GW200129_065458, and preserved anti-aligned spin for GW191109_010717 across raw and deglitched frames).
Significance. If validated, the approach could provide a practical route to more robust spin and precession measurements for events affected by known data-quality issues, complementing existing waveform-model comparisons. The work applies previously developed parametric uncertainty models to real LVK events and reports concrete overlapping credible intervals, which is a positive step toward falsifiable improvements in systematic-error control.
major comments (2)
- [Results section (GW200129_065458 and GW191109_010717 analyses)] The central claim that broad priors on the phase/amplitude uncertainty parameters absorb dominant systematics (including glitches) without biasing astrophysical parameters rests on real-event analyses without ground truth. No injection-recovery tests with simulated signals containing known parameters plus injected glitches or deglitching artifacts are reported; such tests are required to distinguish genuine error mitigation from posterior broadening or new degeneracies introduced by the additional free parameters.
- [Discussion of χ_p results and prior choices] The reported consistency in χ_p for GW200129_065458 (0.60^{+0.31}_{-0.33}, 0.58^{+0.30}_{-0.29}, 0.56^{+0.31}_{-0.28} across the three waveform models) could arise in part from the added model flexibility rather than independent error reduction. The manuscript should quantify the information loss or degeneracy introduced by the uncertainty parameters (e.g., via prior-posterior overlap or effective degrees of freedom) and compare against a baseline without uncertainty parameters.
minor comments (1)
- [Abstract and results] The abstract and results text use inconsistent notation for event names (e.g., GW191109_010717 with escaped underscore); standardize to standard LVK naming convention throughout.
Simulated Author's Rebuttal
We thank the referee for their constructive and insightful comments on our manuscript. We appreciate the positive assessment of the potential utility of our data-driven approach for mitigating systematics in gravitational-wave parameter estimation. We address each major comment below, agreeing where revisions are warranted and outlining specific changes to strengthen the validation of our results.
read point-by-point responses
-
Referee: [Results section (GW200129_065458 and GW191109_010717 analyses)] The central claim that broad priors on the phase/amplitude uncertainty parameters absorb dominant systematics (including glitches) without biasing astrophysical parameters rests on real-event analyses without ground truth. No injection-recovery tests with simulated signals containing known parameters plus injected glitches or deglitching artifacts are reported; such tests are required to distinguish genuine error mitigation from posterior broadening or new degeneracies introduced by the additional free parameters.
Authors: We agree that the absence of controlled injection-recovery tests represents a limitation in fully validating that the uncertainty models mitigate rather than merely broaden posteriors. Our focus has been on demonstrating practical consistency for real events flagged in the literature, where ground truth is unavailable by definition. In the revised manuscript we will add a dedicated subsection presenting injection studies: we will inject simulated binary black hole signals with known parameters into real O3 noise segments, superimpose synthetic glitches near the signal, apply deglitching procedures, and then recover the parameters both with and without the parametric uncertainty models. This will allow quantitative assessment of bias, posterior width changes, and recovery fidelity. revision: yes
-
Referee: [Discussion of χ_p results and prior choices] The reported consistency in χ_p for GW200129_065458 (0.60^{+0.31}_{-0.33}, 0.58^{+0.30}_{-0.29}, 0.56^{+0.31}_{-0.28} across the three waveform models) could arise in part from the added model flexibility rather than independent error reduction. The manuscript should quantify the information loss or degeneracy introduced by the uncertainty parameters (e.g., via prior-posterior overlap or effective degrees of freedom) and compare against a baseline without uncertainty parameters.
Authors: We thank the referee for highlighting the need for quantitative diagnostics. In the revised version we will augment the discussion with explicit metrics: we will report the overlap between prior and posterior distributions for the phase and amplitude uncertainty parameters, compute the Kullback-Leibler divergence to measure information gain, and estimate effective degrees of freedom. We will also include side-by-side posterior comparisons for χ_p and component spins obtained from the standard waveform models versus the augmented models, thereby demonstrating that the improved inter-model agreement is not solely attributable to extra flexibility. revision: yes
Circularity Check
Broad priors on self-cited uncertainty parameters produce consistency by added flexibility
specific steps
-
self citation load bearing
[Abstract]
"we reanalyze these events using a couple of parametric models developed in previous work that incorporate uncertainties in both the phase and amplitude of the GW waveform. In this data-driven approach, we apply sufficiently broad priors on the uncertainty parameters to account for potential systematic errors. Our findings show that the proposed method effectively reduces systematic errors, even those arising from data artifacts, such as glitches occurring near a signal and the deglitching process in GW frame files. Similarly, inconsistent results from different waveform models become much more"
The claim that the method 'effectively reduces systematic errors' and makes 'inconsistent results from different waveform models become much more consistent' rests on the parametric uncertainty models from previous work plus broad priors; the observed consistency (e.g., χ_p ≈ 0.58 across IMRPhenomXPHM, IMRPhenomXO4a, NRSur7dq4) is a direct statistical consequence of the extra degrees of freedom absorbing both waveform discrepancies and glitches, rather than an independent external validation.
full rationale
The paper's central result (consistent χ_p across waveform models for GW200129_065458 and similar for GW191109_010717) is obtained by reanalyzing real events with parametric phase/amplitude uncertainty models taken from prior work, equipped with broad priors. This added flexibility absorbs model differences and data artifacts by construction, yielding the reported consistency without ground-truth injection tests to separate genuine mitigation from degeneracy broadening. The derivation therefore depends on the self-cited framework for its load-bearing claim, producing moderate circularity even though the priors are not explicitly tuned to the target events.
Axiom & Free-Parameter Ledger
free parameters (1)
- phase and amplitude uncertainty parameters
axioms (1)
- domain assumption The parametric models from prior work adequately represent the dominant classes of waveform and data-analysis systematics present in the selected events.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
we model our ignorance by introducing amplitude and phase uncertainties in the waveform models... apply sufficiently broad priors on the uncertainty parameters
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
ΔΘ_total = ΔΘ_WF + ΔΘ_MP + ΔΘ_DA
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
GW150914 095045 Even though the GW150914 095045 runs did not show any signs of systematic errors across studies, we still use it as an example of how the framework behaves when no significant systematic errors are present. The mismatch study (Figures 2 and 3) indicate that this event shows the minimal intrinsic difference across different waveform models....
-
[2]
GW190412 065458 The event GW190412 065458 is one of the events with an asymmetric mass ratio and a positive effec- tive spin parameterχ eff [95]. In GWTC-2.1 analy- sis, the reported values of the effective spin parame- ter areχeff = 0.25+0.10 −0.10 (IMRPhenomXPHM) and 0.14+0.17 −0.07 (SEOBNRv4PHM). The 4OGC analysis recoveredχ eff = 0.25+0.10 −0.09 (IMRP...
-
[3]
The final mass of the remanent black hole is an intermediate mass black hole (IMBH) [99]
GW190521 065458 The event GW190521 065458 [98] is a high mass sig- nal where the significant probability distribution of the primary black hole mass lies in the mass gap produced by the pair-instability supernova process. The final mass of the remanent black hole is an intermediate mass black hole (IMBH) [99]. This event could have formed through various ...
-
[4]
GW191109 010717 The event GW191109 010717 is particularly significant because it exhibits anti-aligned spin properties, which may hint at a dynamical origin [101]. However, this event was affected by a scattered light glitch that occurred very close to the signal in the LIGO Livingston detec- tor. In the GWTC-3 analysis, both waveform models, SEOBNRv4PHMa...
-
[5]
These observations are consistent among baseline runs and RPE parametrization runs. In a previous study [29], the authors noted that the negative effective spin mea- surement is attributed to the segment of the LIGO Liv- ingston data between 30 and 40 Hz, or 0.1-0.004 seconds before the merger, which coincides with the excess glitch power. In figure 7, we...
-
[6]
GW200129 065458 GW200129 065458 is another interesting event that may have originated through a dynamical formation channel, especially considering its reported high pre- 13 35 40 45 50 55 60 65 0.0 0.1Deglitched Baseline Runs 35 40 45 50 55 60 65 0.0 0.1 WF-Error Runs 35 40 45 50 55 60 65 source 0.0 0.1Raw 35 40 45 50 55 60 65 source 0.0 0.1 1.0 0.5 0.0 ...
-
[7]
GW200208 222617 GW200208 222617 is a comparatively low SNR (net- work SNR∼8) event that shows multimodal mass pos- 15 0.1 0.0 0.1 0.2 rel H1 20.0 24.5 30.0 36.7 45.0 55.1 67.5 82.7 101.2 124.0 Frequency [Hz] 0.1 0.0 0.1 0.2 rel L1 GW191109 (Single detector PE runs) Deglitched Raw Prior FIG. 8. The figure shows the 1D marginalized posterior samples of the ...
work page 2000
-
[8]
A. G. Abacet al.(LIGO Scientific, VIRGO, KAGRA), (2025), arXiv:2508.18082 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2025
-
[9]
The LIGO Scientific Collaboration, The Virgo Collab- oration, The KAGRA Collaboration,et al., “GWTC-3: Compact Binary Coalescences Observed by LIGO and Virgo During the Second Part of the Third Observing Run,” (2021)
work page 2021
- [10]
- [11]
- [12]
-
[13]
A. G. Abacet al.(LIGO Scientific, VIRGO, KAGRA), (2025), arXiv:2508.18083 [astro-ph.HE]
work page internal anchor Pith review Pith/arXiv arXiv 2025
-
[14]
The population of merging compact binaries inferred using gravitational waves through GWTC-3
R. Abbottet al.(KAGRA, VIRGO, LIGO Scientific), Phys. Rev. X13, 011048 (2023), arXiv:2111.03634 [astro-ph.HE]
work page internal anchor Pith review Pith/arXiv arXiv 2023
-
[15]
A. G. Abacet al.(LIGO Scientific, Virgo, KAGRA), Phys. Rev. Lett.136, 041403 (2026), arXiv:2509.08099 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[16]
Tests of General Relativity with GWTC-3
R. Abbottet al.(LIGO Scientific, VIRGO, KAGRA), (2021), arXiv:2112.06861 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2021
-
[17]
A. G. Abacet al.(LIGO Scientific, VIRGO, KAGRA), (2026), arXiv:2603.19019 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[18]
A. G. Abacet al.(LIGO Scientific, VIRGO, KAGRA), (2026), arXiv:2603.19020 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[19]
A. G. Abacet al.(LIGO Scientific, VIRGO, KAGRA), (2026), arXiv:2603.19021 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[20]
B. P. Abbottet al.(LIGO Scientific, Virgo), Phys. Rev. Lett.121, 161101 (2018), arXiv:1805.11581 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2018
- [21]
-
[22]
R. Abbottet al.(LIGO Scientific, Virgo, KAGRA), Astrophys. J.949, 76 (2023), arXiv:2111.03604 [astro- 22 ph.CO]
-
[23]
A. G. Abacet al.(LIGO Scientific, VIRGO, KAGRA), (2025), arXiv:2509.04348 [astro-ph.CO]
work page internal anchor Pith review Pith/arXiv arXiv 2025
- [24]
-
[25]
Abbottet al.(LIGO Scientific, KAGRA, VIRGO), Astrophys
R. Abbottet al.(LIGO Scientific, KAGRA, VIRGO), Astrophys. J.970, 191 (2024), arXiv:2304.08393 [gr-qc]
-
[26]
A. G. Abacet al.(LIGO Scientific, VIRGO, KAGRA), Astrophys. J. Lett.993, L25 (2025), arXiv:2507.08219 [astro-ph.HE]
work page internal anchor Pith review Pith/arXiv arXiv 2025
- [27]
- [28]
-
[29]
A. Samajdar and T. Dietrich, Phys. Rev. D98, 124030 (2018), arXiv:1810.03936 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[30]
A. Samajdar and T. Dietrich, Phys. Rev. D100, 024046 (2019), arXiv:1905.03118 [gr-qc]
- [31]
- [32]
- [33]
-
[34]
S. Hourihane and K. Chatziioannou, Phys. Rev. D112, 084006 (2025), arXiv:2506.21869 [gr-qc]
- [35]
- [36]
-
[37]
S. Hourihane, K. Chatziioannou, M. Wijngaarden, D. Davis, T. Littenberg, and N. Cornish, Phys. Rev. D 106, 042006 (2022), arXiv:2205.13580 [gr-qc]
- [38]
- [39]
- [40]
-
[41]
C. Plunkett, S. Hourihane, and K. Chatziioannou, Phys. Rev. D106, 104021 (2022), arXiv:2208.02291 [gr- qc]
-
[42]
S. Bini, K. Kr´ ol, K. Chatziioannou, and M. Isi, (2026), arXiv:2601.09678 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2026
- [43]
- [44]
-
[45]
Spin-Precession: Breaking the Black Hole--Neutron Star Degeneracy
K. Chatziioannou, N. Cornish, A. Klein, and N. Yunes, Astrophys. J. Lett.798, L17 (2015), arXiv:1402.3581 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[46]
Can we measure individual black-hole spins from gravitational-wave observations?
M. P¨ urrer, M. Hannam, and F. Ohme, Phys. Rev. D 93, 084042 (2016), arXiv:1512.04955 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2016
-
[47]
J. Calder´ on Bustillo, S. Husa, A. M. Sintes, and M. P¨ urrer, Phys. Rev. D93, 084019 (2016), arXiv:1511.02060 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2016
-
[48]
First higher-multipole model of gravitational waves from spinning and coalescing black-hole binaries
L. London, S. Khan, E. Fauchon-Jones, C. Garc´ ıa, M. Hannam, S. Husa, X. Jim´ enez-Forteza, C. Kalaghatgi, F. Ohme, and F. Pannarale, Phys. Rev. Lett.120, 161102 (2018), arXiv:1708.00404 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[49]
R. Cotesta, A. Buonanno, A. Boh´ e, A. Taracchini, I. Hinder, and S. Ossokine, Phys. Rev. D98, 084028 (2018), arXiv:1803.10701 [gr-qc]
-
[50]
C. Kalaghatgi, M. Hannam, and V. Raymond, Phys. Rev. D101, 103004 (2020), arXiv:1909.10010 [gr-qc]
-
[51]
G. Pratten, P. Schmidt, R. Buscicchio, and L. M. Thomas, Phys. Rev. Res.2, 043096 (2020), arXiv:2006.16153 [gr-qc]
-
[52]
C. Garc´ ıa-Quir´ os, M. Colleoni, S. Husa, H. Estell´ es, G. Pratten, A. Ramos-Buades, M. Mateu-Lucena, and R. Jaume, Phys. Rev. D102, 064002 (2020), arXiv:2001.10914 [gr-qc]
- [53]
-
[54]
B. P. Abbottet al.(KAGRA, LIGO Scientific, Virgo), Living Rev. Rel.19, 1 (2016), arXiv:1304.0670 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2016
-
[55]
A. Buikemaet al.(aLIGO), Phys. Rev. D102, 062003 (2020), arXiv:2008.01301 [astro-ph.IM]
-
[56]
J. Aasiet al.(LIGO Scientific), Class. Quant. Grav.32, 074001 (2015), arXiv:1411.4547 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[57]
B. P. Abbottet al., Phys. Rev. D93, 112004 (2016), [Addendum: Phys.Rev.D 97, 059901 (2018)], arXiv:1604.00439 [astro-ph.IM]
work page internal anchor Pith review Pith/arXiv arXiv 2016
-
[58]
Model Waveform Accuracy Standards for Gravitational Wave Data Analysis
L. Lindblom, B. J. Owen, and D. A. Brown, Phys. Rev. D78, 124020 (2008), arXiv:0809.3844 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2008
-
[59]
M. P¨ urrer and C.-J. Haster, Phys. Rev. Res.2, 023151 (2020), arXiv:1912.10055 [gr-qc]
- [60]
- [61]
-
[62]
L. Pompili, A. Buonanno, and M. P¨ urrer, Phys. Rev. D112, 044062 (2025), arXiv:2410.16859 [gr-qc]
-
[63]
S. Mezzasoma, C.-J. Haster, and N. Yunes, (2026), arXiv:2603.15765 [gr-qc]
-
[64]
B. P. Abbottet al.(LIGO Scientific, Virgo), Class. Quant. Grav.34, 104002 (2017), arXiv:1611.07531 [gr- qc]
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[65]
R. Abbottet al.(KAGRA, VIRGO, LIGO Scientific), Phys. Rev. X13, 041039 (2023), arXiv:2111.03606 [gr- qc]
work page internal anchor Pith review Pith/arXiv arXiv 2023
-
[66]
G. Ashton and S. Khan, Phys. Rev. D101, 064037 (2020), arXiv:1910.09138 [gr-qc]
- [67]
-
[68]
G. Ashton and T. Dietrich, Nature Astron.6, 961 (2022), arXiv:2111.09214 [gr-qc]
- [69]
-
[70]
A. Puecher, A. Samajdar, G. Ashton, C. Van Den Broeck, and T. Dietrich, Phys. Rev. D109, 023019 (2024), arXiv:2310.03555 [gr-qc]
- [71]
- [72]
- [73]
-
[74]
R. Bachhar, M. P¨ urrer, and S. R. Green, (2024), arXiv:2410.17168 [gr-qc]
-
[75]
S. Mezzasoma, C.-J. Haster, C. B. Owen, N. J. Cor- nish, and N. Yunes, Phys. Rev. D112, 044057 (2025), arXiv:2503.23304 [gr-qc]
- [76]
- [77]
-
[78]
A. Ramos-Buades, S. Husa, G. Pratten, H. Estell´ es, C. Garc´ ıa-Quir´ os, M. Mateu-Lucena, M. Colleoni, and R. Jaume, Phys. Rev. D101, 083015 (2020), arXiv:1909.11011 [gr-qc]
-
[79]
D. A. Brown and P. J. Zimmerman, Phys. Rev. D81, 024007 (2010), arXiv:0909.0066 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2010
-
[80]
P. Narayan, N. K. Johnson-McDaniel, and A. Gupta, Phys. Rev. D108, 064003 (2023), arXiv:2306.04068 [gr- qc]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.