3910 words
20 minutes
【計算化学】自作pythonモジュールで遷移状態構造を求めてみる(Tris(trifluoromethyl)phosphineが2つ配位したPd触媒によるtrifliorotoleneの還元的脱離, NNP使用）
最終更新：2025-08-11
概要
本記事では、自作モジュール(MultiOptPy)で、Tris(trifluoromethyl)phosphineが2つ配位したPd触媒によるtrifliorotoleneの還元的脱離の遷移状態構造を算出してみる。以前の記事で反応を扱ったため、モデル系中の配位子の種類を変えた。使用する計算レベルは、Meta社のFAIR Chemistryが開発したニューラルネットワークポテンシャル(NNP)であるUMA(Meta’s Universal Model for Atoms)を使用する。
MultiOptPyは電子状態計算ソフトウェアを用いた分子構造最適化手法の勉強を目的として作成したpythonモジュールである。
MultiOptPyのレポジトリ:https://github.com/ss0832/MultiOptPy
還元的脱離について：
今回使用したニューラルネットワークポテンシャルについて：
- https://ai.meta.com/blog/meta-fair-science-new-open-source-releases/ (UMAの公開に関する記事)
- https://github.com/facebookresearch/fairchem (FAIR Chemistryの提供するGitHubのレポジトリ)
- https://fair-chem.github.io/ (同上のレポジトリの内容に関して説明したサイト)
- https://huggingface.co/facebook/UMA (NNPの配布サイト, Hugging Faceへのアカウント登録と配布元の使用許諾が必要である。)
使用した自作モジュールMultiOptPyのバージョン
v1.10c環境
Windows 11※Windows 11環境下でAnaconda PowerShell Promptを使用した。
Source codeのダウンロード
wget https://github.com/ss0832/MultiOptPy/archive/refs/tags/v1.10c.zip
unzip v1.10c.zip
cd MultiOptPy-v1.10c環境構築手順
今回は、Windows 11のPower Shellを使用した。初めに、NNPを使用できる環境が整ったAnaconda PowerShell Promptを用意する手順を説明する。
1, https://repo.anaconda.com/archive/ より、Anaconda3-2025.06-1-Windows-x86_64.exeでAnacondaをインストールする。
2, 検索機能を使い、スタートからAnaconda PowerShell Promptを開く。
3, 以下のコマンドを実行し、仮想環境を作成する。
conda create -n (任意の仮想環境名) python=3.12.74, 先ほど作成した仮想環境をconda activate (仮想環境名)で起動させる。
5, 以下のコマンドを実行し、必要なライブラリを導入する。
pip install fairchem-core==2.2.0 ase==3.25.0 torch==2.6.0- fairchem-coreは、FAIR Chemistryが管理しているNNPを動作させるために必要なライブラリである。
- aseはNNPに電子エネルギーを算出したい分子構造を渡すために必要なインターフェイスの役割を果たすために必要なライブラリである。
- torchはPyTorchライブラリを指す。これはニューラルネットワークなどの機械学習を行ったり、学習結果を扱ったりするために必須なライブラリである。
これで、Anaconda PowerShell Promptから仮想環境を立ち上げることで、NNPを使用する準備が整えることが出来る。
次に、NNPを使用するために必要なModelの情報が保存されている.ptファイルのダウンロードおよびNNPの自作モジュールへの導入方法について説明する。
1, 以下のサイトにアクセスして、uma-s-1p1.ptをダウンロードする。（使用許諾が下りていれば可能である。）
https://huggingface.co/facebook/UMA2, ダウンロード後、MultiOptPy-v1.10cディレクトリ内に存在するsoftware_path.confに対して、uma-s-1p1.ptの絶対パスを用いて以下を追記する。
uma-s-1p1::(uma-s-1p1.ptの絶対パス)これで、MultiOptPy-v1.10cがNNPuma-s-1p1を使用できるようになる。
使用するNNPに関する具体的な説明
今回使用するNNPについて具体的に説明する。
- UMAのModel Checkpointは
uma-s-1p1を使用した。
IMPORTANTモデルをダウンロードした時点(2025/8/9)で、
uma-s-1のモデルが存在していたが、https://huggingface.co/facebook/UMA より、” *Note uma-s-1 is now deprecated and will be removed by the next release. uma-s-1p1 should be used instead for a replacement.
とあり、今後のアップデートで削除されるため、使用しなかった。
- 小分子系のトレーニングセットである
Omol25(omol)を使用して学習したニューラルネットワークポテンシャルを使用する。
※自作モジュールでの具体的な使用の仕方に関しては、https://github.com/ss0832/MultiOptPy/blob/main/multioptpy/ase_calculation_tools.py を参照
手順
1. 初期構造の準備
モデル反応系として、以下の構造を用意した。今回はファイルの名前をreductive_elimination_PCF3_3.xyzとした。 初期構造は以下のものを使用した。
42
Pd 0.624705911130 0.010022575618 0.084070998866
C 0.487303578640 -2.000829163119 0.102597066303
C 0.509790270624 -2.691517337181 -1.102471727478
C 0.441750574619 -2.692231805035 1.305395515460
C 0.478103697787 -4.079428238082 -1.100183198149
H 0.565763420690 -2.161065201006 -2.046255557755
C 0.403570438263 -4.080957609155 1.299385933616
H 0.447152766089 -2.162719659892 2.251043081517
C 0.420940053562 -4.775235272887 0.099104199788
H 0.504226884160 -4.616965029999 -2.039575452427
H 0.365284726057 -4.618936121418 2.238147661773
H 0.396565900941 -5.856873773845 0.098529558682
P -1.735430551859 -0.132864095825 0.011170254099
P 1.199800850586 2.375500069289 -0.033560217738
C 2.647677679794 -0.414273856006 0.202062589882
C -2.770739824663 1.097910629827 1.018032795230
C -2.490620636678 0.087831326166 -1.710953500444
C -2.691322149589 -1.702938844224 0.521767547397
C 0.710389011617 3.445951599701 1.450512946094
C 0.375630544452 3.313315432777 -1.456489825819
C 2.948002515938 3.096059529016 -0.285392891784
F -1.817754408572 1.036822477605 -2.366327894092
F -3.778575395105 0.424445261247 -1.703690430788
F -2.356549793721 -1.060016902901 -2.384809636816
F -4.081401526924 0.886417219111 1.035167905663
F -2.561754392383 2.324871456904 0.527607979852
F -2.331183683198 1.055961311887 2.281167214110
F -2.272454693614 -2.790782830446 -0.103735637869
F -3.997931483525 -1.577969795217 0.262553541805
F -2.544558580705 -1.865510676654 1.841431203998
F -0.338274396657 2.895946862214 2.067820084475
F 1.728124520806 3.472693939012 2.316505710855
F 0.390706313251 4.698233210912 1.130400975465
F 0.528931977614 2.584090785580 -2.567547100099
F -0.931676288834 3.434400930275 -1.203760807064
F 0.862456419356 4.527498822813 -1.691018511093
F 3.317077893227 2.890201971011 -1.551134047869
F 2.955804450903 4.416685776887 -0.063355588252
F 3.835378686011 2.550536782870 0.530365731053
F 3.326943915839 0.237547330684 -0.793760032889
F 3.122339736010 0.091740833812 1.377319657217
F 3.105805068059 -1.673569922325 0.151861905226
2. 遷移状態構造最適化のための初期構造の算出
a. 初期構造をカレントディレクトリにreductive_elimination_PCF3_3.xyzとして保存し、以下を実行する。
python .\optmain.py .\reductive_elimination_PCF3_3.xyz -os uma-s-1p1 -opt rsirfo_fsb -ma 300 2 15-opt rsirfo_fsbは準ニュートン法であるRS-I-RFO法を構造最適化に使用することを示す。初期のへシアンに関しては、特にオプションで指定しない限り、単位行列が使われる。-spinはスピン多重度の指定である。PySCFを使用するときは目的とするスピン多重度に1を引いた値を指定する。（デフォルトでは1が指定される。）-elec Mは形式電荷をMとすることを示す。（デフォルトでは0が指定される。）-ma yyy a bはyyykJ/molの活性化障壁を超えうるペア同士を近づける力を原子のラベル番号aとbのペアに構造最適化時に加えることを示す。-os uma-s-1p1は今回使用するNNPを指定している。これを使用する際にASEライブラリが必要である。
これを実行すると、omolのデータセットを使用したuma-s-1p1モデルのNNPで得たエネルギーに対して、指定した人工力ポテンシャルを加えた上で初期構造を構造最適化することができる。
結果はyyyy_mm_dd（今日の年月日）ディレクトリの中に存在するディレクトリを開いて確認できる。
正常終了していれば、このディレクトリ中に、reductive_elimination_PCF3_3_traj.xyzが存在するので、これをコピーして、MultiOptPy-v1.10cディレクトリに置く。
reductive_elimination_PCF3_3_traj.xyzは構造最適化の過程をAvogadro（公式ページ:https://avogadro.cc/ ）等で可視化して確認できるようにしている。このreductive_elimination_PCF3_3_traj.xyzは次のNEB計算に使用する。
※reductive_elimination_PCF3_3_traj.xyzをアニメーションとして表示したい場合は、[https://github.com/ss0832/molecule_movie] を使うと良い。
（この人工力ポテンシャルを加えて行った構造最適化の結果はreductive_elimination_PCF3_3_optimized.xyzで確認できる。構造を可視化して、生成系になっているか確認する。反応系のままであれば、-maの設定を見直してb.をやり直す。今回の場合以下の構造が得られた。生成系の安定構造に近いものが得られていると判断できるため、次のNEB計算を行うことが可能である。）
42
OptimizedStructure
Pd 0.188291979218 -0.011091727738 -0.001771164255
C 1.950147674095 -1.838646573926 0.304412951826
C 1.547622109418 -2.180813028761 -0.993298913061
C 1.294852691360 -2.416251472915 1.400543737599
C 0.565081045479 -3.146304712042 -1.180514113852
H 2.056294212271 -1.741783275648 -1.841176923393
C 0.323047847298 -3.382334694513 1.200346094313
H 1.599658037796 -2.144339853186 2.402185538134
C -0.011083298088 -3.775773692598 -0.088585245295
H 0.274418217133 -3.422899499356 -2.184900503001
H -0.159216664186 -3.842280444737 2.052540193260
H -0.743888129543 -4.555458448179 -0.240294675675
P -2.057937298457 -0.288822731255 -0.022275709757
P 1.087741519614 2.078548689607 -0.108462903341
C 3.173075510363 -1.148940225277 0.518749083041
C -2.969986219803 0.933794142928 1.104988214020
C -2.942991081404 0.090948754327 -1.655271452691
C -3.076262603798 -1.823495239393 0.439159719730
C 0.699521203897 3.251053432559 1.330569580762
C 0.294892184771 3.057386542622 -1.525189183753
C 2.870830038163 2.673762158310 -0.394298168237
F -2.296968500377 1.057531297269 -2.309247123206
F -4.216437747251 0.467226692125 -1.522778702981
F -2.915372592402 -1.014148413385 -2.419532113483
F -4.269519793974 0.694779486359 1.285960572427
F -2.848027088343 2.163109278779 0.588743923451
F -2.387808040865 0.917015237463 2.310688187041
F -2.661165784098 -2.879039812033 -0.274374021992
F -4.392224771613 -1.722886192344 0.251885915919
F -2.860214706424 -2.090913552171 1.738203291843
F -0.454672712795 2.907301075930 1.906144342435
F 1.668812810564 3.130146603416 2.250446529109
F 0.613738495789 4.537451862833 0.985301609121
F 0.301575406035 2.292931591246 -2.626378920471
F -0.979184239635 3.321187355052 -1.204423725932
F 0.885863420906 4.210647330060 -1.835676596292
F 3.240069804744 2.301285085792 -1.629734124156
F 3.034492326505 3.996329930470 -0.296545113967
F 3.698286399644 2.109491028071 0.484997032955
F 3.544306685999 -0.397570734705 -0.535308540653
F 3.142631057491 -0.347227846778 1.602110048754
F 4.187710594502 -2.020905404276 0.732061373702
b. reductive_elimination_PCF3_3_traj.xyzを初期パスとして、NEB法で経路の緩和を行う。
NEB法を用いることで、先ほど得られたreductive_elimination_PCF3_3_traj.xyz全体のエネルギーを下げることができる。これにより、パスのエネルギー極大値を持つ構造を遷移状態構造に近づける。（この時点ではまだ正確な遷移状態構造は求められていない。）
python .\nebmain.py reductive_elimination_PCF3_3_traj.xyz -ns 10 -modelhess -spng -nd 0.5 -os uma-s-1p1-nd Nはノード間の距離をN Åとして初期パスを作成することを示す。-ns nはn回分NEB法による経路の緩和を行うことを示す。-fc MはM回あたりの経路緩和回数に対して1回だけ正確なHessianを計算し、経路緩和に使用する。これを使用すると、Hessianを使用しない場合の経路緩和アルゴリズムとは別のものを使用して、経路緩和を行う。-spngは緩和中のパスのエネルギープロファイルや各ノードの勾配のRMS値をmatplotlibで可視化するオプションである。-os uma-s-1p1は今回使用するNNPを指定している。これを使用する際にASEライブラリが必要である。
c. 初期構造の決定及び遷移状態構造の計算
MultiOptPy-v1.10cと同じディレクトリ内に、NEBという名前を含むディレクトリが生成されている。 そのディレクトリ内のenergy_plot.csvを確認し、緩和後のパスのエネルギー極大値を示す構造を確認する。
パスの緩和後の各ノードのエネルギー一覧(単位
※bias_force_rms.csvにて、各Iterationごとのすべてのノードの勾配のRMS値を確認できる。
私が実行した環境では、4番のノード(グラフでは5番)がエネルギー極大値を示していた。
※こちら[https://ss0832.github.io/molecule_viewer/] を使うことでも可視化は可能である。
reductive_elimination_PCF3_3_traj_4.xyz
42
0 1
Pd 0.452521000220 0.119737362716 0.075081092663
C 0.971655189656 -1.929401927912 0.118010730292
C 0.820702695396 -2.597796521050 -1.097063011003
C 0.753357646877 -2.598108307190 1.324667979064
C 0.551704073324 -3.952376179768 -1.092422865762
H 0.858846255558 -2.076419179134 -2.045658728076
C 0.481705052501 -3.951152650042 1.302879425118
H 0.734093505695 -2.076081849080 2.273700573485
C 0.434978003158 -4.655294976251 0.101947560536
H 0.409573857801 -4.466155737658 -2.035198182454
H 0.282574829368 -4.466066178602 2.234837220375
H 0.289883924920 -5.728039437820 0.097893861456
P -1.811020264823 -0.110568767032 0.013141600302
P 1.177043221825 2.389699282943 -0.040014017853
C 2.415313392845 -0.888728504389 0.193782441658
C -2.833960556817 1.103499666565 1.020099764347
C -2.540208650192 0.097373038472 -1.711329548161
C -2.700381047866 -1.692088775257 0.519517904420
C 0.710329097559 3.466959880826 1.446615043372
C 0.371828533690 3.344621604639 -1.459450830412
C 2.933018671601 3.071474397988 -0.285914785203
F -1.855697628704 1.037357108088 -2.369723339461
F -3.829274927666 0.433068092412 -1.734168272419
F -2.392165288870 -1.057954828866 -2.378958429699
F -4.147358556789 0.896440023217 1.055394479860
F -2.621334499853 2.331469707825 0.530476364444
F -2.378496442145 1.061770594528 2.280622078081
F -2.185523841839 -2.754883772683 -0.099597518696
F -4.008892344512 -1.656357200102 0.248269842850
F -2.552686222419 -1.863106226236 1.840882366502
F -0.338622302799 2.916453718628 2.066101973230
F 1.731968570288 3.479725139017 2.314725441178
F 0.398118923167 4.727380486887 1.146586703437
F 0.526576441650 2.611358993341 -2.571842836597
F -0.936748394711 3.463144259061 -1.209727800664
F 0.855939242918 4.560482237171 -1.704253989815
F 3.316940541812 2.852853312276 -1.548479148897
F 2.993405670422 4.388698035215 -0.055664978926
F 3.801061609742 2.473700786386 0.528961857616
F 2.895365064690 -0.120093756340 -0.840382318771
F 2.696563338903 -0.271113519279 1.374690526246
F 3.267302614419 -1.915479433510 0.170963772336
構造が壊れていないので、これを遷移状態を求めるための初期構造とする。
reductive_elimination_PCF3_3_traj_4.xyzをMultiOptPy-v1.10cと同じディレクトリ内にコピーする。
そして、以下を実行する。
python .\optmain.py reductive_elimination_PCF3_3_traj_4.xyz -freq -tcc -opt rsirfo_bofill -fc 5 -order 1 -os uma-s-1p1-opt rsirfo_bofillは遷移状態構造の最適化向けのoptimizerを指定することを意味する。準ニュートン法であるRS-I-RFO法を使用する。今回は-fcで正確なへシアンを計算するようにしているので、初期へシアンは正確なへシアンを使用するようになっている。-order 1は一次の鞍点を求めることを指定する。（デフォルトだと極小値を求めるようになっている。）-fc 5は5回の反復回数当たり1回正確なへシアンを計算することを指定する。-freqは収束条件を満たした後に基準振動解析を行うことを示す。（自前で実装しているため、あくまで目安として使用することを推奨する。各振動モードをvibration_animation内のxyzファイルで可視化できる。）UMAモデルから算出されるHessianは数値微分により求めているため、原子数Zが多いとZの二乗オーダーで計算コストが急増する。-tccは収束条件を厳しくすることを示す。（Gaussianのtightと同等）
実行して得られた正確な遷移状態構造を以下に示す。
(実行して得られた正確な遷移状態構造はreductive_elimination_PCF3_3_traj_4_optimized.xyzとして保存されている。)
42
OptimizedStructure
Pd 0.202478944346 0.029997291182 0.214184503861
C 1.036833667458 -1.862375258332 0.296182324628
C 1.918399279462 -1.916306992002 -0.775427988218
C 0.354285602192 -2.986002848747 0.731603325021
C 2.036177489529 -3.103765633798 -1.489614571862
H 2.525781851843 -1.059713044045 -1.037916681742
C 0.499623149226 -4.168912632242 0.023589895954
H -0.286475842524 -2.941842773908 1.600590131984
C 1.327523945042 -4.225459859219 -1.091275975141
H 2.703995566745 -3.150566446281 -2.339656854719
H -0.044669331322 -5.048532501750 0.341488587257
H 1.432807147136 -5.153097392934 -1.638002527087
P -2.111860512732 -0.080529024955 0.325915732298
P 1.082515155036 2.084740701212 -0.436653546686
C 1.607007478942 -0.579291071130 1.802938409949
C -2.772050016338 1.209146924080 1.537993485255
C -3.226429564264 0.294913298403 -1.152700720717
C -3.020128148630 -1.602706549750 0.971417183618
C 1.190015597078 3.486640959308 0.827458352209
C 0.259835781195 3.007065024995 -1.861867219827
C 2.859601825503 2.111719284794 -1.068659096922
F -2.678764284731 1.272312980062 -1.880079059643
F -4.466504743773 0.665886075105 -0.833663139796
F -3.301363185882 -0.804039146765 -1.916674844251
F -4.030135701561 1.028492419094 1.937507398254
F -2.687262594083 2.410404917677 0.947453105835
F -1.993921804262 1.218673674812 2.624751079537
F -2.655806213038 -2.673112538663 0.256610996998
F -4.349689138219 -1.522581124011 0.932473590991
F -2.650062029941 -1.802895758734 2.245811106427
F 0.112293709354 3.445009675934 1.615047822207
F 2.277645224012 3.320270694394 1.584483360312
F 1.250533062111 4.700292570906 0.274738754863
F -0.040401286301 2.126868088781 -2.824692375195
F -0.882163666457 3.542826185471 -1.410104664678
F 0.990549831141 3.976547513202 -2.411635799998
F 2.879289437861 1.568034022679 -2.296073360153
F 3.407450193587 3.323653355307 -1.140351037599
F 3.621703728437 1.358084618663 -0.269659499644
F 2.917133035398 -0.804864130351 1.638943555750
F 1.549205357411 0.674040574610 2.354046451899
F 1.155002004014 -1.369026123053 2.789479808771
210回程度の反復計算で遷移状態構造が得られた。-freqオプションにより生成されたnormal_modes.txtやvibration_animationディレクトリ内の振動モードのアニメーションを確認した。
以下に-freqオプションで生成されたnormal_modes.txtの一部を示す。
Mode 0 1 2
Freq [cm^-1] -223.7771 10.5750 15.6855
Reduced mass [au] 10.4650 12.5469 16.7304
Force const [Dyne/A] -0.3088 0.0008 0.0024
Char temp [K] 0.0000 15.2151 22.5680
Normal mode x y z x y z x y z
Pd 0.02299 -0.01528 0.02736 -0.00237 -0.00032 -0.00454 -0.00086 0.00502 -0.00042
C -0.09968 -0.01085 -0.10199 0.00664 0.00543 0.01992 -0.00439 0.00276 -0.00086
C -0.04322 -0.02191 -0.05157 0.03172 0.00264 0.04097 -0.01686 -0.00021 -0.01125
C -0.03952 -0.01908 -0.04785 0.00083 0.00513 0.01046 0.00040 0.00422 0.00991
C 0.01697 -0.03969 -0.00555 0.05313 -0.00128 0.05066 -0.02591 -0.00171 -0.00988
H -0.07877 -0.00328 -0.07660 0.03466 0.00327 0.05036 -0.01944 -0.00113 -0.02054
C 0.01751 -0.03937 -0.00391 0.02143 0.00176 0.02054 -0.00773 0.00221 0.01130
H -0.06078 -0.00869 -0.06403 -0.01851 0.00727 -0.00393 0.00943 0.00681 0.01645
C 0.02796 -0.04410 0.00598 0.04787 -0.00152 0.04022 -0.02123 -0.00077 0.00163
H 0.03971 -0.04830 0.01253 0.07349 -0.00367 0.06678 -0.03621 -0.00394 -0.01784
H 0.04499 -0.04913 0.01685 0.01739 0.00150 0.01292 -0.00423 0.00317 0.01994
H 0.05324 -0.05322 0.02662 0.06422 -0.00428 0.04803 -0.02788 -0.00216 0.00271
P -0.00332 -0.00124 -0.00686 -0.00283 -0.00169 -0.00918 0.00033 0.00038 0.00705
P -0.00170 0.01108 0.00476 0.00319 -0.00479 -0.01327 -0.00059 0.00302 -0.00381
C -0.02500 0.09138 -0.00215 -0.03136 0.02086 0.02527 0.00998 0.00171 -0.00838
C 0.00384 0.00035 -0.00177 0.00828 -0.02258 0.01904 0.00163 -0.02439 0.03414
C -0.00494 0.00006 -0.00201 -0.01239 0.02742 0.00550 -0.00581 0.02799 0.01868
C 0.00260 -0.00289 -0.00138 -0.00209 -0.01149 -0.03148 0.00444 -0.01385 -0.02074
C 0.00206 0.00755 0.00294 0.02338 -0.00978 -0.00960 -0.02441 0.02044 -0.02105
C -0.00070 0.00314 0.00051 -0.00523 0.00576 -0.00157 0.01911 -0.01710 -0.02825
C -0.00048 -0.00463 0.00214 -0.00285 -0.01228 -0.03083 0.00969 -0.00125 0.02536
F -0.00096 0.00043 0.00023 -0.01535 0.03892 0.01885 -0.00885 0.04124 0.03421
F -0.00367 0.00123 -0.00269 -0.00932 0.02463 0.02056 -0.00429 0.02253 0.03112
F -0.00248 0.00053 -0.00193 -0.02000 0.04086 -0.01324 -0.00953 0.04215 -0.00134
F 0.00415 -0.00107 -0.00282 0.01088 -0.02870 0.02463 0.00647 -0.04134 0.04196
F 0.00296 0.00016 -0.00170 0.00548 -0.01184 0.04052 -0.01304 -0.01390 0.05347
F 0.00419 -0.00154 -0.00325 0.01609 -0.04342 0.01351 0.01086 -0.03627 0.02764
F -0.00182 -0.00144 -0.00209 -0.00106 -0.00071 -0.04715 0.01003 -0.00019 -0.03814
F 0.00268 0.00300 -0.00085 -0.00209 -0.01169 -0.03048 0.00421 -0.01795 -0.02254
F 0.00258 0.00024 -0.00205 -0.00242 -0.03015 -0.03417 0.00188 -0.03424 -0.02309
F 0.00064 0.00364 0.00110 0.02892 -0.00345 -0.00183 -0.03681 0.02776 -0.03746
F 0.00144 0.00660 0.00104 0.02802 -0.02114 -0.01867 -0.03615 0.03503 -0.00121
F -0.00044 0.00644 0.00175 0.02916 -0.00849 -0.00619 -0.01968 0.01288 -0.03730
F 0.00110 -0.00148 0.00381 -0.02114 0.01133 -0.00167 0.02375 -0.03214 -0.01606
F -0.00146 0.00161 -0.00052 0.00316 0.01191 0.01237 0.01762 -0.00234 -0.04962
F -0.00139 0.00147 -0.00112 -0.00445 0.00285 -0.00570 0.03079 -0.03110 -0.03748
F -0.00364 -0.00260 0.00185 -0.01738 -0.00954 -0.03208 0.03058 -0.01672 0.03227
F -0.00409 -0.00282 0.00117 0.00299 -0.01512 -0.03301 0.00899 -0.00135 0.01903
F -0.00140 -0.00217 0.00109 0.00115 -0.01810 -0.03962 -0.00242 0.00979 0.04710
F -0.03686 0.02330 -0.03115 -0.02363 0.05641 0.03827 0.00414 -0.03907 0.00074
F 0.02603 0.05853 0.05832 -0.07301 0.01656 0.03001 0.04738 0.01208 -0.02777
F -0.01742 0.03424 -0.03432 -0.02044 0.00468 0.01723 -0.01798 0.02946 0.00111
(...snip...)その結果、虚振動が１つであることが確認できた。
次に、vibration_animation内のmode_1_604i_wave_number.xyzをAvogadroで確認すると、想定される反応系と生成系をつなぐ方向に振動していることを確認できた。
終わりに
自作モジュールで、UMAモデルのニューラルネットワークポテンシャル(NNP)を用いて、Tris(trifluoromethyl)phosphineが2つ配位したPd触媒によるtrifliorotoleneの還元的脱離の遷移状態構造を算出する手順を説明した。
参考
- https://github.com/ss0832/MultiOptPy (自作モジュールMultiOptPyのレポジトリ)
- https://avogadro.cc/ (Avogadro、分子構造可視化ツール)
- https://ai.meta.com/blog/meta-fair-science-new-open-source-releases/ (UMAの公開に関する記事)
- https://github.com/facebookresearch/fairchem (FAIR Chemistryの提供するGitHubのレポジトリ)
- https://fair-chem.github.io/ (同上のレポジトリの内容に関して説明したサイト)
- https://huggingface.co/facebook/UMA (NNPの配布サイト, Hugging Faceへのアカウント登録と配布元の使用許諾が必要である。)
- The Journal of Chemical Physics 2010, 132, 241102.
- The Journal of Chemical Physics 1991, 94, 751–760.
- In Classical and Quantum Dynamics in Condensed Phase Simulations; WORLD SCIENTIFIC: LERICI, Villa Marigola, 1998; pp 385–404.
- The Journal of Chemical Physics, 2020, 153, 024109.
- The Journal of Chemical Physics, 2022, 144, 214108.
【計算化学】自作pythonモジュールで遷移状態構造を求めてみる(Tris(trifluoromethyl)phosphineが2つ配位したPd触媒によるtrifliorotoleneの還元的脱離, NNP使用）
https://ss0832.github.io/posts/20250811_mop_usage_41/