Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2606.24643

The Quantum Split-Step Fourier Algorithm for Nonlinear Optical Waveguides

Fabio Biancalana

physics.opticsnlin.PSquant-phSubmitted Jun 23, 2026• Announced Jun 24, 2026• 0 citations
arXivHTMLPDF

Abstract

We introduce the Quantum Split-Step Fourier (QSSF) algorithm for nonlinear optical waveguides, a numerical framework that combines split-step propagation of the nonlinear Schrödinger equation with a commutator-preserving Bogoliubov evolution of Gaussian quantum fluctuations. The method propagates the classical mean field together with the Bogoliubov matrices $U$ and $V$, from which reduced second moments, covariance matrices, symplectic eigenvalues, and entropic measures are constructed for arbitrary spectral windows. Applied to soliton-driven resonant radiation, QSSF shows that the selected radiation band acquires a steadily increasing von Neumann entropy and a corresponding loss of purity, quantifying its entanglement with the rest of the spectrum in the lossless Gaussian setting. The analysis also reveals a surprisingly pronounced low-dimensional structure: although the radiation occupies many Fourier bins, its reduced Gaussian state is dominated by only a few Williamson modes. QSSF therefore provides a practical information-theoretic diagnostic for quantum correlations in nonlinear frequency conversion, supercontinuum generation, and multimode squeezed-light formation in ultrafast waveguide platforms.

Related papers

Quantum noise dynamics in nonlinear pulse propagation

Edwin Ng, Ryotatsu Yanagimoto, Marc Jankowski +2 more

The propagation of ultrafast pulses in dispersion-engineered waveguides, exhibiting strong field confinement in both space and time, is a promising avenue towards single-photon nonlinearities in an all-optical platform. However, quantum engineering in such systems requires new numerical tools and physical insights to harness their complicated multimode and nonlinear quantum dynamics. In this work, we use a self-consistent, multimode Gaussian-state model to capture the nonlinear dynamics of broadband quantum fluctuations and correlations, including entanglement. Notably, despite its parametrization by Gaussian states, our model exhibits nonlinear dynamics in both the mean field and the quantum correlations, giving it a marked advantage over conventional linearized treatments of quantum noise, especially for systems exhibiting gain saturation and strong nonlinearities. Numerically, our approach takes the form of a Gaussian split-step Fourier (GSSF) method, naturally generalizing highly efficient SSF methods used in classical ultrafast nonlinear optics; the equations for GSSF evaluate in $O(M^2\log M)$ time for an $M$-mode system with $O(M^2)$ quantum correlations. To demonstrate the broad applicability of GSSF, we numerically study quantum noise dynamics and multimode entanglement in several ultrafast systems, from canonical soliton propagation in third-order ($χ^{(3)}$) waveguides to saturated $χ^{(2)}$ broadband parametric generation and supercontinuum generation, e.g., as recently demonstrated in thin-film lithium niobate nanophotonics.

arXiv:2307.05464•Jul 11, 2023•5 citations

Quantum noise in energy-efficient slow light structures for optical computing: sqeezed light from slow light

Ryan Hamerly, Kambiz Jamshidi, Hideo Mabuchi

Due to their strong light confinement, waveguides with optical nonlinearities may be a promising platform for energy-efficient optical computing. Slow light can enhance a waveguide’s effective nonlinearity, which could result in devices that operate in low-power regimes where quantum fluctuations are important, and may also have quantum applications including squeezing and entanglement generation. In this manuscript, slow-light structures based on the Kerr (χ(3)) nonlinearity are analyzed using a semi-classical model to account for the quantum noise. We develop a hybrid split-step / Runge-Kutta numerical model to compute the mean field and squeezing spectrum for pulses propagating down a waveguide, and use this model to study squeezing produced in optical waveguides. Scaling relations are explored, and the benefits and limitations of slow light are discussed in the context of squeezing.

Proceedings of SPIE, the International Society for Optical Engineering/Proceedings of SPIE•Apr 29, 2016•0 citations

Quantum computing of the nonlinear Schrödinger equation via measurement-induced potential reconstruction

Kaiwen Weng, Zhaoyuan Meng, Guohui Hu

The nonlinear Schrödinger equation (NLSE) is a fundamental model that describes diverse complex phenomena in nature. However, simulating the NLSE on a quantum computer is inherently challenging due to the presence of the nonlinear term. We propose a hybrid quantum-classical framework for simulating the NLSE based on the split-step Fourier method. During the linear propagation step, we apply the kinetic evolution operator to generate an intermediate quantum state. Subsequently, the Hadamard test is employed to measure the Fourier components of low-wavenumber modes, enabling the efficient reconstruction of nonlinear potentials. The phase transformation corresponding to the reconstructed potential is then implemented via a quantum circuit using the phase kickback technique. To validate the efficacy of the proposed algorithm, we numerically simulate the evolution of a Gaussian wave packet, a soliton wave, and the wake flow past a cylinder. The simulation results demonstrate excellent agreement with the corresponding classical solutions. This work provide a concrete basis for analyzing accuracy-cost trade-offs in quantum-classical simulations of nonlinear dispersive wave dynamics.

quant-pharXiv:2601.19184•Jan 28, 2026•0 citations

Quantum computing of the nonlinear Schrödinger equation via measurement-assisted potential reconstruction

Kaiwen Weng, Guohui Hu, Zhaoyuan Meng

The nonlinear Schrödinger equation (NLSE) is a fundamental model that describes diverse complex phenomena in nature. However, simulating the NLSE on a quantum computer is inherently challenging due to the presence of the nonlinear term. We propose a hybrid quantum–classical framework for simulating the NLSE based on the split-step Fourier method. During the linear propagation step, we apply the kinetic evolution operator to generate an intermediate quantum state. The Hadamard test is then used to measure the Fourier components of low-wavenumber modes, from which an approximate nonlinear potential is reconstructed by classical post-processing. The phase transformation corresponding to the reconstructed potential is then implemented as a diagonal phase operation. To validate the proposed method, we numerically simulate the evolution of a Gaussian wave packet, a soliton wave, snake instability, and the wake flow past a cylinder. The simulation results show that the proposed method captures the dominant wave-function profiles, phase structures, and characteristic nonlinear dynamics. This work provides an implementation-level study and an error characterization to assess accuracy-cost trade-offs in hybrid quantum–classical simulations of nonlinear dispersive wave dynamics.

AIP Advances•Jun 1, 2026•0 citations

Phase-space open-systems dynamics of second-order nonlinear interactions with pulsed quantum light

Emanuel Hubenschmid, Victor Rueskov Christiansen

The theoretical description of broadband, multimode quantum pulses undergoing a second-order $χ^{(2)}$-nonlinear interaction can be quite intricate, due to the large dimensionality of the underlying phase space. However, in many cases only a few broadband (temporal) modes are relevant before and after the nonlinear interaction. Here we present an efficient framework to calculate the relation between the quantum states at the input and output of a nonlinear element in their respective relevant modes. Since the number of relevant input and output modes may differ, resulting in an open quantum system, we introduce the generalized Bloch-Messiah decomposition (GBMD), reducing the description to an equal number of input and output modes. The GBMD enables us to calculate the multimode Wigner function of the output state by convolving the rescaled Wigner function of the reduced input quantum pulse with a multivariate Gaussian phase-space function. We expand on this result by considering two examples input states: A Fock state in a single broadband mode and a two-mode squeezed vacuum, both in the THz-frequency regime, up-converted to a single output broadband mode of optical frequencies. We investigate the effect, the convolution and thermalization due to entanglement breakage have on the output Wigner function by calculating the von Neumann entropy of the output Wigner function. The methods presented here can be used to optimize the amplification or frequency conversion of broadband quantum states, opening an avenue to the generation and characterization of optical quantum states on ultrafast time scales.

quant-pharXiv:2512.03933•Physical Review A•Dec 4, 2025•0 citations

Numerical Simulations on Quantum Noise Squeezing for Soliton-like Pulses in Optical Fiber

Arseny A. Sorokin, Nikolay I. Salnikov, Anna N. Kosareva +2 more

Generation of light with nonclassical properties, such as squeezed light, has gathered a lot of attention because of possible uses in such areas as quantum state engineering, quantum imaging, continuous variable quantum computing, and detection of gravitational waves. Squeezed light is generally produced by transporting light through a nonlinear medium. One such media is Kerr-nonlinear optical fibres. We aim to find the parameters for a fibre polarization squeezing setup by studying the dependencies of quantum noise suppression on the duration of the pulse and the input power. This is performed by modeling the light field propagation using the split-step Fourier method to numerically solve the stochastic nonlinear Schrödinger equation obtained by using the Wigner representation. This equation includes such physical effects as damping, quantum loss noise, dispersion, Kerr and Raman nonlinearities, and stochastic Raman noise.

International Conference on Numerical Simulation of Optoelectronic Devices (NUSOD)•Sep 12, 2022•0 citations
We introduce the Quantum Split-Step Fourier (QSSF) algorithm for nonlinear optical waveguides, a numerical framework that combines split-step propagation of the nonlinear Schrödinger equation with a commutator-preserving Bogoliubov evolution of Gaussian quantum fluctuations. The method propagates the classical mean field together with the Bogoliubov matrices UUU and VVV, from which reduced second moments, covariance matrices, symplectic eigenvalues, and entropic measures are constructed for arbitrary spectral windows. Applied to soliton-driven resonant radiation, QSSF shows that the selected radiation band acquires a steadily increasing von Neumann entropy and a corresponding loss of purity, quantifying its entanglement with the rest of the spectrum in the lossless Gaussian setting. The analysis also reveals a surprisingly pronounced low-dimensional structure: although the radiation occupies many Fourier bins, its reduced Gaussian state is dominated by only a few Williamson modes. QSSF therefore provides a practical information-theoretic diagnostic for quantum correlations in nonlinear frequency conversion, supercontinuum generation, and multimode squeezed-light formation in ultrafast waveguide platforms.
The propagation of ultrafast pulses in dispersion-engineered waveguides, exhibiting strong field confinement in both space and time, is a promising avenue towards single-photon nonlinearities in an all-optical platform. However, quantum engineering in such systems requires new numerical tools and physical insights to harness their complicated multimode and nonlinear quantum dynamics. In this work, we use a self-consistent, multimode Gaussian-state model to capture the nonlinear dynamics of broadband quantum fluctuations and correlations, including entanglement. Notably, despite its parametrization by Gaussian states, our model exhibits nonlinear dynamics in both the mean field and the quantum correlations, giving it a marked advantage over conventional linearized treatments of quantum noise, especially for systems exhibiting gain saturation and strong nonlinearities. Numerically, our approach takes the form of a Gaussian split-step Fourier (GSSF) method, naturally generalizing highly efficient SSF methods used in classical ultrafast nonlinear optics; the equations for GSSF evaluate in O(M2log⁡M)O(M^2\log M)O(M2logM) time for an MMM-mode system with O(M2)O(M^2)O(M2) quantum correlations. To demonstrate the broad applicability of GSSF, we numerically study quantum noise dynamics and multimode entanglement in several ultrafast systems, from canonical soliton propagation in third-order (χ(3)χ^{(3)}χ(3)) waveguides to saturated χ(2)χ^{(2)}χ(2) broadband parametric generation and supercontinuum generation, e.g., as recently demonstrated in thin-film lithium niobate nanophotonics.
Due to their strong light confinement, waveguides with optical nonlinearities may be a promising platform for energy-efficient optical computing. Slow light can enhance a waveguide’s effective nonlinearity, which could result in devices that operate in low-power regimes where quantum fluctuations are important, and may also have quantum applications including squeezing and entanglement generation. In this manuscript, slow-light structures based on the Kerr (χ(3)) nonlinearity are analyzed using a semi-classical model to account for the quantum noise. We develop a hybrid split-step / Runge-Kutta numerical model to compute the mean field and squeezing spectrum for pulses propagating down a waveguide, and use this model to study squeezing produced in optical waveguides. Scaling relations are explored, and the benefits and limitations of slow light are discussed in the context of squeezing.
The nonlinear Schrödinger equation (NLSE) is a fundamental model that describes diverse complex phenomena in nature. However, simulating the NLSE on a quantum computer is inherently challenging due to the presence of the nonlinear term. We propose a hybrid quantum-classical framework for simulating the NLSE based on the split-step Fourier method. During the linear propagation step, we apply the kinetic evolution operator to generate an intermediate quantum state. Subsequently, the Hadamard test is employed to measure the Fourier components of low-wavenumber modes, enabling the efficient reconstruction of nonlinear potentials. The phase transformation corresponding to the reconstructed potential is then implemented via a quantum circuit using the phase kickback technique. To validate the efficacy of the proposed algorithm, we numerically simulate the evolution of a Gaussian wave packet, a soliton wave, and the wake flow past a cylinder. The simulation results demonstrate excellent agreement with the corresponding classical solutions. This work provide a concrete basis for analyzing accuracy-cost trade-offs in quantum-classical simulations of nonlinear dispersive wave dynamics.
The nonlinear Schrödinger equation (NLSE) is a fundamental model that describes diverse complex phenomena in nature. However, simulating the NLSE on a quantum computer is inherently challenging due to the presence of the nonlinear term. We propose a hybrid quantum–classical framework for simulating the NLSE based on the split-step Fourier method. During the linear propagation step, we apply the kinetic evolution operator to generate an intermediate quantum state. The Hadamard test is then used to measure the Fourier components of low-wavenumber modes, from which an approximate nonlinear potential is reconstructed by classical post-processing. The phase transformation corresponding to the reconstructed potential is then implemented as a diagonal phase operation. To validate the proposed method, we numerically simulate the evolution of a Gaussian wave packet, a soliton wave, snake instability, and the wake flow past a cylinder. The simulation results show that the proposed method captures the dominant wave-function profiles, phase structures, and characteristic nonlinear dynamics. This work provides an implementation-level study and an error characterization to assess accuracy-cost trade-offs in hybrid quantum–classical simulations of nonlinear dispersive wave dynamics.
The theoretical description of broadband, multimode quantum pulses undergoing a second-order χ(2)χ^{(2)}χ(2)-nonlinear interaction can be quite intricate, due to the large dimensionality of the underlying phase space. However, in many cases only a few broadband (temporal) modes are relevant before and after the nonlinear interaction. Here we present an efficient framework to calculate the relation between the quantum states at the input and output of a nonlinear element in their respective relevant modes. Since the number of relevant input and output modes may differ, resulting in an open quantum system, we introduce the generalized Bloch-Messiah decomposition (GBMD), reducing the description to an equal number of input and output modes. The GBMD enables us to calculate the multimode Wigner function of the output state by convolving the rescaled Wigner function of the reduced input quantum pulse with a multivariate Gaussian phase-space function. We expand on this result by considering two examples input states: A Fock state in a single broadband mode and a two-mode squeezed vacuum, both in the THz-frequency regime, up-converted to a single output broadband mode of optical frequencies. We investigate the effect, the convolution and thermalization due to entanglement breakage have on the output Wigner function by calculating the von Neumann entropy of the output Wigner function. The methods presented here can be used to optimize the amplification or frequency conversion of broadband quantum states, opening an avenue to the generation and characterization of optical quantum states on ultrafast time scales.
Generation of light with nonclassical properties, such as squeezed light, has gathered a lot of attention because of possible uses in such areas as quantum state engineering, quantum imaging, continuous variable quantum computing, and detection of gravitational waves. Squeezed light is generally produced by transporting light through a nonlinear medium. One such media is Kerr-nonlinear optical fibres. We aim to find the parameters for a fibre polarization squeezing setup by studying the dependencies of quantum noise suppression on the duration of the pulse and the input power. This is performed by modeling the light field propagation using the split-step Fourier method to numerically solve the stochastic nonlinear Schrödinger equation obtained by using the Wigner representation. This equation includes such physical effects as damping, quantum loss noise, dispersion, Kerr and Raman nonlinearities, and stochastic Raman noise.