[Paper Review] Propagation of the main signal in a dispersive Lorentz medium

A. Ciarkowski|ArXiv.org|Dec 7, 2004

Orbital Angular Momentum in Optics9 references3 citations

TL;DR

This paper derives uniform and non-uniform asymptotic representations for the main signal in a dispersive Lorentz medium excited by a sine-modulated pulse with a hyperbolic tangent envelope that turns on abruptly at finite time. Using the Bleistein-Handelsman method for integrals with nearby saddle points and algebraic singularities, it shows how the uniform representation reduces to the non-uniform one via residues, with numerical validation under varying pulse rise times and medium parameters.

ABSTRACT

Evolution of the main signal in a Lorentz dispersive medium is considered. The signal propagating in the medium is excited by a sine-modulated pulse signal, with its envelope described by a hyperbolic tangent function. Both uniform and non-uniform asymptotic representations for the signal are found. It is shown when the uniform representation can be reduced to the non-uniform one. The results obtained are illustrated with a numerical example.

Motivation & Objective

• To analyze the evolution of the main signal in a dispersive Lorentz medium when excited by a pulse with finite rise time.
• To derive both uniform and non-uniform asymptotic representations for the signal's propagation dynamics.
• To clarify the mathematical conditions under which the uniform asymptotic form reduces to the non-uniform form obtainable via residue calculus.
• To validate the theoretical results with a numerical example using realistic medium parameters and varying pulse rise rates.

Proposed method

• Formulates the electromagnetic signal propagation using a complex phase function derived from the Lorentz model of dielectric response.
• Applies the Laplace transform to the hyperbolic tangent envelope to express the initial signal in the frequency domain.
• Employs the Bleistein-Handelsman method for uniform asymptotic evaluation of integrals with a nearby saddle point and algebraic singularity.
• Deforms the integration contour to the Olver-type path passing through near and distant saddle points to isolate pole contributions.
• Derives the uniform asymptotic representation in terms of the complementary error function, capturing behavior across critical transition points.
• Compares the uniform form with the non-uniform residue-based representation, showing their equivalence under specific parameter regimes.

Experimental results

Research questions

• RQ1 Under what conditions does the uniform asymptotic representation of the main signal reduce to the non-uniform representation derived via residues?
• RQ2 How does the rise time of the excitation pulse, parameterized by β, affect the signal evolution and precursor structure in a Lorentz medium?
• RQ3 What role do the poles of the Laplace-transformed envelope play in the asymptotic behavior of the main signal?
• RQ4 How does the merging of saddle points at θ = θ₁ affect the validity of the asymptotic approximations?
• RQ5 In what way does the choice of medium parameters (b, ω₀, δ) influence the signal's amplitude and temporal evolution?

Key findings

• The uniform asymptotic representation based on the complementary error function remains valid even when the second derivative of the phase function vanishes at the saddle point, unlike the non-uniform form.
• The uniform representation reduces to the non-uniform residue-based form when the saddle point and pole are well-separated, confirming consistency between methods.
• Numerical results show that for large β (fast rise time), the main signal evolves with a sharp peak, while for small β (slow rise time), the signal growth is delayed and attenuated due to the influence of the second pole.
• At θ ≈ 1.5, the uniform representation breaks down due to φωω → 0, indicating a critical transition point where the standard asymptotic approximation fails.
• The inclusion of the second pole (ωc2 = ωc − 2iβ) in the asymptotic expansion leads to a slower signal growth, as confirmed by Fig. 4 in the numerical example.
• The numerical example using Brillouin’s medium parameters (b = √20×10¹⁶ s⁻¹, ω₀ = 4.0×10¹⁶ s⁻¹, δ = 0.28×10¹⁶ s⁻¹) confirms the theoretical predictions across varying β values.