Gaussianの計算について
- 白金錯体の構造最適化をGaussian03Wで行っています。
- 計算中にSyntax Errorが発生し、lanl2tz(f)が使用できないエラーが発生しました。
- Gaussian計算にお詳しい方のご助力をお願いします。
- ベストアンサー
gaussianの計算について。
gaussianの計算について。 白金錯体の構造最適化をgaussian03Wで行っています。 bp86/genを使用して白金以外を6-31+g(d),白金をlanl2tz(f)として計算をサブミットしても 次のようなエラーで計算がすぐにとまってしまいます。 QPERR --- A SYNTAX ERROR WAS DETECTED IN THE INPUT LINE. lanl2tz(f) ' Last state="Top" TCursr= 746 LCursr= 0 Error reading general basis specification. lanl2tz(f)は使用できないということでしょうか。 ご存知の方はお力添えいただけると助かります。宜しくお願いします。
- chackyb
- お礼率81% (18/22)
- 化学
- 回答数2
- ありがとう数4
- みんなの回答 （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)というキーワードの代わりとなっているということでしょうか。。。
関連するQ&A
- Gaussian98W のエラーメッセージ
Gaussian98W で常磁性アニオン性錯体の励起状態計算を行いましたところ，SCF 収束後に以下のようなエラーメッセージを出して中断してしまいます。 **** Fatal Problem: The smallest beta delta epsilon is -0.15543365D+00 この現象は，ub3lyp/6-31+g* rpa，ub3lyp/3-21g* rpa，ub3lyp/sto-3g rpa，ucis/3-21g*，zindo のどの励起状態計算でも起こります。エラーメッセージの意味や，解決方法などについて，どんな些細なことでも結構ですので，何かご存知の方がいらっしゃいましたら，是非ヒントを下さい。 どうかよろしくお願いいたします。
- ベストアンサー
- 化学
- Gaussianでのバンド計算
量子化学計算ソフト Gaussian では周期構造を設定できるようですが、 無機物質のバンド構造図やフォノン分散図の計算はできますか？ 適用可能な元素、精度、一般のバンド計算やDMol3との比較など、 どのような情報でも歓迎です。宜しくお願い致します。
- ベストアンサー
- 化学
- Gaussian x Gaussian
いつもお世話になっています。 ガウシアン分布の処理について質問させてください。 ある広がりを持った無相関な二つの正規分布があるとして、それぞれの中心地が離れているものの、お互いの裾野同士は重なっています。 例えばその重なりが1%程度の領域同士であった時、 その二つの事象が同時に起こる確率とは、 10^-2 x 10^-2 = 10^-4 ？ それとも SQRT(10^-4 + 10^-4) = 1. 4x10^-2 ? さらに具体的な適用例を示すと、 D-FlipFlopの入力Dataは既知のランダムにばらつく成分を持っているのですが、 入力クロックもまた既知のランダムなばらつき成分を持っています。 クロックの立ち上がりがDataのちょうどど真ん中をたたくとして、 結果出力されるDataのBit Error Rateはどう計算するのだろう。。 というのがこの質問のそもそもの発端です。 よろしくお願いできればと思っております。
- ベストアンサー
- 数学・算数
- gaussianを使った配座エネルギーの計算
gaussian 03 package を使って自由回転によって得られる配座ポテンシャルエネルギーを計算したいのですが方法がわかりません。 よろしくお願いします。 例えばCR4のようにCに側鎖が４つついていて特定のCとRの結合の二面角に準ずる自由回転によって得られる配座ポテンシャルエネルギーを計算したいです。
- ベストアンサー
- 化学
- gaussianで蛍光スペクトルの計算
こんにちは。私立大学4年生です。 gaussianで蛍光波長の求め方が知りたいです。 蛍光波長を求めるインプットファイルを教えて下さい。 下記インプットファイルでやっても、すぐエラーとなって終わります。 分子はナフタレンなどの簡単なものです。 # opt td(Nstates=5,direct) B3LYP/6-31G よろしくお願いします。
- ベストアンサー
- 化学
- Gaussianで計算した分子中の結合次数などを調べたい
いつもお世話になります。 GaussianでH2OやAg(NH3)2などの分子中の結合次数や、 ある分子軌道にどの原子軌道がどれだけ含まれているかというのを調べたいのですが これはどのようなキーワードを使うと調べられるのでしょうか？ Gaussianにはあまり慣れてないので使い方がよく分からなくて…… どうぞよろしくお願いします。
- ベストアンサー
- 化学
- intel fortran cimpilerでプログラムファイルがコンパイルできない
こんにちは。 当方プログラムに関しては全くの初心者です。 研究で早急に計算プログラムを使用したく、 現在intel fortran cimpiler8.0を使い、とあるソースファイルをコンパイルしたいのですがうまくいきません。 コンパイルの手順としてあっているのか自信がないのですが、 とりあえずやってみて、以下のようにエラーが出ます。 $ source /opt/intel_fc_80/bin/ifortvars.csh bash: /opt/intel_fc_80/bin/ifortvars.csh: line 28: syntax error: unexpected end of file $ ifort /home/ディレクトリのパス/ファイル名.f90 ifort: error: could not find directory in which g++ resides これはどういうことなのでしょうか？？ どなたかご教授お願いいたします。 ちなみに、使用しているOSはlinuxのFedora.core1です。
- ベストアンサー
- その他（プログラミング・開発）
- python print文のエラーがでます
#!/usr/bin/env python # -*- coding: utf-8 -*- from __future__ import print_function import datetime class TZ(datetime.tzinfo): def __init__(self, name, offset): self.name = name self.offset = offset def utcoffset(self, dt): return datetime.timedelta(hours=self.offset) def tzname(self, dt): return self.name def dst(self, dt): return datetime.timedelta(0) def main(): path = '/sys/bus/w1/devices/28-000*********/w1_slave' with open(path) as f: data = f.read() temp = float(data[data.index('t=')+2:])/1000 FMT = '{},{:.1f}' JST = TZ('JST', 9) now = datetime.datetime.now(JST) print(FMT.format(now.strftime("%Y/%m/%d %H:%M:%S"), temp)) if __name__ == '__main__': main() とうってじっこうすると print(FMT.format(now.strftime("%Y/%m/%d %H:%M:%S"), temp)) の部分に SyntaxError: invalid syntax とエラーがでます。 どういうことなのか教えていただけるとうれしいです！
- 締切済み
- その他（プログラミング・開発）
- glibc2.4がmakeできない
OSはUbuntu14.04でapt-getのUPDATEなどはしていません。glibcは書籍「自分で作るLinux OS」から取得したものです。INTERNET環境がアナログ電話回線なのでINTERNETからDOWNLOADするというのは現実的ではありません。 はじめはbinutilsのVERSION違いでERRORだったので、 ac_prog_version="$ac_prog_version, ok"; ac_verc_fail=no;; 2.2[0-9]*) を追加しました。 その後makeすると以下のようなERRORが出てglibc2.4がmakeできません。どのようにすればこのERRORを解決することが出来るでしょうか。初心者ですが、よろしくおねがいします。字数制限の関係上、ERROR内容を省略したので、解決に必要な部分があったら補足します。 $ sudo make make -r PARALLELMFLAGS="" CVSOPTS="" -C ../glibc-2.4 objdir=`pwd` all make[1]: ディレクトリ `/usr/local/source/glibc-2.4' に入ります ---中略--- mv -f /usr/local/source/build-glibc/sysd-rulesT /usr/local/source/build-glibc/sysd-rules mawk -f scripts/gen-sorted.awk \ -v subdirs='csu assert ctype locale intl catgets math setjmp signal stdlib stdio-common libio malloc string wcsmbs time dirent grp pwd posix io termios resource misc socket sysvipc gmon gnulib iconv iconvdata wctype manual shadow po argp crypt nss localedata timezone rt conform debug dlfcn elf' \ -v srcpfx='' \ nptl/sysdeps/pthread/Subdirs sysdeps/unix/inet/Subdirs sysdeps/unix/Subdirs assert/Depend intl/Depend catgets/Depend stdlib/Depend stdio-common/Depend libio/Depend malloc/Depend string/Depend wcsmbs/Depend time/Depend posix/Depend iconvdata/Depend nss/Depend localedata/Depend rt/Depend debug/Depend > /usr/local/source/build-glibc/sysd-sorted-tmp mawk: scripts/gen-sorted.awk: line 19: regular expression compile failed (bad class -- [], [^] or [) /[^ mawk: scripts/gen-sorted.awk: line 19: syntax error at or near ] mawk: scripts/gen-sorted.awk: line 19: runaway regular expression /, "", subd ... make[1]: ディレクトリ `/usr/local/source/glibc-2.4' から出ます make[1]: ディレクトリ `/usr/local/source/glibc-2.4' に入ります mawk -f scripts/gen-sorted.awk \ -v subdirs='csu assert ctype locale intl catgets math setjmp signal stdlib stdio-common libio malloc string wcsmbs time dirent grp pwd posix io termios resource misc socket sysvipc gmon gnulib iconv iconvdata wctype manual shadow po argp crypt nss localedata timezone rt conform debug dlfcn elf' \ -v srcpfx='' \ nptl/sysdeps/pthread/Subdirs sysdeps/unix/inet/Subdirs sysdeps/unix/Subdirs assert/Depend intl/Depend catgets/Depend stdlib/Depend stdio-common/Depend libio/Depend malloc/Depend string/Depend wcsmbs/Depend time/Depend posix/Depend iconvdata/Depend nss/Depend localedata/Depend rt/Depend debug/Depend > /usr/local/source/build-glibc/sysd-sorted-tmp mawk: scripts/gen-sorted.awk: line 19: regular expression compile failed (bad class -- [], [^] or [) /[^ mawk: scripts/gen-sorted.awk: line 19: syntax error at or near ] mawk: scripts/gen-sorted.awk: line 19: runaway regular expression /, "", subd ... make[1]: *** `/usr/local/source/build-glibc/abi-versions.h' に必要なターゲット `/usr/local/source/build-glibc/Versions.all' を make するルールがありません. 中止. make[1]: ディレクトリ `/usr/local/source/glibc-2.4' から出ます make: *** [all] エラー 2
- 締切済み
- Linux系OS
お礼
ご回答ありがとうございます。とても的確に教えて下さるのでとても理解しやすいです。 テスト計算のことはどうぞお気になさらずに。見比べながらやっていると良い練習になります☆ 無事に白金錯体のインプットファイルも作成でき、サブミット時にはじかれることもなく順調に計算が開始されました。 あとは収束してくれることを祈ります。 色々と教えて下さって本当にありがとうございました！