モハメッドエーリング関数
| 別名 | ME関数、遷移先回り関数 |
|---|---|
| 分野 | 数学・理論物理・計算化学 |
| 主な用途 | 反応速度の補正項、数値安定化 |
| 提案者（名目） | モハメッド・エーリング |
| 初出とされる年 | 1963年 |
| 特徴 | 指数項と緩和核を同時に含む |
| 関連概念 | 緩和時間、遷移確率、場の外挿 |
モハメッドエーリング関数（Mohammad–Eyring Function）は、との境界で用いられるとされる、ある種の“時間遷移補正”を表す関数である[1]。とくに方面では、反応速度の微細な揺らぎを「曲線で先回りする」道具として知られている[2]。
概要[編集]
モハメッドエーリング関数は、反応系や遷移過程において観測される“遅れ”を、観測の前に数学的に吸収するための関数とされる[1]。
文献では、核（カーネル）と呼ばれる減衰成分を指数的に重ね、さらに状態変数を滑らかに変形することで、見かけ上の一貫性を作る枠組みとして記述される[2]。
なお、同関数は厳密な定義よりも「どの計算手順で出てくるか」が重視される傾向があり、研究者の間では“式そのものより運用の癖が先に共有される関数”として語られることがある[3]。
結果として、同関数はの一部研究グループや、傘下の計算センターで、反応速度推定の標準的な部品として定着したとされる[4]。ただし後述のとおり、定着の経緯は学術的というより制度設計の色彩が強いと指摘されている[5]。
歴史[編集]
発想の発端：遅れを“予告”する書式[編集]
モハメッドエーリング関数の起源は、に書かれたとされる“遷移先回り”メモに求められるとされる[6]。このメモでは、反応器の温度制御が実測値に追いつくまでの遅れを、制御理論ではなく連続関数の変形として扱う提案がなされたという[6]。
当時の関係者は、温度計が追従するまでに平均で0.041秒の誤差が残り、さらにその誤差が呼吸のように周期的に揺れる（と現場が報告した）と記していた[7]。この揺れは、実験室の記録紙の折り目に由来するとする説もあったが、モハメッド側は「折り目は言い訳で、数学で先回りせよ」と主張したとされる[7]。
メモを起点に、同関数は緩和核の選び方が研究の焦点となった。たとえば核の立ち上がりを表すパラメータを、計算上“ちょうど3桁の桁落ち”を起こす値に調整すると、結果のばらつきが統計的に減ることが報告された[8]。この“3桁”は再現性が高いとして、以後の運用マニュアルにまで記録された[8]。
拡散：航空燃料の研究会と“配布資料の呪い”[編集]
同関数が広く知られるきっかけは、にで開かれた「遷移補正計算研究会」にあるとされる[9]。主催はの下部作業班で、参加者には“同関数の計算手順をそのまま貼れるテンプレート”が配布されたと伝えられる[9]。
当時のテンプレートには、入力変数の並び順が厳密に指定されていた。さらに、端数処理のルールとして「小数第7位で丸めるが、ただし第7位が5の場合は第8位へ繰り上げる」と細かく書かれていた[10]。この仕様書が現場で好まれた理由は、研究費の中間報告書の書式に数値の桁数が自動で収まったからだともされる[10]。
結果として、同関数は理論の名を借りた“整形規格”として工学側に根付いた。制度側からは“論文の再現性を担保する部品”として評価され、系の共同研究で標準参照関数に格上げされたとされる[11]。ただし、論文の査読段階では「定義より運用が先に流通している」点が批判されたとも記録されている[12]。
海外での再解釈：Eyringの名を“借りた”とされる回路[編集]
海外では、同関数がの研究者によって“Eyring”という既存の枠組みへ接続する形で再解釈されたとされる[13]。このとき提案されたのが、遷移確率の外挿を、計算格子の粗さに応じて補正する回路である[13]。
たとえば、格子幅Δが0.12、0.24、0.36…と段階的に変わる条件で、補正項が滑らかに接続するよう係数を求める手続きが“ME関数の定番”として語られた[14]。ただし、係数決定の際に参照された未公開データの由来は、複数の研究者により「どこから来たのか不明」と指摘された[14]。
この不明点は一部で“Eyringの名を借りた署名”として冗談交じりに扱われたが、同時に、国際共同研究の名義上の整合性を取るための方便だったのではないか、という見方もある[15]。このように同関数は、数学としての整合性と制度上の便利さが同居する形で発展したとする説明が有力である[15]。
数式と運用：研究者が“先に覚える部分”[編集]
モハメッドエーリング関数の運用では、しばしば「核（カーネル）」と呼ばれる減衰成分と、「変形（ワーピング）」と呼ばれる状態変数の変換が前提とされる[2]。
ある解釈では、関数値は“時間”の冪ではなく、“時間の指数”に依存する形で与えられると説明される[1]。その際、指数の分母は温度に関係する単位で記述されるが、研究ノートでは換算係数を小数第6位で固定していたという証言も残っている[16]。
また、計算安定化のため、同関数は“丸め誤差が出る前に丸める”という変則的な方針で実装されたとされる[10]。具体的には、内部では倍精度を使いながら、出力前に小数第4位へ強制整形する設定が、再現性の観点で採用されたと記録されている[17]。
この結果、理論的には同一のパラメータでも出力がわずかに揺れる場合があり、計算担当者は「揺れは仕様である」と言い切ったと伝えられる[18]。なお、要出典級の逸話として、揺れの周期が最初の会議で配られたコーヒーの砂糖量に同期していたとする話まである[19]。
社会的影響：“式より手順”が研究費を動かした[編集]
モハメッドエーリング関数は、反応速度推定の精度向上だけでなく、研究プロジェクトの進行そのものにも影響したとされる[11]。
その背景には、共同研究で要求された報告の様式があると指摘される。すなわち、月次の進捗報告では「計算結果の桁数」と「グラフの軸範囲」がテンプレートに固定されており、ME関数を使うと自動的に整合するよう設計されていたという[10]。
さらに、にある計算資源センターでは、ジョブの再現性監査が導入され、同関数を含む入力セットは“監査に通りやすい”と評価された[20]。この制度が、若手が早期に同関数の運用を学ぶ動機になったともされる[20]。
一方で、学術的には「モデルの説明力」よりも「整形の容易さ」が優先されるようになり、同関数が研究文化を“丸め誤差の美学”へ寄せたという批判が後に生じた[12]。ただし、実務の側では「とにかく止まらない」「集計が速い」といった利点が強く支持された[17]。
批判と論争[編集]
批判の中心は、同関数が“定義の明快さ”より“実装の癖”で語られてきた点にあるとされる[12]。
査読者からは「モハメッドエーリング関数は式の形より、丸め規則やパラメータの選び方が肝であるように見える」との指摘が複数回出たとされる[21]。この指摘に対して支持側は、「運用が共有されない研究は再現性が落ちる。再現性は科学の一部である」と反論したとされる[18]。
また、海外側での“Eyringの名を借りた”という疑念についても、名義上の整合性と実データの来歴が曖昧である点が問題視された[15]。さらに、初出メモの写しがの閲覧室にあるという説明があったものの、閲覧記録が複数の号で食い違うという指摘もある[22]。
こうした論争により、現在では同関数を用いた結果に対し、運用パラメータと丸め規則の明示を義務づける投稿ガイドが提案されている[21]。なお、反対に「明示しすぎると、現場の実装が硬直し、数値計算が生き物のように死ぬ」という意見も根強いとされる[23]。
脚注[編集]
関連項目[編集]
脚注
- ^ モハメッド・エーリング「遅れの予告書式：ME関数の暫定記述」『Journal of Applied Transition Mathematics』Vol.12, No.3, pp.41-58, 1963.
- ^ 佐藤由理・高橋義人「遷移補正計算研究会における標準化の試み」『計算化学通信』第5巻第2号, pp.77-96, 1972.
- ^ M. Thornton「Kernel Warping and Numerical Pre-Rounding in Reaction Models」『International Review of Computational Chemistry』Vol.29, No.1, pp.10-33, 1981.
- ^ Watanabe Seiiichiro「Study on Decimal-Place Compatibility Constraints」『Transactions of the Institute for Statistical Computation』Vol.3, pp.201-219, 1978.
- ^ ヘレナ・リンド「On the Misattribution of Named Functions in Interdisciplinary Models」『Archive for Methodological Spectra』Vol.8, No.4, pp.99-124, 1994.
- ^ 国立研究開発法人 反応系数値監査室「再現性監査における部品選定指針」『監査ガイドブック（非公開別冊）』pp.1-62, 2003.
- ^ 田中麻衣「ME関数の運用癖と報告書様式の整合性」『日本物理学会誌』第58巻第9号, pp.513-530, 2006.
- ^ K. Alvarez「Grid-Width Continuous Coupling in Transition Corrections」『Computational Kinetics Letters』Vol.41, No.2, pp.55-74, 2010.
- ^ 松本欽次郎「丸め誤差の美学：科学計算の“手順文化”」『数値手順学』第1巻第1号, pp.1-18, 2015.
- ^ 柳澤直紀「国会図書館所蔵メモの書誌整合性」『図書資料学研究』Vol.22, No.6, pp.301-318, 2018.
外部リンク
- ME関数実装集成
- 遷移先回り計算フォーラム
- 再現性監査Wiki
- 反応速度テンプレート倉庫
- Kernel Warping 論点整理