The Dirac sea of phase: Unifying phase paradoxes and Talbot revivals in multimode waveguides

Abstract: The quantum mechanical description of phase remains a fundamental challenge, with theoretical efforts tracing from the early works of London and Dirac to discrete formalisms. In this work, we extend the action-angle formalism to the Helmholtz-Schrödinger equation by introducing a phase-dependent wavefunction $φ(θ, t)$ residing in the Hardy space $H<sup>2(\mathbb{D})$. This mathematical structure, defined by functions analytic on the unit disk with square-integrable boundary values, naturally ensures the positivity of the energy spectrum while providing a rigorous framework for wave dynamics in photonic systems. We demonstrate that establishing a self-adjoint phase operator requires extending the Hilbert space to $L<sup>2$, a procedure that necessitates the admission of negative energy states. We interpret these states through an analogy with the Dirac sea, where the existence of antiphase or antiphoton modes provides a conceptual framework for understanding the fundamental limits of phase localization and quantum uncertainty. This formalism is applied to light propagation in multimode waveguides characterized by anharmonic refractive index profiles. By mapping modal dispersion to our phase representation, we show that the deviation of propagation constants from linear spacing governs the spatial evolution of the optical field. This approach offers a clear mechanism for the emergence of periodic self-imaging known as the Talbot effect, the generation of fractional revivals, and the formation of complex fractal interference patterns, providing a robust toolkit for the characterization and design of multimode interference devices.

• The paper develops a Hardy-space phase representation and extends it to an L² “Dirac sea” to address phase-operator self-adjointness and negative-frequency modes.
• The paper models quartically perturbed waveguides numerically, showing that the dispersive coefficient a₂ = 3λ/2 controls revival timing while higher-order spectral terms reduce later-revival fidelity.
• The paper links phase-space shredding, fractional revivals, and fractal Talbot carpets to multimode interference, offering design parameters for self-imaging optical couplers and classical analogues of quantum collapse and revival.

Overview and motivation

This paper by Korneev, Ramos-Prieto, and Moya-Cessa addresses two problems that are usually treated separately: the mathematical definition of quantum phase, and the dynamics of light propagation in multimode waveguides with anharmonic modal spectra. The unifying observation is that paraxial propagation in a graded-index waveguide is governed by a Helmholtz-Schrödinger equation, so that the longitudinal coordinate $z$ plays the role of time. Within this isomorphism, the authors construct a continuous phase representation for the optical field in the Hardy space $H^2(\mathbb{D})$, extend it to $L^2$ via a Dirac-sea-like argument, and then apply the resulting machinery to predict Talbot revivals, fractional revivals, and fractal Talbot carpets in quartically perturbed waveguides (2603.07872).

The phase representation in Hardy space

The central construction is a phase-dependent wavefunction defined by a one-sided Fourier series,

$$\phi(\theta, t) = \frac{1}{\sqrt{2\pi}} \sum_{n=0}^{\infty} C_n(t) e^{in\theta},$$

with $\theta \in [0, 2\pi)$. Because only non-negative frequencies appear, $\phi$ is the boundary value of a function analytic in the open unit disk, i.e., an element of $H^2(\mathbb{D})$. The number operator maps to $-i\,\partial/\partial\theta$, and for the harmonic Hamiltonian $\hat{\mathcal{H}} = \hat{n}$ the phase evolution is a rigid transport equation $i\,\partial_t \phi = -i\,\partial_\theta \phi$.

The Hardy space formulation enforces positivity of the photon-number spectrum geometrically, through the analyticity of the domain, rather than through the truncation or limiting procedures of the Susskind-Glogower and Pegg-Barnett approaches (2603.07872). A notable consequence is that the real and imaginary parts of $H^2(\mathbb{D})$0 are locked by a Hilbert transform: the full state is encoded in either component alone, in a manner the authors compare to Kramers-Kronig relations. This means the real transport equations for $H^2(\mathbb{D})$1 and $H^2(\mathbb{D})$2 are individually equivalent to the full quantum evolution. The authors also invoke the rigged-Hilbert-space program of Bohm and collaborators, in which Hardy spaces encode time-asymmetry and a semigroup (rather than group) evolution.

The Dirac sea of phase

A self-adjoint phase (angle) operator cannot be defined within $H^2(\mathbb{D})$3: multiplication by a real function of $H^2(\mathbb{D})$4 generates negative-frequency components that leave the physical subspace, and while $H^2(\mathbb{D})$5 is admissible, its inverse is not, precluding unitary phase shifts. The authors therefore extend the Hilbert space to the full $H^2(\mathbb{D})$6, where both $H^2(\mathbb{D})$7 and the multiplication operator $H^2(\mathbb{D})$8 are well defined. The cost is that the spectrum of $H^2(\mathbb{D})$9 becomes the full set of integers, admitting negative-energy states.

The paper interprets these negative modes through an explicit analogy with Dirac's hole theory: the $L^2$0 extension constitutes a "Dirac sea of phase," with holes in the negative-energy sea manifesting as positive-energy states with inverted phase evolution—antiphoton modes. The authors state that these virtual modes provide a dynamical basis for the number-phase uncertainty relation, preventing infinite phase localization. It should be noted that this interpretation is presented as a heuristic; the paper does not derive observable consequences of antiphoton modes, and the connection to the uncertainty principle is asserted rather than quantified within the formalism.

Anharmonic waveguides and spectral analysis

The second part of the paper applies the formalism to waveguides whose refractive index profile deviates from the ideal parabola. The effective Hamiltonian is taken as a harmonic oscillator plus a quartic perturbation, $L^2$1. Parity conservation restricts the quartic term to couple Fock states of the same parity ($L^2$2), giving a pentadiagonal matrix in the number basis. Numerical diagonalization yields the anharmonic eigenvalues $L^2$3 and eigenvectors $L^2$4, which can be projected into both the position basis (Hermite-Gauss expansion) and the phase basis. The spectral deformation is non-uniform: higher-index modes shift more strongly, which is the direct mechanism for dephasing and rephasing during propagation.

Propagation dynamics and the Talbot mechanism

Using the numerically exact anharmonic basis, the evolved state is expanded as $L^2$5, with coherent-state input amplitudes $L^2$6 (simulations use $L^2$7). For weak anharmonicity, the Hamiltonian reduces to

$L^2$8

which, projected onto the phase basis, yields an effective dispersive equation

$L^2$9

with $$\phi(\theta, t) = \frac{1}{\sqrt{2\pi}} \sum_{n=0}^{\infty} C_n(t) e^{in\theta},$$0 and $$\phi(\theta, t) = \frac{1}{\sqrt{2\pi}} \sum_{n=0}^{\infty} C_n(t) e^{in\theta},$$1. The first-order term produces rigid rotation; the second produces phase diffusion, which is reversible because of the periodic boundary condition and the discrete spectrum. Revivals occur at $$\phi(\theta, t) = \frac{1}{\sqrt{2\pi}} \sum_{n=0}^{\infty} C_n(t) e^{in\theta},$$2: doubling the anharmonicity from $$\phi(\theta, t) = \frac{1}{\sqrt{2\pi}} \sum_{n=0}^{\infty} C_n(t) e^{in\theta},$$3 to $$\phi(\theta, t) = \frac{1}{\sqrt{2\pi}} \sum_{n=0}^{\infty} C_n(t) e^{in\theta},$$4 halves the revival period (from roughly 300–400 to about half that value in the simulations), at the cost of degraded fidelity in subsequent revivals due to higher-order spectral terms that break exact commensurability.

The simulations identify three propagation regimes: near-field coherent transport, a mid-field fractal regime with high-density Talbot carpets, and a self-imaging window where the field reconstructs. The phase-space representation $$\phi(\theta, t) = \frac{1}{\sqrt{2\pi}} \sum_{n=0}^{\infty} C_n(t) e^{in\theta},$$5, plotted in polar coordinates, exposes the analytic structure directly—the linear spectrum shreds the phase distribution while the compactness of the angular domain guarantees refocusing. The collapse-revival dynamics of the beam centroid $$\phi(\theta, t) = \frac{1}{\sqrt{2\pi}} \sum_{n=0}^{\infty} C_n(t) e^{in\theta},$$6 mirror the Jaynes-Cummings collapse and revival phenomenology, with spectral nonlinearity playing the role of the atom-cavity detuning. The authors argue that multimode interference devices can therefore serve as classical analogues for quantum coherence phenomena, and that the coefficients $$\phi(\theta, t) = \frac{1}{\sqrt{2\pi}} \sum_{n=0}^{\infty} C_n(t) e^{in\theta},$$7, $$\phi(\theta, t) = \frac{1}{\sqrt{2\pi}} \sum_{n=0}^{\infty} C_n(t) e^{in\theta},$$8 provide direct design parameters for self-imaging couplers beyond the parabolic approximation.

Limitations and open questions

Several caveats bear on the results. The Dirac sea of phase is a conceptual analogy; the paper offers no operational prescription for detecting or exploiting antiphoton modes, and the claimed link to phase localization limits is qualitative. The revival analysis relies on a weak-anharmonicity expansion truncated at $$\phi(\theta, t) = \frac{1}{\sqrt{2\pi}} \sum_{n=0}^{\infty} C_n(t) e^{in\theta},$$9, and the observed degradation of later revivals confirms that higher-order terms matter quantitatively. All propagation results are numerical and assume lossless, time-independent, Kerr-free waveguides; fabrication disorder, material chromaticity, and nonlinear effects are cited as sources of anharmonicity but are not modeled explicitly. The authors themselves flag non-Hermitian and parity-time-symmetric waveguides—where gain and loss could modulate the dephasing rate and revival fidelity—as an unexplored extension, and the connection to philophase minimum-uncertainty states in weighted Bergman spaces is invoked rather than developed.

Conclusion

The paper provides a coherent framework in which the spectral positivity of the mode index is a geometric property of the Hardy space, the self-adjointness of phase requires an $\theta \in [0, 2\pi)$0 extension interpretable as a phase Dirac sea, and the deviation of modal spacing from linearity—captured by the dispersive coefficient $\theta \in [0, 2\pi)$1—quantitatively controls Talbot revival periods and carpet formation in multimode waveguides. The main open questions left by the work concern the physical observability of the antiphoton sector, the quantitative role of higher-order dispersion on revival fidelity, and the behavior of the formalism in non-Hermitian waveguide geometries.