| 查看: 3748 | 回复: 9 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】高斯计算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几何优化后，再导出的得到的 想知道如何改如何修正才能得到我想要的结果，现在这样的错误原因是什么，希望各位达人能够不吝赐教。 |
回复此楼» 猜你喜欢
华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员
已经有11人回复
-
华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员
已经有10人回复
-
物理化学论文润色/翻译怎么收费?
已经有156人回复
-
200年来人们一直无法破解的难题： 气液临界点状态方程
已经有7人回复
-
关于引力延迟下双星系统角动量不守恒的疑惑---求高手解答！
已经有10人回复
-
求标准卡片Cu0.4In0.4Zn1.2S2 PDF#32-0340
已经有2人回复
-
量子点电致发光器件（QLED）：蓝光量子点（QDs）和纳米氧化锌电子传输材料
已经有0人回复
-
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建
已经有0人回复
-
求ENDNOTE破解软件安装包
已经有4人回复
-
河海大学 （211.双一流高校）钙钛矿光电实验室丁勇团队欢迎有志青年！
已经有3人回复
-
黑龙江省原子能研究院化学课题组招收2026级化学专业或化学相关专业硕士生
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于化合物超极化率的计算方法
已经有6人回复
- 高斯算紫外误差很大，咋办？
已经有14人回复
- 求高斯td计算紫外的问题
已经有16人回复
- 高斯计算激发态并行问题
已经有6人回复
- 【求助】关于激发态优化及发射光谱的计算
已经有36人回复
- 【求助】高斯计算时无缘无故中断的问题
已经有8人回复
- 【求助】高斯计算出现Fatal Problem: The smallest alpha delta epsilon is -0.785373
已经有7人回复
- 【求助】G09荧光光谱例子 运行出错
已经有25人回复
- 【求助】含时密度泛函计算金团簇的激发态输出出错
已经有8人回复
- 【求助】计算荧光光谱的疑惑，求教！
已经有18人回复
- 【求助】关于TD-DFT模拟Uv-Vis光谱的问题
已经有45人回复
- 【求助】如何在以下高斯03输出文件中得到二阶超极化率
已经有15人回复
- 【求助】单重激发态的高斯算法
已经有29人回复
» 抢金币啦！回帖就可以得到:
-
祈福2026国基面上项目
+5/435
-
散金币，求好运，祝面上顺利！
+2/96
-
物理学 调剂
+1/80
-
博士后出站经验请教
+1/78
-
河南工学院诚聘优秀凝聚态物理和电子科学与技术相关专业硕博研究生
+2/40
-
0854电子信息调剂，闽南师大光电芯片研发实验室，柯少颖教授课题组，集成电路
+1/40
-
华南师范大学（211）博士招生- 电子、自动化、机械、生物学、物理相关专业
+2/38
-
上海交通大学电催化方向博后招聘-有锂电、高分子、燃料电池背景者优先
+1/31
-
2026年工科硕士调剂-上海大学全国重点实验室团队-材料数据挖掘方向-研究生3-5人
+1/30
-
南京理工大学招收2026级硕士生-张光普博导
+2/30
-
医学327分专硕调剂
+1/16
-
山东省植物基蛋白技术创新中心课题组联合培养硕士研究生招募启事
+1/12
-
【2026博士招生/博后招聘】北京航空航天大学潘彪课题组——AI芯片设计方向
+1/11
-
聊城大学化院青椒招收化学相关专业调剂学生
+1/7
-
教育部长江学者和创新团队发展计划”入选团队招收 材料，化学与化工博士研究生
+1/7
-
青岛科技大学可持续高分子团队 考研招生
+1/7
-
香港大学David Srolovitz教授课题组招聘「香江学者」博士后
+1/7
-
浙江大学高纯分离团队诚聘博士后
+1/5
-
东北林业大学（211，双一流）宋老师课题组招收高分子、化学、复合材料等方向研究生
+1/4
-
招收环工、环科、化工、应化、化学、生工、生物、材化等相关专业研究生，要求过国家线
+1/3
3楼2011-03-28 20:09:03
5楼2011-03-29 09:21:37
★
375642546(金币+1): 谢谢！ 2011-04-08 09:57:37
ellemiu(金币+10): 谢谢，多谢指点，，看了你写的明白了很多，不过还是没有发现我错误的地方 2011-04-08 17:15:51
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 對這溶劑效應做適當的處理。 |
6楼2011-04-07 22:03:25
7楼2011-04-08 03:30:55
