pith. sign in

arxiv: 2510.13677 · v2 · submitted 2025-10-15 · 🌀 gr-qc · astro-ph.IM· cs.NE· physics.comp-ph

APRIL: Auxiliary Physically-Redundant Information in Loss -- A physics-informed framework for parameter estimation with a gravitational-wave case study

Pith reviewed 2026-05-18 07:10 UTC · model grok-4.3

classification 🌀 gr-qc astro-ph.IMcs.NEphysics.comp-ph
keywords physics-informed neural networksgravitational wave parameter estimationloss function augmentationcompact binary coalescencechirp massmachine learning for physicsinspiral waveforms
0
0 comments X

The pith

Adding auxiliary physically-redundant terms to the loss preserves the physical minimum while reshaping the landscape for faster convergence in neural network parameter estimation.

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

The paper presents APRIL as a way to augment the standard supervised loss in neural networks with extra terms drawn from exact physical redundancy relations among the outputs. These additions are shown to leave the true physical solution unchanged as the global minimum but alter the surrounding loss surface so that gradient descent reaches physically consistent solutions more readily. In a gravitational-wave case study the method is applied to recover chirp mass, total mass and symmetric mass ratio from noise-free inspiral waveforms using a simple fully-connected network. The resulting test accuracy improves by up to an order of magnitude, particularly for the harder-to-learn parameters. The approach is positioned as a scalable complement to conventional physics-informed networks when many realizations of the same physics must be processed.

Core claim

By including auxiliary physically-redundant information in the loss, the training objective keeps the original physical minimum intact while changing the geometry of the loss surface, which demonstrably improves convergence and yields substantially higher accuracy when estimating the chirp mass, total mass and symmetric mass ratio of compact binary systems from simulated inspiral-frequency waveforms.

What carries the argument

APRIL augments the supervised output-target loss with auxiliary terms that exploit exact physical redundancy relations among the neural-network outputs, such as algebraic identities linking chirp mass, total mass and mass ratio.

If this is right

  • The method scales to large collections of systems that share the same underlying physics while still enforcing physical consistency.
  • Accuracy gains are largest for parameters that standard supervised training learns poorly.
  • The framework remains compatible with future extensions that incorporate realistic noise and broader parameter ranges.
  • It supplies a complementary route to standard PINNs when the task involves many independent realizations rather than a single system.

Where Pith is reading between the lines

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

  • The same redundancy-augmentation idea could be transferred to other inverse problems in which outputs are linked by known algebraic or differential identities.
  • If the network architecture is known to be under-expressive, the auxiliary terms may still regularize training but would require separate validation against injected signals.
  • Pairing APRIL with existing noise-robust training schedules might extend the observed accuracy gains into the realistic-data regime without changing the core construction.

Load-bearing premise

The auxiliary redundancy terms can be added without introducing bias or inconsistency even when the input signals contain realistic noise or when the neural network cannot perfectly represent the underlying waveform model.

What would settle it

A controlled experiment that adds realistic detector noise to the same inspiral waveforms and then checks whether the recovered parameters become systematically biased relative to the noise-free case would directly test whether the physical minimum remains unbiased.

Figures

Figures reproduced from arXiv: 2510.13677 by Francesco Di Clemente, Leigh Smith, Matteo Scialpi, Micha{\l} Bejger.

Figure 1
Figure 1. Figure 1: An example of f(tk) from the 1.5PN CBC GW event, corresponding to m1 ≃ 78.1 M⊙, m2 ≃ 12.6 M⊙ (M ≃ 25.4 M⊙, Mtot ≃ 90.7 M⊙, and η ≃ 0.12). 3.2 Methodology We will now describe the simulated datasets, the network used for the benchmark study and the benchmark itself. The dataset (Sec. 3.2.1) and the algorithm (Sec. 3.2.2) are purposefully simple in order to focus our study on the impact of the loss component… view at source ↗
Figure 2
Figure 2. Figure 2: Training, validation and test datasets for M, Mtot and η. Training and validation datasets are generated by sampling Mtot and η from a uniform distribution. The test dataset is obtained by sampling m1 and q from the mass distribution inferred by LVK collaboration from the GWTC-4 catalog [45, 46], as described in Sec. 3.2.1. Summarizing, an input dataset is composed of D frequency arrays {fk} K k=1 = f(tk),… view at source ↗
Figure 3
Figure 3. Figure 3: Algorithm training flow. The frequency array {fk} K k is given as input to the FCNN architecture, to give {M, Mtot, η}θ outputs. The outputs are then combined in different algebraic quantities that will be substituted in the different loss terms. The total loss is then computed as the sum of all terms. Its value and the gradient meta-data permits to update the NN parameter θ for a new epoch. During the tra… view at source ↗
Figure 4
Figure 4. Figure 4: RL1 results (median and 68% CI) for all the runs. Left panels show the study of the mass parameters individually, upper right panel shows the same study for the sum of the three, while bottom right panel shows the needed epochs to converge. The shape and the color of the markers determines the training dataset and the batch sizes. When APRIL losses are absent the RL1 result for the common sum is worse by a… view at source ↗
Figure 5
Figure 5. Figure 5: Test output relative errors on mass components comparing the runs for {D, B, seed} =  5 × 103 , 16, 1 [PITH_FULL_IMAGE:figures/full_fig_p012_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Ground truth loss components for the different runs. The combination of hyperparameters for these runs are {D, B, seed} =  5 × 103 , 16, 1 [PITH_FULL_IMAGE:figures/full_fig_p013_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: RL1 error comparing Lt + LAPRIL for the same batch size and a different data size. Here we set always the seed value to 1. We can see clearly a shift for M, which is anyway compensated with increasing D. 4.3 Dependence of parameter’s accuracy on the training data size Looking again at [PITH_FULL_IMAGE:figures/full_fig_p013_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: RL1 results (median and 68% CI) for all the runs. In the y axis the runs are labeled following their {αt, αdf, αp, αa} values. The different panels are dividing the study of the mass parameters individually, while the shape and the color of the markers determines the training dataset and the batch sizes. One can clearly notice that when APRIL losses are absent the RL1 result for the common sum is worse tha… view at source ↗
read the original abstract

Physics-Informed Neural Networks (PINNs) embed the partial differential equations (PDEs) governing the system under study directly into the training of Neural Networks, ensuring solutions that respect physical laws. While effective for single-system problems, standard PINNs scale poorly to datasets containing many realizations of the same underlying physics with varying parameters. To address this limitation, we present a complementary approach by including auxiliary physically-redundant information in loss (APRIL), i.e. augment the standard supervised output-target loss with auxiliary terms which exploit exact physical redundancy relations among outputs. We mathematically demonstrate that these terms preserve the true physical minimum while reshaping the loss landscape, improving convergence toward physically consistent solutions. As a proof-of-concept, we benchmark APRIL on a fully-connected neural network for gravitational wave (GW) parameter estimation (PE). We use simulated, noise-free compact binary coalescence (CBC) signals, focusing on inspiral-frequency waveforms to recover the chirp mass $\mathcal{M}$, the total mass $M_\mathrm{tot}$, and symmetric mass ratio $\eta$ of the binary. In this controlled setting, we show that APRIL achieves up to an order-of-magnitude improvement in test accuracy, especially for parameters that are otherwise difficult to learn. This method provides physically consistent learning for large multi-system datasets and is well suited for future GW analyses involving realistic noise and broader parameter ranges.

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 / 2 minor

Summary. The paper introduces APRIL, a framework that augments the standard supervised loss in neural networks with auxiliary terms enforcing exact physical redundancy relations among output parameters (e.g., among chirp mass ℳ, total mass M_tot, and symmetric mass ratio η). It claims a mathematical demonstration that these terms preserve the true physical minimum while reshaping the loss landscape to improve convergence, and reports up to an order-of-magnitude accuracy gain in a proof-of-concept gravitational-wave parameter estimation task using a fully-connected network on noise-free simulated inspiral waveforms.

Significance. If the mathematical preservation result holds beyond the noise-free exact-match case and the accuracy gains are reproducible with baselines and error bars, APRIL could provide a lightweight, scalable complement to standard PINNs for enforcing consistency across large multi-system datasets in gravitational-wave astronomy and similar parameter-estimation problems.

major comments (2)
  1. [Mathematical demonstration] Mathematical demonstration (likely §3 or equivalent): the claim that auxiliary redundancy terms 'preserve the true physical minimum' is shown only under the assumption that network outputs satisfy the relations exactly (true by construction for noise-free simulated waveforms); the manuscript provides no explicit derivation steps or analysis of how the weighted auxiliary terms shift the effective minimum when additive noise or imperfect network representation makes the relations inconsistent with the data.
  2. [Results section] Results on accuracy improvement (likely §4 or Table/Figure reporting test accuracy): the order-of-magnitude gain is stated without baseline comparisons to the unaugmented loss, without error bars or statistical significance on the metrics, and without details on test-set size or cross-validation, undermining verification of the central empirical claim.
minor comments (2)
  1. [Introduction] Notation: the abstract and introduction use both ℳ and M_tot without an early explicit statement of the exact redundancy relation (e.g., ℳ = (M_tot η)^{3/5} (1-η)^{1/5} or equivalent) that the auxiliary terms enforce.
  2. [Discussion] The manuscript states the method is 'well suited for future GW analyses involving realistic noise' but contains no discussion or preliminary test of how the auxiliary weighting should be chosen or annealed when noise is present.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive and detailed report. We address each major comment below and describe the revisions we will make to strengthen the manuscript.

read point-by-point responses
  1. Referee: Mathematical demonstration (likely §3 or equivalent): the claim that auxiliary redundancy terms 'preserve the true physical minimum' is shown only under the assumption that network outputs satisfy the relations exactly (true by construction for noise-free simulated waveforms); the manuscript provides no explicit derivation steps or analysis of how the weighted auxiliary terms shift the effective minimum when additive noise or imperfect network representation makes the relations inconsistent with the data.

    Authors: We agree that the current mathematical demonstration is developed explicitly for the exact-match case that holds by construction in our noise-free proof-of-concept. The derivation establishes that the auxiliary terms evaluate to zero at the ground-truth parameters (where the redundancy relations are satisfied) and are strictly non-negative elsewhere, thereby leaving the location of the true minimum unchanged. We did not supply a full perturbation analysis for small inconsistencies arising from noise or network approximation error. In the revision we will add explicit derivation steps in Section 3 together with a short analytical example showing that, for small deviations, the auxiliary terms act as a restoring force toward consistency without introducing new spurious minima near the physical solution. revision: yes

  2. Referee: Results on accuracy improvement (likely §4 or Table/Figure reporting test accuracy): the order-of-magnitude gain is stated without baseline comparisons to the unaugmented loss, without error bars or statistical significance on the metrics, and without details on test-set size or cross-validation, undermining verification of the central empirical claim.

    Authors: We acknowledge that the results section would benefit from more complete reporting. While the improvement is measured relative to standard supervised training, we will revise the section to include an explicit side-by-side comparison table, report means and standard deviations over five independent runs with different random seeds, state the test-set size (1000 waveforms), and describe the single hold-out evaluation protocol. We will also add a brief note explaining why cross-validation was not performed in this controlled proof-of-concept study. revision: yes

Circularity Check

0 steps flagged

No significant circularity in APRIL derivation or claims

full rationale

The paper introduces auxiliary loss terms exploiting exact physical redundancies among outputs (chirp mass, total mass, symmetric mass ratio) and states a mathematical demonstration that these terms preserve the true physical minimum while reshaping the landscape. This demonstration relies on the algebraic properties of the added terms being zero at the true parameter values satisfying the redundancies, which is an independent property of the loss construction rather than a reduction to fitted inputs or self-citation. The GW parameter estimation results are presented as a controlled proof-of-concept on noise-free simulated waveforms, with reported accuracy gains treated as empirical observations rather than quantities forced by the equations themselves. No load-bearing step reduces by construction to the inputs; the framework remains self-contained against external benchmarks of loss augmentation.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the existence of exact, a-priori known physical redundancy relations among output parameters that can be turned into auxiliary loss terms without shifting the true minimum.

axioms (1)
  • domain assumption Exact physical redundancy relations exist among the network outputs (e.g., algebraic relations linking chirp mass, total mass, and symmetric mass ratio) and can be expressed as auxiliary loss terms.
    The method explicitly relies on these relations to augment the loss while preserving the physical minimum.

pith-pipeline@v0.9.0 · 5802 in / 1303 out tokens · 44936 ms · 2026-05-18T07:10:44.904401+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

59 extracted references · 59 canonical work pages · 1 internal anchor

  1. [1]

    Machine learning:a review,

    C. Isonkobong, “Machine learning:a review,”Semiconductor Science and Information Devices, vol. 2, 11 2020

  2. [2]

    Enhancing gravitational-wave science with machine learning,

    E. Cuoco, J. Powell, M. Cavaglià, K. Ackley, M. Bejger, C. Chatterjee, M. Coughlin, S. Cough- lin, P. Easter, R. Essick, H. Gabbard, T. Gebhard, S. Ghosh, L. Haegel, A. Iess, D. Keitel, Z. Márka, S. Márka, F. Morawski, T. Nguyen, R. Ormiston, M. Pürrer, M. Razzano, K. Staats, G. Vajente, and D. Williams, “Enhancing gravitational-wave science with machine ...

  3. [3]

    Applications of machine learning in gravitational-wave research with current interferometric detectors,

    E. Cuoco, M. Cavaglià, I. S. Heng, D. Keitel, and C. Messenger, “Applications of machine learning in gravitational-wave research with current interferometric detectors,”Living Reviews in Relativity, vol. 28, Feb. 2025

  4. [4]

    Cuoco, ed.,Gravitational Wave Science with Machine Learning, vol

    E. Cuoco, ed.,Gravitational Wave Science with Machine Learning, vol. 1 ofSpringer Series in Astrophysics and Cosmology. Springer Singapore, 1 ed., 2025. 4 b/w illustrations, 104 colour illustrations

  5. [5]

    Interpretable machine learning in physics: A review,

    S. J. Wetzel, S. Ha, R. Iten, M. Klopotek, and Z. Liu, “Interpretable machine learning in physics: A review,” 2025

  6. [6]

    A systematic literature review on the use of machine learning in software engineering,

    N. Fred and I. O. Temkin, “A systematic literature review on the use of machine learning in software engineering,” 2024

  7. [7]

    Machine learning applications in structural engineering - a review,

    P. Haneena Jasmine and S. Arun, “Machine learning applications in structural engineering - a review,”IOP Conference Series: Materials Science and Engineering, vol. 1114, p. 012012, mar 2021

  8. [8]

    A review of multimodal explainable artificial intelligence: Past, present and future,

    S. Sun, W. An, F. Tian, F. Nan, Q. Liu, J. Liu, N. Shah, and P. Chen, “A review of multimodal explainable artificial intelligence: Past, present and future,” 2024

  9. [9]

    Focus on explainable machine learning in sciences

    “Focus on explainable machine learning in sciences.”https://iopscience.iop.org/ collections/mlst-240319-506, 2024. IOPscience Collection. 20

  10. [10]

    Physics-informed neural networks: A deep learn- ing framework for solving forward and inverse problems involving nonlinear partial differential equations,

    M. Raissi, P. Perdikaris, and G. Karniadakis, “Physics-informed neural networks: A deep learn- ing framework for solving forward and inverse problems involving nonlinear partial differential equations,”Journal of Computational Physics, vol. 378, pp. 686–707, 2019

  11. [11]

    Scientific machine learning through physics-informed neural networks: Where we are and what’s next,

    S. Cuomo, V. S. di Cola, F. Giampaolo, G. Rozza, M. Raissi, and F. Piccialli, “Scientific machine learning through physics-informed neural networks: Where we are and what’s next,” 2022

  12. [12]

    Adaptive physics-informed neural networks: A survey,

    E. Torres, J. Schiefer, and M. Niepert, “Adaptive physics-informed neural networks: A survey,” 2025

  13. [13]

    Explainable autoencoder for neutron star dense matter parameter estimation,

    F. Di Clemente, M. Scialpi, and M. Bejger, “Explainable autoencoder for neutron star dense matter parameter estimation,” 2025

  14. [14]

    Maggiore,Gravitational Waves: Volume 1: Theory and Experiments

    M. Maggiore,Gravitational Waves: Volume 1: Theory and Experiments. Oxford University Press, 10 2007

  15. [15]

    Gravitational waves from merging compact binaries: How accurately can one extract the binary’s parameters from the inspiral waveform?,

    C. Cutler and E. E. Flanagan, “Gravitational waves from merging compact binaries: How accurately can one extract the binary’s parameters from the inspiral waveform?,”Physical Review D, vol. 49, p. 2658–2697, Mar. 1994

  16. [16]

    Learning orbital dynamics of binary black hole systems from gravitational wave measurements,

    B. Keith, A. Khadse, and S. E. Field, “Learning orbital dynamics of binary black hole systems from gravitational wave measurements,”Physical Review Research, vol. 3, Nov. 2021

  17. [17]

    Using physics-informed neural networks to compute quasinormal modes,

    A. S. Cornell, A. Ncube, and G. Harmsen, “Using physics-informed neural networks to compute quasinormal modes,”Phys. Rev. D, vol. 106, p. 124047, Dec 2022

  18. [18]

    Calculating quasi-normal modes of schwarzschild black holes with physics informed neural networks,

    N. Patel, A. Aykutalp, and P. Laguna, “Calculating quasi-normal modes of schwarzschild black holes with physics informed neural networks,” 2024

  19. [19]

    Solv- ing the teukolsky equation with physics-informed neural networks,

    R. Luna, J. Calderón Bustillo, J. J. Seoane Martínez, A. Torres-Forné, and J. A. Font, “Solv- ing the teukolsky equation with physics-informed neural networks,”Phys. Rev. D, vol. 107, p. 064025, Mar 2023

  20. [20]

    Machine learning for conservative-to- primitive in relativistic hydrodynamics,

    T. Dieselhorst, W. Cook, S. Bernuzzi, and D. Radice, “Machine learning for conservative-to- primitive in relativistic hydrodynamics,”Symmetry, vol. 13, no. 11, 2021

  21. [21]

    Magnetohydrodynamics with physics informed neural op- erators,

    S. G. Rosofsky and E. A. Huerta, “Magnetohydrodynamics with physics informed neural op- erators,”Machine Learning: Science and Technology, vol. 4, p. 035002, jul 2023

  22. [22]

    Grinn: a physics-informed neural network for solving hydrodynamic systems in the presence of self-gravity,

    S. Auddy, R. Dey, N. J. Turner, and S. Basu, “Grinn: a physics-informed neural network for solving hydrodynamic systems in the presence of self-gravity,”Machine Learning: Science and Technology, vol. 5, p. 025014, apr 2024

  23. [23]

    Modellingforce-freeneutronstarmagne- tospheres using physics-informed neural networks,

    J.F.Urbán, P.Stefanou, C.Dehman, andJ.A.Pons, “Modellingforce-freeneutronstarmagne- tospheres using physics-informed neural networks,”Monthly Notices of the Royal Astronomical Society, vol. 524, pp. 32–42, 06 2023

  24. [24]

    Solving the pulsar equation using physics-informed neural networks,

    P. Stefanou, J. F. Urbán, and J. A. Pons, “Solving the pulsar equation using physics-informed neural networks,”Monthly Notices of the Royal Astronomical Society, vol. 526, pp. 1504–1511, 09 2023

  25. [25]

    Solving einstein equations using deep learning,

    Z.-H. Li, C.-Q. Li, and L.-G. Pang, “Solving einstein equations using deep learning,” 2023

  26. [26]

    April implementation and benchmark results

    M. Scialpi, F. Di Clemente, L. Smith, and M. Bejger, “April implementation and benchmark results.”https://github.com/matteoscialpi/APRIL.git, 2025. accessed 2025-09-08

  27. [27]

    Multilayer feedforward networks are universal approximators,

    K. Hornik, M. Stinchcombe, and H. White, “Multilayer feedforward networks are universal approximators,”Neural Networks, vol. 2, no. 5, pp. 359–366, 1989

  28. [28]

    Gradient-based learning applied to document recognition,

    Y. Lecun, L. Bottou, Y. Bengio, and P. Haffner, “Gradient-based learning applied to document recognition,”Proceedings of the IEEE, vol. 86, no. 11, pp. 2278–2324, 1998

  29. [29]

    Learning representations by back- propagating errors,

    D. E. Rumelhart, G. E. Hinton, and R. J. Williams, “Learning representations by back- propagating errors,”nature, vol. 323, no. 6088, pp. 533–536, 1986

  30. [30]

    Variational inference with normalizing flows,

    D. J. Rezende and S. Mohamed, “Variational inference with normalizing flows,” 2016

  31. [31]

    Attention is all you need,

    A. Vaswani, N. Shazeer, N. Parmar, J. Uszkoreit, L. Jones, A. N. Gomez, L. Kaiser, and I. Polosukhin, “Attention is all you need,” 2023

  32. [32]

    R. M. Neal,Bayesian Learning for Neural Networks, vol. 118 ofLecture Notes in Statistics. Springer New York, NY, 1 ed., 1996

  33. [33]

    Die feldgleichungen der gravitation,

    A. Einstein, “Die feldgleichungen der gravitation,”Sitzungsberichte der Königlich Preussischen Akademie der Wissenschaften (Berlin), pp. 844–847, 1915. Presented on 25 November 1915

  34. [34]

    Näherungsweise integration der feldgleichungen der gravitation,

    A. Einstein, “Näherungsweise integration der feldgleichungen der gravitation,”Sitzungsberichte der Königlich Preussischen Akademie der Wissenschaften (Berlin), pp. 688–696, 1916. First derivation of gravitational waves from linearized field equations. 21

  35. [35]

    Über gravitationswellen,

    A. Einstein, “Über gravitationswellen,”Sitzungsberichte der Königlich Preussischen Akademie der Wissenschaften (Berlin), pp. 154–167, 1918. Corrected derivation of gravitational wave formula

  36. [36]

    Advanced LIGO,

    J. Aasi, B. P. Abbott, R. Abbott,et al., “Advanced LIGO,”Classical and Quantum Gravity, vol. 32, no. 7, p. 074001, 2015

  37. [37]

    Advanced Virgo: a second-generation interfer- ometric gravitational wave detector,

    F. Acernese, M. Agathos, K. Agatsuma,et al., “Advanced Virgo: a second-generation interfer- ometric gravitational wave detector,”Classical and Quantum Gravity, vol. 32, no. 2, p. 024001, 2015

  38. [38]

    Overview of KAGRA: Calibration, detector charac- terization, physical environmental monitors, and the geophysics interferometer,

    T. Akutsu, M. Ando, K. Arai,et al., “Overview of KAGRA: Calibration, detector charac- terization, physical environmental monitors, and the geophysics interferometer,”Progress of Theoretical and Experimental Physics, vol. 2021, no. 5, p. 05A101, 2021

  39. [39]

    Status of GEO-600,

    K. L. Dooley, “Status of GEO-600,”Journal of Physics: Conference Series, vol. 610, p. 012015, May 2015

  40. [40]

    On the gravitational field of a mass point according to Einstein's theory

    K. Schwarzschild, “Über das gravitationsfeld eines massenpunktes nach der einsteinschen the- orie,”Sitzungsberichte der Königlich Preussischen Akademie der Wissenschaften (Berlin), pp. 189–196, 1916. English translation available at arXiv:physics/9905030

  41. [41]

    Gravitational radiation from post-newtonian sources and inspiralling compact binaries,

    L. Blanchet, “Gravitational radiation from post-newtonian sources and inspiralling compact binaries,”Living Reviews in Relativity, vol. 5, Apr. 2002

  42. [42]

    The basic physics of the binary black hole merger GW150914,

    B. P. Abbottet al., “The basic physics of the binary black hole merger GW150914,”Annalen der Physik, vol. 529, Oct. 2016

  43. [43]

    Ueber die numerische auflösung von differentialgleichungen.,

    C. Runge, “Ueber die numerische auflösung von differentialgleichungen.,”Mathematische An- nalen, vol. 46, pp. 167–178, 1895

  44. [44]

    Beitrag zur nähenmgsweisen integration totaler differentialgleichungen.,

    V. W. Kütta, “Beitrag zur nähenmgsweisen integration totaler differentialgleichungen.,” Zeitschrift für mathematik und physik, vol. 46, pp. 435–453, 1901

  45. [45]

    GWTC-4.0: Updating the gravitational-wave transient catalog with observations from the first part of the fourth LIGO-Virgo-KAGRA observing run,

    LIGO Scientific Collaboration, Virgo Collaboration, and KAGRA Collaboration, “GWTC-4.0: Updating the gravitational-wave transient catalog with observations from the first part of the fourth LIGO-Virgo-KAGRA observing run,” 2025

  46. [46]

    GWTC-4.0: Population properties of merging compact binaries,

    LIGO Scientific Collaboration, Virgo Collaboration, and KAGRA Collaboration, “GWTC-4.0: Population properties of merging compact binaries,” 2025

  47. [47]

    Various techniques used in connection with random digits,

    J. von Neumann, “Various techniques used in connection with random digits,” inMonte Carlo Method, vol. 12 ofApplied Mathematics Series, pp. 36–38, Washington, DC: National Bureau of Standards, 1951. Original description of acceptance–rejection sampling

  48. [48]

    C. P. Robert and G. Casella,Monte Carlo Statistical Methods. New York, NY: Springer, 2 ed., 2004

  49. [49]

    Noticesurlaloiquelapopulationsuitdanssonaccroissement,

    P.-F.Verhulst, “Noticesurlaloiquelapopulationsuitdanssonaccroissement,”Correspondance mathématique et physique, vol. 10, pp. 113–121, 1838. Original introduction of the logistic equation

  50. [50]

    Adam: A method for stochastic optimization,

    D. P. Kingma and J. Ba, “Adam: A method for stochastic optimization,”International Con- ference on Learning Representations (ICLR), 2015

  51. [51]

    Decoupled weight decay regularization,

    I. Loshchilov and F. Hutter, “Decoupled weight decay regularization,”International Conference on Learning Representations (ICLR), 2019

  52. [52]

    Stochastic gradient descent tricks,

    L. Bottou, “Stochastic gradient descent tricks,” inNeural Networks: Tricks of the Trade (G. Montavon, G. B. Orr, and K.-R. Müller, eds.), pp. 421–436, Berlin, Heidelberg: Springer, 2012

  53. [53]

    torch.optim.lr_scheduler.reducelronplateau

    PyTorch Contributors, “torch.optim.lr_scheduler.reducelronplateau.”https://pytorch.org/ docs/stable/generated/torch.optim.lr_scheduler.ReduceLROnPlateau.html, 2025. Ac- cessed: 2025-08-13

  54. [54]

    torch.use_deterministic_algorithms

    PyTorch Contributors, “torch.use_deterministic_algorithms.”https://docs.pytorch.org/ docs/stable/generated/torch.use_deterministic_algorithms.html, 2025. Accessed: 2025-09-08

  55. [55]

    Nvidia GeForce RTX 4060 graphics card

    NVIDIA Corporation, “Nvidia GeForce RTX 4060 graphics card.”https://www.nvidia.com/ en-us/geforce/graphics-cards/40-series/rtx-4060/, 2023. Accessed: 2025-08-13

  56. [56]

    Nvidia GeForce RTX 4070 graphics card

    NVIDIA Corporation, “Nvidia GeForce RTX 4070 graphics card.”https://www.nvidia.com/ en-us/geforce/graphics-cards/40-series/rtx-4070/, 2023. Accessed: 2025-08-13

  57. [57]

    GWTC-4.0: Population properties of merging compact binaries

    LIGO Scientific Collaboration, Virgo Collaboration, and KAGRA Collaboration, “GWTC-4.0: Population properties of merging compact binaries.” Zenodo, Aug. 26 2025. 22

  58. [58]

    analy- ses_bbh.tar

    LIGO Scientific Collaboration, Virgo Collaboration, and KAGRA Collaboration, “analy- ses_bbh.tar.” Zenodo, Aug. 26 2025

  59. [59]

    Bbh- massspinredshift_brokenpowerlawtwopeaks_gaussiancomponentspins_powerlawredshift.h5

    LIGO Scientific Collaboration, Virgo Collaboration, and KAGRA Collaboration, “Bbh- massspinredshift_brokenpowerlawtwopeaks_gaussiancomponentspins_powerlawredshift.h5.” Zenodo, Aug. 26 2025. 23