| 24小時(shí)熱門版塊排行榜 |
| 9 | 1/1 | 返回列表 |
| 查看: 4517 | 回復(fù): 8 | |||
shenzhendsh鐵蟲(chóng) (小有名氣)
|
[求助]
Gaussian TDDFT 激發(fā)態(tài) 頻率計(jì)算 輸入文件
|
|
我要計(jì)算單重激發(fā)態(tài)的頻率,輸入文件如下: %chk=chk.chk %nprocshared=8 %mem=3200MB # freq td=(nstates=30,root=1,singlet) b3lyp tzvp title Cartesians coordination 請(qǐng)問(wèn)為何一直算不完,而且有很多SCF Done能量?正確的輸入命令是什么呢? 對(duì)于激發(fā)態(tài)過(guò)渡態(tài)的優(yōu)化,輸入命令又是什么呢? 剛接觸激發(fā)態(tài)計(jì)算,還望大家不吝賜教!非常感謝! |
木蟲(chóng) (著名寫手)
|
一直算不完是什么意思?算了一天??jī)商?而且你算?0個(gè)激發(fā)態(tài),如果體系大的話算個(gè)十幾天都很正常。 SCF的能量不用管它,都會(huì)算的,只需關(guān)注最后一次循環(huán)的輸出結(jié)果即可。 如果要做激發(fā)態(tài)過(guò)渡態(tài)優(yōu)化,加opt關(guān)鍵詞就可以了,但注意如此做的話最好用--Link1--把opt與freq分成兩個(gè)任務(wù)來(lái)算,不然很有可能會(huì)報(bào)錯(cuò)。 |
鐵蟲(chóng) (小有名氣)
|
多謝回復(fù)!體系含三個(gè)V,8個(gè)O,一個(gè)C,四個(gè)H,8個(gè)CPU跑了兩天了,而且SCF Done的能量幾乎不變。 做激發(fā)態(tài)過(guò)渡態(tài)優(yōu)化時(shí),只加opt不行吧,我加了opt(ts,calcfc),結(jié)果報(bào)錯(cuò)如下:No analytic 2nd derivatives for this method. >Error termination via Lnk1e in /home/soft/g09/l1.exe at Tue Nov 12 11:30:43 2013. 優(yōu)化時(shí)應(yīng)加什么關(guān)鍵詞呢?多謝回復(fù)! |
木蟲(chóng) (著名寫手)
鐵桿木蟲(chóng) (正式寫手)
| Gaussian的TDDFT算激發(fā)態(tài)頻率是數(shù)值的,所以需要算6N+1次(N是體系的原子數(shù)目)TD,所以你會(huì)找到6N+1個(gè)"SCF Done",什么時(shí)候"SCF Done"的數(shù)目到達(dá)這個(gè)了,什么時(shí)候你的激發(fā)態(tài)頻率計(jì)算就結(jié)束了。2天時(shí)間太少,TDDFT算個(gè)激發(fā)態(tài)頻率20天到一個(gè)月都是正常的,慢慢等吧。沒(méi)事數(shù)數(shù)"SCF Done"的數(shù)目,大概估算一下需要算多久吧。 |
鐵蟲(chóng) (小有名氣)
木蟲(chóng) (正式寫手)
木蟲(chóng) (著名寫手)
|
大神你好,我在這問(wèn)您個(gè)問(wèn)題,就是我有個(gè)分子,大約含有110額原子,之前優(yōu)化激發(fā)態(tài)一直出錯(cuò),使用的方法是b3lyp/6-31g*。假如說(shuō),我計(jì)算激發(fā)態(tài)的頻率,能使用b3lyp/6-31g*嗎?sob給的解說(shuō)是這樣的:Gaussian沒(méi)有TDDFT解析二階導(dǎo)數(shù),這是幾乎做不動(dòng)的,耗時(shí)我估計(jì)是600~1200次激發(fā)能計(jì)算的耗時(shí)。必須要計(jì)算頻率的話,只能是在CIS下優(yōu)化并計(jì)算頻率,因?yàn)镃IS下有解析二階導(dǎo)數(shù),速度和TDDFT計(jì)算頻率耗時(shí)完全不在一個(gè)數(shù)量級(jí)。 我想,基態(tài)是b3lyp/6-31g*,激發(fā)態(tài)用CIS能行嗎?b3lyp真的做不了我的這個(gè)計(jì)算?我曾經(jīng)計(jì)算過(guò)一個(gè)70個(gè)分子,很長(zhǎng)很長(zhǎng)時(shí)間,我該怎么辦?謝謝pwzhou的指點(diǎn)。謝謝您 |

鐵蟲(chóng) (初入文壇)
|
你好,請(qǐng)問(wèn)激發(fā)態(tài)的頻率計(jì)算怎么做?為什么我計(jì)算老出錯(cuò) 分子的激發(fā)態(tài)優(yōu)化后做頻率計(jì)算 輸入文件是 %chk=b.chk %rwf=b.rwf %mem=200MW %nprocshare=12 # freq td=(singlets,nstates=6,root=1) b3lyp/6-31g(d) 計(jì)算到下邊就停止了 SavETr: write IOETrn= 770 NScale= 10 NData= 16 NLR=1 LETran= 118. IDoAtm=1111111111111111111111111111111111111111 Differentiating once with respect to electric field. with respect to dipole field. Differentiating once with respect to nuclear coordinates. There are 1 degrees of freedom in the 1st order CPHF. IDoFFX=0. LinEq1: Iter= 0 NonCon= 1 RMS=2.02D-03 Max=2.21D-01 AX will form 1 AO Fock derivatives at one time. LinEq1: Iter= 1 NonCon= 1 RMS=4.31D-04 Max=2.25D-02 LinEq1: Iter= 2 NonCon= 1 RMS=2.91D-04 Max=1.36D-02 LinEq1: Iter= 3 NonCon= 1 RMS=2.27D-04 Max=1.60D-02 LinEq1: Iter= 4 NonCon= 1 RMS=1.99D-04 Max=1.26D-02 LinEq1: Iter= 5 NonCon= 1 RMS=1.04D-04 Max=6.16D-03 LinEq1: Iter= 6 NonCon= 1 RMS=5.06D-05 Max=1.92D-03 LinEq1: Iter= 7 NonCon= 1 RMS=2.28D-05 Max=1.17D-03 LinEq1: Iter= 8 NonCon= 1 RMS=1.49D-05 Max=1.22D-03 LinEq1: Iter= 9 NonCon= 1 RMS=1.48D-05 Max=1.31D-03 LinEq1: Iter= 10 NonCon= 1 RMS=7.12D-06 Max=2.68D-04 LinEq1: Iter= 11 NonCon= 1 RMS=3.05D-06 Max=1.36D-04 LinEq1: Iter= 12 NonCon= 1 RMS=2.41D-06 Max=1.46D-04 LinEq1: Iter= 13 NonCon= 1 RMS=1.28D-06 Max=5.58D-05 |
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 303求調(diào)劑 +3 | 睿08 2026-03-17 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 0703化學(xué)調(diào)劑 ,六級(jí)已過(guò),有科研經(jīng)歷 +8 | 曦熙兮 2026-03-15 | 8/400 |
|
|
[考研] 299求調(diào)劑 +4 | △小透明* 2026-03-17 | 4/200 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 5/250 |
|
|
[考研] 【0856】化學(xué)工程(085602)313 分,本科學(xué)科評(píng)估A類院;瘜W(xué)工程與工藝,誠(chéng)求調(diào)劑 +7 | 小劉快快上岸 2026-03-11 | 8/400 |
|
|
[碩博家園] 深圳大學(xué)碩士招生(2026秋,傳感器方向,僅錄取第一志愿) +4 | xujiaoszu 2026-03-11 | 9/450 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[基金申請(qǐng)]
今年的國(guó)基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑,求各位老師收留 +8 | 秋有木北 2026-03-14 | 8/400 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
[考研] 材料080500調(diào)劑求收留 +3 | 一顆meteor 2026-03-13 | 3/150 |
|
|
[考研] 331求調(diào)劑(0703有機(jī)化學(xué) +5 | ZY-05 2026-03-13 | 6/300 |
|
|
[考研] 學(xué)碩285求調(diào)劑 +13 | Wisjxn 2026-03-12 | 46/2300 |
|
|
[考研] [0860]321分求調(diào)劑,ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
[考研] 333求調(diào)劑 +3 | 球球古力 2026-03-11 | 3/150 |
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級(jí)已過(guò) 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-12 | 3/150 |
|