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: Fiber Bragg Grating for Coupling into Higher-order Modes

Here we show how the numerical beam propagation feature of RP Fiber Power can be used for analyzing light propagation in fiber Bragg gratings, as long as no reflections are involved. This takes a little script programming, but you will see that it is actually quite straightforward.

Description of the Model

We first need to define the refractive index profile of the fiber from which the Bragg grating is made:

n_cl := 1.44  { cladding index }
n_co := n_cl + 2e-3  { core index }
r_co := 30 um  { core radius }
r_cl := 40 um  { region where n must be exactly constant }
L_f := 200 mm  { length of the simulated region }
n_f(r):= { index of pure fiber }
  if r <= r_cl
  then n_cl + exp(-ln(2) * (r / r_co)^8) * (n_co - n_cl)
  else n_cl

Instead of a simple step-index profile, we have used a supergaussian function.

Now we use this is an input for the mode solver in order to check how many modes we have:

calc set_n_profile("n_f", r_cl)
for m := 1 to m_max(lambda) do
 show "LP0" , m, ": ", beta_lm(0, m, lambda) * um:f6:" / µm"

We can then also use the mode solver for calculating the difference between the phase constants of the fundamental mode and of the targeted higher-order mode, from which we get the grating period required for effective coupling:

l_target := 0  { targeted mode }
m_target := 3
dbeta := beta_lm(0, 1, lambda) - beta_lm(l_target, m_target, lambda)
show "dbeta:  ", dbeta * um:f6:" / µm"
p := 2pi / dbeta
show "Period: ", p:d3:"m"

Now we define the three-dimensional refractive index profile of the fiber Bragg grating:

n_FBG(x, y, z):=  { index of fiber Bragg grating }
  var r;
  r := sqrt(x^2 + y^2);
  if r <= r_cl
  then n_cl
   + exp(-ln(2) * (r / r_co)^8) * (n_co - n_cl) * (1 + 0.1 * sin(dbeta * z))
  else n_cl;

Now we set up the numerical beam propagation:

; Grid parameters for numerical beam propagation:
r_max := 4 * r_co
N_r := 2^7
dr := 2 * r_max / N_r
dz := 200 um
N_z := L_f / dz
N_s := 20  { number of sub-steps for each dz interval }
  bp_set_grid(r_max, N_r, r_max, N_r, L_f, N_z, N_s);
  bp_set_n_z('n_FBG(x,y,z)', 'z');
  bp_set_A0('A_lm_xy(0, 1, lambda, x, y)');  { fundamental mode as the input field }

Finally, we can just use the function bp_I() for making a color plot of the resulting intensity distribution in the device. In addition, we can add some code for plotting the evolution of optical powers in the two modes:

diagram 1:
x: 0, L_f / mm
"z (mm)", @x
y: -r_max / um, +r_max / um
"y (µm)", @y
y2: 0, 100
"optical powers (%)", @y2
! I_max := 1.5 * bp_I_max(0)
cp: color_I((bp_I(0, y * um, x * mm) / I_max)),
  order := yx
P_lm(l, m, z) := 
  { power in LP_lm mode at position z }
  abs2(int(int(bp_A%(x, y, z) * A_lm_xy(l, m, lambda, x, y),
    x := -r_max to +r_max step dr),
    y := -r_max to +r_max step dr))
f: 100 * P_lm(0, m_target, x * mm), yscale = 2,
  color = red, width = 1, ["LP", l_target, m_target, " power"]
f: 100 * P_lm(0, 1, x * mm), yscale = 2,
  color = blue, width = 1, "LP01 power"

What we get is the following diagram:

coupling to a higher-order mode in a fiber Bragg grating

Figure 1: Intensity distribution inside the fiber Bragg grating. The red and blue curve shows the evolution of optical power in the LP03 mode and LP01 mode, respectively.

For interactively exploring the field profile, you can also use the interactive beam profile viewer:

beam profile viewer of the RP Fiber Power software

It is also interesting to check the evolution of mode powers in the very beginning of the fiber Bragg grating:

coupling to a higher-order mode in a fiber Bragg grating

Figure 2: Evolution of mode powers in the first 5 millimeters of the grating.

One can see that only the power in the LP03 mode rises consistently, while the other mode powers oscillate up and down on very low levels. This is because the coupled contributions to the corresponding mode amplitude from different z positions in the device are normally out of phase – except for the one mode where we have compensated that phase change by using the right grating period.

Of course, it would be easy to change various details of the model or plot various other quantities. For example, you could study the effects of asymmetric grating structures, nearly degenerate modes in strongly multimode fibers, effects of bending, or whatever else may come to your mind!

(back to the list of example cases)