arxiv: 2605.13745 · v1 · submitted 2026-05-13 · ⚛️ physics.optics

Recognition: 2 theorem links

· Lean Theorem

Amplitude Noise Suppression in Frequency-Doubled Lasers: A Lyapunov Mechanism for Intensity Stabilization in Coupled Oscillator Systems

Thomas M. Baer

Authors on Pith no claims yet

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

classification ⚛️ physics.optics

keywords amplitude noise suppressionfrequency-doubled lasersLyapunov functionalchi2 nonlinear opticsintracavity doublingcoupled oscillatorsintensity stabilization

0 comments

The pith

Chi-squared coupled-wave dynamics in frequency-doubled lasers admit a Lyapunov functional that suppresses amplitude noise by driving modes to constant intensity.

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

This paper establishes that the nonlinear chi2 interactions inside the frequency-doubling crystal create a Lyapunov functional whose steady decrease forces the laser modes toward states of constant total intensity. A reader would care because this explains why multimode intracavity frequency-doubled lasers show amplitude noise suppression far exceeding what independent mode statistics predict, as seen in experiments with 100-fold noise contrast. The mechanism arises from the algebraic form of the quadratic dissipative coupling shared among the oscillators, offering a passive stabilization route. It is confirmed experimentally in an Nd:YVO4-LBO setup where full output noise drops dramatically compared to filtered single-mode output.

Core claim

We show that the chi2 coupled-wave dynamics in the doubling crystal admit a Lyapunov functional whose monotone decrease under each crystal pass establishes a constant-intensity manifold as the per-pass descent target of the mode dynamics. This algebraic structure of the coupling, a coherent superposition of oscillators sharing a quadratic dissipative channel, accounts for the observed noise suppression in multimode intracavity frequency-doubled lasers.

What carries the argument

The Lyapunov functional constructed from the chi2 coupled-wave equations, which monotonically decreases to the constant-intensity manifold with each pass through the crystal.

If this is right

Amplitude noise can be suppressed orders of magnitude beyond partition noise limits in multimode lasers.
The effect holds for any coupled oscillator system sharing the same quadratic dissipative channel algebra.
Experimental validation shows 100-fold noise reduction contrast in Nd:YVO4-LBO lasers at fixed detector bandwidth.
The stabilization occurs passively through the crystal dynamics without external intervention.

Where Pith is reading between the lines

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

Similar effects could occur in other chi2-based devices where multiple frequencies interact quadratically.
The passive nature suggests applications in systems where active feedback is impractical.
Generalization to other algebraic forms of coupled oscillators might reveal analogous stabilization in non-optical domains.

Load-bearing premise

The chi2 coupled-wave equations dominate the intracavity dynamics so that the Lyapunov functional decreases monotonically without disruption from other loss channels or cavity effects.

What would settle it

Measuring equal amplitude noise levels in the full multimode output and in the Fabry-Perot-filtered single-mode output at the same bandwidth would falsify the claim, as the mechanism predicts a large contrast due to the multi-mode coupling.

Figures

Figures reproduced from arXiv: 2605.13745 by Thomas M. Baer.

**Figure 2.** Figure 2: FIG. 2. Total green output [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗

**Figure 3.** Figure 3: FIG. 3. Primary measurement: full Millennia Vs green output [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗

**Figure 4.** Figure 4: FIG. 4. Mode-count dependence: filtered green output for [PITH_FULL_IMAGE:figures/full_fig_p010_4.png] view at source ↗

**Figure 5.** Figure 5: FIG. 5. PCA decomposition of the per-mode intensity fluc [PITH_FULL_IMAGE:figures/full_fig_p011_5.png] view at source ↗

read the original abstract

Multimode intracavity frequency-doubled lasers can reach states of amplitude noise suppression orders of magnitude beyond the predictions of independent-mode partition statistics. We show that the chi2 coupled-wave dynamics in the doubling crystal admit a Lyapunov functional whose monotone decrease under each crystal pass establishes a constant-intensity manifold as the per-pass descent target of the mode dynamics. We confirm the mechanism in an intracavity frequency-doubled Nd:YVO4-LBO laser, observing a 100 fold contrast between full and Fabry-Perot-filtered output noise at fixed detector bandwidth, well beyond the statistical-averaging baseline. The mechanism rests on the algebraic structure of the coupling, a coherent superposition of oscillators sharing a quadratic dissipative channel, and is therefore a candidate for analogous noise-suppression effects in other coupled oscillator systems with the same algebraic form.

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.

Desk Editor's Note private letter to a colleague

The paper derives a Lyapunov functional from the algebraic structure of chi2 coupling that drives intensity to a constant manifold per pass, and backs it with a 100-fold experimental noise contrast in an Nd:YVO4-LBO laser.

read the letter

The main takeaway is that the chi2 coupled-wave dynamics alone admit a Lyapunov functional whose monotone decrease per crystal pass targets a constant-intensity state, and the authors report a 100-fold noise suppression in experiment that exceeds partition-statistics expectations. This algebraic framing is new relative to the empirical models cited in the abstract; it treats the modes as a coherent superposition sharing a quadratic dissipative channel rather than independent oscillators. The experimental contrast between full and Fabry-Perot-filtered output at fixed bandwidth is a clear data point in favor of the mechanism operating in a real cavity. The derivation appears to rest on the specific form of the chi2 equations without fitted parameters, which is a point in its favor. The soft spot is the round-trip map. The stress-test concern is valid on the abstract alone: gain saturation in the Nd:YVO4, pump depletion, and linear losses are not obviously included in the per-pass descent argument. If those terms can make the functional increase along some trajectories, the constant-intensity manifold is no longer guaranteed. The full text needs to show the explicit functional and confirm dV/dt remains non-positive once the complete operator is restored. This work is aimed at nonlinear-optics and precision-laser groups who care about first-principles noise control. It is coherent enough on its own terms to merit peer review, though the referee should press on the full-cavity closure.

Referee Report

2 major / 2 minor

Summary. The manuscript claims that the χ² coupled-wave dynamics in the frequency-doubling crystal admit a Lyapunov functional whose monotone per-pass decrease drives the intracavity mode amplitudes to a constant-intensity manifold, thereby suppressing amplitude noise far below independent-mode partition statistics. This algebraic mechanism is asserted to arise from the coherent superposition of oscillators sharing a quadratic dissipative channel. The claim is supported by an experimental demonstration in an intracavity frequency-doubled Nd:YVO4-LBO laser that reports a 100-fold contrast in output noise between the full multimode case and a Fabry-Perot-filtered case at fixed detector bandwidth.

Significance. If the Lyapunov functional can be shown to remain monotonically non-increasing once the full round-trip map (including gain saturation, pump depletion, and linear losses) is restored, the result supplies a concrete, parameter-free stabilization mechanism for coupled-oscillator systems with quadratic dissipation. The reported experimental contrast is large enough to be of practical interest for low-noise frequency-doubled sources, and the algebraic framing suggests possible extensions to other nonlinear optical or classical oscillator networks.

major comments (2)

[Abstract / central claim paragraph] Abstract and the paragraph stating the central claim: the existence of a Lyapunov functional V for the χ² coupled-wave equations is asserted and said to produce monotone decrease to the constant-intensity manifold, but neither the explicit functional nor the calculation establishing dV/dt ≤ 0 (or per-pass ΔV ≤ 0) is exhibited. Because the skeptic correctly notes that non-χ² terms (Nd:YVO4 gain saturation, pump depletion, linear cavity losses) are present in the round-trip operator, the manuscript must demonstrate that these terms do not reverse the sign of the derivative on any physically relevant trajectory; without that step the descent argument is incomplete.
[Round-trip map / composite dynamics] Section describing the round-trip map: the paper isolates the χ² crystal dynamics for the Lyapunov construction but does not provide the composite map that folds in the gain medium and cavity losses, nor does it verify that the functional remains a descent function for the composite operator. This verification is load-bearing for the claim that the mechanism survives in a real laser cavity.

minor comments (2)

[Experimental results] The experimental contrast is quoted as “100 fold” but the precise definition of the noise metric (e.g., integrated RIN over what bandwidth, normalization to shot-noise level) and the statistical-averaging baseline calculation are not stated with sufficient detail to allow direct replication.
[Theory / coupled-wave equations] Notation for the coupled-wave amplitudes and the quadratic dissipative term should be introduced once, with all symbols defined before the Lyapunov functional is introduced.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments. The points raised about the explicit presentation of the Lyapunov functional and its extension to the composite round-trip dynamics are well taken. We have revised the manuscript to supply the missing explicit expressions, calculations, and verifications while preserving the original claims.

read point-by-point responses

Referee: [Abstract / central claim paragraph] Abstract and the paragraph stating the central claim: the existence of a Lyapunov functional V for the χ² coupled-wave equations is asserted and said to produce monotone decrease to the constant-intensity manifold, but neither the explicit functional nor the calculation establishing dV/dt ≤ 0 (or per-pass ΔV ≤ 0) is exhibited. Because the skeptic correctly notes that non-χ² terms (Nd:YVO4 gain saturation, pump depletion, linear cavity losses) are present in the round-trip operator, the manuscript must demonstrate that these terms do not reverse the sign of the derivative on any physically relevant trajectory; without that step the descent argument is incomplete.

Authors: We agree that the original submission did not exhibit the explicit functional or the derivative calculation. The revised manuscript now defines the Lyapunov functional explicitly and provides the algebraic steps establishing dV/dt ≤ 0 (or per-pass ΔV ≤ 0) for the isolated χ² coupled-wave equations. We have also added an analysis of the composite round-trip operator that incorporates gain saturation, pump depletion, and linear losses. Both analytical bounds and numerical integration over the experimental parameter range show that these additional terms do not reverse the sign of the change in V on physically relevant trajectories, thereby completing the descent argument. revision: yes
Referee: [Round-trip map / composite dynamics] Section describing the round-trip map: the paper isolates the χ² crystal dynamics for the Lyapunov construction but does not provide the composite map that folds in the gain medium and cavity losses, nor does it verify that the functional remains a descent function for the composite operator. This verification is load-bearing for the claim that the mechanism survives in a real laser cavity.

Authors: We acknowledge that the original manuscript presented only the isolated χ² dynamics. The revised version now supplies the explicit composite round-trip map that folds in the Nd:YVO4 gain (including saturation and pump depletion) and linear cavity losses. We verify that the Lyapunov functional remains monotonically non-increasing for this full operator by providing both a perturbative analytical argument (showing non-positive contributions from the extra terms) and direct numerical simulations of the complete map using parameters matched to the Nd:YVO4-LBO experiment. This confirms the mechanism is robust in the physical cavity. revision: yes

Circularity Check

0 steps flagged

No load-bearing circularity; Lyapunov functional derived from chi2 algebraic structure

full rationale

The derivation claims that the chi2 coupled-wave equations admit a Lyapunov functional whose monotone decrease establishes the constant-intensity manifold. No quoted step shows the functional defined in terms of the target manifold itself, no fitted parameter renamed as prediction, and no self-citation chain invoked as the sole justification for monotonicity. The mechanism is presented as following from the algebraic form of the quadratic dissipative channel. Experimental observations are reported separately. This qualifies as a self-contained derivation against the stated model with only minor risk of unverified extension to full round-trip dynamics.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the assumption that chi2 coupled-wave dynamics govern the mode evolution and admit a Lyapunov functional; no free parameters or invented entities are mentioned in the abstract.

axioms (1)

domain assumption The intracavity dynamics are governed by chi2 coupled-wave equations without dominant competing loss mechanisms.
Invoked to justify that the Lyapunov functional controls the per-pass descent.

pith-pipeline@v0.9.0 · 5434 in / 1317 out tokens · 52662 ms · 2026-05-14T17:35:27.974561+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

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

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

M4 ≡ ⟨I(t)²⟩ / ⟨I(t)⟩² ... dM4/dz ≤ 0 ... constant-intensity manifold as the per-pass descent target

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

31 extracted references · 1 canonical work pages

[1]

Baer, Large-amplitude fluctuations due to longitudi- nal mode coupling in diode-pumped intracavity-doubled Nd:YAG lasers, J

T. Baer, Large-amplitude fluctuations due to longitudi- nal mode coupling in diode-pumped intracavity-doubled Nd:YAG lasers, J. Opt. Soc. Am. B3, 1175 (1986)

1986
[2]

S. E. Harris and R. Targ, FM oscillation of the He-Ne laser, Appl. Phys. Lett.5, 202 (1964)

1964
[3]

C. R. Wu and L. Mandel, Deterministic chaos in a laser with a saturable absorber, J. Opt. Soc. Am. B4, 1870 (1987)

1987
[4]

Wiesenfeld, C

K. Wiesenfeld, C. Bracikowski, G. James, and R. Roy, Observation of antiphase states in a multimode laser, Phys. Rev. Lett.65, 1749 (1990)

1990
[5]

Chaos in a multimode solid-state laser system,

R. Roy, C. Bracikowski, and G. E. James, Dynamics of a multimode laser with nonlinear, birefringent intracavity elements, inRecent Developments in Quantum Optics, edited by R. Inguva (Plenum Press, New York, 1993) see also C. Bracikowski and R. Roy, “Chaos in a multimode solid-state laser system,” Chaos1, 49 (1991)

1993
[6]

Erneux and P

T. Erneux and P. Mandel, Minimal equations for an- tiphase dynamics in multimode lasers, Phys. Rev. A52, 4137 (1995)

1995
[7]

Kozyreff and P

G. Kozyreff and P. Mandel, Antiphase dynamics and self- pulsing due to a low-frequency spatial population grating in a multimode laser, Phys. Rev. A58, 4946 (1998)

1998
[8]

M. E. Pietrzyk and M. B. Danailov, Elimination of chaos in multimode, intracavity-doubled lasers in the presence of spatial hole-burning, arXiv preprintnlin/0012051 (2000)

work page arXiv 2000
[9]

Adams, P

M. Adams, P. D. Maker, and I. Ferguson, FM operation of Nd:YAG lasers with standing wave and ring cavity configurations, Opt. Commun.76, 127 (1990)

1990
[10]

Tsunekane, N

M. Tsunekane, N. Taguchi, and H. Inaba, Reduction of intensity noise in a diode-pumped intracavity-doubled Nd:YAG laser, Opt. Lett.22, 1000 (1997), see also OSA TOPS Vol. 10, p. 64 (1997)

1997
[11]

D. W. Anthon, Passive FM laser operation and the sta- bility of intracavity-doubled lasers, Appl. Opt.38, 5144 (1999)

1999
[12]

Matos, O

L. Matos, O. D. M¨ ucke, J. Chen, and F. X. K¨ artner, Carrier-envelope phase dynamics and noise analysis in octave-spanning Ti:sapphire lasers, Opt. Express14, 2497 (2006)

2006
[13]

D. V. Sutyrin, N. Poli, N. Beverini, S. V. Chepurov, M. Prevedelli, M. Schioppo, F. Sorrentino, M. G. Tarallo, and G. M. Tino, Frequency noise performances of a Ti:sapphire optical frequency comb stabilized to an opti- cal reference, Opt. Commun.291, 291 (2013)

2013
[14]

Ikeda, Multiple-valued stationary state and its insta- bility of the transmitted light by a ring cavity system, Opt

K. Ikeda, Multiple-valued stationary state and its insta- bility of the transmitted light by a ring cavity system, Opt. Commun.30, 257 (1979)

1979
[15]

Sargent III, M

M. Sargent III, M. O. Scully, and W. E. Lamb, Jr.,Laser Physics(Addison-Wesley, 1974)

1974
[16]

A. E. Siegman,Lasers(University Science Books, 1986)

1986
[17]

K. Y. Tsang, R. E. Mirollo, S. H. Strogatz, and K. Wiesenfeld, Dynamics of a globally coupled oscilla- tor array, Physica D48, 102 (1991)

1991
[18]

Fabiny, P

L. Fabiny, P. Colet, R. Roy, and D. Lenstra, Coher- ence and phase dynamics of spatially coupled solid-state lasers, Phys. Rev. A47, 4287 (1993)

1993
[19]

Silber, L

M. Silber, L. Fabiny, and K. Wiesenfeld, Stability results for in-phase and splay-phase states of solid-state laser arrays, J. Opt. Soc. Am. B10, 1121 (1993)

1993
[20]

S. H. Strogatz, From Kuramoto to Crawford: exploring the onset of synchronization in populations of coupled oscillators, Physica D143, 1 (2000). Appendix A: Model details The simulation treats each longitudinal mode as a complex field amplitudeE i(t) (i= 1, . . . , N), evolved per round trip through three sequential physical stages: gain saturation includ...

2000
[21]

Each operator acts on the full mode vector{E i}

Round-trip structure A single round tripE (n) i →E (n+1) i is the composition E(n+1) i = ˆN ˆC ˆG[E(n) i ] +ξ (n) i ,(A1) where ˆGis the gain stage, ˆCis theχ (2) crystal stage, and ˆNis the linear cavity loss including output coupling. Each operator acts on the full mode vector{E i}. The additive termξ (n) i models spontaneous emission injected once per ...
[22]

Gain stage with spatial-hole-burning The gain stage applies a per-mode amplification Gi = exp " g0/2 1 +P j βij|Ej|2/Is −ℓ/2 # ,(A2) whereg 0 is the small-signal gain,ℓis the round-trip linear loss,I s is the saturation intensity, andβ ij =β(|i−j|) is the standing-wave spatial-hole-burning cross-saturation coefficient between modesiandj. For a gain medium...
[23]

Nonlinear stage: split-step FFT computation The crystal stage ˆCis computed in the time domain. Given the input mode vector{E (n) i }at the crystal en- trance, the time-domain fieldE(t) is constructed by in- 15 verse FFT, the coupled-wave equations (A4) are inte- grated through the forward pass, the boundary condi- tionδψis applied at the high-reflector, ...
[24]

Noise sources Two noise mechanisms are included in the simulation. Spontaneous emission is modeled as additive complex Gaussian noise injected once per mode per round trip: ξ(n) i = q Rsp/2 ξ(n,R) i +iξ (n,I) i ,(A5) whereξ (n,R) i , ξ(n,I) i ∼ N(0,1) are independent standard normal variates andR sp =n sphν/τc is the spontaneous emission rate per mode, wi...
[25]

The cavity length is taken to be 42 cm based on direct measurement of the IR longitu- dinal beat-note frequency at 358 MHz (corresponding to free spectral rangec/2L= 358 MHz)

Parameter values All simulations use parameters matched to the plat- form studied here. The cavity length is taken to be 42 cm based on direct measurement of the IR longitu- dinal beat-note frequency at 358 MHz (corresponding to free spectral rangec/2L= 358 MHz). TABLE I. Simulation parameters used in all production runs. σgain was calibrated against the ...
[26]

Define the instantaneous fundamental intensity entering the crystal, I0(t)≡ |E(t)| 2.(B2) The time average⟨·⟩is taken over one round-trip period T= 2π/Ω

Setup Let the total intracavity field at the entrance of the crystal be E(t) = X j Aj exp i(ωjt+φ j) ,(B1) with mode amplitudesA j and phasesφ j at frequencies ωj =ω 0 +jΩ. Define the instantaneous fundamental intensity entering the crystal, I0(t)≡ |E(t)| 2.(B2) The time average⟨·⟩is taken over one round-trip period T= 2π/Ω. The fourth-moment ratio is M4 ...
[27]

Exact coupled-wave solution Within the crystal, the fundamental fieldEand second-harmonic fieldE grn evolve under the coupled- wave equations of Appendix A, Sec. A.3. For perfect phase matching and on timescales short compared to the round-trip period — so that each time slicetof the input waveformI 0(t) propagates through the crystal indepen- dently — th...
[28]

ComputingdM 4/dz DifferentiatingM 4 =B/A 2 with respect tozand using ∂I/∂z=−S ′, dA dz =−⟨S ′⟩, dB dz =−2⟨IS ′⟩,(B6) so that dM4 dz = 1 A2 dB dz − 2B A3 dA dz = 2 A3 ⟨I2⟩⟨S′⟩ − ⟨I⟩⟨IS ′⟩ .(B7) Equivalently, usingB=⟨I 2⟩andA=⟨I⟩, dM4 dz = 2 A3 B⟨S ′⟩ −A⟨IS ′⟩ .(B8) The sign ofdM 4/dzis determined by the sign ofB⟨S ′⟩ − A⟨IS ′⟩
[29]

We need thatIis monotoni- cally increasing inI 0 within the regime of validity stated below

Co-monotonicity inI 0 At each fixedz >0, bothI(t, z) andS(t, z) are func- tions ofI 0(t) alone, since each time slice evolves indepen- dently through the crystal. We need thatIis monotoni- cally increasing inI 0 within the regime of validity stated below. Co-monotonicity ofS ′/IwithIis established alongside the Chebyshev step in Sec. B 5. Writingf(x) =x 2...
[30]

This is Hardy-Littlewood-P´ olya,Inequalities(Cam- bridge, 1952), Theorem 236

Chebyshev integral inequality For any two real-valued functionsf, gthat are both monotonically increasing (or both monotonically decreas- ing) in a common variable and averaged with respect to a common probability measureµ, Z f g dµ≥ Z f dµ· Z g dµ,(B11) with equality if and only ifforgis constantµ- a.e. This is Hardy-Littlewood-P´ olya,Inequalities(Cam- ...

1952
[31]

No restriction on mode count, no assumption of weak coupling, and no prior assumption on phase organization is needed

Scope and remarks Scope.The result (B18) holds for arbitrary mode am- plitudes and phases at the input to the crystal, at every positionzwithin the forward pass, requiring only that the peak instantaneous single-pass conversion efficiency satisfyη peak = tanh2(κL p I0,max)<0.70. No restriction on mode count, no assumption of weak coupling, and no prior as...

1964