pulse propagation modeling (original) (raw)
Definition: working with physical models describing the propagation of ultrashort pulses e.g. in lasers or optical fibers
Categories: 
light pulses, 
methods, 
physical foundations
Related: Tutorial on Modeling of Pulse AmplificationRaman Scattering in a Fiber AmplifierNonlinear Pulse Compression in a FiberNumerical Experiments With Soliton Pulses in FibersSoliton Pulses in a Fiber AmplifierParabolic Pulses in a Fiber AmplifierCollision of Soliton Pulses in a FiberErbium-doped Fiber Amplifier for Rectangular Nanosecond Pulsesnumerical beam propagationdispersionnonlinearitiesnonlinear pulse distortionpulse compressiondouble pulsesparabolic pulsessupercontinuum generationHaus Master equation
Page views in 12 months: 899
DOI: 10.61835/ywa Cite the article: BibTex BibLaTex plain textHTML Link to this page! LinkedIn
Content quality and neutrality are maintained according to our editorial policy.
📦 For purchasing pulse propagation modeling software, use the RP Photonics Buyer's Guide — an expert-curated directory for finding all relevant suppliers, which also offers advanced purchasing assistance.
Contents
What is Pulse Propagation Modeling?
Relevance of Pulse Propagation Effects
Simulation Models for Ultrashort Pulse Propagation
General Numerical Pulse Propagation
Summary:
This article provides a comprehensive introduction to pulse propagation modeling, a crucial tool for understanding and developing systems involving ultrashort pulses. It explains the key physical effects that alter light pulses during propagation, such as chromatic dispersion and various optical nonlinearities like the Kerr effect.
The text highlights the relevance of such modeling for applications like the design of mode-locked lasers, fiber amplifiers, supercontinuum sources, and systems for optical fiber communications.
It also details different simulation models, from analytical approaches like the Haus Master equation and soliton perturbation theory to powerful numerical techniques. The widely used split-step Fourier method is explained, and the distinction between 1D models (e.g., for single-mode fibers) and more complex 3D models (for effects like Kerr lens mode locking) is clarified.
(This summary was generated with AI based on the article content and has been reviewed by the article’s author.)
What is Pulse Propagation Modeling?
When propagating in transparent optical media, the properties of ultrashort light pulses can undergo complicated changes. Typical physical effects influencing pulses are:
- Chromatic dispersion can lead to dispersive pulse broadening, but also to pulse compression, generation of a chirp, etc.
- Various nonlinearities can become relevant at high peak powers. For example, the Kerr effect can cause self-phase modulation, and Raman scattering may e.g. induce Raman gain within the pulse spectrum (Raman self-frequency shift).
- Optical gain and losses can modify the pulse energy and the spectral shape.
- The spatial properties can change due to linear effects such as diffraction and waveguiding, but also due to nonlinear effects such as self-focusing. In highly nonlinear interactions, filamentation may occur.
Of course, different effects can act simultaneously, and often interact in surprising ways. For example, chromatic dispersion and Kerr nonlinearity can lead to soliton effects.
Raman Scattering in a Fiber Amplifier
We investigate the effects of stimulated Raman scattering in an ytterbium-doped fiber amplifier for ultrashort pulses, considering three very different input pulse duration regimes. Surprisingly, the effect of Raman scattering always gets substantial only on the last meter, although the input peak powers vary by two orders of magnitude.
For longer pulse durations (nanosecond durations or longer), the situation is normally much simpler: chromatic dispersion is quite irrelevant, and nonlinearities can substantially change the pulse spectrum, but not substantially affect the temporal profile — unless extreme processes occur at very high optical intensity, such as laser-induced damage with a dielectric breakdown.
Nonlinear Pulse Compression in a Fiber
We explore how we can spectrally broaden light pulses by self-phase modulation in a fiber and subsequently compress the pulses using a dispersive element. A substantial reduction in pulse duration by more than an order of magnitude is easily achieved, while the pulse quality is often not ideal.
Relevance of Pulse Propagation Effects
Pulse propagation effects as mentioned above are relevant in various kinds of situations. Some examples are:
Figure 1: This animated spectrogram shows how a third-order soliton evolves in a fiber. Solitons of higher orders exhibit even more complicated behavior. The image has been generated with the RP ProPulse software.
- Details of the propagation of ultrashort pulses in a mode-locked laser determine the steady-state pulse properties such as pulse duration, bandwidth and chirp, and the stability of pulse generation, multiple pulsing, etc.
- The propagation in fibers is relevant e.g. for pulse amplification, pulse compression and supercontinuum generation, and in optical fiber communications.
- Nonlinear frequency conversion of ultrashort pulses can lead to complicated changes in pulse shapes. In addition to the nonlinear interaction, there can be influences from effects such as temporal spatial walk-off and dispersive broadening.
A detailed understanding of pulse propagation is thus important, for example, for the development of ultrafast laser and amplifier systems, of supercontinuum sources and other kinds of nonlinear frequency conversion devices.
For gaining such an understanding, the modeling and simulation of pulse propagation is usually indispensable, as experimental exploration is far more cumbersome and limited. For example, there are sophisticated (and expensive) instruments for pulse characterization, but a complete characterization remains difficult, and can usually be applied only to freely accessible pulses — for example, to the output pulses of a laser system, but not to pulses at any location within a laser resonator or a fiber. Further, a computer simulation model can be used to very quickly and easily find out how the pulse propagation would react to various changes. That way, optimal performance can be realized far more efficiently.
Note that the article on laser modeling and simulation explains in detail various general aspects of modeling and simulations, which also apply to pulse propagation modeling. In the following, we treat only aspects which are specifically relevant for ultrashort pulse propagation.
Simulation Models for Ultrashort Pulse Propagation
Depending on the situation, different kinds of physical modeling techniques are required. Some of the most important ones are described in the following:
Haus Master Equation
The Haus Master equation is an analytical tool (a differential equation) mainly for calculating the steady-state pulse properties obtained in mode-locked lasers. It can be seen as a generalization of the nonlinear Schrödinger equation.
This equation is sometimes used for simulating the evolution of pulses over the resonator round trips of an ultrashort pulse in a mode-locked laser. It does not deal with details of the intracavity evolution (i.e., the evolution within a single round trip), but only describes the total effects in one round trip. It is thus mostly applicable to cases where the nonlinear and dispersive effects within a single resonator round trip are relatively weak. That is often the case in mode-locked bulk lasers, but not in mode-locked fiber lasers.
In not too complicated cases, some analytical results can directly be derived from the Haus Master equation, but it is also common to solve it numerically.
Soliton Perturbation Theory
Soliton perturbation theory describes the propagation of soliton pulses which can be subject to gain or loss, spectral filtering, or additional details of nonlinearities such as the delayed nonlinear response. A number of dynamic equations describe the evolution of the basic parameters of solitons (e.g., pulse energy and duration, center wavelength and chirp) under the influence of various effects. Also, the so-called continuum is included, i.e. a temporally broad background radiation with which a soliton can interact. Soliton perturbation theory can be used, e.g., to describe the generation of Kelly sidebands in a mode-locked fiber laser, or propagation effects in optical fiber communications.
Other Dynamic Models
Some models based on second-order moments of the complex electric field of a pulse [5] can also greatly reduce the number of dynamic variables. However, they are applicable only as long as the pulse shapes remain relatively simple.
A difficulty is that it is not always obvious where the parameter region with a reasonable accuracy ends.
The advantage of a significantly faster computation (compared with a full numerical simulation) becomes less important as the power of computers is increasing.
General Numerical Pulse Propagation
Numerical techniques are available for simulating pulse propagation in more general cases. The common approach to describe a short pulse in such models is briefly explained in the following:
- Analytically, one uses a complex amplitude ($A(t)$) in the time domain, or alternatively an amplitude ($A(\nu)$) in the frequency domain. Typically, this is a “slowly varying” amplitude, not containing the fast optical oscillation; the electric field strength is proportional to ($A(t) \exp(-i \omega_0 t)$) with the optical (angular) center frequency ($\omega_0$). The frequency-dependent amplitudes then have their maximum around ($\nu = 0$), rather than around the center frequency. The time- and frequency-dependent functions are related to each other by a Fourier transform.
Simulations on Pulse Propagation
Ultrashort pulses change in complicated ways when propagating through a fiber, for example. A suitable simulator is essential for getting complete insight — not only on the resulting output pulses, but also on the pulses at any location in your system. The RP Fiber Power software is an ideal tool for such work.
- Numerically, we represent the pulse with an array of discrete complex amplitudes in the time or frequency domain. Normally, the number of amplitudes is an integer power of 2, because fast Fourier transform algorithms work best with that.
- In some cases, a moderate number of amplitudes (e.g., 28 = 256) are well sufficient to represent a pulse. A substantially higher number of amplitudes is required in cases where, for example, pulses can be relatively long in time and broadband at the same time. A typical case is supercontinuum generation, where one sometimes uses 215 or more amplitudes.
- The chosen numerical parameters (for example, width of the considered temporal range and number of amplitudes) determine the temporal resolution, which is related to the width of the optical frequency range, and the spectral resolution. Inappropriately chosen parameters can lead to numerical artifacts.
Various aspects need to be considered for the numerical propagation of pulses:
- Linear effects such as chromatic dispersion or frequency-dependent propagation losses are easily treated in the frequency domain, whereas nonlinear interactions are often (but not always) more conveniently handled in the time domain. As required, switching between both domains can be done with a fast Fourier transform algorithm (FFT techniques).
- There are some more challenging cases where both time and frequency dependencies must be considered. An example is the saturation of frequency-dependent gain [9].
- One frequently uses the symmetrized split-step Fourier method, used particularly for pulse propagation in fibers [12]. The (weak) dispersive and nonlinear effects corresponding to short fiber pieces are alternately applied. The numerical errors associated with the finite longitudinal step size can be minimized with a special symmetrization technique, which allows for higher accuracies without excessively increased computation times.
- Note that the split-step Fourier method is similarly used in numerical beam propagation where spatial rather than temporal Fourier transforms are used. It is also possible to apply both spatial and temporal transforms.
- Automatic control of the step size (in the propagation direction) can be very important for ensuring numerical accuracy while also maintaining a high computational efficiency.
One distinguishes between simpler 1D models and more complex 3D models:
- In many cases, it is sufficient to ignore the transverse dimensions. For example, for propagation of pulses in a single-mode fiber, the transverse profile can usually be regarded as fixed (although strictly speaking the mode radius is frequency-dependent, and at extremely high intensities there are self-focusing effects). A pulse at one longitudinal position within a fiber, for example, can then be represented by a single time- or frequency-dependent complex amplitude. A model can propagate such a pulse along a fiber, for example, by using a differential equation containing influences of chromatic dispersion and nonlinearities. In the simplest case, only the simple Kerr effect is considered, but one can also include self-steepening and the delayed nonlinear response for stimulated Raman scattering. Further, in an active fiber or a laser crystal there is gain (often with relevant gain saturation).
- In special situations, such as the detailed investigation of Kerr lens mode locking or filamentation phenomena of terawatt pulses propagating in gases, one requires a full 3D model, also considering the transverse dimensions. This can be realized with sophisticated numerical beam propagation methods, including all the above-mentioned details for ultrashort pulse propagation. Such models are difficult to set up and slow to execute on a computer, but essential in some situations.
- A kind of intermediate approach, for example for pulse propagation in multimode waveguides, is based on a description of the optical field as a superposition of propagation modes, which can be coupled e.g. via nonlinearities. Similarly, one can describe interactions between light components with orthogonal polarization states.
By applying statistical techniques, pulse propagation models can also be used to investigate noise phenomena [7].
Although the construction even of 1D pulse propagation models involves various non-trivial challenges, it is relatively easy to simulate pulse propagation with sufficiently powerful and flexible software such as RP Fiber Power. It makes sense to have such models and software developed by few people and used by many others, who do not need to have programming skills and experience with numerical algorithms. A substantially less detailed understanding of the underlying physics is also sufficient.
Frequently Asked Questions
This FAQ section was generated with AI based on the article content and has been reviewed by the article’s author (RP).
What is pulse propagation modeling?
Pulse propagation modeling involves using physical models and computer simulations to analyze how the properties of ultrashort light pulses, such as their shape, duration, and spectrum, change as they travel through optical media like fibers or laser crystals.
Which physical effects influence the propagation of ultrashort pulses?
Why is the modeling of pulse propagation important?
What is the Haus Master equation?
The Haus Master equation is an analytical tool, specifically a differential equation, used to calculate the steady-state pulse properties in a mode-locked laser. It describes the total effects of dispersion and nonlinearities over one full resonator round trip.
What is the split-step Fourier method?
The split-step Fourier method is a numerical technique for simulating pulse propagation. It works by alternately applying the effects of dispersion in the frequency domain and nonlinearities in the time domain, switching between these domains using fast Fourier transforms.
What is the difference between 1D and 3D pulse propagation models?
1D models consider only the evolution of the pulse over time at a single spatial point, which is sufficient for propagation in single-mode fibers. 3D models also account for the pulse's transverse spatial dimensions, which is necessary for phenomena like Kerr lens mode locking or filamentation.
Suppliers
Bibliography
| [1] | P. V. Mamyshev and S. V. Chernikov, “Ultrashort-pulse propagation in optical fibers”, Opt. Lett. 15 (19), 1076 (1990); doi:10.1364/OL.15.001076 |
|---|---|
| [2] | G. P. Agrawal, “Optical pulse propagation in doped fiber amplifiers”, Phys. Rev. A 44 (11), 7493 (1991); doi:10.1103/PhysRevA.44.7493 |
| [3] | H. A. Haus et al., “Structures for additive pulse mode locking”, J. Opt. Soc. Am. B 8 (10), 2068 (1991); doi:10.1364/JOSAB.8.002068 |
| [4] | P. L. François, “Nonlinear propagation of ultrashort pulses in optical fibers: total field formulation in the frequency domain”, J. Opt. Soc. Am. B 8 (2), 276 (1991); doi:10.1364/JOSAB.8.000276 |
| [5] | M. Potasek et al., “Analytic and numerical study of pulse broadening in nonlinear dispersive fibers”, J. Opt. Soc. Am. B 3 (2), 205 (1992); doi:10.1364/JOSAB.3.000205 |
| [6] | D. Marcuse, “RMS width of pulses in nonlinear dispersive fibers”, IEEE J. Lightwave Technol. 10 (1), 17 (1992); doi:10.1109/50.108730 |
| [7] | R. Paschotta, “Noise of mode-locked lasers. Part I: numerical model”, Appl. Phys. B 79, 153 (2004)“,”http://link.springer.com/article/10.1007%2Fs00340-004-1547-x; R. Paschotta, “Noise of mode-locked lasers. Part II: timing jitter and other fluctuations”, Appl. Phys. B 79, 163 (2004); doi:10.1007/s00340-004-1548-9 |
| [8] | B. Burgoyne et al., “Nonlinear pulse propagation in optical fibers using second order moments”, Opt. Express 15 (16), 10075 (2007); doi:10.1364/OE.15.010075 |
| [9] | R. Paschotta, “Modeling of ultrashort pulse amplification with gain saturation”, Opt. Express 25 (16), 19112 (2017); doi:10.1364/OE.25.019112 |
| [10] | Y.-H. Chen et al., “Accurate modeling of ultrafast nonlinear pulse propagation in multimode gain fiber”, J. Opt. Soc. Am. B 40 (10), 2633 (2023); doi:10.1364/JOSAB.500586 |
| [11] | C. R. Phillips, M. Jankowski, N. Flemens and M. M. Fejer, “General framework for ultrafast nonlinear photonics: unifying single and multi-envelope treatments”, Opt. Express 32 (5), 8284 (2024); doi:10.1364/OE.513856 |
| [12] | G. P. Agrawal, Nonlinear Fiber Optics, 4th edn., Academic Press, New York (2007) |
| [13] | R. Paschotta, tutorial “Passive Fiber Optics” |
| [14] | R. Paschotta, tutorial “Passive Fiber Optics”, Part 12: Ultrashort Pulses and Signals in Fibers |
| [15] | R. Paschotta, tutorial “Modeling of Fiber Amplifiers and Lasers”, part 7 |
| [16] | R. Paschotta, tutorial “Modeling of Pulse Amplification” |
(Suggest additional literature!)
Questions and Comments from Users
Here you can submit questions and comments. As far as they get accepted by the author, they will appear above this paragraph together with the author’s answer. The author will decide on acceptance based on certain criteria. Essentially, the issue must be of sufficiently broad interest.
Please do not enter personal data here. (See also our privacy declaration.) If you wish to receive personal feedback or consultancy from the author, please contact him, e.g. via e-mail.
By submitting the information, you give your consent to the potential publication of your inputs on our website according to our rules. (If you later retract your consent, we will delete those inputs.) As your inputs are first reviewed by the author, they may be published with some delay.