| 查看: 1406 | 回復(fù): 10 | ||
w824194485木蟲 (正式寫手)
|
[求助]
計算出錯 已有1人參與
|
|
|
輸入文件 %NOSAVE %mem=500MB # OPT=gdiis freq CBS-QB3 Title Card Required 0 2 N S 1 B1 O 2 B2 1 A1 B1 1.50866118 B2 1.45832917 A1 124.00909224 輸出文件 ----------------------------------------------------------- #N Geom=AllCheck Guess=TCheck SCRF=Check CCSD(T)/6-31+G(d') ----------------------------------------------------------- 1/18=20,19=11,29=7,38=1,40=1/1; 2/12=2,40=1/2; 3/5=11,6=6,7=11,11=9,14=-4,16=1,25=1,30=1,70=2,116=-2/1,2,3; 4/5=101/1; 5/5=2/2; 8/6=4,9=120000,10=1/1,4; 9/5=7,14=2/13; 6/7=2,8=2,9=2,10=2/1; 99/5=1,9=1/99; ------------------- Title Card Required ------------------- Redundant internal coordinates taken from checkpoint file: /data/tmp/Gau-29637.chk Charge = 0 Multiplicity = 2 N,0,-0.0063063572,0.,-0.0065648995 S,0,0.003093632,0.,1.5051797991 O,0,1.2120929758,0.,2.3343866301 Recover connectivity data from disk. Input orientation: --------------------------------------------------------------------- Center Atomic Atomic Coordinates (Angstroms) Number Number Type X Y Z --------------------------------------------------------------------- 1 7 0 -0.006306 0.000000 -0.006565 2 16 0 0.003094 0.000000 1.505180 3 8 0 1.212093 0.000000 2.334387 --------------------------------------------------------------------- Distance matrix (angstroms): 1 2 3 1 N 0.000000 2 S 1.511774 0.000000 3 O 2.639044 1.466037 0.000000 Stoichiometry NOS(2) Framework group CS[SG(NOS)] Deg. of freedom 3 Full point group CS NOp 2 Largest Abelian subgroup CS NOp 2 Largest concise Abelian subgroup C1 NOp 1 Standard orientation: --------------------------------------------------------------------- Center Atomic Atomic Coordinates (Angstroms) Number Number Type X Y Z --------------------------------------------------------------------- 1 7 0 1.401806 -0.231005 0.000000 2 16 0 0.000000 0.335034 0.000000 3 8 0 -1.226580 -0.467938 0.000000 --------------------------------------------------------------------- Rotational constants (GHZ): 68.7554633 9.7059073 8.5052574 Standard basis: 6-31+(d') (6D, 7F) There are 45 symmetry adapted basis functions of A' symmetry. There are 16 symmetry adapted basis functions of A" symmetry. Integral buffers will be 131072 words long. Raffenetti 2 integral format. Two-electron integral symmetry is turned on. 61 basis functions, 120 primitive gaussians, 61 cartesian basis functions 16 alpha electrons 15 beta electrons nuclear repulsion energy 96.6357984684 Hartrees. NAtoms= 3 NActive= 3 NUniq= 3 SFac= 7.50D-01 NAtFMM= 80 NAOKFM=F Big=F One-electron integrals computed using PRISM. NBasis= 61 RedAO= T NBF= 45 16 NBsUse= 61 1.00D-06 NBFU= 45 16 Initial guess read from the checkpoint file: /data/tmp/Gau-29637.chk B after Tr= 0.000000 0.000000 0.000000 Rot= 1.000000 0.000000 0.000000 0.000000 Ang= 0.00 deg. Initial guess orbital symmetries: Alpha Orbitals: Occupied (A') (A') (A') (A') (A') (A") (A') (A') (A') (A') (A') (A") (A') (A') (A") (A') Virtual (A") (A') (A') (A") (A') (A') (A') (A") (A') (A") (A') (A') (A") (A') (A') (A') (A") (A') (A') (A') (A") (A') (A') (A') (A") (A') (A') (A") (A") (A') (A") (A') (A') (A") (A') (A") (A') (A') (A') (A') (A') (A') (A') (A') (A') Beta Orbitals: Occupied (A') (A') (A') (A') (A') (A") (A') (A') (A') (A') (A") (A') (A') (A") (A') Virtual (A') (A") (A') (A') (A") (A') (A') (A') (A") (A') (A") (A') (A') (A") (A') (A') (A") (A') (A') (A') (A') (A") (A') (A') (A') (A") (A') (A") (A') (A") (A') (A") (A') (A') (A") (A') (A") (A') (A') (A') (A') (A') (A') (A') (A') (A') Initial guess <Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 0.5000 <S**2>= 0.7644 S= 0.5072 Requested convergence on RMS density matrix=1.00D-08 within 128 cycles. Requested convergence on MAX density matrix=1.00D-06. Requested convergence on energy=1.00D-06. No special actions if energy rises. Keep R1 and R2 ints in memory in canonical form, NReq=4427933. Rare condition: small coef for last iteration: 0.000D+00 Rare condition: small coef for last iteration: 0.000D+00 Rare condition: small coef for last iteration: 0.000D+00 Rare condition: small coef for last iteration: 0.000D+00 Rare condition: small coef for last iteration: 0.000D+00 Rare condition: small coef for last iteration: 0.000D+00 Rare condition: small coef for last iteration: 0.000D+00 Rare condition: small coef for last iteration: 0.000D+00 Rare condition: small coef for last iteration: 0.000D+00 Rare condition: small coef for last iteration: 0.000D+00 Rare condition: small coef for last iteration: 0.000D+00 >>>>>>>>>> Convergence criterion not met. SCF Done: E(UHF) = -526.685148166 A.U. after 129 cycles Convg = 0.3454D-05 -V/T = 2.0011 <Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 0.5000 <S**2>= 0.9197 S= 0.5815 <L.S>= 0.000000000000E+00 Annihilation of the first spin contaminant: S**2 before annihilation 0.9197, after 0.7571 Convergence failure -- run terminated. Error termination via Lnk1e in /opt/g09/l502.exe at Tue Oct 6 13:37:29 2015. |

木蟲 (著名寫手)

木蟲 (著名寫手)
|
你看看帖子就知道了,就是你的計算沒有收斂, 這不是你采用了opt=gdiis的方式來調(diào)節(jié)收斂了, 不收斂的問題分為兩個:幾何優(yōu)化不收斂和scf不收斂。 我的建議是:(個人見解,不一定對哈) 1:opt=(maxstep=10,notrust) 這個是因為默認(rèn)補償是0.03bohr,我們調(diào)節(jié)為0.01,減小步長,看看能否成功。 2:可以調(diào)節(jié)對稱性: symm=loose 3:我曾經(jīng)成功的方式:先用小的機組,再用大的機組。 4: freq和opt是否可以分開計算? 5:使用maxcyc來調(diào)節(jié)步數(shù),基本不管用。 以上是個人見解,不一定對,關(guān)于兩個收斂的調(diào)節(jié)方式有很多帖子,樓主可以多看看帖子,都是從不會,到慢慢熟悉的,我的感覺是這樣哈 |

木蟲 (正式寫手)

木蟲 (著名寫手)
論文?苦!主。
木蟲 (正式寫手)

木蟲 (正式寫手)

木蟲 (正式寫手)

木蟲 (正式寫手)

木蟲 (著名寫手)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +6 | tcxiaoxx 2026-03-20 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-18 | 6/300 |
|
|
[考研] 280求調(diào)劑 +11 | 咕嚕曉曉 2026-03-18 | 12/600 |
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
[考研] 材料工程(專)一志愿985 初試335求調(diào)劑 +3 | hiloiy 2026-03-17 | 4/200 |
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
[考研] 初始318分求調(diào)劑(有工作經(jīng)驗) +3 | 1911236844 2026-03-17 | 3/150 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩,總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
[考研] 求調(diào)劑,一志愿:南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
[考研] 081700化工學(xué)碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
[考研]
|
不想起名字112 2026-03-19 | 3/150 |
|
|
[考研] 生物學(xué)調(diào)劑招人。。 +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[考研] 一志愿蘇州大學(xué)材料工程(085601)專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|