Search tips
Search criteria 


Logo of wtpaEurope PMCEurope PMC Funders GroupSubmit a Manuscript
Phys Rev A. Author manuscript; available in PMC 2010 December 10.
Published in final edited form as:
PMCID: PMC3000527

Soliton absorption spectroscopy


We analyze optical soliton propagation in the presence of weak absorption lines with much narrower linewidths as compared to the soliton spectrum width using the novel perturbation analysis technique based on an integral representation in the spectral domain. The stable soliton acquires spectral modulation that follows the associated index of refraction of the absorber. The model can be applied to ordinary soliton propagation and to an absorber inside a passively modelocked laser. In the latter case, a comparison with water vapor absorption in a femtosecond Cr:ZnSe laser yields a very good agreement with experiment. Compared to the conventional absorption measurement in a cell of the same length, the signal is increased by an order of magnitude. The obtained analytical expressions allow further improving of the sensitivity and spectroscopic accuracy making the soliton absorption spectroscopy a promising novel measurement technique.


Light sources based on femtosecond pulse oscillators have now become widely used tools for ultrashort studies, optical metrology, and spectroscopy. Such sources combine broad smooth spectra with diffraction-limited brightness, which is especially important for high-sensitivity spectroscopic applications. Advances in near- and mid-infrared femtosecond oscillators made possible operation in the wavelength ranges of strong molecular absorption, allowing direct measurement of important molecular gases with high resolution and good signal-to-noise ratio [1]. At the same time, it was observed that such oscillators behave quite differently, when the absorb-ing gas fills the laser cavity or introduced after the output mirror [2, 3]. The issue has become especially important with introduction of the mid-IR femtosecond oscillators such as Cr:ZnSe [4], which operate in the 2–3 μm wave-length region with strong atmospheric absorption.

As an example, Figure 1 presents a typical spectrum of a Cr:ZnSe femtosecond oscillator, operating at normal atmospheric conditions. It is clearly seen, that the pulse spectrum acquires strong modulation features which resemble the dispersion signatures of the atmospheric lines. Being undesirable for some applications, such spectral modulation might at the same time open up interesting opportunity of intracavity absorption spectroscopy. Compared with the traditional intracavity laser absorption spectroscopy [5, 6] based on transient processes, this approach would have an advantage of being a well-quantified steady-state technique, that can be immediately coupled to frequency combs and optical frequency standards for extreme accuracy and resolution.

FIG. 1
Output spectrum of a 100-fs Cr:ZnSe oscillator (black solid line) when operated at open air. The atmospheric transmission (gray) is calculated from HITRAN database [7] and corresponds to a full round-trip. The lower graph (b) shows the expanded central ...

In this paper, we present a numerical and analytical treatment of the effect of a narrowband absorption on a femtosecond pulse, considered as a dissipative soliton. Such a treatment covers both, passively modelocked ultrashort pulse oscillators with intracavity absorbers, and soliton propagation in fibers with impurities. The theoretical results are compared with the experiment for a femtosecond Cr:ZnSe oscillator operating at normal atmospheric conditions. We prove that the spectral modulation imposed by a narrowband absorption indeed accurately follows the associated index of refraction when the absorber linewidth is sufficiently narrow.


Our approach is based on the treatment of an ultrashort pulse as one-dimensional dissipative soliton of the nonlinear complex Ginzburg-Landau equation (CGLE) [8, 9]. This equation has such a wide horizon of application that the concept of “the world of the Ginzburg-Landau equation” has become broadly established [10]. In particular, such a model describes a pulse with the duration T0 inside an oscillator or propagating along a nonlinear fiber.

To obey the CGLE, the electromagnetic field with the amplitude A(z, t) should satisfy the slowly-varying amplitude approximation, provided by the relation ω0 >> 1/T0, where ω0 is the field carrier frequency, t is the local time, and z is the propagation coordinate. This approximation is well satisfied even for pulses of nearly single optical cycle duration [11]. When additionally we can neglect the field variation along the cavity round-trip or the variation of material parameters along a fiber, as well as the contribution of higher-order dispersions, the amplitude dynamics can be described on the basis of the generalized CGLE [9, 1214]

[partial differential]A(z,t)[partial differential]z={σ+Σ^[P(z,t)]iγP(z,t)}A(z,t)+(α+iβ22)[partial differential]2[partial differential]t2A(z,t)+Γ^[A(z,t)],

where P [equivalent] |A|2 is the instant field power and α is the inverse gain bandwidth squared. The nonlinear terms in Eq. (1) describe i) saturable self-amplitude modulation (SAM) with nonlinear gain defined by the non-linear operator Σ^, and ii) self-phase modulation (SPM), defined by the parameter γ. For a laser oscillator, γ = 4πnn2lcryst0Aeff. Here λ0 is the wavelength; n and n2 are the linear and nonlinear refractive indexes of an active medium, respectively; lcryst is the length of the active medium; Aeff = πω2 is the effective area of a Gaussian mode with the radius w inside the active medium. The propagation coordinate z is naturally normalized to the cavity length, i.e. z becomes the cavity round-trip number. For a fiber propagation, γ = 2πnn20Aeff, where n and n2 are the linear and nonlinear refractive indexes of a fiber, respectively, and Aeff is the effective mode area of the fiber [15]. Finally, β2 is the round-trip net group delay dispersion (GDD) for an oscillator or the group velocity dispersion parameter for a fiber with β2 < 0 corresponding to anomalous dispersion.

The typical explicit expressions for Σ^[P] in the case, when the SAM response is instantaneous, are i) Σ^[P]=κP (cubic nonlinear gain), ii) Σ^[P]=κ(PζP2) (cubic-quintic nonlinear gain), and iii) Σ^[P]=κP/(1+ζP) (perfectly saturable nonlinear gain) [9, 16]. The second case corresponds to an oscillator mode-locked by the Kerr-lensing [8]. The third case represents, for instance, a response of a semiconductor saturable absorber, when T0 exceeds its excitation relaxation time [17]. However, if the latter condition is not satisfied, one has to add an ordinary differential equation for the SAM and Eq. (1) becomes an integro-differential equation (see below).

The σ-term is the saturated net-loss at the carrier frequency ω0, which is the reference frequency in the model. This term is energy-dependent: the pulse energy E(z)[equivalent]P(z,t)dt can be expanded in the vicinity of threshold value σ = 0 as σ ≈ δ (E/E* – 1) [18], where δ = [ell]2/g0 ([ell] is the frequency-independent loss and g0 is the small-signal gain, both for the round-trip) and E* is the round-trip continuous-wave energy equal to the average power multiplied by the cavity period.

The operator Γ^ describes an effect of the frequency-dependent losses, which can be attributed to an absorp-tion within the dissipative soliton spectrum. That can be caused, for instance, by the gases filling an oscillator cavity or the fiber impurities for a fiber oscillator. Within the framework of this study, we neglect the effects of loss saturation and let Γ^ be linear with respect to A (z, t). The expression for Γ^[A(z,t)] is more convenient to describe in the Fourier domain, Ã(z, ω) being the Fourier image of A(z, t). If the losses result from the l independent homogeneously broadened lines centered at ωl (relative to ω0) with linewidths Ωl and absorption coefficients An external file that holds a picture, illustration, etc.
Object name is ukmss-32699-ig0008.jpgl < 0, then the action of operator Γ^ can be written down in the form of a superposition of causal Lorentz profiles [19, 20]

Γ^[A~]=(l[sm epsilon]l1+i(ωωl)/Ωl)A~(z,ω).

In the more general case the causal Voigt profile has to be used for Γ^[A~] [21]. Causality of the complex profile of Eq. (2) demonstrates itself in the time domain, where one has

Γ^[A(z,t)][proportional, variant]l[sm epsilon]lΩlte(Ωliωl)(tt)A(z,t)dt.

The conventional analysis of perturbed soliton propagation includes approximation of the effective group delay dispersion of the perturbation as Taylor series β(ω)=β2(ωω0)2/2+β3(ωω0)3/6, assuming that the additional terms β2, β3, … are sufficiently small. This approach is absolutely not applicable in our case, because the dispersion, associated with a narrow linewidth absorber can be extremely large. For example, an atmospheric line with a typical width Ω = 3GHz and peak absorption of only 10−3 produces the group delay dispersion modulation of β2=±0.9ns2, far exceeding the typical intracavity values of β2 ~ 102 ..104 fs2. Moreover, decreasing the linewidth Ω (and thus reducing the overall absorption of the line) causes the group delay dispersion term to diverge as β2[proportional, variant]Ω2.

In the following we shall therefore start with a numerical analysis to establish the applicability and stability of the model, and then present a novel analysis technique based on integral representation in the spectral domain.


Without introducing any additional assumptions, we have solved the Eqs. (1,2) numerically by the symmetrized split-step Fourier method. To provide high spectral resolution, the simulation local time window contains 222 points with the mesh interval 2.5 fs. The simulation parameters for the cubic-quintic version of Eq. (1) are presented in Table I. The GDD parameter β2 = −1600 fs2 provides a stable single pulse with FWHM ≈100 fs. The single low-power seeding pulse converges to a steady-state solution during z ≈5000.

Laser simulation parameters. The numbers correspond to a Cr:ZnSe femtosecond oscillator of Fig. 1 with lcryst =0.4 cm, w =80 μm, λ0 =2.5 μm, n = 2.44, n2 = 10−14 cm2/W, [ell] =0.075, g0 =2.5[ell].

The pulse propagation within a linear medium (e.g. an absorbing gas outside an oscillator, a passive fiber containing some impurities, or microstructured fiber filled with a gas) is described by Eqs. (1,2) with zero α, γ, Σ^, and the initial A (0, t) corresponding to output oscillator pulse. The obvious effect of the absorption lines on a pulse spectrum are the dips at ωl (Fig. 4a), that simply follow the Beer’s law. This regime allows using the ultrashort pulse for conventional absorption spectroscopy [1]. The nonzero real part of an absorber permittivity (i.e. I(Γ^)0) does significantly change the pulse in time domain [22], but does not alter the spectrum. The pulse spectrum reveals only imaginary part of an absorber permittivity (i.e. R(Γ^)0).

FIG. 4
(Color online) Central parts of the dissipative soliton spectra in an oscillator with the single absorption line centered at ω =0. (a) An external file that holds a picture, illustration, etc.
Object name is ukmss-32699-ig0008.jpg = −0.05; Ω = 4 GHz (solid curve, open circles and crosses) and 1 GHz (dashed curve, open squares ...

Introducing the nonzero SPM coefficient γ=[mid ]β2[mid ]/A02T02 with zero α and Σ^ transforms Eq. (1) to a perturbed nonlinear Schrödinger equation and results in the true perturbed soliton propagation. In this case, as shown in Fig. 2b, the situation becomes dramatically different. Besides the dips in the spectrum shown by the gray curve corresponding to a contribution of R(Γ^) only, there is a pronounced contribution from the phase change induced by the dispersion of absorption lines (solid curve in Fig. 2b corresponds to the complex profile of Γ^ in Eq. (2)). As a result, the spectral profile has the sharp bends with the maximum on the low-frequency side and the minimum on the high-frequency side of the corresponding absorption line. At the same time, the dips in the spectrum due to absorption are strongly suppressed. In addition to spectral features, the soliton decays, acquires a slight shift towards the higher frequencies, and its spectrum gets narrower due to the energy loss.

FIG. 2
(Color online) Part of the pulse spectrum after: (a) linear propagation for 25 dispersion lengths inside a fiber with two absorption lines and (b) perturbed soliton propagation for 100 dispersion lengths. Gray curve in (b) corresponds to the contribution ...

The soliton spectrum reveals in this case the real part of the absorber permittivity. However, the continuous change of the soliton shape due to energy decay renders the problem as a non-steady state case. The situation becomes different in a laser oscillator, where pumping provides a constant energy flow to compensate the absorption loss.

Let us consider the steady-state intra-cavity narrowband absorption inside a passively modelocked femtosecond oscillator, where the pulse is controlled by the SPM and the SAM, which is described by the cubic-quintic Σ^ in Eq. (1) modeling the Kerr-lens mode-locking mechanism [12]. Such an oscillator can operate both, in the negative dispersion regime [12] with an chirped-free soliton-like pulse, and in the positive dispersion regime [18], where the propagating pulse acquires strong positive chirp. In this study, we consider only the negative dispersion regime, the positive dispersion regime will be a subject of following studies.

The results of the simulation are shown in Figs. Figs.33 and and4,4, and they demonstrate the same dispersion-like modulation of the pulse spectrum. Fig. 3 demonstrates action by three narrow (Ω =2 GHz) absorption lines centered at −10, 0 and 10 GHz in the neighborhood of ω = 0. One can see (Fig. 3a), that the absorption lines do not cause spectral dips at ωl, but produce sharp bends, very much like the case of the true perturbed Schrödinger soliton considered before. One can also clearly see the collective redistribution of spectral power from higher- to lower-frequencies, which enhances local spectral asymmetry (Fig. 3a). Such an asymmetry suggests that the dominating contribution to a soliton perturbation results from the real part of an absorber permittivity, which, in particular, causes the time-asymmetry of perturbation in the time domain. This asymmetry is seen in time domain as a ns-long modulated exponential precursor in Fig. 3b.

FIG. 3
(Color online) Dissipative soliton in an oscillator: (a) central part of the spectrum and (b) power P(t). An oscillator is filled with an absorber described by Eq. (2) with the triplet of lines: An external file that holds a picture, illustration, etc.
Object name is ukmss-32699-ig0008.jpg1 = An external file that holds a picture, illustration, etc.
Object name is ukmss-32699-ig0008.jpg2 = An external file that holds a picture, illustration, etc.
Object name is ukmss-32699-ig0008.jpg3 = −0.005, Ω1 = Ω2 = Ω ...

The simulated effect of a single narrow absorption line centered at ω = 0 is shown in Fig. 4 for different values of peak absorption An external file that holds a picture, illustration, etc.
Object name is ukmss-32699-ig0008.jpg and width Ω. In Fig. 4a, An external file that holds a picture, illustration, etc.
Object name is ukmss-32699-ig0008.jpg = −0.05 and Ω = 4 GHz (solid curve, open circles and crosses) and 1 GHz (dashed curve, open squares and triangles). The solid and dashed curves demonstrate the action of complex profile (2), whereas circles (squares) and crosses (triangles) demonstrate the separate action of I(Γ^) and R(Γ^), respectively. One can see, that the profile of perturbed spectrum traces that formed by only I(Γ^) (i.e. it traces the real part of an absorber permittivity). One can say, that the pure phase effect (I(Γ^), circles and squares in Figs. 4a,b) strongly dominates over the pure absorption (R(Γ^), crosses in Figs. 4a,b), like that for the Schrödinger soliton. Such a domination enhances with a lowered An external file that holds a picture, illustration, etc.
Object name is ukmss-32699-ig0008.jpg (Fig. 4b, crosses), however the |An external file that holds a picture, illustration, etc.
Object name is ukmss-32699-ig0008.jpg| growth increases the relative contribution of R(Γ^) and causes the frequency downshift of the bend (Fig. 4a). Amplitude of the bend traces the An external file that holds a picture, illustration, etc.
Object name is ukmss-32699-ig0008.jpg value, while its width is defined by Ω.

The SAM considered above is modeled by the cubicquintic nonlinear term Σ^ in Eq.(1). Such a SAM is typically realized by using the self-focusing inside an active medium (Kerr-lens modelocking). For reliable self-starting operation of mode-locked oscillator it is often desirable to use a suitable saturable absorber (SA), e.g. a semiconductor-based SESAM [23, 24]. Such an absorber can be described in the simplest case by a single-lifetime two-level model, giving the time-dependent loss coefficient Λ (t) as

[partial differential]Λ(t)[partial differential]t=Λ0Λ(t)TsΛ(t)P(t)JsAeff,

where Λ0 is the loss coefficient for a small signal, Aeff is the effective beam area on the SA, Ts and Js are the SA relaxation time and the saturation energy fluency, respectively. Eq. (4) supplements Eq. (1) and the SAM term Σ^ in the latter has to be replaced by (Λ0 – Λ(t)). When the pulse width is longer than the SA relaxation time, one can replace Eq. (4) by its adiabatic solution so that


where ξ [equivalent] Ts/JsS is the inverse saturation power.

We have simulated Eqs. (1,2,4) in the case of Js =50 μJ/cm2 and Ts =0.5 ps, that correspond to the measurement in Fig. 1. Two cases have been considered: weak focusing (Aeff =4000 μm2 or saturation energy Es=2 nJ) and hard focusing (Aeff =1000 μm2 or saturation energy Es=0.5 nJ). We also considered Eq. (5) for the same peak saturation level as weakly focused SA (i.e. ξ−1=4 kW). In the latter case the SA effectively becomes instantaneous, and the perturbed soliton spectrum is the same as for the Kerr-lens modelocking, i.e. as for the cubic-quintic Σ^. When the saturation energy is sufficiently large, there is no difference between the models expressed by Eqs. (4) and (5). The effect of a narrow absorption line is similar to that of the soliton of the cubic-quintic Eq. (1). Decrease of the saturation energy Es causes down-shift of the pulse spectrum as a whole, but the narrow bend on the soliton spectrum reproduces the real part of the absorber permittivity. One can thus conclude that the type of the SAM is irrelevant for an effect of the narrowband absorption lines on a dissipative soliton spectrum.

Another important conclusion from the numerical simulations is the demonstrated stability of the dissipative soliton against perturbations induced by narrowband absorption. In the following analytical treatment we shall therefore omit the stability analysis.


To study the transformation of dissipative soliton spectrum under action of narrow absorption lines, we apply the perturbation method [9, 25]. Since the basic features already become apparent for the perturbed Schrödinger soliton and do not depend on SAM details, we shall consider the simplest case of the cubic nonlinear gain Σ^[P]=κP. The unperturbed solitonic chirp-free solution of such reduced equation with Γ^=0 is a (z, t) = A0sech (t/T0) exp [(t) + iqz] with dϕ/dt=[var pi]. The unperturbed soliton parameters are [25]

[var pi]=0,A02=2ακT02,q=β22T02,T0=ασ,

where the equation parameters are confined


Hence, the soliton wavenumber is q = −γσ/κ.

Its is reasonable to treat the soliton of the reduced Eq. (1) as the Schrödinger soliton with the parameters constrained by the dissipative terms σ, α and κ (see Eqs. (6,7)). This implies that the equation, which has to be linearized with respect to a small perturbation copropagating with the soliton without beating, decay or growth (i.e. having a wavenumber real and equals to q [9]), is the perturbed nonlinear Schrödinger equation

[partial differential]A(z,t)[partial differential]z=iβ22[partial differential]2[partial differential]t2A(z,t)iγP(z,t)A(z,t)+Γ^[A(z,t)].

Linearization of the latter with respect to a perturbation f (t) exp (iqz) results in

iqf(t)=iβ22d2f(t)dt2iγ[2[mid ]a(t)[mid ]2f(t)+a(t)2f*(t)]+Γ^(a+f).

In the spectral domain, Eq. (9) becomes


where [25]

U(ω)=πγT02A02ωcsch(πT0ω/2),S(ω)=iA0πT0cosh(πT0ω/2)l[sm epsilon]l1i(ωωl)/Ωl1+(ωωl)2/Ω2,k(ω)=β22ω2[sm epsilon]l(ωωl)/Ωl+i1+(ωωl)2/Ω2.

Here k (ω) is the frequency-dependent complex wave number, and S (ω) is the perturbation source term for Γ^ corresponding to Eq. (2).

Further, one may assume the phase matching between the soliton and its perturbation. This assumption in combination with the equality U (ω) = U* (ω), which holds for the Schrödinger soliton, results in dωU(ωω)f~*(ω)=dωU(ωω)f~(ω).

The equation (10) for the Fourier image of perturbation is the Fredholm equation of second kind. Its solution can be obtained by the Neumann series method so that the iterative solution becomes [25]


where f~n(ω) is the n-th iteration and f~0(ω)=S(ω)/[k(ω)q].

The “phase character” of a soliton perturbation (i-multiplier in lhs. of Eq. (9) and the expression for the source term (11)) demonstrate that the real part of absorber permittivity contributes to the real part of soliton spectral amplitude. Simultaneously, the resonant condition k (ω) – q =0, which is responsible for a dispersive wave generation caused by, for instance, the higher-order dispersions [9], is not reachable in our case. The resonance can appear in case of large |An external file that holds a picture, illustration, etc.
Object name is ukmss-32699-ig0008.jpg|, κ/γ, and ΩlT0, but such regimes are beyond the scope of this work.

Eq. (12) can be solved numerically. Fig. 5 shows R(f~1(ω)) (dashed curve) and I(f~1(ω)) (dotted curve). One can see, that the real part of absorber permittivity defines R(f~(ω)) while the imaginary part of absorber permittivity defines I(f~(ω)). That agrees with the simulation results and is opposite to the case of a linear pulse propagation. One can also see a tiny frequency down-shift θ of the of f~(ω) minimum from ωl like that in the simulations.

FIG. 5
(Color online) Central part of the dissipative soliton spectra perturbed by f~1(ω) (solid curve) or f~0(ω) (open squares) as well as the profiles of R(f~1) (dashed curve) and I(f~1) (dotted curve) from Eq. (12). Single ...

The pulse spectrum (solid curve in Fig. 5), results from interference of the perturbation with the soliton. For the chosen parameters of the absorption line the zero-order approximation f~0(ω) (open squares) is very close to the first-order approximation (solid curve) but is slightly down-shifted in the vicinity of the bend maximum and minimum.

With even narrower line width Ω of 1 GHz (Fig. 6) the spectral perturbation gets very close to the real part of absorber permittivity and, simultaneously, the spectral down-shift θ (location of the I(f~1) minimum, dashed curve) vanishes. The f~0(ω) (gray solid curves) now perfectly matches the f~1(ω) (open circles and crosses) within a broad range of An external file that holds a picture, illustration, etc.
Object name is ukmss-32699-ig0008.jpg (gray solid curves 1 and 2 as well as circles and crosses belong to An external file that holds a picture, illustration, etc.
Object name is ukmss-32699-ig0008.jpg =−0.005 and −0.05, respectively). The bend amplitudes are in agreement with Eq. (11).

FIG. 6
(Color online) Central part of the dissipative soliton spectra perturbed by f~1(ω) (open circles and crosses) or f~0(ω) (solid gray curves) as well as the profiles of R(f~1) (solid black curve) and I(f~1) (dashed curve) ...

A superposition of three identical absorption lines, which corresponds to the numerical spectra in Fig. 2, is shown in Fig. 7. One can see, that the lowest-order analytical solution f~0(ω)accurately reproduces the numerical result. It is important, that a cumulative contribution of lines into k (ω) does not distort a superposition contribution of S (ω) into a soliton spectrum (see Eqs. (11,12)). This means that the individual contribution of a single line within a group is easily distinguishable and can be quantitatively assessed, opening way for interesting spectroscopic applications.

FIG. 7
(Color online) Central part of the dissipative soliton spectra perturbed by f~0(ω) (black solid curve) and the profiles of R(f~0) (dashed curve) and I(f~0) (dotted curve) from Eq. (12). Triplet of absorption lines is centered at ...

As Figs. Figs.55 and and66 suggest, the zero-order approximation f~0(ω)=S(ω)/[k(ω)q] is quite accurate for a description of perturbation in the limit of |An external file that holds a picture, illustration, etc.
Object name is ukmss-32699-ig0008.jpg| << 1. This allows expressing the perturbed spectrum of an isolated line (see (11,12)) in analytical form [25]:

P~(ω)[equivalent][mid ]a~(ω)+f~0(ω)[mid ]2A02π2T02sech(πT0ω2)2(β22ω2+q)2(1+(ωωl)2Ωl2)[[sm epsilon]l2+2(ωωl)[sm epsilon]lΩl(β22ω2+q)+(β22ω2+q)2(1+(ωωl)2Ωl2)].

Eq. (13) allows further simplification in the case of |An external file that holds a picture, illustration, etc.
Object name is ukmss-32699-ig0008.jpg| << 1

P~(ω)A02π2sech(ασπω2)2ασ[1+2[sm epsilon]lκγ(αωl2+σ)(ωωl)(1+(ωωl)2Ωl2)Ωl],

where Eqs. (6) and the condition αΩl2[dbl less than]1 have been used.

Eq. (14) demonstrates that the spectral bend follows the real part of absorber permittivity. Spectral down-shift of bend is the effect of O([sm epsilon]2) and not included in (14). The perturbation is represented by the term in square brackets and its relative amplitude is proportional to An external file that holds a picture, illustration, etc.
Object name is ukmss-32699-ig0008.jpg. Furthermore, the aspect ratio of the kink grows with i) the increase of the relative contribution of the SAM κ/γ ; ii) the gain bandwidth 1/α; iii) approaching of the resonance frequency ωl to the center of soliton spectrum (but the ratio of the aspect ratio to the local soliton spectral power increases with |ωl|, because the former decreases as ωl2 while the latter falls faster as cosh (πT0ωl/2)2); and iv) with an approaching to the soliton stability border, which corresponds to vanishing σ. It should be noted, that smaller σ entails the soliton width growth (Eq. (6)).

Since the soliton parameters are interrelated, it is instructive to express σ through the observable parameters such as soliton energy E or the soliton width T0. When αωl2[dbl less than]σ (e.g. ωl ≈0 or/and an oscillator operates far from the stability border σ =0), the perturbation amplitude is inversely proportional to γκE2:

2[sm epsilon]lκγ(αωl2+σ)32[sm epsilon]lαγκE2=2[sm epsilon]lκT02γα

For a fixed gain bandwidth, the amplitude scales with squared pulsewidth T02. Ultimately, the latter equation is equivalent to

2[sm epsilon]lκγ(αωl2+σ)2[sm epsilon]lq,

i.e. the relative perturbation amplitude near soliton central frequency is the ratio of incurred loss coefficient to the soliton wavenumber, regardless of the z coordinate normalization. Therefore, this analytical expression that has been derived for the self-consistent oscillator, should also be valid for the case of a soliton propagation in a long fiber when the conditions of applicability |An external file that holds a picture, illustration, etc.
Object name is ukmss-32699-ig0008.jpg| << 1, Ω << 1/T0 are met. The final form of the soliton spectrum thus becomes

P~(ω)=P~0(ω)(12ql[sm epsilon]l(ωωl)/Ωl1+(ωωl)2Ωl2),

where P~0(ω) is the spectrum of an unperturbed soliton.


In the analysis above we have shown that, for the case of sufficiently sparse, narrow and weak Lorentzian absorber lines, their spectral signatures are equivalent to the dispersion-like modulation with a relative amplitude equal to the peak absorption coefficient over oscillator round-trip (or nonlinear length for passive propagation) divided by the soliton wavenumber. For quantitative comparison with the experiment we recall equation (6) and express the maximum spectrum deviation of a single line at |ωωl|; = Ωl through observable parameters:

[mid ][sm epsilon]lq[mid ]=χlLT02[mid ]β2[mid ]=χlL0.0319[mid ]β2[mid ](Δν)2,

where χl = 2An external file that holds a picture, illustration, etc.
Object name is ukmss-32699-ig0008.jpgl is the peak absorption coefficient of the line, L is the absorber path length and Δν is the full width at half maximum of the soliton spectrum. Substituting the actual values of the setup in Fig. 1 (β2 = −820 ± 40 fs2, Δν =113 cm−1= 3.39 THz, round-trip air path length L = 149 cm, relative humidity 50 ± 1% at 21 ± 0.5 °C) and taking e.g. the line at 4088 cm−1 (122.56 THz, marked with an asterisk), we obtain |An external file that holds a picture, illustration, etc.
Object name is ukmss-32699-ig0008.jpg/q| = 0.33 ± 0.02 for the maximum modulation, which is in perfect agreement with the observed value of 31.5% (Fig. 1b, black line). The agreement is remarkably accurate given the less than optimal resolution of the spectrometer (0.25 cm−1) and significant third-order dispersion of about +104 fs3 that was not accounted for in the presented analysis.

It is important to notice, that the expression (18) includes only the externally observable soliton band-width and relatively stable dispersion parameter. The alignment-sensitive values like saturated losses σ, nonlinearity γ, nonlinearity saturation parameter κ, etc. which are in practice not known with sufficient accuracy, are all accounted for by the self-consistent soliton parameters.

Another important point is the fact, that the signal amplitude 2|An external file that holds a picture, illustration, etc.
Object name is ukmss-32699-ig0008.jpg/q| can be much bigger than that from conventional absorption spectroscopy χlL of the cell with the same length. The signal enhancement factor can be controlled by the pulse parameters and it exceeds an order of magnitude for the presented case (χlL = 5% for selected line and a single-pass cell of a resonator size). For additional sensitivity improvement one can apply the well-developed intracavity multi-pass cell technique [26]. The expression (18) suggests, that ultimate sensitivity can be obtained at the expense of the reduced bandwidth coverage Δν. In this respect, the presented technique has the same quadratic dependence of sensitivity on spectral bandwidth as the conventional intracavity absorption spectroscopy [5].

Further refinement of the presented theory should include demonstration of its applicability to arbitrarily shaped absorption features. The superposition property provides a strong argument for such extension, but it has to be rigorously proven for Doppler- and more general Voigt-shaped lines, and also for the dense line groups in e.g. Q-branches. It would be interesting also to extend the theory to the absorber lines at the soliton wings (Fig. 1a).

With the above issues resolved, the soliton-based spectroscopy may become a powerful tool for high-resolution, high-sensitivity spectroscopy and sensing. Possible implementation include soliton propagation in a gas-filled holey fibers, as well as already presented intracavity spectroscopy with femtosecond oscillators. The latter, being a natural frequency comb source, allows direct locking to optical frequency standards, providing for ultimate resolution and spectral accuracy.


We have been able to derive an analytical solution to the problem of an one-dimensional optical dissipative soliton propagating in a medium with narrowband absorption lines. We predict appearance of spectral modulation that follows the associated index of refraction rather than absorption profile. The novel perturbation analysis technique is based on integral representation in the spectral domain and is insensitive to the diverging differential terms, inherent to the Taylor series representation of the narrow spectral lines.

The model is applicable to a conventional soliton propagation and to a passively modelocked laser with intracavity absorber, the only difference being the characteristic propagation distance (dispersion length and cavity round-trip, respectively). In the latter case the prediction has been confirmed for a case of water vapour absorption lines in a mid-IR Cr:ZnSe oscillator. The model provides very good qualitative and quantitative agreement with experimental observations, opening a way to metrological and spectroscopical applications of the novel technique, which can provide a significant (order of magnitude and more) enhancement of the signal over conventional absorption for the same cell length.


We gratefully acknowledge insightful discussions and experimental advice from N. Picqué, G. Guelachvili (CNRS, Univ. Paris-Sud, France), and I. T. Sorokina (NTNU, Norway). This work has been supported by the Austrian science Fund FWF (projects 17973 and 20293) and the Austrian-French collaboration Amadée.


[1] Sorokin E, Sorokina IT, Mandon J, Guelachvili G, Picqué N. Opt. Express. 2007;15:16540. [PMC free article] [PubMed]
[2] Mandon J, Guelachvili G, Sorokina I, Sorokin E, Kalashnikov V, Picqué N. Europhysics Conference; Abstract Volume 32G. paper WEoB.4 at Europhoton 2008.
[3] Kalashnikov V, Sorokin E, Mandon J, Guelachvili G, Picqué N, Sorokina IT. paper TUoA.3 at Europhoton 2008. Abstract Volume 32G.
[4] Sorokina IT, Sorokin E, Carrig T. paper CMQ2 at CLEO/QELS 2006; Technical Digest on CD.
[5] Baev VM, Latz T, Toschek PE. Appl. Phys. B. 1999;69:171.
[6] Picqué N, Gueye F, Guelachvili G, Sorokin E, Sorokina IT. Opt. Lett. 2005;30:3410. [PMC free article] [PubMed]
[7] The HITRAN. database
[8] Haus HA, Fujimoto JG, Ippen EP. Structures for additive pulse mode locking. JOSA B. 1991;8:2068–2076.
[9] Akhmediev NN, Ankiewicz A. Solitons: Nonlinear Pulses and Beams. Chapman&Hall; 1997.
[10] Aranson IS, Kramer L. Rev. Mod. Phys. 2002;74:99.
[11] Brabec T, Krausz F. Phys. Rev. Lett. 1997;78:3282.
[12] Kärtner FX, editor. Few-cycle Laser Pulse Generation and its Applications. Springer Verlag; Berlin: 2004.
[13] Akhmediev NN, Ankiewicz A, editors. Dissipative Solitons. Springer Verlag; Berlin: 2005.
[14] the signs before i in Eq. (1) correspond to those in [8] and Akhmanov SA, Vysloukh VA, Chirkin AS. Optics of Femtosecond Laser Pulses. AIP; NY: 1992.
[15] Agrawal GP. Nonlinear Fiber Optics. Academic Press; San Diego: 2001.
[16] Biswas A, Konar S. Introduction to non-Kerr Law Optical Solitons. Chapman&Hall; Boca Raton: 2007.
[17] Haus HA, Silberberg Y. J. Opt. Soc. Am. B. 1985;2:1237.
[18] Kalashnikov VL, Podivilov E, Chernykh A, Apolonski A. Applied Physics B. 2006;83:503.
[19] Butylkin VS, Kaplan AE, Khronopulo Yu.G., Yakubovich EI. Resonant Nonlinear Interactions of Light with Matter. Springer Verlag; Berlin: 1989.
[20] Oughstun KE. Electromagmetic and Optical Pulse Propagation 1. Springer; NY: 2006.
[21] De Sousa Meneses D, Gruener G, Malki M, Echegut P. J. Non-Crystalline Solids. 2005;351:124.
[22] Yamaoka Y, Zeng L, Minoshima K, Matsumoto H. Appl. Opt. 2004;43:5523. [PubMed]
[23] Islam MN, Sunderman ER, Soccolich CE, Bar-Joseph I, Sauer N, Chang TY, Miller BI. IEEE J. Quantum Electron. 1989;25:2454.
[24] Keller U, Weingarten KJ, Kärtner FX, Kopf D, Braun B, Jung ID, Fluck R, Honninger C, Matuschek N, Aus der Au J. IEEE J. Selected Topics in Quantum Electron. 1996;2:435.
[25] Kalashnikov VL. Maple worksheet (unpublished)
[26] Kowalevicz AM, Sennaroglu A, Zare AT, Fujimoto JG. J. Opt. Soc. Am. B. 2006;23:760.