DOCSLIB.ORG

(2018) Observation of Asymmetric Electromagnetic Field Profiles In

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
(2018) Observation of Asymmetric Electromagnetic Field Profiles In PHYSICAL REVIEW B 97, 085105 (2018) Observation of asymmetric electromagnetic field profiles in chiral metamaterials Nobuyuki Hisamoto,1 Tetsuya Ueda,1,* Kei Sawada,2 and Satoshi Tomita3 1Department of Electrical Engineering and Electronics, Kyoto Institute of Technology, Matsugasaki, Sakyo-ku, Kyoto 606-8585, Japan 2RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo, Hyogo 679-5148, Japan 3Graduate School of Materials Science, Nara Institute of Science and Technology, Ikoma, Nara 630-0192, Japan (Received 29 June 2017; published 2 February 2018) We experimentally observe asymmetric electromagnetic field profiles along two-dimensional chiral metamate- rials. The asymmetric field profiles depending on the chirality and the operation frequency have been reproduced well by the numerical simulation. Around a chiral meta-atom, distribution of a Poynting vector is found to be shifted asymmetrically. These results are explained in terms of an analogy with the side-jump mechanism in the electronic anomalous Hall systems. DOI: 10.1103/PhysRevB.97.085105 I. INTRODUCTION chiral metamaterials, determine the macroscopic constitutive relations of the effective media including not only permittivity The electromagnetic phenomena such as unidirectional and permeability but also chirality for bi-anisotropic media, by electromagnetic wave propagation have been of much in- appropriately designing meta-atoms [26–34]. So far, various terest. Unidirectional backscattering immune propagation of applications of the chiral metamaterials have been proposed surface waves, referred to as electromagnetic edge modes, and demonstrated [35–39]. However, little is known about the was observed in electromagnetic band gap structures such as edgelike mode propagation in the chiral metamaterials at the photonic crystals [1–3] with topological properties [4,5]. The operation frequencies far below the photonic band gaps due to unidirectional surface wave propagations that appear at the the Bragg scattering. boundary between two different areas with and without broken In this paper we demonstrate electromagnetic edgelike time-reversal symmetry using ferromagnetic medium [6,7] mode propagation in chiral metamaterials at microwave fre- or broken space-inversion symmetry due to chirality [8–13] quencies. We consider two-dimensional (2D) microwave chiral are topologically protected and robust against scattering by metamaterials of triangular lattice that are constructed us- disorder. ing helical meta-atoms [40,41]. Experimental results show Even in nontopological systems without periodicity, asym- asymmetric magnetic field profiles, which are relevant to metric propagation of surface plasmon polaritons due to the edgelike mode propagation in the metamaterials. The field unidirectional excitation [14–17] has been demonstrated. In concentration side is switched by inverting the chirality of these systems, unidirectional wave propagation was observed the meta-atom and by changing the operation frequencies. in various structures with space-inversion symmetry by ap- The asymmetric field profiles are reproduced by numerical propriately setting up incident wave conditions, such as by simulations based on the finite element method. Furthermore, controlling source positions [14], incident angles [15,16], and the numerical simulations visualize asymmetric energy and the polarization [17]. However, these approaches focused on Poynting vector distribution in the metamaterials, which en- wave incidence techniques rather than a design of wave- ables us to understand the behaviors of asymmetric field guiding structures. profiles. At frequencies far below the electromagnetic band gaps, where the wavelength of propagating waves becomes much longer than the lattice constant, the artificial periodic structures II. EXPERIMENTAL AND NUMERICAL SETUPS can be regarded as effective media based on the concept of metamaterials [18]. Metamaterials are artificial electromag- Figure 1 shows schematic illustrations of 2D chiral metama- netic media which are composed of small elements compared terials and the calculated dispersion diagram. A right-handed to the wavelength in order to manipulate electromagnetic wave (RH) chiral meta-atom was composed of a metallic-wire RH propagation and to discover new electromagnetic phenomena. helix. The chiral axis is along the z direction. As illustrated in On the other hand, chirality with broken space-inversion Fig. 1(a), design parameters for the meta-atom are as follows: = = symmetry, such as helices [19,20] or gammadions [21], gives wire diameter dw 0.4 mm, numbers of turns in the helix t = = rise to optical activity [22–25], which can rotate a plane of 3.5, pitch of turns p 1 mm, and helix diameter dh 2.4 mm. the polarization or may convert linear to circular polarizations The resonance frequency of the lowest mode for the and vice versa. Metamaterials with the chirality, referred to as isolated chiral meta-atom is calculated to be 7.0 GHz, which is determined by the total length of metal wire corresponding to a half wavelength of microwaves. Array of the meta-atoms comprises the 2D chiral metamaterial with thickness of *[email protected] a single meta-atom. Figure 1(c) illustrates the dispersion 2469-9950/2018/97(8)/085105(6) 085105-1 ©2018 American Physical Society HISAMOTO, UEDA, SAWADA, AND TOMITA PHYSICAL REVIEW B 97, 085105 (2018) FIG. 1. (a) Side view of a RH chiral meta-atom. (b) Top view of a finite rhombus-shaped 2D RH chiral metamaterial with triangular FIG. 2. Measured [(a)–(c)] and simulated [(d)–(f)] profiles for the lattice. Red arrow at the bottom indicates the magnetic-dipole excita- y component of the magnetic fields. (a) and (d) at 6.85 GHz for the tion position. (c) Dispersion diagram of a 2D RH chiral metamaterial LH chiral metamaterial, (b) and (e) at 6.85 GHz for the RH chiral with infinite triangular lattice. (d) The enlarged dispersion diagram in metamaterial. (c) and (f) also for the RH chiral meta-atoms but at 7.5 the vicinity of the lowest resonant mode in (c). and 8.4 GHz, respectively. will be a future issue. Here we focus our attention on the lowest diagram calculated for the 2D RH chiral metamaterials band at approximately 7.0 GHz where such topological effects consisting of the RH chiral meta-atoms in infinite triangular are negligibly small. lattice under the periodic boundary condition. In the numerical We experimentally and numerically investigate finite simulation setup for the eigenmode analysis, a perfect magnetic rhombus-shaped 2D chiral metamaterials as shown in Fig. 1(b). conductor (PMC) has been set on top and bottom boundaries The meta-atoms with the chiral axis parallel to the z axis were of the calculation region, because we are focusing on the TE placed on each grid of the 2D triangular lattice on the x-y plane mode propagation near the lowest-order magnetic dipolelike with 7 × 7 cells. The meta-atom design is the same as shown resonance of chiral meta-atoms. The lattice spacing among in Fig. 1(a). In microwave experiments, RH and LH chiral the meta-atoms were set to 8 mm. meta-atoms were made of stainless-steel wires. Forty-nine For the eigenmode calculation, a commercially available RH(LH) chiral meta-atoms embedded in a polystyrene foam electromagnetic field solver based on the finite element method slab comprise the 2D RH(LH) chiral metamaterials so that the ANSYS HFSS ver. 16 was utilized. Although Fig. 1(c) shows wire top edges are aligned to the −x direction. the dispersion diagram of the 2D RH chiral metamaterial, the In the measurement of microwave field profiles, two diagram is independent of the sign of chirality and identical identical loop antennas with diameter of 2.4 mm were used; to that of 2D left-handed (LH) chiral metamaterials consisting one for a magnetic dipole excitation source and another for a of the LH chiral meta-atoms. Figure 1(d) shows the enlarged magnetic field detection probe. Both antennas are connected to dispersion diagram in the vicinity of the lowest resonant mode two ports of a vector network analyzer (Agilent PNA N5224A) in Fig. 1(c).InFig.1(d), mutual coupling among neighboring through coaxial cables. While the meta-atom at the bottom meta-atoms under the magnetic-dipolelike lowest-mode reso- vertex of the rhombus designated by a red arrow in Fig. 1(b) nance at approximately 7.0 GHz gives rise to a backward-wave was excited by the excitation source antenna, magnetic field mode propagation. In the present case, the passband width of distribution at the slab top surface was measured by scanning the mode is extremely narrow due to the weak coupling. The the magnetic field detection probe in the x-y plane at a distance coupling between the chiral backward-wave mode and TEM of 1 mm on the opposite side of the excitation. In order to mode in free space results in a narrow gap at approximately reproduce experimental results by the numerical simulation, a 7.2 GHz between the and K points as well as between and time-varying magnetic dipole moment normal to the lattice is M points. placed as an incident wave source just below the meta-atom Additionally, Fig. 1(c) shows a band gap due to the Bragg at the bottom vertex of a rhombus at a distance of 3 mm. scattering at the K and M points around 20 GHz. We have numerically confirmed that this gap is caused by the meta- atom’s chirality. When the meta-atom structure is achiral such III. RESULTS AND DISCUSSION as with finite-length straight wire or ring shape, this gap will disappear. It would be interesting to control the dispersion A. Measured and calculated magnetic field profiles by changing the meta-atom’s chirality among RH, LH, and Figure 2 shows the measured and calculated magnitude achiral, giving topological states to the electromagnetic waves. of the magnetic field distribution for the y component However, such a study is out of the scope of this paper, and around 7.0 GHz corresponding to the magnetic-dipolelike 085105-2 OBSERVATION OF ASYMMETRIC ELECTROMAGNETIC … PHYSICAL REVIEW B 97, 085105 (2018) lowest-mode resonance frequency of the meta-atoms.
Load more

Recommended publications
  • Chapter 3 Electromagnetic Waves & Maxwell's Equations
    Chapter 3 Electromagnetic Waves & Maxwell’s Equations Part I Maxwell’s Equations Maxwell (13 June 1831 – 5 November 1879) was a Scottish physicist. Famous equations published in 1861 Maxwell’s Equations: Integral Form Gauss's law Gauss's law for magnetism Faraday's law of induction (Maxwell–Faraday equation) Ampère's law (with Maxwell's addition) Maxwell’s Equations Relation of the speed of light and electric and magnetic vacuum constants 1 c = 2.99792458 108 [m/s] 0 0 permittivity of free space, also called the electric constant permeability of free space, also called the magnetic constant Gauss’s Law For any closed surface enclosing total charge Qin, the net electric flux through the surface is This result for the electric flux is known as Gauss’s Law. Magnetic Gauss’s Law The net magnetic flux through any closed surface is equal to zero: As of today there is no evidence of magnetic monopoles See: Phys.Rev.Lett.85:5292,2000 Ampère's Law The magnetic field in space around an electric current is proportional to the electric current which serves as its source: B ds 0I I is the total current inside the loop. ds B i1 Direction of integration i 3 i2 Faraday’s Law The change of magnetic flux in a loop will induce emf, i.e., electric field B E ds dA A t Lenz's Law Claim: Direction of induced current must be so as to oppose the change; otherwise conservation of energy would be violated. Problem with Ampère's Law Maxwell realized that Ampere’s law is not valid when the current is discontinuous.
    [Show full text]
  • Poynting Vector and Power Flow in Electromagnetic Fields
    Poynting Vector and Power Flow in Electromagnetic Fields: Electromagnetic waves can transport energy as a result of their travelling or propagating characteristics. Starting from Maxwell's Equations: Together with the vector identity One can write In simple medium where and are constant, and Divergence theorem states, This equation is referred to as Poynting theorem and it states that the net power flowing out of a given volume is equal to the time rate of decrease in the energy stored within the volume minus the conduction losses. In the equation, the following term represents the rate of change of the stored energy in the electric and magnetic fields On the other hand, the power dissipation within the volume appears in the following form Hence the total decrease in power within the volume under consideration: Here (W/mt2) is called the Poynting vector and it represents the power density vector associated with the electromagnetic field. The integration of the Poynting vector over any closed surface gives the net power flowing out of the surface. Poynting vector for the time harmonic case: Using the convention, the instantaneous value of a quantity is the real part of the product of a phasor quantity and when is used as reference. Considering the following phasor: −푗훽푧 퐸⃗ (푧) = 푥̂퐸푥(푧) = 푥̂퐸0푒 The instantaneous field becomes: 푗푤푡 퐸⃗ (푧, 푡) = 푅푒{퐸⃗ (푧)푒 } = 푥̂퐸0푐표푠(휔푡 − 훽푧) when E0 is real. Let us consider two instantaneous quantities A and B such that where A and B are the phasor quantities. Therefore, Since A and B are periodic with period , the time average value of the product form AB, 푇 1 퐴퐵 = ∫ 퐴퐵푑푡 푎푣푒푟푎푔푒 푇 0 푇 1 퐴퐵 = ∫|퐴||퐵|푐표푠(휔푡 + 훼)푐표푠(휔푡 + 훽)푑푡 푎푣푒푟푎푔푒 푇 0 1 퐴퐵 = |퐴||퐵|푐표푠(훼 − 훽) 푎푣푒푟푎푔푒 2 For phasors, and , where * denotes complex conjugate.
    [Show full text]
  • Electromagnetic Energy
    Physics 142 Electromagnetic Enmergy Page !1 Electromagnetic Energy A child of five can understand this; send someone to fetch a child of five. — Groucho Marx Energy in the fields can move from place to place We have discussed the energy in an electrostatic field, such as that stored in a capacitor. We found that this energy can be thought of as distributed in space, with an energy per unit volume (energy density) at each point in space. We obtained a formula for this 1 2 quantity, ! ue = 2 ε0E (if there are no dielectric materials). This formula applies to any E- field. We also found that magnetic fields possess energy distributed in space, described 2 by the magnetic energy density ! um = B /2µ0 . If there are both electric and magnetic fields, the total electromagnetic energy density is the sum of ! ue and ! um . These specify how much electromagnetic field energy there is at any point in space. But we have not yet considered how this energy moves from place to place. Consider the energy flow in a flashlight. We know that energy moves from the battery to the bulb, where it is converted into heat and light. It is tempting to assume that this energy flows through the conductors, like water in a pipe. But if we look carefully we find that the electromagnetic energy density in the conductors is much too small to account for the amount of energy in transit. Nearly all of the energy gets to the bulb by flowing through space near the conductors. In a sense they guide the energy but do not carry much of it.
    [Show full text]
  • Poynting's Theorem and the Wave Equation
    Chapter 18: Poynting’s Theorem and the Wave Equation Chapter Learning Objectives: After completing this chapter the student will be able to: Use Poynting’s theorem to determine the direction and magnitude of power flow in an electromagnetic system. Use Maxwell’s Equations to derive a general homogeneous wave equation for the electric and magnetic field. Derive a simplified wave equation assuming propagation in a vacuum and an electric field polarized in only one direction. Use Maxwell’s Equations to derive the speed of light in a vacuum. You can watch the video associated with this chapter at the following link: Historical Perspective: John Henry Poynting (1852-1914) was an English physicist who did work in electromagnetic energy flow, elasticity, and astronomy. He coined the term “Greenhouse Effect.” Both the Poynting Vector and Poynting’s Theorem are named in his honor. Photo credit: https://upload.wikimedia.org/wikipedia/commons/5/5f/John_Henry_Poynting.jpg, [Public domain], via Wikimedia Commons. 1 18.1 Poynting’s Theorem With Maxwell’s Equations, we now have the tools necessary to derive Poynting’s Theorem, which will allow us to perform many useful calculations involving the direction of power flow in electromagnetic fields. We will begin with Faraday’s Law, and we will take the dot product of H with both sides: (Copy of Equation 16.24) (Equation 18.1) Next, we will start with Ampere’s Law and will take the dot product of E with both sides: (Copy of Equation 17.13) (Equation 18.2) Now, let’s subtract both sides of Equation 18.2 from both sides of Equation 18.1: (Equation 18.3) We can now apply the following mathematical identity to the left side of Equation 18.3: (Equation 18.4) This substitution yields: (Equation 18.5) Distributing the E across the right side gives: (Equation 18.6) 2 Now let’s concentrate on the first time on the right side.
    [Show full text]
  • The Poynting Vector: Power and Energy in Electromagnetic fields
    The Poynting vector: power and energy in electromagnetic fields Kenneth H. Carpenter Department of Electrical and Computer Engineering Kansas State University October 19, 2004 1 Conservation of energy in electromagnetics The concept of conservation of energy (along with conservation of momen- tum) is one of the basic principles of physics – both classical and modern. When dealing with electromagnetic fields a way is needed to relate the con- cept of energy to the fields. This is done by means of the Poynting vector: P = E × H. (1) In eq.(1) E is the electric field intensity, H is the magnetic field intensity, and P is the Poynting vector, which is found to be the power density in the electromagnetic field. The conservation of energy is then established by means of the Poynting theorem. 2 The Poynting theorem By using the Maxwell equations for the curl of the fields along with Gauss’s divergence theorem and an identity from vector analysis, we may prove what is known as the Poynting theorem. 1 EECE557 Poynting vector – supplement to text - Fall 2004 2 The Maxwell’s equations needed are ∂B ∇ × E = − (2) ∂t ∂D ∇ × H = J + (3) ∂t along with the material relationships D = 0E + P (4) B = µ0H + µ0M (5) or for isotropic materials D = E (6) B = µH (7) In addition, the identity from vector analysis, ∇ · (E × H) ≡ −E · (∇ × H) + H · (∇ × E), (8) is needed. 2.1 The derivation If P is to be power density, then its surface integral over the surface of a volume must be the power out of the volume.
    [Show full text]
  • PHYS 110B - HW #4 Fall 2005, Solutions by David Pace Equations Referenced As ”EQ
    PHYS 110B - HW #4 Fall 2005, Solutions by David Pace Equations referenced as ”EQ. #” are from Griffiths Problem statements are paraphrased [1.] Problem 8.2 from Griffiths Reference problem 7.31 (figure 7.43). (a) Let the charge on the ends of the wire be zero at t = 0. Find the electric and magnetic fields in the gap, E~ (s, t) and B~ (s, t). (b) Find the energy density and Poynting vector in the gap. Verify equation 8.14. (c) Solve for the total energy in the gap; it will be time-dependent. Find the total power flowing into the gap by integrating the Poynting vector over the relevant surface. Verify that the input power is equivalent to the rate of increasing energy in the gap. (Griffiths Hint: This verification amounts to proving the validity of equation 8.9 in the case where W = 0.) Solution (a) The electric field between the plates of a parallel plate capacitor is known to be (see example 2.5 in Griffiths), σ E~ = zˆ (1) o where I define zˆ as the direction in which the current is flowing. We may assume that the charge is always spread uniformly over the surfaces of the wire. The resultant charge density on each “plate” is then time-dependent because the flowing current causes charge to pile up. At time zero there is no charge on the plates, so we know that the charge density increases linearly as time progresses. It σ(t) = (2) πa2 where a is the radius of the wire and It is the total charge on the plates at any instant (current is in units of Coulombs/second, so the total charge on the end plate is the current multiplied by the length of time over which the current has been flowing).
    [Show full text]
  • PHYS 332 Homework #1 Maxwell's Equations, Poynting Vector, Stress
    PHYS 332 Homework #1 Maxwell's Equations, Poynting Vector, Stress Energy Tensor, Electromagnetic Momentum, and Angular Momentum 1 : Wangsness 21-1. A parallel plate capacitor consists of two circular plates of area S (an effectively infinite area) with a vacuum between them. It is connected to a battery of constant emf . The plates are then slowly oscillated so that the separation d between ~ them is described by d = d0 + d1 sin !t. Find the magnetic field H between the plates produced by the displacement current. Similarly, find H~ if the capacitor is first disconnected from the battery and then the plates are oscillated in the same manner. 2 : Wangsness 21-7. Find the form of Maxwell's equations in terms of E~ and B~ for a linear isotropic but non-homogeneous medium. Do not assume that Ohm's Law is valid. 3 : Wangsness 21-9. Figure 21-3 shows a charging parallel plate capacitor. The plates are circular with radius a (effectively infinite). Find the Poynting vector on the bounding surface of region 1 of the figure. Find the total rate at which energy is entering region 1 and then show that it equals the rate at which the energy of the capacitor is increasing. 4a: Extra-2. You have a parallel plate capacitor (assumed infinite) with charge densities σ and −σ and plate separation d. The capacitor moves with velocity v in a direction perpendicular to the area of the plates (NOT perpendicular to the area vector for the plates). Determine S~. b: This time the capacitor moves with velocity v parallel to the area of the plates.
    [Show full text]
  • Lecture 16 Waves in Isotropic Media: Dielectrics And
    Lecture 16 Waves in Isotropic Media: Dielectrics and Conductors In this lecture you will learn: • Wave propagation in dielectric media • Waves propagation in conductive media ECE 303 – Fall 2005 – Farhan Rana – Cornell University Review: Plane Waves in Free Space Faraday’s Law: r r r r ∇ × E()r = − j ω µo H ()r Ampere’s Law: r r r v r r ∇ × H()r = J()r + j ω εo E(r ) r Complex Wave Equation: Assume: J(rv) = ρ(rr) = 0 r r r r 2 r r ∇ × ∇ × E()r = − j ω µo ∇ × H(r ) = ω µo εo E(r ) 0 r r 2 r r 2 r r ⇒ ∇()∇ .E()r − ∇ E()r = ω µo εo E ()r 2 r r 2 r r ⇒ ∇ E()r = −ω µo εo E ()r For a plane wave in free space we know the E-field and H-field phasors to be: r r r r ˆ − j k .r ω E r E()r = n Eo e k = ω µ ε = k o o c r E r r H()rr = kˆ × nˆ o e− j k .r µ ( ) η = o ≈ 377 Ω H ηo o εo ECE 303 – Fall 2005 – Farhan Rana – Cornell University 1 Waves in a Dielectric Medium – Wave Equation Suppose we have a plane wave of the form, E r r r E()rr = nˆ E e− j k .r o ε traveling in an infinite dielectric medium with permittivity ε H What is different from wave propagation in free space? Faraday’s Law: r r r r ∇ × E()r = − j ω µo H ()r Ampere’s Law: r r r ∇ × H()rr = J()rv + j ω ε E(rr) r Complex Wave Equation: Assume: J(rv) = ρ(rr) = 0 r r r r 2 r r compare with the ∇ × ∇ × E()r = − j ω µo ∇ × H(r ) = ω µo ε E(r ) 0 complex wave equation r r 2 r r 2 r r in free space ⇒ ∇()∇ .E()r − ∇ E()r = ω µo ε E ()r 2 r r 2 r r 2 r r 2 r r ⇒ ∇ E()r = −ω µo ε E ()r ∇ E(r )()= −ω µo εo E r ECE 303 – Fall 2005 – Farhan Rana – Cornell University Waves in a Dielectric Medium
    [Show full text]
  • Electromagnetism - Lecture 13
    Electromagnetism - Lecture 13 Waves in Insulators Refractive Index & Wave Impedance • Dispersion • Absorption • Models of Dispersion & Absorption • The Ionosphere • Example of Water • 1 Maxwell's Equations in Insulators Maxwell's equations are modified by r and µr - Either put r in front of 0 and µr in front of µ0 - Or remember D = r0E and B = µrµ0H Solutions are wave equations: @2E µ @2E 2E = µ µ = r r r r 0 r 0 @t2 c2 @t2 The effect of r and µr is to change the wave velocity: 1 c v = = pr0µrµ0 prµr 2 Refractive Index & Wave Impedance For non-magnetic materials with µr = 1: c c v = = n = pr pr n The refractive index n is usually slightly larger than 1 Electromagnetic waves travel slower in dielectrics The wave impedance is the ratio of the field amplitudes: Z = E=H in units of Ω = V=A In vacuo the impedance is a constant: Z0 = µ0c = 377Ω In an insulator the impedance is: µrµ0c µr Z = µrµ0v = = Z0 prµr r r For non-magnetic materials with µr = 1: Z = Z0=n 3 Notes: Diagrams: 4 Energy Propagation in Insulators The Poynting vector N = E H measures the energy flux × Energy flux is energy flow per unit time through surface normal to direction of propagation of wave: @U −2 @t = A N:dS Units of N are Wm In vacuo the amplitudeR of the Poynting vector is: 1 1 E2 N = E2 0 = 0 0 2 0 µ 2 Z r 0 0 In an insulator this becomes: 1 E2 N = N r = 0 0 µ 2 Z r r The energy flux is proportional to the square of the amplitude, and inversely proportional to the wave impedance 5 Dispersion Dispersion occurs because the dielectric constant r and refractive index
    [Show full text]
  • 1 Vectors and Vector Analysis 2 Electrostatics
    1 Vectors and Vector Analysis ² Scalar and Vector Products: ¡ ¢ a ¢ b = axbx + ayby + azbz ; a £ b = ay bz ¡ az by ; az bx ¡ axbz; axby ¡ aybx : a ¢ b is a scalar, a £ b is a pseudovector. ² Multiple Products: a ¢ (b £ c) = (a £ b) ¢ c ; a £ (b £ c) = b(a ¢ c) ¡ c(a ¢ b) : ² Scalar and vector operators: 2 2 2 ¡ @@@ ¢ @ @ @ r = ;; ; ¢ = r ¢ r = 2 +2 + 2 : @x @y @z @x @y @z ² Gradient, divergence, and curl of ¯elds: ¡ @f @f @f ¢ grad f(r) = rf(r) = ;; ; @x @y @z @ax @ay @az div a(r) = r ¢ a(r) = + + ; @x @y @z ¡ @az @ay @ax @az @ay @ax ¢ rot a(r) = r £ a(r) = ¡ ; ¡ ; ¡ : @y @z @z @x @x @y ² Second derivatives: ¡ ¢ rot grad f(r) = r £ rf(r) = o ; ¡ ¢ div rot a(r) = r ¢ r £ a(r) = 0 ; ¡ ¢ div grad f(r) = r ¢ rf(r) = ¢f(r) ; ¡ ¢ ¡ ¢ ¡ ¢ rot rot a(r) = r £ r £ a(r) = r r ¢ a(r) ¡ r £ r a(r) = grad div a(r) ¡ ¢a(r) : ² Integral theorems of Gauss and Stokes: If ­ is a ¯nite volume, and @­ its closed surface, I I ZZZ £ ¤ 3 b(r) ¢ da(r) = b(r) ¢ n^(r) da(r) = div b(r) d r : @­ @­ ­ If A is a ¯nite surface, and @A is its closed boundary, I ZZ £ ¤ b(r) ¢ dr = rot b(r) ¢ n^(r) da(r) : @A A 2 Electrostatics ² Elementary charge e and dielectric susceptibility ²0 : ¡19 e = 1:60217733 ¢ 10 As ; ¡12 As ²0 = 8:854187818 ¢ 10 : Vm 0 ² Coulomb's law (point particle, charge q, at position r ): 0 q r ¡ r E(r) = 0 3 : 4¼²0 jr ¡ r j 1 ² Coulomb's law (charge distribution ½(r)): ZZZ 0 0 1 ½(r )(r ¡ r ) 3 0 E(r) = 0 3 d r : 4¼²0 jr ¡ r j ² Coulomb's law in restricted geometries (surface charge density σ(r), line charge density ¸(r)): ZZ 0 0 Z 0 0 1 σ(r )(r ¡ r )
    [Show full text]
  • The Maxwell Stress Tensor and Electromagnetic Momentum
    Progress In Electromagnetics Research Letters, Vol. 94, 151–156, 2020 The Maxwell Stress Tensor and Electromagnetic Momentum Artice Davis1, * and Vladimir Onoochin2 Abstract—In this paper, we discuss two well-known definitions of electromagnetic momentum, ρA and 0[E × B]. We show that the former is preferable to the latter for several reasons which we will discuss. Primarily, we show in detail—and by example—that the usual manipulations used in deriving the expression 0[E×B] have a serious mathematical flaw. We follow this by presenting a succinct derivation for the former expression. We feel that the fundamental definition of electromagnetic momentum should rely upon the interaction of a single particle with the electromagnetic field. Thus, it contrasts with the definition of momentum as 0[E × B] which depends upon a (defective) integral over an entire region, usually all space. 1. INTRODUCTION The correct form of the electromagnetic energy-momentum tensor has been debated for almost a century. A large number of papers have been written [1], some of which concentrate upon a single particle interaction with the electromagnetic field in free space [2, 3]; others discuss the interaction between a field and a dielectric material upon which it impinges [4]. In works which discuss the Minkowski-Abraham controversy, the system considered is usually on the macroscopic level. On this level, singularities caused by point charges cannot be treated. We show that neglecting the effect due the properties of classical charges creates difficulties which cannot be resolved within the framework typically discussed by textbooks on classical electrodynamics [5]. Our present concept of electromagnetic field momentum density has two interpretations: (i) As 0[E × B], the scaled Poynting vector.
    [Show full text]
  • Chapter 5 Energy and Momentum
    Chapter 5 Energy and Momentum The equations established so far describe the behavior of electric and magnetic fields. They are a direct consequence of Maxwell’s equations and the properties of matter. Although the electric and magnetic fields were initially postulated to explain the forces in Coulomb’s and Ampere’s laws, Maxwell’s equations do not provide any information about the energy content of an electromagnetic field. As we shall see, Poynting’s theorem provides a plausible relationship between the electromagnetic field and its energy content. 5.1 Poynting’s Theorem If the scalar product of the field E and Eq. (1.33) is subtracted from the scalar product of the field H and Eq. (1.32) the following equation is obtained ∂B ∂D H ( E) E ( H) = H E j E . (5.1) · ∇ × − · ∇ × − · ∂t − · ∂t − 0 · Noting that the expression on the left is identical to (E H), integrating both ∇· × sides over space and applying Gauss’ theorem the equation above becomes ∂B ∂D (E H) n da = H + E + j0 E dV (5.2) ∂V × · − V · ∂t · ∂t · Z Z h i Although this equation already forms the basis of Poynting’s theorem, more insight is provided when B and D are substituted by the generally valid equations (1.19). 59 60 CHAPTER 5. ENERGY AND MOMENTUM Equation (5.2) then reads as ∂ 1 (E H) n da + D E + B H dV = j0 E dV (5.3) ∂V × · ∂t V 2 · · − V · Z Z h i Z 1 ∂P ∂E µ ∂M ∂H E P dV 0 H M dV.
    [Show full text]