In Silico

マテリアルズインフォマティクス・材料

対称性ゼロの予測可否、パリティ・ビット一つで決まる

2026/8/28

目次
※ 概念図(背景→問い)・作図:AI 【背景】 対称性の形だけで、量がぴったりゼロに固定される場合がある そこでの非ゼロ予測は、不正確ではなく、起こりえない値になる 【問い】 そのゼロを外せるかどうかは、訓練前の何が決めているのか 特徴量が鏡映の区別を持つかどうか——設計の一点で決まる
※ 概念図(背景→問い)・作図:AI。数値は含まない。対称性が量をゼロへ固定するかどうかと、そのゼロを模型が外せてしまうかどうかは、別の話である。
背景

結晶に力を加えると電圧が生じることがある。この現象を圧電効果といい、その大きさは圧電テンソルという量で表す。ただし、どの結晶でも生じるわけではない。たとえばルチル型の二酸化チタン(TiO₂)は圧電効果を示さない。測ると小さい値になる、という意味ではない。ルチル型TiO₂の圧電テンソルは厳密にゼロである。

ゼロになる理由は、測定の精度ではなく結晶の形から決まる。ルチル型TiO₂は中心対称の結晶である——構造の中に反転の中心があり、その点を通して座標を反転させると元の構造に重なる。一方、圧電テンソルは反転操作のもとで符号が反転する量である。したがって中心対称の結晶では、圧電テンソルは自分自身の符号を反転させたものと等しくなければならない。この条件を満たす値はゼロだけである。結晶が持つ対称操作のもとで物性テンソルの形が保たれなければならないというこの規則を、結晶学ではノイマンの原理と呼ぶ。

したがって、中心対称の結晶に非ゼロの圧電応答を予測する模型は、少し外れた答えを出しているのではない。物理的に起こりえない答えを出している。

機械学習で物性を予測する模型がこのゼロを守れるかどうかは、精度の指標には表れない。訓練データにゼロという答えが十分に含まれていれば模型もゼロに近い値を返すようになる、という期待は自然である。だが、ゼロを訓練データから学ぶことと、模型の設計上ゼロ以外を出せないこととは、別のことである。

問い

中心対称の結晶に対して、模型が禁じられた非ゼロの値を出せてしまうかどうかは、何が、いつ決めているのか。

要点

出せてしまうかどうかは、訓練でどう学ばせたかではなく、訓練を始める前の設計一つで決まる。決め手は、模型の特徴量がパリティのラベルを持つかどうかである。 パリティとは、構造を鏡に映す操作(鏡映)に対して量の符号がどう変わるかの区別をいう。回転と鏡映の両方に出力が対応して変わる模型をO(3)等変、回転にだけ対応して鏡映を区別しない模型をSO(3)等変という。

2026年8月19日公開のプレプリントが、この違いを測っている1。対象は、ノイマンの原理で圧電テンソルが厳密にゼロへ強制される中心対称の結晶2,000種である。O(3)等変の模型(NequIP・Allegro・MACEの3核)は、18,000回の結晶ごとの実行すべてで閾値0.01 C m⁻²を下回り、中央値は浮動小数点演算の床に張り付いた。SO(3)等変の模型は、結晶の90〜96%でこの閾値を超え、その中央値は訓練データに含まれる本物の圧電応答と同じ大きさだった。両者の差は6桁である。

しかもO(3)側は、圧電応答が本当にゼロでない結晶——非中心対称の結晶——での予測誤差(MAE=平均絶対誤差。小さいほど良い)も低い。厳密性は精度と引き換えに買ったものではない。

回転の等変性を外して速度を得る設計は、すでに主流である

分子や結晶を扱う機械学習ポテンシャルの多くは、回転についての等変性を模型の構造そのものに組み込んで作られてきた。この設計は精度で有利になりやすい一方、計算コストを押し上げる。だから、この等変性をあえて外し、データからの学習に任せるという設計判断が広がっている。

Orb-v3はその代表例である。著者らは、非等変・非保存力(力をエネルギーの微分として定義しない構成)の構成でも「物理的な性質——ポテンシャルエネルギー面の高階微分を要するものも含めて——を正確に模型化できる」と述べ、性能と速度のトレードオフの前線を広げたと主張する2。レイテンシを10倍を超えて削減し、メモリ使用量を8倍を超えて削減しながら、多くの評価でほぼ最先端の性能を維持したという2。等変性を外すことは、実務者がいま10倍規模の高速化を手にするための主流の手段であり、Orb-v3の主張は、著者ら自身が走らせたベンチマークにもとづく。なおOrb-v3は、力をエネルギーの微分として計算するconservative版と、しないdirect版の両方を公開しており、著者らはどちらを選ぶかはワークフロー次第だと書いている2

より一般の設定でも、対称性は学習で埋め合わせられるという報告がある。剛体どうしの相互作用を汎用の変換器で学習させた研究は、模型規模・学習ステップ数・データ量を振ったうえで、等変性はデータ効率を上げる一方、「十分なエポック数があれば、データ拡張を使った非等変模型の訓練でこの差を埋められる」と報告する(TMLR掲載)3。計算予算を同じだけ振った比較では、検証したどの予算でも等変模型が非等変模型を上回るが、模型規模と学習時間への予算配分の最適点そのものは、等変・非等変の両方で異なるとも書かれている3

Orb-v3の著者らは結論で、保存性と等変性が汎用ポテンシャルの厳格な前提だという枠組みを反証したと書いている2。Orb-v3とBrehmerらの比較が示すのは、等変性を外すという選択が理にかなう場面が現に存在するということだ。回転についての対称性の誤りは、精度の指標に現れる近似の誤差であり、速度と引き換えに許容できる場合がある。問題は、対称性の誤りを近似の誤差として扱うこの考え方が、どこまで通用するかである。

対称性が強制するゼロは、回転の話とは別の問題

2026年8月19日に公開されたプレプリントは、対称性の誤りを近似の誤差として扱う考え方が通用しない場合があると指摘する1。対象は圧電テンソルで、判定に使うのはノイマンの原理である。結晶が反転対称の中心を持つ——中心対称——場合、圧電応答は反転操作のもとで自分自身の符号を反転させたものと一致しなければならず、それを満たす値はゼロしかない。

ノイマンの原理でゼロかどうかを判定するには、回転についての情報だけでは足りない。反転操作は鏡映——鏡に映したときの変化——に属する操作で、回転をどれだけ組み合わせても再現できない。模型の特徴量が回転にしか対応しないSO(3)等変で作られている場合、鏡映によって符号が反転するという性質そのものを表現する手段を持たない。パリティのラベルまで特徴量に持たせるO(3)等変の設計だけが、この符号反転を表現できる。

原典は次のように主張する。回転についての対称性を模型にどこまで組み込むかという議論は、これまで近似の誤差の大小として扱われてきた。だが一部の制約は厳密で、対称性はある種の物性テンソルを厳密にゼロへ追い込み、そこでの非ゼロ予測は不正確なのではなく物理的に起こりえない1。模型がこの種の予測をできてしまうかどうかは、訓練の結果ではなく、訓練前に特徴量がパリティのラベルを持つかどうかという、一つの、めったに報告されない設計上のビットで決まる、と原典は書く1。原典はこの基準を「パリティ・ギャップ」と呼び、群論だけから、どの物性・どの結晶がこの制約に実際に晒されているかを計算する1

6桁の差は、精度と引き換えに買ったものではない

NequIP・Allegro・MACEという3つの核それぞれにO(3)腕とSO(3)腕を用意し、パリティの有無以外は同じ条件に揃える。加えてEqV2(SO(3)腕しか存在しない核)も含め、2,000種の中心対称結晶で圧電テンソルを予測させ、閾値0.01 C m⁻²と比べる。この閾値は訓練データにある本物の圧電テンソルの経験分布から較正した値で、著者らが置いたものではない1。以下の数字は、対称性を理想化した座標での結果である——生の座標での挙動は後述する。O(3)腕は18,000回の結晶ごとの実行のうち、一件もこの閾値を超えなかった(内訳は3核×2,000結晶×3訓練seed=18,000で、本稿の再構成)。中央値は浮動小数点演算の床——NequIPで3.19×10⁻⁷ C m⁻²、Allegroで3.84×10⁻⁷、MACEで2.76×10⁻⁶——に張り付き、最大の単一予測でも閾値の7.6分の1にとどまる。SO(3)腕は結晶の90〜96%でこの閾値を超え、原典はこれを「誤検出(false flags)」と呼ぶ(2,000件のうち166件は m3̄m で、パリティが無くても回転だけでゼロが強制される。それを除いた1,834件に限れば、3核のSO(3)腕は97.6〜99.2%で閾値を超える1)。EqV2にはO(3)腕自体が存在せず、この範囲の上端に位置する。誤検出した予測の中央値は、訓練データに含まれる本物の圧電応答の大きさと同じ桁——たとえばNequIPのSO(3)腕(拡張なし)では、2,000結晶全体での違反の中央値が6.66×10⁻¹ C m⁻²——にあり、数値誤差ではなく本物のサイズの値を出している。

厳密性は精度と引き換えではない。圧電応答が本当にゼロでない結晶——非中心対称の分割——での予測誤差は、NequIPでO(3)腕0.2083に対しSO(3)腕(拡張なし)0.2405、Allegroで0.2140対0.2589、MACEで0.2222対0.2567(単位はC m⁻²)——いずれもO(3)腕のほうが低い、つまり良い。原典は、ゼロを混ぜた拡張がどの核でも回帰誤差を改善したとも書き、それを「不可能な値を出し続ける、より良い回帰器」と呼んでいる1。原典は限界として、O(3)腕3種がいずれもe3nn実装であり、訓練後の精度比較に実装由来の寄与が混じる可能性を排除できないと自ら断っている1。ただし同じ段落は、非e3nn実装でも構造的なゼロ自体は再現したと書いている1。しかもSO(3)腕はO(3)腕よりパラメータ数が0.2〜39.8%多い。容量で有利なほうが、精度でも劣っている。測っている量は誤検出と同じ圧電テンソルで、違うのは分割のほうだ——中心対称の結晶ではゼロを外さないか、非中心対称の結晶では本物の値をどれだけ当てられるか。

この分離は、原典が閾値を10⁻⁴から10⁰ C m⁻²まで4桁振っても崩れない。読み出しを平均プーリングに替え、閾値もその分布から較正し直して訓練をやり直しても、O(3)側の誤検出率は、この理想化した座標ではすべての核・すべてのseedで正確にゼロのままだが、SO(3)側の誤検出率は平均0.65ポイント動く——原典はこれをseed間のばらつきより小さい値だと書き添えている。

同じ形の所見が、別の厳密な制約でも独立に報告されている。Orb-v3の著者らはこの監査を名指しで引き、だからこそ保存版と非保存版の両方を公開したと書いている2。力がエネルギーの微分として定義される——保存力である——という制約を外すと計算は速くなるが、構造最適化の収束が定まらなくなり、種類の異なる分子動力学で不安定化が起きると、ICML 2025でオーラル発表されたBigiらの監査は報告する4。著者らは「エネルギー保存の欠如を監視し補正することの難しさを踏まえると、直接予測した力は細心の注意を払って使うべきだ」と書き、まず直接予測した力で事前学習し、そのあと逆伝播で微調整するという折衷案を勧める4。評価時には両方の力を併用すれば、非物理的な効果を避けながら直接力の計算効率をほぼそのまま享受できる、とも書いている4。原典自身がこの区別を書いている——回転対称性は推論時に監視・補正でき、データ拡張でも学べるが、非保存性の評価にはヤコビアンの明示的な計算が要る、というのが著者らの言い分である4。制約の中身はパリティとはまったく別——エネルギー保存であって鏡映対称ではない——なので、これは同じ現象の再現ではなく、別の厳密な制約が失われたときに同じ形の失敗が起きるという、独立した収束の証拠として引いておく。

ゼロを訓練データに入れても直らない、そして万能の判定でもない

訓練データにゼロという答えを混ぜれば直るのではないか、という直感は自然である。パリティの論文はこれも測っている。SO(3)腕にゼロを明示的に教えるため、別に用意した1,000種の中心対称結晶に厳密なゼロをラベルとして与え、訓練セットへ加えて再学習した「拡張あり」の腕を作った。結果、違反の大きさは縮む——NequIPでは、拡張セットと同じ空間群に属する評価用の結晶(原典の「seen」分割)の中央値が8.26×10⁻¹から4.49×10⁻¹へ下がる——が、誤検出の割合はほとんど動かない。明示的なゼロを教えた当の1,016件(注入した1,000件と、もともと訓練セットにあった実結晶16件)に対してすら、拡張ありの腕はNequIPで89.5%、Allegroで89.9%、MACEで89.6%の割合で、まだ閾値を超えた値を返す。教師データとしてゼロを見せても、鏡映を区別できない特徴量からは、鏡映で符号が反転するはずの応答を表現しようがない——量を小さくすることはできても、その量をゼロに固定する仕組み自体が模型の側に無いためだ。原典は一文でこう書いている。「明示的なゼロで訓練しても厳密性は取り戻せない」1。ここでいう厳密性は精度とは別のものだ——精度は近づける対象だが、厳密性は満たすか満たさないかしかない。

同じ理由が、既存の汎用モデルを土台に使う場面にも及ぶ。凍結した汎用ポテンシャルの上に新しい物性を出す読み出しヘッドだけを載せて訓練する構成は、計算コストの面で広く使われている。原典はこれを凍結バックボーンで実際に測っている。パリティのラベルを持つMACE-MP-0の凍結特徴に圧電ヘッドを載せると、違反はすべての結晶とseedで厳密にゼロになった。全次数を偶と宣言するeSEN-30M-OAMでは、同じ構成が本編と同じ割合で誤検出した1。土台がSO(3)なら、その上に載せたヘッドをどれだけ訓練しても、鏡映を区別する能力は生まれない。

では、手元の模型がO(3)等変とSO(3)等変のどちらなのかは、どうやって確かめるのか。原典が示す手続きは単純で、初期化直後の模型に、ある構造とその鏡映像を入力し、出力が符号を反転させて一致するかどうかを見る——一回の鏡映操作で、ラベルは数秒で検証できる1。訓練を待つ必要はない。

ただし、パリティの有無による差がどの結晶にもどの物性にも現れるわけではない。原典は対照実験で切り分けている——双極子は圧電テンソルと同じくパリティ奇だが、階数1ではギャップがゼロか小さく、QM9の10,831分子のうち制約が働くのは7分子(約0.065%)にすぎないため、O(3)腕とSO(3)腕の成績が一致する1。分けているのは対象のパリティではなく、評価する結晶でギャップが非ゼロかどうかである1。原典自身、O(3)腕も生の座標のうえでは完全に無欠ではないと認めている。もっとも厳しい閾値(0.0001 C m⁻²)では、O(3)腕の誤検出率は2.45〜2.82%まで上がる。原典はこれをモデルの誤りではなく、データそのものの読みだと説明する——2,000種の生の結晶構造は、選定に使った緩い許容誤差(10⁻³)では全件が中心対称と確認されるが、厳しい許容誤差(10⁻⁴。こちらは予測の閾値ではなく、構造を中心対称と判定するときの座標の許容誤差である)では1,956件(97.8%)、さらに厳しい10⁻⁵では1,869件(93.5%)しか中心対称と確認できない。つまりこの2.45〜2.82%は、データに実在するわずかな対称性の破れを模型が正しく検出している分だ、というのが原典の読みである1

もう一つの限定は結晶の種類にかかわる。m3̄mという点群に属す結晶では、回転の対称性だけで既にこの応答がゼロに強制されており、パリティの有無は関係なくなる。原典はこの状態を「パリティ・ギャップが消える」と表現し、SO(3)腕もここでは正しい。パリティ・ギャップという基準は、まさにこの、群論だけから計算できる「どの物性がどの結晶で実際に晒されているか」を求めるためのものである1。原典が導くのは、この種の模型を使うなという禁止ではなく、区別である——回転についての対称性の誤りは近似の誤差として扱えるが、ゼロを強制する種類の制約は近似ではない、という一線だ1

訓練の前に、何を確認するか


出典4件
  1. Can Polat, Mustafa Kurban, Erchin Serpedin, Hasan Kurban, “A single design choice determines whether machine learning models of materials make physically impossible predictions”(arXiv:2608.18714, 2026年8月19日公開・査読前)。所属はTexas A&M・Texas A&M Qatar・Ankara University・Hamad Bin Khalifa University。NequIP・Allegro・MACEの3核それぞれにO(3)腕とSO(3)腕を、加えてEqV2(SO(3)腕のみ)を用意し、2,000種の中心対称結晶の圧電テンソルを閾値0.01 C m⁻²で判定した。O(3)腕は3核×2,000結晶×3seedの18,000回全実行で閾値を超えず、中央値は浮動小数点演算の床、最大の単一予測でも閾値の7.6分の1にとどまった。SO(3)腕は結晶の90〜96%で閾値を超え、その中央値は訓練データ中の本物の圧電応答と同じ大きさだった。訓練データにゼロを混ぜた「拡張あり」の腕でも、当のゼロに対する誤検出率はNequIP89.5%・Allegro89.9%・MACE89.6%とほとんど下がらず(違反の大きさだけ縮む)、凍結した汎用ポテンシャルに読み出しヘッドを載せる構成でもヘッドは土台の対称性群を受け継ぐと報告する。初期化直後の模型に一回の鏡映を通すだけでラベルは数秒で検証できるとも書く。著者ら自身の限定として、O(3)腕も生の座標では最も厳しい閾値(0.0001 C m⁻²)で2.45〜2.82%の誤検出が残ること(これはデータ中のわずかな対称性の破れの検出だと著者らは読む)、m3̄m点群の結晶では回転だけで既にゼロが強制されパリティ・ギャップが消えることを挙げる。この論文が導くのは特定の設計を禁止する主張ではなく、回転についての対称性誤差(近似)と厳密にゼロを強制する制約(非近似)を区別する基準(パリティ・ギャップ)である。https://arxiv.org/abs/2608.18714 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32

  2. Benjamin Rhodes, Sander Vandenhaute, Vaidotas Šimkus, James Gin, Jonathan Godwin, Tim Duignan, Mark Neumann, “Orb-v3: atomistic simulation at scale”(arXiv:2504.06231, 2025年4月8日公開。Orbital Materials 社による自社モデルの技術報告で、査読誌・会議への掲載記録は無い)。非等変・非保存力の汎用原子間ポテンシャルOrb-v3を示し、「非等変・非保存の構成でも、ポテンシャルエネルギー面の高階微分を要する性質を含め、物理的性質を正確に模型化できる」と述べる。レイテンシを10倍超、メモリ使用量を8倍超削減しながら、多くの評価でほぼ最先端の性能を維持したと報告する。本稿では、等変性を外す設計判断が実務で広く選ばれている理由として引いた——ただしこの論文が測るのは精度と速度のトレードオフであり、対称性で厳密にゼロが強制される量を外せてしまうかどうかは扱っていない。https://arxiv.org/abs/2504.06231 2 3 4 5 6 7

  3. Johann Brehmer, Sönke Behrends, Pim de Haan, Taco Cohen, “Does equivariance matter at scale?”(arXiv:2410.23179, 2024年10月30日投稿・v2は2025年7月27日。TMLR掲載)。剛体どうしの相互作用を汎用の変換器で学習させ、模型規模・学習ステップ数・データ量を振って比較した。等変性はデータ効率を上げるが、「十分なエポック数があれば、データ拡張を使った非等変模型の訓練でこの差を埋められる」と報告する。計算予算を揃えた比較では等変模型が各予算で非等変模型を上回り続け、模型規模と学習時間への予算配分の最適点は等変・非等変で異なるとも書く。本稿では、回転についての対称性は学習で埋め合わせうるという、本稿の見立てに対する最も一般的な反証材料として引いた——原典が焦点に据える群は、並進・回転・鏡映からなる E(3) で、等変側のアーキテクチャ GATr は鏡映まで含めて等変である。ただし差を埋めたデータ拡張のほうは、2次元の並進と重力軸まわりの回転だけを掛けており、鏡映は含まれていない。対称性が量を厳密にゼロへ強制する場面も扱っていない。https://arxiv.org/abs/2410.23179 2 3 4 5

  4. Filippo Bigi, Marcel Langer, Michele Ceriotti, “The dark side of the forces: assessing non-conservative force models for atomistic machine learning”(arXiv:2412.11569, 2024年12月16日投稿・v6は2025年10月25日。ICML 2025 Oral採択)。力をエネルギーの微分として定義しない非保存力モデルは計算が速い一方、幾何最適化の収束が不定になり、複数種類の分子動力学で不安定化を起こすと報告する。著者らは「エネルギー保存の欠如を監視し補正することの難しさを踏まえると、非保存力は注意深く使うべきだ」と書き、非保存力で事前学習しその後保存力で微調整する折衷案を勧める。本稿ではパリティとは別の厳密な物理制約(エネルギー保存)が失われたときに同じ形の失敗が起きるという独立した収束の証拠として引いた——同じ結論の再現ではなく、制約の種類が異なる別の監査である。https://arxiv.org/abs/2412.11569 2 3 4

この記事はAIが執筆しています。内容には誤りが含まれる可能性があります。ご注意ください。