Parameters for litharge SF56 glasses.
In this contribution some basic properties of femtosecond laser pulse are summarized. In sections 2.1–2.5 the generation of femtosecond laser pulses via mode locking is described in simple physical terms. In section 2.6 we deal with measurement of ultrashort laser pulses. The characterization of ultrashort pulses with respect to amplitude and phase is therefore based on optical correlation techniques that make of the short pulse itself. In section 3 we start with the linear properties of ultrashort light pulses. However, due to the large bandwidth, the linear dispersion is responsible for dramatic effects. To describe and manage such dispersion effects a mathematical description of an ultrashort laser pulse is given first before we continue with methods how to change the temporal shape via the frequency domain. The chapter ends with a paragraph of the wavelet representation of an ultrashort laser pulse.
- femtosecond laser pulse
- ultrashort pulse
- spectral phase
Propagation of ultrashort optical pulses in a linear optical medium consisting of free space [1, 2, 3, 4, 5], dispersive media [6, 7], diffractive optical elements [8, 9, 10], focusing elements  and apertures [12, 13] has been extensively studied analytically, though only a few isolated attempts have been made on numerical simulation. However, analytical methods have the limitations of not being able to handle arbitrary pulse profiles.
2. Ultrashort laser pulses generations
The central aim of this section is to give a concise introduction to nonlinear optics and to provide basic information about the most-widely used tunable femtosecond laser sources, in particular tunable Ti:sapphire oscillators and Ti:sapphire amplifiers or optical parametric amplifiers.
2.1 Titane sapphire oscillator
In 1982, the first Ti:sapphire laser was built by Mouton . The laser tunes from 680 nm to 1130 nm, which is the widest tuning range of any laser of its class1. Nowadays Ti:sapphire lasers usually deliver several watts of average output power and produce pulses as short as 6.5 fs (Figure 1) .
At high intensities, the refractive index depends nonlinearly on the propagating field. The lowest order of this dependence can be written as follows:
n0: linear index refractive.
where is the nonlinear index coefficient and describes the strength of the coupling between the electric field and the refractive index . The intensity is:
The refractive index changes with intensity along the optical path and it is larger in the center than at the side of the nonlinear crystal. This leads to the beam self-focusing phenomenon, which is known as the Kerr lens effect (see Figure 2).
Consider now a seed beam with a Gaussian profile propagating through a nonlinear medium, e.g. a Ti:sapphire crystal, which is pumped by a cw radiation. For the stronger focused frequencies, the Kerr lens favors a higher amplification. Thus, the self-focusing of the seed beam can be used to suppress the cw operation, because the losses of the cw radiation are higher. Forcing all the modes to have equal phase (mode-locking) implies that all the waves of different frequencies will interfere (add) constructively at one point, resulting in a very intense, short light pulse. The pulsed operation is then favored, and it is said that the laser is mode-locked. Thus, the mode-locking occurs due to the Kerr lens effect induced in the nonlinear medium by the beam itself and the phenomenon is known as Kerr-lens mode-locking.
Auto TPL Tripler for laser oscillator.
Sprite XT: Tunable ultrafast Ti: sapphire laser (Figure 3).
The modes are separated in frequency by , L being the resonator length, which also gives the repetition rate of the mode-locked lasers:
Where L is length of cavity and T is period.
Moreover, the ratio of the resonator length to the pulse duration is a measure of the number of modes oscillating in phase. For example, if and the emerging pulses have time duration, there are 105 modes contributing to the pulse bandwidth. There are two ways of mode-locking a femtosecond laser: passive mode-locking and active mode-locking. In a laser cavity, these modes are equally spaced (with spacing depending on the cavity length). The electric field distribution with N such modes in phase (considered to be zero, for convenience) can be written as:
Where is the central frequency and is the mode spacing, this appears as a carried wave with frequency domain.
n: is integer from 1 to N.
The laser intensity is given by
Where E(t): electrical field.
This is series of pulses with width inversely proportional to the number of modes that are locked in phase of the mode spacing. The concept of mode-locking is easier said than done.
Figure 4. shows how the time distribution of a laser output depends upon the phase relations between the modes. Figure 4a is the resultant intensity of two modes in phase Figure 4b, is the resultant intensity of five modes in phase and a period repetition of a wave packet from the resultant constructive interference can be seen.
2.2 CPA laser system
CPA is the abbreviation of chirped pulse amplification. Chirped pulse amplification is a technique to produce a strong and at the same time ultrashort pulse. The concept behind CPA is a scheme to increase the energy of an ultrashort pulse while avoiding very high peak power in the amplification process.
In the CPA technique, ultrashort pulses are generated typically at low energy ∼10−9 J, with a duration around 10−12 –10−14 seconds and at a high repetition rate of about 108 1/s in an oscillator (Figure 5).
2.3 Multipass and regenerative amplification
Two of the most widely used techniques for amplification of femtosecond laser pulses are the multipass and the regenerative amplification. In the multipass amplification different passes are geometrically separated (see Figure 6a).
The regenerative amplification technique implies trapping of the pulse to be amplified in a laser cavity (see Figure 2b). Here the number of passes is not important. The pulse is kept in the resonator until all the energy stored in the amplification crystal is extracted. Trapping and dumping the pulse in and out of the resonator is done by using a Pockelcell (or pulse-picker) and a broad-band polarizer .
Of all potential amplifier media, titanium- doped sapphire has been the most wide spread used. It has several desirable characteristics which make it ideal as a high-power amplifier medium such as a very high damage threshold (∼8
By using a dispersive line (combination of gratings and/or lenses), the individual frequencies within a femtosecond pulse can be separated (stretched) from each other in time (see Figure 7a).
In normal materials, low frequency components travel faster than high frequency components; in other words, the velocities of large wavelength components are higher than that of shorter ones. These materials induce a positive group velocity dispersion on a propagating pulse. To compensate the positive GVD (Group Velocity Dispersion) and rephase the dephased components a setup which produces negative group velocity dispersion is needed .
2.4.1 Grating compressor
Four identical gratings in a sequence as shown in Figure 8 make up a grating compressor. A pulse impinges on the first gratings with an angle of . Then from the second grating the spectral components in the spectrum travel together in parallel directions but with wavelength dependent position (spatially chirped). The gratings are set in such a way that their wavelength dispersions are reversed which implies that the exiting ray from the second grating is parallel to the incident ray to the first grating.
The group delay induced by the grating compressor is given
where, is the light wavelength,
GDD is Group Delay dispersion.
The third order dispersion produced in the grating compressor will be:
TOD is Third Order Dispersion.
2.4.2 Prism compressor
A prism compressor is built of four sequentially arranged identical prisms used in a geometry similar to Figure 9; often at their minimum deviation (to decrease geometrical (spatial) distortion of prisms on the beam) and in their Brewster angle (to minimize power loss). Because of the symmetry in the arrangement, it is possible to place a mirror after the second prism (as we did in the grating compressor setup) perpendicular to the beam propagation direction. The first prism spreads the pulse spectral components out in space. In the second prism the red frequencies of the spectrum must pass through a longer length in the glass than the blue frequencies.
The optical path of a ray propagating in the compressor is defined as : GDD will result
: is angle.
The third order dispersion can also be evaluated in the same manner to that used above
n: is refractive index.
2.5 Mathematical description of laser pulses
In order to understand the behavior of ultrashort light pulses in the temporal and spectral domain, it is necessary to formulate the relation between the two domains mathematically. It is important to introduce the concept of the amplitude and the phase of the electric field because the generation, measurement, and shaping of ultrashort laser pulses is based on measuring and influencing these properties. The electric field in the time-domain is invariably connected with its counterpart in the frequency-domain via a Fourier transform:
2.5.1 Pulse duration and spectral width
The statistical definitions are usually used in theoretic calculations and given as
In case of Gaussian pulse is easy to determine pulse duration and spectral width by applying FWHM (Full Width Half at Maximum) of intensity. One can show that these quantities are related through the following universal inequality.
Therefore, one defines the pulse duration as the Full Width at Half Maximum (FWHM) of the intensity profile and the spectral width as the FWHM of the spectral intensity. The Fourier inequality is then usually given by
where K is a numerical constant, depending on the assumed shape of the pulse.
22.214.171.124 Gaussian pulse
The Gaussian pulse, which is most commonly used in ultrashort laser pulse characteristics. The pulse is linearly chirped and represented by
A0: amplitude, : chirp,
: pulse duration at FWHM.
:pulse duration at FWHM after propagation.
The instantaneous frequency is given as
W0: central pulsation.
: instantaneous pulse.
is instantaneous phase.
The shortest possible pulse, for a given spectrum, is known as the
2.5.2 Time domain description
Since in this paper the main emphasis is on the temporal dependence, all spatial dependence is neglected, i.e., the electric field , is a real quantity and all measured quantities are real. However, the mathematical description is simplified if a complex representation is used:
where is the complex envelope, usually chosen such that the real physical field is twice the real part of the complex field, and is the carrier frequency, usually chosen to the center of the spectrum. In this way the rapidly varying is separated from the slowly varying envelope . can be further decomposed into:
is often to as the temporal phase of the pulse and the absolute phase, which relates the position of the carrier wave to the temporal envelope of the pulse (see Figure 11). In the strong linear term due to the carrier frequency, is omitted.which means that a nonlinear temporal phase yields a time-dependent frequency modulation- the pulse is said to carry a chirp (illustrated in Figure 11b).
An ultrashort pulse of light will lengthen after it has passed through glass as the index of refraction, which dictates the speed of light in the material, depends nonlinearly on the wavelength of the light. The wavelength of an ultrashort pulse of light is formed from the distribution of wavelengths either side of the center wavelength with the width of this distribution inversely proportional to the pulse duration.
126.96.36.199 Phase and chirp
Instantaneous phase function of can be described as the sum of temporal phase andproduct of carrier frequency with time by the relation
Carrier frequency has been choosen by minimizing of temporal variationof phase . The first deriviation of is defined by temporally-dependentcarrier frequency as the result of applying the derivation we receive relation expended in series. Then carrier frequency time denotes quadratic chirp.
Positive chirp is when leading edge of pulse is red-shifted in relation to central wavelength and trailing edge is blue-shifted. Negative chirp happens in opposite case. Linear chirp, instantaneous frequency varies linearly with time. The presence of chirp results in significant different delays between the spectrally different components of laser pulse causing pulse broadening effect and leading to a duration-bandwidth.
Chirps always appear when ultrashort laser pulses propagate through a medium such as air or glass, where the spectral components of the pulse are subject to a different refractive index. This effect is called
2.5.3 Lens frequency domain description
The frequency representation is obtained from the time domain by a complex Fourier transform,
Just as in the time domain, can be written as:
where now denotes the spectral phase. An inverse transform leads back to the time domain,
From Eq. (20) it is clear that can be seen as a superposition of monochromatic waves. A common procedure is to employ Taylor expansion
is spectral phase, w is pulsation, and is variation of the spectral phase with pulsation.
It can be seen by inserting this Taylor expansion into Eq. (21) that the first, two terms will not change the temporal profile of the pulse.
2.6 Ultrashort pulse measurement techniques
The most commonly used technique of ultrashort laser pulse width measurement is concerned with the study of the temporal intensity profile through its second-order correlation function that is obtained by the second-harmonic generation. Ultrashort - pulse characterization techniques, such as the numerous variants of frequency resolved optical gating (FROG) and spectral phase interferometry for direct electric-field reconstruction (SPIDER), fail to fully determine the relative phases of wellseparated frequency components. If well-separated frequency components are also characterized gate pulses are used .
2.6.1 Non-interferometric techniques
188.8.131.52 Intensity autocorrelation
Complex electric field corresponds to intensity and an intensityautocorrelation function is defined by
is quadrature detection
I(t) is intensity and E(t) is electrical field
Two parallel beams with a variable delay are generated, then focused into a second-harmonic-generation crystal to obtain a signal proportional to . Only the beam propagating on the optical axis, proportional to the cross product is retained. This signal is then recorded by a slow detector, which measures
is exactly the intensity autocorrelation .
This numerical factor, which depends on the shape of the pulse, is sometimes called the deconvolution factor. If this factor is known, or assumed, the time duration (intensity width) of a pulse can be measured using an intensity autocorrelation. However, the phase cannot be measured [17, 18].
The technique of Frequency-Resolved Optical Gating (FROG) has been introduced by Trebino and coworkers. In FROG technique signal has been temporally shifted about τ through time-delay element in respect with signal . Then, two signals have been in nonlinear medium non-interferometrically overlapped. As the result of SFG or DFG process (at the efficient phase matching conversion) one receive the FROG signal.
In construction with process II also collinear geometry is possible.
The spectral intensity
is delay between two temporal electrical field E1(t) and E2(t)
is delay between two spectral electrical field E1(w) and E2(w)
is called FROG-trace. Relations (25), (26) are only the mathematical representation of Eq. (24). The task of receiving unknown complex signals and from measured FROG trace is known as the FROG reconstruction problem.
The resulting trace of intensity versus frequency and delay is related to the pulse’s spectrogram, a visually intuitive transform containing both time and frequency information. Using phase retrieval concepts that the FROG trace yields the full intensity and phase of an arbitrary ultrashort pulse. As has been already mentioned, several schemes and methods exist for frequency-resolved optical gating as a technique for the full characterization of ultrashort optical signals as complex electric fields [14, 19].
Experimental setup. PARS microscopy with 532-nm excitation and 1310-nm integration beams. BC, beam combiner; GM, galvanometer mirror; L, lens; OL, objective lens; PBS, polarized beam splitter; PD, photodiode; QWP, quarter wave plate; SMF, single mode fiber.
2.6.2 Interferometric techniques
184.108.40.206 Interferometric autocorrelation
Setup for an interferometric autocorrelator is similar to the field autocorrelator above, with the following optics added: L: converging lens, SHG: secondharmonic generation crystal, F: spectral filter to block the fundamental wavelength. A nonlinear crystal can be used to generate the second harmonic at the output of a Michelson interferometer in a collinear geometry. In this case, the signal recorded by a slow detector (Figure 12).
is called the interferometric autocorrelation. It contains some information about the phase of the pulse: the fringes in the autocorrelation trace wash out as the spectral phase becomes more complex (Figure 13).
The Spectral Phase Interferometry for Direct Electric-field Reconstruction technique (SPIDER) is based on spectral interferometry and needs no components which has to be shifted over the measurement process. From the signal that should be characterized, the copy-signal is being generated by beam splitter
The time between the signal and copy itself has been established through fixed position at the optical delay-line. Then the copy of signal goes through phase filter (dispersive medium, for instance SF10 glass), so arises the signal
EM electrical field at point M
: phase of electrical field
Through the phase filter electric field get additional spectral phase , which corresponds to temporal extension . From the signals
SFG-Signal can be created
As square law detectors are not sensitive to the phase, the measurement of the intensity (whether it is spatial or spectral) is an easy task but the measurement of the phase needs indirect solutions (Figure 14).
Spectral interferometry allows us obtain difference between two spectral phases. To the spectral interferometry spectrum, one should apply fast Fourier transform and as the product one will achieve in form of one center peak and two sidebands lower peaks in the time domain (Figure 15).
Centered peak contains only spectrum information. One filter out two peaks and from existing one can receive spectral phase difference by applied inverse fast Fourier transform. To the main advantages of the SPIDER method can be counted following properties: pulse retrieval is direct (non-iterative), minimal data are required: only one spectrum yields spectral phase .
3. Propagation of a light pulse in a transparent medium
The frequency Fourier transform of a Gaussian pulse has already been given as
is coefficient of function Gaussian
An ultrashort Fourier limited pulse has a broad spectrum and no chirp; when it propagates a distance through a transparent medium, the medium introduces a dispersion to the pulse inducing an increase in the pulse duration. To investigate and determine the dispersion, we assume a Gaussian shape for the pulse. The electric field of the pulse is given as Eq. (35).
After the pulse has propagated a distance z, its spectrum is modified to
where is now a frequency-dependent propagation factor. In order to allow for a partial analytical calculation of the propagation effects, the propagation factor is rewritten using a Taylor expansion as a function of the angular frequency, assuming that (this condition is only weakly true for the shortest pulses). Applying the Taylor expansion to Eq. (37), the pulse spectrum becomes.
is bandwidth of the pulse and is central pulsation
The time evolution of the electric field in the pulse is then derived from the calculation of the inverse Fourier transform of Eq. (39),
In the first exponential term of Eq. 40, it can be observed that the phase of the central frequency is delayed by an amount after propagation over a distance z. Because the phase is not a measurable quantity, this effect has no observable consequence. The phase velocity measures the propagation speed of the plane wave components of the pulse in the medium. These plane waves do not carry any information, because of their infinite duration.
The second term in Eq. 40 shows that, after propagation over a distance z, the pulse keeps a Gaussian envelope. This envelope is delayed by an amount , being the group velocity. The second term in Eq. 40 also shows that the pulse envelope is distorted during its propagation because its form factor , defined as
Depends on the angular frequency through
This term is called the “Group Velocity Dispersion”. The temporal width of the pulse at point z:
initial pulse before propagation inside the medium
3.1 Application in litharge index SF56 medium
In optical materials, the refractive index is frequency dependent. This dependence can be calculated for a given material using a Sellmeier equation, typically of the form
Bi: is the coefficients of glass see in Table 1
The index of litharge SF56 is given by the following expression (43):
where is the frequency of resonance and is the amplitude of resonance. In the case of optical fibers, the parameters and are obtained experimentally by fitting the measured dispersion curves to Eq. (44) with and depend on the core constituents .
3.2 Parameter of dispersion
An ultrashort Fourier limited pulse has a broad spectrum and no chirp; whenit propagates a distance through a transparent medium, the medium introduces dispersion to the pulse inducing an increase in the pulse duration. We consider dispersions of orders two. The pulse broadens on propagation because of group velocity dispersion (GVD).
In summary, the propagation of a short optical pulse through transparent medium results in a delay of the pulse, a duration broadening and a frequency chirp. All these phenomena are increase with distance of propagation. We shown in Figure 16 that the duration broadening is not linear for ultrashort pulse. Specially under 70 fs or less. The Eq. (43) is not applicable for pulse less than 70 fs. For minimize these parameters we introduce the nonlinear phenomena as Self phase modulation (SPM), Soliton pulse, dispersion compensate fiber.
3.3 Group velocity dispersion
The Group Velocity Dispersion (GVD) is defined as the propagation of different frequency components at different speeds through a dispersive medium. This is due to the wavelength-dependent index of refraction of the dispersive material .
Itseemsto methatwe can write asarecurrence, giving basedonderivativesoforderi, the index of refraction. Matrix form, wecan writet
The various terms of the Taylor expansion to order n can be written in the shape of a matrix [A], which’s we can express various terms A ij.
is Taylor series of phase and is the orders of Taylor series
Analytically known and experimentally observed propagation affects such as spectral shift, pulse broadening and asymmetry in dispersive media can be easily brought out in the simulation using formalism presented here. In addition, such studies can be extended to pulses of arbitrary temporal shape without any further algorithmic complexity by numerical simulation. Higher order dispersion effects can be handled easily in the numerical simulation unlike in the case of analytical calculation (Figure 17) .
4. Time-frequency decomposition
4.1 Wavelet theory
The wavelets are very particular elementary functions, these are the shortest vibrations and most elementary that one can consider. One can say that the wavelet east carries out a zooming on any interesting phenomenon of the signal which place on a small scale in the vicinity of the point considered .
4.2 Wavelet techniques
Starting with a signal , in plane , we define wavelet centered at Ω by (Figure 18):
TF it is the Fourier Transform equation
z is the distance of propagation
In time, the pulse is also Gaussian, of parameter
The maximum of amplitude of the wavelet vary with Ω, center frequency of analysis on Gaussian of parameter .
The signal propagates in the positive direction in a linear dispersive and transparent medium, which fills the half space and whose refractive index is . After propagation, the wavelet may be written as
As already mentioned, is large enough to ensure that analyzing function has only non negligible values over a spectral range lying in the neighborhood of Ω in (Figure 19). Under these circumstances, we have
Neglecting the higher terms in Eq. (49):
We calculates the temporal electric field associated with the wavelet
The amplitude of the incident Ω wavelet is given from Eq. (58) by
This wavelet is characterized by a Gaussian envelope. This decomposition is valid only for the values of much larger than ().
The delay of group of the wavelet is characterized by a Gaussian envelope which is the temporal width.
The delay of group of the wavelet is inversely proportional to the velocity of group its envelope propagates without deformation .
4.3.1 Parameters of the simulations
Pulse of the wavelet:
Longer of the medium:
To describe the propagation of the pulse, we only consider the propagation of the maximum of each wavelet in a three dimensional representation:
Generating ultrashort light pulses requires a laser to operate in a particular regime, called mode-locking which many be illustrated either in the frequency or the time domain. Depending on the particular case, one description is much more intuitive than the other and we have chosen to present the simpler approach. The generation of femtosecond laser pulses via mode locking is described in simple physical terms. As femtosecond laser pulses can be generated directly from a wide variety of lasers with wavelengths ranging from the ultraviolet to the infrared no attempt is made to cover different technical approaches.
The ability to accurately measure ultrashort laser pulses is essential to creating, using, and improving them, but the technology for their measurement has consistently lagged behind that for their generation. The result has been a long and sometimes quite painful history of attempts—and failures—to measure these exotic and ephemeral events. The reason is that many pulse-measurement techniques have suffered from, and continue to suffer from, a wide range of complications, including the presence of ambiguities, insufficient temporal and/or spectral resolution and/or range, an inherent inability to measure the complete pulse intensity and/or phase, an inability to measure complex pulses, and misleading results due to the loss of information due to idiosyncrasies of the technique or multishot averages over different pulses.
Finally, we have demonstrated here the possible decomposition of an ultrashort pulse into an infinite number of longer Fourier transform limited wavelets which propagate without any deformation through a dispersive medium. After propagation through the medium, the pulse may be visualized in a three-dimensional representation by the locus of the wavelet maxima.
Florenţa A.C, Dynamics of Ultra-short Laser Pulse Interaction with Solids at the Origin of Nanoscale Surface Modification” geboren am 16.11.1971 in Berca, Rumänien Cottbus (2006)
Chritov, I.P., Propagation of femtosecond light pulses, Opt. Commun. 53, 364–366.
Sheppard, C.R, Xianosong, G, Free-space propagation of femtosecondligthpulses, Opt. Commun.133, 1–6. (1997)
Agrawal,G.P, Spectrum-induced changes in diffraction of pulsed beams, Opt. Commun.157, 52–56. (1998).
Kaplan, A.E., Diffraction-induced transformation of near cycle and subcycle pulses, J . Opt. Soc. Am. B 15, 951–956. (1998)
Agrawal, G.P, Far-Field diffraction of pulsed optical beams in dispersive media, Opt. Commun.167, 15–22. (1999).
Porras,M.A, Propagation of single-cycle pulse light beams in dispersive media, Phys. Rev. A60, 5069–5073. (1998).
Ichikawa, H, Analysis of femtosecond-order optical pulses diffracted by periodic structure, J. Opt. Soc. Am. A 16, 299–304. (1999).
Pietstun, R., Miller, D.B., Spatiotemporal control of ultrashort optical pulses by refractive-diffractive-dispersive structured optical elements, Opt. let. 26, 1373–1375. (2001).
Kempe, M., Stamm, U., Wilhelmi, B., and Rudolph, W., Spatialand temporaltransformation of femtosecond laser pulses by lenses systems, J . Opt. Soc. 38, 1058–1064. (1999).
Matei,G.O,Gili, M.A., Spherical aberration in spatial and temporal transformation of femtosecond laser pulses by lenses and lens systems, J. Opt. Soc. Am. B9, 1158–1165. (1992).
Fuchss,U., Zeintner, U.D., and Tunnermenn (A.),Ultrashort pulse propagation in complex optical systems , Opt. Express 13, 3852–3861. (2005).
Leonid.S., Ferrari.A, and Bertolotti.M, Diffraction of a time Gaussian-shaped pulsed plane wave from a slit, Pure. Appl. Opt5, 349–353. (1996).
Moulton.P.f, Quant. Electron. Conf., Munich, Germany, June (1982).
Holzwarth,“Measuring the Frequency of Light using Femtosecond Laser Pulses” aus Stuttgart den 21. December (2000).
Rulliere, C, Femtosecond laser pulses: principles and Experiments, Springer,(1998).
Joanna.M, “Seutp of very advanced for phase and amplitude reconstruction of electric field (VAMPIRE)“Master Thesis submitted to the Faculty for the Natural Sciences and for Mathematics of the Rostock University Germany (2007).
Abdolah.M.K, “Manipulation and characterization of femtosecond laser pulses for cluster spectroscopy” vorgelegt von aus Teheran (Iran) July (2007).
Greg.T, Andy. R Margaret.M, Measurement of 10 fs Laser Pulses. 1077-260X,IEEE (1996).
Hofmann.A, Bestimmung des elektrischen Feldesultrakurzer Laserpulsemit SPIDER, Teil1: Theorie, Projektpraktikumsbericht, PhysikalischesInstitut, EP1, Universität Würzburg (2005).
Jong Kook K. Investigation of high nonlinearity glass fibers for potential applications in ultrafast nonlinear fiber devices, Dissertation submitted to the Faculty of the Virginia (2005).
Khelladi.M, Seddiki.O, and Bendimerad.T, Time-Frequency Decomposition of an Ultrashort Pulse:Wavelet Decomposition, Radioengineering ISSN, April (2008) Vol.17, pp.1210–2512, N1
Meyer,Y, Jaffort, S., Riol,O, waveletanalysis, edition Française de scientifique American(1987)