※ 引述《allenc38 (我現在只想被肛)》之銘言： : 用 #pw91pw91/lanl2dz opt 計算含五個Au的分子, 基本上如果你預設用g09沒帶電荷, 5個Au的cluster預設是doublet gaussian程式預設會直接用open-shell去算, 也就是用upw91pw91, 已經是開殼層計算了. : 得到優化的等腰梯型結構,再以 : # ropw91pw91/lanl2dz freq pop=full 計算 請問你是要算open-shell singlet state嗎? 如果是的話才可能會考慮到使用RO 不然如果單純就你的問題, 只想使用開殼層計算的話 就是像上面說的, 使用upw91pw91/lanl2dz即可 : 結果是 error L502 : 加了use=L506 依然失敗 : 有解決計算不了開殼系統的方法嗎? : (編輯:用#RO3LYP/lanl2dz 可以成功計算, 但我需要pw91pw91...) L502一般是做SCF不收斂, 解決方法有很多 常見的像是提高SCF的cycle數(不過真正無法收斂的系統怎麼加都是沒用的) 或者是提高計算精度, 例如SCF=XQC 等等方法 不過你的問題可能需要先釐清是否需要計算open-shell singlet state 如果不用的話我想直接使用upw91pw91應該是可以的 BTW: 不知道你是否有打錯, RO3LYP是利用O3LYP泛函做閉殼層計算喔 除非你是說ROB3LYP 這個才是用B3LYP做RO的計算這樣 --

※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 59.125.27.109 (臺灣) ※ 文章網址: https://webptt.com/m.aspx?n=bbs/Chemistry/M.1587557787.A.934.html