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: Tapered Fiber

Here we investigate how light propagates in tapered fibers, i.e., fibers which have been heated and stretched such that the core dimension is reduced in some taper region.

Description of the Model

The reduction in fiber diameter is described with a simple function:

t_min := 0.5
t(z) := t_min + (1 - t_min) * 0.5 * (1 + cos(1 * 2pi * z / z_max))

When we later define the refractive index profile, we can easily use that function to modify a predefined two-dimensional function:

  bp_set_n_z('n_f(sqrt(x^2 + y^2) / t(z))', 'z'); { index profile }

The numerical resolution needs to be chosen such that cladding modes are realistically modeled, because we must expect that some light is lost from the core in the tapered region.


Initially, we assume a single-mode step-index fiber. The incident light is assumed to be fully in the guided mode (the LP01 mode) of the fiber. Figure 1 shows the amplitude distribution along the fiber. The gray curves illustrate the taper, reducing the fiber diameter in the middle by 50%. The red curve shows about 13% of the optical power are lost. This is because the LP01 becomes quite weakly guided in the taper region, and the taper is a bit too fast for an adiabatic transition.

amplitude distribution along the fiber

Figure 1: Amplitude distribution along the fiber.

Figure 2 shows how various parameters evolve. Three different measures are used for the beam size:

amplitude distribution along the fiber

Figure 2: Evolution of beam parameters along the fiber.

In a second simulation, we assume a fiber with doubled numerical aperture (0.2 instead of 0.1). That more strongly guiding fiber also supports LP11 modes – but not in the tapered region. If we launch the initial beam into the LP01 mode (as above), it can “survive” the tapered region without substantial losses, as it is guided quite strongly – see Figure 3. The beam size (not shown here) undergoes only quite moderate changes.

amplitude distribution along the fiber

Figure 3: Amplitude distribution along the fiber with higher numerical aperture.

If we launch into the LP11 mode, the light is entirely lost in the taper region, as seen in Figure 4:

amplitude distribution along the fiber

Figure 4: Amplitude distribution along the fiber for launching into the LP11 mode.

One sees in Figure 4 that reflections from the edges of the grid – a numerical artifact – are not completely suppressed. This, however, is no problem for the purpose of this investigation; only the detailed curve for the power drop and the calculated 0.1% transmission for the LP11 mode are not accurate – the latter should be essentially zero.

(back to the list of example cases)