量子化学計算でのＳＣＦ収束について
- 量子化学計算において、ＳＣＦが収束しない現象が発生しています。
- 質問者は量子化学計算に興味を持ち、硫酸銅５水和物のエネルギー計算を行いたいと考えています。
- 使用している量子化学計算ソフトはgamessで、計算には次の設定を使用しています。
- ベストアンサー
量子化学計算について質問です！ＳＣＦが収束しない。
どうも始めまして。私は現在量子化学計算に興味を持っていて、硫酸銅５水和物のエネルギー計算をしたいのですが、ＳＣＦが収束しません。量子化学計算ソフトは gamess を使っております。 $CONTRL SCFTYP=ROHF RUNTYP=ENERGY EXETYP=RUN MAXIT=30 ICHARG=2 MULT=2 $END $SYSTEM TIMLIM=525600 MEMORY=1000000 PARALL=.TRUE. $END $BASIS GBASIS=N31 NGAUSS=6 $END $SCF DIRSCF=.TRUE. $END $DATA cuso4*2-5h2o C1 Cu 29.0 0.00032 0.00006 -0.00482 O 8.0 1.87391 0.72893 0.20761 H 1.0 2.30361 0.87851 1.05120 H 1.0 2.45529 0.97384 -0.51367 O 8.0 0.72672 -1.82045 -0.20189 H 1.0 1.65791 -2.03664 -0.17990 H 1.0 0.19403 -2.60517 -0.32474 O 8.0 -1.87905 -0.74390 -0.03441 H 1.0 -2.41031 -0.87040 -0.82288 H 1.0 -2.36790 -1.01005 0.74544 O 8.0 -0.72037 1.83340 0.01802 H 1.0 -1.65133 2.04934 0.00492 H 1.0 -0.18435 2.62459 0.05330 S 16.0 -0.71811 -0.14151 3.54580 O 8.0 0.32389 0.25649 4.53181 O 8.0 -0.10911 -0.10651 2.17980 O 8.0 -1.13011 -1.53351 3.83680 O 8.0 -1.87411 0.76449 3.60880 O 8.0 -0.30103 -0.12963 -4.52492 S 16.0 0.73997 0.26837 -3.53892 O 8.0 0.13197 0.23337 -2.17292 O 8.0 1.15297 1.66137 -3.82992 O 8.0 1.89597 -0.63663 -3.60192 O 8.0 -2.22803 1.96237 -3.57392 H 1.0 -3.16503 1.87236 -3.80092 H 1.0 -1.79303 1.19037 -3.95492 $END これが今回プログラムなのですが、やはり根本的に間違えてるのでしょうか？ それとも初期構造が原因なのでしょうか？ 量子化学計算、化学に詳しい人私にアドバイスをください！ よろしくお願いします！
- kyuyang
- お礼率66% (2/3)
- 化学
- 回答数1
- ありがとう数2
- みんなの回答 （1）
- 専門家の回答
質問者が選んだベストアンサー
■明らかな間違い (1) 電子数が違います。ICHARGは2ではなく-2のはずです。 ■収束させるために (2) 繰り返し回数の上限を大きくします。MAXIT=200 (3) 私の経験からいうとROHFよりUHFの方が収束し易いです。SCFTYP=UHF 上記の三点を変更して計算してみると、収束しました。 FINAL UHF ENERGY IS -3411.7442821010 AFTER 32 ITERATIONS ■より良い計算のために (4) 初期構造をみると、硫酸イオンの酸素と銅イオンの間の距離が短いようです。実験値は2.4オングストロームくらいです。 (5) DFTTYP=PBE0 などのDFT計算も試してみて下さい。常磁性金属イオンを含む系の計算では、ハートリー・フォック法よりも密度汎関数法の方が多用されます。 (6) ECP=SBKJC などの有効内殻ポテンシャルも試してみて下さい(計算するときは、基底関数系を$BASIS GBASIS=SBKJC $ENDのようにECPに合わせる必要があります)。遷移金属錯体の計算をするときに、金属原子の基底関数に6-31Gを使うことは、あまりありません。
関連するQ&A
- He元素のSCF計算
SCF計算について質問があります。 わたしの持っている本（”基礎量子化学” 菊池修著 朝倉書店）にHe原子のSCF計算が具体例としてのっているのですが，1sのスレーター関数χ_iに入れるζ_iの値が ζ_1 = 1.45，ζ_2 = 2.90 となっています。 この数値はどこから出てきたんでしょうか？ ζ=有効核電荷/主量子数だとすると1.70になる気がするのですが…。 わかりにくい文章で済みません。何か補足すべきことがあれば言ってください。
- ベストアンサー
- 化学
- 量子化学計算_対称な分子構造を作るには
量子化学計算用の入力ファイルを、左右対称な形でつくる方法を教えて下さい。 ガウシアン DFT B3LYP/6-31G(d) を用いて、量子化学計算をしたいと思っています。有機分子の構造最適化を行い、HOMO、LUMOエネルギーレベル、分子軌道分布を調べることが目的です。分子量1000程度の、C, H, O, S, B, Fを含む分子を扱っています。 自分の所属する研究室では、chem3Dを使って粗く最適化した入力ファイルの作製をしています。 chemdrawで構造を作って、chem3Dに移します。まずMM2で最適化、さらにMOPACで最適化します。 扱っている分子の構造式は左右対称な形なのですが、chem3Dで最適化して出来上がった入力ファイルの立体構造が左右不対称になっています。 そこで質問です。 1.非対称ということは最適化がうまくできていないと考えていいでしょうか？ 2.計算条件で、対称になるように条件を制限することはできるのでしょうか？できるならその方法（chem3Dなどのソフトの操作方法）を教えて頂けないでしょうか？ 3.実際にどのような形をしているかは計算では分からないものなのでしょうか？入力ファイルの構造をもとに、少しだけねじったり伸ばしたりして、エネルギーが最小化される点を探索しているイメージでしょうか？ なお、分子のX線構造解析のデータは持っていません。試料はあるのですが単結晶が得られていないので。 量子化学のことを殆ど知らないものですが、ご回答よろしくお願いします。
- ベストアンサー
- 化学
- 物理化学実験についてです
物理化学実験の予習で ヘスの法則に従い、硫酸銅五水和物と無水硫酸銅の 生成エンタルピーから溶解熱を計算せよ というものがあったのですが、どうやって計算するのかわかりません。 硫酸銅五水和物の生成エンタルピーは-2279.7で無水硫酸銅のせエンタルピーは-771.36です。 溶解の式CuSO4・5H2O (s) + aq → CuSO4・5H2O (aq) (この式も全然自信ないですが)を使うのですか？でも硫酸銅五水和物の生成エンタルピーしかわかってないし、どのように求めるのかわかりません。お願いします。
- ベストアンサー
- 化学
- 「量子」と「量子化されている」について
古いタイプの有機化学者の質問ですが宜しくお願いいたします。 軌道を説明するとき、s軌道、p軌道、、とエネルギーが不連続に上がっていきます。その時に「量子化されている」との説明がされています。 専門家の皆さんに、「量子」と「量子化されている」について、分かり易く（数式抜きで）ご教示頂きたく宜しくお願いいたします。
- ベストアンサー
- 物理学
- センター化学 過去問 解説お願いします
いつも大変お世話になっております。 化学について質問させてください。 「硫酸銅（II）五水和物62mgを空気中で900°Cで加熱したところ、質量が20mgになった。 加熱して得られた物質の組成として最も適当なものをひとつ選べ。 １．Cu ２．ＣuO ３．Ｃu2O ４．ＣuＳO4 ５．ＣuＳO4・H2O」 に対して、250（硫酸銅（II）五水和物の分子量）×20÷62＝80 その80の分子量に相当する２が答えとする解説があるのですが、 このような単純な比の計算で答えに至る理由がどうしてもわかりません。 よろしければ詳しく解説いただけないでしょうか。よろしくお願いします。
- ベストアンサー
- 化学
- 量子化学 HFの軌道
量子化学ですが、 Hは1ｓに1電子。 Fは1ｓに2で、2ｓに2。2ｐｘ、2ｐｙに2つずつ2ｐｚに1つ。 HとFの結合は、Hの1ｓとFの2ｐｚの軌道で起こる。 ここで、結合性分子軌道σのモデルみたいなのは ●○ （●が少し○より大きい。●は+、○は-） で σに2つ電子が入っているので理解できます。 エネルギー的にその上の非結合性結合にFは2ｐｘ、2ｐｙに2つずつ電子使い。 質問はここからで、 Hの1ｓとFの2ｐｚは結合性分子軌道で1つずつつかって、エネルギー的に非結合性結合の上の反結合性分子軌道σ*に、使える電子がないのに、●○●というモデルが、本に載せられていました。電子がないのに、この反結合性分子軌道は存在するのでしょうか？それとも、自分の考え方が違うのでしょうか？ どなたか教えてください。お願いします。
- ベストアンサー
- 化学
- GAMESSの出力ファイルの内容について
はじめまして。 私は応用物理専攻の院生です。 現在量子化学計算プログラム「GAMESS」を使って、メタンやエチレンなど比較的小さな有機分子の価電子の束縛エネルギーと波動関数を求めようとしています。 基底関数を6-31+G(d)として計算させたときの出力ファイルの内容について質問します。 以下が出力ファイルの一部です。 ATOMIC BASIS SET ---------------- THE CONTRACTED PRIMITIVE FUNCTIONS HAVE BEEN UNNORMALIZED THE CONTRACTED BASIS FUNCTIONS ARE NOW NORMALIZED TO UNITY SHELL TYPE PRIM EXPONENT CONTRACTION COEFFICIENTS C 1___S___1___3047.524880______0.536345______( 0.001835) 1___S___2_____457.369518______0.989452______( 0.014037) 1___S___3_____103.948685______1.597283______( 0.068843) 1___S___4_______29.210155______2.079187______( 0.232184) 1___S___5_________9.286663______1.774174______( 0.467941) 1___S___6_________3.163927______0.612580______( 0.362312) 炭素原子の基底関数6-31+G(d)の1s(n=1)軌道の変数です。 見やすくするために_(アンダーバー)を挿入しています。 ４列目はexponents係数、６列目はexpantion係数だということがわかりました。 これらの値は今勉強している参考書[1]と同じ値でした。 ５行目の値が何を示しているのかがわかりません。 GAMESSのマニュアルや他の参考書、HPなど見ましたがわかりませんでした。 量子化学計算関係の論文は私にはまだ難解でよくわかりませんでした。 量子化学では常識的なことかもしれませんが、宜しくお願い致します。 [1] Physical Chemistry - a molecular approach(原著), A.Donald et al., Univercity science Book, p424
- ベストアンサー
- 化学
- gamessのアウトプットについて質問です。
私は最近量子化学計算に興味を持ち、gamessとゆうプログラムソフトに手を出したのですが、少々わからないことがあるので質問します。 今、水分子(H2O)を、HF/6-31Gでエネルギー計算をして、アウトプットを見ていて疑問に思ったのですが、TOTAL ENERGYとTOTAL POTENTIAL ENERGYの違いはなんですか？ あと、水分子のTOTAL ENERGY = -76hartreeと出たので、単位をkcal/molにするため627.51掛けたのですが、とても大きな値が出て焦っています。 なぜ627.51も掛けるのか？ あと、単位 hartreeがなにを表しているのか教えて下さい。 よろしくお願いします。
- 締切済み
- 化学
お礼
本当にありがとうございます！ たくさん勉強になりましたし、 量子化学計算についてもっと勉強したいと思いました！