SciELO - Scientific Electronic Library Online

 
vol.64 número1Electrones y huecos: perspectivas en torno al debate sobre el realismo científicoSome remarks on the Bel-Robinson tensor and gravitational radiation índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Revista mexicana de física E

versión impresa ISSN 1870-3542

Rev. mex. fís. E vol.64 no.1 México ene./jun. 2018

 

Education

The retarded potential of a non-homogeneous wave equation: introductory analysis through the Green functions

A. Téllez-Quiñonesa  * 

J.C. Valdiviezo-Navarroa 

A. Salazar-Garibaya 

A.A. López-Calocab 

aCONACYT-Centro de Investigación en Geografía y Geomática, Ing. J.L. Tamayo, A.C. (Unidad Mérida), Carretera Sierra Papacal-Chuburna Puerto Km 5, Sierra Papacal-Yucatán, 97302, México.

bCentro de Investigación en Geografía y Geomática, Ing. J.L. Tamayo, A.C., Contoy 137-Lomas de Padierna, Tlalpan-Ciudad de México, 14240, México.


Abstract

The retarded potential, a solution of the non-homogeneous wave equation, is a subject of particular interest in many physics and engineering applications. Examples of such applications may be the problem of solving the wave equation involved in the emission and reception of a signal in a synthetic aperture radar (SAR), scattering and backscattering, and general electrodynamics for media free of magnetic charges. However, the construction of this potential solution is based on the theory of distributions, a topic that requires special care and time to be understood with mathematical rigor. Thus, the goal of this study is to provide an introductory analysis, with a medium level of formalism, on the construction of this potential solution and the handling of Green functions represented by sequences of well-behaved approximating functions.

Keywords: Mathematical methods in physics; mathematics; diffraction theory; backscattering; radar

PACS: 01.30.Rr; 02.30.Em; 41.20.Jb

1. Introduction

Potential theory can be simply understood as the art of solving a linear distributional non-homogeneous differential equation through the Green functions 1,2,3. In the context of this study, our interest resides in the construction of an integral solution that derives from the divergence Gauss theorem, certain Green identities, and the handling of the Green functions. In the construction of this potential solution we also find other results of great importance, such as the integral theorem of Helmholtz and Kirchhoff, which is the main result that supports the scalar diffraction theory in optics 4. However, Green functions are not properly functions in the usual sense, since they are formally defined as distributions. Distributional theory, Green functions, and the use of Green identities have been successfully implemented in many theoretical and applied works, e.g., SAR theory 5,6,7, scattering and wave propagation 8,9,10,11,12, wave diffraction and electrodynamics 13,14,15,16, phase unwrapping 17,18,19,20, etc.

The concept of distribution is not easy to explain and in most of the references where distributions are mentioned, such as basic courses in calculus 21, differential equations 22, linear systems 23, or Fourier analysis 24, a detailed explanation of such abstractions is not usually given. Hence, in order to have a more solid notion of the concept of distributions, or generalized functions, specialized literature should be consulted. This literature is specifically related to a coarse field of mathematics called functional analysis 25,26 and some notions on Lebesgue measure 27. Certainly, in this work we will not provide specific and detailed information about distributions. Instead of that, we are going to use the artifice of working with sequences of well-behaved approximating functions 28, which permit us to talk about distributions concisely and without too much complexity. However, we hope to motivate the reader in the study of distributions through one of their most important applications: the standard solution of the non-homogeneous wave equation, also known as the retarded potential. To understand the construction of this integral solution, we first need to expose the main problem involved with the non-homogeneous wave equation. Thus, in order to establish the context of such a problem, we start in Sec. 2 with a typical deduction of the non-homogeneous wave equation from the Maxwell equations. Motivated in the perspective of SAR theory, in the subsequent sections, we provide some descriptive guidelines for constructing the potential solution of this non-homogeneous wave equation. The potential solution is supported by the divergence Gauss theorem and the Green identities, described in Sec. 3. Of course, the definition of Green functions is also explained in Sec. 4, where we discuss about certain incongruities when working with Green functions, found in certain references (for example 21). These incongruities refer to the question: How can be demonstrated that a function is a Green function with some of formalism? We emphasize that we do not want to criticize the descriptive style of references 21,22,23 or 24, which are indeed quite advisable. We simply want to point out that the treatment of distributions should be made with more care. Additionally to this discussion, a description of the integral theorem of Helmholtz and Kirchhoff, and its relation with the retarded potential is presented in Sec. 5. Finally, Sec. 6 outlines our conclusions.

2. Derivation of the non-homogeneous wave equation

It is well known that, an electromagnetic wave, such as a radar wave, is characterized in each point of the space x^=(x,y,z), and each time t, by the vectorial functions E=E(t,x^), the electric field (EF), and H=H(t,x^), the magnetic field (MF). When assuming as known the variables J=J(t,x^), the current density (vectorial field), σ=σ(t,x^), the charge density (scalar function), q, the permittivity constant, and p, the permeability constant, the fields E and H can be found. These fields are determined by the Maxwell equations: E=σ/q, the Gauss law for EF, H=0, the Gauss law for MF, ×E+p(H/t)=0^, the Faraday law, and ×H-q(E/t)=J, the Ampère law. Here, :=(/x,/y,/z), 0^=(0,0,0), and the media is assumed free of magnetics sources. In this case, operations F and ×F refer, respectively, to the divergence and the rotational of any vectorial field F=F(t,x^). Hence, a method to find E and H can be constructed from the previous Maxwell equations and the next theorems:

Theorem 1: Let F:R3R3 be a real valued vectorial field of class C 1 (continuous function with continuous first order partial derivatives), exceptionally in a finite number of points. Then, F is the gradient of some scalar function f (that is, F=f), if and only if ×F=0^.

Theorem 2: If F:R3R3 is a C 1-vectorial field such that F=0, then, there is a C 1-vectorial field G such that F=×G.

These theorems are demonstrated in 21 for real valued functions dependent on variable x^ (where F=F(x^)), however, they can be generalized to complex valued functions of the form F:R3C3. Moreover, these theorems are independent of variable t, so, they are also true for complex fields 29 of the form F:R4C3 where F=F(t,x^). Now, by assuming E and H as C 2-fields (that means, C 1-fields with continuous second order partial derivatives), from the Gauss law for MF, and Theorem 2, there exists a field A 0 such that

H=×A0. (1)

Additionally, from Faraday law and Eq. (1), we have that

0^=×E+pt(×A0)=×E+pA0t, (2)

which, from Theorem 1, implies the existence of a function μ0 satisfying

E+pA0t=μ0. (3)

On the other hand, from the Ampère law, Eqs.(1) and (3), and the vectorial identity ×(×A0)=(A0)-2A0 (where 2:=(2/x2)+(2/y2)+(2/z2) is the Laplacian operator), it follows that

J =×H-qEt=××A0-qtμ0-pA0t=A0-2A0-qtμ0+pq2A0t2 (4)

From the previous relation, it is obtained that

2A0-pq2A0t2=-J+(A0)-qt(μ0)=-J+(A0)-(qμ0t)=-J+(A0-qμ0t). (5)

Now, by using the Gauss law for EF and Eq. (3), we get

2μ0=σq+pt(A0). (6)

Let us define f as a scalar solution of the equation

2f-pq2ft2=-A0-qμ0t, (7)

and μ, as the function μ:=p(f/t)+μ0. From these definitions we obtain the relation

ft=(μ-μ0)p. (8)

Thus, it can be noticed that

H=×A, (9)

with A=A0+f. The result in Eq. (9) is logical because ×(f)=0^ for all differentiable scalar field f; in other words, ×A=×A0 for such A. Since H can be calculated, as suggested in Eq. (9), there is a function μ1 that satisfies

E+pAt=μ1, (10)

in analogy to the steps given for deriving Eq. (3). However, μ=μ1 due to the fact that

μ1=E+pt(A0+f)=(E+pA0t)+pt(f)=μ0+p(ft)=μ, (11)

a consequence of Eqs. (3), (8), and (10). In this case, we can replace μ1 by μ in Eq. (10) in order to obtain

E+pAt=μ. (12)

Thus, by using Eqs. (9) and (12), this A must satisfy

2A-pq2At2=-J+A-qμt (13)

and

2μ=σq+pt(A), (14)

in analogy to the steps for concluding Eqs. (5) and (6) from Eqs. (1) and (3). However, we have now a simplification in our calculations because

A-qμt=0, (15)

as the reader can confirm.

Since Eq. (15) takes place, then Eq. (13) reduces to

2A-pq2At2=-J. (16)

In a similar fashion, from Eq. (15), the expression in Eq. (14) is rewritten as

2μ-pq2μt2=σq. (17)

Therefore, by considering the wave propagation velocity c0:=1/pq (as defined in 30), and introducing the function ζ:=-σ/q, Eqs. (16) and (17) are correspondingly expressed as

1c022At2-2A=J,                            18a1c022μt2-2μ=ζ.                            (18b)

Then, when assuming J and ζ as known functions, the solutions of Eqs. (18) for A and μ, permit us to find H and E, from Eqs. (9) and (12), respectively. Equation (18) for μ, and analogously for A, is known as the non-homogeneous wave equation, or d’Alembert equation5 in the distributional sense. Due to its undulatory nature, the solution of this equation is called scalar wave in the case of μ, or vectorial wave, in the case of A.

3. The Helmholtz equation, the divergence Gauss theorem, and the Green identities

For constructing a solution of the non-homogeneous wave equation, from the methods applied for solving differential equations 22,31,32, it is typically proposed a general solution of the form μ=μ0+μ1, where μ0 is solution of the homogeneous version of the wave equation (when ζ=0), and μ1 is solution of the original non-homogeneous equation (when ζ0). In radar language 33,34,35, function ζ is interpreted as the source or the electric pulse emitted by an antenna. In this context, when ζ0, it is understood that the pulse induces a propagating wave which reaches an object. This wave is known as emitting or incident field and it is denoted by μ1. In opposite form, the object, assumed to be like a non-emitting electric pulse source (ζ=0), reflects or backscatters the incident wave. So, the reflected wave or backscattered field is denoted by μ0 (details of this conception are explained in 5). Consequently, the general solution μ is known as total field. Since the total field can be found by solving first the homogeneous wave equation, we are going to focus our attention in the construction of μ0. Thus, in this section, variable μ0 will be denote as μ for simplifying notation. Moreover, it is important to remark that any component of the vectorial function A in Eq.(18), with J=0^, can be found just as solving the scalar case for μ in Eq.(18), with ζ=0. In this sense, it will be sufficient to establish the theory for solving this scalar case.

Let us consider a scalar wave with spatial period or wavelength λ0, refractive index media n = 1 (air) and angular-temporal frequency ω0. Such wave can be represented by the function μ̃(t,x^)=a(x^)cos(ω0t+ϕ(x^)), where a(x^) is the amplitude, and ϕ(x^) is the phase. This cosinusoidal form is well known from the classic theory for solving the homogeneous wave equation; specifically speaking, the method called separation of variables, and the Fourier series22. However, in an equivalent way, the form of the wave can be generalized to the complex representation μ(t,x^)=f(x^)eiω0t, where μ̃ is the real part of μ, f(x^)=a(x^)eiϕ(x^) is the spatial part of μ (the so-called phasor or complex perturbation4), and i is the imaginary unit. Since this complex representation must satisfy the homogeneous wave equation, we have that 2μ-(1/c02)(2μ/t2)=0, where c0=λ0/T0 is the wave propagation velocity and T0 denotes the temporal period of such wave. In this case τ0=1/T0 would be called simply as the temporal frequency, where ω0=2πτ0. Thus, when substituting the complex representation of μ in the homogeneous wave equation, we obtain that eiω0t2f(x^)-(i2ω02/c02)f(x^)eiω0t=0, which implies

(2+k02)f(x^)=0. (19)

This last expression is called Helmholtz equation, where k0=ω0/c0=(2πτ0)/(λ0τ0)=(2π)/λ0 is the wave number or the angular-spatial frequency.

On the other hand, a well-known result in vectorial calculus is the divergence Gauss theorem21, which can be written as

Ω3FdV=Ω2Fn^dA, (20)

where F=F(x^) is a C 1-vectorial field on Ω such that F:R3R3, and Ω is an elemental region in R3 with positive parametrized boundary Ω. The surface Ω can be a sphere, an ellipsoid, a parallelepiped, etc. In this theorem, symbol Ω3 denotes the triple integral over the region Ω, Ω2 denotes the double integral over the surface Ω, dV is a volume differential element, dA is an area differential element, and symbol denotes dot product. In addition, vectorial function n^ represents the unitary normal vector with respect to the surface Ω, in such a way that, it points towards the outside of the surface. So, if we take F=fg, where f=f(x^) and g=g(x^) are two differentiable scalar fields from R3 to R, then

F=(fg)=fg+f2g. (21)

When substituting this equation in the divergence Gauss relation, and introducing the directional derivative (g/n):=gn^, we obtain

Ω3fgdV+Ω3f2gdV=Ω2fgndA. (22)

Equation (22) is known as the first Green identity. Analogously, when considering F=gf, the equality

Ω3gfdV+Ω3g2fdV=Ω2gfndA, (23)

is deduced. By subtracting Eq. (23) from Eq. (22), it is concluded that

Ω3(f2g-g2f)dV=Ω2fgn-gfndA, (24)

a relation called second Green identity. It is important to remark that the divergence Gauss theorem is also valid for C 1-complex vectorial fields F:R3C3 that can be expressed in terms of complex scalar fields f,g:R3C. This last includes the possibility of considering fields of the form F=F(t,x^), f=f(t,x^), and g=g(t,x^), which means that F:R4C3, and f,g:R4C. The validity of this theorem is a consequence of the linearity of the integrals for complex valued expressions that can be denoted as F = Re(F) + iIm(F), where Re(F) and Im(F), are the real and the complex parts of F, respectively. Thus, the only required condition is to have Re(F) and Im(F), as C 1-real valued functions.

4. The Green functions: language of distributions

Informally, any function g that satisfies the Helmholtz equation almost everywhere26,27 on R3, could be called Green function4, however, such g requires to satisfy another properties in the context of distributions26. In this sense, expression almost everywhere refers to a property which is satisfied at all points of a domain with the exception of the points in a zero volume subset of the domain. For this particular case, the property would be the fulfilling of the Helmholtz equation. On the other hand, the Green functions and the Green identities are important and useful, specially in optics, for establishing a transcendental theorem: the integral theorem of Helmholtz and Kirchhoff. This theorem is the key result for supporting the scalar diffraction theory and it is related in part with the solution of a particular distributional non-homogeneous wave equation, the so-called retarded potential of the d’Alembert equation5. Now, let L be a linear operator applied to scalar functions depending on x^. Formally, a distribution G(x^,y^) is said to be a Green function with respect to L, if it satisfies: a) G(x^,y^)=G(y^,x^) for all x^=(x,y,z) and y^=(u,v,w) in R3, and b) L[Gx^,y^]=δ(y^-x^), where δ(y^-x^) is the Dirac delta distribution. However, every distribution D is said to be a Dirac delta, if it satisfies: 1) D(y^-x^)=0, for all y^x^, and 2)

R33D(y^-x^)dV(y^)=1.

In this definition dV(y^) refers to a volume differential element with respect to the integration variable y^. Thus, when an arbitrary distribution D fulfills properties 1) and 2), we write δ:=D.

Let us consider G(x^,y^):=gy^(x^), where gy^(x^)=eik0||x^-y^||/||x^-y^||. Here, |||| denotes the Euclidean norm for vectors in R3. This function satisfies gy^(x^)=gx^(y^) by symmetry, then G(x^,y^) fulfills property a). Is this G a Green function? What linear operator could be related to this G to declare it as a Green function? Well, the obvious answer is that such operator must be involved with the Helmholtz equation. So, if we consider L:=[-14π](2+k02), then Lfx^=[-14π](2+k02)f(x^) for any smooth function f:R3C. Moreover, Dy^-x^:=LGx^,y^=[-14π](2+k02)gy^(x^)=0 for all x^y^, from Proposition 1 in Appendix A. This also implies that property 1) is satisfied by D(y^-x^), every time that x^y^. Thus, assuming this D as a simple function depending on y^, which is a discontinuous function in y^=x^, we rigorously have that

R33D(y^-x^)dV(y^)=0

by using improper integrals. Even in the case of more complex integrals, if a function is zero almost everywhere in certain domain, then its Lebesgue integral on such domain should be zero 27. Thus, since property 2) is not achieved by this D, it implies that G(x^,y^)=gy^(x^) is not a Green function. But, why so many references 2,3,4,8 declare gy^(x^) as a Green function? Well, may be the answer is in the interpretation of D, and consequently G, as distributions. In the spirit of considering a sequence of well-behaved approximating functions (see the Remarks in Appendix A), G can be rewritten to the equivalent form

G(x^,y^):=limγ0+eik0rfγ(r), (25)

where r=||x^-y^|| and limγ0+fγ(r)=1/r for r > 0. In this case for each γ>0, fγ is a smooth function for all r > 0 and right continuous at r = 0. In other words, dfγ/dr exists for all r > 0, and limr0+fγ(r)=fγ(0), respectively. Now, it is evident that G in Eq. (25) satisfies properties a) and 1), however, if property 2) is required for D with this new G, we need to construct function fγ conveniently. There are many forms to do this construction, but we are going to propose a particular one. Before starting with this construction, let us think again that G(x^,y^) is gy^(x^), and consider the next argumentation:

From property 1), it is inferred that the integral of D on R3 is the same that the integral of D on any open ball Bε(x^)={y^R3:||y^-x^||<ε}. Then,

R33D(y^-x^)dV(y^)=(-14π)×Bε(x^)3(2+k02)gy^(x^)dV(y^)=(-14π)×[Bεx^32gy^x^dVy^+k02Bεx^3gy^x^dVy^], (26)

for all ε>0 fixed. Thus, when considering a translation to the origin, a change to spherical coordinates, and the use of improper integrals, we have that

Bε(x^)3gy^(x^)dV(y^)=Bε(x^)3gx^(y^)dV(y^)=4πk02[eik0ε(1-ik0ε)-1], (27)

independently of the discontinuity of gx^(y^) at y^=x^. Now,

Bε(x^)32gy^(x^)dV(y^)=Bε(x^)3gy^(x^)dV(y^)=-Bε(x^)3y^gy^(x^)dV(y^)=-Bε(x^)2gy^(x^)n^(y^)dA(y^), (28)

where n^(y^) is the unitary normal vector to the surface Bε(x^) and dA(y^) is an area differential element with respect to y^. This is a consequence of Proposition 2 in Appendix A and the divergence Gauss theorem in variable y^. Following the calculus of the previous triple integral (now a double integral) we have that

=-Bεx^2eik0r[ik0r-1r3](x^-y^)y^-x^y^-x^dA(y^)=Bεx^2eik0r[ik0r-1r3]×(r2r)dA(y^)=4πeik0ε(ik0ε-1), (29)

where the values r=||x^-y^|| in the double integrals are all equal to ε, because y^Bε(x^)={y^R3:||y^-x^||=ε}. When considering Eqs. (27), (28), and (29) in connection with Eq. (26), we get property 2). Therefore, function gy^(x^) results to be a Green function.

The last argumentation, although desirable, is false and its main fail is: function gy^(x^)=eik0r(ik0r-1)(x^-y^)/r3 is not a C 1-function (a smooth function) on Bε(x^), specifically at y^=x^ where r = 0. Then, we have incorrectly applied the divergence Gauss theorem in Eq. (28), which requires smoothness for the vector field gy^(x^). However, this mistake motivates us to think about an appropriate election of function fγ in Eq. (25). The sequence of functions fγ satisfy limγ0+eik0rfγ(r)=eik0r/r=gx^(y^) for all y^x^ (r > 0), but that is not enough. For each γ>0, we also need smoothness for the terms eik0rfγ(r) for all y^Bε(x^). However, since the term eik0r is smooth for all r0, then a sequence of functions of the form

fγ(r):=p0(r)rγpγ(r)0r<γ, (30)

where p 0(r) := 1/r, could be useful. For each γ>0, pγ(t) could be a convenient polynomial function, in such a way that pγ(γ)=p0(γ), pγ'(γ)=p0'(γ), and pγ(γ)=p0(γ). Here, symbols ´ and ´´ denote, correspondingly, the first and the second derivatives with respect to r. Hence, these conditions would warrant the smoothness of fγ for all r > 0. However, if pγ(r) is a linear combination of r-powers and we want to analyze the behavior of eik0rfγ(r) (particularly at r = 0), we first need to calculate the resultant expressions of the partial derivatives of eik0rrn, with respect to x, y, and z. For instance, from the chain rule we have

x[eik0rrn]=eik0r[nrn-1+ik0rn](x-u)r, (31)

for x^y^ (or r > 0). This expression is not necessarily right continuous at r = 0 when considering the definition of continuity by lateral limits. However, if we achieved to avoid divisions by r (removing the discontinuity at r = 0), the formula in Eq. (31) would be better behaved. Thus, by considering integer powers of r such that n2 and defining the vector i^:=(1,0,0), we get

x[eik0rrn]|x^=y^:=limh0eik0rrn|x^=y^+hi^-eik0rrn|x^=y^h=limh0eik0|h||h|nh=0, (32)

due to the fact that

limh0+eik0|h||h|nh=limh0-eik0|h||h|nh=0,

as the reader can confirm. Moreover, from Eq.(31) it follows that

limx^y^x[eik0rrn]=limr0+xu{eik0r[nrn-1+ik0rn](x-u)r}=limr0+xu{eik0r[nrn-2+ik0rn-1](x-u)}=0. (33)

In consequence, Eqs. (31)-(33) warrant continuity for the term &[eik0rrn]/x for all x^y^, and also for x^=y^, when considering this term as function of x^. Nevertheless, due to the radial symmetry of eik0r[nrn-1+ik0rn]/r in Eq. (31), term [eik0rrn]/x, interpreted as function of y^, must be also a continuous function. The same reasoning applies to [eik0rrn]/y and [eik0rrn]/z; therefore, if fγ(r) is a polynomial of powers n2 for r values such that 0r<γ, then eik0rfγ(r) will be continuous for all points y^ corresponding to those r values. Moreover, without loss of generality, we can assume that pγ(r)=ar2+br3+cr4+dr5 with pγ(γ/2)=p0(γ/2), which implies that

pγ(r)=(42/γ3)r2-(111/γ4)r3+(102/γ5)r4-(32/γ6)r5. (34)

So, from Eqs. (30) and (34) we get

eik0rfγ(r)=eik0rA(r)(x^-y^)rγeik0rB(r)(x^-y^)0r<γ, (35)

where Ar=[ik0r-1]/r3 and

B(r)=84γ3-333rγ4+408r2γ5-160r3γ6+ik0[42rγ3-111r2γ4+102r3γ5-32r4γ6], (36)

are two functions such that limrγ+A(r)=A(γ)=B(γ)=limrγ-B(r), limrγ+A'(r)=A'(γ)=B'(γ)=limrγ-B'(r), limr0+B(r)=B(0), and limr0+B'(r)=B'(0). Equation (35) is now a smooth function (C 1-class) for all y^x^ and also for y^=x^. In the same way, term eik0rfγ(r) is another smooth function for all y^ by construction. From this construction and Eq. (25), we can calculate again

R33D(y^-x^)dV(y^)=-14πlimγ0+[Bε(x^)32eik0rfγ(r)dV(y^)+k02Bεx^3eik0rfγ(r)dV(y^)], (37)

for any ε>0 arbitrary and fixed. But this time we have

Bε(x^)3eik0rfγ(r)dV(y^)=Bε(x^)Bγ(x^)3eik0rrdV(y^)+Bγ(x^)3eik0rpγ(r)dV(y^)=4πk02[eik0ε1-ik0ε-eik0γ1-ik0γ]+eik0rpγ(r)4πγ33, (38)

by definition of fγ when considering γ<ε. In the last expression, the integral of (eik0r/r) can be calculated by a translation to the origin and spherical coordinates, while the integral of eik0rpγ(r) is obtained from the mean value theorem for integrals21,36 due to the integrand’s continuity. In this case r=||x^-y^|| with y^Bγ(x^)Bγ(x^), and the factor (4πγ3/3) is the volume of the ball Bγ(x^). On the other hand, when taking the complex module of the integral in Eq. (38), we have that

|Bε(x^)3eik0rfγ(r)dV(y^)|4πk02|eik0ε(1-ik0ε)-eik0γ(1-ik0γ)|+4πγ33|pγ(r)|. (39)

Additionally, from Eq. (34), the triangle inequality, and the fact that rγ, it follows that

4πγ33|pγ(r)|1148πγ23. (40)

In an extreme case (4πγ3/3)|pγ(r)|(4πγ2/3) for r-values close or equal to γ, but the term (4πγ3/3)|pγ(r)| is always dominated by proportional factors to γ2. Then, when considering very small values of γ (the limit context of Eq. (37)), the inequalities in Eqs. (39) and (40) permit us to establish that eik0rpγ(r)(4πγ3/3)0 and

k02Bε(x^)3eik0rfγ(r)dV(y^)4π[eik0ε(1-ik0ε)-1], (41)

from Eq. (38). Now, due to the form of term eik0rfγ(r) in Eq. (35) and Proposition 2, it follows that 2eik0rfγ(r)=eik0rfγ(r)=-y^eik0rfγ(r). Therefore,

Bε(x^)32eik0rfγ(r)dV(y^)=-Bε(x^)3y^eik0rfγ(r)dV(y^)=-Bε(x^)2eik0rA(r)(x^-y^)(y^-x^)||y^-x^||dA(y^)==4πeik0ε(ik0ε-1), (42)

as a consequence of the divergence Gauss theorem in variable y^, applied to the field eik0rfγ(r), and the formula for eik0rfγ(r) when γ<ε=r. Naturally, the reduction of the calculations in the double integral of Eq. (42) is implied from the fact that y^Bε(x^). So, since the smoothness of the field eik0rfγ(r) on Bε(x^) is warranted in this case, then function G in Eq. (25) satisfies property 2), when considering Eqs. (41), (42), and (37). In conclusion, G(x^,y^) in Eq. (25), with fγ in Eq. (30) and pγ in Eq. (34), is a true Green function. And of course, this G is understood as the limit of a sequence of smooth functions that converge to gy^(x^) almost everywhere. In mathematical terms, G(x^,y^)=gy^(x^) for all x^y^, except in x^=y^, where gy^(y^) is not defined and G(y^,y^)=0. Informally speaking, that is why gy^(x^) inherits the name of its equivalent G(x^,y^).

Finally, we could be attempted to believe that G is an ordinary function like

G(x^,y^)=eik0||x^-y^||/||x^-y^||forx^y^0forx^=y^, (43)

and that is true, in some way, with respect to the image sets induced by both expressions. Although formally

R33L[G]dV(y^)=0,

there is no obstacle to call G as a Green function because G is also a generalization of G or, equivalently, G=gy^ almost everywhere. Something similar happens when thinking in the one dimensional step function, also known as the Heaviside step. The derivative of the Heaviside step is zero at all points, except in the jump discontinuity where it is not defined. The integral of this derivative along the real line is identically zero. However, if the Heaviside step is understood as a distribution, then its weak derivative26 is a Dirac delta (see the remarks in Appendix A). As well known, the integral of this weak derivative along the real line is identically one. These affirmations may seem contradictory, but they are only a question of abstract interpretation.

5. The integral theorem of Helmholtz and Kirchhoff and the retarded potential

Let us consider an elemental region Ω in R3 as a closed set, in such a way that its boundary is conformed by two independent surfaces. This means that Ω=SBε(x^0) with SBε(x^0)=, where S is a surface that bounds Bε(x^0) and denotes the empty set. Also, we are going to assume that S is smooth by parts, where SΩ and Bε(x^0)Ω, due to the fact that Ω is closed. So, region Ω can be assumed as a glass ovoid with an inside air sphere, where the elliptical surface of the ovoid corresponds to S, and the surface of the air sphere is Bε(x^0). It should be understood that Ω limits two exterior zones and one inner zone (see Fig. 1): the big exterior zone given by the open set (ΩBε(x^0))c={x^R3:x^ΩBε(x^0)}, the small exterior zone given by Bε(x^0), and the inner zone Ω\Ω=Ω\(SBε(x0))={x^Ω:x^SBε(x0)}.

Figure 1. Sketch of the region Ω. In this example, it is assumed as an ovoid with an inside sphere of radius ε

Therefore, when considering possible parametrizations of the surfaces S and Bε(x^0), we must think that the unitary normals n^(x^) should point towards the big exterior zone when x^S, and point towards the small exterior zone when x^Bε(x^0), respectively. So, if we apply the second Green identity to a C 1-phasor f(x^)=a(x^)eiϕ(x^), corresponding to some μ solution of the homogeneous wave equation on Ω, and the function g(x^)=eik0r/r with r=||x^-x^0||, then we get

Ω2fgn-gfndA=Ω3(f2g-g2f)dV. (44)

However,

Ω3(f2g-g2f)dV=Ω3(-k02fg+k02gf)dV=0, (45)

because both the phasor f and function g satisfy the Helmholtz equation for all points in Ω. This implies that

Ω2(fgn-gfn)dA=S2(fgn-gfn)dA+Bε(x^0)2(fgn-gfn)dA=0, (46)

or, equivalently,

Bε(x^0)2(fgn-gfn)dA=-S2(fgn-gfn)dA. (47)

Since the double integral on the right hand side of Eq.(47) is a constant, independently of the ε value, it follows that

limε0+Bε(x^0)2(fgn-gfn)dA=-S2(fgn-gfn)dA, (48)

including the possibility that Ω=S{x^0} for the limit case. Now, the integral on Bε(x^0) in Eqs.(47) and (48) can be expressed as

Bε(x^0)2(fgn-gfn)dA=4πε2(fgn-gfn)|x^*, (49)

from the mean value theorem of integrals, where 4πε2 is the area of Bε(x^0) and x^* is some fixed point in Bε(x^0). Thus, if x^*Bε(x^0), then

g(x^*)=eik0εε,                                                                    (50a)g(x^*)=eik0εik0ε-1ε3x^*-x^0,                               (50b)

by using the formulas in Proposition 1. Owing to the fact that vector n^ points towards the small exterior zone in this particular case, we have that

gn(x^*)=gn^x^*, (51a)

n^(x^*)=x^0-x^*||x^0-x^*||=-x^*-x^0ε, (51b)

then

gn(x^*)=eik0ε(1-ik0ε)ε2, (52)

from Eqs. (50) and (51). Hence

Bε(x^0)2(fgn-gfn)dA=4π[fx^*eik0ε1-ik0ε-εeik0εfnx^*], (53)

after substituting Eqs. (50) and (52) into Eq. (49). If ε0+, then x^*x^0 and

limε0+Bε(x^0)2fgn-gfndA=4πf(x^0), (54)

from Eq. (53) and the fact that f and f/n are continuous functions (f is C 1-function). The limit in Eq. (54) represents the same to write

f(x^0)=14πS2gfn-fgndA, (55)

by considering Eqs. (48) and (54). Equation (55) is valid for any C 1-phasor f:R3C, where g(x^)=eik0r/r and r=||x^-x^0||. In this case S{x^0}=Ω is the union of a smooth-by-parts surface S with the frontier point x^0. Examples of smooth-by parts surfaces could be a sphere, an ellipsoid, the three faces of triangular pyramid, the six faces of a parallelepiped, etc. So, without loss of generality and in the context of distributions, the result in Eq. (55) is also valid when simply assuming S=Ω in such a way that x^0Ω\S. Equation (55) corresponds to the integral theorem of Helmholtz and Kirchhoff, which can be generalized to any g function satisfying the Helmholtz equation in an equivalent sense. In other words, g in Eq. (55) could be any Green function. In this argumentation we have assumed that L[f] = 0 for all x^Ω, but such hypothesis can be modified to an equivalent case as follows: Let G be a Green function with respect to the linear operator L, and let h be a distribution similar to h0δ(x^-x^1) with a constant factor (an independent term with respect to variable x^) h 0, and x^1Ω\Ω. Thus, if f is a function such that L[fx^]=h(x^), then

f(y^)=14πΩ2(G(x^,y^)fn(x^)-f(x^)Gn(x^,y^))dA(x^), (56)

for all y^Ω\Ω. Equation (56) can be simply deduced, for instance, from the second Green identity we have

Ω2(Gfn-fGn)dA(x^)=Ω3(G2f-f2G)dV(x^)=Ω3G(2+k02)fdV(x^)-Ω3f(2+k02)GdV(x^), (57)

but this is the same that

=-4πΩ3GL[f]dV(x)+4πΩ3fL[G]dV(x^)=-4πΩ3G(y^,x^)h(x^)dV(x^)+4πΩ3f(x^)δ(x^-y^)dV(x^)=-4πΩ3G(y^,x^)h(x^)dV(x^)+4πf(y^), (58)

from definition of G. Then, from Eqs.(57) and (58) we get

f(y^)=14πΩ2(Gfn-fGn)dA(x^)+Ω3G(y^,x^)h(x^)dV(x^). (59)

So, if h(x^)=h0δ(x^-x^1) with x^1Ω\Ω, then the triple integral on Ω in Eq. (59) is identically zero from the properties of the Dirac delta. This is because

Ω3G(y^,x^)δ(x^-x^1)dV(x^)=G(y^,x^1),

only if x^1Ω\Ω. Moreover, for all x^Ω\Ω, the integrand G(x^,y^)h(x^)=h0G(x^,y^)δ(x^-x^1)=0 does not have any discontinuity in Ω\Ω when considering G in Eq. (25). Thus, by interchanging x^ by y^ and vice-versa in Eq. (59), we get a potential solution for the problem of solving L[fx^]=h(x^) for all x^Ω\Ω. Such a problem would be solved with the next hypothesis: the knowledge of functions f and f/n on Ω, a given function h, preferable smooth and defined on Ω, and a possible and convenient Green function G, defined with respect to the operator L. The inferred solution from Eq. (59) can be expressed as f(x^)=f0(x^)+f1(x^) with

f0(x^):=14πΩ2Gfn-fGndA(y^) (60)

and

f1(x^):=Ω3G(x^,y^)h(y^)dV(y^), (61)

where

L[f1x^]=Ω3L[Gx^,y^]h(y^^)dV(y^)=Ω3δ(y^-x^)h(y^)dV(y^)=h(x^). (62)

Since L[fx^]=h(x^) by hypothesis, it follows that L[f0x^]=0. In other words, f 0 and f 1 are particular solutions of the homogeneous and the non-homogeneous Helmholtz equation, respectively. We say that f 0 and f 1 are “particular” functions, because both of them depend on the G chosen. Eventually, beyond of considering a bounded set as Ω, if we think in solving L[f] = h on R3, the result in Eq. (59) would be equivalent to consider

f(x^)=R33G(x^,y^)h(y^)dV(y^), (63)

as the potential solution.

On the other hand, the theory of Green functions, exposed in a previous section, was described with respect to the operator -14π[2+k02]. However, this theory does not change if we consider other similar operators as (2+k02) or [2+(k0c0)2]. Thus, when considering L:=-14π[2+(k0c0)2], function

g0(x^):=-ei(-k0/c0)||x^-y^||||x^-y^||, (64)

satisfies that

2+k02c02g0(x^)=0, (65)

when assuming y^ as a constant vector and x^y^. To verify Eq. (65), it is only required to replace k 0 and x^0 in Proposition 1 by -k 0/c 0 and y^, respectively. Now, we want to solve

1c022t2-2μ(t,x^)=ζ(t,x^), (66)

on a domain Ω, where ζ (considered only as function of x^) behaves as a Dirac delta translated to some point outside of Ω\Ω. So, by proposing a solution of the form μ(t,x^)=f(x^)eik0t and substituting this solution in Eq. (66), we get

eik0t-2-k02c02f(x^)=ζ(t,x^). (67)

Therefore, when Eq.(67) is evaluated at (t-(||x^-y^||/c0),x^), it is obtained that

eik0(t-(||x^-y^||/c0))(-2-k02c02)f(x^)=ζ(t-(||x^-y^||/c0),x^). (68)

In some way, the equality in Eq. (67) suggests that ζ could be interpreted as a function with separable variables, time t, and position x^. This is similar to conveniently think that ζ(t,x^)=P(t)Q(x^), then Eq. (67) can be expressed as L[fx^]=h(x^), where h(x^)=h0Q(x^) and h0=h0(t)=P(t)e-ik0t/(4π) is a constant term with respect to variable x^. Consequently, if the spatial part of ζ is a function like Q(x^)=δ(x^-x^1) with x^1Ω\Ω, then the phasor f can be calculated by the integral theorem of Helmholtz and Kirchhoff, the second Green identity, and any convenient Green function G, as

f(y^)=14πΩ2(Gfn-fGn)dA(x^)=14πΩ3(G2f-f2G)dV(x^), (69)

for all y^Ω\Ω. In addition, when taking G=g0, it is concluded that

f(y^)=14πΩ3e-ik0||x^-y^||/c0||x^-y^||×(-2-k02c02)f(x^)dV(x^), (70)

from Eqs.(64), (65), and (69). Furthermore, Eq.(70) implies that

μ(t,y^)=14πΩ31||x^-y^||ζt-||x^-y^||c0,x^dV(x^), (71)

from the form assumed for μ, and Eqs. (68) and (70). The expression in Eq. (71) is known as the retarded potential of ζ, as mentioned in 5. Its name reveals that the signal μ is recovered from the source ζ with a delay in time. This potential represents a standard solution of the d’Alembert equation displayed in Eq. (66). Moreover, the formula in Eq. (71) is independent of the fact that ζ behaves as a Dirac delta. For instance, if we simply assume that ζ is a smooth function on Ω, then Eq. (68) can be expressed as

L[fx^]=h(x^):=Ht-||x^-y^||c0,x^, (72)

where H(t,x^)=(e-ik0t/(4π))ζ(t,x^). From Eq. (61), by interchanging x^ by y^, and considering G=-g0, it follows that

f(y^)=Ω3e-i(k0/c0)||x^-y^||||x^-y^||H(t-||x^-y^||c0,x^)dV(x^)=e-ik0t4πΩ31||x^-y^||ζ(t-||x^-y^||c0,x^)dV(x^), (73)

is a potential solution of the equation L [f] = h given in Eq. (72). Indeed, from the form assumed for μ, Eq. (73) induces again the formula given in Eq. (71). Finally, if we want to solve L [f] = h on R3, we could use Eq.(63) and conclude that

μ(t,y^)=14πR331||x^-y^||ζt-||x^-y^||c0,x^dV(x^), (74)

in analogy to the previous formulas that preserve the name of retarded potential. Nevertheless, Eq. (74) is also valid for the case when ζ is a Dirac delta δ, specially in the case of approximating this δ with a sequence of smooth functions like the Gaussians.

6. Discussion and Conclusions

We have derived the retarded potential of a non-homogeneous wave equation by considering certain subtle mathematical details. These details refer to the use of distributions, in our own and simplified interpretation of these generalized functions, and in what sense it is said that a given function is a Green function. According to our analysis, we obtained a distributional solution f for Eq. (72), which permits us to construct a solution μ for Eq. (66). When considering Eq. (72) in a bounded set Ω, the solution can take place by establishing boundary conditions in the frontier Ω, as exposed by Eq. (59). These conditions may be imposed in f, or in f/n, depending on the problem. For instance, in the simple case of h= 0 for all points in Ω, it follows that f = f 0 from Eq. (60). In this particular case, if we only know f on the frontier Ω (Dirichlet problem), then a desirable election of G could be a Green function that vanishes on this boundary.

In a general perspective, the base idea of the potential solutions is manifested by Eqs. (59)-(61) and (63), where the election of an appropriate Green function is crucial to obtain specific results. For example, the formula given by Eq. (71) in the bounded case, or by Eq. (74) in the unbounded case, respectively. The retarded potentials allow us to build solutions of classic problems in electrodynamics, wave propagators, radar, among others. For instance, potential theory can be applied for modeling equations related with the emission and detection of a SAR signal 5,6. In a SAR configuration, the main equation that involves the recovered values of the signal and the scattering object density is a consequence of using Eqs. (71) or (74), in connection with the first Born approximation. Such an approximation reduces the ill-posed problem of recovering the scattering object to a simple convolution-filtering problem.

About the mathematical rigor found in some references, an explicit and formal explanation about Green functions, by using operator theory, is found in 1. Although in this reference there is no mention on that gy^(x^)=eik0r/r is Green function with respect to the Helmholtz operator, an extensive analysis to derive Green functions from many different linear operators is given. Nevertheless, such analysis is out of the scope of these notes. On the other hand, it is interesting to notice that expression gy^(x^) is declared as a Green function in many books (for example 2, 3, 4, 8), without any formal proof of that fact. Whereas in some other books and for some other kind of Green functions, a proof is provided but with drawbacks. For instance, the argumentation in 21 when justifying f(x^,y^)=-1/(4π||x^-y^||) as a Green function with respect to a distributional Poisson equation. In that reference, the drawback is exactly the same that was exposed in Sec. 4 on the illegal use of the divergence Gauss theorem. Of course, we do not pretend to criticize that books because they are actually excellent references and we are far from exposing a formal proof. However, we expect at least to motivate the reader on the importance of considering sequences of approximating functions 28, to have a more clear idea about the handling of distributions. Such as made in Eq. (25), by considering these kind of sequences, an equivalence relation could be established between G(x^,y^), defined in Eq. (25), and function gy^(x^). Since G(x^,y^)=gy^(x^) is valid almost everywhere with respect to the variable y^R3 (or x^R3) 26,27, there is no doubt now of calling gy^(x^) as a Green function. In a similar manner, the argumentation in 21 when justifying f(x^,y^) as a Green function could be improved when considering sequences of smooth approximating functions. Just as Jackson explains in his analysis on Poisson and Laplace equations in reference 3, we want to emphasize a very good footnote which refers to the volume integral of Eq. (1.36) in that reference: “The reader may complain that (1.36) has been obtained in an illegal fashion since 1/|x - x´| is not well-behaved inside the volume V. Rigor can be restored by using a limiting process…” Well, such limiting process has been exemplified in this work.

Acknowledgment

The authors thank Dr. Gerardo García-Almeida for their valuable comments on this paper.

References

1. B. Friedman, Principles and Techniques of Applied Mathematics (John Wiley & Sons, Inc., NY, 1956). pp. 134-163. [ Links ]

2. G. Arfken, Mathematical Methods for Physicists, (Academic Press, Inc., NY, 1985). pp. 897-924. [ Links ]

3. J.D. Jackson, Classical Electrodynamics, (John Wiley & Sons, Inc., NJ, 1999), pp.34-37. [ Links ]

4. J.W. Goodman, Introduction to Fourier Optics, (McGraw-Hill, Inc., NY, 1968). pp. 30-40. [ Links ]

5. M. Gilman and S. Tsynkov, SIAM Journal on Imaging Sciences 8 (2015) 186. [ Links ]

6. M. Cheney, SIAM Review 43 (2001) 301. [ Links ]

7. S.V. Tsynkov, SIAM Journal on Imaging Sciences 2 (2009) 646. [ Links ]

8. M. Bertero and P. Boccacci, Introduction to inverse problems in imaging, (IOP Publishing Ltd, London, 1998). pp. 191-219. [ Links ]

9. A. Karlsson, H. Otterheim, and R. Stewart, J. Opt. Soc. Am. A 10 (1993) 886. [ Links ]

10. K.F. Ren, G. Gréhan, and G. Gouesbet, J. Opt. Soc. Am. A 14 (1997) 3014. [ Links ]

11. M.A. Alonso and G.W. Forbes, J. Opt. Soc. Am. A 14 (1997) 1076. [ Links ]

12. S.A. Prahl, D.G. Fischer, and D.D. Duncan, J. Opt. Soc. Am. A 26 (2009) 1533. [ Links ]

13. J.A. Stratton and L.J. Chu, Phys. Rev. 56 (1939) 99. [ Links ]

14. X. Wang, Z. Fan, and T. Tang, J. Opt. Soc. Am. A 23 (2006) 872. [ Links ]

15. V.I. Tatarskii, J. Opt. Soc. Am. A 12 (1995) 1254. [ Links ]

16. H. Ardavan, A. Ardavan, J. Singleton, J. Fasel, and A. Schmidt, J. Opt. Soc. Am. A 25 (2008) 543. [ Links ]

17. I. Lyuboshenko and H. Maître, J. Opt. Soc. Am. A 16 (1999) 378. [ Links ]

18. G. Fornaro, G. Franceschetti, and R. Lanari, IEEE Transactions on Geoscience and Remote Sensing 34 (1996) 720. [ Links ]

19. G. Fornaro, G. Franceschetti, R. Lanari, and E. Sansosti, J. Opt. Soc. Am. A 13 (1996) 2355. [ Links ]

20. S. Marano, F. Palmieri, and G. Franceschetti, J. Opt. Soc. Am. A 19 (2002) 1319. [ Links ]

21. J.E. Marsden and M.J. Hoffman, Vector Calculus, (W.H. Freeman and Company, NY, 1996). pp. 522-531. [ Links ]

22. D.G. Zill and M.R. Cullen, Differential Equations with Boundary-Value Problems, (Cengage Learning, Belmont, CA, 2009). pp. 292-295. [ Links ]

23. A.V. Oppenheim, A.S. Willsky, and S.H. Nawab, Signals and Systems, (Prentice-Hall, Inc., NY, 1997). pp. 32-38. [ Links ]

24. H.P. Hsu, Fourier Analysis, (Simon & Schuster, Inc., NY, 1970). pp. 37-51. [ Links ]

25. E. Kreyszig, Introductory Functional Analysis with Applications, (John Wiley & Sons, NY, 1978). pp. 50-110. [ Links ]

26. E.H. Lieb and M. Loss, Analysis, (AMS, Rhode Island, 2001). pp. 135-158. [ Links ]

27. R.G. Bartle, The Elements of Integration and Lebesgue Measure, (John Wiley & Sons, NY , 1995). pp. 6-51. [ Links ]

28. R.N. Bracewell, The Fourier Transform and its Applications, (McGraw-Hill, Boston, 2000). pp. 74-104. [ Links ]

29. J.E. Marsden and M.J. Hoffman, Basic Complex Analysis, (W.H. Freeman and Company, NY , 1999). pp. 95-110. [ Links ]

30. H. Maître, Processing of Synthetic Aperture Radar Images, (John Wiley & Sons, Inc., NJ, 2008). pp. 1-6. [ Links ]

31. L.C. Evans, Partial Differential Equations, (AMS, Providence, Rhode Island, 2010). pp. 65-85. [ Links ]

32. M.K. Keane, A very Applied First Course in Partial Differential Equations, (Pearson Education Limited, Harlow, UK, 2001). pp. 139-163. [ Links ]

33. F.M. Henderson, Principles and Applications of Imaging Radar, (John Wiley & Sons, Inc., NY, 1998). pp. 9-47. [ Links ]

34. R. Bamler and P. Hartl, Inverse Problems 14 (1998) R1. [ Links ]

35. W.M. Brown, IEEE Transactions on Aerospace and Electronic Systems AES-3 (1967) 217. [ Links ]

36. M. Spivak, Calculus, (W.A. Benjamin Inc., NY, 1967). p. 383. [ Links ]

Appendix A

Remarks:

I) Based on 28, definition of distribution for the general case of functions h:R3C is given as follows: Let D be the set of C-functions ϕ such that, ϕ, with all its derivatives vanish at infinity as fast as ||x^||-N when ||x^||, and independently of how long is the positive integer N. Any function ϕ in the set D is said to be particularly well-behaved and such set would correspond to the space of test functions defined in 26. According to 28, a sequence of functions {hγ} (considering values γ>0) is said to be regular, if and only if limγ0+R33hγϕdV exists, for all ϕD. Thus, a distribution H is a regular sequence of functions in D given by {hγ}, where symbol

R33HϕdV

means

limγ0+R33hγϕdV. (A1)

Of course, this symbol is not a true integral, because H is a sequence. Even in the case when H is interpreted as the limit of this sequence, such limit could not be an ordinary function and that is why H is declared as a symbolic function28. In a more strict sense, a distribution, is a continuous linear functional T H defined on the space of functions D 1,26, as the symbolic notation

TH(ϕ):=R33HϕdV

suggests.

II) In this discussion, we have relaxed the hypothesis by considering D, as the set conformed by functions that are at least of C 1-class and that vanish at infinity in the next weak sense 26: function h vanishes at infinity, if and only if V[S t,h ] is finite, for all constant t > 0. Here, St,h:={x^:|h(x^)|>t} and V[S t,h ] is the volume (or Borel measure) of the set S t,h . Moreover, in this work, a distribution H is understood as the limit of a sequence {hγ}, conformed by elements in the set D, our particular set of well-behaved or test functions. However, for practical purposes, our test functions do not require to satisfy more properties, like the regularity condition imposed in 28. Conceptually, we are assuming that

H(x^):=limγ0+hγ(x^). (A2)

Therefore, when an ordinary function h is such that h(x^)=H(x^) almost everywhere, we say that h is the distribution H.

III) When considering {eik0rfγ(r)} with r=||x^|| in Eq. (25) (the simple case when y^=0^), we clearly have that functions hγ(x^):=eik0rfγ(r) are of C 1-class, as described in Sec. 4. Now, since |hγ(^)|=|fγ(r)|=fγ(r), then St,hγ=St,γ:={x^:fγ(r)>t}. Hence, given a pair of constants t,γ>0, it is not difficult to observe that V[St,γ] is finite in any case: a) For 2/γt, there is no x^R3 such that fγ(r)>t because the maximum value of fγ is achieved at r=γ/2, then St,γ= and V[]=0. b) For 1/γt<2/γ, we have St,γ=Br2(0^)\Br1(0^)¯ for some fixed values r 1 and r 2 such that 0<r1<r2γ. Here, Br1(0^)¯=Br1(0^)Br1(0^) and that implies St,γBr2(0^)Bγ(0^), which means that V[St,γ]V[Br20^]V[Bγ0^]=(4πγ3/3). c) For t<1/γ, we have St,γ=B1/t(0^)\Br1(0^)B1/t(0^) for some fixed value r 1 such that 0 < r 1 < 1/t, then V[St,γ]V[B1t0^]=(4π/3t3).

Therefore, remark III) establishes that the sequence of functions used in Eq. (25) is a sequence of well-behaved functions, a sequence in the set D defined in remark II).

IV) According to our definition of D, any function ϕ in this set has at least a continuous partial derivative ϕ/x (it could be also with respect to y, or z). Then, a sequence of C 2-class functions in D given by {hγ}, could induce a sequence {hγ/x} contained in D. If so, we can talk about wH/x=limγ0+hγ/x and H=limγ0+hγ, independently if these limits represent functions or not. Therefore, if the sequences {hγ/x} and {hγ} are such that

R33(wH/x)ϕdV=-R33H(ϕ/x)dV, (A3)

for all ϕD (by considering the symbolic representation in Eq. (A1)), then wH/x is said to be a weak derivative of H in a distributional sense. Of course, we have used the notation wH/x, to distinguish this limit (or generalized function) from the usual partial derivative of H, which is H/x every time that H is interpreted as an ordinary function (H = h).

Proposition 1 The function g(x^)=eik0r/r, where r=||x^-x^0||>0 and x^0 is a constant position, satisfies

g(x^)=eik0r(ik0r-1r3)(x^-x^0). (A4)

Moreover, such function satisfies the Helmholtz equation (2+k02)g(x^)=0 for all x^x^0.

Proposition 2 Let us consider the variables x^=(x,y,z), y^=(u,v,w), and the operator y^:=(/u,/v,/w). Thus, for any vectorial field of the form H(r)(x^-y^), with r=||x^-y^|| and H(r) as a smooth function for all r > 0, we have that

y^H(r)(x^-y^)=-H(r)(x^-y^). (A5)

In particular, function gy^(x^)=eik0||x^-y^||/||x^-y^|| satisfies that gy^(x^)=eik0r[ik0r-1r3](x^-y^), as a direct consequence of Proposition 1 (when x^0=y^). Therefore,

y^gy^(x^)=-gy^(x^), (A6)

for x^y^.

The last two propositions can be demonstrated by a careful calculation of partial derivatives and an adequate use of the chain rule.

Received: August 17, 2017; Accepted: September 11, 2017

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License