We give an integral expression for the vector potential of a time-independent, steady azimuthal current density. Our derivation is substantially simpler and somewhat more general than others given in the literature. As an illustration, we recover the results for the vector potential of a circular current loop as an orthogonal expansion in spherical and cylindrical coordinates. Additionally, we obtain closed analytical expressions for the vector potential and the magnetic induction of a circular current loop in terms of Legendre functions of the second kind, that are simpler than the results in terms of complete elliptic integrals given in textbooks.
論文ID : 2510.11663タイトル : The vector potential of a steady azimuthal current density. Once again著者 : Antonio O. Bouzas (CINVESTAV-IPN, Mexico)分類 : physics.class-ph (古典物理学)発表日 : 2025年10月13日論文リンク : https://arxiv.org/abs/2510.11663 本論文は、時間に無関係な定常方位角電流密度のベクトルポテンシャルの積分表現を提示する。既存文献の他の導出と比較して、本論文の導出方法はより簡潔で、より普遍的である。応用例として、本論文は円形電流ループのベクトルポテンシャルの球座標系および円筒座標系における直交展開結果を復元する。さらに、本論文は円形電流ループのベクトルポテンシャルと磁束密度の閉形式解析表現を得た。この表現は第2種ルジャンドル関数に基づいており、教科書の完全楕円積分に基づく結果よりも簡潔である。
本研究が解決する中核的な問題は、方位角対称性を有する定常電流密度分布に対して、より簡潔にそのベクトルポテンシャルを求解する方法である。具体的には、電流密度が以下の形式を有する場合である:
J ( r ) = J ϕ ( r ) ϕ ^ \mathbf{J}(\mathbf{r}) = J_\phi(\mathbf{r})\hat{\phi} J ( r ) = J ϕ ( r ) ϕ ^
ここで ϕ ^ = − sin ( ϕ ) x ^ + cos ( ϕ ) y ^ \hat{\phi} = -\sin(\phi)\hat{\mathbf{x}} + \cos(\phi)\hat{\mathbf{y}} ϕ ^ = − sin ( ϕ ) x ^ + cos ( ϕ ) y ^ は方位角単位ベクトルである。
理論的意義 :方位角対称電流分布は電磁気学において重要な地位を占め、より複雑な電磁場分布を理解するための基礎である教育的価値 :大学院生および高学年の学部生に対して、より明確な磁気静力学の理解を提供する実用的価値 :円形電流ループは多くの実際の電磁デバイスの基本構成要素である複雑性 :既存文献の導出過程は比較的複雑であり、特に球座標系での処理が複雑である限定性 :ほとんどの導出は特定の座標系にのみ適用可能である表現の複雑性 :従来の結果は通常、完全楕円積分で表現され、形式が複雑である著者は、より簡潔で普遍的な導出方法を提供し、より簡単な解析表現を得ることを目指している。特に、円形電流ループというこの古典的問題に対してである。
普遍的な積分表現の提示 :方位角を含む任意の直交曲線座標系に適用可能なベクトルポテンシャルの積分公式を提示した導出過程の簡略化 :文献3 の球座標導出および文献2 の円筒座標導出と比較して、本論文の方法はより簡潔である新しい閉形式解析解の獲得 :第2種ルジャンドル関数を用いた円形電流ループのベクトルポテンシャルの表現を得た。これは楕円積分形式よりも簡潔である関数関係の確立 :第2種ルジャンドル関数と完全楕円積分の間の新しい関係式を導出した完全な磁場解の提供 :円形電流ループの磁束密度の閉形式解析表現を提示した定常条件を満たす方位角電流密度が与えられた場合:
∇ ⋅ J ( r ) = 0 , ∂ J ϕ ∂ ϕ ( r ) = 0 \nabla \cdot \mathbf{J}(\mathbf{r}) = 0, \quad \frac{\partial J_\phi}{\partial \phi}(\mathbf{r}) = 0 ∇ ⋅ J ( r ) = 0 , ∂ ϕ ∂ J ϕ ( r ) = 0
クーロンゲージ下でのベクトルポテンシャルを求解する:
∇ ⋅ A ( r ) = 0 , ∇ 2 A ( r ) = − μ 0 J ( r ) \nabla \cdot \mathbf{A}(\mathbf{r}) = 0, \quad \nabla^2\mathbf{A}(\mathbf{r}) = -\mu_0\mathbf{J}(\mathbf{r}) ∇ ⋅ A ( r ) = 0 , ∇ 2 A ( r ) = − μ 0 J ( r )
ベクトルポアソン方程式の解から出発する:
A ( r ) = μ 0 4 π ∫ d 3 r ′ J ( r ′ ) ∣ r − r ′ ∣ \mathbf{A}(\mathbf{r}) = \frac{\mu_0}{4\pi}\int d^3\mathbf{r}' \frac{\mathbf{J}(\mathbf{r}')}{|\mathbf{r} - \mathbf{r}'|} A ( r ) = 4 π μ 0 ∫ d 3 r ′ ∣ r − r ′ ∣ J ( r ′ )
電流密度を代入し、方位角積分を分離する:
A ( r ) = μ 0 4 π ∫ d V 2 ′ J ϕ ( r ′ ) ∫ − π π d ϕ ′ ϕ ^ ′ ∣ r − r ′ ∣ \mathbf{A}(\mathbf{r}) = \frac{\mu_0}{4\pi}\int dV'_2 J_\phi(\mathbf{r}') \int_{-\pi}^{\pi} d\phi' \frac{\hat{\boldsymbol{\phi}}'}{|\mathbf{r} - \mathbf{r}'|} A ( r ) = 4 π μ 0 ∫ d V 2 ′ J ϕ ( r ′ ) ∫ − π π d ϕ ′ ∣ r − r ′ ∣ ϕ ^ ′
ここで d V 2 ′ dV'_2 d V 2 ′ は方位角を除く体積要素を表す。
距離関数を導入する:
R ( α ) = ∣ r − r ′ ∣ ϕ ′ = ϕ + α R(\alpha) = |\mathbf{r} - \mathbf{r}'|_{\phi' = \phi + \alpha} R ( α ) = ∣ r − r ′ ∣ ϕ ′ = ϕ + α
変数置換 ϕ ′ ′ = ϕ ′ − ϕ \phi'' = \phi' - \phi ϕ ′′ = ϕ ′ − ϕ により、被積分関数の周期性と奇偶性を利用する:
A ( r ) = μ 0 4 π ϕ ^ ∫ d V 2 ′ J ϕ ( r ′ ) ∫ − π π d ϕ ′ cos ( ϕ ′ − ϕ ) R ( ϕ ′ − ϕ ) \mathbf{A}(\mathbf{r}) = \frac{\mu_0}{4\pi}\hat{\boldsymbol{\phi}}\int dV'_2 J_\phi(\mathbf{r}') \int_{-\pi}^{\pi} d\phi' \frac{\cos(\phi' - \phi)}{R(\phi' - \phi)} A ( r ) = 4 π μ 0 ϕ ^ ∫ d V 2 ′ J ϕ ( r ′ ) ∫ − π π d ϕ ′ R ( ϕ ′ − ϕ ) c o s ( ϕ ′ − ϕ )
単位ベクトル n ^ \hat{\mathbf{n}} n ^ を巧妙に定義することにより:
球座標系:n ^ = sin ( θ ) r ^ + cos ( θ ) θ ^ \hat{\mathbf{n}} = \sin(\theta)\hat{\mathbf{r}} + \cos(\theta)\hat{\boldsymbol{\theta}} n ^ = sin ( θ ) r ^ + cos ( θ ) θ ^ 円筒座標系:n ^ = ρ ^ \hat{\mathbf{n}} = \hat{\boldsymbol{\rho}} n ^ = ρ ^ 異なる座標系の統一的な処理を実現した。
被積分関数の奇偶性を識別することにより、奇関数項を自動的に消去し、積分過程を大幅に簡略化した。
三角関数と距離関数の 2 π 2\pi 2 π 周期性を十分に活用し、積分区間の処理を簡略化した。
本論文は円形電流ループを検証ケースとして選択した。理由は以下の通りである:
比較のための完全な解析解を有する 球座標系および円筒座標系の両方で標準的な結果を有する 電磁気学教育における古典的問題である 球座標系 :J ϕ = I a sin θ ′ δ ( cos θ ′ ) δ ( r ′ − a ) J_\phi = \frac{I}{a}\sin\theta' \delta(\cos\theta')\delta(r' - a) J ϕ = a I sin θ ′ δ ( cos θ ′ ) δ ( r ′ − a ) 円筒座標系 :J ϕ = I δ ( z ′ ) δ ( ρ ′ − a ) J_\phi = I\delta(z')\delta(\rho' - a) J ϕ = I δ ( z ′ ) δ ( ρ ′ − a ) ここで I I I はループ電流強度、a a a はループ半径である。
Mathematica ソフトウェアを使用して数値計算と特殊関数の評価を実施し、異なる表現の数学的等価性を確保した。
Jackson教科書と一致する球面調和関数展開を得た:
A ϕ = μ 0 a I 2 ∑ ℓ = 1 ∞ 1 ( ℓ + 1 ) ℓ r < ℓ r > ℓ + 1 P ℓ 1 ( 0 ) P ℓ 1 ( cos θ ) A_\phi = \frac{\mu_0aI}{2}\sum_{\ell=1}^{\infty} \frac{1}{(\ell+1)\ell} \frac{r_<^\ell}{r_>^{\ell+1}} P_\ell^1(0)P_\ell^1(\cos\theta) A ϕ = 2 μ 0 a I ∑ ℓ = 1 ∞ ( ℓ + 1 ) ℓ 1 r > ℓ + 1 r < ℓ P ℓ 1 ( 0 ) P ℓ 1 ( cos θ )
Zangwill教科書のベッセル関数積分形式を復元した:
A ϕ ( r ) = μ 0 I a 2 ∫ 0 ∞ d k J 1 ( k a ) J 1 ( k ρ ) e − k ∣ z ∣ A_\phi(\mathbf{r}) = \frac{\mu_0Ia}{2}\int_0^{\infty} dk J_1(ka)J_1(k\rho)e^{-k|z|} A ϕ ( r ) = 2 μ 0 I a ∫ 0 ∞ d k J 1 ( ka ) J 1 ( k ρ ) e − k ∣ z ∣
ベクトルポテンシャルのルジャンドル関数表現 :
A ϕ ( r ) = μ 0 I 2 π a ρ Q 1 / 2 ( ξ ) , ξ = ρ 2 + z 2 + a 2 2 a ρ A_\phi(\mathbf{r}) = \frac{\mu_0I}{2\pi}\sqrt{\frac{a}{\rho}}Q_{1/2}(\xi), \quad \xi = \frac{\rho^2 + z^2 + a^2}{2a\rho} A ϕ ( r ) = 2 π μ 0 I ρ a Q 1/2 ( ξ ) , ξ = 2 a ρ ρ 2 + z 2 + a 2
磁束密度の閉形式表現 :
B ρ = − μ 0 I 2 π 2 a / ρ z ( ρ − a ) 2 + z 2 ( ρ + a ) 2 + z 2 Q 1 / 2 1 ( ξ ) B_\rho = -\frac{\mu_0I}{2\pi}\frac{2\sqrt{a/\rho} z}{\sqrt{(\rho-a)^2 + z^2}\sqrt{(\rho+a)^2 + z^2}}Q_{1/2}^1(\xi) B ρ = − 2 π μ 0 I ( ρ − a ) 2 + z 2 ( ρ + a ) 2 + z 2 2 a / ρ z Q 1/2 1 ( ξ )
B z = μ 0 I 2 π 1 ρ a ρ ( 1 2 Q 1 / 2 ( ξ ) + ρ 2 − a 2 − z 2 ( ρ − a ) 2 + z 2 ( ρ + a ) 2 + z 2 Q 1 / 2 1 ( ξ ) ) B_z = \frac{\mu_0I}{2\pi}\frac{1}{\rho}\sqrt{\frac{a}{\rho}}\left(\frac{1}{2}Q_{1/2}(\xi) + \frac{\rho^2 - a^2 - z^2}{\sqrt{(\rho-a)^2 + z^2}\sqrt{(\rho+a)^2 + z^2}}Q_{1/2}^1(\xi)\right) B z = 2 π μ 0 I ρ 1 ρ a ( 2 1 Q 1/2 ( ξ ) + ( ρ − a ) 2 + z 2 ( ρ + a ) 2 + z 2 ρ 2 − a 2 − z 2 Q 1/2 1 ( ξ ) )
ρ / a \rho/a ρ / a z / a z/a z / a 積分形式(30) 楕円積分(31) ルジャンドル関数(33) 0.35 0.0 0.577159 0.577159 0.577159 3.50 0.0 0.132367 0.132367 0.132367 0.35 1.74 0.065912 0.065912 0.065912 3.50 1.74 0.091995 0.091995 0.091995
数値検証により、3つの表現が6桁の有効数字内で完全に一致することが示された。
第2種ルジャンドル関数と完全楕円積分の新しい関係を確立した:
Q 1 / 2 ( z ) = 2 z + 1 ( z K ( 2 z + 1 ) − ( z + 1 ) E ( 2 z + 1 ) ) Q_{1/2}(z) = \sqrt{\frac{2}{z+1}}\left(zK\left(\sqrt{\frac{2}{z+1}}\right) - (z+1)E\left(\sqrt{\frac{2}{z+1}}\right)\right) Q 1/2 ( z ) = z + 1 2 ( zK ( z + 1 2 ) − ( z + 1 ) E ( z + 1 2 ) )
Jackson『古典電動力学』 :円形電流ループの楕円積分解と球面調和関数展開を提供Zangwill『現代電動力学』 :円筒座標系におけるベッセル関数積分表現を提示Vasilyev (2012) :球座標系における軸対称電流のベクトルポテンシャルを導出導出の簡潔性 :Vasilyevの球座標導出と比較して、本論文の方法はより直接的である普遍性 :Zangwillの円筒座標法と比較して、本論文は任意の直交曲線座標系に適用可能である表現の簡潔性 :ルジャンドル関数表現は楕円積分形式よりも簡潔である完全性 :ベクトルポテンシャルと磁場の両方の閉形式表現を同時に提供する普遍的公式の確立 :方位角電流密度のベクトルポテンシャルが統一された積分表現を有することを証明した古典的問題の簡略化 :円形電流ループ問題に対してより簡潔な解析解を提供した数学的関連性の解明 :異なる特殊関数間の新しい関係を確立した適用範囲の制限 :方位角対称性を有する電流分布にのみ適用可能である空間局所性の仮定 :電流密度が空間内で局所的に分布することを要求する定常条件 :時間に無関係な静的な場合のみを考慮する時変の場合への推広 より複雑な電流分布の考慮 他の特殊関数表現の探索 方法の革新性が強い :既存文献よりも簡潔な導出方法を提供する数学的厳密性が高い :導出過程は論理的に明確で、数学的処理は規範的である実用的価値が顕著 :新しい解析表現は計算上より便利である教育的価値が大きい :電磁気学教育に対してより明確な理解の途を提供する検証が充分 :複数の方法による数値比較により結果の信頼性を確保する応用範囲が限定的 :特定の対称性を有する電流分布にのみ適用可能である革新の程度が相対的に限定的 :主に既知の結果の再導出と表現方法の改善である新しい物理的洞察の欠如 :物理概念の突破というより、数学技法の改善である学術的貢献 :電磁気学理論に新しい数学的ツールと表現方法を提供する教育的価値 :電磁気学コースの教育品質の改善に有利である実用的意義 :簡略された表現は工学計算と数値シミュレーションに有利である再現可能性 :導出過程が明確で、結果は検証と応用が容易である理論電磁気学研究 :より複雑な問題の基礎構成要素として電磁気学教育 :学生に対してより直感的な理解方法を提供する工学計算 :円形コイルの磁場を正確に計算する必要がある応用において数値シミュレーション :数値方法の精度を検証するための解析基準としてJ. D. Jackson, "Classical Electrodynamics," Third Edition (1999) A. Zangwill, "Modern Electrodynamics" (2012) A. Vasilyev, "Vector Potential and Magnetic Field of Axially Symmetric Currents," arXiv:1208.4983 A. Erdelyi et al., "Higher Transcendental Functions" (1953) A. Erdelyi et al., "Tables of Integral Transforms" (1954) 総合評価 :これは高品質の理論物理学論文である。古典的な問題を扱っているが、巧妙な数学技法により、より簡潔で優雅な解決方法を提供している。本論文の主要な価値は方法論の改善と新しい数学表現形式の発見にあり、電磁気学理論と教育の両方に積極的な貢献をしている。