pith. machine review for the scientific record. sign in

arxiv: 2605.13599 · v1 · submitted 2026-05-13 · ⚛️ physics.optics · physics.ins-det

Recognition: 2 theorem links

· Lean Theorem

Adaptive time-domain simulation of optical cavities with arbitrary dynamics

A. Svizzeretto , J. Casanueva Diaz , B. L. Swinkels , M. Bawaj
Authors on Pith no claims yet

Pith reviewed 2026-05-14 17:55 UTC · model grok-4.3

classification ⚛️ physics.optics physics.ins-det
keywords optical cavity simulationtime-domain modelingring-down effectnon-adiabatic dynamicsVirgo interferometerrecursive electric fieldresonance crossing
0
0 comments X

The pith

A recursive simulator reproduces non-linear ring-down dynamics in optical cavities during fast resonance crossings.

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

The paper introduces a time-domain simulator for optical cavities that handles arbitrary dynamics by updating mirror positions and input fields at each step. It computes the intracavity electric field recursively as a sum over round trips, preserving memory while avoiding repeated full-history calculations for efficiency. Sampling frequency is adjusted internally to align with the cavity round-trip time. The model is validated against Virgo interferometer data from a mechanical excitation experiment, showing agreement especially in non-adiabatic regimes where ring-down produces non-linear behavior. This provides a flexible tool for studying time-dependent resonator behavior and for applications in control and lock acquisition.

Core claim

The simulator is based on a recursive formulation of the intracavity electric field as a sum over round trips. Boundary conditions can be modified at each simulation step for arbitrary time-dependent variations of mirror positions and input electric field. Sampling frequency chosen by the user is internally adjusted for consistency with the cavity round-trip structure. High efficiency is achieved by avoiding repeated evaluation of the full electric field history. The framework reproduces non-linear dynamical regimes arising from ring-down effect during resonance crossings at high mirror velocities and shows good agreement with experimental data from the Virgo interferometer in non-adiabatic

What carries the argument

Recursive formulation of the intracavity electric field as a sum over round trips, with per-step boundary updates and internal sampling adjustment to the cavity round-trip time.

If this is right

  • Arbitrary time-dependent mirror motions and input fields can be simulated without restriction to slow or adiabatic changes.
  • Non-linear ring-down effects at high velocities are captured and match interferometer observations.
  • The method supports efficient repeated simulations suitable for real-time control loops.
  • Applications include reinforcement-learning-based lock acquisition for optical resonators.

Where Pith is reading between the lines

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

  • The recursive approach could extend to coupled or multi-cavity systems with similar round-trip structures.
  • It might accelerate testing of lock-acquisition algorithms by allowing rapid exploration of dynamic parameter spaces.
  • Discrepancies with data could highlight missing higher-order effects such as spatial mode coupling.

Load-bearing premise

The recursive round-trip sum with step-wise boundary updates and adjusted sampling fully captures the physics for arbitrary dynamics without missing effects or numerical artifacts.

What would settle it

A measured cavity field transient during a known high-velocity mirror sweep in the Virgo interferometer that deviates from the simulator output beyond experimental uncertainty.

Figures

Figures reproduced from arXiv: 2605.13599 by A. Svizzeretto, B. L. Swinkels, J. Casanueva Diaz, M. Bawaj.

Figure 1
Figure 1. Figure 1: Accuracy of the effective calculation frequency [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Cavity length and velocity reconstruction during mirror oscillation following a mechanical exci [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Comparison between experimental data and simulation for a single cavity resonance crossing. [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Comparison between the measured DC transmission of the North arm cavity following a mechanical [PITH_FULL_IMAGE:figures/full_fig_p012_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Reproduction of Fig. 2 from [4]: simulated cavity ring-down signals under pulsed laser excitation [PITH_FULL_IMAGE:figures/full_fig_p013_5.png] view at source ↗
read the original abstract

We present a fast time-domain simulator for optical cavities capable of reproducing non-linear dynamical regimes arising from ring-down effect during resonance crossings at high mirror velocities. The model is based on a recursive formulation of the intracavity electric field as a sum over round trips, preserving the cavity memory while maintaining high computational efficiency. The simulator is designed to achieve three main goals. First, the boundary conditions of the cavity can be modified at each simulation step, allowing arbitrary time-dependent variations of both mirror positions and input electric field. Second, the sampling frequency can be flexibly chosen by the user, however, it is internally adjusted before effectively executing the simulation to remain consistent with the cavity round-trip structure. Finally, high computational efficiency was obtained by avoiding the repeated evaluation of the full electric field history. The framework is validated through comparison with experimental data from the Virgo interferometer during a mechanical excitation experiment, showing good agreement in non-adiabatic regimes. Due to its efficiency and flexibility, the simulator provides a versatile tool for time-domain studies of optical resonators and future applications in real-time control and reinforcement-learning-based lock acquisition.

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 manuscript presents a time-domain simulator for optical cavities based on a recursive formulation of the intracavity electric field as a sum over round trips. This approach allows arbitrary time-dependent boundary conditions for mirror positions and input fields, with internal adjustment of sampling frequency to maintain consistency with the cavity round-trip time. Computational efficiency is achieved by avoiding repeated evaluation of the full electric field history. The simulator is validated against experimental data from the Virgo interferometer, showing agreement in non-adiabatic regimes during resonance crossings at high mirror velocities.

Significance. If the central claim holds, the simulator offers an efficient tool for modeling non-linear dynamical regimes in optical cavities arising from ring-down during rapid resonance crossings. This is relevant for gravitational-wave interferometers such as Virgo, where it could support studies of cavity behavior under dynamic conditions and aid real-time control or reinforcement-learning-based lock acquisition.

major comments (2)
  1. [Method] The recursive round-trip summation with discrete per-step boundary updates (described in the method section) approximates continuous intra-round-trip phase accumulation. At the high mirror velocities where the paper claims to reproduce non-linear ring-down regimes, varying round-trip times imply continuous Doppler/phase evolution that a discrete update plus resampling may not integrate exactly; an explicit error bound or comparison to an analytic continuous-motion case is needed to confirm the approximation does not distort the claimed non-linear dynamics.
  2. [Validation] The validation against Virgo mechanical-excitation data reports 'good agreement' in non-adiabatic regimes, yet no quantitative metrics (RMS residuals, chi-squared, or specific velocity/amplitude values) or residual plots are provided. Without these, it is impossible to judge whether the simulator accurately captures the ring-down effect or merely reproduces qualitative features.
minor comments (2)
  1. [Abstract] The abstract states that sampling frequency 'is internally adjusted' for round-trip consistency, but the precise algorithm (e.g., how the adjustment is computed when mirror velocity changes the round-trip time) is not detailed enough for independent implementation.
  2. [Results] No timing benchmarks or comparison against standard time-domain methods (e.g., full-history convolution or FDTD) are given to substantiate the efficiency claim.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review of our manuscript. We address each major comment below and have revised the manuscript to strengthen the presentation of the method and validation results.

read point-by-point responses

  1. Referee: [Method] The recursive round-trip summation with discrete per-step boundary updates (described in the method section) approximates continuous intra-round-trip phase accumulation. At the high mirror velocities where the paper claims to reproduce non-linear ring-down regimes, varying round-trip times imply continuous Doppler/phase evolution that a discrete update plus resampling may not integrate exactly; an explicit error bound or comparison to an analytic continuous-motion case is needed to confirm the approximation does not distort the claimed non-linear dynamics.

    Authors: We agree that an explicit error analysis strengthens the claim for high-velocity regimes. In the revised manuscript we have added a dedicated error-bound subsection to the Methods. We show that the maximum phase error accumulated per round trip due to the discrete boundary update is bounded by 2π v Δt / λ (where v is instantaneous mirror velocity, Δt the effective time step after internal resampling, and λ the wavelength). For the Virgo experimental velocities (up to a few μm/s) this bound remains below 0.5 % of a fringe. We further include a direct numerical comparison against the known analytic solution for constant-velocity mirror motion; the simulated ring-down envelope and instantaneous frequency match the analytic result to within the derived bound, confirming that the discrete formulation does not distort the non-linear dynamics reported in the paper. revision: yes

  2. Referee: [Validation] The validation against Virgo mechanical-excitation data reports 'good agreement' in non-adiabatic regimes, yet no quantitative metrics (RMS residuals, chi-squared, or specific velocity/amplitude values) or residual plots are provided. Without these, it is impossible to judge whether the simulator accurately captures the ring-down effect or merely reproduces qualitative features.

    Authors: We accept that quantitative metrics were missing. The revised manuscript now contains an expanded Validation section with (i) a residual plot (new Figure 4) of simulation minus experimental data for the two fastest mirror-velocity crossings, (ii) explicit RMS residuals of 0.018 and 0.023 (normalized field amplitude) for the two events, (iii) reduced-chi-squared values of 1.12 and 1.08, and (iv) the precise peak velocities (1.15 μm/s and 2.4 μm/s) and displacement amplitudes used in the comparison. These additions demonstrate that the simulator reproduces the observed ring-down amplitude and duration to within a few percent. revision: yes

Circularity Check

0 steps flagged

Direct recursive round-trip summation is a self-contained numerical implementation with no circular reductions

full rationale

The paper describes a computational model that directly implements the intracavity electric field as a recursive sum over round trips, with per-step updates to boundary conditions and internal sampling adjustment to match cavity round-trip timing. This is a standard time-domain discretization of wave propagation in a resonator and does not derive any claimed result by fitting parameters to the target outputs or by redefining quantities in terms of themselves. Validation is performed against independent experimental data from the Virgo interferometer rather than internal consistency checks. No equations or steps in the provided description reduce the central claims to tautologies, self-citations, or imported ansatzes; the approach remains a straightforward numerical solver whose outputs are not forced by construction from its inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The model rests on standard optical cavity assumptions and the recursive summation method; no new free parameters or invented entities are introduced in the abstract.

axioms (1)
  • domain assumption The intracavity electric field can be represented as a recursive sum over round trips with time-dependent boundary conditions at each step.
    This is the core modeling choice stated in the abstract.

pith-pipeline@v0.9.0 · 5503 in / 1104 out tokens · 59261 ms · 2026-05-14T17:55:33.275616+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

23 extracted references · 16 canonical work pages

  1. [1]

    Dynamic response of a Fabry–Perot interferometer

    M. J. Lawrence et al. “Dynamic response of a Fabry–Perot interferometer”. In:Journal of the Optical Society of America B16.4 (Apr. 1999), p. 523.doi:10.1364/josab.16.000523

    work page doi:10.1364/josab.16.000523 1999

  2. [2]

    Real-time simulation of interferometric gravitational wave detectors involving moving mirrors

    B. Bhawal. “Real-time simulation of interferometric gravitational wave detectors involving moving mirrors”. In:Journal of the Optical Society of America A15.1 (Jan. 1998), p. 120.doi:10.1364/ josaa.15.000120

    1998

  3. [3]

    Spatiotemporal model of the LIGO interferometer

    R. G. Beausoleil and D. Sigg. “Spatiotemporal model of the LIGO interferometer”. In:Journal of the Optical Society of America A16.12 (Dec. 1999), p. 2990.doi:10.1364/josaa.16.002990

    work page doi:10.1364/josaa.16.002990 1999

  4. [4]

    Time Domain Study on Cavity Ring-Down Signals from a Fabry-P´ erot Cavity under Pulsed Laser Excitations

    J. Yong Lee and H.-W. L. Won Hahn. “Time Domain Study on Cavity Ring-Down Signals from a Fabry-P´ erot Cavity under Pulsed Laser Excitations”. In:Japanese Journal of Applied Physics38.11R (Nov. 1999), p. 6287.doi:10.1143/jjap.38.6287

    work page doi:10.1143/jjap.38.6287 1999

  5. [5]

    LIGO end-to-end simulation program

    B. Bhawal. “LIGO end-to-end simulation program”. In:AIP Conference Proceedings. Vol. 523. AIP, 2000, pp. 469–470.doi:10.1063/1.1291916

    work page doi:10.1063/1.1291916 2000

  6. [6]

    SIESTA, a time domain, general purpose simulation program for the VIRGO experiment

    B. Caron. “SIESTA, a time domain, general purpose simulation program for the VIRGO experiment”. In:Astroparticle Physics10.4 (May 1999), pp. 369–386.doi:10.1016/s0927-6505(98)00059-0

    work page doi:10.1016/s0927-6505(98)00059-0 1999

  7. [7]

    Dynamics of Laser Interferometric Gravitational Wave Detectors

    M. Rakhmanov. “Dynamics of Laser Interferometric Gravitational Wave Detectors”.https://resolver. caltech.edu/CaltechETD:etd- 05242004- 154845. PhD thesis. Pasadena, California: California In- stitute of Technology, 2000

    2000

  8. [8]

    Optical cavity locking in non linear regime: A Reinforcement Learning strategy for future applications in Gravitational Wave Detectors

    A. Svizzeretto. “Optical cavity locking in non linear regime: A Reinforcement Learning strategy for future applications in Gravitational Wave Detectors”. Master’s thesis. Perugia: Universit` a di Perugia, 2024

    2024

  9. [9]

    Partial Observability and Domain Randomization in RL-Based Strategy for Optical Cavity Locking Optimization

    A. Svizzeretto and M. Bawaj. “Partial Observability and Domain Randomization in RL-Based Strategy for Optical Cavity Locking Optimization”. In:Proceedings of Fifth MODE Workshop on Differentiable Programming for Experiment Design — PoS(MODE2025). MODE2025. Sissa Medialab, Jan. 2026, p. 035.doi:10.22323/1.491.0035

    work page doi:10.22323/1.491.0035 2026

  10. [10]

    Bawaj and A

    M. Bawaj and A. Svizzeretto.Applying reinforcement learning to optical cavity locking tasks: con- siderations on actor-critic architectures and real-time hardware implementation. Jan. 2026. arXiv: 2509.14884 [physics.ins-det]

    work page arXiv 2026

  11. [11]

    Frequency Locking of an Optical Cavity Using a Time-Varying Kalman Filtering Approach

    S. Z. S. Hassen and I. R. Petersen. “Frequency Locking of an Optical Cavity Using a Time-Varying Kalman Filtering Approach”. In:IEEE Transactions on Control Systems Technology22.3 (May 2014), pp. 1143–1150.doi:10.1109/tcst.2013.2266693

    work page doi:10.1109/tcst.2013.2266693 2014

  12. [12]

    Sim-to-Real Transfer in Deep Reinforcement Learning for Robotics: a Survey

    W. Zhao, J. P. Queralta, and T. Westerlund. “Sim-to-Real Transfer in Deep Reinforcement Learning for Robotics: a Survey”. In:2020 IEEE Symposium Series on Computational Intelligence (SSCI). 2020, pp. 737–744.doi:10.1109/SSCI47803.2020.9308468

    work page doi:10.1109/ssci47803.2020.9308468 2020

  13. [13]

    Svizzeretto and M

    A. Svizzeretto and M. Bawaj.oreonspy: Optical REsONator Simulation in PYthon. Version v4.2.2. 2026.doi:10.5281/zenodo.19567205

    work page doi:10.5281/zenodo.19567205 2026

  14. [14]

    An introduction to Pound–Drever–Hall laser frequency stabilization

    E. D. Black. “An introduction to Pound–Drever–Hall laser frequency stabilization”. In:American Journal of Physics69.1 (2001), pp. 79–87.doi:10.1119/1.1286663. 10

    work page doi:10.1119/1.1286663 2001

  15. [15]

    The control of the Virgo interferometer for gravitational wave detection

    L. Barsotti. “The control of the Virgo interferometer for gravitational wave detection”.https://gwic. ligo.org/assets/docs/theses/Barsotti_Thesis.pdf. PhD thesis. Pisa: Universit` a degli Studi di Pisa, 2006

    2006

  16. [16]

    Swinkels and P

    B. Swinkels and P. Ruggi.Lock acquisition of the Advanced Virgo arm-cavities with reduced force. en. https://tds.virgo-gw.eu/ql/?c=8821. 2012

    2012

  17. [17]

    Measurement of the optical parameters of the Virgo interferometer

    F. Acernese et al. “Measurement of the optical parameters of the Virgo interferometer”. In:Applied Optics46.17 (May 2007), p. 3466.doi:10.1364/ao.46.003466

    work page doi:10.1364/ao.46.003466 2007

  18. [18]

    Virgo detector characterization and data quality: results from the O3 run

    F. Acernese et al. “Virgo detector characterization and data quality: results from the O3 run”. In: Classical and Quantum Gravity40.18 (Aug. 2023), p. 185006.doi:10.1088/1361-6382/acd92d

    work page doi:10.1088/1361-6382/acd92d 2023

  19. [19]

    Optical and noise studies for Advanced Virgo and filter cavities for quantum noise reduction in gravitational-wave interferometric detectors

    E. Capocasa. “Optical and noise studies for Advanced Virgo and filter cavities for quantum noise reduction in gravitational-wave interferometric detectors”.https : / / theses . hal . science / tel - 03141038/file/capocasa_thesis_2017.pdf. PhD thesis. Paris: Universit´ e Paris Diderot

  20. [20]

    Microring resonator channel dropping filters

    B. Little et al. “Microring resonator channel dropping filters”. In:Journal of Lightwave Technology 15.6 (June 1997), pp. 998–1005.doi:10.1109/50.588673

    work page doi:10.1109/50.588673 1997

  21. [21]

    Toward the end-to-end optimization of particle physics instruments with differentiable programming

    T. Dorigo et al. “Toward the end-to-end optimization of particle physics instruments with differentiable programming”. In:Reviews in Physics10 (June 2023), p. 100085.doi:10 . 1016 / j . revip . 2023 . 100085

    2023

  22. [22]

    The third generation of gravitational wave observatories and their science reach

    M. Punturo et al. “The third generation of gravitational wave observatories and their science reach”. In: Classical and Quantum Gravity27.8 (Apr. 2010), p. 084007.doi:10.1088/0264-9381/27/8/084007

    work page doi:10.1088/0264-9381/27/8/084007 2010

  23. [23]

    Sensitivity studies for third-generation gravitational wave observatories

    S. Hild et al. “Sensitivity studies for third-generation gravitational wave observatories”. In:Classical and Quantum Gravity28.9 (Apr. 2011), p. 094013.doi:10.1088/0264-9381/28/9/094013. 11 A Auxiliary materials for validation Fig. 4 shows the comparison between the measured DC transmission during the mode-matching measure- ments and the corresponding sim...

    work page doi:10.1088/0264-9381/28/9/094013 2011