# RP Fiber Power – Simulation and Design Software

for Fiber Optics, Amplifiers and Fiber Lasers

Start | Features | Model | Fiber Data | Versions | Interface | Demos |

## Mode Solver

**RP Fiber Power** contains a powerful mode solver, which can very quickly calculate all guided modes (LP modes) for a fiber from a given radially symmetric refractive index profile.
This works even for fibers with several hundred modes.
The mode solver also calculates propagation constants and the fraction of power guided in the core.

The refractive index profile can be specified with a formula or from tabulated values as result from measurements. Any wavelength dependence can also be included, which also allows one to calculate the chromatic dispersion (material and waveguide dispersion) with predefined functions.

The resulting intensity profiles can be very easily used in the definition of optical channels for the calculation of the power evolution in the fiber. (Of course, for that purpose you may also use profiles from other sources, e.g. based on the manufacturer's specifications.)

See also our demo video for using the mode solver.

## Beam Propagation

**RP Fiber Power** can numerically simulate the propagation of arbitrary complex field distributions in fibers (or other waveguides) with arbitrary (but weakly guiding) refractive index profiles.
A split-step Fourier algorithm is used.
Multiple beams at different wavelengths can be propagated, and can interact with each other through laser-active ions and nonlinear effects (XPM, SRS).

The user has full control over the numerical grid parameters, refractive index and dopant distributions, input fields, etc.
As all other features, all this can also be controlled via script commands in order to generate arbitrary plots, feed data from one beam propagation device into another (possibly with interpolation), import and export data to files with any format, etc.
This makes **RP Fiber Power** a powerful tool e.g. for research on tapered fibers, double-clad fibers and various other fiber components.

See also our demo video for numerical beam propagation.

## Level Systems and Transitions

**RP Fiber Power** can be used to model cases where two or even three different types of laser-active ions participate in an amplification or lasing process.
Each ion may have a substantial number of metastable excited levels.

The user can freely define various kinds of transitions between electronic levels:

- spontaneous transitions (radiative or non-radiative)
- optically induced transitions (absorption or stimulated emission)
- quenching processes (including pair-induced quenching and quenching processes in larger ion clusters)
- energy transfers (between ions of the same type, or of different types, e.g. between Er
^{3+}and Yb^{3+}ions)

These inputs determine the rate equation system, which is then very efficiently solved with an advanced algorithm.

The figure to the right shows an example case in the context of a thulium-doped upconversion fiber laser. Here, one is dealing with three metastable levels in addition to the ground state. There are four absorbing transitions (including the quasi-three-level laser transition) and two transitions with stimulated emission (pump and signal), apart from various spontaneous transitions.

## Optical Channels

**RP Fiber Power** allows one to take into account a large number of so-called optical channels, which may represent single propagation modes or some (possibly large) collection of propagation modes.
For modeling amplified spontaneous emission (ASE), channels can have a finite optical bandwidth.
For each channel, the transverse intensity distribution can be freely defined.
For fiber amplifiers, the noise figure of any signal channel can be calculated.

The figure below illustrates the optical channels in an erbium-doped fiber amplifier: a pump channel, two signal channels, and 16 ASE channels.

See also our demo video on fiber amplifiers and lasers.

## Dynamic Simulations

The software also allows you to simulate the temporal evolution of your system. A highly flexible approach has been implemented, allowing both the precise simulation of the dynamics of lasers, fully taking into account the propagation time in the fiber, and very fast simulations of single-pass or double-pass pulse amplification.

As an example, the figure on the right side shows how a pulse shape gets distorted by gain saturation during amplification in a high-gain fiber amplifier, and how much of the stored energy is extracted by such a pulse. Another example corresponds to a Q-switched fiber laser model, which is described on a separate page.

You may automatically perform multiple dynamic simulations. For example, you can simulate the pumping of a pulse amplifier with relatively coarse time steps, and later then signal amplification with a much higher temporal resolution. Then put both simulations into a loop to simulate multiple pumping/amplification cycles.

## Ultrashort Pulse Propagation

An ultrashort pulse can be numerically represented with an array of complex amplitudes in the time or frequency domain. Fourier transforms are used to relate time and frequency to each other. Wavelength-dependent amplification and chromatic dispersion are modeled in the frequency domain, whereas nonlinear effects are treated in the time domain. A refined numerical split-step algorithm with automatic step size control ensures a high accuracy combined with high speed.

Arbitrary chromatic dispersion can be included by specifying a wavelength-dependent function for the GVD (group velocity dispersion) or the effective refractive index. If the integrated mode solver is used, this can also calculate the mode dispersion.

Not only a simple Kerr nonlinearity can be modeled, but also a nonlinearity with a delayed nonlinear response, leading to stimulated Raman scattering. The user can specify an arbitrary response function, not only parameters of a given functional form, so that Raman gain with arbitrary gain spectrum can be treated.

Well over 100 functions of the script language are available for defining model inputs and retrieving all sorts of pulse properties. The results can also be inspected using the convenient interactive pulse display window.

Although various competing products can also simulate ultrashort pulse propagation, **RP Fiber Power** is probably unique in terms of its high flexibility, as results from the full control with a powerful script language.
Therefore, you may use it for a wide range of simulations – even going far beyond what the developer has anticipated.

## Speed and Reliability

The numerical solution of the equations for the steady state in such a model is a challenging task. It requires an iterative procedure. If this procedure is not very well worked out, its convergence properties can strongly depend on the situation. (See the Spotlight article of 2008-11-08 with a discussion of the validation of such simulations.)

**RP Fiber Power** has been very carefully designed, using an iterative procedure with automatically adapted parameters.
As a result, the software achieves a surprisingly high speed and reliable convergence in a very wide range of cases.
This includes cases with laser operation very close to threshold or far above threshold, amplification with backward pumping and strong ASE, etc., where other software may have extreme problems in finding the correct solution, or doing this within a reasonable time.

The diagram above shows with a color scale how the number of iterations required in a fiber laser model varies when the pump power and the output coupler transmission are systematically varied, and the configuration is each time calculated from scratch.
There is a slight increase of computation time near threshold, but the computation for one configuration (with 3 optical channels in that case) typically requires only a tiny fraction second on an ordinary PC, despite the good spatial resolution (with 100 points in *z* direction).

Ultrashort pulse propagation is also simulated with a very robust and efficient algorithm, automatically taking finer numerical steps only where needed, based on a monitoring of numerical errors.

See also the **frequently asked questions**.