General method for solving nonlinear optical scattering problems using fix point iterations

This paper introduces a new fixed-point iteration scheme based on the Bidirectional Pulse Propagation Equations to solve nonlinear electromagnetic scattering problems in slab-like geometries with arbitrary material responses, demonstrating its convergence and accuracy through a case study involving mixed fast electronic and slow molecular vibrational responses.

The Gist

Imagine you are trying to predict how a beam of light behaves when it hits a special piece of glass. This isn't just any glass; it's "smart" glass. When light hits it, the glass doesn't just sit there; it changes its own properties based on how bright the light is. This is called nonlinear optics.

The problem is that predicting this behavior is incredibly hard. It's like trying to predict the path of a ball bouncing through a room full of trampolines that change their bounciness the moment the ball touches them.

Here is a breakdown of what Per Kristen Jakobsen's paper does, using simple analogies.

1. The Old Ways vs. The New Way

Scientists have tried to solve this problem for years using two main methods:

• The "Time-Travel" Method (FDTD): Imagine trying to simulate a movie frame-by-frame. You calculate what happens at 1 second, then 2 seconds, then 3. This works well for simple materials, but for "smart" glass that reacts in complex ways, you need to know the future to calculate the present. It's like trying to drive a car while only looking at the rearview mirror. It gets computationally expensive and slow.
• The "Frequency" Method (UPPE): Instead of looking at time, this method looks at the "colors" (frequencies) of the light. It's like analyzing a song by looking at its sheet music rather than listening to it play. This is fast and handles complex materials well. However, it has a fatal flaw: it assumes the light only moves forward. If the light hits a wall and bounces back (reflects), this method gets confused because it doesn't know the "future" reflection is coming back to hit it.

The Paper's Solution:
The author introduces a new method called Bidirectional Pulse Propagation Equations (BPPE). Think of this as a two-way street. It tracks light moving forward and light bouncing backward simultaneously.

2. The "Guess and Check" Strategy (Fixed-Point Iteration)

The core of the paper is a new mathematical trick to solve the equations for this two-way street.

Imagine you are trying to find the perfect temperature for a shower, but the faucet is broken. You can't just turn the knob to the right spot. Instead, you have to:

1. Guess a temperature.
2. Turn the knob.
3. Feel the water.
4. Adjust your guess based on how hot or cold it feels.
5. Repeat until the water feels "just right."

In math, this is called Fixed-Point Iteration.

• The Problem: The light hitting the back of the glass creates a reflection that travels back to the front. But the front is where the light started. So, the reflection depends on the light, but the light depends on the reflection. It's a "chicken and egg" problem.
• The Trick: The author rewrites the problem so you start with a "linear" guess (assuming the glass is boring and doesn't change). Then, you use a mathematical "map" to see how the "smart" glass would actually react to that guess. You feed the result back into the map, get a new guess, and repeat.
• The Result: The author shows that for most materials, this "guess and check" process converges very quickly. You don't need super-complex math; you just need to keep iterating the map, and it naturally settles on the correct answer.

3. The "Ghost" Problem (Causality)

One of the most interesting parts of the paper is a "scare" the author had.

When they ran the simulation, the computer produced a picture of the light bouncing inside the glass. At first glance, it looked like time travel.

• It looked like a pulse of light was appearing at the back of the glass and traveling backward in time to the front.
• This is impossible (it violates causality—the rule that cause must come before effect).

The Resolution:
The author realized the computer wasn't broken; the interpretation was wrong.

• The Analogy: Imagine you are looking at a shadow on a wall. The shadow looks like a monster, but it's actually just a hand holding a puppet.
• The "backward traveling" pulse wasn't actually moving backward in time. It was a mathematical artifact. The "left-moving" wave the computer calculated was actually a mix of a real backward wave and a weirdly shaped forward wave. When you separate them correctly, everything makes sense: the light moves forward, hits the back, bounces back, and moves forward again. No time travel required.

4. Why This Matters

This paper is a "how-to" guide for a new, efficient way to simulate light.

• Speed: It's much faster than the old "frame-by-frame" methods.
• Accuracy: It handles complex materials (like the glass that changes with light intensity) and reflections better than the old "frequency-only" methods.
• Versatility: While the paper tests it on a simple slab of glass, the math can be applied to any shape of material where light bounces around, including acoustic waves (sound) and ocean waves.

Summary

The author built a new mathematical engine to simulate how light interacts with "smart" materials. They solved a tricky "chicken and egg" problem by using a simple "guess and check" loop that converges quickly. They also solved a mystery where the simulation looked like it was breaking the laws of physics (time travel), proving that the laws of physics were safe all along—it was just a misunderstanding of the data.

It's a significant step forward in our ability to design better lasers, fiber optics, and optical computers.

Technical Summary

1. Problem Statement

The paper addresses the computational challenge of solving nonlinear electromagnetic scattering problems, specifically for pulse propagation in materials with arbitrary temporal responses (both linear dispersion and nonlinear effects).

• Limitations of Existing Methods:
  • Time-Domain Methods (e.g., FDTD): Struggle with materials having complex spectral responses or broad-band pulses because they require modeling non-local time responses, which is computationally expensive and difficult to scale for macroscopic materials.
  • Unidirectional Pulse Propagation Equation (UPPE): Solves problems in the spectral domain and handles complex material responses efficiently. However, it assumes no backward-propagating waves (reflections). If reflections occur (due to interfaces or nonlinearities), UPPE fails because it requires knowledge of the full temporal spectrum at a specific point $z=a$ for all time (past, present, and future), which is physically impossible to know a priori when reflections return to that point.
  • Bidirectional Pulse Propagation Equation (BPPE): A spectral method that accounts for both forward and backward moving waves. However, it leads to a coupled system where the backward component at the input ($z=a$) is unknown because it is generated by interactions deeper in the material.

The Core Problem: How to solve the coupled BPPE equations for a slab geometry where the backward-propagating spectral amplitude at the input is unknown and depends on the nonlinear dynamics within the slab.

2. Methodology

The author proposes a novel reformulation of the BPPE scattering problem into a nonlinear mapping problem solvable via fixed-point iteration.

A. Mathematical Formulation

1. Spectral Decomposition: The electric ($E$) and magnetic ($H$) fields are expressed as superpositions of right-moving ($A^+$) and left-moving ($A^-$) spectral modes within three regions: before the slab ($z<a$), inside the slab ($a<z<b$), and after the slab ($z>b$).
2. BPPE Equations: The propagation of these spectral amplitudes inside the nonlinear slab is governed by first-order differential equations driven by the nonlinear polarization $\hat{p}(z, \omega)$.
3. Mapping to Boundary Conditions:
   • The problem is reduced to finding the unknown reflected spectral amplitude $A^-(a, \omega)$ at the input interface.
   • A nonlinear map $G$ is defined such that $G(A^-(a, \omega)) = S_R(\omega)$, where $S_R$ is the known source at the right boundary.
   • $G$ is constructed by composing:
     • The boundary condition matrix at the input ($M_{12}$).
     • The propagation operator through the slab ($U(a, b)$).
     • The boundary condition matrix at the output ($M_{23}$).
     • A projection operator ($P$) that extracts the backward component.

B. Fixed-Point Iteration Scheme

Directly solving $G(u) = S_R$ is computationally expensive in high dimensions. The author exploits the fact that in most optical applications, nonlinear effects are weak ($\epsilon \ll 1$).

1. Linearization: The map $U(a, b)$ is decomposed into the Identity map ($I$) plus a deviation $V$ (where $V$ is proportional to the nonlinearity).
2. Deviation Variable: Let $A_0$ be the solution to the linear scattering problem. The unknown is redefined as the deviation $a^- = A^-(a) - A_0$.
3. Iterative Map: The equation is transformed into a fixed-point problem:
   $$a^- = H(a^-)$$
   where $H$ is a nonlinear map derived from the deviation $V$.
4. Algorithm: The solution is found by simple iteration: $a^-_{n+1} = H(a^-_n)$. This is computationally cheaper than quasi-Newton methods and avoids the need for Jacobian calculations.

C. Validation via Exact Solution Construction

To verify that the iteration converges to the correct physical solution (and not a spurious one), the author develops a method in Appendix B to generate exact solutions for the nonlinear mapping problem without iteration. This involves:

• Solving the BPPE equations backwards from the output with a free parameter.
• Reflecting and flipping the solution to define a new input.
• Using this constructed "exact" solution to benchmark the convergence of the fixed-point iteration.

3. Key Contributions

1. Reduction to Fixed-Point Problem: The paper successfully reformulates the complex bidirectional scattering problem into a fixed-point iteration for the deviation from the linear solution. This avoids the "future knowledge" paradox of UPPE while handling reflections.
2. Computational Efficiency: By reducing the problem to simple function iteration (rather than solving large nonlinear systems via Newton-Raphson), the method becomes significantly more computationally efficient per step, making it feasible for high-dimensional spectral grids.
3. Causality Resolution: The paper addresses a critical issue where the numerical solution appeared to violate causality (showing "acausal" backward waves). The author proves that this is an artifact of interpretation. The spectral amplitude $A^-(z, \omega)$ (labeled "left-moving") mathematically contains components of both left- and right-moving physical fields. When reconstructed in the time domain, the resulting fields are strictly causal.
4. Generalizability: The method is not limited to optical pulses; it applies to any wave scattering problem (acoustic, ocean waves) where the geometry is slab-like and the material response is arbitrary.

4. Results

The method was tested on a purely nonlinear slab with a polarization response consisting of:

• A fast electronic (Kerr-like) response.
• A slow molecular (Raman-like) response.

Simulation Findings:

• Convergence: For a slab of 50 wavelengths, the fixed-point iteration converged to a solution satisfying the scattering equation to 10 digits of accuracy within 30 iterations. For a 150-wavelength slab, it reached 12 digits of accuracy within 40 iterations.
• Causality Check: Space-time plots of the electric field initially showed "acausal" features (waves appearing to travel backward in time). However, analysis of a linearized test case proved that these features are actually right-moving waves misinterpreted as left-moving due to the spectral decomposition. The final physical fields are strictly causal.
• Exact Solution Verification: Using the exact solution construction method, the authors confirmed that the fixed-point iteration converges to the unique physical solution of the scattering problem.

5. Significance

• Bridging the Gap: This method fills a critical gap between time-domain methods (limited by dispersion complexity) and unidirectional spectral methods (limited by reflections). It allows for the simulation of strong reflections in broadband, dispersive, nonlinear media.
• Handling Complex Materials: It can handle arbitrary material responses (sums of fast and slow vibrational modes) without simplifying assumptions about the spectral width of the pulse.
• Robustness: The demonstration that the method converges to a causal, physically meaningful solution, even when the spectral decomposition suggests otherwise, validates the robustness of the approach for practical engineering and physics applications.
• Future Applications: The framework is presented as a general tool for wave scattering, potentially applicable to acoustic and oceanographic problems, and offers a pathway to generalize to compound slabs and complex transverse geometries.

In summary, Jakobsen presents a mathematically rigorous and computationally efficient framework for solving nonlinear optical scattering with reflections, resolving previous causality ambiguities and providing a verified numerical scheme that outperforms standard approaches in specific high-complexity regimes.