古典電磁気学の定理
古典電磁気学 において 、 相互性は、一定の制約の下で時間不変の線形媒体に対する マクスウェル方程式 における時間 調和 電流 密度 (発生源)と結果として生じる 電磁場 の交換を含む 、さまざまな関連定理を指します。相互性は、電磁気学に適用される
線型代数 の対称 演算子 の概念と密接に関連しています。
おそらく最も一般的で一般的な定理は ローレンツの相反性 (および レイリー・カーソンの相反性 などのさまざまな特殊なケース)であり、 レイリー卿 による 音 と ヘルムホルツ による 光 に関する類似の結果を受けて1896年に ヘンドリック・ローレンツ が行った研究にちなんで名付けられました(Potton 2004)。大まかに言えば、電流が流される点と電場が測定される点を入れ替えても、 振動電流と結果として生じる 電場の関係は変わらないと述べています。 電気ネットワークの特定のケースでは、ネットワーク内の異なるポイントでの 電圧 と 電流は 交換可能であるというステートメントとして表現されることがあります 。より技術的には、2番目の回路による1番目の回路の相互 インピーダンス は、1番目の回路による2番目の回路の相互インピーダンスと同じになります。
相互性は光学 において有用であり 、(量子効果とは別に)古典的電磁気学だけでなく 放射測定 法でも表現できます。
静電気学 にも、 電位 と 電荷密度 の交換に関する グリーンの相反性 として知られる類似の定理があります 。
相互性定理の形式は、電気ネットワークや アンテナ システムの解析など、多くの電磁気学の応用で使用されています。 [1]
たとえば、相互性は、アンテナが送信機としても受信機としても同様に機能すること、具体的には、アンテナの 放射パターンと受信パターン が同一であることを意味します。相互性は、 インピーダンス行列 と 散乱行列の対称性、 境界要素および転送行列の計算方法で使用する グリーン関数 の対称性、 導波路 システムにおける 高調波モード の 直交性特性(これらの特性を 固有演算子 の対称性から直接証明する代わりに) など、電磁気システム に関する他の定理を証明するために使用される基本的な補題でもあります 。
ローレンツの相互性
具体的には、電流密度があり、 電場 と 磁場を 生成するとします 。この 3 つはすべて、 角周波数 ω を持つ時間の周期関数であり、特に時間依存性があります。同様に 、同じ周波数 ω で2 番目の電流があり 、それ自体で電場 と磁場を生成するとします。ローレンツの相反定理は、以下で説明する媒体の材料に関する特定の単純な条件下で、 体積 V を 囲む任意の表面 S について、次のことを述べています。
J
1
{\displaystyle \mathbf {J} _{1}}
E
1
{\displaystyle \mathbf {E} _{1}}
H
1
,
{\displaystyle \mathbf {H} _{1}\,,}
exp
(
−
i
ω
t
)
.
{\displaystyle \exp(-i\omega t)\,.}
J
2
{\displaystyle \mathbf {J} _{2}}
E
2
{\displaystyle \mathbf {E} _{2}}
H
2
.
{\displaystyle \mathbf {H} _{2}\,.}
∫
V
[
J
1
⋅
E
2
−
E
1
⋅
J
2
]
d
V
=
∮
S
[
E
1
×
H
2
−
E
2
×
H
1
]
⋅
d
S
.
{\displaystyle \int _{V}\left[\mathbf {J} _{1}\cdot \mathbf {E} _{2}-\mathbf {E} _{1}\cdot \mathbf {J} _{2}\right]\mathrm {d} V=\oint _{S}\left[\mathbf {E} _{1}\times \mathbf {H} _{2}-\mathbf {E} _{2}\times \mathbf {H} _{1}\right]\cdot \mathbf {\mathrm {d} S} \ .}
同様に、微分形式では( 発散定理 により)次のようになります。
J
1
⋅
E
2
−
E
1
⋅
J
2
=
∇
⋅
[
E
1
×
H
2
−
E
2
×
H
1
]
.
{\displaystyle \mathbf {J} _{1}\cdot \mathbf {E} _{2}-\mathbf {E} _{1}\cdot \mathbf {J} _{2}=\nabla \cdot \left[\mathbf {E} _{1}\times \mathbf {H} _{2}-\mathbf {E} _{2}\times \mathbf {H} _{1}\right]\ .}
この一般形は、多くの特殊なケースでは簡略化されることが多い。特に、通常は と が 局所化されている(つまり、 コンパクトなサポート を持つ)こと、また無限遠方から入射する波がないことが想定される。この場合、空間全体で積分すると、面積分項が打ち消され(以下を参照)、次式が得られる。
J
1
{\displaystyle \ \mathbf {J} _{1}\ }
J
2
{\displaystyle \mathbf {J} _{2}}
∫
J
1
⋅
E
2
d
V
=
∫
E
1
⋅
J
2
d
V
.
{\displaystyle \int \mathbf {J} _{1}\cdot \mathbf {E} _{2}\,\mathrm {d} V=\int \mathbf {E} _{1}\cdot \mathbf {J} _{2}\,\mathrm {d} V\ .}
この結果(および以下の簡略化)は、 レイリー卿の音波に関する研究と カーソン (1924年、1930年)による無線周波数アンテナへの応用への拡張にちなんで、 レイリー-カーソン相反定理 と呼ばれることもあります。多くの場合、この関係は点状の 双極子源を考慮することでさらに簡略化されます。その場合、積分は消え、単に電界と電流の対応する 双極子 モーメントの積が得られます。または、無視できるほど薄い電線の場合は、1本の電線に印加された電流と別の電線にかかる電圧を掛け合わせたものが得られ、その逆も同様です。以下も参照してください。
ローレンツの相互定理の別の特殊なケースは、体積 V が 両方 の局所的なソースを完全に含む場合 (または V が どちら のソースとも交差しない場合 ) に適用されます。この場合、
∮
S
(
E
1
×
H
2
)
⋅
d
S
=
∮
S
(
E
2
×
H
1
)
⋅
d
S
.
{\displaystyle \ \oint _{S}(\mathbf {E} _{1}\times \mathbf {H} _{2})\cdot \mathbf {\mathrm {d} S} =\oint _{S}(\mathbf {E} _{2}\times \mathbf {H} _{1})\cdot \mathbf {\mathrm {d} S} \ .}
実際の問題では、ローレンツ関係やその他の相互関係のより一般化された形式があり、 電流密度 に加えて 磁気電流密度も使用されます。これらのタイプの相互関係は通常、 電気工学の 文献 で議論されています。 [2] [3] [4] [5] [6] [7]
J
{\displaystyle \ \mathbf {J} \ }
M
{\displaystyle \ \mathbf {M} \ }
電気ネットワークの相互性
上記では、ローレンツの相反性は、外部から印加される電流源とその結果生じる場という観点から説明されました。多くの場合、特に電気ネットワークでは、代わりに外部から印加される電圧とその結果生じる電流について考えることが好まれます。ローレンツの相反性定理は、このケースも説明しており、 対称 であることが求められる3×3 導電率 マトリックス σを持つ オーム性材料 (つまり、印加された場に対して線形に反応する電流)を想定しています 。これは、以下の他の条件によって暗示されています。この状況を適切に説明するには、外部から印加される場 (駆動電圧から) と、結果として生じる全場とを注意深く区別する必要があります ( King 、 1963)。
より具体的には、 上記はマクスウェル方程式に導入された外部「ソース」項のみで構成されていました。これを、外部ソースと材料内の結果として生じる電界の両方によって生成される 合計 電流と区別するために、 と表記します。この外部電流が導電率 σ の材料内にある場合 、それは外部から印加される電界に対応し 、 σ の定義により次のようになります。
J
{\displaystyle \ \mathbf {J} \ }
J
(
e
)
{\displaystyle \ \mathbf {J} ^{(e)}\ }
E
(
e
)
{\displaystyle \ \mathbf {E} ^{(e)}\ }
J
(
e
)
=
σ
E
(
e
)
.
{\displaystyle \ \mathbf {J} ^{(e)}=\sigma \mathbf {E} ^{(e)}\ .}
さらに、上記の電場はこの電流に対する 応答 のみで構成されており 、「外部」の場は含まれていません。 したがって、前の場を次のように表します。 ここで、 全体の 場は次のように表されます。
E
{\displaystyle \mathbf {E} }
E
(
e
)
.
{\displaystyle \ \mathbf {E} ^{(e)}\ .}
E
(
r
)
,
{\displaystyle \ \mathbf {E} ^{(r)}\ ,}
E
=
E
(
e
)
+
E
(
r
)
.
{\displaystyle \ \mathbf {E} =\mathbf {E} ^{(e)}+\mathbf {E} ^{(r)}\ .}
さて、ローレンツの相互定理の左側の式は、 σ を外部電流項から 応答場項に移動し 、項を加算および減算することで書き直すことができ 、外部場と 全 電流の積が得られる。
J
(
e
)
{\displaystyle \mathbf {J} ^{(e)}}
E
(
r
)
,
{\displaystyle \ \mathbf {E} ^{(r)}\ ,}
σ
E
1
(
e
)
E
2
(
e
)
{\displaystyle \ \sigma \mathbf {E} _{1}^{(e)}\mathbf {E} _{2}^{(e)}\ }
J
=
σ
E
:
{\displaystyle \ \mathbf {J} =\sigma \mathbf {E} \ :}
∫
V
[
J
1
(
e
)
⋅
E
2
(
r
)
−
E
1
(
r
)
⋅
J
2
(
e
)
]
d
V
=
∫
V
[
σ
E
1
(
e
)
⋅
(
E
2
(
r
)
+
E
2
(
e
)
)
−
(
E
1
(
r
)
+
E
1
(
e
)
)
⋅
σ
E
2
(
e
)
]
d
V
=
∫
V
[
E
1
(
e
)
⋅
J
2
−
J
1
⋅
E
2
(
e
)
]
d
V
.
{\displaystyle {\begin{aligned}&\int _{V}\left[\mathbf {J} _{1}^{(e)}\cdot \mathbf {E} _{2}^{(r)}-\mathbf {E} _{1}^{(r)}\cdot \mathbf {J} _{2}^{(e)}\right]\operatorname {d} V\\={}&\int _{V}\left[\sigma \mathbf {E} _{1}^{(e)}\cdot \left(\mathbf {E} _{2}^{(r)}+\mathbf {E} _{2}^{(e)}\right)-\left(\mathbf {E} _{1}^{(r)}+\mathbf {E} _{1}^{(e)}\right)\cdot \sigma \mathbf {E} _{2}^{(e)}\right]\operatorname {d} V\\={}&\int _{V}\left[\mathbf {E} _{1}^{(e)}\cdot \mathbf {J} _{2}-\mathbf {J} _{1}\cdot \mathbf {E} _{2}^{(e)}\right]\operatorname {d} V\ .\end{aligned}}}
細いワイヤの極限では、これは外部から印加された電圧 (1) と結果として生じる全電流 (2) の積となり、その逆も同様である。特に、レイリー・カーソンの相反定理は単純な合計となる。
∑
n
V
1
(
n
)
I
2
(
n
)
=
∑
n
V
2
(
n
)
I
1
(
n
)
{\displaystyle \ \sum _{n}{\mathcal {V}}_{1}^{(n)}I_{2}^{(n)}=\sum _{n}{\mathcal {V}}_{2}^{(n)}I_{1}^{(n)}}
ここで 、および Iは、 2つの可能な電圧セット と電圧セットに対する 回路要素のセット( nでインデックス付け)における AC 印加電圧と結果として生じる電流の 複素振幅を それぞれ表します。
V
{\displaystyle \ {\mathcal {V}}\ }
V
1
{\displaystyle \ {\mathcal {V}}_{1}\ }
V
2
.
{\displaystyle \ {\mathcal {V}}_{2}\ .}
最も一般的には、これは各システムが単一の 電圧源 を持つ場合にさらに簡略化され 、 その 場合定理は単純に
V
s
,
{\displaystyle \ {\mathcal {V}}_{\text{s}}\ ,}
V
1
(
1
)
=
V
s
{\displaystyle \ {\mathcal {V}}_{1}^{(1)}={\mathcal {V}}_{\text{s}}\ }
V
2
(
2
)
=
V
s
.
{\displaystyle \ {\mathcal {V}}_{2}^{(2)}={\mathcal {V}}_{\text{s}}\ .}
I
1
(
2
)
=
I
2
(
1
)
{\displaystyle I_{1}^{(2)}=I_{2}^{(1)}}
または言葉で言うと:
(2)の電圧からの位置(1)の電流は、(1)の同じ電圧からの位置(2)の電流と同一である。
ローレンツの相互性の条件と証明
ローレンツの相反定理は、固定周波数 (線型媒体)
で と を 関連付ける 線型演算子が次の式を満たすという事実を反映しているだけです。
ここで は
、
通常、 ベクトル場 と [8] の 「 内積 」の下で 対称演算子 です
(技術的には、この 非共役 形式は、複素数値場に対しては実数値ではないため、真の内積ではありませんが、ここでは問題になりません。この意味で、演算子は真のエルミート型ではなく、複素対称型です)。これは、 与えられた ωにおける 誘電率 ε と 透磁率 μが 対称 3×3 行列(対称階数 2 テンソル)である場合に常に当てはまります。もちろん、これらが スカラー である一般的なケース (等方性媒体の場合)も含まれます。これらは実数である必要は ありません。複素数値は、有限の導電率 σ を持つ導体( を介して ε に含まれる) などの損失のある材料に対応しており、 このため、相互定理は 時間反転不変性 を必要としませ ん。対称 ε 行列と μ 行列の条件は 、ほぼ常に満たされます。例外については以下を参照してください。
O
^
{\displaystyle \operatorname {\hat {O}} }
J
{\displaystyle \mathbf {J} }
E
{\displaystyle \mathbf {E} }
ω
{\displaystyle \omega }
J
=
O
^
E
{\displaystyle \mathbf {J} =\operatorname {\hat {O}} \mathbf {E} }
O
^
E
≡
1
i
ω
[
1
μ
(
∇
×
∇
×
)
−
ω
2
ε
]
E
{\displaystyle \operatorname {\hat {O}} \mathbf {E} \equiv {\frac {1}{i\omega }}\left[{\frac {1}{\mu }}\left(\nabla \times \nabla \times \right)-\;\omega ^{2}\varepsilon \right]\mathbf {E} }
(
F
,
G
)
=
∫
F
⋅
G
d
V
{\textstyle (\mathbf {F} ,\mathbf {G} )=\int \mathbf {F} \cdot \mathbf {G} \,\mathrm {d} V}
F
{\displaystyle \mathbf {F} }
G
.
{\displaystyle \mathbf {G} \ .}
ε
→
ε
+
i
σ
/
ω
{\displaystyle \varepsilon \rightarrow \varepsilon +i\sigma /\omega \ }
内積 の下での任意の エルミート演算子に対して 、定義により が成り立ち 、レイリー・カーソンの相互定理は、この特定の演算子に対するこのステートメントのベクトル版に過ぎず、 となります 。 ここでの演算子のエルミート特性は、 部分積分 によって導くことができます。有限積分体積に対して、この部分積分からの表面項は、上記のより一般的な表面積分定理をもたらします。特に、重要な事実は、ベクトル場 と、 表面 S で囲まれた 体積 V 上の部分積分 (または 発散定理 ) に対して、次の 恒等式を与えることです。
O
^
{\displaystyle \operatorname {\hat {O}} }
(
f
,
g
)
{\displaystyle (f,g)\!}
(
f
,
O
^
g
)
=
(
O
^
f
,
g
)
{\displaystyle (f,\operatorname {\hat {O}} g)=(\operatorname {\hat {O}} f,g)}
J
=
O
^
E
:
{\displaystyle \mathbf {J} =\operatorname {\hat {O}} \mathbf {E} \ :}
(
E
1
,
O
^
E
2
)
=
(
O
^
E
1
,
E
2
)
.
{\displaystyle (\mathbf {E} _{1},\operatorname {\hat {O}} \mathbf {E} _{2})=(\operatorname {\hat {O}} \mathbf {E} _{1},\mathbf {E} _{2})\ .}
F
{\displaystyle \mathbf {F} }
G
,
{\displaystyle \mathbf {G} \ ,}
∫
V
F
⋅
(
∇
×
G
)
d
V
≡
∫
V
(
∇
×
F
)
⋅
G
d
V
−
∮
S
(
F
×
G
)
⋅
d
A
.
{\displaystyle \int _{V}\mathbf {F} \cdot (\nabla \times \mathbf {G} )\,\mathrm {d} V\equiv \int _{V}(\nabla \times \mathbf {F} )\cdot \mathbf {G} \,\mathrm {d} V-\oint _{S}(\mathbf {F} \times \mathbf {G} )\cdot \mathrm {d} \mathbf {A} \ .}
この恒等式を 2 回適用すると、 プラス面項
が得られ、ローレンツの相互関係が得られます。
(
E
1
,
O
^
E
2
)
{\displaystyle (\mathbf {E} _{1},\operatorname {\hat {O}} \mathbf {E} _{2})}
(
O
^
E
1
,
E
2
)
{\displaystyle (\operatorname {\hat {O}} \mathbf {E} _{1},\mathbf {E} _{2})}
マクスウェル方程式とベクトル演算を用いたローレンツの相互性の条件と証明 [9]
ローレンツの電磁相互作用定理の一般形を証明します。これは 、 2つの異なる正弦波電流密度によってそれぞれ生成される 同じ周波数の
場 と場が、次の条件を満たすというものです。
E
1
,
H
1
{\displaystyle \mathbf {E} _{1},\mathbf {H} _{1}}
E
2
,
H
2
{\displaystyle \mathbf {E} _{2},\mathbf {H} _{2}}
J
1
{\displaystyle \mathbf {J} _{1}}
J
2
{\displaystyle \mathbf {J} _{2}}
∫
V
[
J
1
⋅
E
2
−
E
1
⋅
J
2
]
d
V
=
∮
S
[
E
1
×
H
2
−
E
2
×
H
1
]
⋅
d
S
.
{\displaystyle \int _{V}\left[\mathbf {J} _{1}\cdot \mathbf {E} _{2}-\mathbf {E} _{1}\cdot \mathbf {J} _{2}\right]\mathrm {d} V=\oint _{S}\left[\mathbf {E} _{1}\times \mathbf {H} _{2}-\mathbf {E} _{2}\times \mathbf {H} _{1}\right]\cdot \mathbf {\mathrm {d} S} .}
誘電率と透磁率が位置の関数ではあるが時間の関数ではない領域を考えてみましょう。領域の全フィールド、電流、電荷で記述されるマクスウェル方程式は、領域の電磁気的挙動を記述します。2 つの回転方程式は次のとおりです。
∇
×
E
=
−
∂
∂
t
B
,
∇
×
H
=
J
+
∂
∂
t
D
.
{\displaystyle {\begin{array}{ccc}\nabla \times \mathbf {E} &=&-{\frac {\partial }{\partial t}}\mathbf {B} \ ,\\\nabla \times \mathbf {H} &=&\mathbf {J} +{\frac {\partial }{\partial t}}\mathbf {D} \ .\end{array}}}
安定した一定周波数条件下では、2 つの回転方程式から時間周期の場合の Maxwell 方程式が得られます。
∇
×
E
=
−
j
ω
B
,
∇
×
H
=
J
+
j
ω
D
.
{\displaystyle {\begin{array}{ccc}\nabla \times \mathbf {E} &=&-j\omega \mathbf {B} \ ,\\\nabla \times \mathbf {H} &=&\mathbf {J} +j\omega \mathbf {D} \ .\end{array}}}
この記事の方程式の記号は の複素乗数を表し 、選択された基準に対する同相部分と逆相部分を表すことを認識する必要があります。 の複素ベクトル乗数は、一般に 位相器 と呼ばれる複素スカラー量と同様に、 ベクトル位相器 と呼ばれることがあります 。
e
j
ω
t
{\displaystyle e^{j\omega t}}
e
j
ω
t
{\displaystyle e^{j\omega t}}
ベクトル演算の同値性は、
すべてのベクトル と
H
⋅
(
∇
×
E
)
−
E
⋅
(
∇
×
H
)
=
∇
⋅
(
E
×
H
)
{\displaystyle \mathbf {H} \cdot (\nabla \times \mathbf {E} )-\mathbf {E} \cdot (\nabla \times \mathbf {H} )=\nabla \cdot (\mathbf {E} \times \mathbf {H} )}
E
{\displaystyle \mathbf {E} }
H
.
{\displaystyle \mathbf {H} \ .}
この等価性を に適用する と 次のようになります。
E
1
{\displaystyle \mathbf {E} _{1}}
H
2
{\displaystyle \mathbf {H} _{2}}
H
2
⋅
(
∇
×
E
1
)
−
E
1
⋅
(
∇
×
H
2
)
=
∇
⋅
(
E
1
×
H
2
)
.
{\displaystyle \mathbf {H} _{2}\cdot (\nabla \times \mathbf {E} _{1})-\mathbf {E} _{1}\cdot (\nabla \times \mathbf {H} _{2})=\nabla \cdot (\mathbf {E} _{1}\times \mathbf {H} _{2})\ .}
時間周期方程式の積をこの最後の同値で示されるように取って足すと、
−
H
2
⋅
j
ω
B
1
−
E
1
⋅
j
ω
D
2
−
E
1
⋅
J
2
=
∇
⋅
(
E
1
×
H
2
)
.
{\displaystyle -\mathbf {H} _{2}\cdot j\omega \mathbf {B} _{1}-\mathbf {E} _{1}\cdot j\omega \mathbf {D} _{2}-\mathbf {E} _{1}\cdot \mathbf {J} _{2}=\nabla \cdot (\mathbf {E} _{1}\times \mathbf {H} _{2})\ .}
これは現在、懸念される量にわたって統合される可能性がある。
∫
V
(
H
2
⋅
j
ω
B
1
+
E
1
⋅
j
ω
D
2
+
E
1
J
2
)
d
V
=
−
∫
V
∇
⋅
(
E
1
×
H
2
)
d
V
.
{\displaystyle \int _{V}\left(\mathbf {H} _{2}\cdot j\omega \mathbf {B} _{1}+\mathbf {E} _{1}\cdot j\omega \mathbf {D} _{2}+\mathbf {E} _{1}\mathbf {J} _{2}\right)\mathrm {d} V=-\int _{V}\nabla \cdot (\mathbf {E} _{1}\times \mathbf {H} _{2})\mathrm {d} V\ .}
発散定理から、 の体積積分は 境界上の
の表面積積分に等しくなります。
div
(
E
1
×
H
2
)
{\displaystyle \operatorname {div} (\mathbf {E} _{1}\times \mathbf {H} _{2})}
E
1
×
H
2
{\displaystyle \mathbf {E} _{1}\times \mathbf {H} _{2}}
∫
V
(
H
2
⋅
j
ω
B
1
+
E
1
⋅
j
ω
D
2
+
E
1
⋅
J
2
)
d
V
=
−
∮
S
(
E
1
×
H
2
)
⋅
d
S
^
.
{\displaystyle \int _{V}\left(\mathbf {H} _{2}\cdot j\omega \mathbf {B} _{1}+\mathbf {E} _{1}\cdot j\omega \mathbf {D} _{2}+\mathbf {E} _{1}\cdot \mathbf {J} _{2}\right)\mathrm {d} V=-\oint _{S}(\mathbf {E} _{1}\times \mathbf {H} _{2})\cdot {\widehat {\mathrm {d} S}}\ .}
この形式は一般的な媒体では有効ですが、線形、等方性、時間不変の物質の一般的なケースでは、 εは 時間に依存しないスカラーです。したがって、一般に物理的な大きさ と
D
=
ϵ
E
{\displaystyle \mathbf {D} =\epsilon \mathbf {E} }
B
=
μ
H
.
{\displaystyle \mathbf {B} =\mu \mathbf {H} \ .}
最後の方程式は
∫
V
(
H
2
⋅
j
ω
μ
H
1
+
E
1
⋅
j
ω
ϵ
E
2
+
E
1
⋅
J
2
)
d
V
=
−
∮
S
(
E
1
×
H
2
)
⋅
d
S
^
.
{\displaystyle \int _{V}\left(\mathbf {H} _{2}\cdot j\omega \mu \mathbf {H} _{1}+\mathbf {E} _{1}\cdot j\omega \epsilon \mathbf {E} _{2}+\mathbf {E} _{1}\cdot \mathbf {J} _{2}\right)\mathrm {d} V=-\oint _{S}(\mathbf {E} _{1}\times \mathbf {H} _{2})\cdot {\widehat {\mathrm {d} S}}\ .}
まったく同様の方法で 、ベクトルについては 次の式が得られます。
E
2
{\displaystyle \mathbf {E} _{2}}
H
1
{\displaystyle \mathbf {H} _{1}}
∫
V
(
H
1
⋅
j
ω
μ
H
2
+
E
2
⋅
j
ω
ϵ
E
1
+
E
2
⋅
J
1
)
d
V
=
−
∮
S
(
E
2
×
H
1
)
⋅
d
S
^
.
{\displaystyle \int _{V}\left(\mathbf {H} _{1}\cdot j\omega \mu \mathbf {H} _{2}+\mathbf {E} _{2}\cdot j\omega \epsilon \mathbf {E} _{1}+\mathbf {E} _{2}\cdot \mathbf {J} _{1}\right)\operatorname {d} V=-\oint _{S}(\mathbf {E} _{2}\times \mathbf {H} _{1})\cdot {\widehat {\mathrm {d} S}}\ .}
最後の2つの方程式を要素で減算すると
、微分形式
QEDで同等になります。
∫
V
[
J
1
⋅
E
2
−
E
1
⋅
J
2
]
d
V
=
∮
S
[
E
1
×
H
2
−
E
2
×
H
1
]
⋅
d
S
.
{\displaystyle \int _{V}\left[\mathbf {J} _{1}\cdot \mathbf {E} _{2}-\mathbf {E} _{1}\cdot \mathbf {J} _{2}\right]\operatorname {d} V=\oint _{S}\left[\mathbf {E} _{1}\times \mathbf {H} _{2}-\mathbf {E} _{2}\times \mathbf {H} _{1}\right]\cdot \mathbf {\mathrm {d} S} \ .}
J
1
⋅
E
2
−
E
1
⋅
J
2
=
∇
⋅
[
E
1
×
H
2
−
E
2
×
H
1
]
{\displaystyle \ \mathbf {J} _{1}\cdot \mathbf {E} _{2}-\mathbf {E} _{1}\cdot \mathbf {J} _{2}=\nabla \cdot \left[\mathbf {E} _{1}\times \mathbf {H} _{2}-\mathbf {E} _{2}\times \mathbf {H} _{1}\right]\ }
表面条件のキャンセル
全空間にわたる積分に対するローレンツの相互定理の右辺の表面積項の打ち消しは、完全に自明ではないが、いくつかの方法で導くことができる。表面積分の厳密な取り扱いは、相互作用する波場状態の 因果関係 を考慮に入れる:無限遠における表面積分の寄与は、2つの因果波場の時間畳み込み相互作用に対してのみゼロになる(時間相関相互作用は非ゼロの寄与につながる)。 [10]
もう 1 つの単純な議論は、局所的なソースの場合、フィールドは無限遠でゼロになるというものですが、この議論は損失のない媒体の場合には当てはまりません。吸収がない場合、放射フィールドは距離に反比例して減衰しますが、積分の表面積は距離の 2 乗で増加するため、積分では 2 つの率が互いにバランスします。
代わりに、媒体は十分に離れたところから均質かつ等方性であると仮定するのが一般的です (例: King、1963)。この場合、放射場は、 および とともに、放射状に外側 ( 方向) に伝播する平面波の形を漸近的に取ります 。 ここ で 、 Z は 周囲の媒体のスカラー インピーダンス です。次に、 単純な ベクトル恒等式 により となり、同様に、 2 つの項は互いに打ち消し合います。
O
^
r
{\displaystyle \operatorname {\hat {O}} {\mathbf {r} }}
O
^
r
⋅
E
=
0
{\displaystyle \operatorname {\hat {O}} {\mathbf {r} }\cdot \mathbf {E} =0}
H
=
r
^
×
E
/
Z
{\displaystyle \mathbf {H} ={\hat {\mathbf {r} }}\times \mathbf {E} /Z}
μ
/
ϵ
{\textstyle {\sqrt {\mu /\epsilon }}}
E
1
×
H
2
=
E
1
×
r
^
×
E
2
Z
,
{\displaystyle \ \mathbf {E} _{1}\times \mathbf {H} _{2}={\frac {\mathbf {E} _{1}\times {\hat {\mathbf {r} }}\times \mathbf {E} _{2}}{Z}}\ ,}
E
1
⋅
E
2
Z
r
^
.
{\displaystyle {\frac {\mathbf {E} _{1}\cdot \mathbf {E} _{2}}{Z}}\ {\hat {\mathbf {r} }}\ .}
E
2
×
H
1
=
E
2
⋅
E
1
Z
r
^
{\displaystyle \mathbf {E} _{2}\times \mathbf {H} _{1}={\frac {\mathbf {E} _{2}\cdot \mathbf {E} _{1}}{Z}}\ {\hat {\mathbf {r} }}}
上記の議論は、表面項が打ち消される理由を明示的に示していますが、一般性に欠けています。あるいは、 極限吸収原理(LAP) によって課せられた放射境界条件で、損失のない周囲媒体の場合を扱うことができます。損失 ( ε の虚数部) がゼロになる極限を取ります 。ゼロ以外の損失については、媒体が均質であるかどうかに関係なく、フィールドは距離とともに指数関数的に減衰し、表面積分はゼロになります。ローレンツの相反定理の左辺は、ゼロ以外の損失を持つすべての空間での積分に対してゼロになるため、損失がゼロになる極限でもゼロにならなければなりません。(LAP は、無限遠からの入射波がゼロであるという ゾンマーフェルトの放射条件 を暗黙的に課していることに注意してください。そうでないと、任意の小さな損失でも入射波が除去され、極限では損失のない解が得られないためです。)
相互性とグリーン関数
演算子 の逆 、すなわち (これは、損失のないシステムでの無限遠での境界条件の指定を必要とする)は、 と同じ対称性を持ち 、本質的に グリーン関数の 畳み込みである。したがって、ローレンツの相互性に関する別の観点は、電磁グリーン関数との畳み込みが ε および μ の適切な条件下で複素対称(または以下では反エルミート)線形演算であるという事実を反映しているというものである 。より具体的には、グリーン関数は、における m 番目の方向 の点双極子電流から におけるの n 番目の成分 を与えると書くことができ (本質的に、 の行列要素を与える )、レイリー・カーソンの相互性は、 というステートメントと同等である。 とは異なり 、グリーン関数の明示的な式を与えることは一般には不可能であるが(均質媒体などの特殊な場合を除く)、数値的方法によって日常的に計算される。
O
^
,
{\displaystyle \operatorname {\hat {O}} \ ,}
E
=
O
^
−
1
J
{\displaystyle \mathbf {E} =\operatorname {\hat {O}} ^{-1}\mathbf {J} }
O
^
{\displaystyle \operatorname {\hat {O}} }
G
n
m
(
x
′
,
x
)
{\displaystyle G_{nm}(\mathbf {x} ',\mathbf {x} )}
E
{\displaystyle \mathbf {E} }
x
′
{\displaystyle \mathbf {x} '}
x
{\displaystyle \mathbf {x} }
G
{\displaystyle G}
O
^
−
1
{\displaystyle \operatorname {\hat {O}} ^{-1}}
G
n
m
(
x
′
,
x
)
=
G
m
n
(
x
,
x
′
)
.
{\displaystyle G_{nm}(\mathbf {x} ',\mathbf {x} )=G_{mn}(\mathbf {x} ,\mathbf {x} ')\ .}
O
^
,
{\displaystyle \operatorname {\hat {O}} \ ,}
ロスレス磁気光学材料
ε が 対称行列で は ない 場合の 1 つは、 磁気光学 材料の場合です。この場合は、通常のローレンツの相反性は成り立ちません (ただし、一般化については以下を参照)。磁気光学材料を許容し、材料の 吸収が無視できる 状況に限定すると、 ε と μ は一般に 3×3 複素 エルミート行列 です。この場合、演算子は 共役 内積 の下でエルミートであり 、相反定理の変形 [ 要出典 ] が依然として成立します。
ここで、符号が変化する理由は、 上の式の によるもので、これにより演算子は 反エルミート になります(表面項を無視した場合)。の特殊なケースでは、 エネルギー保存 則または ポインティング の定理が再述されます (ここでは上記とは異なり、損失のない材料を想定しているため)。電流による仕事の時間平均率 (の実数部で与えられる ) は、電力の時間平均外向きフラックス ( ポインティング ベクトル の積分) に等しくなります。ただし、同様に、この相互性バリアントについて空間全体で積分すると、表面項は一般には消えないため、追加の仮定なしではレイリー カーソン形式は成立しません。
1
μ
(
∇
×
∇
×
)
−
ω
2
c
2
ε
{\textstyle \ {\frac {1}{\mu }}\left(\nabla \times \nabla \times \right)-{\frac {\omega ^{2}}{c^{2}}}\varepsilon }
(
F
,
G
)
=
∫
F
∗
⋅
G
d
V
,
{\textstyle (\mathbf {F} ,\mathbf {G} )=\int \mathbf {F} ^{*}\cdot \mathbf {G} \,\mathrm {d} V\ ,}
−
∫
V
[
J
1
∗
⋅
E
2
+
E
1
∗
⋅
J
2
]
d
V
=
∮
S
[
E
1
∗
×
H
2
+
E
2
×
H
1
∗
]
⋅
d
A
{\displaystyle -\int _{V}\left[\mathbf {J} _{1}^{*}\cdot \mathbf {E} _{2}+\mathbf {E} _{1}^{*}\cdot \mathbf {J} _{2}\right]\mathrm {d} V=\oint _{S}\left[\mathbf {E} _{1}^{*}\times \mathbf {H} _{2}+\mathbf {E} _{2}\times \mathbf {H} _{1}^{*}\right]\cdot \mathbf {\mathrm {d} A} }
1
i
ω
{\displaystyle {\frac {1}{i\omega }}}
O
^
{\displaystyle \operatorname {\hat {O}} }
J
1
=
J
2
,
{\displaystyle \mathbf {J} _{1}=\mathbf {J} _{2}\ ,}
−
J
∗
⋅
E
{\displaystyle -\mathbf {J} ^{*}\cdot \mathbf {E} }
磁気光学材料がレイリー・カーソンの相反性を破壊するという事実は、 ファラデーアイソレータ や サーキュレータ などのデバイスにとって重要です。ファラデーアイソレータの片側の電流は反対側に磁場を生成しますが、 その逆は
起こりません。
非対称材料への一般化
損失性材料と磁気光学材料の組み合わせの場合、一般に ε テンソルと μ テンソルが対称行列でもエルミート行列でもない場合でも、 異なるシステム に と が存在する と見なすことで、ローレンツの相互関係の一般化バージョンを取得できます 。
(
J
1
,
E
1
)
{\displaystyle (\mathbf {J} _{1},\mathbf {E} _{1})}
(
J
2
,
E
2
)
{\displaystyle (\mathbf {J} _{2},\mathbf {E} _{2})}
特に、 が ω におけるマクスウェル方程式を物質を持つシステムに対して満たし 、 が転置を 表す物質を 持つ システムに対して ω におけるマクスウェル方程式を満たす場合、ローレンツの相互法則の式が成り立ちます。これは、 6×6の完全な磁化率テンソルを転置することで、 双異方性物質 にさらに一般化できます。 [11]
(
J
1
,
E
1
)
{\displaystyle (\mathbf {J} _{1},\mathbf {E} _{1})}
(
ε
1
,
μ
1
)
,
{\displaystyle (\varepsilon _{1},\mu _{1})\ ,}
(
J
2
,
E
2
)
{\displaystyle (\mathbf {J} _{2},\mathbf {E} _{2})}
(
ε
1
T
,
μ
1
T
)
,
{\displaystyle \left(\varepsilon _{1}^{\mathsf {T}},\mu _{1}^{\mathsf {T}}\right)\ ,}
T
{\displaystyle {}^{\mathsf {T}}}
相互主義の例外
非線形媒体 の場合 、一般的に相互性定理は成立しません。また、相互性は、時間によって変化する(「アクティブ」)媒体にも一般的には適用されません。たとえば、 ε が 何らかの外部プロセスによって時間的に変調される場合などです。(どちらの場合も、周波数 ω は一般的に保存量ではありません。)
フェルド・タイ相互主義
1992年に、密接に関連した相互性定理がYA Feld [12] とCT Tai [13]によって独立に表現され、 Feld-Tai相互性 または Feld-Tai補題
として知られています 。これは、2つの時間調和局所電流源とその結果生じる 磁場 を関連付けます。
∫
J
1
⋅
H
2
d
V
=
∫
H
1
⋅
J
2
d
V
.
{\displaystyle \int \mathbf {J} _{1}\cdot \mathbf {H} _{2}\,\operatorname {d} V=\int \mathbf {H} _{1}\cdot \mathbf {J} _{2}\,\operatorname {d} V\ .}
しかし、フェルド・タイの補題は、ローレンツの相反性よりもはるかに制限された条件下でのみ有効です。一般的に、この補題では、完全導電性材料の領域を除いて、等方性均質 インピーダンス 、つまり一定の スカラー μ / ε 比を持つ時間不変の線形媒体が必要です。
より正確には、フェルド-タイの相互性は、上記のように電磁演算子のエルミート対称性(またはむしろ複素対称性)を必要とするが、 および に関連する演算子が およびに 関連する演算子の定数スカラー倍であるという仮定にも依存しており 、 これは ε が μ の定数スカラー倍であるときに真である(2 つの演算子は一般に ε と μ の入れ替えによって異なる )。 上記のように、有限体積上の積分のより一般的な定式化を構築することもできます。
E
{\displaystyle \ \mathbf {E} \ }
i
ω
J
{\displaystyle \ i\omega \mathbf {J} \ }
H
{\displaystyle \ \mathbf {H} \ }
∇
×
(
J
/
ε
)
,
{\displaystyle \ \nabla \times (\mathbf {J} /\varepsilon )\ ,}
放射測定における光学的相反性
量子効果とは別に、古典理論は、任意の時間経過を伴う近距離場、中距離場、遠距離場の電気および磁気現象をカバーします。光学は、遠距離場のほぼ正弦波振動電磁気効果を指します。電気変数と磁気変数のペアの代わりに、光学的相互性を含む光学は、従来 比強度 と呼ばれていた スペクトル放射輝度 などの偏 光 ペア放射変数で表現できます。
1856年、 ヘルマン・フォン・ヘルムホルツは 次のように書いています。
「点A から発せられた光線は、 屈折、反射などを経て 点 Bに到達します。点 A において、光線の方向に2 つの垂直な平面 a 1 、 a 2 をとり、光線の振動をこれらの平面のそれぞれで 1 つずつ 2 つの部分に分割します。点 B の光線に同様の平面 b 1 、 b 2 を とります。すると、次の命題が証明されます。平面 a 1 で偏光された光量 J が A から 与えられた光線の方向に発せられるとき、その b 1 で偏光された光の部分 Kが B に到達します。逆に、 b 1 で偏光された 光量 J が B から発せられると 、同じ量の a 1 で偏光された光量 K が A に到達します 。」 [14]
これはヘルムホルツの相反性 (または反転)原理 と呼ばれることもあります。 [15] [16] [17] [18] [19] [20] 波が磁場の作用を受けた物質を伝わる場合、相反性は破れる可能性があり、この原理は適用されません。 [14] 同様に、光線の進路上に移動する物体がある場合、この原理はまったく適用できない場合があります。歴史的には、1849年に ジョージ・ストークス卿が 偏光を考慮せずに光学反転原理を述べました。 [21] [22] [23]
熱力学の原理と同様に、この原理は、実験が提案された法則のテストであるという通常の状況とは対照的に、実験が正しく実行されているかどうかをチェックするのに十分信頼できる。 [24] [25]
この原理を最も簡単に述べると、「 私があなたを見ることができれば、あなたも私を見ることができる」となります。この原理は 、グスタフ・キルヒホッフが 熱放射の法則 を導出する際に 使用され 、 マックス・プランクが 熱放射の法則 を分析する際にも使用されました 。
レイトレーシング グローバルイルミネーションアルゴリズムでは、入射光と出射光は 、双方向反射分布関数 (BRDF)の結果に影響を与えることなく、互いに反転したものとして考えることができます 。 [25]
グリーンの相互関係
上記の相互定理は振動場に対するものでしたが、 グリーンの相互定理は 電荷の 分布が固定された静電気に対する類似の定理です (Panofsky および Phillips、1962)。
特に、 全電荷密度 から生じる電位を とします 。電位は ポアソン方程式 を 満たし、 は 真空の誘電率 です 。同様に、 全電荷密度 から生じる電位を とします 。 を満たします 。どちらの場合も、電荷分布は局所的であると仮定し、電位は無限遠でゼロになるように選択できます。すると、グリーンの相反定理は、全空間の積分について次のようになります。
ϕ
1
{\displaystyle \phi _{1}}
ρ
1
{\displaystyle \rho _{1}}
−
∇
2
ϕ
1
=
ρ
1
/
ε
0
{\displaystyle -\nabla ^{2}\phi _{1}=\rho _{1}/\varepsilon _{0}}
ε
0
{\displaystyle \varepsilon _{0}}
ϕ
2
{\displaystyle \phi _{2}}
ρ
2
{\displaystyle \rho _{2}}
−
∇
2
ϕ
2
=
ρ
2
/
ε
0
{\displaystyle -\nabla ^{2}\phi _{2}=\rho _{2}/\varepsilon _{0}}
∫
ρ
1
ϕ
2
d
V
=
∫
ρ
2
ϕ
1
d
V
.
{\displaystyle \int \rho _{1}\phi _{2}dV=\int \rho _{2}\phi _{1}\operatorname {d} V~.}
この定理はグリーンの第二恒等式 から簡単に証明できる 。つまり、
∫
ϕ
2
(
∇
2
ϕ
1
)
d
V
=
∫
ϕ
1
(
∇
2
ϕ
2
)
d
V
,
{\displaystyle \int \phi _{2}(\nabla ^{2}\phi _{1})\operatorname {d} V=\int \phi _{1}(\nabla ^{2}\phi _{2})\operatorname {d} V\ ,}
つまり、それは エルミート演算子です(次のように部分を 2 回積分するとわかります)。
∇
2
{\displaystyle \nabla ^{2}}
参照
参考文献
^ Stumpf, M. (2018). アンテナ 理論における電磁的相互関係 。ピスカタウェイ、ニュージャージー:Wiley-IEEE プレス。Bibcode :2017erat.book.....S.
^ Balanis, CA (2024). Advanced Engineering Electromagnetics (第3版). Wiley . pp. 335–337. doi :10.1002/9781394180042. ISBN 978-1-394-18001-1 . 2024年 2月29日 閲覧 。
^ Harrington, RF (2001). 時間調和電磁場. Wiley - IEEE . pp. 116–120. doi :10.1109/9780470546710. ISBN 978-0-471-20806-8 . 2024年 2月29日 閲覧 。
^ Van Bladel、J. (2007)。電磁場 (第 2 版)。 ワイリー - IEEE 。 397–402ページ。 土井 :10.1002/047012458X。 ISBN 978-0-471-26388-3 . 2024年 2月29日 閲覧 。
^ Felsen, LB ; Marcuvitz, N. (2003). 波の放射と散乱. Wiley - IEEE . pp. 90–93. doi :10.1109/9780470546307. ISBN 978-0-780-31088-9 . 2024年 2月29日 閲覧 。
^ Kong, JA (2008). 電磁波理論 (第3版). EMW. pp. 697–702. ISBN 978-0-9668143-9-2 。
^ Collin, RE (1991). 導波の場の理論 (第 2 版). Wiley - IEEE . pp. 49–50. doi :10.1109/9780470544648. ISBN 978-0-879-42237-0 . 2024年 2月29日 閲覧 。
^ Chew, Wen Cho (2008 年 4 月)。「電磁気学における相互性とエネルギー保存定理の新た な 考察」 IEEE Transactions on Antennas and Propagation 56 ( 4): 970–975。Bibcode : 2008ITAP ...56..970C。doi :10.1109/TAP.2008.919189。S2CID 13615400 。
^ ラモ、サイモン、ウィナリー、ジョン、ヴァン・デューザー、セオドア (1965)。 通信エレクトロニクスにおける場と波 (国際版)。ジョン・ワイリー・アンド・サンズ 。ISBN
978-047170720-2 。
0471707201 出版社
^ Stumpf, M. (2018). アンテナ 理論における電磁的相反性 。IEEE Press / Wiley。§1.4.3。Bibcode :2017erat.book.....S.
^ Kong, Jin Au (1972 ) . 「 二 等方性媒体の定理」。IEEE 紀要 。60 (9): 1036–1046。doi :10.1109/PROC.1972.8851。
^ Feld, Ya.N. (1992). 「電気力学における二次方程式について」. ソビエト物理学-ドクラディ . 37 : 235–236.
^ Tai, C.-T. (1992). 「電磁理論における相補的相互性定理」. IEEE Transactions on Antennas and Propagation . 40 (6): 675–681. Bibcode :1992ITAP...40..675T. doi :10.1109/8.144602. hdl : 2027.42/21036 .
^ ab von Helmholtz, H. (1856)。 Handbuch der physiologischen Optik [生理光学 ハンドブック]。 Vol. 1(第1版)。ライプツィヒ:レオポルド・フォス。 p. 169; プランクが引用。ここで引用した英語版は「ヘルムホルツの翻訳」 に基づいています。 哲学雑誌 (第2刷)。シリーズ4。20 。 ガスリー、F.による翻訳:2–21 1867。
^ Minnaert, M. (1941). 「月の光測定における相互性原理」. 天体物理学ジャーナル . 93 : 403–410. Bibcode :1941ApJ....93..403M. doi :10.1086/144279.
^ Chandrasekhar, S. (1950). 放射伝達 . オックスフォード、イギリス: オックスフォード大学出版局. pp. 20–21, 171–177, 182.
^ ティングワルト、CP (1952)。 「Über das Helmholtzsche Reziprozitätsgesetz in der Optik」[光学におけるヘルムホルツの相反則について]。 オプティック 。 9 (6): 248–253。
^ Levi, L. (1968). 応用光学: 光学システム設計ガイド . 第1巻. ニューヨーク、NY: Wiley. p. 84. (全2巻)
^ Clarke, FJJ; Parry, DJ (1985). 「ヘルムホルツの相互性: その有効性と反射測定への応用」. 照明研究と技術 . 17 (1): 1–11. doi :10.1177/14771535850170010301. S2CID 123394330.
^ Born, M.; Wolf, E. (1999). 光学原理 :光の伝播、干渉、回折の電磁理論 (第7版)。ケンブリッジ大学出版局。p. 423。ISBN 0-521-64222-1 。
^ ストークス、GG (1849)。「ニュートンリングの中心点の完全な黒さについて、および反射光線と屈折光線の強度に関するフレネルの公式の検証について」。 ケンブリッジおよびダブリン数学ジャーナル 。新シリーズ。4 : 1–14。
^ Mahan, AI (1943). 「ストークスの可逆性原理の数学的証明」 アメリカ光学会誌 33 (11): 621–626. doi : 10.1364/JOSA.33.000621.
^ Lekner, J. (1987). 電磁波と粒子波の反射理論. ドルドレヒト: Martinus Nijhoff. pp. 33–37. ISBN 90-247-3418-5 – Google ブックス経由。
^ レイリー、JW ストラット、男爵 (1900)。「拡散反射における相互法則について」。 哲学雑誌 。シリーズ 5。49 :324-325。 {{cite journal}}: CS1 maint: multiple names: authors list (link)
^ ab Hapke, B. (1993). 反射および放射分光法の理論 。ケンブリッジ、イギリス: ケンブリッジ大学出版局。セクション 10C、263-264 ページ 。ISBN 0-521-30789-9 。
出典
Landau, LD ; Lifshitz, EM (1960). 連続媒体の電気力学 . Reading, MA: Addison-Wesley. §89.
King, RWP (1963). 基礎電磁気理論 (再版). ニューヨーク: ドーバー. §IV.21.
Altman, C.; Such, K. (1991). 電磁気学における相互性、空間マッピングおよび時間反転 。ドルドレヒト、オランダ: Kluwer。
ハワイ州ローレンツ (1895–1896)。 「電磁場のエネルギーに関するポインティングの定理と光の伝播に関する 2 つの一般命題」 (PDF) 。 Verhandelingen der Koninklijke Akademie van Wetenschappen 、アムステルダム[アムステルダム王立科学アカデミーの論文 ] 。 4 : 176–187 – ライデン大学 (leidenuniv.nl) 経由。
Potton, RJ (2004). 「光学における相互関係」. 物理学の進歩に関する報告書 . 67 (5): 717–754. Bibcode :2004RPPh...67..717P. doi :10.1088/0034-4885/67/5/R03. — このトピックの歴史に関するレビュー記事。
カーソン、JR (1924)。「逆数定理の一般化」。 ベルシステム技術ジャーナル 。3 (3): 393–399。doi :10.1002/ j.1538-7305.1924.tb00009.x – インターネットアーカイブ ( archive.org)経由 。
カーソン 、JR (1930)。「逆エネルギー定理」。 ベルシステム技術ジャーナル 。9 (4): 325–331。doi : 10.1002/j.1538-7305.1930.tb00369.x – インターネットアーカイブ (archive.org) 経由。
Feld, Ya. N. (1992). 「電気力学における二次補題について」. ソビエト物理学-ドクラディ . 37 : 235–236.
Tai, C.-T. (1992). 「電磁理論における相補的相互性定理」. IEEE Transactions on Antennas and Propagation . 40 (6): 675–681. Bibcode :1992ITAP...40..675T. doi :10.1109/8.144602. hdl : 2027.42/21036 .
Panofsky, Wolfgang KH; Phillips, Melba (1962). 「古典電気と磁気」 マサチューセッツ州レディング: Addison-Wesley.
Stumpf, M. (2018). アンテナ理論における電磁相互作用 。ピスカタウェイ、ニュージャージー:Wiley-IEEE プレス。Bibcode : 2017erat.book.....S.
Stumpf, M. (2020). アンテナモデリングにおける時間領域電磁相互作用 。ピスカタウェイ、ニュージャージー州: Wiley-IEEE プレス。
アサチー、VS;ミシシッピ州ミルムーサ。ディアス・ルビオ、A.ファン、S。南京州トレチャコフ(2020年10月)。 「電磁非相反性とその起源に関するチュートリアル」。 IEEE の議事録 。 108 (10): 1684 ~ 1727 年。 arXiv : 2001.04848 。 土井 :10.1109/JPROC.2020.3012381。