 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 1471 | 回復(fù): 0 | ||
zhlu929銅蟲 (初入文壇)
|
[求助]
用高斯搜索過渡態(tài)反應(yīng)物和產(chǎn)物自旋多重度不一樣怎么辦？
|
|
如題，用高斯搜索過渡態(tài)，反應(yīng)物自旋多重度為1，產(chǎn)物自旋多重度為5，過渡態(tài)自旋多重度應(yīng)該怎么設(shè)？我設(shè)自旋多重度為1搜到過渡態(tài)后，跑IRC很容易搜到反應(yīng)物，但是怎么跑到自旋多重度為5的產(chǎn)物很迷惑，我可以用過渡態(tài)的結(jié)構(gòu)把自旋多重度改為5來進(jìn)行計(jì)算嗎？這是否意味著我的結(jié)構(gòu)在過渡態(tài)瞬間電子結(jié)構(gòu)改變？這樣做對(duì)嗎？ 附上計(jì)算腳本 1.先掃描鍵長 #p opt=modredundant b3lyp/gen geom=connectivity pseudo=read scfcyc=200 symm=veryloose Fe(H2B(pz)2)2phen_2 6-31g 0 1 Fe 0.00000000 0.00000000 2.62521400 S -1.24849800 1.31269500 -6.46666200 S 1.24849800 -1.31269500 -6.46666200 N -1.25380300 0.39111200 1.08503100 分子坐標(biāo)好長，不寫了 2.找到可能的過渡態(tài)，計(jì)算頻率 #p freq ub3lyp/gen pseudo=read scfcyc=200 symm=veryloose Fe_phen_2freq 0 1 Fe 2.54767100 0.00002800 0.00002100 S -6.28746100 -2.17776300 0.42803200 S -6.28746900 2.17775800 -0.42805800 這里還是用自旋多重度=1計(jì)算的，結(jié)果發(fā)現(xiàn)只有1個(gè)-200的虛頻，很開心 3.人為把自旋多重度改為5進(jìn)行優(yōu)化計(jì)算，跑出來的結(jié)果跟產(chǎn)物很相像，但是有一個(gè)-8.71的虛頻，震動(dòng)方向如fig1紅色箭頭所示，我是否需要把它重新調(diào)關(guān)鍵詞優(yōu)化？要怎么調(diào)關(guān)鍵詞呢？可以直接沿震動(dòng)方向調(diào)坐標(biāo)優(yōu)化嗎？ #p opt freq=noraman ub3lyp/gen iop(3/24=1,5/13=1) pseudo=read symm=veryloose ts1 6-31+g* optfreq 0 5 Fe 0.00000000 0.00000000 2.53767100 S 2.27317600 -0.66823400 -6.22746100 S -2.27317600 0.66823400 -6.22746100 4.分別跑IRC找產(chǎn)物和反應(yīng)物，問題來了 a）自旋多重度設(shè)為1找反應(yīng)物 #p irc=(CalcFC,forward,maxpoints=10,stepsize=20) b3lyp/gen pseudo=read ts1_irc forward 0 1 Fe 2.54767100 0.00002800 0.00002100 S -6.28746100 -2.17776300 0.42803200 S -6.28746900 2.17775800 -0.42805800 很順利結(jié)束 b）自旋多重度設(shè)為5找產(chǎn)物 #p irc=(CalcFC,reverse,maxpoints=30,stepsize=20,maxcycle=400) IOp(5/13=1) b3lyp/gen pseudo=read ts1_irc reverse 0 5 Fe 2.54767100 0.00002800 0.00002100 S -6.28746100 -2.17776300 0.42803200 S -6.28746900 2.17775800 -0.42805800 各種報(bào)錯(cuò)，各種調(diào)參數(shù),報(bào)錯(cuò)信息如下： Cartesian Forces: Max 0.027937985 RMS 0.005666714 Z-matrix is all fixed cartesians, so copy forces. Leave Link 716 at Wed Dec 9 11:03:35 2015, MaxMem= 4026531840 cpu: 1.0 (Enter /public/software/g09/l123.exe) IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC SCF convergence failure! Path Direction Number 1 Point Number 2 Unrecoverable SCF failure Error termination via Lnk1e in /public/software/g09/l123.exe at Wed Dec 9 11:03:35 2015. Job cpu time: 24 days 21 hours 43 minutes 58.5 seconds. File lengths (MBytes): RWF= 17086 Int= 0 D2E= 0 Chk= 106 Scr= 2 求助各位大神  fig1.JPG |
回復(fù)此樓» 猜你喜歡
可生物降解聚酯正在重塑現(xiàn)代醫(yī)療器械
已經(jīng)有0人回復(fù)
-
26屆博士申請(qǐng)
已經(jīng)有4人回復(fù)
-
有機(jī)高分子材料論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有72人回復(fù)
-
PLLA 微球在生物醫(yī)學(xué)領(lǐng)域的研究與應(yīng)用
已經(jīng)有1人回復(fù)
-
PLLA 微球的材料特性與制備工藝解析
已經(jīng)有0人回復(fù)
-
醫(yī)用聚己內(nèi)酯（PCL）全解析：從兩段式降解到組織工程應(yīng)用
已經(jīng)有0人回復(fù)
-
反相硅膠如何再生？
已經(jīng)有1人回復(fù)
-
哈爾濱理工大學(xué)2026年研究生調(diào)劑，材料科學(xué)與化學(xué)工程學(xué)院研究生調(diào)劑
已經(jīng)有3人回復(fù)
-
武漢一本高校招收材料、化學(xué)、高分子、紡織、生物相關(guān)專業(yè)碩士生
已經(jīng)有0人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 zhlu929 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 281求調(diào)劑（0805） +6 | 煙汐憶海 2026-03-16 | 16/800 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
|
[考研] 294求調(diào)劑材料與化工專碩 +5 | 陌の森林 2026-03-18 | 5/250 |
|
|
|
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +3 | xbxudjdn 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +4 | 1孫悟空 2026-03-17 | 4/200 |
|
|
|
[考研] 293求調(diào)劑 +11 | zjl的號(hào) 2026-03-16 | 16/800 |
|
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
|
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
|
[基金申請(qǐng)]
今年的國基金是打分制嗎？
 50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-12 | 6/300 |
|
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油！? 2026-03-15 | 5/250 |
|
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
|
[考研] 招收0805（材料）調(diào)劑 +3 | 18595523086 2026-03-13 | 3/150 |
|
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +5 | 一定有學(xué)上- 2026-03-12 | 5/250 |
|
|
|
[考研] 085600材料與化工 309分請(qǐng)求調(diào)劑 +7 | dtdxzxx 2026-03-12 | 8/400 |
|