- ベストアンサー
Gaussianの計算について
- 白金錯体の構造最適化をGaussian03Wで行っています。
- 計算中にSyntax Errorが発生し、lanl2tz(f)が使用できないエラーが発生しました。
- Gaussian計算にお詳しい方のご助力をお願いします。
- みんなの回答 （2）
- 専門家の回答
質問者が選んだベストアンサー
> lanl2tz(f)というキーワードをインプットファイルに記すことなく > 計算が行われているのは、「Au o」から続く数字をインプットファイルに示すことで > lanl2tz(f)というキーワードの代わりとなっているということでしょうか。。。 はい。おおまかには、そんな感じです。 詳しく言うと、むしろ話は逆で、lanl2tz(f)という名前の基底関数がもしGaussianの中に入っていたならば、「lanl2tz(f)というキーワードをインプットファイルに示すことで、『Au 0』から続く数字の代わりになる」ということです。 「Au 0」から続く数字は、計算に使われるAuの基底関数そのものと考えてください。基底関数の開発者がこの基底関数につけた名前が、lanl2tz(f)になります。lanl2tz(f)が発表されたのはGaussian03のリリース後なので、インプットファイルにlanl2tz(f)と書いても、gaussian03はそれが基底関数の名前ということを知りませんから、エラーになります。ですので、基底関数の名前でなく、基底関数そのものをインプットファイルに示す必要があります。 あと、いま気づいたのですけど、テスト計算はAuでなくてPtにしておけばよかったですね。白金錯体を金錯体と勘違いしてました。すみません。でも#1のテスト計算がうまく走ったのなら、あとはPtの基底関数を入手するだけで白金錯体の計算ができると思います。
その他の回答 (1)
- 101325
- ベストアンサー率80% (495/617)
> lanl2tz(f)は使用できないということでしょうか。 いいえ。使用できます。注意点が二つあります。 (1) lanl2tz(f)はgaussian03に入っていません。まず EMSL Basis Set Exchange などから基底関数を入手する必要があります。 https://bse.pnl.gov/bse/portal フォーマットはGaussian94でいいです。 (2) lanl2tz(f)などの有効内殻ポテンシャル(ECP, 擬ポテンシャル)を使う基底関数をgenで指定したときには、キーワード pseudo=read が必要です。 テスト計算用の入力ファイル -----ここから----- # bp86/gen pseudo=read opt （空行） AuH （空行） 0 1 Au H 1 1.524 （空行） H 0 6-31+g(d) **** Au 0 S 1 1.00 2.8090000 1.0000000 （２４行省略） F 1 1.00 1.0500000 1.0000000 **** （空行;ここの空行が２行あるとうまくいかない） AU 0 AU-ECP 4 60 g potential 5 1 622.6287956 -60.0000000 （３１行省略） 2 2.8965118 15.3424188 （空行） （空行） -----ここまで----- 計算して結合距離が1.5456オングストロームになれば大丈夫だと思います。
補足
とても丁寧に解説してくださって本当にありがとうございます。 早速AuHのテスト計算を行ってみました。無事に結合距離1.5456?の構造を得ることができました！ 計算初心者でもトライ出来るように回答してくださってとても感謝しています。ありがとうございます☆ 一つ、お伺いしたいのですが、、、lanl2tz(f)というキーワードをインプットファイルに記すことなく 計算が行われているのは、「Au o」から続く数字をインプットファイルに示すことで lanl2tz(f)というキーワードの代わりとなっているということでしょうか。。。
お礼
ご回答ありがとうございます。とても的確に教えて下さるのでとても理解しやすいです。 テスト計算のことはどうぞお気になさらずに。見比べながらやっていると良い練習になります☆ 無事に白金錯体のインプットファイルも作成でき、サブミット時にはじかれることもなく順調に計算が開始されました。 あとは収束してくれることを祈ります。 色々と教えて下さって本当にありがとうございました！