β版はβ版だからこそβ版と呼ばれる
何だこのタイトル。
ORCA 2.7.0 betaを使っていて、並列MP2計算がおかしいことにやっと気づいた話。
NORIでのみ、並列計算したときのエネルギー値がずれる。RI-MP2ならずれない。2.6.35でどうだったかを検証するのが面倒なので(OpenMPI版無いし)、実は前からだったのかは分からないが…
具体的には、1,1'4',1'',4'',1'''-Quaterphenyl(ベンゼン環がパラ位で4つつながった分子)の一点計算をしていて気がついた。RI-MP2(PAL1～4)とMP2(PAL1)ではE(MP2corr)=3.11 Ehになるところ、MP2(PAL2～4)ではE(MP2corr)=2.71 Ehになるという罠。ちなみにbiphenylぐらいのサイズの分子なら大丈夫(そういうところが気づきにくい原因かなぁ？)。
ま、β版なんで、正式版で修正されている類の話でしょう。