应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】高斯计算td的问题
101/1返回列表
查看: 3455  |  回复: 9
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

ellemiu

铁杆木虫 (正式写手)

[交流] 【求助】高斯计算td的问题

我想算一个物质在乙腈溶液中的紫外图谱
输入文件如下

#p rb3lyp/genecp td=(singlet,nstates=15) scrf=(iefpcm,solvent=acetonitrile)

0 1
分子坐标

C H N 0
6-31G(d)
****
Mn Cl 0
6-31G+(d,p)

运算时总是提示下面的错误

Leave Link  202 at Mon Mar 28 10:31:45 2011, MaxMem=  268435456 cpu:       0.3
(Enter /opt/g03/l301.exe)
General basis read from cards:  (5D, 7F)
Centers:   1  2  3  5  6  7  8 10 11 12 13 14 16 17 18 19 20 22 23 24
Centers:  25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44
Centers:  45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64
Centers:  65 66 67 68 69 70 71 72 73 74 75 76  4  9 15 21
6-31G(d)
****
Centers:  77 78
6-31g+(d,p)

  QPERR --- A SYNTAX ERROR WAS DETECTED IN THE INPUT LINE.

'
Last state="Top"
TCursr=  746 LCursr=    0
Error reading general basis specification.
Error termination via Lnk1e in /opt/g03/l301.exe at Mon Mar 28 10:31:47 2011.
Job cpu time:  0 days  0 hours  0 minutes  7.2 seconds.
File lengths (MBytes):  RWF=     13 Int=      0 D2E=      0 Chk=      1 Scr=      1

注：分子坐标是先通过b3lyp/6-31g几何优化后，再导出的得到的

想知道如何改如何修正才能得到我想要的结果，现在这样的错误原因是什么，希望各位达人能够不吝赐教。
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼2011-03-28 11:11:47
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

abbott

金虫 (著名写手)
★
ellemiu(金币+2): 哪里的语法错了，能不能说的详细一点 2011-03-28 18:26:31
375642546(金币+1): 谢谢! 2011-04-08 09:58:16
语法错误
回复此楼
2楼2011-03-28 12:00:36
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

daiyulan85

木虫 (正式写手)
★
ellemiu(金币+1): 6-31g都是常用的基组，不太明白你说的错误是指？ 2011-03-29 00:36:21
375642546(金币+1): 谢谢！ 2011-04-08 09:58:06
Error reading general basis specification

会不会是基组错误？
赞一下(1人) 回复此楼
3楼2011-03-28 20:09:03
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

blueidea

禁虫 (小有名气)
ellemiu(金币+1): 多谢来访。。 2011-03-29 00:37:40
本帖内容被屏蔽
4楼2011-03-29 00:35:26
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

shangji2

金虫 (正式写手)
ellemiu(金币+1): 哥，指点下迷津，我加你 2011-03-29 17:12:13
这个小问题，QQ谈，多给点金币啊，83842088
赞一下 回复此楼
5楼2011-03-29 09:21:37
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

woxiangfei

木虫 (文坛精英)
★
375642546(金币+1): 谢谢！ 2011-04-08 09:57:37
ellemiu(金币+10): 谢谢，多谢指点，，看了你写的明白了很多，不过还是没有发现我错误的地方 2011-04-08 17:15:51
给你一个例子：

螢光範例： Acetaldehyde ?牡谝患ぐl態（n→π*） 放射的螢光
現在?硌芯窟@個迴圈：

Acetaldehyde 激發和放射迴圈
這裡感興趣的過程是放射（emission），這個範例顯示如何完整激發和放射迴圈的溶液效應進行研究。
步驟 1: 基態幾何結構優選和頻率計算（平衡溶解）。這是分子包含 PCM 平衡溶解的基態標準的 Opt Freq 計算。
%chk=01-ac
# B3LYP/6-31+G(d,p) Opt Freq SCRF=(Solvent=Ethanol)

Acetaldehyde ground state

0 1
C
C,1,RA
X,2,1.,1,A
O,2,RB,3,A,1,180.,0
X,1,1.,2,90.,3,0.,0
H,1,R1,2,A1,5,0.,0
H,1,R23,2,A23,5,B23,0
H,1,R23,2,A23,5,-B23,0
H,2,R4,1,A4,3,180.,0

RA=1.53643
RB=1.21718
R1=1.08516
R23=1.08688
R4=1.10433
A=62.1511
A1=110.51212
A23=109.88119
A4=114.26114
B23=120.56468
步驟 2: 垂直激發，線性回應溶解模型。這是垂直激發的 TD-DFT 計算，因此採用基態平衡幾結構，加上預設的溶解模型：線性回應，非平衡狀態。進行單結構 TD-DFT 計算，預設選項採用非平衡溶解模型。這計算工作的結果用?碚页瞿男顟B是感興趣的狀態，以及這些狀態的次序。這些結果提供激發態溶解的合理描述，但還是不及特定狀態的溶解計算所得的結果。在這個範例中， 越遷 n->π* 狀態是第一個激發態。下一步，使用特定狀態方法產生垂直激發步驟較佳的描述。
%chk=02-ac
# B3LYP/6-31+G(d,p) TD=NStates=6 SCRF=(Solvent=Ethanol)
Geom=Check Guess=Read

Acetaldehyde: linear response vertical excited states

0 1
步驟 3: 垂直激發的特定狀態溶解模型。需要兩個工作步驟：第一先做分子基態計算，在 PCM 輸入段落中指定NonEq=write，將基態非平衡溶解模型的計算結果儲存起?怼５诙糠诌M行真正的指定狀態計算，用 NonEq=read ?臋z查檔讀回所需的資料進行非平衡溶解模型計算。
%chk=03-ac
# B3LYP/6-31+G(d,p) SCRF=(Solvent=Ethanol,Read)
Geom=Check Guess=Read

Acetaldehyde: prepare for state-specific non-eq solvation
by saving the solvent reaction field from the ground state

0 1

NonEq=write

--link1--
%chk=03-ac
# B3LYP/6-31+G(d,p) TD(NStates=6,Root=1)
SCRF=(Solvent=Ethanol,StateSpecific,Read)
Geom=Check Guess=Read

Acetaldehyde: read non-eq solvation from ground state and
compute energy of the first excited with the state-specific method

0 1

NonEq=read
步驟 4: 激發態幾何結構優選計算。 接下?恚M行一次 TD-DFT 幾何優選計算，採用平衡線性回應溶解模型，找出激發態位能曲面的最低能量結構。因為這是一個 TD-DFT 邀選計算，程式預設的選項是平衡溶解模型。像這種典型的情況，基態的分子有平面對稱性，但在激發態時，這對稱性就被破壞了，因此在優選計算開始的時候就將基態幾何結構做稍微的擾動，將對稱性破壞，往後的優選計算才能順利進行。
%chk=04-ac
# B3LYP/6-31+G(d,p) TD=(Read,NStates=6,Root=1) SCRF=(Solvent=Ethanol)
Geom=Modify Guess=Read Opt=RCFC

Acetaldehyde: excited state opt
Modify geometry to break Cs symmetry
since first excited state is A"

0 1

4 1 2 3 10.0
5 1 2 7 -50.0
步驟 5: 激發態分子結構的振動頻率。現在執行頻率計算，確認步驟 4 找到的幾何結構是在能量曲面的最低點。若有需要，計算結果也可被用?磉M行 Franck-Condon 計算（見後面的討論）。這是一個數值頻率計算。
%chk=05-ac
# B3LYP/6-31+G(d,p) TD=(Read,NStates=6,Root=1) Freq
SCRF=(Solvent=Ethanol) Geom=Check Guess=Read

Acetaldehyde excited state freq

0 1
步驟 6: 溶液中特定狀態的電子放射(部分 1)。這步驟進行分子激發態在其平衡幾何結構上的特定狀態平衡溶解計算，用 PCM NonEq=write 將溶解數據寫道檢查檔，給下一步計算工作使用。
%chk=06-ac
# B3LYP/6-31+G(d,p) TD=(Read,NStates=6,Root=1)
SCRF=(Solvent=Ethanol,StateSpecific,Read)
Geom=Check Guess=Read

Acetaldehyde emission state-specific solvation
at first excited state optimized geometry

0 1

NonEq=write
步驟 7: 電子放射到分子基態（部分 2）。最後，用非平衡溶解模型計算基態能量，使用激發態的幾何結構和激發態的靜態溶解模型。
%chk=07-ac
# B3LYP/6-31+G(d,p) SCRF=(Solvent=Ethanol,Read) Geom=Check Guess=Read

Acetaldehyde: ground state non-equilibrium
at excited state geometry.

0 1

NonEq=read
步驟 1，2，和 4 足夠用?碛嬎銡庀喾肿拥募ぐl和放射能量（再用步驟 5 確認穩定點的特性）。當要考慮溶液效應時，這幾步驟的計算就不夠了，因為步驟 4 計算的能量相當於溶劑反應場的基態，而回應激發態電荷分?旬a生在反應場內發生的放射。步驟 6 和 7 對這溶劑效應做適當的處理。
赞一下(1人) 回复此楼
6楼2011-04-07 22:03:25
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

yongma2008

木虫 (著名写手)
375642546: 其实也有这种基组，只是不常用而已！ 2011-04-08 09:43:58
ellemiu(金币+1): 谢谢来访，其实是我写错，不过我刚研究了下说明书，这个基组也是存在的 2011-04-08 17:20:02
6-31G+(d,p)，第一次见这种写法
赞一下(1人) 回复此楼
7楼2011-04-08 03:30:55
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

鱼妃

金虫 (著名写手)
★
375642546(金币+1): 谢谢！ 2011-04-08 09:56:44
ellemiu(金币+2): 谢谢指点，其实是我写错了，不过现在长见识了。 2011-04-08 17:24:53
引用回帖:
Originally posted by yongma2008 at 2011-04-08 03:30:55:
6-31G+(d,p)，第一次见这种写法

楼主想用的机组 不是波普尔基组6-31+G(d,p) 而是含George Petersson等人发明的基组
这种基组写作6-31G+，6-31H+等形式
这种基组在高斯中写作6-31G(d',p')
赞一下(1人) 回复此楼
8楼2011-04-08 09:50:38
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

375642546

木虫 (正式写手)
ellemiu(金币+3): 谢谢指教，我查了说明书确实如你所说，不过我后来将基组改成了Lanl2dz计算Mn原子，用6-31g计算C、H、N、Cl原子，依然出现同样的错误 2011-04-08 17:23:50
我来说说我的看法吧，首先lz的出错原因：语法错误，基组关键词写错了。
lz应该是想写6-31+g（d,p）吧？结果写成了6-31g+（d,p）了？
其实想用6-31g+（d,p）也是可以的，这是George Petersson等人开发的基组，高斯中也包含这种基组，只是这种基组不这样写罢了。如6-31g+的关键词是6-31g(d')。

参考：高斯说明书。
赞一下 回复此楼
9楼2011-04-08 09:55:50
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

yongma2008

木虫 (著名写手)
引用回帖:
Originally posted by 鱼妃 at 2011-04-08 02:50:38:
楼主想用的机组 不是波普尔基组6-31+G(d,p) 而是含George Petersson等人发明的基组
这种基组写作6-31G+，6-31H+等形式
这种基组在高斯中写作6-31G(d',p')

谢谢指教。。
回复此楼
10楼2011-04-08 12:41:15
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ellemiu 的主题更新
101/1返回列表
信息提示
关闭
请填处理意见
关闭