Pith. sign in

REVIEW 3 major objections 5 minor 52 references

Periodic orbits underlying spatiotemporal chaos in the Lugiato-Lefever model

T0 review · 3 major / 5 minor · reviewed 2026-08-04 · deepseek-v4-flash

Pith's one-line read Unstable periodic orbits, computed as exact solutions of the Lugiato–Lefever equation, organize the chaotic dynamics of optical cavities: chaos is a walk that intermittently shadows these oscillating Turing rolls.

desk verdict Solid new OTR existence and bifurcation results, but the large-domain shadowing claim leans on fitted domain sizes and a 2D projection; worth serious refereeing with a focus on that evidence. read the letter →

arxiv 2509.10283 v1 pith:Z2IILGEQ submitted 2025-09-12 nlin.PS nlin.CDphysics.optics

classification nlin.PSnlin.CDphysics.optics MSC 35B3235Q5537C27
keywords Lugiato-LefeverequationoscillatingTuringrollsperiodicorbitsspatiotemporalchaosmodulationalinstabilitybifurcationanalysisopticalfrequencycombsexactcoherentstates
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper seeks to show that the chaotic light patterns produced in a driven optical cavity are not structureless randomness but are organized by a family of exact, time-periodic solutions of the underlying Lugiato–Lefever equation. The authors compute stable and unstable oscillating Turing rolls (OTRs) in the modulation-instability regime and trace them through parameter space, revealing how they emerge from Hopf bifurcations of stationary Turing rolls and reconnect to soliton and breather branches. They then show, in direct numerical simulations, that the chaotic trajectory intermittently adopts the spatial and temporal signature of these OTR solutions before switching to another pattern or to disorganized dynamics. If the claim holds, optical turbulence in this system becomes describable within the dynamical-systems picture of chaos as a walk among exact invariant solutions, the same picture used for wall-bounded fluid turbulence.

What carries the argument

The central object is the OTR (oscillating Turing roll): a time-periodic solution of the LLE in which adjacent rolls oscillate in anti-phase, defined as a fixed point x* of the return map F^T(x*) − x* = 0 over the temporal period T. The authors compute these solutions with a Jacobian-free Newton–Krylov solver in a truncated Fourier basis, continue the branches in detuning and pump parameters, and characterize stability both on the full cavity domain L and on the reduced domain L0 = 2L/N (or 3L/N for the relative orbit OTR II). The shadowing mechanism is demonstrated by comparing spatiotemporal intensity portraits and two-dimensional phase-space projections of chaotic simulations with the exa

What would settle it

Directly measure the state-space distance between the chaotic LLE trajectory and the OTR I/II solutions over time in the regime of Figure 5 (f^2=16, ζ0=0.45). If the minimal distance does not approach the solver tolerance for intervals of order one OTR period, or if shifting the fitted domain by ±1% destroys the agreement in the spatiotemporal portraits, the shadowing claim fails. A complementary check is to perform a recurrence analysis: count near-returns of the chaotic trajectory to the same point in the Poincaré section and compare the associated periods with the OTR periods.

Watch

Extended reading notes

Core claim

The central claim is that unstable periodic orbit solutions, called oscillating Turing rolls (OTR), are exact invariant solutions of the Lugiato–Lefever equation and that they serve as the elementary building blocks of both temporal and spatiotemporal chaos in the modulation-instability regime. The authors report three OTR families, follow their solution branches by numerical continuation, and connect them to continuous-wave, Turing-roll, soliton, and breather branches, completing the classical phase diagram. In small cavities, the chaotic field spontaneously switches between OTR I and OTR III patterns. In large cavities, localized patches of anti-phase and drifting roll oscillations match O

Load-bearing premise

The load-bearing premise is that the OTRs computed on the fitted periods 0.988 L0,I and 0.970 L0,II actually represent the coherent structures observed in the chaotic field; if the pattern size in the cavity is merely an artifact of measuring roll-maxima distances while interacting waves distort the structure, the shadowing match in Figure 5 is a fitting product rather than evidence that the chaotic trajectory visits exact solutions.

Editorial extensions

If this is right

  • The OTR solution branches complete the phase diagram of the LLE by connecting continuous waves, Turing rolls, solitons, and breathers through bifurcations.
  • The L0-stable OTRs (stable on their own period but not on the full cavity) can be observed directly in microresonators as spontaneous temporal switching of intracavity power.
  • The unstable OTRs provide a language for describing optical turbulence: a chaotic state is a concatenation of visits to different exact periodic solutions.
  • Additional OTR families are expected, which would further refine the description of spatiotemporal patterns.
  • The approach extends to dispersion-engineered and coupled resonators, since the OTRs are solutions of the LLE, the universal mean-field model for these systems.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the shadowing picture holds, statistical properties of the chaotic field (such as residence times and transition rates between patterns) could be computed from the Floquet multipliers and stable/unstable manifolds of the OTRs using periodic-orbit theory, though the paper does not attempt this.
  • The sensitivity of the OTR profile to a few percent change in domain size, shown in the supplementary materials, suggests that in larger cavities pattern interactions will distort the embedded periodic orbits; a quantitative theory of this distortion would be needed to make shadowing predictive rather than post hoc.
  • The relative periodic orbit OTR II (periodic in a moving frame) hints that drift-type symmetries matter in optical cavities; tracking such relative orbits in other driven-dissipative systems may reveal similar organizing states.
  • Because the OTRs exist on domains that are rational fractions of the cavity, they may correspond to rational harmonic combs in the frequency domain; an experimental check is whether the corresponding comb teeth show the predicted phase and amplitude oscillations.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper computes unstable time-periodic invariant solutions—'oscillating Turing rolls' (OTRs)—of the Lugiato–Lefever equation (LLE) using Newton–Krylov continuation and arclength continuation. It traces OTR branches from Turing-roll Hopf bifurcations, identifies three families (OTR I, II, III), characterizes their stability on the full cavity domain and on reduced L0 domains, and connects them to continuous waves, Turing rolls, solitons, and breathers. The central claim is that chaotic LLE dynamics in small and large cavities transiently shadow these OTR periodic orbits, making them the 'fundamental building blocks' of both temporal and spatiotemporal chaos.

Significance. If substantiated, the result would be a significant transfer of the exact-coherent-states framework from fluid dynamics to nonlinear optics: optical turbulence in a Kerr cavity would be described as a chaotic walk among exact invariant solutions. The existence, continuation, and stability analysis of the OTR branches is a credible numerical contribution, and the methods are modern and appear internally consistent. The paper also offers a potentially testable prediction in the form of L0-stable OTRs that could be observed in microresonators. However, the manuscript's strongest claim—that chaotic trajectories shadow OTRs in large domains—rests on qualitative visual and low-dimensional evidence, and one part of that evidence is obtained by fitting the OTR spatial period to the chaotic pattern size. The significance therefore depends on closing this gap.

major comments (3)
  1. [Fig. 5 and surrounding text (p. 5)] The large-domain shadowing evidence is partly circular. The text states that the OTR solutions in Fig. 5(c1,c2) are computed on domains 0.988L0,I and 0.970L0,II, 'chosen according to the size of the MI patterns from (b1,b2)'. Supplementary Fig. 7 explicitly shows that even a few percent change in domain size significantly modifies the OTR profile and amplitude. Thus matching the chaotic pattern's wavelength and amplitude is to a considerable degree enforced by construction. The authors should either determine the OTR periods independently (e.g., from the continuation branch structure or from a symmetry argument) or quantify how the shadowing comparison degrades as the domain size is varied within the plausible range.
  2. [Fig. 5(d), phase-space projection] No quantitative shadowing statistic is provided. The confirmation of shadowing in Fig. 5(d) is a two-dimensional projection of |ψ| at two fixed θ positions. A chaotic trajectory can appear recurrent in such a low-dimensional projection even when its full state-space distance to the OTR is not small. The authors should report a normalized L2 (or equivalent) distance between the chaotic segments and the corresponding OTR solutions, minimized over allowed phase shifts, time shifts, and, for OTR II, frame drift. Without this, the assertion that 'spatiotemporal chaos as a chaotic trajectory in state space transiently shadows the OTR periodic orbits' is not quantitatively supported.
  3. [Abstract and Fig. 1] The claim that the results 'complete the classical phase diagram of the optical cavity' is stronger than the evidence presented. The continuation covers sixteen OTR I branches with N=8–12 and specific pump values (f^2=4, 8, 16 are used in figures), plus OTR II and OTR III for selected N. This is not a systematic exploration of the full parameter space, and the phrase 'completes' overreaches. I recommend either tempering the claim to 'connects the main dynamical regimes' or providing a systematic coverage map, e.g., in the (ζ0, f) plane, showing which branches have been tracked.
minor comments (5)
  1. [Introduction, sentence after Eq. (1)] The sentence 'where the two control parameters are the pump-cavity detuning ζ0 is and pump laser strength f' contains a grammatical error ('ζ0 is').
  2. [Abstract and main text] The word 'intermittancy' should be 'intermittency'.
  3. [Fig. 4 caption] The caption says panels (d1,d2) show OTR I and OTR III solutions 'which fit the observed MI patterns'; this wording assumes the conclusion. It would be clearer to say 'computed for comparison'.
  4. [Supplementary Materials, Fig. 8] Fig. 8(b) shows the OTR III bifurcation diagram for N=10 at f^2=8, while Fig. 4 uses N=11 at f^2=8. Please clarify whether the OTR III solution used in Fig. 4 is on the N=11 branch and whether the bifurcation structure is representative.
  5. [General] The manuscript would benefit from a data-availability statement indicating whether the continuation code and solution branches will be made publicly available, given the reproducibility-oriented methods.

Circularity Check

1 steps flagged · score 6.0 of 10

Large-domain shadowing claim is partly built from fitted OTR domain sizes; OTR existence and small-domain switching remain independent.

  1. fitted input called prediction [Main text, paragraph describing Figure 5 (Section 'Spatiotemporal chaos emerges in large cavities'), p. 4–5; see also Figure 5 caption (c1,c2)]
    "In order to provide a fair comparison of the observed MI patterns shown in Figure 5(b1,b2) with the OTR periodic orbits, we measure the size of the patterns based on the distance between roll maxima, that gives us values 0.988L0,I and 0.970L0,II respectively. Then we compute OTR I and OTR II solutions on these domains, which results in an accurate fit of the MI patterns by the periodic orbits in terms of both spatiotemporal periods and intracavity field amplitudes, see Figure 5(b,c)."

    The spatial periods of the OTRs used for the shadowing comparison are not taken from the OTR branch itself; they are measured from the chaotic MI patterns (0.988 L0,I and 0.970 L0,II). Computing an OTR on a domain equal to the measured pattern size forces its spatial period to match the chaotic pattern's wavelength by construction. The subsequent agreement in panel (b)/(c) is therefore partly guaranteed by the fitting step, not an independent confirmation of shadowing. This is amplified by Supplementary Figure 7, which states that 'even a few percent variations of the domain size significantly modify the spatio-temporal profile' — so the fitted domain is a strong tuning parameter. To demonstrate shadowing independently one would need a quantitative state-space distance comparison, ideally

full rationale

The OTR solutions themselves are computed by solving Eq. (2) with Newton-Krylov methods and continuation; their existence and bifurcation structure do not depend on the chaotic trajectories and are credible. The small-domain temporal switching demonstrated in Figure 4 also does not involve domain-size fitting, and the self-citation to [37] for numerical methods is not load-bearing. The only genuinely circular step is in the large-domain shadowing evidence of Figure 5: the OTR spatial periods are fitted to the measured chaotic pattern sizes, and this fitted agreement is then presented as confirmation that the chaotic trajectory 'transiently shadows the OTR periodic orbits.' Since Supplementary Figure 7 shows the strong sensitivity of OTR profiles to few-percent domain-size changes, the fitted domain choice has substantial leverage over the match. Thus part of the central claim — that OTRs are fundamental building blocks of spatiotemporal chaos — reduces by construction for the large-domain case. Score 6 reflects this partial circularity; it is not higher because the OTR existence, bifurcation analysis, and the small-domain temporal switching remain independent, and because the paper does not rely on a self-citation chain to force its conclusions.

Assumptions & free parameters 1 free parameters · 4 assumptions · 3 invented entities

The central computations rest on the LLE as the governing model, on the validity of the Newton-Krylov/continuation methods for this PDE, on the sufficiency of the Fourier truncation, and on the assumption that stability within a symmetry subspace (L0-stability) is relevant for observability. The only numbers fitted to the chaotic data are the domain-size multipliers used in the shadowing comparison.

free parameters (1)
  • Pattern-domain size multipliers = 0.988 L0,I and 0.970 L0,II
    Used in Fig. 5 to set the spatial period of OTR solutions so that they match the observed chaotic patterns; this is a fit to the data, not a parameter-free prediction.
assumptions (4)
  • domain assumption The Lugiato-Lefever equation (1) correctly models dissipative Kerr cavity dynamics
    The entire claim concerns solutions of this PDE; if the model is not faithful, the results would not transfer to experiments.
  • standard math The flow map F^T of the LLE is smooth and the return map fixed-point condition (2) is solvable by Newton-Krylov methods
    Required for the numerical computation of periodic orbits; standard for globally well-posed dissipative PDEs.
  • domain assumption The truncated Fourier basis (512 modes for L=8π, 128 modes for L0) is sufficient for converged solutions
    No grid-convergence study is presented; the continuation and stability results depend on this resolution.
  • domain assumption L0-stability, i.e. stability in a discrete N-fold translational symmetry subspace, is the relevant stability for observability in the full cavity
    The paper claims L0-stable OTRs could be directly observed; this assumes long-wavelength perturbations that are excluded by the symmetry subspace do not destroy the state in practice.
invented entities (3)
  • OTR I (oscillating Turing roll) solution family independent evidence
    purpose: Anti-phase time-periodic Turing roll states used to explain temporal and spatiotemporal intermittency
    Exact solutions of Eq. (2) with mapped stability regions (Fig. 1); the paper suggests they can be observed in microresonators as periodic power oscillations, giving a falsifiable experimental handle.
  • OTR II relative periodic orbit family
    purpose: Periodic orbits in a moving frame that underlie wavy oscillating roll patterns in Fig. 5(b2)
    Reported at a single parameter set (f^2=16, N=12) with no parameter map or experimental prediction beyond the numerical pattern match.
  • OTR III periodic orbit family
    purpose: Additional periodic orbit involved in temporal switching in Fig. 4
    Introduced to fit an intermittency event; its bifurcation structure is only partially described in the supplement.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Periodic orbits underlying spatiotemporal chaos in the Lugiato-Lefever model." pith.science (2026). https://pith.science/paper/Z2IILGEQ

@misc{pith2026250910283,
  author       = {Pith},
  title        = {Pith review of: Periodic orbits underlying spatiotemporal chaos in the Lugiato-Lefever model},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/Z2IILGEQ}},
  note         = {Machine review of arXiv:2509.10283}
}
read the original abstract

We obtain and investigate theoretically a broad family of stable and unstable time-periodic orbits-oscillating Turing rolls (OTR)-in the Lugiato-Lefever model of optical cavities. Using the dynamical systems tools developed in fluid dynamics, we access the OTR solution branches in parameter space and elucidate their bifurcation structure. By tracking these exact invariant solutions deeply into the chaotic region of the modulation instability, we connect the main dynamical regimes of the Lugiato-Lefever model: continuous waves, Turing rolls, solitons, and breathers, which completes the classical phase diagram of the optical cavity. We then demonstrate that the OTR periodic orbits play a fundamental role as elementary building blocks in the regime of the intracavity field transition from stable Turing rolls to fully developed turbulent regimes. Depending on the cavity size, we observe that the chaotic intracavity field driven by modulation instability displays either spatiotemporal or purely temporal intermittancy between chaotic dynamics and different families of the OTR solutions, exhibiting locally the distinctive wave patterns and large amplitude peaks. This opens avenues for a theoretical description of optical turbulence within the dynamical systems framework.

Figures

Figures reproduced from arXiv: 2509.10283 by the authors.

Figure 1
Figure 1. FIG. 1: Existence and stability range of OTR [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 3
Figure 3. FIG. 3: First row: spatiotemporal portraits of the normalized intensity [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 4
Figure 4. FIG. 4: Spontaneous temporal switching between OTR [PITH_FULL_IMAGE:figures/full_fig_p003_4.png] view at source ↗
Figures from the paper (4 more)
Figure 5
Figure 5. Figure 5: FIG. 5: Spontaneous spatiotemporal emergence of the oscillating roll patterns and the underlying OTR solutions in the regime [PITH_FULL_IMAGE:figures/full_fig_p004_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6: Schematic description of the numerical algorithms for computation and parametric continuation of periodic orbit [PITH_FULL_IMAGE:figures/full_fig_p008_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7: First row: spatiotemporal portraits of the normalized intensity [PITH_FULL_IMAGE:figures/full_fig_p008_7.png]
Figure 8
Figure 8. Figure 8: FIG. 8: (a) Bifurcation diagram in [PITH_FULL_IMAGE:figures/full_fig_p008_8.png]

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

52 extracted references · 2 linked inside Pith

  1. [1]

    L. A. Lugiato and R. Lefever, Physical review letters58, 2209 (1987)

  2. [2]

    Within the whole domainL, theL 0-stable OTRI solu- tions correspond to the OTR solutions that are stable in a discreteN−fold translational symmetry subspace of the intracavity fields. Since the long wavelength instabilities are suppressed within this symmetry subspace, the OTRs in the symmetry subspace lose stability at large detuning compared to theL-sta...

  3. [3]

    Grelu,Nonlinear optical cavity dynamics: from mi- croresonators to fiber lasers(John Wiley & Sons, 2015)

    P. Grelu,Nonlinear optical cavity dynamics: from mi- croresonators to fiber lasers(John Wiley & Sons, 2015)

  4. [4]

    Matsko, A

    A. Matsko, A. Savchenkov, W. Liang, V. Ilchenko, D. Sei- del, and L. Maleki, Optics letters36, 2845 (2011)

  5. [5]

    Barashenkov and Y

    I. Barashenkov and Y. S. Smirnov, Physical Review E 54, 5707 (1996)

  6. [6]

    optical turbulence

    in combination with perturbation theoretical argu- ments valid for small damping and forcing [7–10]. In gen- eral, the nonlinear partial differential equation (1) can, however, only be solved numerically. Depending on the values of the control parameters (ζ0, f), several dynamical regimes have been character- ized and are reliably observed in many optical...

  7. [7]

    Kippenberg, A

    T. Kippenberg, A. Gaeta, M. Lipson, and M. Gorodet- sky, Science361, eaan8083 (2018)

  8. [8]

    V. E. Zakharov and A. B. Shabat, Sov. Phys. JETP34, 62 (1972)

Show all 52 references
  1. [9]

    D. J. Kaup and A. C. Newell, Proceedings of the Royal Society of London. A. Mathematical and Physical Sci- ences361, 413 (1978)

  2. [10]

    Nozaki and N

    K. Nozaki and N. Bekki, Physica D: Nonlinear Phenom- ena21, 381 (1986)

  3. [11]

    Wabnitz, Journal of the Optical Society of America B 13, 2739 (1996)

    S. Wabnitz, Journal of the Optical Society of America B 13, 2739 (1996)

  4. [12]

    Coppini and P

    F. Coppini and P. Santini, Physical Review E102, 062207 (2020)

  5. [13]

    T. Herr, V. Brasch, J. D. Jost, C. Y. Wang, N. M. Kon- dratiev, M. L. Gorodetsky, and T. J. Kippenberg, Nature Photonics8, 145 (2014)

  6. [14]

    F. Leo, S. Coen, P. Kockaert, S.-P. Gorza, P. Emplit, and M. Haelterman, Nature Photonics4, 471 (2010)

  7. [15]

    F. Leo, L. Gelens, P. Emplit, M. Haelterman, and S. Coen, Optics Express21, 9180 (2013)

  8. [16]

    Lucas, M

    E. Lucas, M. Karpov, H. Guo, M. Gorodetsky, and T. J. Kippenberg, Nature communications8, 736 (2017)

  9. [17]

    Godey, I

    C. Godey, I. V. Balakireva, A. Coillet, and Y. K. 6 Chembo, Physical Review A89, 063814 (2014)

  10. [18]

    Parra-Rivas, D

    P. Parra-Rivas, D. Gomila, L. Gelens, and E. Knobloch, Physical Review E98, 042212 (2018)

  11. [19]

    Z. Qi, S. Wang, J. Jaramillo-Villegas, M. Qi, A. M. Weiner, G. D’Aguanno, T. F. Carruthers, and C. R. Menyuk, Optica6, 1220 (2019)

  12. [20]

    Parra-Rivas, D

    P. Parra-Rivas, D. Gomila, L. Gelens, and E. Knobloch, Physical Review E97, 042204 (2018)

  13. [21]

    Haelterman, S

    M. Haelterman, S. Trillo, and S. Wabnitz, Optics com- munications91, 401 (1992)

  14. [22]

    Haelterman, S

    M. Haelterman, S. Trillo, and S. Wabnitz, Optics com- munications93, 343 (1992)

  15. [23]

    T. Herr, K. Hartinger, J. Riemensberger, C. Y. Wang, E. Gavartin, R. Holzwarth, M. L. Gorodetsky, and T. J. Kippenberg, Nature photonics6, 480 (2012)

  16. [24]

    Coillet, Z

    A. Coillet, Z. Qi, I. V. Balakireva, G. Lin, C. R. Menyuk, and Y. K. Chembo, Optics letters44, 3078 (2019)

  17. [25]

    Gomila, P

    D. Gomila, P. Parra-Rivas, P. Colet, A. Coillet, G. Lin, T. Daugey, S. Diallo, J.-M. Merolla, and Y. K. Chembo, Physical Review A106, 053518 (2022)

  18. [26]

    Coillet and Y

    A. Coillet and Y. K. Chembo, Chaos: An Interdisci- plinary Journal of Nonlinear Science24(2014)

  19. [27]

    Panajotov, M

    K. Panajotov, M. G. Clerc, and M. Tlidi, The European Physical Journal D71, 1 (2017)

  20. [28]

    Cheng and X

    G. Cheng and X. Yao, Optics Communications , 132170 (2025)

  21. [29]

    Coulibaly, M

    S. Coulibaly, M. Taki, A. Bendahmane, G. Millot, B. Ki- bler, and M. G. Clerc, Physical Review X9, 011054 (2019)

  22. [30]

    Lukashchuk, J

    A. Lukashchuk, J. Riemensberger, A. Tusnin, J. Liu, and T. J. Kippenberg, Nature Photonics17, 814 (2023)

  23. [31]

    Z. Wang, K. Zhu, F. Zhang, Q. Gong, and Q.-F. Yang, Science Advances11, 9422 (2025)

  24. [32]

    Hopf, Communications on Pure and Applied Mathe- matics1, 303 (1948)

    E. Hopf, Communications on Pure and Applied Mathe- matics1, 303 (1948)

  25. [33]

    Nagata, Journal of Fluid Mechanics217, 519 (1990)

    M. Nagata, Journal of Fluid Mechanics217, 519 (1990)

  26. [34]

    R. R. Kerswell, Nonlinearity18, R17 (2005)

  27. [35]

    J. F. Gibson, J. Halcrow, and P. Cvitanovi´ c, Journal of Fluid Mechanics611, 107 (2008)

  28. [36]

    B. Suri, J. Tithof, R. O. Grigoriev, and M. F. Schatz, Physical review letters118, 114501 (2017)

  29. [37]

    Reetz, T

    F. Reetz, T. Kreilos, and T. M. Schneider, Nature Com- munications10, 2277 (2019)

  30. [38]

    M. D. Graham and D. Floryan, Annu. Rev. Fluid Mech. 53, 227 (2021)

  31. [39]

    Deshmukh, A

    S. Deshmukh, A. Tusnin, A. Tikan, T. J. Kippenberg, and T. M. Schneider, (2025),arXiv:2508.09921

  32. [40]

    Skryabin, Z

    D. Skryabin, Z. Fan, A. Villois, and D. Puzyrev, Physical Review A103, L011502 (2021)

  33. [41]

    D. C. Cole, E. S. Lamb, P. Del’Haye, S. A. Diddams, and S. B. Papp, Nature Photonics11, 671 (2017)

  34. [42]

    Avila, F

    M. Avila, F. Mellibovsky, N. Roland, and B. Hof, Phys. Rev. Lett.110, 224502 (2013)

  35. [43]

    Haelterman, S

    M. Haelterman, S. Trillo, and S. Wabnitz, Physical Re- view A47, 2344 (1993)

  36. [44]

    Villois, D

    A. Villois, D. Puzyrev, D. Skryabin, and M. Onorato, Physical Review A108, 013520 (2023)

  37. [45]

    Tikan, J

    A. Tikan, J. Riemensberger, K. Komagata, S. H¨ onl, M. Churaev, C. Skehan, H. Guo, R. N. Wang, J. Liu, P. Seidler,et al., Nature Physics17, 604 (2021)

  38. [46]

    M. H. Anderson, A. Tikan, A. Tusnin, J. Riemensberger, A. Davydova, R. N. Wang, and T. J. Kippenberg, Phys- ical Review X13, 011040 (2023)

  39. [47]

    Lucas, S.-P

    E. Lucas, S.-P. Yu, T. C. Briles, D. R. Carlson, and S. B. Papp, Nature Photonics17, 943 (2023)

  40. [48]

    A. K. Tusnin, A. M. Tikan, and T. J. Kippenberg, Phys- ical Review A102, 023518 (2020)

  41. [49]

    Englebert, N

    N. Englebert, N. Goldman, M. Erkintalo, N. Mostaan, S.-P. Gorza, F. Leo, and J. Fatome, Nature Physics19, 1014 (2023)

  42. [50]

    Sun, W.-C

    W.-R. Sun, W.-C. Ma, and B. A. Malomed, Physical Review A112, 013501 (2025)

  43. [51]

    Saad and M

    Y. Saad and M. H. Schultz, SIAM Journal on scientific and statistical computing7, 856 (1986)

  44. [52]

    Viswanath, Journal of Fluid Mechanics580, 339 (2007)

    D. Viswanath, Journal of Fluid Mechanics580, 339 (2007). 7 SUPPLEMENT AR Y MA TERIALS NUMERICAL METHODS Newton-Krylov Solver Periodic orbit solutions underlying OTRs, correspond- ing to fixed pointsx ∗ of the return map (2) are com- puted using Jacobian-free Newton-Krylov algo...

Pith tools

Reviewed August 4, 2026 · model on record in the stance chip above.