GAMESSの概要
GAMESS(US)は、Iowa州立大のGordon，Moscow州立大のGranovsky，大阪府立大の小関先生らによって開発及び配布が行われている、無償提供の非経験的分子軌道法/密度汎関数理論計算プログラムです。無償で提供されているにも関わらず、詳細なマニュアルが配布され、またありとあらゆる OS，プラットフォームへ移植されています(Windows,Mac,Unix,Linux,DOS,OS/2,Solaris...)。
GAMESS(US)は、Iowa州立大のGordon，Moscow州立大のGranovsky，大阪府立大の小関先生らによって開発及び配布が行われている、無償提供の非経験的分子軌道法/密度汎関数理論計算プログラムです。無償で提供されているにも関わらず、詳細なマニュアルが配布され、またありとあらゆる OS，プラットフォームへ移植されています(Windows,Mac,Unix,Linux,DOS,OS/2,Solaris...)。
GAMESSの入力ファイルは、その内容の全てが「$XXXXX～$END」というブロック(グループという)に分けて記述されています。各グループ内に「keyword=property」の形で具体的な計算指示が格納されます。プロパティの多くにはデフォルト値が設定されていて入力が不要なケースも多いですが、ユーザー側で正確に指定する必要のある項目も多数存在します。
GAMESSのキーワードはMOPAC以上に膨大な数に上り、とても全てを紹介することはできません (当然筆者は全てのキーワードを知っているわけではありません)。ここでは、当サイトで使用頻度の高いキーワードに絞り、その意味と利用方法を説明していきたいと思います。
GAMESSには、基本的な基底関数はプログラムに含まれていて、特定のキーワードを使って呼び出して使います。$BASISグループが、基底関数に関する指定を行うセクションになります。既に「キーワードの説明」で、よく使うものについて説明していますが、ここではさらに詳しい入力について説明を加えます。
PC GAMESSはWindows環境でも簡単に並列化することができます。最近はデュアルコアCPU(人によってはクアッドコアも？)を利用している方も多いと思いますので、SMP環境に適応した簡単な並列化の方法をお示しします。
非常に大きな分子を計算するとき、その計算時間を削減する方法としてQM/MM法が現在有力な方法の一つとなっています。大きな分子においては全ての原子が重要ではなく、多くの原子は立体的な影響のみで、電子状態が重要な寄与を果たすのは一部であるという考え方は、直感的に受け入れやすいものです。そういう場合、分子を「精密に計算する部分:高精度ab initio法やDFT」と「粗く計算する部分:低精度ab initio法や半経験法，分子力学」という風に区分して計算できれば、計算資源の節約につながります。
GaussianではONIOM法が実装されており、分子を最大3つのレイヤーに切り分け、各層に別々のモデル化学(非経験的分子軌道法および密度汎関数理論，半経験的分子軌道法，分子力学)を割り当てることができます。
一方、GAMESSではONIOM法の基となったIMOMM/SIMOMM法が実装されています。これは、分子を2つのレイヤーに切り分け、一方に量子力学計算を、もう一方に分子力学計算を割り当てます。分子力学のエンジンとしてはTINKERが使われています。ここでは、IMOMM法を使ったQM/MM計算の簡単な例を紹介します。
分子間の相互作用は、計算化学が扱う主要な問題の一つです。単純に考えると、例えばAという分子とBという分子が相互作用してABという錯体を形成する場合、A,Bが単独で存在したときのエネルギーと錯体ABのエネルギーを比較して、その相互作用がどれぐらい強いかを議論するということになります。しかし、ここには思いがけない落とし穴が存在します。それが本稿で紹介するBSSE(Basis Set Superposition Error:基底関数重なり誤差)です。以下、GAMESSでのBSSEの計算方法とモデル化学による差について見てみたいと思います。
量子化学計算により出力されるデータの中で、必ずと言っていいほど利用されるのがエネルギー値です。その値と出力パターンは計算の方法(プログラム・計算目的・モデル化学)の違いで変わってきます。
ここでは、GAMESS(WinGAMESS/PC GAMESS)の出力ファイルのどこにエネルギーが記載され、その値はどのように処理すべきかについて、ごく簡単に述べたいと思います。
PC GAMESSの並列化については以前に書きましたが、ここではWinGAMESSの簡潔な並列化について紹介します。PC GAMESSとWinGAMESSはどちらも「GAMESS」とつきますが、その中身はかなり違っていて、例えばPC GAMESSではMP3/4計算やRSURFACE計算，Cube出力の実装，そして何より計算が高速なのが特長で、WinGAMESSではCCやNEO，EFP/PCM，NMR，TinkerによるQM/MM，UMP2 GradientなどPC GAMESSより幅広い計算が可能になっています。
分子には基底状態と励起状態があります。多くの有機反応は基底状態で起こりますが、光照射下では励起状態からも反応が起こります。また、紫外・可視スペクトルは共役系の状態を知るために有用な情報を与えますが、これは電子の励起エネルギーを観測していることに他なりません。当然、物質の色もこの電子励起に伴う光吸収に起因しています。
励起状態は、分子軌道計算でもちろん扱うことができます。ですが、励起状態は基底状態と異なり無数にあるので、その取り扱いは時に非常に難しくなります。ここでは、最もよく励起状態が調べられている分子の一つ・ホルムアルデヒドを例に、励起エネルギーの計算と励起状態の構造最適化について概観します。
(C) 2002-2006 s2k (W. Sasaki) All Rights Reserved.