<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">JMP</journal-id><journal-title-group><journal-title>Journal of Modern Physics</journal-title></journal-title-group><issn pub-type="epub">2153-1196</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jmp.2017.810100</article-id><article-id pub-id-type="publisher-id">JMP-78934</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Sobolev Spaces, Schwartz Spaces, and a Definition of the Electromagnetic and Gravitational Coupling
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Jean-Philippe</surname><given-names>Montillet</given-names></name><xref ref-type="aff" rid="aff1"><sub>1</sub></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><label>1</label><addr-line>ESPlab, Ecole Polytechnique de Lausanne, Neuchatel, Switzerland</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>jean-philippe.montillet@epfl.ch</email></corresp></author-notes><pub-date pub-type="epub"><day>01</day><month>09</month><year>2017</year></pub-date><volume>08</volume><issue>10</issue><fpage>1700</fpage><lpage>1722</lpage><history><date date-type="received"><day>August</day>	<month>3,</month>	<year>2017</year></date><date date-type="rev-recd"><day>Accepted:</day>	<month>September</month>	<year>3,</year>	</date><date date-type="accepted"><day>September</day>	<month>6,</month>	<year>2017</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  The concept of 
  <em>multiplicity of solutions</em> was developed in [1] which is based on the theory of energy operators in the Schwartz space 
  S<sup>-</sup>(R) and some subspaces called energy spaces first defined in [2] and [3]. The main idea is to look for solutions of a given linear PDE in those subspaces. Here, this work extends previous developments in 
  S
  <sup>-</sup>(R
  <sup><em>m</em></sup>) (
  <em>m</em>
  &amp;isin;Z<sup>+</sup>) using the theory of Sobolev spaces. Furthermore, we also define the concept of 
  <em>Energy Parallax</em>, which is the inclusion of additional solutions when varying the energy of a predefined system locally by taking into account additional smaller quantities. We show that it is equivalent to take into account solutions in other energy subspaces. To illustrate the theory, one of our examples is based on the variation of Electro Magnetic (EM) energy density within the skin depth of a conductive material, leading to take into account derivatives of EM evanescent waves, particular solutions of the wave equation. The last example is the derivation of the Woodward effect [4] with the variations of the EM energy density under strict assumptions in general relativity. It finally leads to a theoretical definition of an electromagnetic and gravitational (EMG) coupling.
 
</p></abstract><kwd-group><kwd>Electromagnetism</kwd><kwd> General Relativity</kwd><kwd> Schwartz Space</kwd><kwd> Sobolev Spaces</kwd><kwd>  Multiplicity of Solutions</kwd><kwd> Energy Operators</kwd><kwd> Woodward Effect</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Overview</title><p>Teager-Kaiser energy operator was defined in [<xref ref-type="bibr" rid="scirp.78934-ref5">5</xref>] and the family of Teager- Kaiser energy operators in [<xref ref-type="bibr" rid="scirp.78934-ref6">6</xref>] . Many applications in signal processing were found over the past 25 years such as detecting transient signals [<xref ref-type="bibr" rid="scirp.78934-ref7">7</xref>] , filtering modulated signals [<xref ref-type="bibr" rid="scirp.78934-ref8">8</xref>] , image processing [<xref ref-type="bibr" rid="scirp.78934-ref9">9</xref>] . However, [<xref ref-type="bibr" rid="scirp.78934-ref2">2</xref>] and [<xref ref-type="bibr" rid="scirp.78934-ref10">10</xref>] introduced</p><p>the conjugate Teager-Kaiser energy operator and associated family ( Ψ k + ) k ∈ ℤ .</p><p>Subsequently using iterations of the Lie Bracket, [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] defined the generalized con-</p><p>jugate Teager-Kaiser energy operators ( [ [   .   ] p ] k + ) k ∈ ℤ ( p ∈ ℤ + ). To abbreviate</p><p>the notation, we sometimes use the generic name energy operator in order to refer to the conjugate Teager-Kaiser energy operators and the generalized con- jugate Teager-Kaiser energy operators. Precision is made in the denomination when it is required. Furthermore, the purpose of the energy operators and generalized energy operators was the decomposition of the successive derivatives of a finite energy function f n ( n in ℤ + − { 0,1 } ) in the Schwartz space S − ( ℝ ) . The generalized energy operators were introduced when decomposing</p><p>the successive derivatives of a finite energy function of the form ( [ [ f ] p ] 1 + ) n</p><p>( n in ℤ + − { 0,1 } ) in the Schwartz space. It then follows in [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] and [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] the definition of Energy Spaces, which are subspaces of the Schwartz Space S − ( ℝ ) associated with energy operators and generalized energy operators. This definition was used to define the concept of multiplicity of solutions in [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] (Theorem 2 and Corollary 1). The idea is to consider those energy spaces and functions associated with them when solving linear PDEs. More precisely, we look for solutions of a nominated linear PDE within those energy spaces (including the space reduced to { 0 } ). The concept was further developed using the Taylor series of the energy of a solution S − ( ℝ ) for a nominated PDE. The work was based on finding when the successive derivatives, defined through the Taylor series coefficients, are also solutions of this particular PDE (see Section 4 in [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] ).</p><p>This work first generalizes in S − ( ℝ m ) ( m ∈ ℤ + ) the theorems and lemmas established in [<xref ref-type="bibr" rid="scirp.78934-ref2">2</xref>] and [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] stated for S − ( ℝ ) using the properties of the L 2 space called here L 2 ( ℝ m ) ( m ∈ ℤ + ) together with the general property of the Schwartz space S − ( ℝ m ) ⊂ L 2 ( ℝ m ) ( m ∈ ℤ + ) [<xref ref-type="bibr" rid="scirp.78934-ref11">11</xref>] . However, this work imposes the condition of the stability by Fourier transform for any functions in S − ( ℝ ) in order to use the Sobolev space(see Appendix I, Definition I.1). Thus, in this work we consider S − ( ℝ ) together with its dual: the tempered distributions S * , − ( ℝ ) . Secondly, the energy spaces M p k ( p ∈ ℤ + , k ∈ ℤ + ) are also redefined as subspaces of S − ( ℝ m ) . Furthermore, with the definition of the Sobolev spaces, and in particular the Hilbert spaces H k ( ℝ m ) , it allows to show the inclusion M p k ⊂ H k ( ℝ m ) ⊂ L 2 ( ℝ m ) . Then, we finally redefine in ℝ m the Theorem 3 established in [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] and the concept of multiplicity of solutions.</p><p>The next section together with Appendix I are reminders about some important definitions and properties for the Sobolev spaces, the Schwartz space and the L2-norm. Section 3 deals with the generalization of the work exposed in [<xref ref-type="bibr" rid="scirp.78934-ref2">2</xref>] and [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] in S − ( ℝ m ) and the redefinition of the energy spaces. Section 4 recalls the concept of multiplicity of solutions defined in [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] and generalized in S − ( ℝ m ) with Theorem 4. The last section focuses on some applications of this theory. The first application is the wave equation and the discussion of taking into account more solutions from other energy spaces. We then define another concept called energy parallax (i.e. mathematically in Definition 4, see discussion on the physical interpretation in Appendix II) which is directly related to multiplicity of solutions. In order to illustrate this concept, a second example is the variation of energy density in the skin depth of a conductor material. The idea is to show that the variation of energy density can lead to consider multiple derivatives of evanescent waves resulting from the electromagnetic field. The last section is dedicated to the derivation of the Woodward effect [<xref ref-type="bibr" rid="scirp.78934-ref4">4</xref>] from the Hoyle-Narlikar theory [<xref ref-type="bibr" rid="scirp.78934-ref12">12</xref>] [<xref ref-type="bibr" rid="scirp.78934-ref13">13</xref>] using the EM energy density and a discussion takes place about the relationship to the presented theory of energy spaces. It leads to a theoretical definition of an Electromagnetic and Gravitational cou- pling (EMG).</p></sec><sec id="s2"><title>2. Definition of L-2 Norm and Schwartz Space</title><sec id="s2_1"><title>2.1. Notation and Symbols</title><p>In this work, several symbols are used. The set of integer numbers ℤ is some- times called only for the positive integer such as ℤ + or ℤ + m (for a space with dimension m). When the integer 0 is not included, it is explicitly mentioned such as ℤ + − { 0 } . The set of natural numbers is ℕ , with only the positive numbers defined as ℕ + . ℝ is the set of real numbers. Also, the Schwartz space is here called S − ( ℝ m ) which is the notation used in previous works such as [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] and [<xref ref-type="bibr" rid="scirp.78934-ref2">2</xref>] . Several notations describe the relationship between spaces such as intersection ( ∩ ), union ( ∪ ), inclusion ( ⊂ , inclusion without the equality ⊊ , inclusion with equality ⊆ ). Reader can refer to [<xref ref-type="bibr" rid="scirp.78934-ref14">14</xref>] or advanced mathematical textbooks for more explanations.</p></sec><sec id="s2_2"><title>2.2. L-2 Norm and Schwartz Space</title><p>With the difference in Appendix I and the generalities with the Sobolev spaces, here the analysis focuses on the L-2 norm (p equal to 2 for the L p norm). It</p><p>allows to state the Plancherel identity ∀ f ∈ L 2 ( ℝ m ) :</p><p>∫ R m f 2 d t = ∫ R m F ( f ) ( ξ ) 2 d ξ (1)</p><p>We are here interested in the functions belonging to the Schwartz space f ∈ S − ( ℝ m ) ⊂ L 2 ( ℝ m ) . The Schwartz space consists of smooth functions whose derivatives (including the function) are rapidly decreasing (e.g., the space of all bump functions [<xref ref-type="bibr" rid="scirp.78934-ref15">15</xref>] ). The Schwartz space S − ( ℝ m ) is defined as (for m ∈ [ 1,2 ] in [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] , for m ∈ ℤ + in [<xref ref-type="bibr" rid="scirp.78934-ref14">14</xref>] [<xref ref-type="bibr" rid="scirp.78934-ref16">16</xref>] ):</p><p>S − ( ℝ m ) = { f ∈ C ∞ ( ℝ m ) | f α , β &lt; ∞ , ∀ α , β ∈ ℤ + m } (2)</p><p>where α , β are multi-indices and</p><p>f α , β = sup t ∈ ℝ m | t β D α f ( t ) | (3)</p><p>Note that one can define S − ( ℝ m ) with ∀ α , β ∈ ℝ m according to [<xref ref-type="bibr" rid="scirp.78934-ref17">17</xref>] , but we decide to use ℤ + m following the development in the next sections. It is useful for the remainder of the work to remember some properties of the Schwartz functions in S − ( ℝ m ) .</p><p>Properties 1. [<xref ref-type="bibr" rid="scirp.78934-ref18">18</xref>] Some Properties of S − ( ℝ m ) .</p><p>・ If 1 ≤ p ≤ ∞ , then S − ( ℝ m ) ⊂ L p ( ℝ m )</p><p>・ S − ( ℝ m ) is a dense subspace of H k ,2 ( ℝ m ) ( k ∈ ℕ ).</p><p>・ (Stability with Fourier transform) The Fourier transform is a linear isomorphism S − ( ℝ m ) → S − ( ℝ m ) .</p><p>・ If f ∈ S − ( ℝ m ) , then f is uniformly continuous on ℝ m .</p><p>The proof of those properties are standard results with Schwartz spaces established in many harmonic analysis books (e.g., [<xref ref-type="bibr" rid="scirp.78934-ref17">17</xref>] [<xref ref-type="bibr" rid="scirp.78934-ref18">18</xref>] ).</p><p>Remark (1) Note that in [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] [<xref ref-type="bibr" rid="scirp.78934-ref2">2</xref>] [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] , the author used the general term of finite energy functions for Schwartz functions in S − ( ℝ m ) , with m restricted to [ 1,2 ] . It is a common definition in signal processing for the functions in L 2 ( ℝ m ) and generally associated with the Plancherel identity.</p><p>Remark (2) One way to interpret the property that S − ( ℝ m ) is stable by Fourier transform is:</p><p>for f ∈ S − ( ℝ m ) , k ∈ ℕ</p><p>sup ξ ∈ ℝ m { | ( 1 + ξ 2 ) k / 2 F ( f ) ( ξ ) | } &lt; ∞</p><p>↔ ∃ a ∈ ℝ ,     | ( 1 + ξ 2 ) k / 2 F ( f ) ( ξ ) | ≤ a 1 + ξ 2 (4)</p><p>Now, let us recall the definition of the Hilbert spaces H k , p ( ℝ m ) (Sobolev spaces W k , p ( ℝ m ) for p = 2 , see Appendix I, Definition I.1) from (35) and drop the sup-script p in the remainder of this work:</p><p>W k ( ℝ m ) = H k ( ℝ m )                       : = { f ∈ S * , − ( ℝ m ) | ( 1 + ξ 2 ) k / 2 F ( f ) ∈ L 2 ( ℝ m ) } (5)</p><p>Note that S * , − ( ℝ m ) is the space of tempered distributions, dual of S − ( ℝ m ) via the Fourier transform. A function belongs to L 2 ( ℝ m ) if and only if its Fourier transform belongs to L 2 ( ℝ m ) and the Fourier transform preserves the L<sup>2</sup>-norm. As a result, the Fourier transform provides a simple way to define L<sup>2</sup>- Sobolev spaces on ℝ m (including ones of fractional and negative order m [<xref ref-type="bibr" rid="scirp.78934-ref18">18</xref>] ). Finally, the stability via Fourier transform is the key for S − ( ℝ m ) ⊊ H k ( ℝ m ) .</p><p>Remark (3) Following the remark (Remark 3.4 in [<xref ref-type="bibr" rid="scirp.78934-ref19">19</xref>] ) and the general properties of the Fourier transform, one can state the equivalence relationship in L 2 ( ℝ m )</p><p>f ∈ H k ( ℝ m ) ↔ D α f ∈ L 2 ( ℝ m )   ∀ α ≤ k</p><p>↔ F ( D α f ) ∈ L 2 ( ℝ m )   ∀ α ≤ k</p><p>↔ ξ α F ( f ) ∈ L 2 ( ℝ m )   ∀ α ≤ k</p><p>↔ ( 1 + ξ 2 ) α / 2 F ( D α f ) ∈ L 2 ( ℝ m )   ∀ α ≤ k (6)</p><p>Using the definition of H k ( ℝ m ) and the properties of the Fourier transform, it is also possible to show that for k &gt; k ′ , H k ( ℝ m ) ⊂ H k ′ ( ℝ m ) [<xref ref-type="bibr" rid="scirp.78934-ref20">20</xref>] , and the relationship H 0 ( ℝ m ) = L 2 ( ℝ m ) . It is also possible to define</p><p>H ∞ ( ℝ m ) = ∩ k ∈ ℕ H k ( ℝ m ) with S − ( ℝ m ) ⊂ H ∞ ( ℝ m ) , and to extend this equality to k ∈ ℝ following [<xref ref-type="bibr" rid="scirp.78934-ref19">19</xref>] .</p></sec></sec><sec id="s3"><title>3. On Some Subsets of Schwartz Spaces: Energy Spaces</title><p>This section first recalls generalities on the Teager-Kaiser energy operator and its conjugate operator with the application to decompose Schwartz functions from the work developed in [<xref ref-type="bibr" rid="scirp.78934-ref2">2</xref>] and [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] . We call in this work Energy operators the families of operators based on the Teager-Kaiser energy operator. The definitions and theorems are here stated for the Schwartz space S − ( ℝ m ) ( m ∈ ℕ ) whereas the preliminary work in [<xref ref-type="bibr" rid="scirp.78934-ref2">2</xref>] and [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] stated the definitions and main theorems for m ∈ [ 1,2 ] . For m = 2 in Section 6 in [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] , a discussion takes place during the application of the theory to linear partial differential equations. Secondly, the energy spaces defined in [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] and [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] are here generalized on S − ( ℝ m ) with novel relationships with Sobolev spaces H k ( ℝ m ) ( k ∈ ℕ ).</p><sec id="s3_1"><title>3.1. Definition and Properties of the Energy Operators in S − ( ℝ m )</title><p>Let us call the set F ( S − ( ℝ m ) , S − ( ℝ m ) ) all Schwartz functions (or operators) defined such as γ : S − ( ℝ m ) → S − ( ℝ m ) . For f ∈ S − ( ℝ m ) , let us define ∂ i k f</p><p>( k ∈ ℤ , i ∈ [ 1, ⋯ , m ] ), with f defined with the vector parameter T = [ t 1 , t 2 , ⋯ , t m ] ∈ ℝ m such as</p><p>{ ∂ i k f = ∂ k f ∂ t i k ,   ∀ i ∈ [ 1, ⋯ , m ] ,   ∀ k ∈ ℤ + − { 0 } ∂ i k f = ∫ − ∞ t i ( ⋯ ( ∫ − ∞ τ 1 f ( t 1 , t 2 , ⋯ , τ 1 , t k + 1 , t m ) d τ 1 ) ⋯ ) d τ k ,   ∀ i ∈ [ 1, ⋯ , m ] ,   ∀ k ∈ ℤ − − { 0 } ∂ i 0 f = f ,   ∀ i ∈ [ 1, ⋯ , m ] (7)</p><p>Combining multiple integrals and derivatives justify the use of the Schwartz space S − ( ℝ m ) and echoes the choice made previously in [<xref ref-type="bibr" rid="scirp.78934-ref2">2</xref>] (see equation (10)). The definitions and results given in [<xref ref-type="bibr" rid="scirp.78934-ref2">2</xref>] and [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] in the case S − ( ℝ ) are now formulated for S − ( ℝ m ) . Section 2 in [<xref ref-type="bibr" rid="scirp.78934-ref2">2</xref>] and Section 4 in [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] defined the energy operators Ψ k + , Ψ k − ( k in ℤ ) and the generalized energy operators</p><p>[ [   .   ] p ] k + and [ [   .   ] p ] k − ( p in ℤ + ). Following [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] , let us define the energy</p><p>operators with multi-index derivative in (7):</p><p>Ψ k + ( . ) = ∑ i = 1 m     ∂ i 1 . ∂ i k − 1 . + ∂ i 0 . ∂ i k . Ψ k + ( . ) = ∑ i = 1 m     ψ k , i + ( . ) [ . , . ] k + = Ψ k + ( . ) [ . , . ] k , i + = ψ k , i + ( . ) (8)</p><p>Further more, we also use the short notation [ . , . ] k + = [   .   ] k + in the remainder of this work. Note that Ψ k − is the conjugate operator of Ψ k + and ψ k , i − respec- tively to ψ k , i + .</p><p>Remark (4) The families of (generalized) energy operators ( [ [   .   ] p ] k + ) k ∈ ℤ and ( [ [   .   ] p ] k − ) k ∈ ℤ ( p in ℤ + ) are also called families of differential energy operator</p><p>(DEO) [<xref ref-type="bibr" rid="scirp.78934-ref2">2</xref>] [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] .</p><p>Furthermore, [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] defined the generalized energy operators [ [   .   ] 1 ] k + and [ [   .   ] 1 ] k − ( k ∈ ℤ ):</p><p>[ [ . , . ] k , i + , [ . , . ] k , i + ] k , i + = ∂ i 1 ψ k , i + ( . ) ∂ i k − 1 ψ k , i + ( . ) + ∂ i 0 ψ k , i + ( . ) ∂ i k ψ k , i + ( . ) [ [ . , . ] k , i + , [ . , . ] k , i + ] k , i + = ∂ i 1 [ [   .   ] 0 ] k , i + ∂ i k − 1 [ [   .   ] 0 ] k , i +                                                   + ∂ i 0 [ [   .   ] 0 ] k , i + ∂ i k [ [   .   ] 0 ] k , i + [ [ . , . ] k + , [ . , . ] k + ] k + = ∑ i = 1 m [ [ . , . ] k , i + , [ . , . ] k , i + ] k , i + [ [ . , . ] k + , [ . , . ] k + ] k + = ∑ i = 1 m [ [   .   ] 1 ] k , i + = [ [   .   ] 1 ] k + (9)</p><p>By iterating the bracket [   .   ] , [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] defined the generalized operator [ [   .   ] p ] k , i − and the conjugate [ [   .   ] p ] k , i + with p in ℤ + . Note that [ [ f ] p ] 1, i − = 0 ∀ p in ℤ + and i in ℤ .</p><p>Now, the derivative chain rule property and bilinearity of the energy operators and generalized operators (for i in [ 1,2 ] ) are shown respectively in [<xref ref-type="bibr" rid="scirp.78934-ref2">2</xref>] , Section 2 and [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] , Proposition 3. The generalisation of this property to i in [ 1, ⋯ , m ] for</p><p>the operators ψ k , i + ( . ) , ψ k , i − ( . ) , [ [   .   ] p ] k , i − and [ [   .   ] p ] k , i + ( k ∈ ℤ , p ∈ ℤ + ) is</p><p>trivial due to the linearity of the derivatives and integrals when defining ∂ i k in (7). Due to the linearity of the sum, the bilinearity property is also generalized to</p><p>Ψ k + ( . ) , Ψ k − ( . ) , [ [   .   ] 1 ] k + and [ [   .   ] 1 ] k − ( k ∈ ℤ , p ∈ ℤ + ).</p><p>Definition 1. [<xref ref-type="bibr" rid="scirp.78934-ref2">2</xref>] ∀ f in S − ( ℝ m ) , ∀ v ∈ ℤ + − { 0 } , ∀ n ∈ ℤ + and n &gt; 1 , the family of operators ( G k ) k ∈ ℤ (with ( G k ) k ∈ ℤ ⊆ F ( S − ( ℝ m ) , S − ( ℝ m ) ) ) decom- poses ∂ i v f n in ℝ m ( i ∈ [ 1, ⋯ , m ] ), if it exists ( N j ) j ∈ ℤ + ∪ { 0 } ⊆ ℤ + , ( C l ) l = − N j N j ⊆ ℝ , and it exists ( α j ) and r in ℤ + ∪ { 0 } (with r &lt; v ) such as ∂ i v f n = ∑ j = 0 v − 1 ( v − 1 j ) ∂ i v − 1 − j f n − r ∑ u = − N j N j C u G u ( ∂ i α u f ) .</p><p>In addition, one has to define s − ( ℝ m ) as:</p><p>s − ( ℝ m ) = { f ∈ S − ( ℝ m ) | f ∉ ( ∪ k ∈ ℤ K e r ( Ψ k + ) ) ∪ ( ∪ k ∈ ℤ − { 1 } K e r ( Ψ k − ) ) } (10)</p><p>or with the energy operators ψ k , i + and ψ k , i − defined in (8)</p><p>s − ( ℝ m ) = { f ∈ S − ( ℝ m ) | f ∉ ∪ i ∈ [ 1, ⋯ , m ] ( ∪ k ∈ ℤ K e r ( ψ k , i + ( f ) ) ∪ ( ∪ k ∈ ℤ − { 1 } K e r ( ψ k , i − ( f ) ) ) ) } (11)</p><p>K e r ( . ) is the notation for the kernel associated here with the operators Ψ k + , Ψ k − , ψ k , i + and ψ k , i − ( k in ℤ ) (see [<xref ref-type="bibr" rid="scirp.78934-ref2">2</xref>] , Properties 1 and 2). By definition, one can state that s − ( ℝ m ) ⊊ S − ( ℝ m ) . Following Definition 1, the uniqueness of the decomposition in s − ( ℝ m ) with the families of differential operators can be stated as:</p><p>Definition 2. [<xref ref-type="bibr" rid="scirp.78934-ref2">2</xref>] ∀ f in s − ( ℝ m ) , ∀ v ∈ ℤ + − { 0 } , ∀ n ∈ ℤ + and n &gt; 1 , the</p><p>families of operators ( Ψ k + ) k ∈ ℤ and ( Ψ k − ) k ∈ ℤ ( ( Ψ k + ) k ∈ ℤ and ( Ψ k − ) k ∈ ℤ ⊆ F ( s − ( ℝ m ) , S − ( ℝ m ) ) decompose uniquely ∂ i v f n in ℝ m , if for any</p><p>family of operators ( S k ) k ∈ ℤ ⊆ F ( S − ( ℝ m ) , S − ( ℝ m ) ) decomposing ∂ i v f n in ℝ m , there exists a unique couple ( β 1 , β 2 ) in ℝ 2 m such as:</p><p>S k ( f ) = β 1 Ψ k + ( f ) + β 2 Ψ k − ( f ) ,   ∀ k ∈ ℤ (12)</p><p>Two important results shown in [<xref ref-type="bibr" rid="scirp.78934-ref2">2</xref>] are:</p><p>Lemma 1 For f in S − ( ℝ m ) , the family of DEO Ψ k + ( k ∈ ℤ ) decomposes ∂ i v f n , ∀ v ∈ ℤ + − { 0 } , n ∈ ℤ + − { 0,1 } and i ∈ [ 1, ⋯ , m ] .</p><p>Theorem 1. For f in s − ( ℝ m ) , the families of DEO Ψ k + and Ψ k − ( k ∈ ℤ ) decompose uniquely ∂ i v f n , ∀ v ∈ ℤ + − { 0 } , n ∈ ℤ + − { 0,1 } and i ∈ [ 1, ⋯ , m ] .</p><p>The Lemma 1 and Theorem 1 were then extended in [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] to the family of generalized operator with :</p><p>Lemma 2. For f in S p − ( ℝ m ) , p in ℤ + , the families of generalized energy operators [ [   .   ] p ] k + ( k ∈ ℤ ) decompose ∂ i v ( [ [ f ] p − 1 ] 1 + ) n ∀ v ∈ ℤ + − { 0 } , n ∈ ℤ + − { 0,1 } and i ∈ [ 1, ⋯ , m ] .</p><p>Theorem 2. For f in s p − ( ℝ m ) , for p in ℤ + , the families of generalized operators [ [   .   ] p ] k + and [ [   .   ] p ] k − ( k ∈ ℤ ) decompose uniquely ∂ i v ( [ [ f ] p − 1 ] 1 + ) n ∀ v ∈ ℤ + − { 0 } , n ∈ ℤ + − { 0,1 } and i ∈ [ 1, ⋯ , m ] .</p><p>S p − ( ℝ m ) and s p − ( ℝ m ) ( p in ℤ + ) are energy spaces in S − ( ℝ m ) defined in the next section.</p><p>Remark (5) One can extend the Theorem 1, Theorem 2, Lemma 1 and Lemma 2 for f n with n in ℤ following previous discussions in [<xref ref-type="bibr" rid="scirp.78934-ref2">2</xref>] (Section 3, p.74) and [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] (Section 4). n is here restricted to ℤ + − { 0,1 } in order to easy the whole mathematical development.</p></sec><sec id="s3_2"><title>3.2. Energy Spaces in S − ( ℝ m )</title><p>Let us introduce the energy spaces and some properties.</p><p>Definition 3. ( [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] , Definition 3) The energy space E p ⊊ S − ( ℝ m ) , with p in ℤ + , is equal to E p = ∪ v ∈ ℤ + ∪ { 0 } M p v .</p><p>With M p v ⊊ S − ( ℝ m ) for v in ℤ + defined as</p><p>M p v = { g ∈ S − ( ℝ m ) | g = ∂ i k ( [ [ f ] p ] 1 + ) n , [ [ f ] p ] 1 + ∈ S − ( ℝ m ) , k ∈ ℤ + ,                 ∀ k ≤ v , n ∈ ℤ + − { 0 } , i ∈ [ 1, ⋯ , m ] } (13)</p><p>The energy spaces, S p − ( ℝ m ) and s p − ( ℝ m ) ( p ∈ ℤ + ) , cited in Lemma 2 and Theorem 2 are defined:</p><p>S p − ( ℝ m ) = { E p = ∪ i ∈ ℤ + ∪ { 0 } M p i } s p − ( ℝ m ) = { f ∈ S p − ( ℝ m ) | f ∉ ∪ i ∈ [ 1, ⋯ , m ] ( ∪ k ∈ ℤ K e r ( [ [ f ] p ] k , i + )                                 ∪ ( ∪ k ∈ ℤ − { 1 } K e r ( [ [ f ] p ] k , i − ) ) ) } (14)</p><p>Remark (6) Definition 3 does not follow completely Definition 3 in [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] , because the energy space M p v is defined here ∀ k ≤ v , and only for k = v in [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] .</p><p>Remark (7) In the previous definition, M p ∞ = { 0 } ( ∀ p ∈ ℤ + ). Also, M p ∞ ⊂ S p − ( ℝ m ) , whereas M p ∞ ⊂ S p − ( ℝ m ) in [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] and [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] . The inclusion does not change Lemma 2 and Theorem 2 (i.e. M p ∞ ⊂ s p − ( ℝ m ) ). The justification of not</p><p>including this space was only based on the applications of the theory in [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] and [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] which is not justified in this work.</p><p>We can now state some properties associated with the energy spaces on S − ( ℝ m ) .</p><p>Properties 2. ∀ v in ℤ + , and in particular v 1 , v 2 in ℤ + (with v 1 &lt; v 2 ), p in ℤ + , we have the following inclusions:</p><p>・ M p v ⊊ H v ( ℝ m )</p><p>・ M p v 2 ⊊ M p v 1</p><p>・ E p = ∪ v ∈ ℤ + ∪ { 0 }   M p v ⊊ H 0 ( ℝ m )</p><p>Proof.</p><p>1) Let us recall the definition of the Hilbert space on ℝ m according to Appendix I, Definition I.1 and Definition 1.</p><p>H v ( ℝ m ) = { f ∈ L 2 ( ℝ m ) | D α f ∈ L 2 ( ℝ m ) , ∀ α ≤ v } (15)</p><p>Looking at the definition of the energy space M p v and H v ( ℝ m ) , one can notice the similitude. However , the multi-index derivative D α ( [<xref ref-type="bibr" rid="scirp.78934-ref21">21</xref>] , chap. 1.1)</p><p>contains also the cross-derivatives (e.g., ∂ 2 ∂ t 1 ∂ t 2 ), whereas there are no cross-</p><p>derivatives in the definition of ∂ i v at the beginning of Appendix I. Thus, the energy spaces M p v ( p ∈ ℤ + , v ∈ ℤ + − { 0 } ) is defined without the cross- derivatives. In addition with Properties 1, S − ( ℝ m ) ⊊ L 2 ( ℝ m ) . Thus, by definition we have the relationship M p v ⊊ H v ( ℝ m ) :</p><p>2) With Remark (3), we know that for v 1 &lt; v 2 , H v 2 ( ℝ m ) ⊂ H v 1 ( ℝ m ) . Now, with 1), M p v 1 ⊊ H v 1 ( ℝ m ) and M p v 2 ⊊ H v 2 ( ℝ m ) . Now by definition of M p v 1 and M p v 2 , M p v 1 H v 2 ( ℝ m ) = M p v 2 . Finally, M p v 2 ⊊ M p v 1 .</p><p>3) From Remark (3), H 0 ( ℝ m ) = L 2 ( ℝ m ) , S − ( ℝ m ) ⊂ L 2 ( ℝ m ) and (by defi-</p><p>nition of the energy space) E p = ∪ v ∈ ℤ + − { 0 } M p v ⊊ S − ( ℝ m ) ( p ∈ ℤ + ). Thus, E p = ∪ v ∈ ℤ + − { 0 } M p v ⊊ H 0 ( ℝ m ) ( p ∈ ℤ + ).</p><p>Furthermore, Appendix III discusses the relationship between the subspaces M p v and M p − 1 v ( p ∈ ℤ + ). Finally, because we are studying functions and operators in subspaces of S − ( ℝ m ) with S − ( � m ) ⊂ L 2 ( � m ) , one need to extend Proposition 1 in [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] and [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] .</p><p>Proposition 1. If for n ∈ ℤ + , f n ∈ S − ( ℝ m ) and analytic; for any ( p i , q i ) ∈ ℝ 2 and τ i ∈ [ q i , p i ] ( ∀ i ∈ [ 1, ⋯ , m ] ), and E ( f n ) is analytic, where</p><p>E ( f n ( τ i ) ) = ∫ q i τ i f n ( t i ) 2 d t i &lt; ∞ (16)</p><p>then</p><p>E ( f n ( p i ) ) = E ( f n ( q i ) ) + ∑ k = 0 ∞     ∂ t i k ( f n ( q i ) ) 2 ( p i − q i ) k k ! &lt; ∞ (17)</p><p>is a convergent series.</p><p>Proof. The proof of Proposition 1 for i equal 1 is given in [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] (p.4). The extension of the proof for the case i equal m is straightforward with the general definition for any ( p i , q i ) ∈ ℝ 2 and τ i ∈ [ q i , p i ] ( ∀ i ∈ [ 1, ⋯ , m ] ).</p></sec></sec><sec id="s4"><title>4. Multiplicity of the Solutions in S − ( ℝ m )</title><p>To recall [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] , a possible application of the theory of the energy operators is to look at solutions of a given partial differential equation for solutions in S − ( ℝ m ) of the form ∂ i v ( f n ) . Instead of solving the equation for specific values (e.g., boundary conditions), the work in [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] ( [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] , Theorem 1 and corollary)defines the concept of multiplicity of solutions in S − ( ℝ m ) ( m ∈ [ 1,2 ] ) such as the study of the multiple solutions of a PDE based on the definition of the energy spaces E p ( p ∈ ℤ + ). One way to understand this concept, is to study the convergence of the development in Taylor series of the energy function associated to a nominated energy space. It was shown in [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] that taking into account additional terms of the Taylor series leads to define additional solutions of the wave equation (see Section 4 [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] ). In this section, we extend this concept to S − ( ℝ m ) ( m ∈ ℕ 0 ) and we reformulate the results from [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] for the solutions in the subspaces M p v ⊂ E p ⊂ S − ( ℝ m ) ⊂ L 2 ( ℝ m ) ( p ∈ ℤ + , v ∈ ℤ + ).</p><p>Let us define any PDEs of the form:</p><p>{ ∑ j ∈ ℤ + ∑ i ∈ [ 0 , ⋯ , m ] a i j ∂ i v j g = 0, ∀ g ∈ A ( ℝ m ) ⊆ S − ( ℝ m ) , ∀ a i j ∈ ℝ ,   v j ∈ ℤ + (18)</p><p>Thus, all the solutions are here defined in A ( ℝ m ) ⊆ S − ( ℝ m ) . Now, we are interested in the solutions which can be defined on the energy spaces E p</p><p>( p ∈ ℤ + ). In other words, A ( ℝ m ) ∩ p ∈ ℤ + E p ≠ { ∅ } . In particular, we choose the solution g = 0 ∈ A ( ℝ m ) ∩ p ∈ ℤ + E p . Furthermore, one can define g ∈ A ( ℝ m ) ∩ p ∈ ℤ + E p , such as ∃ v ∈ ℤ + for g ∈ A ( ℝ m ) ∩ p ∈ ℤ + M p v . In other words, ∃ f ∈ S − ( ℝ m ) and n ∈ ℤ + − { 0 } , such as g = ( [ [ f ] p ] 1 + ) n . Now, one can then state a general theorem of multiplicty of solutions based on [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] . It follows:</p><p>Theorem 3. (Multiplicity of Solutions in ℝ m ) If A ( ℝ m ) ⊆ S − ( ℝ m ) is a subspace of all the solutions of a nominated linear PDE. For p ∈ ℤ + , g is in E p . Then, g is a solution for this linear PDE if and only if:</p><p>1) (General condition to be a solution) A ( ℝ m ) ∩ E p ≠ { ∅ } .</p><p>2) (Solutions in S − ( ℝ m ) ) g ∈ A ( ℝ m ) ∩ E p , ∃ m ∈ ℝ + such as m = sup ( E ( g ) ) .</p><p>3) (Multiplicity of the solutions) If g ∈ M p v ( v ∈ ℤ + ), ∃ f ∈ S − ( ℝ m ) and</p><p>n ∈ ℤ + − { 0 } , such as g = ∂ i v ( [ [ f ] p ] 1 + ) n ( i ∈ [ 0, ⋯ , m ] ) and ∀ k ≥ v , k ∈ ℤ + , ∂ i k ( [ [ f ] p ] 1 + ) n ∈ A ( ℝ m ) ∩ E p .</p><p>4) (Superposition of solutions and energy conservation ) If F ∈ A ( ℝ m ) ∩ E p , with F = ∑ k ∈ ℤ + , ∀ k ≥ v ∂ i k ( [ [ f ] p ] 1 + ) n such as ∂ i k ( [ [ f ] p ] 1 + ) n ∈ M p k ( i ∈ [ 0, ⋯ , m ] ), then E ( F ) &lt; ∞ .</p><p>Proof. The proof is the generalization of what was already written in [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] (see Theorem 2 in [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] ) for the case m equal 1. Here is the generalization to m.</p><p>1) This is the definition of a solution for a nominated PDE with solutions in A ( ℝ m ) and in the energy space E p .</p><p>2) g ∈ A ( ℝ m ) ∩     E p ⊂ L 2 ( ℝ m ) , thus E ( g ) &lt; ∞ . With Proposition 1, it</p><p>means that for any ( p i , q i ) ∈ ℝ 2 and τ i ∈ [ q i , p i ] ( ∀ i ∈ [ 1, ⋯ , m ] )</p><p>E ( g ( τ i ) ) = ∫ q i τ i g ( t i ) 2 d t i &lt; ∞ (19)</p><p>Thus, following [<xref ref-type="bibr" rid="scirp.78934-ref18">18</xref>] , one can define m i ∈ ℝ such as m i = sup τ i ∈ [ q i , p i ] E ( g ( τ i ) ) and then we define m = max i ∈ [ 1 , ⋯ , m ] m i . With our notation, it is equivalent to write m = sup ( E ( g ) ) .</p><p>3) It is sufficient to show that for v ∈ ℤ + , ∀ k ≥ v , A ( ℝ m ) ∩ M p k ≠ { ∅ } .</p><p>Now, with the definition A ( ℝ m ) ∩ E p ≠ { ∅ } , and A ( ℝ m ) ∩ M p k ≠ { ∅ } . In addition, M p ∞ = { 0 } , M p ∞ ⊂ M p k ( ∀ k ≥ v ) and 0 ∈ A ( ℝ m ) ∩ M p k . The interest</p><p>of this statement is the function ∂ i v h ∈ S − ( ℝ m ) such as ∃ k ∈ ℤ + with k ≥ v</p><p>and ∂ i k h = 0 . In particular, if we introduce a numerical approximation in order to get the condition ∂ i k h ~ 0 . In other words,</p><p>∂ i k h ~ 0 ↔ { ∃ ϵ ∈ ℝ + , ϵ ≪ 1 , ∀ k ∈ ℤ + , k &gt; 0 ,   such   as   | ∂ i k h | ≤ ϵ } (20)</p><p>In some examples in Section 4 in [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] and Section 6 in [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] , it is shown that the evanescent waves when solving the wave equation for specific solutions, is a particular example of those functions.</p><p>4) The proof follows [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] (Theorem 2). This statement is to guarantee that there is a finite sum of energy with the superposition of multiple solutions. Thus with the development in statement (2.), one can use the Minkowski inequality (e.g, [<xref ref-type="bibr" rid="scirp.78934-ref22">22</xref>] , Theorem 202) for τ i in [ p i , q i ] ( ∀ i ∈ [ 1, ⋯ , m ] )</p><p>E ( F ( τ i ) ) = ∫ p i τ i | ∑ k ∈ ℤ + , ∀ k ≥ v ∂ i k ( [ [ f ( t i ) ] p ] 1 + ) n | 2 d t i</p><p>( E ( F ( τ i ) ) ) 0.5 ≤ ∑ k ∈ ℤ + , ∀ k ≥ v ( ∫ p i τ i | ∂ i k ( [ [ f ( t i ) ] p ] 1 + ) n | 2 d t i ) 0.5</p><p>( E ( F ( τ i ) ) ) 0.5 ≤ ∑ k ∈ ℤ + , ∀ k ≥ v ( E ( ∂ i k ( [ [ f ( τ i ) ] p ] 1 + ) n ) ) 0.5</p><p>( E ( F ( τ i ) ) ) 0.5 ≤ ∑ k ∈ ℤ + , ∀ k ≥ v m k 0.5 (21)</p><p>with m k = sup τ i ∈ [ q i , p i ] ( E ( ∂ i k ( [ [ f ( τ i ) ] p ] 1 + ) n ) ) . Thus, (4.) stands if and only if ∑ k ∈ ℤ + , ∀ k ≥ v m k 0.5 &lt; ∞ . As ∀ k ∈ ℤ + , ∀ k ≥ v , m k is in ℝ + , it then exists M = sup ∑ k ∈ ℤ + , ∀ k ≥ v m k 0.5 . One possibility is ∃ k o in ℤ + such as ∀ k &gt; k o , then</p><p>m k = 0 .</p></sec><sec id="s5"><title>5. Some Applications</title><p>This section focuses on the application of the energy space theory. The first section is the study of the concept of multiplicity of solutions with a simple mathematical example using the wave equation. Then, the second section is discussing the application of this concept within the Woodward effect [<xref ref-type="bibr" rid="scirp.78934-ref4">4</xref>] [<xref ref-type="bibr" rid="scirp.78934-ref13">13</xref>]</p><sec id="s5_1"><title>5.1. Energy Variation and Wave Equation</title><p>As a simple case of linear PDE, the wave equation with the particular solutions of the form of evanescent waves, was already discussed in Section 6 of [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] and [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] . However, it is an interesting example to apply and understand the concept of multiplicity stated in Theorem 3. From [<xref ref-type="bibr" rid="scirp.78934-ref23">23</xref>] , the wave equation can be formulated in ℝ 2 (with t and r the time and space variables):</p><p>{ ∂ r 2 g ( r , t ) − 1 c 2 ∂ t 2 g ( r , t ) = 0 , t ∈ [ 0 , T ] ,   r ∈ [ r 1 , r 2 ] ,         ( r 1 , r 2 , T ) ∈ ℝ 3 ,     r 1 &lt; r 2 t 0 ∈ [ 0 , T ] ,   r 0 ∈ [ r 1 , r 2 ] (22)</p><p>c is the speed of light. Note that the values of t and r are restricted to some interval, because it is conventional to solve the equation for a restricted time interval in ℝ + and a specific region in space. According to the previous section, we are here interested in the solutions in the energy (sub)space M p k , of the kind</p><p>g ( r , t ) = ∂ t k ( [ [ f ] p ] 1 + ) n ( r , t ) ( n in ℤ + − { 0 } , p in ℤ + , k in ℤ + ). Furthermore, the relationship M p k ⊂ S − ( ℝ 2 ) ⊂ L ( ℝ 2 ) imposes that the solu-</p><p>tions should be finite energy functions, decaying for large values of r and t. It was previously underlined in [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] and [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] that planar waves should be rejected, because this type of solutions does not belong to L ( ℝ 2 ) . However, evanescent waves are a type of solutions included in S − ( ℝ 2 ) and considered in this work. They are here defined such as:</p><p>{ f ( r , t ) = R e a l { A e x p ( u 2 r ) e x p ( i ( ω t − u 1 r ) ) } , t ∈ [ 0, T ] ,   r ∈ [ r 1 , r 2 ] ,   ( r 1 , r 2 ) ∈ ℝ 2 ,   r 1 &lt; r 2 (23</p><p>i 2 = − 1 , u 1 and u 2 are the wave numbers, ω is the angular frequency and A is the amplitude of this wave [<xref ref-type="bibr" rid="scirp.78934-ref23">23</xref>] . Assuming ω and ( u 1 , u 2 ) known, one can add some boundary conditions in order to estimate u 1 , u 2 and A . Furthermore, a traveling wave solution of (19) should satisfy the dispersion relationship between u 1 , u 2 and ω [<xref ref-type="bibr" rid="scirp.78934-ref23">23</xref>] . However, our interest is just the general form assuming that all the parameters are known. For p = 0 , the type of solutions in M 0 k are:</p><p>{ ∂ t k f n ( r 0 , t ) = ( i ω n ) k f n ( r 0 , t ) , ∂ r k f n ( r , t 0 ) = ( n ( u 2 − i u 1 ) ) k f n ( r , t 0 ) , t ∈ [ 0 , T ] ,   r ∈ [ r 1 , r 2 ] ,   ( r 1 , r 2 , T ) ∈ ℝ 3 ,   r 1 &lt; r 2 t 0 ∈ [ 0 , T ] ,   r 0 ∈ [ r 1 , r 2 ] ,   n ∈ ℤ + − { 0 } ,   k ∈ ℤ + − { 0 } (24)</p><p>In M 1 k , one can then write the type of solutions</p><p>∂ t k Ψ 1 , t + ( f ) ( r 0 , t ) = ∂ t k ( [ [ f ( r 0 , t ) ] 1 ] 1 , t + ) ∂ t k Ψ 1 , t + ( f ) ( r 0 , t ) = R e a l { ( i 2 k ω ) f 2 ( r 0 , t ) } ∂ r k Ψ 1 , r + ( f ) ( r , t 0 ) = R e a l { ( 2 k ( u 2 − i u 1 ) ) f 2 ( r , t 0 ) } t ∈ [ 0 , T ] ,   r ∈ [ r 1 , r 2 ] ,   ( r 1 , r 2 , T ) ∈ ℝ 3 ,   r 1 &lt; r 2 t 0 ∈ [ 0 , T ] ,   r 0 ∈ [ r 1 , r 2 ] ,   n ∈ ℤ + − { 0 } ,   k ∈ ℤ + − { 0 } (25)</p><p>Let us consider the form of solutions which propagates in a closed cavity (e.g., closed wave guide [<xref ref-type="bibr" rid="scirp.78934-ref23">23</xref>] ). One possible solution is the evanescent wave described in (20). Now, if f and E ( f ) are analytic in ℝ 2 , with Proprsition 1 we can assume that f is finite energy (and more generally in S − ( ℝ 2 ) ) with a wise choice on the parameters A , u 1 , u 2 and ω . One can estimate the difference of energy in time over d t inside the cavity at a specific location r 0 ( r 0 in [ r 1 , r 2 ] ) such as</p><p>E ( f ( r 0 , T ) ) = ∫ 0 T ( f ( r 0 , h ) ) 2 d h &lt; ∞</p><p>E ( f ( r 0 , T + d t ) ) = E ( f ( r 0 , T ) ) + ∑ k = 0 ∞     ∂ t k ( f 2 ( r 0 , T ) ) ( d t ) k k ! &lt; ∞</p><p>E ( f ( r 0 , T + d t ) ) = E ( f ( r 0 , T ) ) + f 2 ( r 0 , T ) d t + ∑ k = 1 ∞     ∂ t k − 1 Ψ 1 , t + ( f ) ( r 0 , T )</p><p>E ( f ( r 0 , T + d t ) ) ≃ E ( f ( r 0 , T ) ) + f 2 ( r 0 , T ) d t (26)</p><p>Here the symbol ‘ ≃ ‘ means that</p><p>∃ ϵ ∈ ℝ + ,   ϵ ≪ 1 ,   ∀ k ∈ ℤ + ,   k &gt; 0 |   | ∂ t k − 1 ( Ψ 1 , t + ( f ) ( r 0 , T ) ) | &lt; ϵ | f 2 ( r 0 , T ) | (27)</p><p>Now, let us do a hypothesis that E ( f ( r 0 , T + d t ) ) increases significantly over d t modifying the approximation in (24)</p><p>∃ ϵ ∈ ℝ + ,   ϵ ≪ 1 ,   ∀ k ∈ ℤ + ,   k &gt; 1 |   | ∂ t k − 1 Ψ 1 , t + ( f ) ( r 0 , T ) | &lt; ϵ | Ψ 1 , t + ( f ) ( r 0 , T ) | (28)</p><p>and then,</p><p>E ( f ( r 0 , T + d t ) ) ≃ E ( f ( r 0 , T ) ) + f 2 ( r 0 , T ) d t + Ψ 1, t + ( f ) ( r 0 , T ) d t 2 2 (29)</p><p>To recall that f 2 ( r , t ) ∈ M 0 0 , ∂ t f 2 ( r , t ) ∈ M 0 1 and Ψ 1, t + ( f ) ( r , t ) ∈ M 1 0 , and using Theorem 3, one can take into account solutions in those subspaces. The multiplicity of the solutions due to the variation of energy can be formulated as an approximation for taking into account additional solutions produced by the wave equation.</p><p>Remark (8): In [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] , the general idea was to look for the solutions of linear PDEs in S − ( ℝ ) associated with energy subspaces s p − ( ℝ ) ( p ∈ ℤ + ) in order to apply Theorem 1 in [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] , which is here generalized in Theorem 3 for ℝ m ( m ∈ ℕ + ). The purpose was to find the subspaces reduced to { 0 } when studying the convergence of the Taylor series of the energy functions. However, the redefinition of the energy subspaces M p k within the Sobolev spaces defined in Section 3 allows us to look for solutions in S p − ( ℝ m ) in order to use Lemma 2. Because of the inclusion of the energy spaces shown in Properties 2 using the Sobolev embedding (e.g., Theorem I.1 in Appendix I) such as M p k + 1 k ⊂ M p k ( k ∈ ℤ + , p ∈ ℤ + ), M p ∞ = { 0 } ⊂ M p k .</p><p>Definition 4. (Energy Parallax) Considering a linear PDE with some solutions in A ( ℝ m ) such as A ( ℝ m ) ∩ S − ( ℝ m ) ≠ { ∅ } . Furthermore, if it exists p and v ∈ ℤ + such as A ( ℝ m ) ∩ M p v ≠ { ∅ } , then we associate the energy E ( f ) for f ∈ A ( ℝ m ) ∩ M p v , such as one can estimate the variation d E ( f ) = ( E ( f ( q + d q ) ) − E ( f ( q ) ) ) d q over an elementary quantity d q (e.g., space or time). If d E ( f ) is not negligible ( ∃ ϵ ∈ ℝ m such as ϵ ≫ 1 and d E ( f ) &gt; ϵ ), then one can consider additional solutions in A ( ℝ m ) ∩ M p v + 1 .</p></sec><sec id="s5_2"><title>5.2. Variation of EM Energy Density and the Woodward Effect</title><p>In this section, the theory of energy space is applied to the possible variations of electromagnetic energy density due to, for example, skin depth effect [<xref ref-type="bibr" rid="scirp.78934-ref23">23</xref>] inside some conductive material. Beyond this application, the interest is to give a physical meaning of taking into account those additional solutions in various energy spaces. The second part is dedicated to the Woodward effect and the possible relationship with the variation of EM energy density in some specific settings.</p><sec id="s5_2_1"><title>5.2.1. Variation of EM Energy Density</title><p>Thus, let us formulate the variation in time of energy density (u) at the second order with a Taylor series development such as:</p><p>d u = ∂ t u d t + ∂ t 2 u d t 2 2 + o ( d t 2 ) (30)</p><p>o is the Landau notation to omit higher order quantities. Note that at the</p><p>first order d u d t = ∂ t u . The higher orders term are based on the assumptions that</p><p>the EM waves inside the skin layer of the copper plate are evanescent waves and thus functions in the Schwartz space ( S − ( ℝ 4 ) -with 3 dimension variables and considering also the time ) [<xref ref-type="bibr" rid="scirp.78934-ref23">23</xref>] . As discussed before, those solutions are finite energy functions and in L ( ℝ 4 ) (i.e. following [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] and [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] , u = E ( f ( x 0 , y 0 , z 0 , T ) ) &lt; ∞ at some given point in the skin layer defined by the coordinates x 0 , y 0 , z 0 ). Now, using the Lemma 1 and the space M 0 k in Section 3, we can state in S − ( ℝ 4 )</p><p>M 0 k = { g ∈ S − ( ℝ 4 ) | g = ∂ t k f n ( x 0 , y 0 , z 0 , t )                     = α n ∂ t k − 1 f n − 2 ( x 0 , y 0 , z 0 , t ) ( Ψ 1 + ( f ( x 0 , y 0 , z 0 , t ) ) ) ,                     f ∈ S − ( ℝ 4 ) , n ∈ ℤ + − { 0 } , α n ∈ ℝ , z 0 ∈ [ 0, L ] , ( x 0 , y 0 ) ∈ [ 0, a ] 2 } (31)</p><p>Here f is either the electric or magnetic field (i.e. the absolute norm of E → and B → respectively). With the concept of multiplicity of solutions (e.g., Theorem 3). If g is a general solution of some linear PDEs, then f n can be identified as a special form of the solution (conditionally to its existence ).</p><p>Now considering the wave equation, the electric field and magnetic field are solutions and belong to the subspace M 0 k and associated with the variation of energy density ∂ t u . Furthermore, we can consider the solutions in M 0 1 associated with the variation of energy density ∂ t 2 u , which can be explained with the concept of multiplicity of solutions. The solutions of interest in M 0 1 are for the electric field g = ∂ t E and the magnetic field g = ∂ t B . The Taylor Series development of the energy of (for example) the electric field on a nominated position in space (i.e., x 0 , y 0 , z 0 ) and in an increment of time d t :</p><p>E ( E ( x 0 , y 0 , z 0 , T + d t ) ) = E ( E ( x 0 , y 0 , z 0 , T ) ) + ∑ k = 0 ∞     ∂ t k ( E 2 ( x 0 , y 0 , z 0 , T ) ) ( d t ) k k ! &lt; ∞</p><p>d E ( E ( x 0 , y 0 , z 0 , T + d t ) ) = ∑ k = 0 ∞     ∂ t k ( E 2 ( x 0 , y 0 , z 0 , T ) ) ( d t ) k k ! (32)</p><p>Finally one can write the relationship with the energy density following (26) and the previous Taylor series development for the electric and magnetic field:</p><p>0.5 ( ϵ 0 d E ( E ( x 0 , y 0 , z 0 , T + d t ) ) d t + 1 μ 0 d E ( B ( x 0 , y 0 , z 0 , T + d t ) ) d t ) = 0.5 ( ϵ 0 E 2 ( x 0 , y 0 , z 0 , T ) + 1 μ 0 B 2 ( x 0 , y 0 , z 0 , T ) ) + ∂ t u d t 2 + ∂ t 2 u d t 2 6 + o ( d t 2 ) (33)</p><p>Therefore, taking into account the second order term of the energy density ∂ t 2 u means that additional solutions should also be considered in the EM modeling. Note that in Appendix IV, we are taking an example of evanescent waves inside a copper wall (i.e. skin depth effect [<xref ref-type="bibr" rid="scirp.78934-ref23">23</xref>] ) and try to give further meaning to the consideration of higher order derivatives of the EM energy density where the additional solutions are defined with the energy spaces (e.g., ∂ t E and ∂ t B in M 0 1 ).</p></sec><sec id="s5_2_2"><title>5.2.2. Derivation of the Woodward Effect Using the Electromagnetic Energy Density</title><p>This section focuses on the derivation of the Woodward effect created in a asymmetric EM cavity (i.e. frustum) due to EM waves reflected on the cavity’s wall. Thus, the assumption is that the EM energy density variation results from the evanescent waves taking place in the skin depth of the asymmetric EM cavity’s walls.</p><p>1) Assumptions with the energy momentum relationship</p><p>When the Woodward effect was established in [<xref ref-type="bibr" rid="scirp.78934-ref4">4</xref>] [<xref ref-type="bibr" rid="scirp.78934-ref13">13</xref>] , the authors implicitly assumed the rest mass of the piezoelectric material via the famous Einstein’s relation in special relativity E = m c 2 ( E the rest energy associated with the rest mass m ) and its variation via electrostrictive effect.</p><p>Here, the system is the asymmetric EM cavity. The rest mass is all the particles within it at time t 0 when no charges are on the cavity’s walls. It excludes the photons considered with a null mass. Thus, the main assumption is that the EM excitation on the walls creates electric charges (i.e. electrons) which makes the rest mass varying with time. This assumption is the same as the mass variation of a capacitor between the charge and discharge times [<xref ref-type="bibr" rid="scirp.78934-ref24">24</xref>] . It allows us to state the variation of rest energy such as:</p><p>d E = E ( t + d t ) − E ( t ) = ( m ( t + d t ) − m ( t ) ) c 2 = d m c 2 (34)</p><p>Finally, the variation of rest energy d E is assumed to be equal to the variation of EM energy density ( d u ) resulting from the charges within the skin depth of the walls. We neglect any electrostrictive effects compared to the variation of EM energy density.</p><p>Note that at the particle level, the rest mass should satisfy the energy momen- tum relationship ( u e ) for a free body in special relativity [<xref ref-type="bibr" rid="scirp.78934-ref25">25</xref>] :</p><p>u e 2 = ( p c ) 2 + ( m e c 2 ) 2 p = v u e c 2 (35)</p><p>with p the momentum and m e the rest mass of the particle associated with the total energy u e . The particle is accelerated via the Lorentz force applied to the whole cavity with obviously v ≪ c . Thus, we have also the relationship p 2 &lt; ( u e / c ) 2 . In the remainder, we also use the elementary variation δ which becomes d for an infinitesimally small variation.</p><p>2) Derivation of the Woodward effect and relationship with EM energy density</p><p>If we define the mass density such as ρ = m / V , then from [<xref ref-type="bibr" rid="scirp.78934-ref13">13</xref>] , one can write the elementary mass variation per unit of volume</p><p>δ ρ = δ m V ~ d ρ   ( infinitesimally   small   variation )</p><p>d ρ = 1 4 π G [ 1 ρ ∂ t 2 ρ − 1 ρ 2 ( ∂ t ρ ) 2 ] (36)</p><p>Let us define the the rest energy E = ρ c 2 , then</p><p>d ρ = 1 4 π G [ 1 ρ c 2 ∂ t 2 E − 1 ( ρ c 2 ) 2 ( ∂ t E ) 2 ]</p><p>d ρ = 1 4 π G [ 1 E ∂ t 2 E − 1 ( E ) 2 ( ∂ t E ) 2 ] (37)</p><p>In some particular cases such as an EM cavity, we assume that the variation in time of the rest energy is equal to the variation of EM energy density u (i.e. ∂ t E ≃ ∂ t u ), but the rest energy is much bigger than the EM energy density E ≫ u . It allows then to state the relationship between the Woodward effect and the EM energy density</p><p>d ρ = 1 4 π G [ 1 E ∂ t 2 u − 1 ( E ) 2 ( ∂ t u ) 2 ] (38)</p><p>The EM energy density u follows the general definition of the sum of energy density from the electric ( u E ) and magnetic ( u B ) fields [<xref ref-type="bibr" rid="scirp.78934-ref23">23</xref>] . Finally, (38) can be seen as the definition of the EMG coupling.</p></sec></sec></sec><sec id="s6"><title>6. Conclusions</title><p>This work generalizes in the Schwartz space S − ( ℝ m ) , the framework on conjugate Teager-Kaiser energy operators established in [<xref ref-type="bibr" rid="scirp.78934-ref2">2</xref>] and [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] for the case m in [ 1,2 ] . The concept of multiplicity of solutions defined in [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] is also redefined here in Theorem 3. However, this concept uses the notion of energy spaces ( M p v ( p ∈ ℤ + , v ∈ ℤ + ), subspaces of S − ( ℝ ) defined previously in [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] and [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] . In order to generalize their definition as subspaces of S − ( ℝ m ) , the theory has been extended to some properties on the Hilbert spaces ( H 1 v ( � m ) ) on L 2 ( � m ) . In particular, we show in Properties 2 that M p 2 v ⊊ M p 1 v ( v 1 &lt; v 2 ) and the inclusion M p 1 v ⊊ H 1 v ( ℝ m ) ⊊ L 2 ( ℝ m ) .</p><p>The concept of multiplicity of solutions focuses on, generally speaking, looking for solutions of a given linear PDE specifically in the energy spaces. In this way, it is not following the classical way of solving a linear PDE with boundary conditions. Three examples illustrate this concept. The first one investigates some type of solutions (e.g., evanescent waves) of the wave equation when analysing the Taylor series development of the energy function associated with an evanescent wave. We then formulate another concept: the energy parallax. It is defined mathematically in Definition 4. Under some specific circumstances (e.g., the energy function exists), we show that the variations of energy locally in a predefined system, should lead to include additional solutions in the energy spaces with higher order v (in ℤ + ). The second example is based on the local variations of EM energy density, which allows to define waves which are first order derivative of the EM field. This example is further explored in Appendix IV. Finally, the last example is the derivation of the Woodward effect with some strong hypothesis in order to include the EM energy density in the specific case of asymmetric EM cavity. We introduce in the Woodward effect, the first and second order derivative of the EM energy density, which can be interpreted such as a theoretical definition of an Electromagnetic and Gravita- tional coupling.</p></sec><sec id="s7"><title>Acknowledgements</title><p>The author would like to acknowledge the important discussions with Dr. Jos&#233; Rodal and Prof. Heidi Fearn (California State University Fullerton, physics department) on the Woodward effect and its derivation from general relativity. In addition, we thank Prof. Paul Jolissaint (Universit&#233; de Neuch&#226;tel, Institute of Mathematics) for his kind advices on the Sobolev spaces, and Prof. Marc Troyanov together with Dr. Luigi Provenzano (Ecole Polytechnique F&#233;d&#233;rale de Lausanne, Institute of Mathematics) for the corrections on the Sobolev spaces and the Schwartz space. Finally, we acknowledge the constructive comments from the anonymous reviewers improving this work.</p></sec><sec id="s8"><title>Cite this paper</title><p>Montillet, J.-P. (2017) Sobolev Spaces, Schwartz Spaces, and a Definition of the Electromagnetic and Gravitational Coupling. Journal of Modern Physics, 8, 1700-1722. https://doi.org/10.4236/jmp.2017.810100</p></sec><sec id="s9"><title>Appendix I: Generalities on Sobolev Spaces</title><p>A Sobolev space is a vector space of functions equipped with a norm that is a combination of L<sup>p</sup>-norms of the function itself and its derivatives up to a given order. Intuitively, a Sobolev space is a space of functions with sufficiently many derivatives for some application domain, such as partial differential equations, and equipped with a norm that measures both the size and regularity of a function. Sobolev spaces are named after the Russian mathematician Sergei Sobolev.</p><p>Definition I.1. [<xref ref-type="bibr" rid="scirp.78934-ref14">14</xref>] Let Ω ⊆ ℝ m ( m ∈ ℕ + ) be open. The Sobolev space W k , p ( Ω ) ( k ∈ ℕ , p ∈ [ 1, ∞ ] ) is defined as:</p><p>W k , p ( Ω ) = { f ∈ L p ( Ω ) | D α f ∈ L p ( Ω ) , ∀ α ≤ k } (39)</p><p>with D α f the a-th partial derivative in multi index notation,</p><p>D α f = ∂ | α | f ∂ t 1 α 1 ⋯ ∂ t n α n . The Sobolev space W k , p ( Ω ) is the space of all locally</p><p>integrable functions f in Ω such as their partial derivatives D α f exist in the weak sense for all multi index α ≤ k and belongs to L p ( Ω ) (i.e. f L p &lt; ∞ ) ( [<xref ref-type="bibr" rid="scirp.78934-ref26">26</xref>] , chap. 5). If f lies in W k , p ( Ω ) , we define the W k , p norm of f by the formula</p><p>f W k , p ( Ω ) = ∑ | α | ≤ k D α f L p ( Ω ) (40)</p><p>Now, let us introduce the Fourier transform F : L 1 ( ℝ m ) → C b ( ℝ m ) as in [<xref ref-type="bibr" rid="scirp.78934-ref19">19</xref>]</p><p>F ( f ) = ∫ Ω f ( x ) e − i x ⋅ ξ d x = F ( f ) ( ξ ) (41)</p><p>Here C b ( ℝ m ) is the space of bounded and continuous functions in ℝ m [<xref ref-type="bibr" rid="scirp.78934-ref17">17</xref>] . Note that ⋅ is the scalar product (with x and ξ in ℝ m ). One can then define the Sobolev spaces for Ω = ℝ m , W k , p ( ℝ m ) using the Bessel potentials and the Fourier transform such as [<xref ref-type="bibr" rid="scirp.78934-ref14">14</xref>] or [<xref ref-type="bibr" rid="scirp.78934-ref17">17</xref>] (chap. 9) :</p><p>W k , p ( ℝ m ) = H k , p ( ℝ m ) : = { f ∈ L p ( ℝ m ) | F − 1 [ ( 1 + ξ 2 ) k / 2 F ( f ) ] ∈ L p ( ℝ m ) } (42)</p><p>The Bessel potential spaces are defined when replacing k by any real number s . They are Banach spaces and, for the special case p = 2 , Hilbert spaces . Now, one can state an important result with Sobolev spaces [<xref ref-type="bibr" rid="scirp.78934-ref14">14</xref>]</p><p>Theorem I.1.: W k , p ( ℝ m ) ⊆ W l , q ( ℝ m ) , whenever k &gt; l ≥ 0 and 1 ≤ p &lt; q &lt; ∞ are such that ( k − l ) p &lt; m</p><p>Proof. The proof of this theorem is rather long and technically delicate which is not our focus. Readers interested in this matter should refer to [<xref ref-type="bibr" rid="scirp.78934-ref14">14</xref>] [<xref ref-type="bibr" rid="scirp.78934-ref26">26</xref>] (chap. 5)</p></sec><sec id="s10"><title>Appendix II: Possible Interpretation of the Energy Parallax in Modern Physics</title><p>In Section 4, we define mathematically the notion of multiplicity of solutions for a given PDE. Through the various examples in Section 5, we define the concept of energy parallax. The general meaning is that additional solutions should be taken into account when varying the amount of energy. Those solutions should be defined based on the associated energy spaces (e.g., E p , p ∈ ℤ + ). Now, if we replace this concept in modern physics, what is the meaning behind it?</p><p>In modern physics, Energy is a global concept across the whole science. The definition varies with for example kinetic energy and potential energy in classical mechanics. It relates respectively to the object’s movement through space and function of its position within a field [<xref ref-type="bibr" rid="scirp.78934-ref27">27</xref>] . Chemical energy can be defined broadly such as the electrical potential energy among atoms and molecules. In quantum mechanics, energy is defined in terms of energy operators (e.g., Hamil- tonian) as a time derivative of the work function. It allows to define particles at nominated energy levels associated with an EM waves emitted at frequencies defined by the Planck’s relation. In General Relativity, energy results from the product of a varying mass and the square of the speed of light. Energy can describe the behavior of a system of two particles (and more). For example, the electron-positron annihilation in which rest mass (invariant mass) is destroyed. At the opposite, the inverse process (creator) in which the rest mass of the particle is created from energy of two (or more) annihilating photons [<xref ref-type="bibr" rid="scirp.78934-ref28">28</xref>] .</p><p>Energy parallax is here defined such as the concept of using additional wave functions. For example in Section 5.2.2 increasing the higher order derivatives of the EM energy density leads to the consideration of additional waves. The energy parallax concept can then help us to state that those additional waves are additional excited photons that we must take into account to vary the EM energy density.</p></sec><sec id="s11"><title>Appendix III: Discussion on the Possible Relationship between the Energy Spaces M p v and M p + 1 v − 1</title><p>This section follows the development in Section 3.2 and especially Properties 2. First, ∀ p ∈ ℤ + , M p + 1 v ∩ M p v ≠ { ∅ } , because 0 ∈ S − ( ℝ m ) , and ( ∀ p ∈ ℤ + , v ∈ ℤ + ) 0 ∈ M p v . Thus, 0 ∈ M p + 1 v ∩ M p v .</p><p>To recall Definition 2 and Lemma 2, ∂ i v ( [ [ f ] p ] 1 + ) n can be decomposed with the family of energy operators ( [ [   .   ] p + 1 ] k + ) k ∈ ℤ ( ∀ n ∈ ℤ + − { 0,1 } , i ∈ [ 1, ⋯ , m ] , p ∈ ℤ + , v ∈ ℤ + − { 0 } ). Thus, one can write ( l &lt; v ):</p><p>∂ i v ( [ [ f ] p ] 1 + ) n = ∑ j = 0 v − 1 ( v − 1 j ) ∂ i v − 1 − j ( [ [ f ] p ] 1 + ) n − l ∑ u = − N j N j C u [ [ ∂ i α u f ] p + 1 ] 1 + (43)</p><p>Thus, for n &gt; 1 , Lemma 2 allows to state that M p + 1 v ∩ M p v = { M p v , n &gt; 1 } , with M p v , n &gt; 1 the subspace of E p , but restricted for n ∈ ℤ + and n &gt; 1 .</p><p>Furthermore, let us define the space s p − , * ( ℝ m ) :</p><p>s p − , * ( ℝ m ) = { f ∈ S p − ( ℝ m ) | f ∉ ∪ i ∈ [ 1, ⋯ , m ] ( ∪ k ∈ ℤ K e r ( [ [ f ] p ] k , i + ) ) } (44)</p><p>Note that M p ∞ ⊂ s p − , * ( ℝ m ) , but the bump functions [<xref ref-type="bibr" rid="scirp.78934-ref15">15</xref>] are included in s p − , * ( ℝ m ) . We can also recall the discussion on n = 1 in [<xref ref-type="bibr" rid="scirp.78934-ref1">1</xref>] and [<xref ref-type="bibr" rid="scirp.78934-ref3">3</xref>] , with the definition</p><p>{ ∀ f ∈ s p − , * ( ℝ m ) , p ∈ ℤ + | ∃ g ∈ S − ( ℝ m ) , g = 1 ( [ [ f ] p ] 1 + ) n , ∀ n ∈ ℤ + , n &gt; 1 } (45)</p><p>On can also state that ∂ i k ( [ [ f ] p ] 1 + ) = ∂ i k ( ( [ [ f ] p ] 1 + ) 3 ( [ [ f ] p ] 1 + ) 2 ) ( k ∈ ℤ + ) and use the</p><p>Leibniz’s rule for derivations in order to expand the multiple derivatives or the decomposition stated in Lemma 2. If we call M p v , * ( p ∈ ℤ + ), the subspaces of s p − , * ( ℝ m ) . For all g 1 in M p v , * can be written as a non linear sum of g 2 in M p + 1 v , * . Finally, we can conclude that M p v , * ⊊ M p + 1 v , * . With the specific extension of M p v to the case n = 1 , we can also conclude M p v ⊊ M p + 1 v . In addition,</p><p>s p − , * ( ℝ m ) ⊊ s p + 1 − , * ( ℝ m ) by definition.</p></sec><sec id="s12"><title>Appendix IV: Consequences in Terms of EM theory</title><p>We are taking the example of the variation of EM energy density inside a copper wall due to planar waves reflecting and refracting on it [<xref ref-type="bibr" rid="scirp.78934-ref23">23</xref>] . To recall Section 5.2, the EM field is now including ( E , δ E ) and ( B , δ B ), contribution of the subspaces M 0 0 and M 0 1 respectively when using the concept of multiplicity of the solutions (i.e. Theprem 3) for the higher order derivatives of the energy density (see (26)). We call the total EM field E t o t and B t o t inside the copper plate (skin layer) with associated permittivity ϵ and permeability μ . They are solutions of the Maxwell equations:</p><p>∇ ⋅ E t o t = ρ t o t ϵ , ∇ &#215; E t o t = − ∂ t B t o t , ∇ ⋅ B t o t = 0 , ∇ &#215; B t o t = μ ϵ ∂ t E t o t + μ j , } (46)</p><p>with the principle of charge conservation:</p><p>∂ t ρ t o t + ∇ ⋅ j = 0 (47)</p><p>Now, the variation of energy density (26) together with the equation of charge conservation is formulated such as:</p><p>d u d t + ∇ ⋅ P t o t = ∇ ⋅ E t o t (48)</p><p>P t o t = E t o t &#215; B t o t μ is the Poynting vector. Now, writing E t o t = E + δ E , B t o t = B + δ B and δ is the first derivative in time ( ∂ t ) (i.e. solutions in M 0 1 ),</p><p>then following [<xref ref-type="bibr" rid="scirp.78934-ref23">23</xref>]</p><p>( E + ∂ t E ) ⋅ j = ( E + ∂ t E ) ⋅ [ 1 μ ∇ &#215; ( B + ∂ t B ) − ϵ ∂ t ( E + ∂ t E ) ] (49)</p><p>using the equalities ∇ ⋅ ( E &#215; B ) = B ⋅ ∇ &#215; E − E ⋅ ∇ &#215; B and the Maxwell equation ∇ &#215; E = − ∂ t B , ∇ &#215; ∂ t E = − ∂ t 2 B the previous equation reduces to:</p><p>E ⋅ j + ∇ ⋅ ( E &#215; B μ ) + ∂ t u + ∂ t E ⋅ j + ∇ ⋅ ( ∂ t E &#215; ∂ t B μ ) + ∂ t 2 u + ∇ ⋅ ( ∂ t E &#215; B μ ) + ∇ ⋅ ( E &#215; ∂ t B μ ) + ∂ B ⋅ ∂ B μ + ε ∂ t E ⋅ ∂ t E = 0 (50)</p><p>We can separate in three groups,</p><p>∂ t u + ∇ ⋅ ( E &#215; B μ ) = − j ⋅ E ∂ t 2 u + ∇ ⋅ ( ∂ t E &#215; B μ ) + ∇ ⋅ ( E &#215; ∂ t B μ ) = − j ⋅ ∂ t E ∇ ⋅ ( ∂ t E &#215; ∂ t B μ ) = − ∂ t B ⋅ ∂ t B μ − ε ∂ t E ⋅ ∂ t E }</p><p>The Poynting vector is defined as P = E &#215; B μ and its derivative ∂ t P = ∂ t E &#215; B μ + E &#215; ∂ t B μ . Thus, the second order term of the energy density is</p><p>the contribution of the EM field generated by ∂ t E and ∂ t B is:</p><p>∂ t u + ∇ ⋅ P = − j ⋅ E ∂ t 2 u + ∇ ⋅ ( ∂ t P ) = − j ⋅ ∂ t E ∇ ⋅ ( ∂ t E &#215; ∂ t B μ ) = − ∂ t B ⋅ ∂ t B μ − ε 0 ∂ t E ⋅ ∂ t E }</p><p>The last line is the contribution from only the fields ∂ t E and ∂ t B .</p><p>Finally, the creation of the wave defined by the EM field ( ∂ t E , ∂ t B ) means that some material properties may allow to create two type of EM waves namely ( E , B ) and ( ∂ t E , ∂ t B ).</p></sec></body><back><ref-list><title>References</title><ref id="scirp.78934-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Montillet, J.P. (2017) Bulletin of Mathematical Analysis and Applications, 9, 134-150. http://bmathaa.org/repository/docs/BMAA9-1-13.pdf</mixed-citation></ref><ref id="scirp.78934-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Montillet, J.P. (2013) Acta Applicandae Mathematicae, 129, 61-80.</mixed-citation></ref><ref id="scirp.78934-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Montillet, J.P. (2015) Acta Applicandae Mathematicae, 140.https://doi.org/10.1007/s10440-014-9978-9</mixed-citation></ref><ref id="scirp.78934-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Woodward, J.F. (2001) Foundations of Physics, 31, 819. https://doi.org/10.1023/A:1017500513005</mixed-citation></ref><ref id="scirp.78934-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Kaiser, J.F. (1990) On a Simple Algorithm to Calculate the “Energy” of a Signal. Proceedings of IEEE Conference on Acoustics, Speech, and Signal Processing, Vol. 1. https://doi.org/10.1109/ICASSP.1990.115702</mixed-citation></ref><ref id="scirp.78934-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Maragos, J.F. and Potamianos, A. (1995) IEEE Signal Processing Letters, 2, 152-154. https://doi.org/10.1109/97.404130</mixed-citation></ref><ref id="scirp.78934-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Dunn, R.B., Quatieri, T.F. and Kaiser, J.F. (1993) Detection of Transient Signals using the Energy Operator. Proceedings of IEEE Conference on Acoustics, Speech, and Signal Processing, Minneapolis, 27-30 April 1993, Vol. 3, 145-148.</mixed-citation></ref><ref id="scirp.78934-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Bovik, A.C., Havlicek, J.P. and Desai, M.D. (1993) Theorems for Discrete Filtered Modulated Signals. Proceedings of IEEE Conference on Acoustics, Speech, and Signal Processing, Minneapolis, 27-30 April 1993, Vol. 3, 153-156.</mixed-citation></ref><ref id="scirp.78934-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Cexus, J.C., Boudraa, A., Ardeyeh, F.H. and Diop, E.H.S. (2010) 2D Cross-ΨB-Energy Operator for Images Analysis. Proceedings of International Symposium on Communications, Control and Signal Processing, Limassol, 3-5 March 2010, 1-4.</mixed-citation></ref><ref id="scirp.78934-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Montillet, J.P. (2010) International Mathematical Forum, 5, 2387-2400. http://www.m-hikari.com/imf-2010/45-48-2010/montilletIMF45-48-2010.pdf</mixed-citation></ref><ref id="scirp.78934-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Hunter, J.K. (2009) Lecture Notes, Chapter 5. https://www.math.ucdavis.edu/~hunter/m218a_09/ch5A.pdf</mixed-citation></ref><ref id="scirp.78934-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Fearn, H. (2015) Journal of Modern Physics, 6, 260-272.</mixed-citation></ref><ref id="scirp.78934-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Woodward, J.F. (2004) AIP Conference Proceedings, 699, 1127. https://doi.org/10.1063/1.1649682</mixed-citation></ref><ref id="scirp.78934-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Adams, R. and Fournier, J. (2003) Sobolev Spaces. Elsevier Academic Press.</mixed-citation></ref><ref id="scirp.78934-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Johnson, S.G. (2015) Saddle-Point Integration of C “Bump” Functions. https://arxiv.org/abs/1508.04376</mixed-citation></ref><ref id="scirp.78934-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Keyl, M., Kiukas, J. and Werner, R.F. (2016) Reviews in Mathematical Physics, 28, Article ID: 1630001. https://doi.org/10.1142/S0129055X16300016</mixed-citation></ref><ref id="scirp.78934-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Folland, G.B. (1999) Real Analysis: Modern Techniques and Their Applications. Wiley-Interscience Series of Texts, Monographs and Tracts.</mixed-citation></ref><ref id="scirp.78934-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Hunter, J.K. (2012) Introduction to Real Analysis. Lecture Notes, Chapter 7. https://www.math.ucdavis.edu/~hunter/m125a/intro_analysis_ch7.pdf</mixed-citation></ref><ref id="scirp.78934-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Reiter, M. and Schuster, A. (2008) Fourier Transform and Sobolev Spaces. Lecture Notes, Summer Term. http://www.mat.univie.ac.at/~stein/lehre/SoSem08/sobolev_fourier.pdf</mixed-citation></ref><ref id="scirp.78934-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">Van den Ban, E. and Crainic, M. (2013) Analysis on Manifolds, Lecture 4: Fourier Transforms. https://www.staff.science.uu.nl/~ban00101/anman2013/anman2013_ch4.pdf</mixed-citation></ref><ref id="scirp.78934-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">Saint Raymond, X. (1991) Elementary Introduction to the Theory of Pseudo Differential Operators. CRC Press, Boca Raton.</mixed-citation></ref><ref id="scirp.78934-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">Hardy, G.H., Littlewood, J.E. and Plya, G. (1952) Inequalities. Cambridge University Press, Cambridge.</mixed-citation></ref><ref id="scirp.78934-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">Petit, R. (1993) Ondes Electromagnetiques en Radioélectricité et en Optique, Masson.</mixed-citation></ref><ref id="scirp.78934-ref24"><label>24</label><mixed-citation publication-type="other" xlink:type="simple">Porcelli, E.B. and Filho, V.S. (2016) Physics Essays, 29.https://doi.org/10.4006/0836-1398-29.1.002</mixed-citation></ref><ref id="scirp.78934-ref25"><label>25</label><mixed-citation publication-type="other" xlink:type="simple">Moller, C. (2011) The Theory of Relativity. Nabu Press.</mixed-citation></ref><ref id="scirp.78934-ref26"><label>26</label><mixed-citation publication-type="other" xlink:type="simple">Evans, L.C. (2010) Partial Differential Equations. American Mathematical Society. https://doi.org/10.1090/gsm/019</mixed-citation></ref><ref id="scirp.78934-ref27"><label>27</label><mixed-citation publication-type="other" xlink:type="simple">Feynman, R. (1964) The Feynman Lectures on Physics. http://www.feynmanlectures.caltech.edu/</mixed-citation></ref><ref id="scirp.78934-ref28"><label>28</label><mixed-citation publication-type="other" xlink:type="simple">Vinogradov, I.M. (1977) Encyclopedia of Mathematics. Kluwer Academic Publishers.</mixed-citation></ref></ref-list></back></article>