.
Our site needs your help!
Site categories

Feynman Path Integrals in the Problems of Wave Propagation in Random Media

Join Our Telegram (Seaman Community)

In this section we introduce an alternative representation for a Green function based on Feynman path integrals. This approach is most suitable for the analysis of the wave propagation and scattering in a line-of-sight region. This section provides some basics for study in an unbounded medium, while later in Section 3*** we introduce a path integral representation in the presence of a boundary.

The basic Feynman postulate states that the probability amplitude for the transition of a particle-wave from point
r0
to point
r
can be represented by contributions made by individual trajectories over which the particle-wave can propagate between points
r0
and
r
. The contribution of each trajectory is proportional to exp {jS/ħ}, where S is a classical action and ħ is Planck’s constant. As is well known, the quantum-mechanical probability amplitude of transition is equivalent to the Green function
Gr, r0
of the Helmholtz equation written in a parabolic approximation, after the time has been replaced by the x-coordinate along the selected direction of propagation and m/ħ has been replaced by the wavenumber k = 2π/λ, where m is a particle’s mass, k is the wave number of the admitted wave, and λ is the wavelength. Turning to the trajectory continuum, we can write:
Gr, r0=Dρx exp jkSρx,            Eq. 1

where:

  • r=x, ρ, r0=0, ρ0
    are the coordinates of the observation point and source, respectively;
  • ρ=γ, z, ρ=γ0, z0
    are the vectors in the plane orthogonal to the direction of propagation x, and;
  • Dρx
    is the differential in the space of continuous trajectories.

The action S is given by:

Sρx=0xdxL0x, ρx+12δεx, ρx,               Eq. 2

where:

  • L0x, ρx=12dρdx2+U0x, ρx
    is an unperturbed (when δε = 0) Lagrangian in the small-angle approximation;
  • U0x, px
    is an unperturbed potential;
  • U0x, px=1/2εx, ρx1
    , the angle brackets < … > here and later denote averaging over the ensemble of dielectric permittivity;
  • εx, ρ, δεx, ρ
    is a random component of
    εx, ρ, δε=0.
Consider calculation of the first two moments of the Green function (1) in the case of an unbounded stationary medium filled with random inhomogeneities of dielectric permeability
εx, ρ
. In continuous notation the average Green function
Gr, r0
and coherence function:
Γx, ρ1, ρ2, ρ0, ρ0=Gr1, r0 G*r2, r0
Gr, r0=Dρx exp jkS0ρx expjk20xdxδεx, ρx,             Eq. 3
Γx, ρ1, ρ2, ρ0, ρ0=Gr1, r0G*r2, r0=Dρx  DRx exp jk0xdxdRdxdρdx·exp jk20xdxδεx, Rx+ρx2δεx, Rxρx2.           Eq. 4
In Eq. (4) we have introduced the sum
Rx=1/2ρ1+ρ2x
and the difference
ρx=ρ1xρ2x
of the trajectories with the following boundary conditions:
Rx=0=12ρ0+ρ0 ρx=0=ρ0ρ0,
Rx=x=12ρ1+ρ2 ρx=x=ρ1ρ2.

We can assume that the fluctuations in δε are statistically uniform:

δεx1, ρ1δεx2, ρ2=Bεx1x2, ρ1ρ2

and have different correlation scales:

  • Lx in direction x and;
  • L in the plane (γ, z).
If x ≫ L the integral
0xdxδεx, ρx
represents a Gaussian random value in accordance with the central-limit theorem. Hence:
exp jk20xdxδεx, ρx=expγs,           Eq. 5
exp jk20xdxδεx, Rx+ρx2δεx, Rxρx2=expDs,            Eq. 6

where:

γS=k280xdx0xdxδεx, ρxδεx, ρx              Eq. 7
is the variance of the phase fluctuations along the trajectory
ρx
and:
DS=k280xdx0xdxδεx, Rx+ρx2δεx, Rxρx2·δεx, Rx+ρx2δεx, Rxρx2                 Eq. 8
is a structure function of the phase difference along the trajectories
Rx+ρx2
and
Rxρx2.

Introduce a two-dimensional spatial spectrum of the fluctuations de in a plane (γ, z):

Fεxx, κ=14π2d2ρBεxx, ρejκρ.          Eq. 9

Then for γS and Ds we obtain:

γS=k240x/2dη2η2ηdξd2κFεξ, κ exp jκξdρdη,                Eq. 10
DSR, ρ=k240x/2dη2η2ηdξdξd2κFεξ, κ exp jκξdRdη·cos κξdρdηcos κ, ρ.                Eq. 11
To analyse Eqs. (10) and (11) let us introduce the parameter
θR=maxdRdx,
the characteristic angle of the ray’s trajectory along the x-axis. The value of θR is contributed to by the three characteristic parameters:

Therefore θR = max{θd, θF, θs}.

Let us also introduce α parameter of anisotropy α in the fluctuations of δε: α = L / Lx. The effective width of function
Fεξ, κ
over ξ does not exceed a correlation scale Lx. When inequality:
θRα                Eq. 12
holds, we can assume that
κξdρ/dη, κξdR/dηdη1
in a significant region over ξ in Eqs. (10) and (11). As a result, we obtain the following approximations:
γSx=πk2x4d2κΦε0, κ;             Eq. 13

and:

DSR, ρDSρ=πk220xdxd2κΦε0, κ1cos κρx                Eq. 14

where:

2πΦε0, κ=dξFεξ, κ                Eq. 15
and
Φεκ
is a three-dimensional spatial spectrum of fluctuations in δε. We can reasonably assume that θR ~ θs, and the inequality (12) takes the form:
kL2Lx1                Eq. 16
which means that the longitudinal correlation scale has to be small compared with the distance
kL2
where the diffraction on the irregularities of the scale L become significant. The expressions (5) and (14) provide the fundamental solution to the equation of the coherence function Γ in the Markov approximation:
Γxjk2ΓRρπk24HρΓ=0,             Eq. 17
Hρ=2d2κΦε0, κ1cos κρx

and:

DSρ=πk2/40xdxHρx.

Local” conditions of the Markov approximation are bounded by inequalities (3) and (4), for an average field we also need smallness of the attenuation of the average field over the distance of LxS(Lx) ≪ 1. “Non-local“ conditions lead to the inequality:

Rcλ,              Eq. 18

where Rc is a coherence radius of the scattered field in a plane (γ, z), which depends on the distance x, Rc ≡ Rc(x). The parameter Rc can be found from:

DSRc.              Eq. 19
To illustrate the path integral approach, we obtain a known solution to Eq. (17) for the unbounded medium filled with statistically uniform fluctuations of
δεx, ρ.
Using Eqs. (4), (16) and (14) we obtain:
Γx, ρ1, ρ2, ρ0, ρ0=DρxDRx·expjk0xdxdRdxdρdxDSρx.            Eq. 20
The integral over trajectories R(x) represents a continuous Fourier transform of the delta-function
δd2ρ/dx2.
A Lagrangian in the exponent can be expanded into a functional series up to the second order terms in the vicinity of the trajectory
ρ0x=ρx/x+ρ01x/x.
Integrating then over
ρx
, we obtain the well-known solution for an unbounded medium:
Γx, ρ1, ρ1, ρ0, ρ0=k24π2x2expjkRR0ρρ0xπk240xdxHρxx+ρ01xx.              Eq. 21
The relationship between the average of the path integrals (3), (4) and the average over the Fermat paths was analysed in Refs. The asymptotic evaluation of the path integral representation for the Green function (3) may be performed using the method of stationary phase. The asymptotic method consists of finding the extreme path
ρ *x
which renders the minimum value of the phase
Sρx
given by Eq. (2). It can be noted that in obtaining the solution for the coherence function in an unbounded medium (21) we actually did not face any difficulties with integration over R(x) due to an isotropy of the medium,
εr=const.
However, this is not the case for more complex media. The variational problem for determining
ρ*x
gives rise to the Euler equation for the Fermat paths:
d2ρdx2+δεx, ρx=0,                 Eq. 22
where
=γ, z.

The rigorous mathematical solution to the problem encountered difficulties in a general case of the multi-scale inhomogeneities of δε and, therefore, the basic physics has to be involved for a qualitative analysis.

As mentioned above, there are two characteristic transverse scales of the phase fluctuations in
Sρ.
The first is the Fresnel zone size resulting from the first term in Eq. (3). The second term is due to a random component of the dielectric permittivity de and is of the order of the coherence scale Rc introduced by Eq. (19).

We assume a small value of the mean-square variation of the fluctuations in a phase difference of the fields approaching any points separated in space by a distance of the order of the Fresnel zone size. Hence, this assumption is rendered by the inequality:

δS1=k2π40xdxHxκ1.             Eq. 23
As known, the meaning of this inequality is that the Fresnel zone size is the characteristic region in the plane x = constant, from which the rays arrive at the receiving point in phase, even in the presence of the δε fluctuations. It also means that the integration in Eq. (3) can be fulfilled along a single canal-ray tube bounded in lateral cross-section by the Fresnel volume. We assume the extreme trajectory
ρ*x
to deviate slightly from the unperturbed path
ρ0x
, defined by the unperturbed Lagrangian
L0ρx,
we should require the fluctuations in the arrival angle of the wave defined by
ρ*x
to be small compared with angular size of the Fresnel zone ~(kx)-1/2. In a functional representation this requirement is equivalent to:
δS2=k2π80xdxδρ2x 2Hρxk2π8x2d2κ1.            Eq. 24
In Eq. (7) we evaluated the deviation
δρ=ρ*ρ0
to be of the order of the Fresnel zone size.
Having satisfied inequalities (22) and (23), the integration in Eq. (3) can be performed along the unperturbed ray trajectories
ρ0x
which are the Fermat paths for a free wave particle propagation (in the particular case of an unbounded medium):
ρ0x=ρρ0xx.                 Eq. 25

Such a straight-line approximation of trajectories in a path integral is similar to the extended Huygens–Fresnel principle. The applicability of extended Huygens–Fresnel principle shown that this approximation besides yielding the exact solution for the first two moments, provides a qualitatively correct solution for the high-order moments of the wave field.

Numerical Methods of Parabolic Equations

Among many numerical methods used in the problems of applied electromagnetics and wave propagation one may distinguish two methods most widely used in a VHF/UHF propagation in the atmospheric boundary layer, namely: split-step-Fourier and split-step Pade methods. Both methods are based on the parabolic approximation to a wave equation and differ in the method of obtaining the approximation of the exponential operator. Initially, the split-step parabolic approximation was introduced in a problem of applied acoustics where the wave propagation can be described by a scalar wave equation of elliptical type:

2ux2+2uz2+k2n2u=0                  Eq. 26

where u is the spectral amplitude of the wave component with frequency ω, k = ω/c is the wavenumber, c is the group velocity of the propagation and n is the refractive index of the medium and may vary with the coordinates, we assume, for now, a constant value of n, n = 1. Equation (26) is written for a two-dimensional case of the propagation since this case is the one most widely considered in applications of the split-step approximation.

Equation (26) can be factorised as follows:

x+jk2+2z2·xjk2+2z2 u=0,                 Eq. 27

representing the forward and backward propagating waves respectively. We may also remove fast phase variations with distance by introducing the slow-varying amplitude ũ, u = ũ exp (jkx). The truncated and factorised equation for amplitude ũ takes the form:

x+jk1+1+1k22z2·x+jk11+1k22z2 ũ=0.                 Eq. 28

By retaining the only forward propagating wave we obtain the forward parabolic equation:

x+jk11+1k22z2 ũ=0.                 Eq. 29

which has a formal solution:

ũx+x=expjkx1+1k22z21ũx                 Eq. 30

where Δx is a range increment. As apparent from Eq. (30) the forward propagating wave at range x + Δx can be obtained from values of the wave amplitude at the previous distance x by applying the exponential operator in Eq. (30). In order to actually calculate the wave field, the exponential operator in Eq. (30) should be approximated in a form suitable for computations. Let us introduce a notation Z = 1/k22/∂z2 and, depending on the type of approximation, we obtain three known types of the split-step approximation to parabolic equations:

a The standard or narrow angle approximation, obtained by expanding the square root in the exponent into a Taylor series and retaining the linear term:

1+Z=1+Z2.                 Eq. 31

This approximation is used in a split-step Fourier method to be discussed later, and is normally valid for narrow angles of propagation relative to axis x, not exceeding 15°–20°.

b The Claerbout approximation in the form:

1+Z1+34Z1+Z4                 Eq. 32

which is reported to be valid for wider angles up to 30–40°. And finally:

c The split-step Pade approximation, when exponential is expanded into series:

expjk1+Z11+m=1MαmZ1+bmZ              Eq. 33

where the coefficients αm, bm are determined numerically in the complex plane. The split-step Pade approximation is reported to be valid for angles of propagation relative to the x-axis of up to 90°.

All the above methods have been realized in a very powerful computational tech- nique well suited for parallel computing and widely used in numerous applications. With regard to the problem of radio wave propagation in the earth’s troposphere the most established computational approach is based on the split-step Fourier method equivalent to a standard narrow angle approximation of the square root in the expo- nential operator. Whilst the wide angle approximations, such as Claerbout and Pade approximations have been used in the problems of scattering on objects submerged in a free-space, the author does not know of any systematic derivation of the PadØ approximation in the case of radio wave propagation in a stratified troposphere over terrain or the sea boundary surface. We may notice that while the fundamental elliptic equation of type (22) can be obtained for potentials (Hertz functions) under free-space propagation conditions, the equations for the Debye’s potentials take a different form in the case of a stratified troposphere. In the presence of a random component of the refractive index the main equation for the slow varying envelope is given by Understanding Parabolic Approximation in Wave Propagation: Analytical Methods and Applicationsthis equation.

1r22W1θ2
on the right-hand side of the equation. While this is not very consistent with the arguments used in Understanding Parabolic Approximation in Wave Propagation: Analytical Methods and Applications“Parabolic Approximation to the Wave Propagation”, we retain this term in order to get in line with the approach of other authors to obtaining the parabolic equation. We may notice that to strictly follow the procedure we have to leave the other terms in the right-hand side of Understanding Parabolic Approximation in Wave Propagation: Analytical Methods and Applicationsthis equation in order to be consistent in the accuracy of the approximation. Introducing the same coordinates x = αθ, z = r – α and γ = ϕ · α sin ϑ we obtain:
2W1x2+j2kW1x+W1+k2εmz1+δεrW1=0.           Eq. 34

Then we may introduce the operator:

L=1k2+k2εmz1+δεr           Eq. 35

and factorize Eq. (34) in a form similar to Eq. (25):

x+jk1+1+L·x+jk11+LW1=0.

Then, following the above described procedure, we may retain the forward propagating wave obeying the truncated equation:

W1x=jk1+L1W1           Eq. 36

with formal solution at the marching step in the form:

W1x+x=W1x exp jkx1+L1.          Eq. 37

The exponential can then be expanded into a series similar to Eq. (33) and the next range step solution is given by:

W1x+x, γ, z=W1x+m=1MαmL1+bmLW1x, γ, z.          Eq. 38

We may introduce the auxiliary function:

fmx+x, γ, z=αmL1+bmLW1x, γ, z.          Eq. 39

which can be found by multiplying both parts of Eq. (40) by 1+bmL and solving a second-order differential equation for fm. Substituting fm into Eq. (39) we have:

W1x+x, γ, z=W1x+m=1Mfmx+x, γ, z.          Eq. 40
It is important to notice that all the functions fm can be solved independently and in parallel at each step along the distance x. The marching equation (40), as well as the forward propagating wave equation (30), shall be appended by both boundary conditions (at z = 0 and
r
) and initial conditions at x = 0.

A somewhat modified approach to the implementation of the Pade method is to substitute a Pade series with a Pade product:

m=1MαmL1+bmL=m=1M1+λmL1+μmL.                Eq. 41

Now, we concentrate on the split-step Fourier method. Consider Understanding Parabolic Approximation in Wave Propagation: Analytical Methods and Applicationsthis equation and assume the presence of deterministic and stratified refractivity in the medium, the random component δε is set to 0. Assume first that the propagation takes place in an unbounded medium, i. e., we have only requirements on the appropriate decay of the field at infinity. The next assumption normally employed in the split-step method applied to radio wave propagation in the troposphere is to assume the medium to be stratified over the z-coordinate, i. e. a stratified troposphere, thus removing the dependence on the y-coordinate and truncating Understanding Parabolic Approximation in Wave Propagation: Analytical Methods and Applicationsthis equation to a two-dimensional (x, z) parabolic equation.

Understanding Parabolic Approximation in Wave Propagation: Analytical Methods and ApplicationsEquation in this case can be presented in the form:

Wx=jk2εmz1+12k2z W.                Eq. 42

In a vertically stratified medium the formal solution to Eq. (42) is then given by:

Wx+x=expjxk2εmz1+12k2z2 Wx                Eq. 43

where Δx is a range increment. Equation (43) forms the basis for a split-step Fourier implementation. Let us define a Fourier transform of the envelope W(x, z) as follows:

W~x, p=JWx, z=dzWx, z exp jpz.                Eq. 44
Apparently p is a vertical component of the wave vector of the incident field W(x, z). The fundamental assumption in a split-step solution to Eq. (43) is that a Fourier transform of the envelope W(x, z) is performed while treating the term
k2εmz1
as a constant, i. e., no dependence on the z-coordinate. With that assumption, we obtain a marching solution:
W~x+x, p=W~x, p exp k2emz1p2x2jk.             Eq. 45

The inverse Fourier transform is given by:

Wx+x, z=expjk2εmz1x·J1W~x, p exp jp2x2k.             Eq. 46
In the above equation the term
k2εmz1
is no longer a constant, its variations actually results in variation of the angular spectrum p of the marching field W(x+Δx, z) with the distance.
The boundary condition in the case of ideal reflection (a perfectly conducting boundary surface) can be realized by adding the mirror image of the field W(x, z) in the upper half-space (z > 0) into the lower half space (z < 0) with an appropriate phase to reproduce the odd or even image of the field W relative to z = 0. The odd composition satisfies a condition W(x, z = 0) = 0 for a vertically polarized field, while the even composition results in the boundary condition
W/z|z=0
for horizontal polarization. In these cases of perfectly conducting boundary the Fourier transform reduces to a one-sided sine or cosine transform, respectively:
W~ox, p=2jJWox, z=2j0dzWox, z sin pz,
W~ex, p=2JWex, z=20dzWex, z cos pz                Eq. 47

where the subscripts e and o indicate that the field is even or odd, which, in turn, corresponds to horizontal or vertical polarization, respectively. The impedance boundary condition (Understanding Parabolic Approximation in Wave Propagation: Analytical Methods and Applicationsthis equation) was also treated in Ref. where it was shown that such a boundary condition can be realized using a mixed Fourier transform.

W~qx, p=0dzWx, zq sin pzp cos pz
where
q=jmη+1
Wx+x, z=exp jk2εmz1x·exp jq2x2kqz Kx+2π·0dpq sin pzp cos pzq2+p2 exp jp22kx·0dzWx, zq sin pzp cos pz

where the term K(x) is defined as:

Kx=2q0dzWx, z exp q, z;        Req>0,         Eq. 48
Kx=0;        Req>0.

It is apparent that the implementation of the split-step Fourier algorithm employing the formulas (47) and (48) will be more complicated than the rather simple implementation for a perfectly conducting interface using formula (46). As shown in many other publications (we may refer to a practical implementation of the numerical methods in radio coverage prediction systems), the value of the implementation of the impedance boundary conditions becomes pronounced at vertical polarization at smaller distances and at lower VHF frequencies for long-range propagation, where the impact of the surface wave produced by the distributed images of the source provides a substantial contribution to the received field strength.

The important development in an implementation of the split step Fourier transform method is the capability to treat an irregular terrain profile, as reported in numerous publications. The approach was initially developed in underwater acoustics and then implemented in the studies of radio propagation in the troposphere.

Let us consider the field of horizontal polarization over an irregular terrain with the terrain profile given by the function of the distance T(x). (We may note that with |q| ≫ 1 the boundary condition is W(x, z = 0) = 0 for both polarizations.) The range dependent boundary condition is then given by W(x, z = T(x)) = 0. The approach is then to map the irregular terrain profile to a smooth or more flat one with consequent modifications of the wave equation. The mapping is made by introducing a change of variables. Let the new height and range variables be represented by:

χ=x,  ς=zTx.       Eq. 49

In the new coordinate system the slow varying envelope W(x, z) can be written in the form:

Wx, z=wχ, ς exp (jθχ, ς       Eq. 50

where θ(χ, ς) is a phase correction term due to the irregular terrain. Substituting the new variables into Understanding Parabolic Approximation in Wave Propagation: Analytical Methods and Applicationsthis equation and following the procedure we end up with the following equation for a new envelope:

2wς2+j2kwχ+k2εmς+Tχ1+2ςTχ w=0       Eq. 51

and:

θχ, ς=kςTχk3/20χTα2 dα.       Eq. 52

As observed, the implementation of the split-step Fourier transform can be done in a way similar to the case of a smooth terrain, the difference from Understanding Parabolic Approximation in Wave Propagation: Analytical Methods and Applicationsthis equation is the presence of the second derivative of the boundary interface, 2ςT′′(χ), in the exponential operator. Nonetheless, given that the split-step Fourier transform algorithm treats this term as a constant, the addition of the second derivative does not make a difference in the implementation though ensuring the important correction of the wave front due to scattering on the irregular terrain.

We may notice apparent similarities between the split-step Fourier method and the discrete implementation of the path integral method. In the case of propagation over the boundary interface, considered in Section 3***, implementation of the boundary conditions into the path integral approach will result in the appearance of the mirrored source, similar to the above formulas (47) and (48). In fact, numerical calculation of the discrete version of the path integral is practically realised by applying the fast-fourier transform algorithm to the exponential propagator at given step along the range x.

Finally, we can state that the above numerical methods constitute a powerful framework for quantitative analysis of radiowave propagation through the tropospheric boundary layer in the very general case of the refractivity condition. In particular, the methods described above can incorporate either non-uniformity of the refractivity profiles in the horizontal plane as well as in an irregular terrain. The limitation of the above methods, as in fact of all numerical methods, is that while they are quite helpful in a quantitative estimation of the propagation phenomena when the physical mechanism is clearly understood they lack the capability to provide a framework for qualitative analysis, especially in the case of combined effects of refraction and scattering on random irregularities of refractive index or a randomly rough sea surface.

Author
Author photo - Olga Nesvetailova
Freelancer
Literature
  1. Fock, V.A. Electromagnetic Diffraction and Propagation Problems, Pergamon Press, Oxford, 1965.
  2. Fock, V.A. Distribution of currents excited by a plane wave on the surface of conductor, J. Exp. Theor. Phys., (Nauka), 1945, 15 (12), 693–702.
  3. Leontovitch, M.A. Method of solution to the problems of em wave propagation over the boundary surface, Izv. Acad. Nauk, Ser. Fiz., 1944, 8 (1), 16–22.
  4. Leontovitch, M.A., Fock, V.A. Solution to a problem of em wave propagation over the Earth’s surface using the method of parabolic equation, J. Exp.Theor. Phys., (Nauka), 1946, 16 (7), 557–573.
  5. Fock, V.A. Theory of wave propagation in a non-uniform atmosphere for elevated source, Izv. Acad. Nauk, Ser. Fiz., 1960, 14 (1), 70–94.
  6. Krasnushkin, A.E. The expansion over normal waves in a spherically-stratified medium, Dokl. Acad. Nauk SSSR, 1969, 185 (6), 1262–1265.
  7. Kravtsov, Yu.A., Orlov, Yu.I., Geometric Optics of Non-uniform Media, Nauka, Moscow, 1980, 340 pp.
  8. Rotheram, S. Radiowave propagation in the evaporation duct, The Marconi Rev., 1974, 67 (12), 18–40.
  9. Booker, H.G., Walkinshaw W. The mode theory of tropospheric refraction and its relation to waveguides and diffraction, in Meteorological Factors in Radiowave Propagation, The Royal Society, London, 1946, pp. 80–127.
  10. Andrianov, V.A. Diffraction of UHF in a bilinear model of the troposphere over the Earth’s surface, Radiophys. Quantum Electron., 1977, 22 (2), 212–222.
Footnotes
Sea-Man

Did you find mistake? Highlight and press CTRL+Enter

Март, 03, 2025 55 0
Add a comment


Notes
Text copied