請(qǐng)問(wèn)G4方法能不能算的動(dòng)三四十個(gè)原子的分子
請(qǐng)問(wèn)G4方法能算的動(dòng)三四十個(gè)原子的分子結(jié)構(gòu)嗎。我看別人的文獻(xiàn)中用G4算乙基咔唑的能量，可是我跑了15個(gè)小時(shí)了，還沒(méi)有算完，目前的狀態(tài)是opt結(jié)束了，后面就一直卡著沒(méi)結(jié)果，實(shí)際上只算了7個(gè)小時(shí)。用的服務(wù)器，36核算的。我的輸入文件如下：
# G4 geom=connectivity
necz
0 1
C -1.14579300 0.49395000 -0.00001300
C -0.96842200 -0.91454000 0.00001300
C 0.45961100 -1.15076200 0.00000700
C 1.10413000 0.12488100 -0.00002000
N 0.10752000 1.12339000 -0.00002400
C 1.21791600 -2.32663400 0.00002200
C 2.60766900 -2.24443400 0.00000500
C 3.23702200 -0.98827700 -0.00003000
C 2.50542400 0.19988000 -0.00004100
C -2.42249000 1.06833400 -0.00002600
C -3.52355000 0.21048700 -0.00000800
C -3.36228800 -1.18688000 0.00002100
C -2.08817400 -1.75394600 0.00003100
C 0.18141300 2.59551500 -0.00007000
C 1.57214900 3.22560800 0.00011600
H 0.72562900 -3.29068700 0.00004600
H 3.20564100 -3.14553000 0.00001300
H 4.31787600 -0.93431500 -0.00004600
H 3.03074100 1.13873800 -0.00007000
H -2.56474100 2.13969000 -0.00005200
H -4.52066400 0.63070900 -0.00001900
H -4.23598900 -1.82446300 0.00003300
H -1.96889900 -2.82973300 0.00005100
H -0.36611100 2.95292500 -0.87792300
H -0.36635300 2.95297800 0.87761000
H 1.45424700 4.31072200 0.00011700
H 2.14403100 2.95805300 0.88734700
H 2.14426000 2.95809000 -0.88698100
1 2 1.5 5 1.0 10 1.5
2 3 1.5 13 1.5
3 4 1.5 6 1.5
4 5 1.0 9 1.5
5 14 1.0
6 7 1.5 16 1.0
7 8 1.5 17 1.0
8 9 1.5 18 1.0
9 19 1.0
10 11 1.5 20 1.0
11 12 1.5 21 1.0
12 13 1.5 22 1.0
13 23 1.0
14 15 1.0 24 1.0 25 1.0
15 26 1.0 27 1.0 28 1.0
16
17
18
19
20
21
22
23
24
25
26
27
28
返回小木蟲查看更多
京公網(wǎng)安備 11010802022153號(hào)
G4 在熱力學(xué)組合方法里都是高端的存在。20多個(gè)原子左右基本已經(jīng)是差不多的上限了 ，里面涉及到了MP4 以及 CCSD(T)的步驟。快是肯定不快的。慢慢等著就行。一般30核心左右的，15,6個(gè)原子左右的 可能跑上一個(gè)小時(shí)差不多。原子數(shù)目翻一倍那計(jì)算時(shí)間翻到50， 60倍差不多。40個(gè)原子 恐怕就很難算的動(dòng)了。
至于cpu使用率低是因?yàn)?post-HF方法 需要大量的讀取 寫入硬盤。cpu在那里就是干等著的 硬盤越快 cpu的效率就越高。
原來(lái)如此，難怪卡在昨晚凌晨四點(diǎn)就再也不動(dòng)了，非常感謝，請(qǐng)問(wèn)前輩，像我這種大分子，有沒(méi)有推薦的熱力學(xué)組合方法，G3方法在現(xiàn)在來(lái)看是不是有些過(guò)時(shí)了？ps. 我暫時(shí)用的36核，學(xué)校服務(wù)器一般一個(gè)人最多調(diào)用72核，不然會(huì)被別的人噴
，
G3確實(shí)過(guò)時(shí)了。高斯內(nèi)置的方法里 G4MP2 和 CBS-QB3 性價(jià)比比較高。比G4快不少。 不過(guò)即使是這兩種方法 30個(gè)原子也基本是上限了。
好的好的，謝謝前輩，我用G4MP2和CBS-QB3提交兩個(gè)任務(wù)試試
如果還算不動(dòng)，那就試試 G3MP2B3 這個(gè)還能再快些 但是精度也低些。CBS-4M 基本是高斯支持的方法里最快的�？梢运阍訑�(shù)能多一點(diǎn)點(diǎn) 但是精度比前幾個(gè)可要低不少了。
好的好的，非常感謝前輩，那我就再提交兩個(gè)任務(wù)，我現(xiàn)在就是在篩選一個(gè)合適的方法