トムソン問題

トムソン問題（トムソンもんだい、英: Thomson problem） は、静電エネルギーが最小になるよう単位球面上に N 個の電子を配置する問題である（電子間力はクーロンの法則に従うものとする）。1904年、物理学者J. J. トムソンが、後年ブドウパンモデルと呼ばれる原子模型の発表後、「（電気的な）中性原子中には、負に帯電した電子が存在している」という知見を基に提出した。 この問題に関連して、最小エネルギー配置の幾何の研究や、N を大きくしたときの最小エネルギーの研究がある。 トムソン問題は、数学者スティーヴン・スメイルによる18個の未解決問題の1つ — "2次元球面上の点の分布" — の特別な場合である。冒頭述べたように単位球面（ r = 1 {\displaystyle r=1} ）の場合だけ考えればよい。 電気量 e i = e j = e {\displaystyle e_{i}=e_{j}=e} （ e {\displaystyle e} は電気素量）の電子ペアがもつ静電エネルギーは、クーロンの法則より U i j ( N ) = k e e i e j r i j {\displaystyle U_{ij}(N)=k_{e}{e_{i}e_{j} \over r_{ij}}} ここで k e {\displaystyle k_{e}} はクーロン定数、 r i j = | r i − r j | {\displaystyle r_{ij}=|\mathbf {r} _{i}-\mathbf {r} _{j}|} は位置ベクトル r i {\displaystyle \mathbf {r} _{i}} 、 r j {\displaystyle \mathbf {r} _{j}} で表した電子間距離。 e = 1 {\displaystyle e=1} 、 k e = 1 {\displaystyle k_{e}=1} であるとしても一般性は失わない。このとき、 U i j ( N ) = 1 r i j {\displaystyle U_{ij}(N)={1 \over r_{ij}}} N-電子配置のエネルギーの総和は、ペアごとの値を合算して次のように書ける。 U ( N ) = ∑ i < j 1 r i j {\displaystyle U(N)=\sum _{i<j}{\frac {1}{r_{ij}}}} 大局的な U ( N ) {\displaystyle U(N)} の最小化配置は、大体において数値計算で得られる。 電子2個の場合のトムソン問題の解は、電子が原点を挟んで最も遠ざかったときが解となる： r i j = 2 r = 2 , U ( 2 ) = 1 2 {\displaystyle r_{ij}=2r=2,\quad U(2)={1 \over 2}} 最小エネルギー配置が厳密に分かっているのはほんの限られた場合だけである。 N =1のとき：解は自明であり、電子は球面上のどこにあってもよい。このときのエネルギーは、電子が別の電場源からの影響を一切受けないことからゼロと定義することにする。 N =2のとき：2電子が対蹠点に位置するときが解。 N =3のとき：3電子が大円上の正三角形の頂点になるときが解。 N =4のとき：4電子が正四面体の頂点になるときが解。 N =5のとき：コンピュータを用いた数学的厳密解が2010年に報告された。5電子が双三角錐の頂点になるときが解。 N =6のとき：6電子が正八面体の頂点になるときが解。 N =12のとき：12電子が正二十面体の頂点になるときが解。 注目すべきことに、 N =4,6,12のときのトムソン問題の解は、全ての面が合同な正三角形であるプラトンの立体を形作っている。 N =8,20のときの数値計算による解は残り2つのプラトンの立体（面が正方形、正五角形）にはなっていない。 粒子間のポテンシャルが任意の形で与えられているときの基底状態を考えることもできる。数学的に正式に述べると、 f を実数値減少関数とし、全エネルギーを ∑ i < j f ( | r i − r j | ) {\displaystyle \sum _{i<j}f(|\mathbf {r} _{i}-\mathbf {r} _{j}|)} とする。 伝統的に f ( x ) = x − α {\displaystyle f(x)=x^{-\alpha }} （リースの α {\displaystyle \alpha } -カーネル）を考えることが多い。可積分なリースカーネルについては文献を参照のこと。また可積分でないリースカーネルについてはPoppy-seed bagel theorem（英語版）が成り立っている。 有名なケースは α → ∞ としたとき：「球面上の2点間の距離の最大値が最小となる配置を求める問題」と解釈する。Tammes problem（英語版）として知られている。 α = 1のとき：トムソン問題。 α = 0のとき： f ( x ) = ln ⁡ ( 1 | x | ) {\displaystyle f(x)=\ln \left({\frac {1}{|x|}}\right)} と解釈する。Whyte's problemとして知られている。 高次元空間での N 個の点の配置を考えることもできる。球デザイン（英語版）を参照。 トムソン問題は、トムソンのブドウパンモデルにおいて一様に広がる正電荷がないものとすれば自然に導かれる。 『自然哲学の大部分は原子の構造及び機構が生み出しているのだから、原子に関する発見に取るに足らないものはないし、物理学の進歩を加速させないものもない。』 —J. J.トムソン卿 実験によってトムソンのブドウパンモデルは原子の完全なモデルにはなり得ないことがわかったが、数値計算によるトムソン問題の変則的な解の中には、周期表上の一連の元素における自然な殻内電子配置と対応が見つかっているものもある。 トムソン問題はまた、電子バブル（英語版）や、四重極イオントラップ内の金属液滴の表面配向といった他の物理モデルの研究でも有用である。 一般化されたトムソン問題は、例えばウイルスの球状殻を構成するタンパク質のサブユニット配置を決定するときに現れる。この場合の「粒子」とは殻上のサブユニットのことである。他の具体例として、コロイドソーム（colloidosome；薬物、栄養物、生細胞等、活性成分のカプセル封じとして提案されている）におけるコロイド粒子の均整な配置、炭素原子のフラーレン構造、原子価殻電子対反発則の理論等が挙げられる。低温状態で、巨大単極子を中心としたシェル構造超伝導金属が形成するとされるアブリコソフ-ボルテックスは、対数ポテンシャルによる長距離相互作用の例である。 下記の表で、 N {\displaystyle N} は点（電荷）の個数、 E 1 {\displaystyle E_{1}} はエネルギー、対称性の記述はシェーンフリース記号（三次元の点群を参照）によるもの、また r i {\displaystyle r_{i}} は点電荷の位置を表す。ほとんどの対称性のタイプにおいてベクトルの総和はゼロになり、従って電気双極子モーメントもゼロになる。 点集合の凸包による多面体も併せて考えることが慣習となっている。ここで、 v i {\displaystyle v_{i}} は添字に等しい辺が集まる頂点の数、 e {\displaystyle e} は辺の数、 f 3 {\displaystyle f_{3}} は三角形の面の数、 f 4 {\displaystyle f_{4}} は四角形の面の数、 θ 1 {\displaystyle \theta _{1}} は隣接する2辺が成す角の最小値である。辺の長さは一定とは限らないことに注意。このため（ N =4, 6, 12, 24を除いて） 凸包は表の最終列に記した一様多面体またはジョンソンの立体とはグラフ同型であるにとどまる。 ^ Thomson, Joseph John (March 1904). “On the Structure of the Atom: an Investigation of the Stability and Periods of Oscillation of a number of Corpuscles arranged at equal intervals around the Circumference of a Circle; with Application of the Results to the Theory of Atomic Structure”. Philosophical Magazine. Series 6 7 (39): 237–265. doi:10.1080/14786440409463107. オリジナルの13 December 2013時点におけるアーカイブ。. https://web.archive.org/web/20131213172104/http://www.cond-mat.physik.uni-mainz.de/~oettel/ws10/thomson_PhilMag_7_237_1904.pdf.  ^ Smale, S. (1998). “Mathematical Problems for the Next Century”. Mathematical Intelligencer 20 (2): 7–15. doi:10.1007/bf03025291.  ^ Föppl, L. (1912). “Stabile Anordnungen von Elektronen im Atom”. J. Reine Angew. Math. (141): 251–301. http://eudml.org/doc/149380.  ^ Schwartz, Richard (2010). "The 5 electron case of Thomson's Problem". arXiv:1001.3702。 ^ Yudin, V.A. (1992). “The minimum of potential energy of a system of point charges”. Discretnaya Matematika 4 (2): 115–121 (in Russian). ; Yudin, V. A. (1993). “The minimum of potential energy of a system of point charges”. Discrete Math. Appl. 3 (1): 75–81. doi:10.1515/dma.1993.3.1.75.  ^ Andreev, N.N. (1996). “An extremal property of the icosahedron”. East J. Approximation 2 (4): 459–462.  MR1426716, Zbl 0877.51021 ^ Landkof, N. S. Foundations of modern potential theory. Translated from the Russian by A. P. Doohovskoy. Die Grundlehren der mathematischen Wissenschaften, Band 180. Springer-Verlag, New York-Heidelberg, 1972. x+424 pp. ^ Hardin, D. P.; Saff, E. B. Discretizing manifolds via minimum energy points. Notices Amer. Math. Soc. 51 (2004), no. 10, 1186–1194 ^ Y. Levin and J. J. Arenzon, ``Why charges go to the Surface: A generalized Thomson Problem Europhys. Lett. Vol. 63 p. 415 (2003) ^ Sir J.J. Thomson, The Romanes Lecture, 1914 (The Atomic Theory) ^ LaFave Jr, Tim (December 2013). “Correspondences between the classical electrostatic Thomson problem and atomic electronic structure”. Journal of Electrostatics 71 (6): 1029–1035. arXiv:1403.2591. doi:10.1016/j.elstat.2013.10.001. オリジナルの22 February 2014時点におけるアーカイブ。. https://web.archive.org/web/20140222100807/http://www.pagesofmind.com/FullTextPubs/La13-LaFave-2013-Correspondences-between-the-Thomson-Problem-and-Atomic-Structure.pdf 10 Feb 2014閲覧。.  ^ Kevin Brown. "Min-Energy Configurations of Electrons On A Sphere". Retrieved 2014-05-01. Whyte, L.L. (1952). “Unique arrangements of points on a sphere”. Amer. Math. Monthly 59 (9): 606–611. doi:10.2307/2306764.  Cohn, Harvey (1956). “Stability configurations of electrons on a sphere”. Math. Comput.: 117–120. doi:10.1090/S0025-5718-1956-0081133-0.  Goldberg, Michael (1969). “Stability configurations of electrons on a sphere”. Math. Comp. 23: 785–786. doi:10.1090/S0025-5718-69-99642-2.  Erber, T.; Hockney, G. M. (1991). “equilibrium configurations of N equal charges on a sphere”. J. Phys. A: Math. Gen. 24 (23): L1369. Bibcode: 1991JPhA...24L1369E. doi:10.1088/0305-4470/24/23/008.  Morris, J. R.; Deaven, D. M.; Ho, K. M. (1996). “Genetic-algorithm energy minimization for point charges on a sphere”. Phys. Rev. B 53: R1740. Bibcode: 1996PhRvB..53.1740M. doi:10.1103/PhysRevB.53.R1740.  Erber, T.; Hockney, G. M. (1997). “Complex Systems: Equilibrium Configurations of N {\displaystyle N} Equal Charges on a Sphere ( 2 ≤ N ≤ 112 ) {\displaystyle (2\leq N\leq 112)} ”. Advances in Chemical Physics 98: 495–594. doi:10.1002/9780470141571.ch5. . Altschuler, E. L.; Williams, T. J.; Ratner, E. R.; Tipton, R.; Stong, R.; Dowla, F.; Wooten, F. (1997). “Possible global minimum lattice configurations for Thomson's problem of charges on a sphere”. Phys. Rev. Lett. 78: 2681. Bibcode: 1997PhRvL..78.2681A. doi:10.1103/PhysRevLett.78.2681.  Bowick, M.; Cacciuto, A.; Nelson, D. R.; Travesset, A. (2002). “Crystalline order on a sphere and the generalized Thomson Problem”. Phys. Rev. Lett. 89 (18): 249902. arXiv:cond-mat/0206144. Bibcode: 2002PhRvL..89r5502B. doi:10.1103/PhysRevLett.89.185502.  Dragnev, P. D.; Legg, D. A.; Townsend, D. W. (2002). “Discrete logarithmic energy on the sphere”. Pacific J. Math. 207 (2): 345–358. doi:10.2140/pjm.2002.207.345. . Katanforoush, A.; Shahshahani, M. (2003). “Distributing points on the sphere. I”. Exper. Math. 12 (2): 199–209. doi:10.1080/10586458.2003.10504492.  Wales, David J.; Ulker, Sidika (2006). “Structure and dynamics of spherical crystals characterized for the Thomson problem”. Phys. Rev. B 74 (21): 212101. Bibcode: 2006PhRvB..74u2101W. doi:10.1103/PhysRevB.74.212101.  Configurations reprinted in “The Cambridge cluster database”. 2018年6月14日閲覧。 Slosar, A.; Podgornik, R. (2006). “On the connected-charges Thomson problem”. Europhys. Lett. 75 (4): 631. arXiv:cond-mat/0606765. Bibcode: 2006EL.....75..631S. doi:10.1209/epl/i2006-10146-1.  Cohn, Henry; Kumar, Abhinav (2007). “Universally optimal distribution of points on spheres”. J. Amer. Math. Soc. 20 (1): 99–148. arXiv:math/0607446. Bibcode: 2007JAMS...20...99C. doi:10.1090/S0894-0347-06-00546-7.  Wales, D. J.; McKay, H.; Altschuler, E. L. (2009). “Defect motifs for spherical topologies”. Phys. Rev. B 79 (22): 224115. Bibcode: 2009PhRvB..79v4115W. doi:10.1103/PhysRevB.79.224115. . Configurations reproduced in “The Cambridge cluster database”. 2018年6月14日閲覧。 Ridgway, W. J. M.; Cheviakov, A. F. (2018). “An iterative procedure for finding locally and globally optimal arrangements of particles on the unit sphere”. Comp. Phys. Commun. 233: 84-109. doi:10.1016/j.cpc.2018.03.029.  “Thomson Problem @ S.U.”. 2009年11月24日閲覧。