pith. sign in

arxiv: 2605.13745 · v2 · pith:GETEZMIXnew · submitted 2026-05-13 · ⚛️ physics.optics

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

Thomas M. Baer This is my paper

Pith reviewed 2026-05-20 20:39 UTC · model grok-4.3

classification ⚛️ physics.optics
keywords frequency-doubled lasersLyapunov functionalamplitude noise suppressionchi2 couplingintracavity frequency doublingcoupled oscillatorsintensity stabilizationmultimode lasers
0
0 comments X

The pith

Chi-squared coupling in frequency-doubling crystals creates a Lyapunov functional that drives multimode lasers to constant total intensity.

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

Multimode intracavity frequency-doubled lasers exhibit amplitude noise suppression far beyond what independent-mode statistics predict. The paper traces this to the algebraic structure of the chi2 nonlinear wave mixing inside the doubling crystal. That structure defines a Lyapunov functional that decreases with every pass through the crystal. The decrease steers the combined modes toward a manifold of constant total intensity even while individual modes continue to exchange energy. Direct measurements in an Nd:YVO4-LBO laser confirm a 100-fold noise contrast between the unfiltered output and a Fabry-Perot-filtered single-mode slice at fixed bandwidth.

Core claim

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 mechanism explains amplitude noise suppression orders of magnitude beyond the predictions of independent-mode partition statistics.

What carries the argument

The Lyapunov functional constructed from the algebraic form of the chi-squared nonlinear coupling, which functions as a shared quadratic dissipative channel among the coupled oscillators.

If this is right

  • The stabilization should appear in any coupled-oscillator system whose equations share the same algebraic quadratic-coupling form.
  • Filtering individual modes or groups of modes should reveal substantially higher noise levels than the unfiltered combined output.
  • The effect is local to each pass through the nonlinear medium and does not require specific cavity losses or gain dynamics.
  • Analogous intensity stabilization should be observable in other physical systems exhibiting quadratic dissipative coupling between oscillators.

Where Pith is reading between the lines

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

  • Laser engineers could use the mechanism to achieve low-noise output without active feedback electronics.
  • The same algebraic descent might suppress fluctuations in non-optical systems such as certain parametric amplifiers or coupled mechanical resonators.
  • Systematic variation of crystal length or mode number would map the rate of descent toward the constant-intensity manifold.

Load-bearing premise

The observed noise suppression is dominated by the algebraic structure of the chi2 coupling and the associated Lyapunov descent rather than by unmodeled cavity losses, gain saturation, or detector effects.

What would settle it

A demonstration that comparable noise suppression occurs in an otherwise identical cavity when the doubling crystal is replaced by a linear absorber or beam splitter that lacks the chi2 interaction would falsify the claim that the Lyapunov mechanism is responsible.

Figures

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

Figure 1
Figure 1. Figure 1: FIG. 1. Per-mode IR amplitudes [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 2
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. 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. 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. 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.

Referee Report

1 major / 2 minor

Summary. The manuscript claims that the algebraic structure of the chi-squared coupled-wave equations in a frequency-doubling crystal admits a Lyapunov functional whose monotone decrease on each discrete crystal pass drives multimode dynamics toward a constant-intensity manifold. This mechanism is invoked to explain amplitude noise suppression far beyond independent-mode partition statistics. The claim is supported by an experimental demonstration in a CW Nd:YVO4-LBO intracavity-doubled laser that reports a 100-fold contrast between full-output and Fabry-Perot-filtered noise spectra at fixed detector bandwidth.

Significance. If the Lyapunov descent survives the transition from the isolated discrete map to the full continuous cavity dynamics, the result would identify a structurally generic route to intensity stabilization in quadratically coupled oscillator systems. The reported experimental contrast is large and the derivation is presented as parameter-free, which would be a notable strength if the mapping to the physical cavity is made rigorous.

major comments (1)
  1. [Lyapunov functional derivation and cavity model] The Lyapunov functional and its monotone decrease are constructed for a discrete per-pass map of the chi2 equations. The actual system is a continuous-wave laser in which the nonlinear interaction is distributed along the crystal length and embedded in a round-trip operator that includes linear losses, gain saturation, and possible dispersion. The manuscript does not show that the descent property is preserved under the continuous spatial derivatives or the full cavity map; without this step the 100-fold contrast cannot be unambiguously attributed to the Lyapunov mechanism rather than conventional nonlinear loss or filtering. (See the section deriving the functional and the experimental comparison.)
minor comments (2)
  1. [Experimental section] The abstract and main text should explicitly state the detector bandwidth and the Fabry-Perot filter parameters used to obtain the 100-fold contrast so that the statistical-averaging baseline can be recalculated by readers.
  2. [Throughout] Notation for the mode amplitudes and the quadratic coupling coefficients should be introduced once and used consistently; several symbols appear to be redefined between the model and the experimental discussion.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on the manuscript. The principal concern is the connection between the discrete Lyapunov analysis and the continuous cavity dynamics; we address this point below and outline the revisions that will be incorporated.

read point-by-point responses

  1. Referee: The Lyapunov functional and its monotone decrease are constructed for a discrete per-pass map of the chi2 equations. The actual system is a continuous-wave laser in which the nonlinear interaction is distributed along the crystal length and embedded in a round-trip operator that includes linear losses, gain saturation, and possible dispersion. The manuscript does not show that the descent property is preserved under the continuous spatial derivatives or the full cavity map; without this step the 100-fold contrast cannot be unambiguously attributed to the Lyapunov mechanism rather than conventional nonlinear loss or filtering. (See the section deriving the functional and the experimental comparison.)

    Authors: We agree that the Lyapunov functional and its monotone decrease are derived for the discrete per-pass map of the chi-squared coupled-wave equations, consistent with the modeling choice stated in the manuscript. This discrete formulation is the conventional approach for describing round-trip evolution in CW intracavity-doubled lasers. We acknowledge that an explicit demonstration of the descent property under continuous spatial propagation within the crystal, together with its compatibility with the full round-trip operator, would strengthen the attribution of the observed noise suppression. In the revised manuscript we will add a dedicated subsection that (i) computes the z-derivative of the Lyapunov functional along the continuous coupled-wave equations inside the crystal and shows it is non-positive, and (ii) discusses how intensity-dependent gain saturation and linear losses, being common to all modes, preserve the descent toward the constant-intensity manifold over successive round trips. These additions will clarify why the 100-fold experimental contrast is consistent with the Lyapunov mechanism rather than solely with conventional nonlinear loss or filtering. revision: yes

Circularity Check

1 steps flagged

Lyapunov functional constructed from chi2 equations makes descent a restatement of the model

specific steps
  1. self definitional [Abstract]
    "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."

    The Lyapunov functional is admitted by and constructed from the same chi2 coupled-wave equations that define the system dynamics. Its monotone decrease is therefore guaranteed by the algebraic structure of those equations, rendering the descent to the constant-intensity manifold a direct consequence of the model assumptions rather than an independent derivation or prediction.

full rationale

The paper derives a Lyapunov functional directly from the chi2 coupled-wave equations and shows its monotone decrease along the discrete-pass trajectories. This establishes the constant-intensity manifold by algebraic construction within the isolated model. While the experimental 100-fold noise contrast is reported, the central mechanism reduces to properties already encoded in the starting equations without an independent benchmark or falsifiable prediction outside the fitted data. The continuous intracavity embedding and other cavity effects are not shown to preserve the same descent property, leaving the explanatory power tied to the discrete model assumptions.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the existence and monotonicity of a Lyapunov functional constructed from the chi2 coupled-wave equations; no free parameters, additional axioms, or invented entities are stated in the abstract.

axioms (1)
  • domain assumption The intracavity dynamics are governed by the standard chi2 coupled-wave equations without additional loss or gain terms that would violate the Lyapunov descent.
    Invoked to establish the monotone decrease under each crystal pass.

pith-pipeline@v0.9.0 · 5665 in / 1233 out tokens · 60281 ms · 2026-05-20T20:39:55.509202+00:00 · methodology

Review history (2 revisions) →

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.

  • IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel echoes
    ?
    echoes

    ECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.

    dM4/dz ≤ 0, z∈[0, Lcrystal]. ... The rigorous proof ... uses the exact sech²/tanh² coupled-wave solution and the Chebyshev integral inequality ... with equality only when I(t) is constant.

  • IndisputableMonolith/Cost.lean Jcost_pos_of_ne_one echoes
    ?
    echoes

    ECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.

    The intensity-channel mechanism is intuitive. The crystal depletes intensity as I², so brighter time-slices lose proportionally more energy than dim time-slices, flattening I(t) toward its mean.

  • IndisputableMonolith/Foundation/BranchSelection.lean branch_selection echoes
    ?
    echoes

    ECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.

    The mechanism rests on the algebraic structure of the coupling — a coherent superposition of oscillators sharing a quadratic dissipative channel

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

29 extracted references · 29 canonical work pages

  1. [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)

    work page 1986

  2. [2]

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

    work page 1964

  3. [3]

    Wu and P

    X.-G. Wu and P. Mandel, Second-harmonic generation in a multimode laser cavity, J. Opt. Soc. Am. B4, 1870 (1987)

    work page 1987

  4. [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)

    work page 1990

  5. [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)

    work page 1993

  6. [6]

    Erneux and P

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

    work page 1995

  7. [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)

    work page 1998

  8. [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. [9]

    C. S. Adams, G. T. Maker, and A. I. Ferguson, FM op- eration of Nd:YAG lasers with standing wave and ring cavity configurations, Opt. Commun.76, 127 (1990)

    work page 1990

  10. [10]

    Tsunekane, N

    M. Tsunekane, N. Taguchi, and H. Inaba, Elimination of chaos in a multilongitudinal-mode, diode-pumped, 6- W continuous-wave, intracavity-doubled Nd:YAG laser, Opt. Lett.22, 1000 (1997), see also OSA TOPS Vol. 10, p. 64 (1997)

    work page 1997

  11. [11]

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

    work page 1999

  12. [12]

    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)

    work page 1979

  13. [13]

    Sargent III, M

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

    work page 1974

  14. [14]

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

    work page 1986

  15. [15]

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

    work page 1991

  16. [16]

    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)

    work page 1993

  17. [17]

    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)

    work page 1993

  18. [18]

    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...

    work page 2000

  19. [19]

    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 ...

    work page

  20. [20]

    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...

    work page

  21. [21]

    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- verse FFT, the coupled-wave equations (A4) are inte- grated through the forward pass, the boundary condi- tionδψis applied at the high-reflector, the...

    work page

  22. [22]

    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...

    work page

  23. [23]

    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 ...

    work page

  24. [24]

    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 ...

    work page

  25. [25]

    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...

    work page

  26. [26]

    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 ′⟩

    work page

  27. [27]

    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...

    work page

  28. [28]

    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- ...

    work page 1952

  29. [29]

    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...

    work page 1964