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: Launching Light into a Single-mode Fiber

Description of the Model

We simulate the following case:

Some technical details of the simulation:

Note that while the mode solver is limited to cases with radial symmetry of the refractive index profile, them beam propagation could be calculated for arbitrary index profiles, as long as they are weakly guiding (which is essentially the case for all all-glass fibers).

Setting up the beam propagation with a few lines of script code is simple, after some parameters have been defined:

lambda := 1 um
; Define the refractive index profile:
n_cl := 1.45 { cladding index }
NA := 0.08
n_co := sqrt(n_cl^2 + NA^2) { core index }
r_co := 4 um { core radius }
n_f(r) := if r <= r_co then n_co else n_cl
; Grid parameters for beam propagation:
r_max := 30 um
N := 2^6
dr := 2 * r_max / N
z_max := 10 mm
dz := 10 um
N_z := z_max / dz
; Incident beam profile: Gaussian beam from a laser
w0 := 3 um { beam radius }
d := 2 um { vertical position error }
theta_in := 0 deg { angle error }
A0%(x, y) := exp(-(x^2 + (y - d)^2) / w0^2) / sqrt(0.5 * pi * w0^2)
  * expi((2pi / lambda) * x * sin(theta_in))
loss(x, y) := 1 * ((x^2 + y^2) / (10 um)^2)^6
    bp_set_grid(r_max, N, r_max, N, z_max, N_z, 2);
    bp_set_n('n_f(sqrt(x^2 + y^2))'); { index profile }
    bp_set_loss('loss(x, y)'); { loss profile }
    bp_set_A0('A0%(x, y)'); { initial amplitude }


Figure 1 shows how the amplitude distribution evolves in the fiber. We see some wiggles, which result from a vertical offset of the incident beam position. Also, one can see how some of the light escapes from the core; this is the non-guided part which is lost after some length of fiber. The red curve shows how the power in the fiber drops. The final value is in good agreement with the analytical estimate.

intensity profiles at the end of a multimode fiber

Figure 1: Field amplitude profile in the yz plane, showing how the radiation evolves in the fiber for a slightly offset input beam.

We can also systematically vary the initial beam radius and see how the launch efficiency varies. This is shown by Figure 2. Here, the beam is assumed to be perfectly aligned to the core.

launch efficiency as a function of the initial beam radius

Figure 2: Launch efficiency as a function of the initial beam radius.

The numerical beam propagation results nicely agree with the analytical calculation. Note that they could also easily be carried out in more complicated cases. For example, one might have a fiber core which has no radial symmetry, any other kind of input beam profile, or a bent fiber.

(back to the list of example cases)