〉〉 先端技術

Fiber Laser and Amplifier Simulations in FETI
Zoltán Várallyay*1, Gábor Gajdátsy*1, András Cserteg*1, Gábor Varga*2 and Gyula Besztercey*3〈概要〉
ファイバレーザは，出力レベルや堅牢性が固体レーザの特性に匹敵するものになってきて，レーザ市
場における存在感が増してきているとともに，ますます要求が高まってきている。その理由は，ひとつにはフォトニック結晶ファイバに代表されるような最新の光ファイバ関連の技術開発成果を使えること，またこれまでに蓄積されてきた多種の光ファイバ技術を用いて洗練されたシステム設計ができることなどが考えられる。ある特定の応用先が決まった時に，必要となるパラメータを最適化したり特性値を最大化したりするために，システムの特性をシミュレーションすることは必要不可欠である。FETI では，連続光（CW）発振器としてだけではなく，パルス発振器や増幅器としてのファイバレーザのシステ
ム特性を計算するアルゴリズムを有している。このアルゴリズムには，ある思想で定義された評価関数を最大化するように系を最適化する際の，「しらみつぶし法」や「単体法」も取り入れることが可能である。
分散や非線形，その他損失や利得を考慮することに伴う様々な効果も取り入れることができる。具体的には，非線形のシュレディンガー方程式とレート方程式を出発点として，様々なソルバーを組み合わせることで，希土類をドープした光ファイバ中のパルス伝搬を計算できるアルゴリズムを構築している。
このアルゴリズムを用いると，直線ファイバでも，リングファイバ発振器でも，パワー増幅器でも，それぞれ境界条件を満足させるような計算を行うことが可能である。本報告書では，開発したアルゴリズ
ムをいくつかのサブシステムに適用させた例を示している。アルゴリズムの確認のみではなく，計算結
果はより安定な超高出力ファイバレーザシステムを開発する道筋を示すものとなっている。

INTRODUCTION
1.1. INTRODUCTION
At the end of the 90s the output power from a diffraction
limited or nearly diffraction limited fiber laser (FL) was
restricted to a few multiple of 10 W power level while
recently this value is increased by three orders of magnitude as in the case of the Yb-doped fiber technology1).
This power level could be achieved, on one hand, by the
introduction of double-clad rare-earth doped fibers 2)
which made possible the use of high power multi-modelaser diodes (LD) as pump sources injected into the large
area cladding, because single mode fiber lasers pumped
with single mode pump diodes are featuring low output
power. On the other hand, the work on such fiber designs
which present an increased core area in order to reduce
the nonlinear distortions of the propagating signal while
keeping still a diffraction limited output field distribution
contributed to the extension of the achievable power levels too3),4). A big leap in the high quality, high energy laser
outputs was achieved by the introduction of the large
mode area (LMA) or large pitch (LP) fibers having photonic crystal cladding where the cladding filters out the high*1 FETI, Simulation group*2 Budapest University of Technology and Economics, Physics
department*3 FETI, General Manager

er order modes providing a nearly diffraction limited output5). Aside from the mode instability6), these type offibers are theoretically infinitely single mode fibers7). Thismode scaling of the fiber core8) and the application of the
chirped pulse amplification (CPA) technique9) can result inenormous output power levels from fiber lasers10). The
output level is comparable to the power levels of solidstate laser systems.
At FETI, we have been modeling nonlinear wave propagation in optical fibers for more than 10 years. We used
our model successfully to predict the nonlinear compression of broad and ultra-short laser pulses in a small corearea fiber11). We also developed and verified amplifiermodels which are able to treat the rate equations inYtterbium (Yb) or Erbium (Er) doped fiber amplifiers andconsequently calculate the gain for the signal, pump andamplified spontaneous emission (ASE) fields in a wide
wavelength range12). We note also that these simulations
can be extended to include the differential gain of the different core modes in case of few modes or multimodeamplifiers13). If one has an amplifier model that provide
reliable calculation results in a certain parameter extent
that model can be easily transfer to modeling fiber oscillators since the gain medium in both devices are identicaland only the boundary conditions are distinct from thephysical point of view. Of course, a fiber oscillator which
古河電工時報第 135 号（平成 28 年 2 月） 57

一般論文

Fiber Laser and Amplifier Simulations in FETI〉〉 先端技術

can provide short laser pulses via Q-switching or modelocking14) is a bit more complex since additional opticalelements such as saturable absorbers and also components with spectral filtering have to be modeled in thesame time. We managed to concatenate these elements
in a fiber ring oscillator and successfully modeled themode-locked, output pulse properties of the oscillator atdifferent parameter selections for the saturable absorber15). The simulation of continuous wave (CW) oscillators
provide a different challenge at high power levels. We had
to include the nonlinear effects in the calculations: they
are crucial because a real CW laser has also finite bandwidth but the Fourier-transform of a finite bandwidth withconstant phase will not be CW in the time domain. By
introducing a phase-diffusion model, the connection
between the spectral domain and the time domain can beestablished and the split-step Fourier method can beapplied to add nonlinear and linear effects to the signal
evaluation16).
In this paper, we describe our physical and numericalmodels how we combine the rate equations associatedwith the gain effects with the nonlinear Schrödinger equation taking into account the dispersive and nonlinear
effects during the light propagation. Using this extended
numerical model, we present simulation results which
intent is to optimize a high power CPA system using LMA
or LP fibers as gain media. The aim of the optimization
will be to determine the optimum applied chirp on the
short pulses at the input end of the fiber to avoid any
nonlinear distortion during the amplification process till
the output end of the fiber.

2.2.THEORY
THEORY
To model the gain properties of the doped optical fiber,
we solve the coupled power evaluation equations of thedifferent signals along with the steady-state, two-levelrate equation17)
(1)

(2)

where + and - in the superscript of P denote the forward
and backward propagating signals, respectively and consequently u is 1 or -1 for forward and backward propagating signals. Γ is the overlap factor of the propagating
mode with the doped region, Cd is the doping concentration, α is the so-called background loss of the fiber without the doping ions, σe and σa are the emission andabsorption cross-sections of the doping ions in the particular host material. These cross-section values have to bemeasured in order to be used in this model12). N2 and N1
are the populations of the metastable and fundamental
states, respectively. In Equation. (2), ξ=πReffh /τ whereτ is the fluorescence lifetime of the doping ion, h is the
Planck constant, Reff is the effective radius of the doped
region and f k is the kth frequency component. One has to
add an amplified spontaneous emission (ASE) term toEquation (1) if the calculated signal is the forward or the
backward ASE signal. This term has a form ofuσeΓCdN2nhf k∆f where n is the number of the propagating modes and ∆f is the ASE frequency resolution.
The dispersion and nonlinearity related effects are governed by the nonlinear Schrödinger equation18):(3)
where A is the complex envelop function of the investigated pulse. The first term at the right-hand side of theequation takes into account the dispersion effects fromthe second order to the higher orders where ß m is the mthorder dispersion contribution of the fiber to the pulse
evaluation and T is the time space in a frame of reference
travelling with the pulse. The second term at the righthand side is the term taking into account the nonlinear
contribution to the pulse evaluation. Here, R(t) is the nonlinear response function and γ is the nonlinear coefficientand they are discussed in details in References8). The
third term could be the loss term but that is alreadyincluded in Equation (1) therefore this term is kept for thegain in the fiber (negative loss). The frequency dependent
G can be calculated from Equation (1) and Equation (2)using a small segment (∆z ) of the fiber:(4)
When solving Equation (1)-(3) with the help of Equation(4), one has to consider that the width of the temporalwindow, where we generate the complex envelope function, is not identical to the time range determined by therepetition rate of the pulse. Therefore, the spectral intensity during the calculations must be corrected accordingly.
The boundary conditions for Equation (1) in case of
amplifier modeling are given in References12) and for a linear fiber oscillator, it is described in References 16). In
order to treat the bidirectional propagation in the amplifier
with great stability we apply a modified shooting method
which is detailed in References12).

RESULTS
3.3.RESULTS
We model an Yb-doped, LMA amplifier to amplify 1μJ
Gaussian pulses at around 1070 nm originating from a 1
MHz fiber oscillator and a pre-amplifier system which
amplified the pulse to the 1 W power level (1μJ·1 MHz).
This pulse is amplified further in the LMA fiber having aneffective core area (ECA) of 4000μm2. Although this core
size is large enough to present low nonlinearity, the mentioned power level will be too large to propagate without
significant spectral broadening if the transform limited
古河電工時報第 135 号（平成 28 年 2 月） 58

一般論文

Fiber Laser and Amplifier Simulations in FETI〉〉 先端技術

pulse duration with its high peak power is the input signal. Fig. 1 shows that even after 10 cm of propagation the
spectral bandwidth is broadened more than 10 times
compared to the original full width at half maximum(FWHM). The used dispersion parameters are the dispersion values of the silica glass at around 1070 nm (nowaveguide contribution at large core sizes). Linear dispersion: D = -28.4366 ps/(nm km), dispersion slope: S =0.3146 ps/(nm 2 km), third order dispersion: T =-7.2853·10-4 ps/(nm3 km) and fourth order dispersion: F =4.6057·10-6 ps/(nm4 km).
Though the doping ions will modify the dispersion of asilica fiber, its small contribution to the above dispersion
values will not alter the calculation results significantly.
Therefore, we calculate with these dispersion values
throughout this paper. Conversely, the nonlinear refractive
index of the fiber is the nonlinear refractive index of thesilica glass which is n2 =2.6·10-20 m2/W.
Figure 1 (b) shows not only the changes of the spectral
FWHM but also the broadening of the temporal width
which extends to 465 fs from the initial 170 fs during 10
cm of propagation due to the dispersion contribution of
the fiber. Due to the higher order dispersions, the modulated spectra becomes slightly asymmetric (Figure 1 (a)).

(a) SPECTRA
Spectral intensity (arb. Units)

8.00E-04
7.00E-04
Input

6.00E-04

Output5.00E-04
4.00E-04
3.00E-04
2.00E-04
1.00E-04

0.00E+00
9.50E+02 1.00E+03 1.05E+03 1.10E+03 1.15E+03 1.20E+03 1.25E+03

Wavelength (nm)

0.6

120

0.5

100

0.4

80

0.3

60

0.2

40
Temporal0.1
0-0.005Spectral

20

Spectral FWHM (nm)Temporal FWHM (ps)(b) TEMPORAL AND SPECTRA FWHM ALONG THE FIBER

0
0.015

0.035

0.055

0.075

0.095

Distance (m)Figure 1 (a) Input and output spectra before and after a 10 cm
LMA fiber amplifier and (b) Temporal and spectral
F WHM evaluations along the fiber without any
amplification.

If we switch on the gain lunching cladding pumps inforward and backward directions in the LMA fiber and weuse at least one meter length of it in order to achieve a
noticeable gain, the temporal and spectral distortions will
be more significant. Therefore, we will use a brute-force
optimization on the system by adding a linear chirp to the
input pulse to broaden the pulse-width and decrease thepeak power sufficiently. At the output, we will compress
the amplified pulse applying linear and second-orderchirps to compensate the phase on the pulse. The meritfunction (MF) of the CPA system is to achieve maximum
compression at the output but simultaneously the highest
pulse quality after the compression. We define the qualityof the pulse as a ratio of stored energy in the main peak
and the total pulse energy (References 11)). Since wewish to obtain the best possible quality pulses at the
shortest possible duration, we found that a proper merit
function considering the pulse quality to a greater weightshould have the form of(5)
where ∆τ is the FWHM of the compressed pulse and
QF is the obtained quality factor of the compressed pulse(a number between 0 and 1) and x is an exponential factor. We found that the exponential factor should be larger
than 3 to obtain a reliable optimum (high quality compressed outputs).
We will minimize MF in Equation (5) and to do this, we
use the following setup of the arrangement: the gain fiberis 1 m long having a core diameter of 88μm and a cladding diameter for the pump of 200μm8). This fiber has an
ECA of 4000μm2. The ratio of the core and cladding area
determines the overlap factor of the cladding pump withthe doping ions in the core12),13) and this way we have
Γ=0.1936. The doping concentration is set to 1026 1/m3
and the fluorescence lifetime of Ytterbium is set to 2.3
ms. Both pumps (forward, backward) have 4 nm spectral
bandwidth and 80 W CW output power injected into thefiber cladding. The input signal has a 10 nm initial bandwidth and 1 W power level which corresponds to 170 fs
transform limited pulse duration at around 1070 nm and5.6 MW peak power. To avoid nonlinear distortions, weare looking for the optimum input chirp on the input signal that way to obtain a minimum for Equation (5).
The shortest pulse width and the possible highest quality factor is found to be at 6.65 ps2 linear, input chirp
scanning the input chirp values in hundred steps between
3 and 8 ps2 repeating the calculations with each pulses.
The corresponding FWHM and QF using a range of compression chirps can be seen in Figure 2. where that is
obvious that the shortest pulse duration does not always
meet with the highest quality factor of the pulse whichmakes necessary to introduce Equation (5) the way we
did. Fig. 3 shows the steady-state inversion level in thegain fiber along with the gain of the signal at the center
frequency. Figure 4 (a) and (b) show the power evolution
古河電工時報第 135 号（平成 28 年 2 月） 59

Fiber Laser and Amplifier Simulations in FETI(a) FWHM VS COMPRESSION CHIRPSINVERSION AND GAIN ALONG THE FIBER
Inversion
0.5

25

Gain

0.4

20

0.3

15

0.2

10

0.1

5

0
-0.05

Gain at 1070 nm (dB/m)

30

0.6

Inversion levelof the propagating signals as well as the temporal andspectral FWHM along the fiber during this amplification
process. Finally, the temporal and spectral shape of thepulse are shown in Figure 5 (a) and (b).
One can see that the first step to amplify high power
pulses is stretching them and the added 6.65 ps2 input
chirp will broaden the pulse to an FWHM of 109.4 ps (SeeFigure 4 (b)) corresponding to a reduced peak power of
8.6 kW. During the amplification and propagation, thepulse width becomes 111.3 ps and due to the gain thepeak power will reach a little more than 1 MW. During thisprocess, the spectral width is increased only by 0.6 nm(Figure 4(b)) which small broadening adumbrates a goodquality compression. The best compression is achieved
at -6.46 ps2 compression chirp and 0.004 ps3 third orderchirp (Figure 2 and Figure 5 (a)). The obtained FWHM ofthe compressed pulse we got is 522 fs with a pulse peak
power of 148.1 MW (Figure 5 (a)).

〉〉 先端技術

0
0.15

0.35

0.55

0.75

0.95

Distance (m)Figure 3 Inversion level and gain along the 1m long amplifier.

(a) POWER EVOLUTION ALONG THE FIBER
100

Total power (W) /logscale/

一般論文

10
1
0.1
0.01

SignalForward pump
Backward pump
Forward ASE
Backward ASE

0.001
0.00010.00001
0.000001
-0.05

0.15

0.35

0.55

0.75

0.95

Distance (m)(b) TEMPORAL AND SPECTRAL FWHM ALONG THE FIBER

Temporal FWHM (ps)(b) QUALITY FACTOR VS COMPRESSION CHIRPS

10.6
111

10.5
10.4

110.5

10.3
110

10.2
Temporal

109.5

Spectral
109
-0.05

10.1

Spectral FWHM (nm)

10.7

111.5

10
9.9

0.15

0.35

0.55

0.75

0.95

Distance (m)Figure 4 (a) Power evolution of the ASE signals, cladding
pumps and the amplified pulse.
(b) Temporal and spectral FWHM of the amplified
pulse along the fiber due to the dispersion andnonlinear effects.

Figure 2 (a) FWHM. (b) QF as functions of compression groupdelay dispersion (GDD) and compression third-orderdispersion (TOD). Color is showing the magnitudeof FWHM and QF in pico-seconds and in proportion,
respectively.

古河電工時報第 135 号（平成 28 年 2 月） 60

Fiber Laser and Amplifier Simulations in FETI〉〉 先端技術(b) SPECTRA

Pulse power (W) /logscale/(a) TEMPORAL SHAPES
1.00E+06
1.00E+05
1.00E+04
1.00E+03
1.00E+02
1.00E+01
1.00E+00
1.00E-01
1.00E-02
1.00E-03
1.00E-04
1.00E-05
1.00E-06
1.00E-07
1.00E-08-270

1.00E-02

Input
OutputCompressed

-170

-70

30

130

230

Duration (ps)Spectral power (arb. units) /logscale/

一般論文

We showed in this paper that our developed software that
calculates the pulse amplification and the nonlinear, dispersive propagation can be used in connection with highpower LMA fibers in CPA systems. Our calculations are
based on a signal having 10 nm bandwidth as an input
for a 1 m long LMA amplifier with 4000μm2 ECA. The signal had 1 W input power and 1 MHz repetition rate. Using
sufficient amount of chirp on the input pulse, we managed to amplify it without significant spectral and temporal distortion to 126.4 W average power level that corresponds to 148.1 MW peak power and 522 fs pulse duration after ideal chirp compensation at the output of theamplifier.

REFERENCES
1) D. J. Richardson, J. Nilsson, and W. A. Clarkson, “High power fiber
lasers: current status and future perspectives [Invited]”, JOSA B,
27, B63-B92 (2010).
2) H. Po, E. Snitzer, L. Tumminelli, F. Hakimi, N. M. Chu, T. Haw“Double clad high brightness Nd fiber laser pumped by GaAlAs
phased array”, Proceedings of the Optical Fiber Communication
Conference, 1989, PD7.
3) J. C. Jasapara, M. J. Andrejco, A. DeSantolo, A. D. Yablon, Z.
Várallyay, J. Fini, D. J. DiGiovanni, and C. Headley, “Diffraction limited fundamental mode operation of core pumped very-large-modearea Er fiber amplifiers”, IEEE Journal of Selected Topics inQuantum Electronics 15, 3-11 (2009).
4) J. C. Jasapara, A. DeSantolo , J. W. Nicholson, A. D. Yablon, and Z.
Várallyay, “Diffraction limited amplification of picosecond pulses in
1170 μm2 effective area erbium fiber”, Optics Express 16, 1886918874 (2008).
5) J. Limpert, N. Deguil-Robin, I. Manek-Hönninger, F. Salin, F. Röser,A. Liem, T. Schreiber, S. Nolte, H. Zellmer, A. Tünnermann, J.
Broeng, A. Petersson, and C. Jakobsen, “High-power, rod-type photonic crystal fiber laser”, Opt. Express, 13, 1055-1058 (2005).
6) C. Jauregui, T. Eidam, J. Limpert, and A. Tünnermann, “The impact
of modal interference on the beam quality of high-power fiber
amplifiers”, Opt. Express 19, 3258–3271 (2011).
7) P. St. J. Russell, “Photonic crystal fibers”, Science 299, 358–362(2003).
8) J. Limpert, F. Stutzki, F. Jansen, H.-J. Otto, T. Eidam, C. Jauregui

1.00E-08

Input

1.00E-11

Output

1.00E-14
1.00E-17
1.00E-20
1.00E-23
1.00E-26
1.00E-29
1.00E-32
1.00E-35
9.50E+02

Figure 5 (a) Input, output and optimally compressed pulse shapes.
(b) The slightly broadened and significantly amplified spectrum.

CONCLUSION
4.4.CONCLUSION

1.00E-05

1.00E+03

1.05E+03

1.10E+03

1.15E+03

1.20E+03

Wavelength (nm)and A. Tünnermann, “Yb-doped large-pitch fibres: effective singlemode operation based on higher-order mode delocalisation”, Light:
Science & Applications 1, 1–5 (2012).
9) D. Strickland and G. Mourou, “Compression of amplified chirped
optical pulses”, Opt. Commun. 56, 219–221 (1985).
10) J Limpert , T. Eidam, M. Baumgartl, F. Röser, M. Plötner, B. Ortaç,S. Nolte, A. Tünnermann, “Compact Ultrafast Oscillators and High
Performance Ultrafast Amplifiers Based on Ytterbium-Doped
Fibers”, Springer Series in Optical Sciences, 195, 75-91 (2015).
11) Z. Várallyay, J. Fekete, Á. Bányász and R. Szipőcs, “Optimizing
input and output chirps up to the third-order for sub-nanojoule,ultra-short pulse compression in small core area PCF”, Applied
Physics B 86, 567-572 (2007).
12) Z. Várallyay, Á. Szabó, A. Rosales, E. Gonzales, H. Tobioka, and C.
Headley, “Accurate modeling of cladding pumped, star-shaped,
Yb-doped fiber amplifiers”, Optical Fiber Technology 21, 180–186(2015).
13) Z. Várallyay, and J. C. Jasapara, “Comparison of amplification inlarge area fibers using cladding-pump and fundamental-mode corepump schemes”, Optics Express 17, 17242-17252 (2009).
14) A. Cserteg, V. Sági, A. Drozdy, Z. Várallyay and G. Gajdátsy,"Characterization of mode locking in an all fiber, all normal dispersion ytterbium based fiber oscillator", SPIE LASE 93442C (2015).
15) Á. Szabó, and Z. Várallyay, “Numerical study on the saturable
absorber parameter selection in an Er fiber ring oscillator”,
Photonics Technology Letters, 24, 122-124 (2012).
16) Á. Szabó, and Z. Várallyay, A. Rosales, and C. Headley, “A fast algorithm to predict spectral broadening in CW bidirectionally pumped
high power Yb-doped fiber lasers”, Applied Physics B 0946-2171,
1-7 (2015).
17) C. R. Giles, and E. Desurvire, “Modeling Erbium-Doped Fiber
Amplifiers. Journal of Lightwave Technology”, 9, 271–283 (1991).
18) G. P. Agrawal, “Nonlinear Fiber Optics”, Third Edition, Academic
Press 2007.

古河電工時報第 135 号（平成 28 年 2 月） 61