`Mathematical Problems in EngineeringVolume 2010 (2010), Article ID 695208, 15 pageshttp://dx.doi.org/10.1155/2010/695208`
Research Article

## Mathematical Transform of Traveling-Wave Equations and Phase Aspects of Quantum Interaction

1Department of Electrical and Computer Engineering, University of West Florida, 11000 University Parkway, Pensacola, FL. 32514, USA
2Faculty of Applied Sciences, Politechnica University, Hagi-Ghita 81, 060032 Bucharest, Romania

Received 28 July 2009; Accepted 30 August 2009

Copyright © 2010 Ezzat G. Bakhoum and Cristian Toma. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The traveling wave equation is an essential tool in the study of vibrations and oscillating systems. This paper introduces an important extension to the Fourier/Laplace transform that is needed for the analysis of signals that are represented by traveling wave equations. Another objective of the paper is to present a mathematical technique for the simulation of the behavior of large systems of optical oscillators.

#### 1. Introduction

The Fourier and the Laplace transforms are great mathematical tools that have served science well for nearly two centuries. Recently, however, the Fourier/Laplace transform has become increasingly inadequate for treating one specific class of signals, namely, signals represented by traveling wave equations. The Fourier/Laplace transform is essentially a tool for time domain to frequency domain translation, and vice versa. In that transformation, the dependence of each component of the original signal on time is described by a time-related constant: the angular frequency . The transformation then maps the time domain into the domain, essentially by finding correlations between the original signal and an infinite set of pure sinusoidal signals. Those correlations collectively become a “frequency-domain" representation of the signal. In the class of signals represented by traveling-wave equations, however, the signal has existence in two domains: the time domain and the spatial domain . More specifically, the signal has dependence on two propagation constants: the temporal propagation constant , and the spatial propagation constant [1, 2]. Hence, a transformation is needed to map the behavior of the signal in the domain into the domain. The Fourier/Laplace transform only maps the domain into the domain and is incapable of performing the mapping that is required in this case. Here, it is important to point out that the transform known as the “multidimensional Fourier transform" [3], used extensively in image processing applications, is still a mapping process. By its definition, the “multidimensional Fourier transform" cannot perform a transformation of a signal that has simultaneous existence in both the time and the spatial domains.

The first objective of this paper is to present an important extension to the Fourier/Laplace transform that will allow the 2 2 mapping that is required in the case of traveling-wave signals to be performed. Another objective of the paper is to demonstrate that certain partial fraction decomposition and differentiation techniques can obviate the need for complex computations to simulate the behavior of large systems of optical oscillators. More particularly, some optical phenomena, such as the phenomenon of optical echo, are very hard to simulate without using physical variables corresponding to the internal states of a great number of atoms. Typically, such simulations require complex (order , where is the number of atoms) computations. We demonstrate that this problem can be fundamentally solved in a much simpler manner by using derivative-taking procedures for symbolic polynomials.

#### 2. A Dual-Domain Extension for the Fourier/Laplace Transform

Signals that physically propagate in various media are usually represented by traveling-wave equations, the most basic of which is the simple sinusoidal traveling wave [1, 2]:

where is the usual temporal propagation constant and where is the spatial propagation constant. The usual Fourier transform for a signal that exists only in the time domain is given by [4, 5]

We first define a companion to the Fourier transform for signals that are defined only in the spatial domain:

This transform maps the domain into the propagation constant () domain. We now introduce a dual-domain transform as follows:

Notice that the Fourier transform measures the correlations between the signal and an infinite set of pure sinusoidal signals in the time domain, while the transform of (2.4) measures the correlations between the signal and an infinite set of pure sinusoidal traveling waves. This is the main difference between the transform of (2.4) and the other known extensions to the Fourier transform.

To write the one-sided version of the transform in (2.4) (which will be equivalent to the Laplace transform), we first decompose the exponential term as follows:

where and . The one-sided transform will now be written as follows:

We can now clearly see that this dual-domain transform differs from the basic Laplace transform in that an extra exponential term is introduced, with being an imaginary number representation of the spatial propagation constant . The basic Laplace transform contains the imaginary number only, which is a representation of the temporal propagation constant .

#### 3. Dual-Domain Transforms of Some Well-Known Propagating Signals

We will now demonstrate how the dual-domain transform is calculated for various propagating signals. The transforms of two well-known signals, specifically, the propagating sinusoid and the propagating unit-step function are calculated here as examples.

##### 3.1. Transform of

By using the well-known decomposition

and carrying out the integration in (2.6) in a straightforward manner, the dual-domain transform of is found to be the following:

We can clearly see that if this transform reverts to the Laplace transform, with the presence of in the denominator. Essentially, what we get in that case will be the Laplace transform of , mapped into the domain.

##### 3.2. Transform of a Propagating Unit Step Function

A propagating unit step function is shown in Figure 1. The velocity of propagation of the wave is .

Figure 1: A propagating unit step function.

Here, due to the positive exponent of in (2.6), we will reverse the order of the limits on the integrals in order to obtain a meaningful result. Specifically, the time domain will be assumed to extend from to 0 and the domain will be assumed to extend from to . The transform of (2.6) will now be written as

The result is

We can again easily observe that if the propagation velocity , the Laplace transform of a stationary unit step function is obtained, with the presence of an extra factor of that represents the mapping into the domain.

#### 4. Dual-Domain Transform Theorems

We will now present two important theorems related to the dual-domain transform. Those theorems, namely, the phase-shift theorem and the convolution theorem, are straightforward extensions to the corresponding theorems of the Fourier/Laplace transforms.

##### 4.1. The Phase-Shift Theorem

The dual-domain transform of a signal that is phase shifted by a phase , where and are shifts along the and dimensions, respectively, is given by

where is the dual-domain transform of .

Proof. One has Let and , so that and . Hence

##### 4.2. The Convolution Theorem

A signal the dual-domain transform of which is given by , where and are the transforms of another signal and a transfer function, respectively, is given by

Proof. will be given by the inverse dual-domain transform of , that is, But Hence By interchanging the order of the double integrals, we have

#### 5. Example Application: Transfer Function of a Wireless Communication Channel with Multiple Phase Shifts

The transfer function characteristics of urban wireless communications channels are very important for cellular phone and wireless multimedia applications [68]. The high-frequency signals that travel through such channels are all represented by traveling wave equations, and each signal suffers from various phase shifts that collectively result in “fading effects" as the signal travels through the channel. Characterizing the behavior of a channel by means of a transfer function has always been a challenging problem, mainly because the Fourier/Laplace transform can only provide a transfer function representation in the frequency domain, that is, the domain, while the behavior of the channel in reality must ultimately be described in both the and the domains. We will now proceed to obtain such a dual-domain transfer function for a wireless communication channel by using the dual-domain transform, and important new conclusions about the behavior of such channels will be reached.

According to the phase-shift theorem (4.1), a signal that encounters multiple phase shifts in a wireless communication channel will be represented by the sum

at the output of the channel, where are the various phase shifts and where are the amplitudes of the phase-shifted components. The transfer function () of the channel is therefore given by

By expanding each exponential term as a power series and rearranging the equation, we get

For comparison, the single-domain Laplace transform of a rectangular pulse in the time domain is given by [5]

where are constants. It is clear that while (5.4) is a power series in , (5.3) is a power series in . This means that (5.3) is the transform of pulse-like characteristics in the time domain and pulse-like characteristics in the spatial domain. This is a very important conclusion, since it was previously assumed that multipath wireless channels exhibit pulse-like characteristics in the time domain only [6]. The pulse-like behavior of the channel in the spatial domain will be then described by a function in the domain, as shown in Figure 2.

Figure 2: Behavior of the multipath wireless channel in the domain.

The conclusion therefore is that the channel acts as a band-pass filter in both the frequency domain and the propagation constant domain. This conclusion has substantial implications for the analysis of multipath wireless channels and will be further developed and investigated in a subsequent report in a specialized communications journal.

#### 6. Aspects Regarding the Propagation of Traveling Waves Corresponding to Wave Functions. Phase Loss due to Multiple Quantum Interaction

We can extend our study at quantum phenomena by taking into account the fact that the expression corresponds also to the phase of a wave function describing the propagation of a free quantum particle.

According to the statistical interpretation of quantum theory, the wave function is representing just a mathematical model suitable for statistical interpretation. In the David Bohm interpretation of quantum theory, the wave function is a real wave in configuration space. This interpretation is also in accordance with measurements aspects of interaction of wave-trains with observer's material medium [9]. Due to these reasons, the phase of the traveling quantum wave requires a more detailed analysis, the dynamical aspect having to be taken into consideration.

Let us consider the propagation of an electron in a central symmetric electrostatic (Coulomb) field. According to the standard quantum analysis (based on second quantification theory, describing creation and annihilation of particles) the central symmetric electrostatic field should be decomposed (using the Fourier transformation) in a set of waves with a certain angular frequency and a certain wave vector :

where the product stands for

with representing the time coordinate, representing to the vector of position, representing the angular frequency ,and representing the wave vector. The measurement system is chosen so as for performing a better correspondence from the angular frequency to the energy and from the wave vector to the momentum of the quantum particle (see [10] for more details).

The quantities and do not correspond to standard photons, and for this reason they cannot be substituted by operators as required by second quantification theory. However, experimental facts have shown that the electromagnetic field effect can be studied using the perturbations method from quantum theory.

Using this method [10] the matrix element (for the first-order of perturbation) corresponding to the electromagnetic interaction between the electron and the electromagnetic field is represented by,

which can be also written as,

In both previous equations the expression corresponds to the sum:

where represents the cuadrivector of the electromagnetic field and represents the Dirac matrices (its argument corresponds to the energy and momentum conservation laws).

If the electromagnetic field does not depend on time, the cuadripotential can be presented as

which shows that the virtual photons composing the electromagnetic field have a nonzero value just for the momentum (the energy corresponding to quantity being equal to zero). As a consequence, the matrix element for the first-order of perturbation can be written as

By performing the integration on , it results

For the coulombian electrostatic field with central symmetry ; it results

where represents the electrical negative charge. The previous equation can be also written (in representation) as

As a consequence, the matrix element for the first-order of perturbation in case of coulombian field with central symmetry can be written as:

which can be also written as

According to standard interpretation of quantum theory, this first-order element from perturbation method is connected to the probability of an interaction between an electron with initial momentum and energy and a virtual photon with momentum so as to result an electron with momentum and energy . This interaction can be also represented using Feynman diagrams. By analyzing higher orders of matrix element (presented in an intutive manner using the same Feynman diagrams) some supplementary phenomena can be also explained, as elastic diffusion of photons by electrons, the annihilation of an electron-positron pair in two photons, the generation of an electron-positron pair by photons, the emission of two photons by an electron, the emission of brehmstrahlung radiation by electrons and positrons and so on.

However, for describing the whole interaction of an electron in a central symmetrical electrostatic field (so as to determine the trajectory of the associated wave-train) we can consider that the electron has to undergo multiple interactions with such virtual photons until the action of the exterior field vanishes (such multiple interactions being allowed by quantum mechanics). According to the quantum laws, each interaction transform the initial electron (corresponding to a certain wave train) into a new wave-train (the final electron) with different characteristics (another momentum and energy), and so on.

However, we should take also into account the fact that the wave-train associated to a quantum particle has not just a statistical interpretation. The Bohm-Aharonov effect has shown that a certain phase should also be associated to such a wave-train. Considering that the wave-train corresponding to an electron undergoes a set of interactions with virtual photons, it results that a certain transient time is required by each interaction, so as physical quantities as wavelength and angular fequency to be defined on a certain time interval (these quantities being used by the formula of matrix elements of interaction). This transient time causes a phase-loss for the initial wave-train, while the timelength of stabilized oscillations (corresponding to the new-generated wave-train associated to the final photon) is represented (in a classical dynamical model) by the difference between the timelength of the initial wave-train and the transient time. Considering that no supplementary stabilized oscillations are added by each interaction (this means that quantum interaction does not generate any output wave-train after the input wave-train has vanished) it results that a very great number of interactions of the initial electron in a central symmetrical electrostatic (coulombian) field would cause a great number of phase loss (each interaction substracting a certain transient time), and, finally, the timelength of the wave-train associated to the final electron tends to zero. This corresponds to a vanishing phenomena for the initial electron in contradiction with experimental facts.

#### 7. Nonlinear Phenomena of Phase-Loss Compensation for Propagating Wave-Trains in Quantum Physics: The Need of a Certain Material Reference System

At an extended analysis of propagating phenomena in quantum physics, we can notice that an approach based on second quantification laws leads to the same aspects of phase-loss (or shortened timelength) of a propagating wave-train due to multiple interactions along its trajectory in different situations, such as the interaction of a propagating free electron with its own electromagnetic field (implying renormalization of electron mass using higher orders of perturbation theory) and the generation of virtual electron-positron pairs due to the interaction of an electron with a weak exterior electromagnetic field (analyzed by taking into account higher order of perturbation theory, implying renormalization of electron charge) [10].

For justifying the fact that the phase-loss does not appear in experimental phenomena of such propagating beams (as presented above), we can take into consideration a possible phenomenon of phase-addition, consisting in an addition of supplementary oscillations after the initial wave-train has vanished. According to mathematical dynamical models, this corresponds to a dead-time transient phenomenon. As a consequence, a kind of internal memory should be associated to the space interval where the interaction occurs, generating an output signal even after the input signal has disappeared.

Another possibility consists in taking into account (from the same dynamical point of view) a resonance phenomena, considering that the initial wave-train can generate (as in a resonant phenomenon) oscillations with amplitude above a certain threshold. When the initial wave-train (similar to an input signal) disappears, these oscillations will become damped oscillations able to generate a final wave-train (similar to an output signal) as long as their amplitude is above a certain threshold. As in the case of the previous dynamical model, this mathematical model requires also an internal memory of the threshold level and of certain state variables able to generate resonant oscillations.

Both possibilities previously presented have a common feature: a certain delay time appears at each interaction, and thus the overall delay time should be connected with the wave-train velocity.

However, the phase-loss phenomenon can be avoided (as a theoretical result) if the whole trajectory is considered to be determined by just one Lagrangean function for the whole interaction of the electron with the electromagnetic field, implying the dynamical equations:

[10]. Thus the wave-function corresponding to the electron and the electromagnetic field is modified in a continuous manner along its trajectory. However, this possible solution is based on quantum laws of first quantification theory, all phenomena corresponding to second quantification (involving creation and annihilation of particles) being neglected. Yet we have to take into account the fact that the exterior electromagnetic field has usually a higher magnitude and a slow time variation as compared to the energy and velocity of the interacting electron. So it is quite naturally to consider that a certain exterior field (more important at a macrosopic scale) is determining the evolution of the electron wave-train, and thus appears the need of a certain reference system which acts upon the wave-trains corresponding to quantum particles. This assumption is supported also by the lack of reversibility for diffraction phenomena in quantum physics. We can notice an electron diffraction phenomenon when an electron beam interacts with a motionless crystal lattice, but we cannot imagine a diffraction phenomenon for the atoms of a crystal lattice when a moving lattice interacts with a motionless spatial distribution of electrons.

#### 8. Phase Aspects and Quantum Dynamics of Cooperative Phenomena: Formal Logic Problems for Cooperative Quantum Phenomena

Our study upon quantum dynamics and phase aspects should be extended by taking into consideration multiparticle cooperative phenomena in quantum physics involving emission of electromagnetic field, while key issues as transition from a set of single particle quantum states to a global cooperative state can be noticed. For example, the photonic echo phenomenon is based on the rephasing that can occur in an inhomogeneously broadened atomic system, when the motion of the Bloch vectors is reversed following the initial dephasing of the atomic dipoles. For the beginning, we consider a group of atoms whose Bloch vectors in the rotating frame have the direction along the axis at the zero moment of time. We consider also that the rotating frame has the angular frequency (the central frequency of the inhomogeneous line). In time, the Bloch vectors of the atoms will precess clockwise or counter-clockwise around the axis, depending whether their natural angular frequencies are either greater or less than . In time a dephasing phenomenon occurs, the atomic Bloch vectors being distributed almost uniformly around the plane. The initial macroscopic moment becomes zero, and the system does not radiate.

By applying a short, intense pulse at time moment a reversal of the coordinate appears, the motion of the atomic Bloch vectors is reversed, and the vectors would be aligned again at a certain time moment . As a consequence, the atomic system generates a high-intensity echo pulse at time moment . This is the basic principle of both spin echo phenomenon and photon echo phenomenon [11].

Usually, a preparatory pulse (from the superradiance states point of view) is applied to the active medium at the initial moment of time. Mathematically, the property of the reference pulse of being a pulse is represented by the relation

where stands for the amplitude of the electric field, stands for its timelength and represents a parameter depending on frequency. The condition for the main optical pulse of being a pulse is represented by the relation

where stands for the electric field of this signal, stands for its timelength, and represents the same parameter depending on frequency.

A basic feature of the high-intensity pulse (the photonic echo pulse) is represented by its power. Unlike amplified optical pulses by stimulated emission (with the power of emergent pulse proportional to ), the power of the photonic echo pulse is proportional to where represents the number of atoms involved in the emission of this high-intensity pulse. This feature can be derived using the formalism of superradiant states , where stands for the cooperative quantum number and is proportional to the unperturbed Hamiltonian function for the ensemble of atoms (due to the the similarities with quantic angular momentum theory, ). According to this similarities, when the state of the atomic system is described by the state , the intensity of the emitted pulse is

If all atoms involved in light emission are in excited state, and the intensity of emitted pulse is equal to

where stands for the intensity of an optical pulse emitted by a single atom ( ). It can be shown that the maximum for the emission rate is obtained for and greater values of ; while is also equal to ( representing the number of atoms in excited state and representing the number of atoms in ground state) it results that or and . In this case the intensity of emitted pulse is

being proportional (for large ) with .

These aspects can be also explained (in an intuitive manner) using classical electromagnetics. When all atoms radiate independently, it can be considered that the phase and angular frequency of emitted electric field are not the same and thus the average value of the intensity corresponds to a sum of terms (all terms having the form vanish while they correspond to alternating functions). It can be also argued that different orientation of Bloch vectors for a great number of atoms imply a null result for the sum of these vectors from the very beginning. On the contrary, when all atoms radiate in a cooperative manner, it can be considered that their emitted electric field has the same phase and angular frequency and thus the electric field generated by all atoms is represented by a sum and the intensity of the emitted electric field is represented by . For it results that the intensity of the emitted electric field is proportional to (in the first case) and to (in the second case). It can be also argued that the situation when all Bloch vectors are aligned corresponds to a great global dipole momentum for the entire system of atoms, and thus the emitted radiation has the same phase.

However, an important question regarding phase appears for these cases: which is the angular frequency of the emitted field for the photonic echo phenomenon? Is there a certain angular frequency or should we consider the emitted field as a mixture of waves with different characteristics? If we consider that we deal with a cooperative phenomenon, it seems logical to search for a certain angular frequency and a certain phase for the whole emitted field (as in the case of correlated photon emission, when two photons act as a single entity even when they are situated at large distances). Thus the emission would be similar to a single photon emission and interference phenomena could be noticed further (according to the Dirac principle that a photon interacts just with itself). Yet this affirmative answer generates another question: under which circumstances an ensemble of atoms can be studied using the quantum formalism of superradiant states, instead of analyzing a superposition of atoms with independent emission? According to formal logic in physics, an interference phenomenon can either occur or not. However, it is quite difficult to imagine that a certain parameter of the ensemble could suddenly generate the transition to a single quantum state when it reaches a specific value (so as interference phenomena to occur). This aspect is similar to the problem of defining a certain material reference system which performs the Lorentz transformation upon a received wave-train [9]: which are the requirements for a certain material reference system (as related to the parameters of the received wave-train)? (see [12] for the theoretical case when the energy of the received wave is comparable to the rest energy of the atoms composing the material medium). More probably, some correlation techniques with past events are required, and thus phase correlation aspects could appear. The phase correlation represents also the main technique for creating and controlling cooperative phenomena (when some spatially distributed particles or atoms can be described by a single quantum function at a certain moment of time, and possible transitions are controlled by imposing phase coincidence).

#### 9. Simulating High-Intensity Superradiant Pulses Using Derivative Procedures and Partial Fraction Decomposition

A step forward in looking for an adequate mathematical model should consist in finding a mathematical quantity proportional to ; this step must be performed in a similar manner used for passing from coordinates representation to momentum representation in quantum theory, by the aid of Fourrier transformation (the basic variables being changed) [13]. We can easily notice that the derivatives of different order for are direct proportional to according to

It can be noticed that the coefficient of final monomial equals which represents the ratio of the intensity of the photonic echo pulse emitted from the state to the intensity of the pulse emitted independently by a single atom (similar to derivative procedures presented in [14]). In an intuitive manner, can be considered as an operator of creation (as in second quantification formula) and the exponent symbolizes the creation of the quantum superradiant state with quantum number as a result of successive application of operator upon an initial preliminary state (with ). The same derivative procedure could be used for obtaining the eigenvalue of the square of the angular momentum according to

and considering an eigenvalue the coefficient of the resulting monomial.

For the case of nonsuperradiant states the ensemble of atoms could be described using a function as

where stands for the specific phase difference for each atom. While the power of equals unity, a second-order differentiation yields zero, and no superradiant pulse appears.

#### 10. Conclusions

This paper has presented an important extension to the Fourier/Laplace transform that is needed for the analysis of signals that are represented by traveling wave equations. Specific mathematical technique for the simulation of the behavior of large systems of optical oscillators was obtained. Signals described by traveling wave equations are typically characterized by a spatial propagation constant, , in addition to the usual temporal propagation constant . The Fourier/Laplace transform can only perform mapping of the time domain into the domain (or frequency domain), and is incapable of performing the more general mapping of the domain into the domain. The new extension introduced here allows for that general mapping to be performed. This study has also emphasized the importance of phase dynamics aspects for traveling wave-trains in multiple interaction phenomena in quantum physics. It was shown that a phase-loss phenomenon should appear due to transient time, implying an annihilation of the quantum particle due to multiple interactions. The absence of this phenomenon in experimental results is explained by a possible correlation between the trajectory and phase of the quantum wave-train before interaction and the trajectory and phase of the quantum wave-train after the corresponding quantum particle has undergone the interaction. The need of a certain exterior reference system acting upon the quantum wave-train is also presented (based on some diffraction aspects). Finally it is shown that some interference aspects of the photonic echo phenomenon imply the need of considering that superradiant quantum states for an ensemble of atoms can either exist or not. Using a specific differentiation technique for symbolic polynomials corresponding to phase difference of each atom, a useful formalism for determining the condition for generating a superradiant pulse and its intensity is obtained.

#### References

1. W. H. Hayt and J. A. Buck, Engineering Electromagnetics, McGraw Hill, New York, NY, USA, 2006.
2. R. Feynmam, R. Leighton, and M. Sands, The Feynman Lectures on Physics, vol. 2, Addison Wesley, Reading, Mass, USA, 1977.
3. J. S. Lim, Two-Dimensional Signal and Image Processing, Prentice Hall, Upper Saddle River, NJ, USA, 1990.
4. T. W. Korner, Fourier Analysis, Cambridge University Press, Cambridge, UK, 2nd edition, 1989.
5. J. E. Marsden and M. J. Hoffman, Basic Complex Analysis, W. H. Freeman, New York, NY, USA, 2nd edition, 1987.
6. T. S. Rappaport, Wireless Communications, Prentice Hall, Upper Saddle River, NJ, USA, 1996.
7. V. K. Garg and J. E. Wilkes, Wireless and Personal Communications Systems, Prentice Hall, Upper Saddle River, NJ, USA, 1996.
8. V. K. Garg, Wireless Network Evolution, Prentice Hall, Upper Saddle River, NJ, USA, 2002.
9. E. Bakhoum and C. Toma, “Relativistic short range phenomena and space-time aspects of pulse measurements,” Mathematical Problems in Engineering, vol. 2008, Article ID 410156, 20 pages, 2008.
10. N. Nelipa, Physique des Particules Elementaires, MIR, Moscou, Russia, 1981.
11. L. Mandel and E. Wolf, Optical Coherence and Quantum Optics, Cambridge University Press, Cambridge, UK, 1995.
12. A. Toma, S. Pusca, and C. Morarescu, “Spatial aspects of interaction between high-energy pulses and waves considered as suddenly emerging phenomena,” in Proceedings of the 1st International Conference on Computational Science and Its Applications (ICCSA '06), vol. 3980 of Lecture Notes Computer Science, pp. 839–846, Glasgow, UK, May 2006.
13. P. Sterian and C. Toma, “Methods for presenting key concepts in physics for MS students by photon-MD program,” Bulgarian Journal of Physics, vol. 27, no. 4, pp. 27–30, 2000.
14. Th. Toma, S. Pusca, and C. Morarescu, “Simulating superradiant laser pulses using partial fraction decomposition and derivative procedures,” in Proceedings of the 1st International Conference on Computational Science and Its Applications (ICCSA '06), vol. 3980 of Lecture Notes Computer Science, pp. 779–784, Glasgow, UK, May 2006.