Simulating multimode nonlinearities at very high mode count via parallel factor decomposition
Abstract
A variety of problems at the frontier of nonlinear optics involve strongly multimode interactions. Simulating these dynamics is often computationally challenging due to the complexity of evaluating a large number of nonlinear overlaps. This is especially acute in ultrafast Kerr effect physics, where simulations in a modal representation have a generic complexity scaling quartically with the number of modes. This has limited the largest number of modes retainable in simulations in practice to a few tens of modes (often fewer), even when run on modern GPUs. Here, I show how by taking advantage of parallel factor decompositions of the nonlinear overlap tensor, a far larger number of modes can be accommodated. Realistic values for generic index profiles without special symmetries or separability can be in the thousands, and with special symmetries, far larger. I show that simulations with >1,000 modes can be done on modern hardware, illustrating this with the example of a multimode soliton propagating over a meter of GRIN fiber. I further show that the improved scaling immediately makes it practical to consider statistical effects, illustrating with the example of noise amplification in soliton fission. These results, implemented in public code with examples, should enable access to simulating the most demanding regimes of multimode nonlinear optics. Compression of the nonlinear interaction tensor may provide a route for more efficient simulation of nonlinear wave dynamics across a variety of other systems including acoustics, water waves, phononics, magnonics, and Bose condensates.
Related
- Thermodynamic origins of multimode spinodal decomposition in multicomponent alloys
- The exact universal exchange-only local field factor of the uniform electron gas and its near-cancellation by parallel-spin correlation near 2kF
- Quantum Vacuum Nonlinearities in Laser Interferometers
- Second- and third-order nonlinearities in silicon nitride waveguides via heterogeneous integration of lithium niobate membranes
- Synchronously-pumped pulsed multimode Raman fiber laser
BibTeX
@misc{rivera2026simulating,
title = {{Simulating multimode nonlinearities at very high mode count via parallel factor decomposition}},
author = {Nicholas Rivera},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.16502},
eprint = {2609.16502},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.16502},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Simulating multimode nonlinearities at very high mode count via parallel factor decomposition
AU - Nicholas Rivera
PY - 2026
DA - 2026/09/15
PB - arXiv
DO - 10.48550/arxiv.2609.16502
UR - https://arxiv.org/abs/2609.16502
AB - A variety of problems at the frontier of nonlinear optics involve strongly multimode interactions. Simulating these dynamics is often computationally challenging due to the complexity of evaluating a large number of nonlinear overlaps. This is especially acute in ultrafast Kerr effect physics, where simulations in a modal representation have a generic complexity scaling quartically with the number of modes. This has limited the largest number of modes retainable in simulations in practice to a few tens of modes (often fewer), even when run on modern GPUs. Here, I show how by taking advantage of parallel factor decompositions of the nonlinear overlap tensor, a far larger number of modes can be accommodated. Realistic values for generic index profiles without special symmetries or separability can be in the thousands, and with special symmetries, far larger. I show that simulations with >1,000 modes can be done on modern hardware, illustrating this with the example of a multimode soliton propagating over a meter of GRIN fiber. I further show that the improved scaling immediately makes it practical to consider statistical effects, illustrating with the example of noise amplification in soliton fission. These results, implemented in public code with examples, should enable access to simulating the most demanding regimes of multimode nonlinear optics. Compression of the nonlinear interaction tensor may provide a route for more efficient simulation of nonlinear wave dynamics across a variety of other systems including acoustics, water waves, phononics, magnonics, and Bose condensates.
ER -