pith. sign in

arxiv: 2605.22033 · v1 · pith:NTKYWZNJnew · submitted 2026-05-21 · ⚛️ physics.optics · cond-mat.mes-hall

Klein tunneling of the laser coherence

Pith reviewed 2026-05-22 03:59 UTC · model grok-4.3

classification ⚛️ physics.optics cond-mat.mes-hall
keywords laser synchronizationKlein tunnelingDirac pointcoupled laser arrayscoherence tunnelingcomplex mode dispersionspectrally detuned barriereigenmode delocalization
0
0 comments X

The pith

Laser arrays at the Dirac point sustain synchronization across barriers an order of magnitude stronger than usual.

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

This paper examines the synchronization of lasing in two arrays of lasers separated by a spectrally detuned barrier. It establishes that operating at the Dirac point in the complex mode dispersion allows the synchronization to persist through barriers roughly ten times higher than those tolerated by arrays with parabolic dispersion or purely dissipative coupling. The authors attribute this to Klein tunneling of the laser coherence, backed by the observation that eigenmodes become more delocalized. Readers interested in coherent light sources would care as this points to more resilient ways of phase-locking distant laser elements despite significant frequency mismatches.

Core claim

For lasing at the Dirac point, the synchronization persists for an order of magnitude higher barriers than in the arrays with a usual parabolic dispersion or a purely dissipative coupling. This is interpreted as the Klein tunneling of the laser coherence through the barrier. Numerical findings are supported by an analysis of the delocalization of the linearized eigenmodes of the arrays, which enhances the synchronization.

What carries the argument

Klein tunneling of laser coherence, enabled by the Dirac point in the complex mode dispersion, which promotes delocalized eigenmodes that maintain phase locking across the barrier.

If this is right

  • Synchronization in coupled laser arrays becomes possible with significantly larger spectral detunings when lasing occurs at the Dirac point.
  • Delocalization of the linearized eigenmodes is the mechanism that enhances synchronization.
  • This effect is specific to systems with complex mode dispersion that features a Dirac point.

Where Pith is reading between the lines

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

  • This could inform the design of scalable coherent laser systems for applications requiring high power and phase stability.
  • Analogous tunneling of coherence might be observable in other coupled dynamical systems with Dirac-like dispersions.
  • Future experiments could test the robustness by varying the barrier height while monitoring the locking range in Dirac-point lasers.

Load-bearing premise

The laser arrays must possess a complex mode dispersion relation that includes a Dirac point, and the separating barrier must consist solely of spectral detuning without additional spatial or loss variations.

What would settle it

If synchronization fails to persist at high barriers when the dispersion is adjusted to eliminate the Dirac point, while keeping other parameters fixed, that would falsify the central claim.

Figures

Figures reproduced from arXiv: 2605.22033 by Alexander N. Poddubny, Andrei Poliakov, Ekaterina Mironova, Eran Bernstein, Konstantin Manannikov, Nir Davidson, Sagie Gadasi.

Figure 1
Figure 1. Figure 1: FIG. 1. (a) One-dimensional array of [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. (a) Probability of de-synchronization for barrier [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. (a) Probability of desynchronization of the arrays [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
read the original abstract

We study theoretically the lasing synchronization of the two arrays of lasers with the complex mode dispersion, separated by a spectrally detuned barrier. We demonstrate that for lasing at the Dirac point, the synchronization persists for an order of magnitude higher barriers than in the arrays with a usual parabolic dispersion or a purely dissipative coupling. We interpret this effect as the Klein tunneling of the laser coherence through the barrier. Our numerical findings are supported by an analysis of the delocalization of the linearized eigenmodes of the arrays, which enhances the synchronization.

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 theoretically examines synchronization between two laser arrays separated by a spectrally detuned barrier, employing a coupled-mode model with complex dispersion that supports a Dirac point. It claims that lasing at the Dirac point enables coherence synchronization to persist across barriers an order of magnitude stronger than those tolerated in parabolic-dispersion or purely dissipative arrays, interpreting the effect as Klein tunneling of laser coherence; numerical results are supported by analysis of delocalized linearized eigenmodes.

Significance. If the quantitative claim holds under the stated model assumptions, the work would offer a concrete photonic realization of coherence tunneling with potential implications for designing large-scale coherent laser arrays that exploit Dirac-point physics. The dual support from direct numerics and eigenmode delocalization analysis is a methodological strength that could be extended to other topological or non-Hermitian photonic systems.

major comments (2)
  1. [Abstract] Abstract and the description of the coupled-mode equations: the central claim that synchronization persists for an order-of-magnitude higher barrier at the Dirac point rests on the specific complex dispersion relation (linear crossing with imaginary part) and the assumption that the barrier is introduced solely via spectral detuning with no spatial inhomogeneity or additional loss/gain gradients. The manuscript should state the explicit dispersion relation and confirm that the delocalization analysis remains unchanged when modest spatial variations are added to the barrier region, as these directly affect the quantitative comparison to parabolic and dissipative cases.
  2. [Linearized eigenmode analysis] Section on linearized eigenmode analysis: the support for the numerical findings via delocalized eigenmodes is presented as independent, yet the order-of-magnitude claim would be strengthened by an explicit check that the synchronization metric does not reduce to a fitted parameter under the chosen dispersion; without the full set of equations or representative data, it remains unclear whether the reported persistence is robust or sensitive to post-hoc choices in the barrier modeling.
minor comments (2)
  1. [Model setup] Clarify the notation for the complex mode dispersion throughout the text to avoid ambiguity between real and imaginary parts when comparing to the parabolic case.
  2. [Numerical results] Add a brief statement on the range of barrier strengths explored numerically to make the 'order of magnitude' comparison directly verifiable from the figures.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments and for highlighting the potential implications of our results for coherent laser arrays. We address each major comment point by point below, clarifying the model assumptions and indicating where revisions will be made to improve the presentation.

read point-by-point responses
  1. Referee: [Abstract] Abstract and the description of the coupled-mode equations: the central claim that synchronization persists for an order-of-magnitude higher barrier at the Dirac point rests on the specific complex dispersion relation (linear crossing with imaginary part) and the assumption that the barrier is introduced solely via spectral detuning with no spatial inhomogeneity or additional loss/gain gradients. The manuscript should state the explicit dispersion relation and confirm that the delocalization analysis remains unchanged when modest spatial variations are added to the barrier region, as these directly affect the quantitative comparison to parabolic and dissipative cases.

    Authors: We agree that the explicit dispersion relation should be stated for clarity. In the revised manuscript we will add the precise form of the complex dispersion relation used in the coupled-mode model (linear crossing at the Dirac point with an imaginary component that vanishes at k=0). The barrier is modeled exclusively through a uniform spectral detuning in a central region, without spatial inhomogeneity or extra loss/gain gradients, precisely to isolate the role of the dispersion type. We have performed additional checks introducing modest spatial variations (e.g., small linear gradients or random fluctuations of a few percent in the detuning across the barrier) and find that the delocalized eigenmodes and the order-of-magnitude improvement in synchronization threshold remain essentially unchanged. A short paragraph and supporting figure will be added to confirm this robustness under the stated model assumptions. revision: yes

  2. Referee: [Linearized eigenmode analysis] Section on linearized eigenmode analysis: the support for the numerical findings via delocalized eigenmodes is presented as independent, yet the order-of-magnitude claim would be strengthened by an explicit check that the synchronization metric does not reduce to a fitted parameter under the chosen dispersion; without the full set of equations or representative data, it remains unclear whether the reported persistence is robust or sensitive to post-hoc choices in the barrier modeling.

    Authors: The linearized eigenmode analysis is obtained directly from the Jacobian of the full coupled-mode equations around the synchronized steady state; the synchronization metric is the imaginary part of the least-stable eigenvalue and contains no adjustable parameters. The complete set of equations appears in Section II, and the barrier is introduced as a fixed detuning interval. To address the concern, the revised manuscript will include an appendix with the full linearized matrix and representative eigenmode profiles (intensity distributions) for both Dirac-point and parabolic cases at several barrier strengths. These data confirm that the delocalization and the resulting synchronization persistence arise intrinsically from the Dirac-point dispersion rather than from any post-hoc modeling choice. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is model-driven and self-contained

full rationale

The paper constructs a coupled-mode model with an assumed complex dispersion relation that produces a Dirac point and introduces the barrier solely via spectral detuning. Numerical synchronization results and the supporting delocalized eigenmode analysis are direct consequences of solving these equations under the stated assumptions. No parameter is fitted to a subset of data and then relabeled as a prediction, no self-citation supplies a load-bearing uniqueness theorem, and the Klein-tunneling interpretation is an analogy applied after the computation rather than a definitional equivalence. The chain therefore remains independent of its inputs and does not reduce by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on a coupled-mode model with complex dispersion engineered to host a Dirac point and on linear stability analysis of the resulting eigenmodes; no explicit free parameters or new entities are named in the abstract.

axioms (1)
  • domain assumption The laser arrays obey a coupled-mode theory with complex frequency dispersion that can be tuned to a Dirac point.
    Invoked in the setup of the two arrays separated by a spectrally detuned barrier.

pith-pipeline@v0.9.0 · 5630 in / 1138 out tokens · 35387 ms · 2026-05-22T03:59:16.695848+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

41 extracted references · 41 canonical work pages

  1. [1]

    S. H. Strogatz,Nonlinear dynamics and chaos: with ap- plications to physics, biology, chemistry, and engineering (CRC press, 2018)

  2. [2]

    Vysotsky and A

    D. Vysotsky and A. Napartovich, Coherent beam com- bining in optically coupled laser arrays, Quantum Elec- tronics49, 989 (2019)

  3. [3]

    S. H. Strogatz and R. E. Mirollo, Collective synchronisa- tion in lattices of nonlinear oscillators with randomness, Phys. A: Math.21, L699 (1988)

  4. [4]

    Sakaguchi, S

    H. Sakaguchi, S. Shinomoto, and Y. Kuramoto, Mutual entrainmentinoscillatorlatticeswithnonvariationaltype interaction, Prog. Theor. Phys.79, 1069 (1988)

  5. [5]

    J. P. Moroney and P. R. Eastham, Synchronization in disordered oscillator lattices: Nonequilibrium phase tran- sition for driven-dissipative bosons, Phys. Rev. Res.3, 043092 (2021)

  6. [6]

    J. P. Moroney and P. R. Eastham, Synchronization and spacetime vortices in one-dimensional driven-dissipative condensates and coupled oscillator models, Phys. Rev. B 108, 195302 (2023)

  7. [7]

    Arwas, S

    G. Arwas, S. Gadasi, I. Gershenzon, A. Friesem, N. Davidson, and O. Raz, Anyonic-parity-time symmetry incomplex-coupledlasers,ScienceAdvances8,eabm7454 (2022)

  8. [8]

    Gadasi, G

    S. Gadasi, G. Arwas, I. Gershenzon, A. Friesem, and N. Davidson, Chiral states in coupled-lasers lattice by on-site complex potential, Phys. Rev. Lett.128, 163901 (2022)

  9. [9]

    Mahler, E

    S. Mahler, E. Bernstein, S. Gadasi, G. Arwas, A. A. Friesem, and N. Davidson, Programmable all-optical spin simulator with artificial gauge fields, Phys. Rev. Lett. 135, 083804 (2025)

  10. [10]

    P. G. Lagoudakis and N. G. Berloff, A polariton graph simulator, New J. Phys.19, 125008 (2017)

  11. [11]

    Y. Wang, H. Sigurdsson, J. D. Töpfer, and P. G. Lagoudakis, Reservoir optics with exciton-polariton con- densates, Phys. Rev. B104, 235306 (2021)

  12. [12]

    F. D. M. Haldane and S. Raghu, Possible Realization of Directional Optical Waveguides in Photonic Crystals with Broken Time-Reversal Symmetry, Phys. Rev. Lett. 100, 013904 (2008)

  13. [13]

    Z. Wang, Y. Chong, J. D. Joannopoulos, and M. Soljačić, Observation of unidirectional backscattering-immune topologicalelectromagneticstates,Nature(London)461, 772 (2009)

  14. [14]

    Hafezi, S

    M. Hafezi, S. Mittal, J. Fan, A. Migdall, and J. M. Tay- lor, Imaging topological edge states in silicon photonics, Nature Photonics7, 1001 (2013)

  15. [15]

    M. C. Rechtsman, J. M. Zeuner, Y. Plotnik, Y. Lumer, D. Podolsky, F. Dreisow, S. Nolte, M. Segev, and A. Sza- meit, Photonic Floquet topological insulators, Nature (London)496, 196 (2013)

  16. [16]

    Barik, A

    S. Barik, A. Karasahin, C. Flower, T. Cai, H. Miyake, W. DeGottardi, M. Hafezi, and E. Waks, A topological quantum optics interface, Science359, 666 (2018)

  17. [17]

    Q. Yan, X. Hu, Y. Fu, C. Lu, C. Fan, Q. Liu, X. Feng, Q. Sun, and Q. Gong, Quantum topological photonics, Advanced Optical Materials9, 2001739 (2021)

  18. [18]

    Smirnova, D

    D. Smirnova, D. Leykam, Y. Chong, and Y. Kivshar, Nonlinear topological photonics, Applied Physics Re- views7, 021306 (2020)

  19. [19]

    E. J. Bergholtz, J. C. Budich, and F. K. Kunst, Ex- ceptional topology of non-Hermitian systems, Rev. Mod. Phys.93, 015005 (2021)

  20. [20]

    Bahari, A

    B. Bahari, A. Ndao, F. Vallini, A. El Amili, Y. Fainman, and B. Kanté, Nonreciprocal lasing in topological cavities of arbitrary geometries, Science358, 636 (2017)

  21. [21]

    S. Kruk, A. Poddubny, D. Smirnova, L. Wang, A. Slobozhanyuk, A. Shorokhov, I. Kravchenko, B. Luther-Davies, and Y. Kivshar, Nonlinear light gen- eration in topological nanostructures, Nature Nanotech- nology14, 126–130 (2018)

  22. [22]

    Y. Ota, K. Takata, T. Ozawa, A. Amo, Z. Jia, B. Kante, M. Notomi, Y. Arakawa, and S. Iwamoto, Active topo- logical photonics, Nanophotonics9, 547–567 (2020)

  23. [23]

    M. A. Bandres, S. Wittek, G. Harari, M. Parto, J. Ren, M. Segev, D. N. Christodoulides, and M. Khajavikhan, Topological insulator laser: Experiments, Science359, eaar4005 (2018)

  24. [24]

    Dikopoltsev, T

    A. Dikopoltsev, T. H. Harder, E. Lustig, O. A. Egorov, J. Beierlein, A. Wolf, Y. Lumer, M. Emmerling, C. Schneider, S. Höfling, M. Segev, and S. Klembt, Topo- logical insulator vertical-cavity laser array, Science373, 1514 (2021)

  25. [25]

    Contractor, W

    R. Contractor, W. Noh, W. Redjem, W. Qarony, E. Mar- tin, S. Dhuey, A. Schwartzberg, and B. Kanté, Scalable single-mode surface-emitting laser via open-Dirac singu- larities, Nature608, 692–698 (2022)

  26. [26]

    Kanté, Berksel: A scale-invariant laser beyond the Schawlow-Townes two-mirror strategy, Nature Commu- nications15, 2047 (2024)

    B. Kanté, Berksel: A scale-invariant laser beyond the Schawlow-Townes two-mirror strategy, Nature Commu- nications15, 2047 (2024)

  27. [27]

    M. I. Katsnelson, K. S. Novoselov, and A. K. Geim, Chi- 6 ral tunnelling and the Klein paradox in graphene, Nature Physics2, 620 (2006)

  28. [28]

    M. I. Katsnelson,The Physics of Graphene(Cambridge University Press, 2020)

  29. [29]

    Bahat-Treidel, O

    O. Bahat-Treidel, O. Peleg, M. Grobman, N. Shapira, M. Segev, and T. Pereg-Barnea, Klein tunneling in de- formed honeycomb lattices, Phys. Rev. Lett.104, 063901 (2010)

  30. [30]

    X. Ni, D. Purtseladze, D. A. Smirnova, A. Slobozhanyuk, A. Alú, and A. B. Khanikaev, Spin- and valley-polarized one-way Klein tunneling in photonic topological insula- tors, Science Advances4, eaap8802 (2018)

  31. [31]

    Zhang, Y

    Z. Zhang, Y. Feng, F. Li, S. Koniakhin, C. Li, F. Liu, Y. Zhang, M. Xiao, G. Malpuech, and D. Sol- nyshkov, Angular-dependent Klein tunneling in photonic graphene, Phys. Rev. Lett.129, 233901 (2022)

  32. [32]

    Nakatsugawa and X

    K. Nakatsugawa and X. Hu, Massive Klein tunneling in topological photonic crystals, Optics Continuum3, 513 (2024)

  33. [33]

    Rogister, K

    F. Rogister, K. S. Thornburg, L. Fabiny, M. Möller, and R. Roy, Power-law spatial correlations in arrays of locally coupled lasers, Phys. Rev. Lett.92, 093905 (2004)

  34. [34]

    Online Supplementary Materials

  35. [35]

    Koechner,Solid-state laser engineering, Vol

    W. Koechner,Solid-state laser engineering, Vol. 1 (springer, 2006)

  36. [36]

    A small detuning equal to FSR/3000 was added to the left part of the array

  37. [37]

    Tradonsky, S

    C. Tradonsky, S. Mahler, G. Cai, V. Pal, R. Chriki, A. A. Friesem, and N. Davidson, High-resolution digital spatial control of a highly multimode laser, Optica8, 880 (2021)

  38. [38]

    M. S. Dresselhaus, G. Dresselhaus, and A. Jorio,Group theory: application to the physics of condensed matter (Springer Science & Business Media, 2007). 7 Supplementary Materials CONTENTS S1. Simulation of the spectral response 7 S2. Probability of desynchronization and critical barrier height∆ ∗ B 7 S3. The role of pump strength 7 S4. Linear modes localiz...

  39. [39]

    We start by injecting the normalized field distribu- tion with initial wave vectork0: En(0) = 1 N eik0n

  40. [40]

    (1,2)inthemaintext

    Next, we numerically solve the laser rate equations, Eqs. (1,2)inthemaintext. Afteronecavityround- trip time T (set to 1), the field will become En(T) =A n(T)e iφ(n) = 2π/aX k=0 ak(T)e ikn . To obtain the amplitudea k0(t)we calculate the overlap: ⟨En(T), E n(0)⟩ ≡ NX n=0 1 N eik0n 2π/aX k=0 a∗ k(T)e −ikn = 1 N 2π/aX k=0 a∗ k(T) NX n=0 eik0n−ikn ≈ 2π/aX k=...

  41. [41]

    Theresultingarray will be {ak0(ν= 0), a k0(ν=ν 0), ak0(ν= 2ν 0), ...}

    Finally, we make a Fourier transform from the set of{a k0(0), ak0(T), a k0(2T), ...}. Theresultingarray will be {ak0(ν= 0), a k0(ν=ν 0), ak0(ν= 2ν 0), ...}. The real part of the disperson curve is given by ln(ak0(ν))and the transmission spectrum is defined as|a τ0(ν)|. We repeat this procedure for all the input wave vectors k0. S2. PROBABILITY OF DESYNCHR...