RP Photonics logo
RP Photonics
Modeling & Design Software
Technical consulting services on lasers, nonlinear optics, fiber optics etc.
Profit from the knowledge and experience of a top expert!
Powerful simulation and design software.
Make computer models in order to get a comprehensive understanding of your devices!
Success comes from understanding – be it in science or in industrial development.
The famous Encyclopedia of Laser Physics and Technology – available online for free!
The ideal place for finding suppliers for many photonics products.
Advertisers: Make sure to have your products displayed here!
powerful tools for efficient laser development and laser science
This page gives you an overview on our software products.
RP Fiber Calculator is a convenient tool for calculations on optical fibers.
RP Fiber Power is an extremely flexible tool for designing and optimizing fiber devices.
RP Resonator is a particularly flexible tool for laser resonator design.
RP ProPulse can simulate the pulse evolution e.g. in mode-locked lasers and sync-pumped OPOs.
RP Coating is a particularly flexible design tool for dielectric multilayer systems.
RP Q-switch can simulate the power evolution in Q-switched lasers.
Most of our software products support a powerful script language, which gives you an extraordinary degree of flexibility.
Here you learn about software license conditions, updates and upgrades, etc.
Competent technical support is a key quality associated with software from RP Photonics.
RP Photonics has distributors in various countries.
The RP Photonics Software News keep you updated on our developments and give the users additional interesting hints.
Here you can make inquiries e.g. concerning technical details, prices and quotations.
en | de

RP Fiber Power – Simulation and Design Software
for Fiber Optics, Amplifiers and Fiber Lasers

Example Case: Nonlinear Self-focusing in a Fiber

Description of the Model

Here, we investigate details of nonlinear self-focusing in a fiber. First, we calculate how the fundamental mode of a large mode area fiber shrinks as the effect of nonlinear self-focusing.

The mode solver actually ignores nonlinear effects. However, with a few lines of script code we can store the refractive index profile including its nonlinear changes and then recalculate the fiber modes. This is repeated until we get a self-consistent solution:

dr := 0.05 um
defarray I[0,200 um,dr]
n_f_nl(r) := n_f(r) + n2 * (if r <= r_max then I~[r])
  { nonlinear refractive index profile }
store_I(P) := 
  for r := 0 to 2 * r_core step dr do I[r] := P * I_lm(0,1,lambda,r)
  { ignore index changes outside 2 * r_core, where the intensity is small }
CalcNonlinearMode(P) := 
  { Calculate the lowest-order mode with self-focusing for the power P. }
    var A, A_l;
    A := 0;
      A_l := A;
      A := A_eff_lm(0,1,lambda);
    until abs(A_l / A - 1) < 1e-6;

We can also numerically simulate the beam propagation, considering the fiber nonlinearity. For that, we need to define a numerical grid and set various other inputs for the beam propagation:

x_max := 30 um { maximum x or y value }
N := 2^5 { number of grid points in x and y direction }
dx := 2 * x_max / N { transverse resolution }
z_max := 30 mm { fiber length }
dz := 100 um { longitudinal resolution }
N_z := z_max / dz { number of z steps }
N_s := 100 { number of sub-steps per dz step }
P_11 := 4 MW
A0%(x,y) := sqrt(P_11) * A_lm_xy(1, 1, lambda, x, y)  { initial field }
    bp_set_grid(x_max, N, x_max, N, z_max, N_z, N_s);
    bp_set_n('n_f(sqrt(x^2 + y^2))'); { index profile }
    bp_set_loss('10e2 * ((x^2 + y^2) / (20 um)^2)^3');  { simulate loss for cladding modes }
    bp_set_A0('A0%(x,y)'); { initial amplitude }
    bp_set_interpol(2); { quadratic interpolation }


Figure 1 shows the mode profile for an optical power of 5 MW (not far below the power for catastrophic self-focusing), and the corresponding refractive index profile.

mode profile with strong self-focusing

Figure 1: Calculated normalized mode intensity profiles with and without self-focusing. The refractive index profiles are also shown. One sees that the index profile is substantially modified by the nonlinear effect.

Figure 2 shows the mode area as a function of the optical power. The mode area shrinks dramatically as the critical power is approached.

mode area vs. optical power

Figure 2: Mode area versus optical power. The red line indicates the critical power for catastrophic self-focusing.

Figure 3 shows the maximum power as a function of the core radius. For each core radius, one has to calculate the optical power for which the on-axis intensity reaches the damage threshold. Of course, the modes need to be recalculated for each power value.

maximum power in fiber vs. core radius

Figure 3: Maximum optical power in the fiber as a function of the core radius.

Initially, the maximum power scales with the core area. For large cores, however, the rise becomes substantially slower, as the mode area is reduced by self-focusing.

Now we investigate what happens if we inject light into the LP11 mode (the first higher-order mode) of the fiber, as calculated without the nonlinearity. For that, we can use numerical beam propagation. Figure 4 shows the result if we inject an optical power of 4 MW, not far below the critical power for self-focusing. Here, the higher-order mode becomes unstable. Even the slightest asymmetry, here resulting from small numerical errors, causes that mode to be transferred into a superposition of the LP01 mode and the LP11 mode after about 10 mm of propagation distance:

numerical beam propagation with Kerr nonlinearity

Figure 4: Amplitude profile in the x-z plane for injection of the LP11 mode as calculated without the nonlinearity.

We can also show the evolution of optical powers in the guided modes:

numerical beam propagation with Kerr nonlinearity

Figure 5: Evolution of optical powers in the LP11 and LP01 mode.

It may seem surprising that the total power undergoes some oscillations: how can the power increase at some locations, even though we only have some losses (for cladding modes)? One can understand this as an energy exchange with cladding modes, enabled by the fiber's nonlinearity. Note also that the nonlinear interaction couples light into the cladding mode, which would not happen at low optical powers.

(back to the list of example cases)