磁單極子存在時的Maxwell方程組的初邊值問題

王治蓉, 楊丕文*, 李又超

(1. 四川師范大學 數學與軟件科學學院, 四川 成都 610066; 2. 東方電氣集團 東方電機有限公司, 四川 德陽 618000)

Maxwell方程組是電磁理論的核心,描述了宏觀電磁現象,方程組的微分形式為

(1)

方程組(1)的4個方程表示的物理意義在文獻[1]中有詳細的解釋.此方程組表明自然界中存在著作為場源的電荷和電流,不存在磁荷與磁流,即磁單極子不存在.理論和實驗卻說明磁單極子存在,國內外都有一些關于磁單極子的研究[2-9].

如果磁單極子存在,設磁荷密度為g0,磁流密度為K,與處理電荷的方法類似,磁流密度K與磁荷密度g0滿足連續性方程

(2)

方程組(1)中的第一與第二式可改寫為

▽·B=g0,

(3)

(4)

方程組(1)中第三與第四式在磁荷存在時仍然成立.

因此,磁單極子存在時,得到一組新的Maxwell方程組為

(5)

設C=1,ε0=μ0=1,E=Φ,H=Ψ, Maxwell方程組(1)與方程組(5)即可簡化為[10-11]

(6)

(7)

引入擬四元數空間[12],即記二階矩陣

則有

e2e3=-e3e2=-ie1,e3e1=-e1e3=-ie2.

引入一階微分算子

▽,

在文獻[12]中討論了某些一階雙曲方程組的初邊值問題,用類似的方法,在文獻[13]中討論了Maxwell方程組(6)的初邊值問題,而在文獻[12]中引入了擬四元數空間后,Maxwell方程組(6)可表示為

D(ψ+iΦ)=if2,

(8)

而方程組(7)可表示為

D(ψ+iΦ)=f1+if2,

(9)

而方程組(7)是磁單極子存在時的Maxwell方程組,下面就討論此方程組的初邊值問題.為了簡化計算,先討論方程

▽)(ψ+iΦ)=-K+i(ρ-J)

2 方程▽)(ψ+iΦ)=-K+i(ρ-J)的初邊值問題

設G是R4中的一個區域,考慮G上關于(ψ+iΦ)的方程

▽)(ψ+iΦ)=-K+i(ρ-J), (10)

(11)

記Ω={(t,x)||x|<1-t,0

(φ2e0+iφ3e1)(0,x)=f2e0+if3e1,

(12)

φ1(0,0,x2,x3)=f1(x2,x3),

(13)

ψ1(0,x)=τ(x),x∈?B3,

(14)

Re(x2e0-ix3e1)κ[ψ2(0,0,x2,x3)e0+

iψ3(0,0,x2,x3)e1]=r(x2,x3)e0,

(15)

其中

ω(x)=ω1+(ω2e0+iω3e1)e2,

(16)

(17)

2)當κ< 0時,問題可解,當且僅當r′(ζ)滿足條件

m=1,2,3,…,-κ-1.

(18)

當此條件滿足時,其解仍如結論(1)中式子表示,但其中ψ(t,x)中的ω(x)=ω1e0+(ω2e0+iω3e1)e2,ω1與結論(1)中表示一致,而

(19)

從而可以推出方程(10)及方程組(11)的相容性條件

▽·J=0,

(20)

方程(20)在電動力學中被稱為電荷守恒方程[14],以下恒假設ρ及J滿足此方程.

定理1對錐形區域ω上的滿足方程組

(21)

的向量函數Φ(t,x)的初始條件(12)和(13)的Cauchy問題,有且僅有唯一的解

(22)

其中

(24)

(25)

(26)

證明由(23)和(24)式知,φ2,φ3滿足方程及初始條件[15]

(27)

(28)

而C(t,x2,x3)滿足方程

(29)

由(25)式與(27)、(28)式,可得

則

從而

即得φ1(t,x)滿足方程(30).

將方程(10)化成關于向量函數ψ的方程組

(31)

由此可得定理2.

(32)

其中ω(x)是滿足方程▽ω=0的任一向量函數.

τ″(x),x∈?B,

Re(x2e0-ix3e1)κ[ω2(0,0,x2,x3)e0+iω3(0,0,x2,x3)e1]=

r(x2,x3)e0-im(x2e0-ix3e1)κ×

r″(x2,x3)e0, (x2,x3)∈Г .

定理3(1)當κ≥ 0時,Ω上的方程(10)的問題F有解ψ+iΦ,其中Φ(t,x)如(22)式所示,ψ(t,x)如式(32)所示,ω(x)類似于式(16)和式(17)表示,僅分別用τ″(x)和r″(x2,x3)代替其中的τ′(x)和r′(x2,x3).

(2)當κ< 0時,Ω上的方程(10)的問題F可解,當且僅當r″(ζ)如r′(ζ)滿足條件(18),當這些條件滿足時,其解ψ+iΦ仍由(22)式與(32)式表示,ψ(t,x)中的ω(x)由τ″(x)和r″(x2,x3)代替(19)式中的τ′(x)和r′(x2,x3).

3 含磁單極子的Maxwell方程組的初邊值問題

設G是R4中的一個區域,考慮G上關于(ψ+iΦ)的方程

▽)(ψ+iΦ)=f1+if2,

(33)

(34)

(34)式即為磁單極子存在時的Maxwell方程組.

從而可以推出方程(33)及方程組(34)中的g0、K、ρ、J滿足條件

▽·K=0,

(35)

(36)

方程(35)和(36)在電動力學中被稱為磁荷守恒方程和電荷守恒方程,以下恒假設g0、K、ρ、J滿足方程(35)和(36).對錐形區域Ω上的滿足方程組

的向量函數Φ(t,x)的初始條件(12)和(13)的Cauchy問題,有且僅有唯一的解

(38)

其中φ2、φ3、φ1如定理1的(23)～(26)式.

將方程(33)化成關于向量函數ψ的方程組

由此可得定理4.

定理4關于Ω上的向量函數ψ(t,x)的方程組(39)可解的充要條件是Φ(t,x)滿足(21)式中的

當此條件滿足時,其一般解可表示為

▽×Φ)(ζ,x)dζ+T3[g0-

(40)

其中ω(x)是滿足方程▽ω=0的任一向量函數.

τ?(x),x∈?B,

Re(x2e0-ix3e1)κ[ω2(0,0,x2,x3)e0+

iω3(0,0,x2,x3)e1]=r(x2,x3)e0-im(x2e0-ix3e1)κ×

r?(x2,x3)e0, (x2,x3)∈Г.

定理51) 當κ≥ 0時,Ω上的方程(33)的問題F有解ψ+iΦ,其中Φ(t,x)如(38)式所示,ψ(t,x)如(40)式所示,ω(x)類似于(16)和(17)式表示,僅分別用τ?(x)和r?(x2,x3)代替其中的τ′(x)和r′(x2,x3).

2) 當κ<0時,Ω上的方程(33)的問題F可解,當且僅當r?(ζ)如r′(ζ)滿足條件(18),當這些條件滿足時,其解ψ+iΦ仍由(38)與(40)式表示,ψ(t,x)中的ω(x)由τ?(x)和r?(x2,x3)代替 (19)式中的τ′(x)和r′(x2,x3).

磁單極子存在時,電磁關系發生改變,定理5的結論即為磁單極子存在時的Maxwell方程組的初邊值問題在不同情況下的解,其中的Φ、ψ即為電場強度和磁場強度函數.

[1] Maxwell, Clerk J. A dynamical theory of the electromagnetic field[J]. Philosophical Transactions of the Royal Society of London,1865,155:459-512.

[2] 劉光華,鄧小燕. 磁單極子進展概述[J]. 大學物理,2010,29(11):1-5.

[3] 胡永紅,毛彩霞. 旋轉磁單極子場方程的研究[J]. 華中師范大學學報:自然科學版,2012,46(2):152-154.

[4] 胡陳果. 推導包括磁荷的麥氏方程組的新方法[J]. 四川師范大學學報:自然科學版,1994,15(1):89-93.

[5] 宮蒂. 含磁荷磁流麥氏方程組解的特性[J]. 中山大學學報:自然科學版,1995,34(1):106-109.

[6] Morris D J P, Tennant D A, Grigera S A, et al. Dirac strings and magnetic monopoles in spin Ice Dy2Ti2O7[J]. Science,2009,326:411-414.

[7] Giblin S R, Bramwell S T, Holdsworth P C W, et al. Creation and measurement of long-lived magnetic monopole currents in spin ice[J]. Nature Physics,2011,7:252-258.

[8] Bonitz M. Spin ice:a plasma of magnetic monopoles[J]. Nature Physics,2011,7:192-194.

[9] Mengotti E, Heyderman L J, Rodríguez A F, et al. Real-space observation of emergent magnetic monopoles and associated Dirac strings in artificial kagome spin ice[J]. Nature Physics,2010,7:68-74.

[10] 陳池一,李康. Maxwell理論的電磁對偶性[J]. 浙江大學學報:自然科學版,2001,28(1):27-34.

[11] 盧里·D. 粒子和場[M]. 董明德,譯. 北京:科學出版社,1981.

[12] 楊丕文,李曼荔,陳穎. 某些一階雙曲方程組的初邊值問題[J]. 應用數學學報,2008,31(1):61-70.

[13] Yang P W, Yang S, Li M L. An initial-boundary value problem for the maxwell equations[J]. J Diff Eqns,2010,249:3003-3023.

[14] 郭碩鴻. 電動力學[M]. 北京:高等教育出版社,1984.

[15] 柯朗, 希爾伯特. 數學物理方程[M]. 北京:科學出版社,1977

[16] Yang P W. The Riemann-Hilbert boundary value problem for the Moisil-Theodorsco system[J]. Acta Math Sci,2006,26 (7):1057-1063.

[17] 楊丕文,楊碩. 可換四元數空間中某些雙曲型方程的Riemann-Hilbert邊值問題[J]. 數學學報,2008,51(1):171-180.

[18] 楊丕文,李曼荔,楊碩. 在可交換四元數空間中的雙曲型方程的特征邊值問題[J]. 數學學報,2007,50(6):1249-1256.

[19] 楊丕文. 四元數正則函數的某些函數論性質[J]. 應用數學學報,2010,33(6):1095-1112.

[20] 楊丕文,李曼荔.n階Moisil-Theodoresc型方程的Riemann-Hilbert邊值問題[J]. 數學學報,2010,53(3):585-596.

[21] 蒲松,楊丕文. 雙解析函數的非正則型及非齊次二階方程的某些邊值問題[J]. 四川師范大學學報:自然科學版,2007,30(3):295-299.

[22] 楊丕文. 雙圓柱區域上的高階四元數函數的Dirichlet邊值問題[J]. 中國科學:數學,2011,41(6);485-496.

[23] 李長江,楊丕文. 擬四元數空間中的一些邊值問題[J]. 四川師范大學學報:自然科學版,2010,33(03):286-291.

[24] Yang P W, Li D.Hλ-regular vector functions and their boundary value problems[J]. Boundary Value Problems,2012,2012:75.

[25] 李覺友. Moisil-Theodorsco方程組的一個非線性邊值問題[J]. 重慶師范大學學報:自然科學版,2008,25(1):42-46.

[26] 李覺友. 非齊次Moisil-Theodorsco方程組的Riemann邊值問題[J]. 重慶師范大學學報:自然科學版,2007,24(4):27-29.