Next Article in Journal
Active Optics–Freeform Segment Mirror Replications from a Deformable Matrix
Previous Article in Journal
Adaptive Detection of Wave Aberrations Based on the Multichannel Filter
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Modeling for Generating Femtosecond Pulses in an Er-Doped Fiber Using Externally Controlled Spectral Broadening and Compression Mechanisms

by
Mohamed Hemdan Abo-elenein
1,*,
Salah Eldeen Ibrahim Hassab Elnaby
1,
Amin Fahim Hassan
2 and
Afaf Mahmoud Abd-Rabou
2
1
Department of Engineering Applications of Lasers, National Institute of Laser Enhanced Sciences, Cairo University, Giza 12613, Egypt
2
Department of Physics, Faculty of Science, Helwan University, Ain Helwan, Cairo 11792, Egypt
*
Author to whom correspondence should be addressed.
Photonics 2022, 9(4), 205; https://doi.org/10.3390/photonics9040205
Submission received: 18 January 2022 / Revised: 14 March 2022 / Accepted: 16 March 2022 / Published: 22 March 2022
(This article belongs to the Topic Optical and Optoelectronic Materials and Applications)

Abstract

:
A model for generating femtosecond laser pulses from a low-power mode-locked laser of moderate temporal half-width was proposed. This was achieved by injecting the pulse into a single-mode inverted-populated Er-doped fiber where self-focusing and absorption were avoided. To initiate spectral broadening, the pulse was phase-modulated by subjecting a part of the fiber to an electric field of suitable intensity and frequency generated into a circular capacitor. To stimulate temporal compression, the phase-modulated pulse was introduced into a combination of two prism sets located symmetrically with respect to the x-axis. After passing the pulse through the first prism set, its spectral components were spatially separated in the y-axis. The spectral phases were manipulated by redirecting the spectral components through a slab cross-section that was subjected to a spatially modulated DC electric field. After passing the slab, the pulse is directed into the second prism set, where the spectral components were spatially overlapped and propagated outside the compressor with the same slope and dimension as before entering the compressor. Constructive super positioning of the phase-manipulated spectral components gave maximum intensity only at a specified location.

1. Introduction

Since the development of mode-locked lasers, shorter pulse durations have been required for many applications in science, medicine, and engineering. Such applications are given in [1,2,3,4,5,6,7,8,9]. To direct a conventional laser pulse into a femtosecond region, it is mathematically necessary to broaden its spectral bandwidth. Because this is achievable by applying nonlinear phenomena to the pulse, different methods for its realization were presented in [10,11,12,13,14]. According to the dispersion relations of the refractive index, the frequency components comprising the spectral bandwidth of the pulse have different velocities. Each frequency will reach an observer at different times, leading either to compressing or expanding the pulse duration. Since most of the materials have a positive-group-delay dispersion (PGDD) that expands the pulse duration, a material of negative-group-delay dispersion (NGDD) is required to reduce its duration. References [15,16,17,18,19,20,21,22] have illustrated different setups for its practical realization. From the above, it is evident that mutual combinations between broadening and compression methods will theoretically yield a temporal compressed pulse. So, for instance, the authors in [23] used nonlinear phenomena to increase the spectral range of a 90 fs laser pulse in a fiber and compressed the broadened pulse by a grating compressor to 30 fs. A pulse compression ratio of 12 (from 5.4 ps to 450 fs) was reported in [24] by using self-phase modulation in a 30 m optical fiber to broaden the spectrum and a grism compressor. The result of a theoretical analysis based on the solution of the nonlinear Schrödinger equation concerning pulse compression ratio and compression pulse quality as a function of fiber length and input pulse intensity was given in [25]. The study was carried out on an optical pulse chirped in a single-mode fiber by self-phase modulation and compressed with a grating compressor. Self-phase modulation in bulk material for spectral broadening the laser pulse and chirped mirror combined with two prisms for its compression was given in [26].
From the previous cited references, it is evident that most of the researchers used:
(1)
Powerful laser pulses for internally driving nonlinear phenomena to initiate the spectral broadening of the laser pulse;
(2)
Compressors with moveable compressor elements for creating compressed pulses at a desired location.

Aim of the Work

To generate a femtosecond laser pulse from a mode-locked laser pulse of low power and moderate halfwidth, it is necessary to externally stimulate nonlinear phenomena. The present article developed a model for their realization. Section 2.1 is devoted for finding a mathematical expression for an electric field that externally phase-modulates the laser pulse during its propagation in the fiber and a setup for its experimental realization. Section 2.2 is devoted to finding a method for external variation of the spectral optical path lengths. Therefore, the spectral phases are controlled such that the pulse compression occurs at a desired location without movement any of the compressor elements. A suggested setup for the experimental realization is also included.

2. Theory

Erbium-doped fibers (EDF) were chosen for generating femtosecond laser pulses for several reasons. They show high gain, polarization-independent amplification, low intrinsic optical noise [27], stable long-term performance, and possibility of being unidirectionally and bidirectionally pumped [28]. To avoid absorption and self-focusing of the laser pulse, the selected fiber was a single-mode and inverted populated. To achieve the cited goal, the pulses had to be externally spectrally broadened during propagation in the EDF. Two prism sets and a slab (for external manipulation of the spectral phases) were responsible for temporally compressing the spectrally broadened laser pulse.

2.1. Broadening of the Spectral Bandwidth of the Laser Pulse Using an Externally Triggered Phase Modulation

A laser pulse can be coherently spectrally broadened when it is subjected to nonlinear effects either internally through self-phase modulation or externally by amplitude, frequency, or phase modulation. For a small intensity, where self-phase modulation would be ineffective, external induction of wide-range spectral broadening would be required. This was achieved by an external high-frequency electric field that modulated the refractive index of the fiber according to the following equation [29]:
n ( Ω ) =   n 0 ( Ω ) +   n 2 E 2 ( t ) ,
where Ω is the angular frequency of the electric field of the laser pulse,   n 0 ( Ω ) and n 2 are the linear and nonlinear refractive indices of the fiber, respectively, and E(t) is the time-dependent external applied electric field. To create a femtosecond laser pulse, the induced spectral bandwidth (Δν), according to the Heisenberg uncertainty principle, must be in the order of 1/Δt, where Δt is the desired laser-pulse duration. To fulfill this criterion, the nonlinear part of Equation (1), n 2 E 2 ( t ) , has to be in the order of the linear part n 0 ( Ω ) . This can be achieved by subjecting the fiber to a high frequency and high electric-field strength. This electric field was created by solving the Maxwell’s equations within a circular-shaped capacitor in which a part of the fiber was wrapped. The required value of the electric-field strength was obtained for certain conditions concerning the radius of the capacitor and the frequency of the generator. Due to the considered symmetry, the Maxwell’s equations have been written in cylindrical coordinates as: (detailed solution see Appendix A)
Ε ( r , t ) r =   μ H ( r , t ) t ,
H ( r , t ) r + 1 r H ( r , t ) =   ε Ε ( r , t ) t .
With
H ( r , t ) =   Re [ H ( r ) exp ( i ω t ) ] ,
Ε ( r , t ) =   Re [ E ( r ) exp ( i ω t ) ] ,
it follows after setting Equations (4) and (5) into Equations (2) and (3) and applying some mathematical manipulations,
2 H ( r ) r 2 + 1 r H ( r ) r 1 r 2 H ( r ) +   ε μ ω 2 H ( r ) = 0 ,
1 μ i ω 2 E ( r ) r 2 + 1 r 1 μ i ω E ( r ) r   ε i ω E ( r ) = 0 ,
where Ε ( r , t ) and H ( r , t ) are the spatial and temporal electric and magnetic field strengths, respectively, E(r) and H(r) are the spatial components of the corresponding fields, respectively, and ε and μ are the electric permittivity and magnetic permeability, respectively.
With ε μ ω 2 =   k 2 , the solution of Equations (6) and (7) are given according to [30], by
E ( r ) =   AJ 0 ( kr ) + A Y 0 ( kr ) ,
H ( r ) =   BJ 1 ( kr ) + B Y 1 ( kr ) ,  
where J 0 and J 1 are the first kind Bessel functions of zero and first order, respectively,   Y 0   and   Y 1 are the second kind Bessel functions of zero and first order, respectively, and A and A   as   well   as   B ,   and   B are integration constants.
To avoid the singularity of Y 0   and   Y 1 at r   = 0 ,   A   and   B were set as equal to zero. To acquire B, Ampere’s law was applied:
H ( r , t ) rd φ   =   Bcos ( ω t ) J 1 ( ka ) ( 2 π a ) = t E ( r , t ) ( 2 π r ) dr   =   I ( t ) ,
where rd φ is an element of the circular path along H ( r ) with radius r, a is the radius of the capacitor, and I ( t ) is the feeding current of the capacitor with amplitude I o .
Accordingly, B and H ( r ) are determined, respectively, as:
B   =   I o / [ J 1 ( ka ) 2 π a ] ,
H ( r ) = [ I o / ( 2 π a ) ] [ J 1 ( kr ) / J 1 ( ka ) ] .
From Equation (3) and considering Equation (12), E ( r ) is calculated as follows:
E ( r )   =   I o kJ 0 ( kr ) / [ ( i ω ε )   ( 2 π a ) J 1 ( ka ) ] ,
With A =   I o k / [ ( i ω ε ) ( 2 π a ) J 1 ( ka ) ] , the relation between A and B is as follows:
B   =   iA ε / μ ,
and the real parts of E ( r , t ) and H ( r , t ) can be written, respectively as
Re [ E ( r , t ) ] =   Re [ AJ 0 ( kr ) exp ( i ω t ) ] =   I o kJ 0 ( kr ) sin ( ω t ) / [ ( ω ε   ( 2 π a ) J 1 ( ka ) ) ] ,
Re [ H ( r , t ) ] =   Re [ iA ε / μ J 1 ( kr ) exp ( i ω t ) ] =   I o k ε / μ J 1 ( kr ) cos ω t / [ ω ε   ( 2 π a ) J 1 ( ka ) ] .
From Equations (15) and (16), the time average of the Poynting vector within the capacitor was equal to zero. This result was in accordance with the assumption of the neglected polarization, absorption losses within the fiber, and the radiation into the surrounding. The latter can be achieved by setting the distance between the capacitor plates to be in the order of 1/100 of the wavelength of the modulating electric field. From Equation (15), it was evident that the electric field could attain large values as J 1 ( ka ) approached small values. Since this was possible for a discrete infinite number of ( ka ) values, we had to choose between those provided easily available frequencies and suitable capacitor dimensions. From the obtained results, the electric-field strength of the laser pulse was phase-modulated according to the following relation:
E laser ( z , t ) =   c 1 exp [ c 2 ( ( t t s ) / Δ t s ) 2 ] × exp ( i Ω 0 t ) exp [ ( 2 π i λ 0 )   ( n 0 ( Ω ) Z + n 2 E 2 ( t ) Z e ) ] ,
where exp [ c 2 ( ( t t s ) / Δ t s ) 2 ] is the time-dependent amplitude of the electric-field strength of the laser pulse to be phase-modulated; ts and Δ t s are the time of maximum amplitude and temporal width of the electric-field strength of the laser pulse, respectively, Ω 0 and λ 0 are the central vacuum angular frequency and the wavelength of the laser pulse, respectively, E 2 ( t ) is the square of the phase-modulating electric field obtained in Equation (15), which is multiplied by n 2 to be in the order of n 0 ( Ω ) , n 0 and n 2 are the linear and nonlinear refractive indices, respectively, Z is the total length of the fiber; and Z e is the length of the fiber subjected to the alternating electric field in the capacitor. To avoid sparking within the capacitor, it had to be placed in a grounded container that had been highly evacuated. This container consists of getter material to adsorb gases that were adsorbed on the surface of the material and might be free during the operation. To maintain a constant fiber temperature, the housing has to be subjected to an air-cooled Peltier cooler.

2.1.1. Computations

By keeping the denominator of Equation (15) constant, the fast Fourier transform (FFT) algorithm was applied on Equation (17) at Z e = 2 m for an amplified laser pulse of a temporal maximum amplitude at 1.8 W, duration Δ t s = 1 × 10 7 s, t s = 4 × 10 7 s, Ω 0 = 1.2 × 10 15 1/s,   λ 0 = 1.531   μ m , phase-modulated with electric-field strength 6.8 × 10 7 V/m, frequencies of 100 MHz, 1.2 GHz, n 2 =   3.2 × 10 16   m 2 / V 2 [31], and n 0 ( λ ) given by Sellmeier’s equation [32].
n 0 2 ( λ ) = [ 1 + B 1 λ 2 / ( λ 2 C 1 ) + B 2 λ 2 / ( λ 2 C 2 ) + B 3 λ 2 / ( λ 2 C 3 ) ] with B 1 = 0.6968 , B 2 = 0.40817 , B 3 = 0.8949 , C 1 = 6.853 × 10 2   μ m 2 ,   C 2 = 0.1161 × 10 2   μ m 2 and C 3 = 9.914   μ m 2 .

2.1.2. Results and Discussion

For computing the Fourier spectrum, the limited computer capacity made it impossible to account for more than 8 × 10 6 samples within the frequency range of 7 × 10 14     1.7 × 10 15   1/s for each of the modulating frequencies.
Figure 1 illustrates the modulus of the amplitude spectrum of the electric field of a laser pulse with a maximum power of 1.8 W, duration of 1 × 10 7 s, and subjected along two meters to a modulating frequency of 1.2 GHz. The computation was performed at Z e   = 2 m using (FFT) algorithm.
Figure 2 shows that according to Equation (17), the normalized temporal distribution of the amplitude of the electric field of the phase-modulated laser pulse had a maximum power of 1.8 W, a duration of 1 × 10 7 s, and a modulating frequency of 1.2 GHz, when computed at Z = 10 m.

2.2. Pulse Compression

To compress the spectral broadened laser pulse, it was directed into a combination of four prisms divided into two identical parallel sets (P1-P2 and P1′-P2′), as shown in Figure 3. The corresponding prisms of the two sets were located at the same x-positions and displaced in the y-axis. The first set acted as a Fourier transformer and was responsible for separating the frequency components comprising the pulse into parallel beams located at different y-positions. The second prism set acted as an inverse Fourier transformer and was responsible for the spatial overlap of the spectral components throughout the environment outside the compressor.
The frequency components comprising the pulse that was emitted from the first prism set were directed via 45° mirror M1 onto the cross section of a slab. The slab, which had a low absorption coefficient as well as a well-defined frequency-dependent Kerr coefficient, was sandwiched between two electrodes. To change the spectral optical paths, each location across the width of the slab was subjected to a DC electric field with a constant value and a variable length in y-axis. This was realized, as shown in Figure 4, by connecting one of the two electrodes to the ground and the other to a DC electric source. The latter electrode was made from a photo-conducting material and was homogeneously illuminated with an array of LEDs. To change the phases of the spectral components, a stretchable computer-controlled opaque screen was used to control the illuminated area across the width of the slab. The emitted radiation from the slab was redirected by the mirror M2 into the lower set of prisms. The location of M2 was arranged such that the spectral components intersected the third prism P2′ at the same locations with respect to the base from which they were emitted from the second prism P2 of the upper set. Due to the symmetry of the configuration, the spectral components were overlapped at the output edge of the fourth prism and had the same slope. The slab’s length was chosen so that the compressor would be approximately quadratic-shaped and the electric field across its thickness would have a moderate value. For all spectral components to be constructively superimposed at a desired location, the spectral phases had to be equal or shifted by multiple of 2 π at that location. As shown in Figure 3, to fulfill this condition, the spectral optical path lengths across the setup and outside it had to be determined without subjecting the slab to an electric field. By dividing each of the obtained spectral optical path lengths by its corresponding vacuum wavelength ( λ ) and ignoring the integer number, a dimensionless fraction of the spectral wavelength was obtained. The corresponding phase of this fraction was obtained as follows:
φ ( Ω ) = 2 π × Fraction ( λ ) ,
and the spectral angle Δ φ ( Ω ) for completing φ ( Ω ) to 2 π was:
Δ φ ( Ω ) = 2 π φ ( Ω ) .
To account for Δ φ ( Ω ) , as found in Equation (19), the nonlinear part of the refractive index of the slab’s material, given in Equation (1), was activated. This was realized by subjecting a particular length l ( Ω ) along the slab at the location corresponding to the frequency across the slab’s width to a DC electric field, according to:
Δ φ ( Ω ) = 2 π λ n 2 E 2 l ( Ω )
Solving Equation (20) for l ( Ω ) , one obtains
l ( Ω ) = Δ φ ( Ω ) λ 2 π   n 2 E 2 = Δ φ ( Ω )   c Ω n 2 E 2 ,
where n 2 is the nonlinear Kerr coefficient of the slab’s material that was chosen to be the same as that of the fiber’s material,   l ( Ω ) is the length along the slab corresponding to the frequency Ω that was subjected to the DC electric field, and c is the vacuum velocity of light.
To simulate the temporal distribution resulting from the compressor at a desired location, an inverse fast Fourier transform (IFFT) algorithm was applied on the manipulated spectral phases.

2.2.1. Computations

The computation for determining l ( Ω ) given in (21) was performed for equilateral prisms of side length l 1 = 0.22 m, height as h = 0.19 m, length between the end of the base of the first prism and the apex of the second one as l 2 = 0.05 m, slab length as 0.2 m, slab thickness as 0.003 m, and electric field influencing the slab as 100 V/mm.

2.2.2. Results and Discussion

Figure 5 shows the fitted curve of the spectral length l ( Ω ) subjected to the DC electric field that influenced the slab to stimulate the frequency components to be in-phase at the desired location. The computation was carried out for a bandwidth of 10 15 1/s induced by modulating the pulsed laser phase by a frequency 1.2 GHz.
Figure 6 represents the compressed laser pulses computed for the modulating frequencies of 100 MHz and 1.2 GHz. For each of the modulating frequencies, the computation was performed for a fixed number of frequency samples 8 × 10 6 . The results showed that the compression for the modulation frequency at 100 MHz was better than that found at 1.2 GHz. This behavior was due to the larger ratio of the selected frequencies relative to the actual frequencies at 100 MHz compared to that at 1.2 GHz. This may be attributed to the fact that 1.2 GHz were more under-sampled or strongly aliased. In practical application, all frequencies would be in-phase if the slab is subjected to the desired distribution of the electric field given by l ( Ω ) . Therefore, in this case, it was suggested that much better compression would be achieved.

3. Conclusions

Although this proposal allowed an ordinary femtosecond laser to be generated, it differed from previous studies in the following respects. The pulses to be compressed were obtained from a mode-locked laser of moderate halfwidth and low power, such that self-phase-modulation could not be initiated. The compression was achieved using a one-stage compressor. The compressed pulses were generated at any desired location without the movement of the compressor elements. Moreover, from the above, it is evident that the spectral broadening of the laser pulse as well as the spectral phase manipulation depend only on n 2 E 2 that must be in the order of n 0 . This condition can be easily achieved by the external variation of the electric field. Consequently, this model can be used for any frequency that is not markedly attenuated during propagation in the proposed setup.

Author Contributions

Conceptualization, M.H.A.-e., S.E.I.H.E. and A.M.A.-R.; methodology, M.H.A.-e. and S.E.I.H.E.; software, M.H.A.-e. and A.M.A.-R.; validation, M.H.A.-e., S.E.I.H.E., A.F.H. and A.M.A.-R.; writing—original draft preparation, M.H.A.-e. and A.M.A.-R.; writing—review and editing, M.H.A.-e. and A.F.H.; visualization, supervision, S.E.I.H.E. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to thank Jala El-Azab, head of the Department of Engineering Applications of Lasers at the National Institute of Laser Enhanced Sciences, Cairo University, for her fruitful discussion.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

To uncover the desired electric-field strength distribution necessary for spectral broadening the laser pulse, Maxwell’s equations were solved:
× H   = D t + J ,
Since the field is in the insulator, J will be set equal to zero to get with D = Ε e
× H   = D t =   ε E t
× E   = B   t
with B = µ H
×   E   = µ B   t
where H = magnetic field strength
Considering E to be parallel to the normal of the area of the capacitor plates, and its radial component to be omitted, H will build concentric circles about the center of the capacitor. Due to the considered geometry, the Maxwell equations have to be rewritten in cylindrical coordinates as:
E ( r , t ) r =   µ H ( r , t ) t
H ( r , t ) r + 1 r H ( r , tr , t ) =   ε E ( r , t ) t
considering the fields to be sinusoidal in time such that H ( r , t ) =   Re [ H ( r ) e i ω t ] , E ( r , t ) =   Re [ E ( r ) e i ω t ] .
The time derivative of the fields gives
H ( r , t ) t =   i ω H ( r ) e i ω t
E ( r , t ) t =   i ω E ( r ) e i ω t
Substituting the obtained expression in Equations (A6) and (A7) in Equations (A4) and (A5), one obtains
E ( r ) r =   i μ ω H ( r )
H ( r ) r + 1 r H ( r ) =   i ε ω E ( r )
Since the above equations are coupled, to decouple them one has to differentiate Equation (A9) w.r.t. (r), to obtain
2 H ( r ) r 2 + 1 r H ( r ) r 1 r 2 H ( r ) =   i ε ω E ( r ) r
substituting for E ( r ) r in Equation(A10) from Equation (A8) one obtains,
2 H ( r ) r 2 + 1 r H ( r ) r 1 r 2 H ( r ) + ε μ ω 2 H ( r ) = 0
Differentiating Equation (A8) w.r.t. (r) one obtains
2 E ( r ) r 2 =   i μ ω H ( r ) r
Substituting for H ( r ) r and H ( r ) in Equation (A9) from Equations (A12) and (A8), respectively, gives,
1 i μ ω 2 E ( r ) r 2 + 1 r 1 i μ ω E ( r ) r i ε ω E ( r ) = 0
with k 2 = ω 2 ε μ , Equations (A11) and (A13) after normalization can be, respectively, written as
2 H ( r ) r 2 + 1 r H ( r ) r + ( k 2 1 r 2 ) H ( r ) = 0
2 E ( r ) r 2 + 1 r E ( r ) r + k 2 E ( r ) = 0
With kr   =   x , k r   = x , so k   = x r and Equations (A14) and (A15) can be written as:
2 E ( r ) x 2 + 1 x E ( r ) x + E ( r ) = 0
2 H ( r ) x 2 + 1 x H ( r ) x + ( 1 1 x 2 ) H ( r ) = 0
The solution of Equations (A16) and (A17) are according to [30] given by:
E ( r ) =   AJ 0 ( kr ) +   A   Y 0 ( kr )
H ( r ) =   BJ 1 ( kr ) +   B   Y 1 ( kr )
Since Y 0 , Y 1 are second-kind Bessel functions of zero and first order, respectively, and have singularity at kr = 0, thus   A   and   B must vanish, giving:
E ( r ) =   AJ 0 ( kr ) H ( r ) =   BJ 1 ( kr )
since J 1 ( 0 ) = 0 ,   then   H ( 0 ) = 0 .
To determine B, the Ampere’s law has to be applied on the magnetic field distribution H ( r ) .
From the relation
H ( r   =   a ) .2 π a   = D t dS   =   I   with   D   =   ε 0 ε r . E
where I is complex amplitude of the current through the wire feeding the capacitor.
Setting r   =   a in the right-hand side of Equation (A19) one obtains, after substituting in the left-hand side for H ( r   =   a ) from the above relation
B   = I 2 π a 1 J 1 ( ka ) ,
and thus H ( r ) can be written as
H ( r ) = I 0 2 π a J 1 ( kr ) J 1 ( ka )
To determine E ( r ) , the obtained equation of H ( r ) and its derivative w.r.t (r) have to be substituted into Equation (A9). This gives:
E ( r ) = 1 i ω ε ( I 0 2 π a 1 J 1 ( ka ) ( k kr J 1 ( kr ) + kJ 0 ( kr ) ) + 1 r I 2 π a J 1 ( kr ) J 1 ( ka ) )
E ( r ) = I 0 k i ω ε ( 2 π a ) J 1 ( k a ) J 0 ( k r )
With ,   A   = I 0 k i ω ε ( 2 π a ) J 1 ( ka )
The following relation is found
B A = i ω ε k
Substituting for k from the relation k   =   ω ε μ it follows that
B A =   i 1 ε μ ε   =   i ε μ
B   =   i ε μ A
From the obtained spatial amplitude distributions of the electrical and magnetic fields, their spatial and time distributions are given by:
E ( r , t ) =   Re [ AJ 0 ( kr ) e i ω t ] ,
H ( r , t ) =   Re [ BJ 1 ( kr ) e i ω t ] ,
E ( r , t ) = | A | J 0 ( kr ) sin ω t ,
and   H ( r , t ) = ε μ | A | J 1 ( kr ) cos ω t   ,
with
| A | =   Ĩ k ω ε ( 2 π a ) J 1 ( ka )
Setting | A | =   A , it follows that
E ( r , t ) =   AJ 0 ( kr ) sin ω t
H ( r , t ) = ε μ AJ 1 ( kr ) cos ω t

References

  1. Dey, D.; Tiwari, A.K. Controlling Chemical Reactions with Laser Pulses. ACS Omega 2020, 5, 17857–17867. [Google Scholar] [CrossRef] [PubMed]
  2. Ikeda, K.; Kotaki, H.; Nakajima, K. X-ray generation by femtosecond laser pulses and its application to soft X-ray imaging microscope. AIP Conf. Proc. 2002, 634, 268–275. [Google Scholar] [CrossRef]
  3. Afshari, M.; Krumey, P.; Menn, D.; Nicoul, M.; Brinks, F.; Tarasevitch, A.; Sokolowski-Tinten, K. Time-resolved diffraction with an optimized short pulse laser plasma X-ray source. Struct. Dyn. 2020, 7, 014301. [Google Scholar] [CrossRef]
  4. Abbas, A.M.; Pan, Q.; Mandon, J.; Cristescu, S.M.; Harren Frans, J.M.; Khodabakhsh, A. Time-resolved mid-infrared dual-comb spectroscopy. Sci. Rep. 2019, 9, 17247. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Borgwardt, M.; Omelchenko, S.T.; Favaro, M.; Plate, P.; Hohn, C.; Abou-Ras, D.; Schwarzburg, K.; van de Krol, R.; Atwater, H.A.; Lewis, N.S.; et al. Femtosecond time-resolved two-photon photoemission studies of ultrafast carrier relaxation in Cu2O photoelectrodes. Nat. Commun. 2019, 10, 2106. [Google Scholar] [CrossRef]
  6. Subramaniam, T.K. Erbium Doped Fiber Lasers for Long Distance Communication Using Network of Fiber Optics. Am. J. Opt. Photonics 2015, 3, 34–37. [Google Scholar] [CrossRef]
  7. Ono, M.; Hata, M.; Tsunekawa, M.; Nozaki, K.; Sumikura, H.; Chiba, H.; Notomi, M. Ultrafast and energy-efficient all-optical switching with graphene-loaded deep-subwavelength plasmonic waveguides. Nat. Photonics 2019, 14, 37–43. [Google Scholar] [CrossRef] [Green Version]
  8. Mottay, E.; Liu, X.; Zhang, H.; Mazur, E.; Sanatinia, R.; Pfleging, W. Industrial applications of ultrafast laser processing. MRS Bull. 2016, 41, 984–992. [Google Scholar] [CrossRef]
  9. Alio, J.L.; AAbdou, A.A.; Puente, A.A.; Zato, M.A.; Nagy, Z. Femtosecond Laser Cataract Surgery: Updates on Technologies and Outcomes. J. Refract. Surg. 2014, 30, 420–427. [Google Scholar] [CrossRef] [PubMed]
  10. Karasawa, N.; Morita, R.; Shigekawa, H.; Yamashita, M. Generation of intense ultrabroadband optical pulses by induced phase modulation in an argon-filled single-mode hollow waveguide. Opt. Lett. 2000, 25, 183–185. [Google Scholar] [CrossRef] [Green Version]
  11. Hanna, M.; Delen, X.; Lavenu, L.; Guichard, F.; Zaouter, Y.; Druon, F.; Georges, P. Nonlinear temporal compression in multipass cells: Theory. J. Opt. Soc. Am. 2017, 34, 1340–1347. [Google Scholar] [CrossRef]
  12. Weidner, P.; Penzkofer, A. Spectral broadening of picosecond laser pulses in optical fibres. Opt. Quantum Electron. 1993, 25, 1–25. [Google Scholar] [CrossRef] [Green Version]
  13. Lu, C.-H.; Witting, T.; Husakou, A.; Vrakking, M.; Kung, A.H.; Furch, F.J. Sub-4 fs laser pulses at high average power and high repetition rate from an all-solid-state setup. Opt. Express 2018, 26, 8941–8956. [Google Scholar] [CrossRef] [PubMed]
  14. Saleh, M.F.; Biancalana, F. Soliton dynamics in gas-filled hollow-core photonic crystal fibers. J. Opt. 2016, 18, 013002. [Google Scholar] [CrossRef] [Green Version]
  15. Schmidt, B.E.; Béjot, P.; Giguère, M.; Shiner, A.; Trallero, C.; Bisson, E.; Kasparian, J.; Wolf, J.-P.; Villeneuve, D.; Kieffer, J.-C.; et al. Compression of 1.8 μm laser pulses to sub two optical cycles with bulk material. Appl. Phys. Lett. 2010, 96, 121109. [Google Scholar] [CrossRef] [Green Version]
  16. Vuong, L.T.; Lopez-Martens, R.B.; Hauri, C.P.; Ruchon, T.; L’Huillier, A.; Foster, M.A.; Gaeta, A.L. Optimal pulse compression via sequential filamentation. In Proceedings of the Quantum Electronics and Laser Science Conference, Baltimore, ML, USA, 6–11 May 2007; p. JWE1. [Google Scholar]
  17. Liao, K.-H.; Chen, M.-Y.; Flecher, E.; Smirnov, V.I.; Glebov, L.B.; Galvanauskas, A. Large-aperture chirped volume Bragg grating based fiber CPA system. Opt. Express 2007, 15, 4876–4882. [Google Scholar] [CrossRef] [PubMed]
  18. Chauhan, V.; Bowlan, P.; Cohen, J.; Trebino, R. Single-diffraction-grating and grism pulse compressors. Opt. Soc. Am. 2010, 27, 619–624. [Google Scholar] [CrossRef]
  19. Akturk, S.; Gu, X.; Kimmel, M.; Trebino, R. Extremely simple single-prism ultrashort-pulse compressor. Opt. Express 2006, 14, 10101–10108. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Dombi, P.; Yakovlev, V.S.; O’Kee, K.; Fuji, T.; Lezius, M.; Tempea, G. Pulse compression with time-domain optimized chirped mirrors. Opt. Express 2005, 13, 10888–10894. [Google Scholar] [CrossRef] [Green Version]
  21. Schulte, J.; Sartorius, T.; Weitenberg, J.; Vernaleken, A.; Russbueldt, P. Nonlinear pulse compression in a multi-pass cell. Opt. Lett. 2016, 41, 4511–4514. [Google Scholar] [CrossRef] [PubMed]
  22. Nada, Y.S.; El-azab, J.M.; Othman, H.A.; Mohamad, T.; Maize, S.M.A. Interaction between Self Phase Modulation and Positive Group Velocity Dispersion in PMMA Polymer for Simplified Thin Film Compressor of High-Intensity Ultrashort Laser Pulses. Nonlinear Opt. Quantum Opt. 2020, 52, 299–311. [Google Scholar]
  23. Shank, C.V.; Fork, R.L.; Yen, R.; Stolen, R.H. Compression of femtosecond optical pulses. Appl. Phys. Lett. 1982, 40, 761–763. [Google Scholar] [CrossRef]
  24. Nikolaus, B.; Grischkowsky, D. 12× pulse compression using optical fibers. Appl. Phys. Lett. 1982, 42, 117–121. [Google Scholar] [CrossRef]
  25. Tomlinson, W.J.; Stolen, R.H.; Shank, C.V. Compression of optical pulses chirped by self-phase modulation in fibers. J. Opt. Soc. Am. B 1984, 1, 139–149. [Google Scholar] [CrossRef]
  26. Mével, E.; Tcherbakoff, O.; Salin, F.; Constant, E. Extracavity compression technique for high-energy femtosecond pulses. J. Opt. Soc. Am. B 2003, 20, 105–108. [Google Scholar] [CrossRef]
  27. Jiménez, M.A.B.; Kuzin, E.A.; Ibarra-Escamilla, B.; Flores-Rosas, A. Optimization of the two-stage single-pump erbium-doped fiber amplifier with high amplification for low frequency nanoscale pulses. Opt. Eng. 2007, 46, 125007. [Google Scholar] [CrossRef]
  28. Quintela, M.A.; Quintela, C.; Lomer, M.; Madruga, F.J.; Conde, O.M.; Lopez-Higuera, J.M. Comparison between a symmetric bidirectional-pumping and a unidrectional-pumping configurations in an erbium fiber ring laser. Int. Soc. Opt. Eng. 2007, 6619, 440–443. [Google Scholar]
  29. Milonni, P.W.; Eberly, J.H. Introduction to Nonlinear Optics, 2nd ed.; John Wiley Sons, Inc.: Hoboken, NJ, USA, 2010; Chapter 10. [Google Scholar]
  30. Bronshtein, I.N.; Semendyayev, K.A.; Musiol, G.; Muehlig, H. Handbook of Mathematics, 5th ed.; Springer: Berlin/Heidelberg, Germany, 2007. [Google Scholar]
  31. Qian, S.; Li, J.; Zhang, Y. Absorptive Nonlinearity in Er-Doped Optical Fiber. Wuhan Univ. J. Nat. Sci. 1999, 4, 175–178. [Google Scholar] [CrossRef]
  32. Fleming, J.W.; Wood, D.L. Refractive index dispersion and related properties in fluorine doped silica. Appl. Opt. 1983, 22, 3102–3104. [Google Scholar] [CrossRef]
Figure 1. The modulus given in Au of the amplitude spectrum resulting from the (FFT) algorithm computed at Ze = 2 m for the electric-field amplitude of a laser pulse of input maximum power 1.8 W, duration of 1 × 10 7 s, and a modulating frequency of 1.2 GHz.
Figure 1. The modulus given in Au of the amplitude spectrum resulting from the (FFT) algorithm computed at Ze = 2 m for the electric-field amplitude of a laser pulse of input maximum power 1.8 W, duration of 1 × 10 7 s, and a modulating frequency of 1.2 GHz.
Photonics 09 00205 g001
Figure 2. The normalized amplitude of the phase−modulated laser pulse of maximum power 1.8 W, duration 1 × 10 7 s computed at the end of the fiber z = 10 m and a modulating frequency 1.2 GHz.
Figure 2. The normalized amplitude of the phase−modulated laser pulse of maximum power 1.8 W, duration 1 × 10 7 s computed at the end of the fiber z = 10 m and a modulating frequency 1.2 GHz.
Photonics 09 00205 g002
Figure 3. Overview of the compressor setup P1, P2, P1′, and P2′ are four prisms, L is a pulsed laser, P is a pump source, C is a capacitor M1 and M2 are mirrors and S is a slab subjected to a DC electric field.
Figure 3. Overview of the compressor setup P1, P2, P1′, and P2′ are four prisms, L is a pulsed laser, P is a pump source, C is a capacitor M1 and M2 are mirrors and S is a slab subjected to a DC electric field.
Photonics 09 00205 g003
Figure 4. Overview of cross-sections in the slab responsible for stimulating the spectral phases to be in phase at any desired location.
Figure 4. Overview of cross-sections in the slab responsible for stimulating the spectral phases to be in phase at any desired location.
Photonics 09 00205 g004
Figure 5. The fitted curve of the slab’s length l ( Ω ) subjected to a DC electric-field distribution across the slab’s width to equalize the phases of the spectral components of the broadened laser pulse at any desired location considering a laser with a pulse duration of 10 7   s, maximum power of 1.8 W, and phase modulated with a frequency of 1.2 GHz.
Figure 5. The fitted curve of the slab’s length l ( Ω ) subjected to a DC electric-field distribution across the slab’s width to equalize the phases of the spectral components of the broadened laser pulse at any desired location considering a laser with a pulse duration of 10 7   s, maximum power of 1.8 W, and phase modulated with a frequency of 1.2 GHz.
Photonics 09 00205 g005
Figure 6. The compressed laser pulses for modulating frequencies (a) 100 MHz and (b) 1.2 GHz.
Figure 6. The compressed laser pulses for modulating frequencies (a) 100 MHz and (b) 1.2 GHz.
Photonics 09 00205 g006aPhotonics 09 00205 g006b
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Abo-elenein, M.H.; Hassab Elnaby, S.E.I.; Hassan, A.F.; Abd-Rabou, A.M. Modeling for Generating Femtosecond Pulses in an Er-Doped Fiber Using Externally Controlled Spectral Broadening and Compression Mechanisms. Photonics 2022, 9, 205. https://doi.org/10.3390/photonics9040205

AMA Style

Abo-elenein MH, Hassab Elnaby SEI, Hassan AF, Abd-Rabou AM. Modeling for Generating Femtosecond Pulses in an Er-Doped Fiber Using Externally Controlled Spectral Broadening and Compression Mechanisms. Photonics. 2022; 9(4):205. https://doi.org/10.3390/photonics9040205

Chicago/Turabian Style

Abo-elenein, Mohamed Hemdan, Salah Eldeen Ibrahim Hassab Elnaby, Amin Fahim Hassan, and Afaf Mahmoud Abd-Rabou. 2022. "Modeling for Generating Femtosecond Pulses in an Er-Doped Fiber Using Externally Controlled Spectral Broadening and Compression Mechanisms" Photonics 9, no. 4: 205. https://doi.org/10.3390/photonics9040205

APA Style

Abo-elenein, M. H., Hassab Elnaby, S. E. I., Hassan, A. F., & Abd-Rabou, A. M. (2022). Modeling for Generating Femtosecond Pulses in an Er-Doped Fiber Using Externally Controlled Spectral Broadening and Compression Mechanisms. Photonics, 9(4), 205. https://doi.org/10.3390/photonics9040205

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop