Next Article in Journal
Real-Time Musical Conducting Gesture Recognition Based on a Dynamic Time Warping Classifier Using a Single-Depth Camera
Next Article in Special Issue
Step-Pulse Modulation of Gain-Switched Semiconductor Pulsed Laser
Previous Article in Journal
Effect of Height Ratio and Mass Ratio on Structure—Soil—Structure Interaction of Two Structures Using Centrifugal Experiment
Previous Article in Special Issue
Spin-ARPES EUV Beamline for Ultrafast Materials Research and Development
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analytical Solutions for the Propagation of UltraShort and UltraSharp Pulses in Dispersive Media

Department of Electrical and Electronic Engineering, Ariel University, Ariel 40700, Israel
Appl. Sci. 2019, 9(3), 527; https://doi.org/10.3390/app9030527
Submission received: 13 January 2019 / Revised: 28 January 2019 / Accepted: 29 January 2019 / Published: 4 February 2019
(This article belongs to the Special Issue Ultrafast Laser Pulses)

Abstract

:
Ultrashort pulses are severely distorted even by low dispersive media. While the mathematical analysis of dispersion is well known, the technical literature focuses on pulses, Gaussian and Airy pulses, which keep their shape. However, the cases where the shape of the pulse is unaffected by dispersion is the exception rather than the norm. It is the objective of this paper to present a variety of pulse profiles, which have analytical expressions but can simulate real-physical pulses with great accuracy. In particular, the dynamics of smooth rectangular pulses, physical Nyquist-Sinc pulses, and slowly rising but sharply decaying ones (and vice versa) are presented. Besides the usage of this paper as a handbook of analytical expressions for pulse propagations in a dispersive medium, there are several new findings. The main findings are the analytical expressions for the propagation of chirped rectangular pulses, which converge to extremely short pulses; an analytical approximation for the propagation of super-Gaussian pulses; the propagation of the Nyquist-Sinc Pulse with smooth spectral boundaries; and an analytical expression for a physical realization of an attenuation compensating Airy pulse.

1. Introduction

Ultrashort pulses are ubiquitous in numerous applications including in bio- and medical-imaging (see, for example, [1,2,3]), optical communications (see, for example, [4,5]), and microscopy (see, for example, Ref. [6,7]), etc.
One of the clear shortcomings of short pulses is the consequences of their wide spectrum. When the spectrum is spectrally wide, the absorption variations within the spectrum increases. However, even if the absorption is approximately homogenous within the pulse spectral range, the refraction index is wavelength-dependent and therefore the pulse experiences dispersion effects.
Dispersion occurs whenever each one of the pulse’s spectral components propagates with a different velocity. Consequently, the pulse experiences shape distortions. The temporally shorter the pulse, the larger the distortions. Moreover, while the dispersion effects are proportional to the medium’s length, it is inversely proportional to the square of the pulse’s temporal width; therefore, ultra-short pulses are prone to being affected by dispersion even in relatively short media [8,9,10].
Clearly, most of the dispersion mitigation techniques were implemented in the optical communication arena, where dispersion can have a destructive effect on the transmitting signals. It has been shown that via ordinary smf28 fiber (the most ubiquitous fiber type) it is impossible to decode data, which is transmitted at a 10 Gb/s rate, beyond a distance of 160 km (see Ref. [11] and in the quantum analogy see Ref. [12]). Clearly, in the presence of noise, this distance shrinks substantially. If the data is carried by picosecond pulses (which is only 1% of the temporal width of the 10 Gb/s signal’s individual pulses), this maximum distance is reduced to merely 16 m, and the distortion effects will be felt within less than a meter. In the femtosecond domain, the dispersion effects are, using the terminology and example of Ref. [13], more dramatic: A 10 fs laser pulse in the 800 nm spectral regime, which propagates through 4 mm of BK7 glass, will be temporally broadened to 50 fs.
It is well known that dispersion is governed by the dispersion equation, which can easily be solved numerically. However, to evaluate the effects of dispersion and to predict them, some intuition is required. In most textbooks, Gaussian pulses are used as a litmus tool to predict the dispersion effect [8,13,14]. Gaussian pulses are very useful since, on the one hand, they are localized pulses with a well-defined width and energy, and on the other hand, the shape of their spectrum is mathematically similar to the dispersion effect. Therefore, the propagation of Gaussian pulses in dispersion medium has a relatively simple analytical solution. This simple pulse model can be used to qualitatively evaluate some fundamental limitations. However, it fails in quantifying their values (see Ref. [11] and the references therein). In fact, the Gaussian model predicts that the pulse keeps its Gaussian shape, and its width is proportional to the square root of the elapsed period. This conduct resembles the spatial spreading of a Gaussian distribution of particles in a diffusion process [15]. Since the “particle” realization of the dispersion process is common in the popular description of chromatic dispersion, the similarity to diffusion may be misleading, since, in fact, this similarity appears only in the Gaussian scenario. Whenever the initial pulse profile deviates from the Gaussian’s one, the differences between diffusion and dispersion are clearly shown. In general, the two processes tend to smooth out the signal. However, the processes are fundamentally different, and the main differences appear after short distances. This discrepancy clearly appears in most optical communication protocols, where the pulses never have Gaussian profiles. Consequently, the signal distortion is completely different from the Gaussian analysis predictions. To overcome this problem, super-Gaussian pulses are commonly used to model real communications scenarios [16,17,18]. Despite the fact that super-Gaussian pulses are much better approximation to real pulses, they have two weaknesses: they cannot be implemented in Non-Return-to-Zero (NRZ) protocols and they do not have a close-form analytical solution.
It is the objective of this paper to assemble important pulses profile, whose propagation in dispersive medium can be expressed in analytical terms. It is shown that there are profiles, which are much better candidates to simulate real pulses, but their dynamics can be formulated in a closed-form analytical expressions.
Moreover, it is shown that the dynamics of almost any practical pulse profile in dispersive media can be simulated and approximated by analytical expressions.
In the first part of the paper, after the introduction of the general theory (Section 2) and the presentation of fundamental theorems (Section 3), the propagation of a Gaussian pulse is presented (Section 4). In the second part, it is shown that it is possible to present, analytically, the propagation of singular pulses (Section 5), which can be written initially with exponentials and step functions (like rectangular pulses and exponential-step pulses). In the third part, it is shown that every step function can be replaced with a smooth step function, and therefore any singular pulse, which is presented in the previous part of the paper, can be replaced with its counterpart smooth version (Section 6). In particular, smooth rectangular pulses can replace super-Gaussian ones. The fourth part is dedicated to the important pulse, which is singular in the spectral domain, i.e., the Nyquist-Sinc pulse (Section 7). First, the analytical expression for its propagation in dispersive medium is presented, and then the analytical expression of the propagation of its spectrally-smooth, i.e., physical, counterpart is presented as well. In the last part of the paper, accelerated-Airy pulses and attenuation compensating pulses are presented, as well as their physical counterparts (Section 8).

2. Generic Dispersion Analysis

Let A ( t , z ) represent the envelope of an electromagnetic (EM) pulse, which propagates in the z-direction in a dispersive medium with the dispersion coefficient β 2 , and thus A ( t , z ) obeys the following equation [14]:
i A ( t , z ) z = β 2 2 2 A ( t , z ) t 2
where t = t z / v is the time in a frame of reference that travels at the velocity of light in the medium ( v = c / n , where c is the vacuum light’s velocity, and n is the mean refractive index), i.e., t is the real time. This equation is valid provided all spectral components of the pulse are governed by the following dispersion relation:
n c ω + β 2 2 ω 2 = k
If A ( t , 0 ) is the initial pulse profile, then the pulse profile at the end of the medium A ( t , z ) can be evaluated by the following convolution:
A ( t , z ) = K ( t t , z ) A ( t , 0 ) d t
where the dispersion Kernel (see Appendix A) is
K ( t t , z ) = ( 2 π i β 2 z ) 1 / 2 exp ( i ( t t ) 2 2 β 2 z )

3. Fundamental Dispersion Theorems

To analyze the dispersion dynamics of the different pulses, we need some analytical tools, which will be formulated as the following theorems.

3.1. Pulse Boosting and Decaying

Let A ( t , 0 ) be an initial pulse whose dispersive dynamics is known, i.e., at the end of the dispersive medium its profile is also given as A ( t , z ) . Then, if the initial pulse is boosted by shifting the pulse’s spectrum, i.e., by multiplying the profile by a harmonic function:
A ( t , z = 0 ) = exp ( i ω 0 t ) A ( t , 0 )
then the pulse’s profile at the end of the medium (after a distance z ) is (see Appendix B)
A ( t , z > 0 ) = exp ( i β 2 2 ω 0 2 z + i ω 0 t ) A ( t + β 2 ω 0 z , z )
Clearly, this theorem can be generalized to the case where the initial pulse is multiplied by an exponent:
A ( t , z = 0 ) = exp ( a t ) A ( t , 0 )
where a can be any real, imaginary, or complex number. In which case, if the initial function does not diverge, the solution at the end of the medium is
A ( t , z > 0 ) = exp ( i β 2 2 a 2 z + a t ) A ( t i β 2 a z , z )

3.2. Pulse Chirping

Let A ( t , 0 ) be an initial pulse whose dispersive dynamics is known, i.e., at the end of the dispersive medium, its profile is A ( t , z ) ; then, if the initial pulse is chirped, i.e.,:
A ( t , 0 ) = exp ( i q t 2 ) A ( t , 0 )
for any real q, the pulse’s profile at the end of the medium (after a distance z) is (see Appendix C)
A ( t , z ) = ( 1 + 2 β 2 z q ) 1 / 2 exp ( i t 2 2 β 2 z + 1 / q ) A ( t 1 + 2 β 2 z q , z 1 + 2 β 2 z q )
Again, this theorem can be generalized to the case where the initial profile is multiplied by a generic Gaussian. That is, if the initial pulse is
A ( t , 0 ) = exp ( q t 2 ) A ( t , 0 )
where q can be any complex number, then as long as q > 0 , the solution at the end of the medium is
A ( t , z ) = ( 1 i 2 β 2 z q ) 1 / 2 exp ( i t 2 2 β 2 z + i / q ) A ( t 1 i 2 β 2 z q , z 1 i 2 β 2 z q )

4. Gaussian Pulse

The simplest pulse, which is ubiquitous in the literature, is the Gaussian one. This is one of the few cases where the pulse has the same shape as its Fourier Transform. Therefore, since the dispersion dynamics can be formulated as a convolution with a complex Gaussian (Equation (4)), then in general, the pulse profile remains a complex Gaussian. In this case, the initial pulse profile is
A ( t , 0 ) = A 0 exp ( ( t / θ ) 2 )
and thus the pulse’s profile at the end of the dispersive medium is (following Equation (12))
A ( t , z ) = A 0 θ θ 2 2 i β 2 z exp ( 1 θ 2 2 i β 2 z t 2 )
The profile of Equation (14) is plotted in Figure 1 and a spatiotemporal presentation of its intensity | A ( t , z ) | 2 is presented in Figure 2. Hereinafter, we adopt the following dimensionless parameters: τ t / θ and ζ β 2 z / θ 2 .
Sample Application: The pulses that many mode locked lasers emit can be simulated by Gaussian pulses (see Ref. [14]). When such a pulse with a Full Width at High Maximum (FWHM) of 1 ps (which corresponds to θ = F W H M / ln 4 0.85   ps ) and a carrier wavelength of λ = 1.06   μ m enters a BK7 glass, then | β 2 | 24   ps 2 / Km , and ζ = 1 in Figure 2 corresponds to z = ζ θ 2 / β 2 30   m .

4.1. Boosted Gaussian

When the initial pulse is a boosted Gaussian,
A ( t , z = 0 ) = A 0 exp ( t 2 / θ 2 + i ω 0 t )
then using Relation (6), the final (after a distance z) profile is
A ( t , z ) = A 0 θ θ 2 2 i β 2 z exp ( i β 2 2 ω 0 2 z + i ω 0 t ) exp ( ( t + β 2 ω 0 z ) 2 θ 2 2 i β 2 z )
which means that beside the pulse spreading, the pulse’s peak propagates with reciprocal “velocity” β 2 ω 0 . The larger the spectral shift ω 0 , the larger the temporal shift of the pulse. The pulse’s profile and its intensity’s dynamics are presented in Figure 3 and Figure 4, respectively.

4.2. Chirped Gaussian

When the initial pulse is chirped, i.e., the pulse’s phase has an additional quadratic dependence on time, i.e.,
A ( t , 0 ) = A 0 exp ( t 2 θ 2 i q t 2 )
then, using Equation (10) the final pulse is
A ( t , z ) = A 0 1 + 2 β 2 z ( q i / θ 2 ) exp ( t 2 θ 2 [ i q θ 2 + 1 1 + 2 β 2 z ( q i θ 2 ) ] )
In Figure 5, the temporal dynamics of the pulse (Equation (18)) is presented for the same parameters as Figure 2 but with additional q = 2 / θ 2 . The shrinkage of the pulse’s width is clearly shown, which is a consequence of the fact that the chirping has widened the spectral width of the initial pulse.

5. Singular Pulses

Before going into the propagation of profiles, which can simulate real pulses, it is instructive to investigate a family of pulses with singularity points. The singularity points are points in which either the pulse’s amplitude or its derivatives are discontinuous. Despite their singularities, these pulses can simulate pulses which have sharp boundaries. Moreover, it will be shown that these generalized, i.e., singular, pulses can be modified very simply into smooth ones and therefore can be implemented in real, physical scenarios.

5.1. The Step Function

The step function response of a dispersive medium is related to the complex error function. Let the initial pulse’s profile be a step function, i.e.,
A ( z = 0 , t ) = A 0 u ( t )
where u ( t ) = { 1 t > 0 0 t 0 is the Heaviside step-function, then the pulse at the end of the medium is
A ( t , z ) = A 0 0 ( 2 π i β 2 z ) 1 / 2 exp ( i ( t t ) 2 2 β 2 z ) d t = 1 2 A 0 e r f c ( t 2 i β 2 z )
Clearly, if initially
A ( z = 0 , t ) = A 0 u ( t + T )
then
A ( t , z ) = 1 2 A 0 e r f c ( t T 2 i β 2 z )
Similarly, using Equation (6), the initial boosted step function pulse
A ( t , z = 0 ) = A 0 exp ( i ω 0 t ) u ( t + T )
will propagate to read
A ( t , z > 0 ) = 1 2 A 0 exp ( i β 2 2 ω 0 2 z + i ω 0 t ) e r f c ( t T + β 2 ω 0 z 2 i β 2 z )
This result is known in the literature as Moshinsky’s function. In quantum mechanics, this function was used to simulate an electrons’ beam, which was blocked by a beam shutter [19,20]. In the quantum analogy, the wavefunction exhibits interference in time.

5.2. Rectangular Pulses

The fundamental generalized pulse is the rectangular one. Unlike the delta function, the rectangular pulse has a finite energy and therefore can simulate real pulses in many scenarios despite its singular edges (a real rectangular pulse), and therefore it is ubiquitous in the literature.
The rectangular pulse can be written with the aid of the step function or the rectangular function:
A ( z = 0 , t ) = A 0 r e c t 1 ( t / θ ) = A 0 [ u ( t / θ + 0.5 ) u ( t / θ 0.5 ) ]
where
r e c t L ( τ ) { 1 | τ | L / 2 0 | τ | > L / 2
Equation (25) has two singular points ( t = ± 0.5 θ ) in which the pulse’s amplitude is discontinuous. The dynamics of this rectangular pulse can be formulated analytically:
A ( t , z ) = 1 2 A 0 [ e r f c ( t θ / 2 2 i β 2 z ) e r f c ( t + θ / 2 2 i β 2 z ) ]
In Figure 6, the real and imaginary components of (27) are plotted for two different medium’s length. The ripples in the pulse’s profile appear after a very short distance.
In Figure 7, a spatiotemporal intensity plot (as a function of τ and ζ ) is presented.
The intricate structure of this function is most pronounced when the incident signal is a train of rectangular pulses:
A ( z = 0 , t ) = A 0 n = r e c t 1 ( t / θ 2 n )
This case is equivalent to the well-known quantum problem of a particle in a 1D box when initially the particle’s probability density is uniform over the box [21].
In this case; the problem’s dynamics can be formulated (using Equation (27)] by
A ( t , z ) = 1 2 A 0 n = [ e r f c ( t 2 n θ / 2 2 i β 2 z ) e r f c ( t 2 n + θ / 2 2 i β 2 z ) ]
However, it should be emphasized that in Ref. [21] a different numerical approach was taken in which no error function analysis was used.
Sample Application: In 10 Gb/s optical communications channel (over smf28 fiber) with a carrier wavelength of λ = 1.55   μ m , the pulses can be simulated by rectangular ones with θ = 100   ps and β 2 20   ps 2 / K m , in which case ζ = 0.3 in Figure 7 corresponds to z = ζ θ 2 / β 2 150   K m . However, as Figure 6 illustrates, considerable distortions appear even for ζ = 0.01 , which corresponds to a fiber length of z 5   Km .

5.3. Chirped Rectangular Pulses

One of the methods to combat dispersion in optical communication is to chirp the pulses, and since most pulses in optical communications are rectangular (in fact they are smooth rectangular, see below) it is instructive to investigate, analytically, the dynamics of a chirped rectangular pulse.
The initial chirped rectangular pulse has the following form:
A ( 0 , t ) = A 0 exp ( i q t 2 ) r e c t 1 ( t / θ ) = A 0 exp ( i q t 2 ) [ u ( t / θ + 0.5 ) u ( t / θ 0.5 ) ]
Therefore, using (10) and (24), the pulse profile after a distance z is
A ( z > 0 , t ) = exp ( i q t 2 2 β 2 z q + 1 ) A 0 2 1 + 2 β 2 z q [ e r f c ( t θ ( 1 + 2 β 2 z q ) / 2 2 i β 2 z ( 1 + 2 β 2 z q ) ) e r f c ( t + θ ( 1 + 2 β 2 z q ) / 2 2 i β 2 z ( 1 + 2 β 2 z q ) ) ]
Besides the widening of the pulse’s boundaries, the boundaries also move at a constant “velocity”, i.e., the two boundaries obey the following relationship:
t B = ± θ 2 ( 1 + 2 β 2 q z )
Clearly, when q is positive, the two boundaries collide at a distance:
z f = 1 / 2 β 2 q ,
which is independent of the pulse’s width θ . In Figure 8, the propagation of this pulse for both positive and negative q is presented.
While the pulse indeed shrinks for negative q , the pulse’s minimum width is reached before the distance z f , i.e., Equation (33). However, if the two pulses are combined, i.e., if the initial pulse is
A ( z = 0 , t ) = A 0 cos ( q t 2 ) r e c t 1 ( t / θ ) = A 0 cos ( q t 2 ) [ u ( t / θ + 0.5 ) u ( t / θ 0.5 ) ]
then the minimum width occurs exactly at the distance z f . In Figure 9, the temporal dynamics of this pulse is presented, and the point z f is clearly shown. Clearly, a smoother, i.e., more realistic, pulse will not converge to an infinitely temporally narrow pulse but instead will converge to a finite temporal width.

5.4. Exponential Pulse

The product of the step function with the exponential function allows the investigation of the dynamics of the simplest exponential function. Utilizing Equation (8), the initial exponential-step function
A ( t , z = 0 ) = A 0 exp ( a t ) A ( t , 0 )
will propagate in a dispersive medium according to
A ( z > 0 , t ) = A 0 2 exp ( i β 2 2 a 2 z + a t ) e r f c ( t i β 2 a z 2 i β 2 z )
In Figure 10, the dynamics of the exponential-step function pulse is plotted.
A spatiotemporal presentation of the pulse’s intensity is presented in Figure 11. This pulse can simulate a pulse where its rise time is considerably shorter than its decay time. In both Figures, a = 1 / θ .

5.5. Cosine Pulse

A bounded cosine pulse is an example of a continuous pulse (unlike the rectangular one), where the discontinuity at the edges of the pulse occurs on the amplitude’s derivative level.
It should be stressed that in terms of the pulse’s intensity, the derivative is continuous as well. Therefore, this is an important pulse because, for many practical purposes, it mimics a continuous pulse, but at the same time, it is (initially) temporally bounded.
In this case, if the initial pulse is
A ( z = 0 , t ) = { A cos ( π t / θ ) | t | < θ / 2 0 | t | θ / 2
which can be rewritten as
A ( z = 0 , t ) = A 2 { [ exp ( i π θ t ) + exp ( i π θ t ) ] u ( t + θ / 2 ) [ exp ( i π θ t ) + exp ( i π θ t ) ] u ( t θ / 2 ) }
then, using Equation (24), the pulse profile at the medium’s end is
A ( z , t ) = A 0 4 exp ( i β 2 z 2 ( π T ) 2 ) [ exp ( i π t T ) e r f c ( t T / 2 + π β 2 z / T 2 i β 2 z ) + exp ( i π t T ) e r f c ( t T / 2 π β 2 z / T 2 i β 2 z ) exp ( i π t T ) e r f c ( t + T / 2 + π β 2 z / T 2 i β 2 z ) exp ( i π t T ) e r f c ( t + T / 2 π β 2 z / T 2 i β 2 z ) ]
Figure 12 illustrates the dynamics of the pulse presented by Equation (39).
A spatiotemporal presentation of the pulse’s intensity is presented in Figure 13.

5.6. Square Cosine Pulse

A smoother pulse profile is the square cosine one. In this case, the discontinuity occurs in the second derivative level, but the function itself and its derivative are both continuous. Therefore, this pulse can simulate pulses, which were generated by discontinuous dielectric media. Let
A ( z = 0 , t ) = { A 0 cos 2 ( π t / θ ) | t | < θ / 2 0 | t | θ / 2
be the signal’s profile at one end of the medium; then, since it can be rewritten as
A ( z = 0 , t ) = 1 4 A 0 [ exp ( 2 i π θ t ) + exp ( 2 i π θ t ) + 2 ] [ u ( t + θ / 2 ) u ( t θ / 2 ) ]
its dynamics can be derived, again using Equation (24), to yield the result at the second end
A ( t , z ) = 1 4 A 0 [ e r f c ( t T / 2 2 i β 2 z ) e r f c ( t + T / 2 2 i β 2 z ) exp ( i β 2 z 2 ( 2 π T ) 2 ) [ exp ( i 2 π t T ) e r f c ( t T / 2 + 2 π β 2 z / T 2 i β 2 z ) + exp ( i 2 π t T ) e r f c ( t T / 2 2 π β 2 z / T 2 i β 2 z ) exp ( i π t T ) e r f c ( t + T / 2 + π β 2 z / T 2 i β 2 z ) exp ( i π t T ) e r f c ( t + T / 2 π β 2 z / T 2 i β 2 z ) ]
Figure 14 presents the dynamics of the square cosine pulse.
A comparison between Equations (27), (39), and (42) is presented in Figure 15. The three pulses represents three levels of singularities. Equation (27) represents a case of amplitude discontinuity, Equation (39) represents discontinuity in the pulse’s amplitude’s derivative, and Equation (42) represents discontinuity in the amplitude’s second derivative.
As can clearly be seen, the higher is the level of the singularity (i.e., the discontinuity occurs at higher derivatives), the faster the pulse decays in accordance with Ref. [22,23].

5.7. Generalization and Applicable Examples

The mathematical tools that were presented in the previous sections can be implemented in numerous practical examples. In what follows we will present some examples.
First example:
Generalization to any power of the bounded cosine pulse, i.e.,
A ( z = 0 , t ) = { A 0 cos n ( π t / T ) | t | < T / 2 0 | t | T / 2
which can be written as
A ( z = 0 , t ) = 1 2 n A 0 [ exp ( i π T t ) + exp ( i π T t ) ] n [ u ( t + T / 2 ) u ( t T / 2 ) ]
and since the brackets can be expanded
A ( z = 0 , t ) = 1 2 n A 0 k = 0 n ( n k ) exp [ i ( 2 k n ) π T t ] [ u ( t + T / 2 ) u ( t T / 2 ) ]
then, using Equation (24), the pulse’s profile at the end of the medium is
A ( z = 0 , t ) = 1 2 n A 0 × k = 0 n ( n k ) exp ( i β 2 z 2 ( 2 π T ) 2 ) ×          [ exp ( i ( 2 k n ) π t T ) e r f c ( t T / 2 + ( 2 k n ) π β 2 z / T 2 i β 2 z )          exp ( i ( 2 k n ) π t T ) e r f c ( t T / 2 ( 2 k n ) 2 π β 2 z / T 2 i β 2 z ) ]
Second example:
The initial singular pulse
A ( z = 0 , t ) = A 0 exp ( a | t | )
consists of two exponential pulses, and therefore its propagation can be expressed in a relatively simple form, namely:
A ( z > 0 , t ) = A 0 2 exp ( i β 2 2 a 2 z ) { exp ( a t ) e r f c ( t i β 2 a z 2 i β 2 z ) + exp ( a t ) e r f c ( t i β 2 a z 2 i β 2 z ) }
Third example:
In the previous example, we sewed two singular pulses; in general, however, one can sew several pulses together to create a combined pulse, which is much smoother than its singular components. For example, take the following pulse profile:
A ( z = 0 , t ) = { A 0 cos ( k t ) | t | < T / 2 A 0 cos ( k T / 2 ) exp [ k ( | t | T / 2 ) tan ( k T / 2 ) ] | t | T / 2
which consist of three singular and discontinuous pulses; however, the combined pulse and its derivative are both continuous. The discontinuity occurs only in the second derivative. In this case the pulse’s amplitude profile obeys
A z , t = A 0 4 exp i β 2 z 2 k 2 exp i k t e r f c t T / 2 + k β 2 z 2 i β 2 z + exp i k t e r f c t T / 2 k β 2 z 2 i β 2 z exp i k t e r f c t + T / 2 + k β 2 z 2 i β 2 z exp i k t e r f c t + T / 2 k β 2 z 2 i β 2 z + 1 2 A 0 cos k T / 2 exp i β 2 2 k tan k T / 2 2 z × exp k tan k T / 2 t e r f c t T / 2 i β 2 k tan k T / 2 z 2 i β 2 z exp k tan k T / 2 t e r f c t T / 2 i β 2 k tan k T / 2 z 2 i β 2 z
Fourth example:
The discontinuity of the exponential pulse can be reduced by multiplying the pulse by a sine function:
A ( z = 0 , t ) = A 0 exp ( a t ) sin ( b t ) u ( t )
Since the pulse can be written as
A ( z = 0 , t ) = i A 0 2 [ exp ( ( a + i b ) t ) exp ( ( a i b ) t ) ] u ( t )
then the pulse’s amplitude after a distance z reads
A ( t , z ) = i A 0 4 exp ( i β 2 2 ( a 2 b 2 ) z + a t ) × [ exp ( β 2 a b z + i b t ) e r f c ( t i β 2 z ( a + i b ) 2 i β 2 z ) exp ( β 2 a b z i b t ) e r f c ( t i β 2 z ( a i b ) 2 i β 2 z ) ]
In Figure 16, this pulse is presented for two cases with different parameters. On the right figure, the frequency is equal to the decaying coefficient a ; however, on the left figure, the frequency is higher than the decaying exponent a , and therefore the profile initially oscillates.

6. Smooth Pulses

Thus far, we have analyzed, except for the Gaussian pulses, only singular pulses, i.e., pulse profiles which have at least one singularity point. Most of them can simulate real physical pulse profiles; however, one may argue that these pulses fundamentally cannot behave like physical pulses. In some respects, it is indeed a justified claim, especially for a short medium. It has been shown that since the dispersion equation is not a causal equation, any singularity has nonlocal effects, which do not exist in real pulses [19,20,22,24]. However, it comes out that there is a very elegant method to convert each one of these singular pulses to a smooth one, which can simulate, in all respect, a real physical profile pulse. This method was first presented in Ref. [23,24] for the rectangular pulse.

6.1. Smooth Step Function

In principle, any smooth step function can replace the Heaviside step function; however, there is a clear advantage to choosing the complementary error function for that purpose. The function
1 2 e r f c ( t Δ )
is an approximation of the step function in the sense that
lim | t / Δ | 1 2 e r f c ( t Δ ) = u ( t )
However, unlike the step function, the transition in Equation (54) occurs within the time scale Δ .
Equivalently, a step function signal that passes through a Gaussian low-pass filter whose spectral FWHM is δ will have the following form [11]:
1 2 e r f c ( π δ t 2 ln 2 )
Therefore, the relation between the transition time scale Δ and the spectral FWHM δ is
δ = 2 ln ( 2 ) θ / π Δ .
Moreover, and this is the reason that Equation (54) is very useful in replacing step functions, the smooth step (Equation (54)) can be written as a convolution between the step function and a Gaussian:
1 2 e r f c ( t Δ ) = 1 π t / Δ exp ( x 2 ) d x = 1 Δ π u ( τ t ) exp ( τ 2 / Δ 2 ) d τ = u ( t ) 1 Δ π exp ( t 2 / Δ 2 )
where the asterisk represents temporal convolution. However, since the effect of dispersion can also be represented as a convolution with the kernel (Equation (4)), the propagation of the initial pulse’ profile
A ( t , 0 ) = A 0 2 e r f c ( t Δ )
can be written simply as
A ( t , 0 ) = A 0 u ( t ) 1 Δ π exp ( t 2 / Δ 2 ) K ( t , z ) = 1 2 A 0 e r f c ( t Δ 2 2 i β 2 z )
Therefore, every step function from the previous section can be replaced with the smooth step function. In particular, the boosted smooth step function
A ( z = 0 , t ) = A 2 exp ( i ω 0 t ) e r f c ( t Δ )
will have the following profile at the end of the medium:
A ( z > 0 , t ) = A 2 exp ( i β 2 2 ω 0 2 z i ω 0 t ) e r f c ( t β 2 ω 0 z Δ 2 2 i β 2 z ) .
Sample Applications: The rise time of pulses in optical communications are usually defined as the time period it takes to rise from 10% to 90% maximum intensity value. Therefore, provided this rise-time ( Δ t 10 90 ) is given, Δ can be derived by Δ 0.67 Δ t 10 90 . If the rise time is defined as the period in which the signal rises from 25% to 75% of its maximum intensity value ( Δ t 25 75 ), then Δ 1.28 Δ t 25 75 .

6.2. Smooth Rectangular Pulse

Similarly, the smooth rectangular pulse
A ( z = 0 , t ) = A 0 2 [ e r f c ( t θ / 2 Δ ) e r f c ( t + θ / 2 Δ ) ]
will appear at the end of the medium in the following form:
A ( z > 0 , t ) = A 0 srect 1 ( t θ , 2 ln ( 2 ) θ π Δ , β 2 z θ 2 ) = A 0 2 [ e r f c ( t θ / 2 2 i β 2 z + Δ 2 ) e r f c ( t + θ / 2 2 i β 2 z + Δ 2 ) ]
where
srect ξ ( τ , δ , ζ ) { e r f c [ ( τ ξ / 2 ) / i 2 ζ + 2 ln ( 2 ) / π 2 δ 2 ] e r f c [ ( τ + ξ / 2 ) / i 2 ζ + 2 ln ( 2 ) / π 2 δ 2 ] } / 2
is the smooth rectangular function (see Ref. [23,24,25]), which describes the dispersion dynamics of a smooth rectangular pulse.
In Figure 17, this smooth rectangular pulse’s profile is presented for the same values as in Figure 6, i.e., the perfect rectangular pulse. However, the difference is that in Figure 17, the pulse’s boundaries are smooth instead of discontinuous. As can be seen, due to the smooth transition, the oscillations, which appear in Figure 6, disappear in Figure 17.
This pulse has numerous applications in optical communications. For example, if the transmitted data is encoded in the amplitudes of these smooth rectangular pulses, then the transmitted signal profile can be written as
A ( t , z > 0 ) = n a n srect ξ ( t / θ n , 2 ln ( 2 ) θ π Δ , β 2 z / θ 2 )
where a n are coefficients which carry the data.
These pulses have two properties that make them especially suitable for stream of pulses. The ξ parameter can determine the signal’s protocol. The Return-to-Zero (RZ) protocol occurs for ξ < 1 , when ξ determines the duty cycle. However, when ξ = 1 , then this is a perfect Non-Return-to-Zero (NRZ) protocol. The second property, which is more important, emerges due to symmetry of the erfc function. Since
e r f c ( t Δ ) + e r f c ( t Δ ) = 2
then two adjacent smooth rectangular pulses look exactly like a single but twice as wide a pulse:
srect 1 ( τ 0.5 , δ , ζ ) + srect 1 ( τ + 0.5 , δ , ζ ) = srect 2 ( τ , δ , ζ )
Sample Application: In 40 Gb/s optical communication channels (over smf28 fiber) with a carrier wavelength of λ = 1.55   μ m , the pulses can be simulated by rectangular ones with θ = 25   ps , β 2 20   ps 2 / Km , in which case ζ = 0.3 and ζ = 0.01 in Figure 17 correspond to z = ζ θ 2 / β 2 9.4   Km and z 0.31   Km , respectively. Since Δ = 0.2 θ , then the 10–90% rise time corresponds to Δ t 10 90 Δ / 0.67 = 0.2 θ / 0.67 7.5   ps .

6.3. Relations to Super-Gaussian Pulses

The smooth rectangular pulses can be an excellent approximation for super-Gaussian ones. Let
A ( z = 0 , t ) = A 0 exp [ ( t / T ) n ]
be a super-Gaussian pulse. This pulse can be approximated by a smooth rectangular one (Equation (63)), A ( z = 0 , t ) = A 0 [ e r f c ( ( t θ / 2 ) / Δ ) e r f c ( ( t + θ / 2 ) / Δ ) ] / 2 , by choosing the following parameters:
θ = 2 T ( ln 2 ) 1 / n
and
Δ = 2 T π n ( ln 2 ) ( n 1 ) / n
In Figure 18, such a comparison is presented for two different values of the super-Gaussian order n . This figure clearly shows that the smooth rectangular is an excellent approximation to super-Gaussian, while, unlike super-Gaussian, the smooth rectangular pulse does have an exact analytical expression for its propagation in dispersive medium.

6.4. Chirped Smooth Rectangular Pulse

The case of the chirped rectangular pulse can be generalized to the chirped smooth rectangular one. That is, if the initial pulse’s profile is
A ( z = 0 , t ) = A 0 2 exp ( i q 2 t ) [ e r f c ( t θ / 2 Δ ) e r f c ( t + θ / 2 Δ ) ]
then, after a certain distance in the dispersive medium, the pulse’s profile will be
A ( z > 0 , t ) = ( 1 + 2 β 2 z q ) 1 / 2 exp ( i t 2 2 β 2 z + 1 / q ) A 0 2 [ e r f c ( t 1 + 2 β 2 z q θ / 2 2 i β 2 z 1 + 2 β 2 z q + Δ 2 ) e r f c ( t 1 + 2 β 2 z q + θ / 2 2 i β 2 z 1 + 2 β 2 z q + Δ 2 ) ]
A spatiotemporal presentation of the pulse’s intensity is presented in Figure 19. The fact that the initial pulse is smooth eliminates the ripples that appear in Figure 8.

6.5. Smooth Cosine Pulse

As was said above, the same procedure can be applied to every singular pulse, i.e., to replace the singular step function with the erfc functions. For example, the singular cosine pulse (Equation (37)) can be replaced with the following smooth cosine one:
A ( z = 0 , t ) = A 0 cos ( π θ t ) [ e r f c ( t θ / 2 Δ ) e r f c ( t + θ / 2 Δ ) ]
which has the following solution for any distance z :
A ( z , t ) = A 0 4 exp ( i β 2 z 2 ( π θ ) 2 ) [ exp ( i π t θ ) e r f c ( t θ / 2 + π β 2 z / θ Δ 2 2 i β 2 z ) + exp ( i π t θ ) e r f c ( t θ / 2 π β 2 z / θ Δ 2 2 i β 2 z ) exp ( i π t θ ) e r f c ( t + θ / 2 + π β 2 z / θ Δ 2 2 i β 2 z ) exp ( i π t θ ) e r f c ( t + θ / 2 π β 2 z / θ Δ 2 2 i β 2 z ) ]
Clearly, when Δ 0 , Equation (75) converges to Equation (39). In Figure 20, the amplitude profile of Equation (75) is presented.
Clearly, the same procedure can be implemented on the square cosine pulse or on any other singular pulse from the previous sections.

6.6. Smooth Exponential Pulse

By following the same procedure, one can replace the singular exponential with the following smooth exponential pulse:
A ( z = 0 , t ) = 1 2 A 0 exp ( a t ) e r f c ( t Δ )
which propagates in dispersive medium according to
A ( z > 0 , t ) = A 0 2 exp ( i β 2 2 a 2 z + a t ) e r f c ( t i β 2 a z Δ 2 2 i β 2 z ) .
This function is illustrated in Figure 21 for two different distances.
Sample Application: A pulse with a 10% to 90% rise time of Δ t 10 90 = 1   ps and a decay time of ~3.7 ps with a carrier wavelength of λ = 2.0   μ m can be simulated by the pulse presented by Equation (76) with Δ = 0.67   ps , a = 0.3   ps 1 . Since a Δ 0.2 , Figure 21 is an illustration of this pulse. If the pulse penetrates a BK7 glass, then for this wavelength, β 2 98   ps 2 / Km , and therefore ζ = 0.1 and ζ = 0.01 in the figure correspond to z = ζ θ 2 / β 2 11   m and z 1.1   m , respectively.

7. Singular Pulses in the Spectral Domain

7.1. The ideal Nyquist-Sinc Pulse

One of the interesting and important pulses, which belongs to this category and received a lot of attention recently in the literature, is the Nyquist-Sinc pulse (see, for example, Ref. [26,27,28,29,30,31]).
These pulses are useful in optical communications for several reasons, one of which is that they are resilient against chromatic dispersion. Moreover, these pulses, like their Fourier counterpart, the rectangular pulses, are orthogonal both in the time and in the frequency domains.
The initial Nyquist-Sinc pulse is singular in the spectral domain; its Fourier transform is:
A ( z = 0 , ω ) = θ A 0 r e c t 1 ( ω θ / 2 π )
In the temporal regime the “sinc” profile appears:
A ( z = 0 , t ) = A 0 sinc ( t θ ) A 0 sin ( π t / θ ) π t / θ
Since the pulse’s spectrum has a simple rectangular shape, the temporal dynamics is straightforward, namely:
A ( z , t ) = A 0 dsinc ( τ , ζ ) A 0 2 i θ 2 2 π β 2 z exp ( i t 2 2 β 2 z ) [ erf ( t | π β 2 z / θ | i 2 β 2 z ) erf ( t + | π β 2 z / θ | i 2 β 2 z ) ]
where
dsinc ( τ , ζ ) 1 2 i 2 π ζ exp ( i τ 2 2 ζ ) [ erf ( τ π | ζ | i 2 ζ ) erf ( τ + π | ζ | i 2 ζ ) ]
is the “dynamic” sinc function [12]. Clearly, dsinc ( τ , ζ = 0 ) = sinc ( τ ) sin ( π τ ) π τ .
In Figure 22, the amplitude profile of this pulse is presented.
The spatiotemporal presentation of the pulse’s intensity is presented in Figure 23.

7.2. Nyquist Sinc Pulse with Smooth Spectrum

The singular spectrum may simulate real pulses, but in real scenarios, singularities cannot occur even in the spectral domain. We can therefore utilize the theorems and conclusions from the previous sections to spectrally smooth the Nyquist-Sinc pulse.
Thus, the singular rectangular spectrum of the Nyquist-Sinc pulse
A ( z = 0 , ω ) = θ A 0 r e c t 1 ( ω θ / 2 π )
can be replaced with its smooth counterpart:
A ( z = 0 , ω ) = θ A 0 1 2 [ e r f c ( ω π / θ Δ ) e r f c ( ω + π / θ Δ ) ]
However, this expression can be written as a convolution and is the spectral domain:
A ( z = 0 , ω ) = θ A 0 r e c t 1 ( ω θ / 2 π ) 1 Δ π exp [ ( ω / δ ) 2 ]
and therefore the initial pulse profile can be written as a product of a sinc function and a Gaussian one:
A ( z = 0 , t ) = A 0 sin ( π t / θ ) π t / θ 1 2 π exp [ ( t δ / 2 ) 2 ]
Using Equations (12) and (80), we finally obtain
A ( t , z ) = A 0 2 exp ( i t 2 2 β 2 z ) i θ 2 2 π β 2 z [ erf ( t | π β 2 z / θ | i 2 β 2 z ( 1 i β 2 z δ 2 / 2 ) ) erf ( t + | π β 2 z / θ | i 2 β 2 z ( 1 i β 2 z δ 2 / 2 ) ) ]
In Figure 24, the amplitude of Equation (86) is plotted. Unlike ordinary Nyquist-Sinc pulses, this pulse decays exponentially.

8. Undistorted Airy Pulses

8.1. Undistorted Ideal Accelerating Pulses

One of the interesting properties of the dispersion equation is that the Airy function is its accelerating solution [32,33,34]. However, despite dispersion, the shape of this accelerating solution remains intact during propagation. This property can be useful in optical communications where the encoded data propagates in the dispersive medium without experiencing distortions.
That is, if the initial pulse is an Airy function pulse:
A ( t , z ) = A 0 Ai ( t / θ )
where Ai ( τ ) stands for the Airy function [35], then the pulse at the end of the dispersive medium keeps its initial profile despite being accelerated; after a distance z, it reads:
A ( t , z ) = A 0 A i [ t / θ ( β 2 z / 2 θ 2 ) 2 ] exp [ i t β 2 z / 2 θ 3 + i ( β 2 z / θ 2 ) 3 / 12 ]
In Figure 25, the dynamics of this pulse is presented for different distances.
As can be seen in Figure 25, the pulse’s amplitude shape | A ( t , z ) / A 0 | = A i [ t / θ ( β 2 z / 2 θ 2 ) 2 ] remains intact, i.e., it keeps the same Airy shape. The only difference is that the pulse’s temporal coordinate [ t / θ ] has shifted with the distance to t / θ ( β 2 z / 2 θ 2 ) 2 , and therefore the pulse’s temporal location accelerates according to the equation
t = ( β 2 z ) 2 / 4 θ 3 .
Therefore, the square of the dispersion coefficient is proportional to the pulse’s “acceleration”.

8.2. Physical Accelerating Pulse

Clearly, the main problem with the Airy function is that it is not a normalizable function, i.e., it has infinite energy and can have no physical realization.
One option to normalize the airy function is to multiply the pulse by an exponent [36]:
A ( z = 0 , t ) = A 0 exp ( a t ) Ai ( t / θ )
In which case, the pulse dynamics can be derived with the aid of Equations (8) and (88), yielding:
A ( z , t ) = A 0 A i [ t i β 2 a z θ ( β 2 z 2 θ 2 ) 2 ] exp [ i β 2 2 a 2 z + a t i β 2 z ( t i β 2 a z ) 2 θ 3 + i 12 ( β 2 z θ 2 ) 3 ]
In Figure 26, the dynamics of this pulse is presented. A spatiotemporal presentation of the pulse’s intensity is presented in Figure 27. This pulse has a finite energy (it can be normalized), but as a consequence, the pulse decays in the z-direction as well. Since
| A ( z , t ) | 2 = | A 0 | 2 | A i [ t i β 2 a z θ ( β 2 z 2 θ 2 ) 2 ] | 2 exp ( 2 a t a ( β 2 z ) 2 θ 3 )
the pulse’s confinement within the time-scale a 1 “costs” resultant confinement in the z-direction, i.e., the pulse decays within a distance of z θ 3 / a / β 2 .

8.3. Attenuation Compensating Airy Pulse

It has been realized that the Airy beams can be implemented to compensate medium’s attenuation. Energy can flow from minus infinity to compensate for the pulse’s attenuation. Such a pulse can be constructed by imaginary period displacement [37,38]:
A ( t , z ) = A 0 Ai [ ( t i η ) / θ ]
In this case, the pulse at the end of the dispersive medium can be directly derived:
A ( t , z ) = A 0 Ai [ ( t i η ) / θ ( β 2 z / 2 θ 2 ) 2 ] exp [ i ( t i η ) β 2 z / 2 θ 3 i ( β 2 z / θ 2 ) 3 / 12 ]
This expression is presented in Figure 28.
The introduction of η is responsible to the additional exponential term exp [ η β 2 z / 2 θ 3 ] , which can compensate for the medium’s absorption.
The problem with these pulses is, again, that they have infinite energy and cannot be normalized. Therefore, only approximations of these pulses can be implemented in real experiments (see, for example, Ref. [38]).

8.4. Physical Attenuation Compensating Airy Pulse

In the previous section, we have seen that a complex Airy pulse can compensate pulse attenuation. However, these pulses cannot be normalized. One option to turn these pulses to a finite energy ones, and therefore normalizable, is to multiply them by an exponential function, i.e., to choose the following initial pulse profile:
A ( z , t ) = A 0 Ai ( t i η θ ) exp ( a t )
Then, following Equations (8) and (93), the pulse after a distance z has the following form:
A ( z , t ) = A 0 Ai [ t i η + i β 2 a z θ ( β 2 z 2 θ 2 ) 2 ] exp [ i β 2 2 a 2 z + a t + i β 2 z ( t i η + i β 2 a z ) 2 θ 3 i 12 ( β 2 z θ 2 ) 3 ]
The intensity of this pulse
| A ( z , t ) | 2 = | A 0 Ai [ t i η + i β 2 a z θ ( β 2 z 2 θ 2 ) 2 ] | 2 exp [ 2 a t β 2 z ( β 2 a z η ) θ 3 ]
is presented in Figure 29 and Figure 30. One can easily see that unlike Figure 25, the pulse’s amplitude increases in short distances, but unlike Figure 28, it has finite energy, and therefore, unlike Equation (94), can be implemented in real physical scenarios.
Clearly, however, the pulse’s amplitude will eventually decay. It will reach its maximum intensity at z = η / β 2 a , beyond which it will decay to zero. Therefore, this pulse may be implemented to compensate for the medium’s losses when the width of the medium is a given parameter.

9. Pulse Broadening Comparison

To illustrate the large differences between equivalent pulses, the FWHM of four different pulses is presented in Figure 31. To keep the initial FWHM of all pulses equal, the pulses’ parameters were adjusted accordingly. For an arbitrary initial FWHM of T, the FWHM as a function of the normalized distance ζ = β 2 z / T 2 is presented for Gaussian pulse, i.e., Equation (14) with the parameter θ = 0.85 T ; for smooth rectangular pulse, i.e., Equation (64) with θ = 1.05 T and Δ = 0.1 T ; for smooth exponential pulse, i.e., Equation (77) with Δ = 0.1 T and a = 0.4 T 1 ; and finally for Nyquist-Sinc pulse, i.e., Equation (80), with θ = 1.1 T . Clearly, the differences will increase if the signals are chirped.

10. Discussion and Conclusions

Dispersive media can have a destructive effect on ultrashort pulses even at relatively short distances.
This paper shows that almost any pulse shape in most practical scenarios can be formulated in a pulse form, whose propagation in dispersive medium has an exact analytical expression.
Therefore, there is no need to restrict the analysis of pulse propagation in a dispersive medium to Gaussian pulses or to numerical analysis. Thus, almost any shape and form of pulses can be analyzed analytically. This is the main result and main conclusion of this paper.
The experimentalist can choose from the multiple examples presented in this paper to predict the distortions that the relevant pulse would experience.
In particular, if the pulses were generated by a mode-locked laser, then a Gaussian pulse is a good approximation. However, if there is an asymmetry between the pulse’s rise and fall times (as in Q-switch lasers), then the experimentalist can use the smooth exponential pulses (Section 6.6) to predict the pulse’s distortions.
Similarly, in modern optical communications channels, the data signals consist of multiple pulses which mostly belong to two categories: rectangular and Nyquist pulses. Both types of pulses are presented in this paper; however, these pulses are both discontinuous: The former is discontinuous in the time domain while the latter is discontinuous in the frequency domain. In practical cases, the signals are always continuous (in both time and frequency domains). Therefore, the experimentalist or the engineer can measure the pulse’s rise time and apply the smooth rectangular pulses (Section 6.2) to the relevant system to predict the pulse’s distortion in the fiber. Similarly, Nyquist-Sinc pulses with a smooth spectrum (Section 7.2) can be used to simulate real systems, which are based on Nyquist-Sinc pulses. The exact dispersion’s distortion can be evaluated without the need for approximations or for numerical analysis.
Moreover, this paper demonstrates that any singular pulse can be replaced with a smooth pulse whose evolution in dispersive medium can be formulated in an exact analytical form.
Besides this generic recipe, the main new findings are the exact analytical derivations of the following pulses: (A) chirped rectangular pulses, which converge to extremely narrow pulses; (B) excellent analytical approximation of super-Gaussian pulses; (C) Nyquist-Sinc Pulses with a smooth spectrum; and (D) a physical realization of attenuation-compensating Airy pulses.

Funding

This research received no external funding.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A. Proof of Equation (3)

The spectral transfer function of the medium is
H ( ω ) = exp ( i β 2 z ω 2 / 2 )
i.e., for the incident signal A ( t , z = 0 ) = A 0 exp ( ± i ω t ) the signal after a distance of z is simply A ( t , z > 0 ) = A 0 H ( ω ) exp ( ± i ω t ) . In general, let A ( ω , z ) = F { A ( t , z ) } be the Fourier transform of A ( t , z ) then
A ( ω , z > 0 ) = H ( ω ) A ( ω , 0 )
Since the Fourier transform of a product is a convolution of the Fourier transforms, then the application of the Fourier transform on Equation (A2) yields
A ( t , z ) = K ( t t , z ) A ( t , 0 ) d t
where
K ( t , z ) = F { H ( ω ) } = ( 2 π i β 2 z ) 1 / 2 exp ( i t 2 2 β 2 z )
is the Fourier transform of H ( ω ) .

Appendix B. Proof of Equation (6)

Let
A ( t , z = 0 ) = exp ( i ω 0 t ) A ( t , 0 )
and let us further assume that
A ( t , z ) = K ( t t , z ) A ( t , 0 ) d t
is given. Then
A ( t , z ) = K ( t t , z ) exp ( i ω 0 t ) A ( t , 0 ) d t
which, after some calculations, can be rewritten as
A ( t , z ) = exp ( i β 2 z ω 0 2 2 + i ω 0 t ) ( 2 π i β 2 z ) 1 / 2 exp ( i [ t ( t + β 2 z ω 0 ) ] 2 2 β 2 z ) A ( t , 0 ) d t
Using Equations (A6) and (A4), Equation (A8) is simply
A ( t , z > 0 ) = exp ( i β 2 2 ω 0 2 z + i ω 0 t ) A ( t + β 2 ω 0 z , z ) .

Appendix C. Proof of Equation (10)

Let
A ( t , 0 ) = exp ( i q t 2 ) A ( t , 0 )
and let us further assume that
A ( t , z ) = K ( t t , z ) A ( t , 0 ) d t
is given. Then
A ( t , z ) = ( 2 π i β 2 z ) 1 / 2 exp ( i ( t t ) 2 2 β 2 z ) A ( t , 0 ) exp ( i q t 2 ) d t
Which, after some algebra, can be written as
A ( t , z ) = ( 1 + 2 β 2 z q 1 + 2 β 2 z q ) 1 / 2 exp ( i t 2 2 β 2 z + 1 / q ) ( 2 π i β 2 z ) 1 / 2 exp ( i ( 1 + 2 β 2 z q 2 β 2 z ) ( t t 1 + 2 β 2 z q ) 2 ) A ( t , 0 ) d t
Using Equation (A11), Equation (A13) can simply be rewritten as
A ( t , z ) = ( 1 + 2 β 2 z q ) 1 / 2 exp ( i t 2 2 β 2 z + 1 / q ) A ( t 1 + 2 β 2 z q , z 1 + 2 β 2 z q )

References

  1. Zevallos, M.E.; Gayen, S.K.; Das, B.B.; Alrubaiee, M.; Alfano, R.R. Picosecond Electronic Time-Gated Imaging of Bones in Tissues. IEEE J. Sel. Top. Quantum Electron. 1999, 5, 916–922. [Google Scholar] [CrossRef]
  2. Gayen, S.K.; Alfano, R.R. Emerging optical biomedical imaging techniques. Opt. Photon. News 1996, 7, 17–22. [Google Scholar] [CrossRef]
  3. Das, B.B.; Yoo, K.M.; Alfano, R.R. Ultrafast time-gated imaging in thick tissues: A step toward optical mammography. Opt. Lett. 1993, 18, 1092–1094. [Google Scholar] [CrossRef] [PubMed]
  4. Marom, D.M.; Sun, P.C.; Fainman, Y. Communication with ultrashort pulses and parallel-to-serial and serial-to-parallel converters. In Proceedings of the LEOS ‘97, 10th Annual Meeting IEEE Lasers and Electro-Optics Society, San Francisco, CA, USA, 10–13 November 1997. [Google Scholar]
  5. Amiri, I.S.; Ahmad, H. Optical Soliton Communication Using Ultra-Short Pulses; Springer: Singapore, 2015. [Google Scholar]
  6. Yamaoka, Y.; Harada, Y.; Sakakura, M.; Minamikawa, T.; Nishino, S.; Maehara, S.; Hamano, S.; Tanaka, H.; Takamatsu, T. Photoacoustic microscopy using ultrashort pulses with two different pulse durations. Opt. Express 2014, 22, 17063–17072. [Google Scholar] [CrossRef] [PubMed]
  7. Gibbs, H.C.; Arne, Y.B.; Alvin, C.L.; Yeh, T. Imaging embryonic development with ultrashort pulse microscopy. Opt. Eng. 2014, 53, 051506. [Google Scholar] [CrossRef]
  8. Technical Note: The Effect of Dispersion on Ultrashort Pulses. Newport Corporation (2018). Available online: https://www.newport.com/n/the-effect-of-dispersion-on-ultrashort-pulses (accessed on 10 January 2018).
  9. Sindhu, T.G.; Bisht, P.B.; Rajesh, R.J.; Satyanarayana, M.V. Effect of higher order nonlinear dispersion on ultrashort pulse evolution in a fiber laser. Microw. Opt. Technol. Lett. 2001, 28, 196–198. [Google Scholar] [CrossRef]
  10. Wang, W.; Liu, Y.; Xi, P.; Ren, Q. Origin and effect of high-order dispersion in ultrashort pulse multiphoton microscopy in the 10 fs regime. Appl. Opt. 2010, 49, 6703–6709. [Google Scholar] [CrossRef] [PubMed]
  11. Granot, E. Fundamental dispersion limit for spectrally bounded On-Off-Keying communication channels and its implications to Quantum Mechanics and the Paraxial Approximation. Europhys. Lett. 2012, 100, 44004. [Google Scholar] [CrossRef]
  12. Granot, E. Information Loss in Quantum Dynamics. In Advanced Technologies of Quantum Key Distribution; INTECH: Rijeka, Croatia, 2017. [Google Scholar]
  13. Wollenhaupt, M.; Assion, A.; Baumert, T. Femtosecond Laser Pulses: Linear Properties, Manipulation, Generation and Measurement. Chap. 12. In Hanbookd of Laser and Optics; Träger, F., Ed.; Springer: New York, NY, USA, 2007. [Google Scholar]
  14. Agrawal, G.P. Fiber-Optic Communications Systems, 3rd ed.; John Wiley & Sons, Inc.: New York, NY, USA, 2002. [Google Scholar]
  15. Crank, J. The Mathematics of Diffusion; Clarendon Press: Oxford, UK, 1975. [Google Scholar]
  16. Ali, R.; Hamza, M.Y. Propagation behavior of super-Gaussian pulse in dispersive and nonlinear regimes of optical communication systems. In Proceedings of the International Conference on Emerging Technologies (ICET 2016), Islamabad, Pakistan, 18–19 October 2016. [Google Scholar]
  17. Anderson, D.; Lisak, M. Propagation characteristics of frequency-chirped super-Gaussian optical pulses. Opt. Lett. 1986, 11, 569–571. [Google Scholar] [CrossRef]
  18. Zhang, L.; Li, C.; Zhong, H.; Xu, C.; Lei, D.; Li, Y.; Fan, D. Propagation dynamics of super-Gaussian beams in fractional Schrödinger equation: From linear to nonlinear regimes. Opt. Express 2016, 24, 14406–14418. [Google Scholar] [CrossRef]
  19. Moshinsky, M. Diffraction in time. Phys. Rev. 1952, 88, 625–631. [Google Scholar] [CrossRef]
  20. Del Campo, A.; Garcia-Calderon, G.; Muga, J.G. Quantum Transients. Phys. Rep. 2009, 476, 1–50. [Google Scholar] [CrossRef]
  21. Berry, M.V. Quantum fractals in boxes. J. Phys. A Math. Gen. 1996, 29, 6617–6629. [Google Scholar] [CrossRef]
  22. Granot, E.; Marchewka, A. Generic Short-Time Propagation of Sharp-Boundaries Wave Packets. Europhys. Lett. 2005, 72, 341–347. [Google Scholar] [CrossRef]
  23. Granot, E.; Luz, E.; Marchewka, A. Generic pattern formation of sharp-boundaries pulses propagation in dispersive media. J. Opt. Soc. Am. B 2012, 29, 763–768. [Google Scholar] [CrossRef]
  24. Granot, E.; Marchewka, A. Emergence of currents as a transient quantum effect in nonequilibrium systems. Phys. Rev. A 2011, 84, 032110–032115. [Google Scholar] [CrossRef]
  25. Marciano, S.; Ben-Ezra, S.; Granot, E. Eavesdropping and Network Analyzing Using Network Dispersion. Appl. Phys. Res. 2015, 7, 27. [Google Scholar] [CrossRef]
  26. Soto, M.A.; Alem, M.; Shoaie, M.A.; Vedadi, A.; Brès, C.-S.; Thévenaz, L.; Schneider, T. Optical sinc-shaped Nyquist pulses of exceptional quality. Nat. Commun. 2013, 4, 2898. [Google Scholar] [CrossRef]
  27. Schmogrow, R.; Bouziane, R.; Meyer, M.; Milder, P.A.; Schindler, P.C.; Killey, R.I.; Bayvel, P.; Koos, C.; Freude, W.; Leuthold, J. Real-time OFDM or Nyquist pulse generation—Which performs better with limited resources? Opt. Express 2012, 20, B543. [Google Scholar] [CrossRef]
  28. Hirooka, T.; Ruan, P.; Guan, P.; Nakazawa, M. Highly dispersion-tolerant 160 Gbaud optical Nyquist pulse TDM transmission over 525 km. Opt. Express 2012, 20, 15001–15007. [Google Scholar] [CrossRef]
  29. Hirooka, T.; Nakazawa, M. Linear and nonlinear propagation of optical Nyquist pulses in fibers. Opt. Express 2012, 20, 19836–19849. [Google Scholar] [CrossRef] [PubMed]
  30. Schmogrow, R.; Hillerkuss, D.; Wolf, S.; Bäuerle, B.; Winter, M.; Kleinow, P.; Nebendahl, B.; Dippon, T.; Schindler, P.C.; Koos, C.; et al. 512QAM Nyquist sinc-pulse transmission at 54 Gbit/s in an optical bandwidth of 3 GHz. Opt. Express 2012, 20, 6439–6447. [Google Scholar] [CrossRef] [PubMed]
  31. Bosco, G.; Carena, A.; Curri, V.; Poggiolini, P.; Forghieri, F. Performance limits of Nyquist-WDM and CO-OFDM in high-speed PM-QPSK systems. IEEE Phot. Technol. Lett. 2010, 22, 1129–1131. [Google Scholar] [CrossRef]
  32. Berry, M.V.; Balázs, N.L. Nonspreading wave packets. Am. J. Phys. 1979, 47, 264–267. [Google Scholar] [CrossRef]
  33. Siviloglou, G.A.; Broky, J.; Dogariu, A.; Christodoulides, D.N. Observation of Accelerating Airy Beams. Phys. Rev. Lett. 2007, 99, 213901. [Google Scholar] [CrossRef] [PubMed]
  34. Bandres, M.A. Accelerating beams. Opt. Lett. 2009, 34, 3791–3793. [Google Scholar] [CrossRef]
  35. Abramowitz, M.; Stegun, A. Handbook of Mathematical Functions; Dover Publications: New York, NY, USA, 1965. [Google Scholar]
  36. Siviloglou, G.A.; Christodoulides, D.N. Accelerating finite energy Airy beams. Opt. Lett. 2007, 32, 979. [Google Scholar] [CrossRef] [PubMed]
  37. Preciado, M.A.; Dholakia, K.; Mazilu, M. Generation of attenuation-compensating Airy beams. Opt. Lett. 2014, 39, 4950–4953. [Google Scholar] [CrossRef]
  38. Preciado, M.A.; Sugden, K. Proposal and design of airy-based rocket pulses for invariant propagation in lossy dispersive media. Opt. Lett. 2012, 37, 4970–4972. [Google Scholar] [CrossRef]
Figure 1. The real (upper panel) and imaginary (lower panel) components of the pulse (Equation (14)). The dashed curve represents the initial profile (for ζ = β 2 z / θ 2 = 0 ), while the solid curve represents the signal after a distance, which corresponds to ζ = β 2 z / θ 2 = 2 . Time is measured in units of θ , which corresponds to the pulse’s temporal width.
Figure 1. The real (upper panel) and imaginary (lower panel) components of the pulse (Equation (14)). The dashed curve represents the initial profile (for ζ = β 2 z / θ 2 = 0 ), while the solid curve represents the signal after a distance, which corresponds to ζ = β 2 z / θ 2 = 2 . Time is measured in units of θ , which corresponds to the pulse’s temporal width.
Applsci 09 00527 g001
Figure 2. A false-color presentation of the pulse’s intensity | A ( t , z ) | 2 (Equation (14)) as a function of the normalized time τ t / θ and the normalized distance ζ β 2 z / θ 2 .
Figure 2. A false-color presentation of the pulse’s intensity | A ( t , z ) | 2 (Equation (14)) as a function of the normalized time τ t / θ and the normalized distance ζ β 2 z / θ 2 .
Applsci 09 00527 g002
Figure 3. The real (upper panel) and imaginary (lower panel) components of the Gaussian pulse (Equation (16)). The dashed curve represents the initial profile (for ζ = β 2 z / θ 2 = 0 ), while the solid curve represents its final shape (for ζ = β 2 z / θ 2 = 2 ). In this example, the carrier frequency is ω 0 = 4 / θ .
Figure 3. The real (upper panel) and imaginary (lower panel) components of the Gaussian pulse (Equation (16)). The dashed curve represents the initial profile (for ζ = β 2 z / θ 2 = 0 ), while the solid curve represents its final shape (for ζ = β 2 z / θ 2 = 2 ). In this example, the carrier frequency is ω 0 = 4 / θ .
Applsci 09 00527 g003
Figure 4. Same as Figure 2, but for the intensity ( | A ( t , z ) | 2 ) of the pulse presented by Equation (16). (In this example, ω 0 = 4 / θ ).
Figure 4. Same as Figure 2, but for the intensity ( | A ( t , z ) | 2 ) of the pulse presented by Equation (16). (In this example, ω 0 = 4 / θ ).
Applsci 09 00527 g004
Figure 5. Same as Figure 2, but for the intensity ( | A ( t , z ) | 2 ) of the pulse presented by Equation (18) (in this example, q = 2 / θ 2 ).
Figure 5. Same as Figure 2, but for the intensity ( | A ( t , z ) | 2 ) of the pulse presented by Equation (18) (in this example, q = 2 / θ 2 ).
Applsci 09 00527 g005
Figure 6. The real (upper panel) and imaginary (lower panel) components of the rectangular pulse (Equation (27)). The dashed curve represents the initial profile (for ζ = β 2 z / θ 2 = 0 ), while the solid curve represents its final shape (for ζ = β 2 z / θ 2 = 0.01 on the right and ζ = β 2 z / θ 2 = 0.3 on the left).
Figure 6. The real (upper panel) and imaginary (lower panel) components of the rectangular pulse (Equation (27)). The dashed curve represents the initial profile (for ζ = β 2 z / θ 2 = 0 ), while the solid curve represents its final shape (for ζ = β 2 z / θ 2 = 0.01 on the right and ζ = β 2 z / θ 2 = 0.3 on the left).
Applsci 09 00527 g006
Figure 7. Same as Figure 2, but for the intensity ( | A ( t , z ) | 2 ) of the pulse presented by Equation (27).
Figure 7. Same as Figure 2, but for the intensity ( | A ( t , z ) | 2 ) of the pulse presented by Equation (27).
Applsci 09 00527 g007
Figure 8. Same as Figure 2, but for the intensity ( | A ( t , z ) | 2 ) of the pulse presented by Equation (31). On the left, q = 4 / θ 2 ; and on the right, q = 4 / θ 2 . The dashed lines correspond to the pulse’s boundaries t B = ± θ ( 1 + 2 β 2 q z ) / 2 .
Figure 8. Same as Figure 2, but for the intensity ( | A ( t , z ) | 2 ) of the pulse presented by Equation (31). On the left, q = 4 / θ 2 ; and on the right, q = 4 / θ 2 . The dashed lines correspond to the pulse’s boundaries t B = ± θ ( 1 + 2 β 2 q z ) / 2 .
Applsci 09 00527 g008
Figure 9. Same as Figure 8 but with the initial pulse Equation (34) for the parameter q = 4 / θ 2 .
Figure 9. Same as Figure 8 but with the initial pulse Equation (34) for the parameter q = 4 / θ 2 .
Applsci 09 00527 g009
Figure 10. The real (upper panel) and imaginary (lower panel) components of the exponential-step function pulse (Equation (36)). The dashed curve represents the initial profile (for ζ = β 2 z / θ 2 = 0 ), while the solid curve represents its final shape (for ζ = β 2 z / θ 2 = 0.01 ).
Figure 10. The real (upper panel) and imaginary (lower panel) components of the exponential-step function pulse (Equation (36)). The dashed curve represents the initial profile (for ζ = β 2 z / θ 2 = 0 ), while the solid curve represents its final shape (for ζ = β 2 z / θ 2 = 0.01 ).
Applsci 09 00527 g010
Figure 11. Same as Figure 2, but for the intensity ( | A ( t , z ) | 2 ) of the pulse presented by Equation (36).
Figure 11. Same as Figure 2, but for the intensity ( | A ( t , z ) | 2 ) of the pulse presented by Equation (36).
Applsci 09 00527 g011
Figure 12. Similar to Figure 10, but for the bounded cosine pulse (Equation (39)) and for the final distance of ζ = β 2 z / θ 2 = 0.3 .
Figure 12. Similar to Figure 10, but for the bounded cosine pulse (Equation (39)) and for the final distance of ζ = β 2 z / θ 2 = 0.3 .
Applsci 09 00527 g012
Figure 13. Same as Figure 2, but for the intensity ( | A ( t , z ) | 2 ) of the pulse presented by Equation (39).
Figure 13. Same as Figure 2, but for the intensity ( | A ( t , z ) | 2 ) of the pulse presented by Equation (39).
Applsci 09 00527 g013
Figure 14. Similar to Figure 12, but for the square cosine pulse (Equation (42)).
Figure 14. Similar to Figure 12, but for the square cosine pulse (Equation (42)).
Applsci 09 00527 g014
Figure 15. Comparison between the intensities of the three singular pulses represented by Equation (27)—dashed curve, Equation (39)—solid curve, and Equation (42)—dotted curve in a logarithmic scale.
Figure 15. Comparison between the intensities of the three singular pulses represented by Equation (27)—dashed curve, Equation (39)—solid curve, and Equation (42)—dotted curve in a logarithmic scale.
Applsci 09 00527 g015
Figure 16. Similar to Figure 10, but for the pulse presented by Equation (53). In these plots, a = 1 / θ on both, but b = 4 / θ on the left figure (final distance corresponds to ζ = β 2 z / θ 2 = 0.1 ) and b = 1 / θ on the right one (final distance corresponds to ζ = β 2 z / θ 2 = 0.2 ).
Figure 16. Similar to Figure 10, but for the pulse presented by Equation (53). In these plots, a = 1 / θ on both, but b = 4 / θ on the left figure (final distance corresponds to ζ = β 2 z / θ 2 = 0.1 ) and b = 1 / θ on the right one (final distance corresponds to ζ = β 2 z / θ 2 = 0.2 ).
Applsci 09 00527 g016
Figure 17. Same as Figure 12, but for Equation (64) with the transition width Δ = 0.2 θ .
Figure 17. Same as Figure 12, but for Equation (64) with the transition width Δ = 0.2 θ .
Applsci 09 00527 g017
Figure 18. A comparison between smooth rectangular pulses (Equation (63), solid curves) and super-Gaussian pulses (Equation (69), dashed curves) for n = 4 (right) and n = 14 (left). In these plots, only the real part of the fields is presented. The imaginary part is zero.
Figure 18. A comparison between smooth rectangular pulses (Equation (63), solid curves) and super-Gaussian pulses (Equation (69), dashed curves) for n = 4 (right) and n = 14 (left). In these plots, only the real part of the fields is presented. The imaginary part is zero.
Applsci 09 00527 g018
Figure 19. Same as Figure 8 but for a smooth rectangular pulse with Δ = 0.25 θ .
Figure 19. Same as Figure 8 but for a smooth rectangular pulse with Δ = 0.25 θ .
Applsci 09 00527 g019
Figure 20. Same as Figure 12, but for Equation (75) with the transition width Δ = 0.4 θ .
Figure 20. Same as Figure 12, but for Equation (75) with the transition width Δ = 0.4 θ .
Applsci 09 00527 g020
Figure 21. Same as Figure 10, but for Equation (77) with the transition width a = 1 / θ and Δ = 0.2 θ for two final distances ζ = β 2 z / θ 2 = 0.01 on the left and ζ = β 2 z / θ 2 = 0.1 on the right.
Figure 21. Same as Figure 10, but for Equation (77) with the transition width a = 1 / θ and Δ = 0.2 θ for two final distances ζ = β 2 z / θ 2 = 0.01 on the left and ζ = β 2 z / θ 2 = 0.1 on the right.
Applsci 09 00527 g021
Figure 22. Same as Figure 10, but for Equation (80) and for the final distance of ζ = β 2 z / θ 2 = 0.3 .
Figure 22. Same as Figure 10, but for Equation (80) and for the final distance of ζ = β 2 z / θ 2 = 0.3 .
Applsci 09 00527 g022
Figure 23. Same as Figure 2, but for the intensity ( | A ( t , z ) | 2 ) of the pulse presented by Equation (80).
Figure 23. Same as Figure 2, but for the intensity ( | A ( t , z ) | 2 ) of the pulse presented by Equation (80).
Applsci 09 00527 g023
Figure 24. Same as Figure 22, but for Equation (86) with δ = 0.6 / θ .
Figure 24. Same as Figure 22, but for Equation (86) with δ = 0.6 / θ .
Applsci 09 00527 g024
Figure 25. Presentation of the temporal dynamics of the accelerating Airy pulse (Equation (88)). On the left, the real and imaginary components of the pulse’s field are presented (for the final distance of ζ = β 2 z / θ 2 = 1 ), and on the right, the pulse’s intensity ( | A ( t , z ) | 2 ) is presented (for the final distance of ζ = β 2 z / θ 2 = 2 ).
Figure 25. Presentation of the temporal dynamics of the accelerating Airy pulse (Equation (88)). On the left, the real and imaginary components of the pulse’s field are presented (for the final distance of ζ = β 2 z / θ 2 = 1 ), and on the right, the pulse’s intensity ( | A ( t , z ) | 2 ) is presented (for the final distance of ζ = β 2 z / θ 2 = 2 ).
Applsci 09 00527 g025
Figure 26. Same as Figure 25, but for Equation (90) with a = 0.1 / θ (for the final distance of ζ = β 2 z / θ 2 = 2 ).
Figure 26. Same as Figure 25, but for Equation (90) with a = 0.1 / θ (for the final distance of ζ = β 2 z / θ 2 = 2 ).
Applsci 09 00527 g026
Figure 27. Same as Figure 2, but for the intensity ( | A ( t , z ) | 2 ) of the pulse presented by Equation (92).
Figure 27. Same as Figure 2, but for the intensity ( | A ( t , z ) | 2 ) of the pulse presented by Equation (92).
Applsci 09 00527 g027
Figure 28. Same as the right figure in Figure 24, but with η / θ = 0.2 .
Figure 28. Same as the right figure in Figure 24, but with η / θ = 0.2 .
Applsci 09 00527 g028
Figure 29. Same as Figure 28, but for the pulse presented by Equation (97), with η / θ = 0.2 and a = 0.1 / θ .
Figure 29. Same as Figure 28, but for the pulse presented by Equation (97), with η / θ = 0.2 and a = 0.1 / θ .
Applsci 09 00527 g029
Figure 30. Same as Figure 2, but for the intensity ( | A ( t , z ) | 2 ) of the pulse presented by Equation (97) for a = 0.2 / θ and η / θ = 1 . The horizontal line corresponds for the maximum intensity distance z = η / β 2 a .
Figure 30. Same as Figure 2, but for the intensity ( | A ( t , z ) | 2 ) of the pulse presented by Equation (97) for a = 0.2 / θ and η / θ = 1 . The horizontal line corresponds for the maximum intensity distance z = η / β 2 a .
Applsci 09 00527 g030
Figure 31. FWHM Comparison between four different pulses: Gaussian (dashed curve), Rectangular (dotted curve), Exponential (solid curve), and Sinc (dot-dash).
Figure 31. FWHM Comparison between four different pulses: Gaussian (dashed curve), Rectangular (dotted curve), Exponential (solid curve), and Sinc (dot-dash).
Applsci 09 00527 g031

Share and Cite

MDPI and ACS Style

Granot, E. Analytical Solutions for the Propagation of UltraShort and UltraSharp Pulses in Dispersive Media. Appl. Sci. 2019, 9, 527. https://doi.org/10.3390/app9030527

AMA Style

Granot E. Analytical Solutions for the Propagation of UltraShort and UltraSharp Pulses in Dispersive Media. Applied Sciences. 2019; 9(3):527. https://doi.org/10.3390/app9030527

Chicago/Turabian Style

Granot, Er’el. 2019. "Analytical Solutions for the Propagation of UltraShort and UltraSharp Pulses in Dispersive Media" Applied Sciences 9, no. 3: 527. https://doi.org/10.3390/app9030527

APA Style

Granot, E. (2019). Analytical Solutions for the Propagation of UltraShort and UltraSharp Pulses in Dispersive Media. Applied Sciences, 9(3), 527. https://doi.org/10.3390/app9030527

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