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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 高斯09寻找过渡态，急需帮助，十分感谢
51/1返回列表
查看: 2534  |  回复: 6
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

浅秋0321

新虫 (初入文坛)

[求助] 高斯09寻找过渡态，急需帮助，十分感谢

各位前辈，第一次发帖，恳请大家多多帮助。如果表述有什么不清楚的地方，请随时告知哈。谢谢
下面的结构是我的基础结构（complex），其中1号Mn是3价的，2号Mn是4价的，每个Mn旁边都一个配体水。我想找一个醋酸根取代一个Mn旁边的配体水的过渡态，可是一个都找不到。尝试了ts的方法，我自己设计的ts结构是一个六元环的结构，如下的ts结构。不知道这种六元环的结构是合理的吗？（我比较傻的对这个猜测的ts结构没有优化，就算了一下频率，只有一个虚频，而且振动方向对应着取代反应的成断键，感觉应该还是有些靠谱的。可是一用ts方法进行过渡态优化的时候醋酸根就会和Mn原子连上，虚频也变了）。因为之前从没算过过渡态，后来尝试用qst3的方法，感觉也是一头雾水。恳请各位前辈给指调明路，十分感谢！

（complex）
3 2
Mn                -1.12334000    0.00185400    0.61431900
Mn                 1.31636300    0.00351200   -0.58131700
O                 -0.35867700    0.00586600   -1.09143600
O                  0.67071600    0.00011400    1.04360500
O                 -1.46711200    0.00675300    2.71027200
O                  1.85023300    0.00828800   -2.62270800
N                 -1.83313200   -2.16576500    0.73532600
N                 -3.09747800   -0.00502100   -0.12647800
N                 -1.84747100    2.16631200    0.72834700
N                  1.76520000   -1.99452800   -0.62427100
N                  3.23842700    0.00022500    0.02852700
N                  1.76870600    2.00088200   -0.61379300
C                 -3.02462700   -2.39354900    0.11122500
C                 -3.67334000   -1.19334900   -0.46131500
C                 -4.83832500   -1.22202800   -1.23766300
C                 -5.41336500   -0.01541000   -1.64643400
C                 -4.84427800    1.19650300   -1.24508700
C                 -3.67927500    1.17836600   -0.46841500
C                 -3.03738000    2.38513400    0.09805800
C                  3.00421100   -2.33677300   -0.13358100
C                  3.81673900   -1.19200100    0.31398900
C                  5.06881000   -1.22232900    0.93815800
C                  5.68887700   -0.00514500    1.25447900
C                  5.07081300    1.21476900    0.94483200
C                  3.81871000    1.18991600    0.32048500
C                  3.00816500    2.33843700   -0.12099300
C                  0.96399800   -2.93660900   -1.16622800
C                  3.43256000   -3.66490800   -0.15515500
C                  1.35548400   -4.27699500   -1.22880900
C                  2.59915700   -4.64708400   -0.70747700
C                  0.96926200    2.94713200   -1.15104600
C                  3.43874100    3.66593500   -0.13545300
C                  1.36302500    4.28717400   -1.20646300
C                  2.60714300    4.65241800   -0.68279800
C                 -1.21982800   -3.17389000    1.39385700
C                 -3.60790800   -3.66715300    0.11247900
C                 -2.96650200   -4.71208300    0.78776000
C                 -1.75835400   -4.46205800    1.44850200
C                 -3.62585200    3.65632400    0.08856300
C                 -1.24138000    3.18107000    1.38338500
C                 -1.78548400    4.46731200    1.42757000
C                 -2.99161200    4.70822700    0.75979400
H                 -5.29409900   -2.16097300   -1.52222200
H                 -5.30479900    2.13138100   -1.53530000
H                  5.55412300   -2.15901600    1.18082200
H                  5.55767000    2.14930700    1.19263500
H                  2.72759700    0.00921300   -3.04467500
H                  1.12298100    0.01014200   -3.27143000
H                 -0.71977700   -0.00240200    3.33440800
H                 -2.32793000    0.01626100    3.16463900
H                 -6.30841300   -0.01948300   -2.25755500
H                  6.65670900   -0.00727700    1.74244900
H                  0.69919000   -5.00868100   -1.68325300
H                  2.92722700   -5.67954400   -0.74357700
H                  4.40995100   -3.93281100    0.22657600
H                  0.01223200   -2.59270100   -1.54716500
H                  0.70809800    5.02234500   -1.65723700
H                  2.93695000    5.68450300   -0.71331800
H                  4.41646900    3.93020600    0.24794400
H                  0.01705400    2.60677200   -1.53405500
H                 -3.41465200   -5.69875100    0.81223400
H                 -1.25165300   -5.24285100    2.00260300
H                 -1.28454200    5.25360700    1.97912900
H                 -3.44379900    5.69322300    0.77585200
H                 -4.55864600   -3.84619000   -0.37298300
H                 -0.28778400   -2.92774000    1.88816500
H                 -4.57511700    3.82838700   -0.40224200
H                 -0.31075600    2.94187900    1.88373000

猜测的ts结构：
#p freq=noraman ub3lyp/3-21g  scf(maxcycle=300)

2 8
Mn                 0.99600100   -0.00142200    0.28129700
Mn                -1.51136400   -0.00076300   -0.80619500
O                  0.13277600    0.02053700   -1.40518600
O                 -0.82216400   -0.06990700    0.76588000
O                  3.45206400    0.05322900    1.58323100
O                 -2.14740300    0.01468500   -2.84578900
N                  1.70186100    2.15975700    0.27282700
N                  2.83793000    0.00370700   -0.74613000
N                  1.70310300   -2.14573100    0.27451100
N                 -1.93571600    1.96670000   -0.80117400
N                 -3.39106700   -0.00164100   -0.10053300
N                 -1.92725200   -1.96225100   -0.83350100
C                  2.79667200    2.38207400   -0.49846500
C                  3.35767600    1.17995900   -1.15668800
C                  4.39974600    1.21195100   -2.09009100
C                  4.91029200    0.00273700   -2.57972000
C                  4.39655300   -1.20333100   -2.09862800
C                  3.35976700   -1.17746400   -1.16566700
C                  2.79374600   -2.37714300   -0.49956400
C                 -3.14381200    2.32565800   -0.25274000
C                 -3.94478300    1.18348700    0.22913700
C                 -5.15548100    1.20656800    0.92150100
C                 -5.75247900   -0.01587600    1.25949000
C                 -5.15139700   -1.23101100    0.90645700
C                 -3.93875200   -1.19355500    0.21422000
C                 -3.13408700   -2.32615900   -0.28332800
C                 -1.14257700    2.89923700   -1.37624300
C                 -3.54976200    3.65549900   -0.24768700
C                 -1.51837700    4.24254700   -1.41409800
C                 -2.72997100    4.62502400   -0.83232300
C                 -1.13177800   -2.87790500   -1.41183300
C                 -3.53225000   -3.66106900   -0.28653700
C                 -1.49753400   -4.22863800   -1.45469000
C                 -2.70583900   -4.62225500   -0.88177600
C                  1.17124700    3.16588200    0.99565200
C                  3.37292900    3.65370500   -0.57961900
C                  2.82186100    4.70230300    0.16145800
C                  1.70515700    4.45565800    0.96774900
C                  3.36436200   -3.65121300   -0.58506000
C                  1.16849400   -3.15311600    1.00215100
C                  1.69981900   -4.44091300    0.97602700
C                  2.80926600   -4.69400400    0.16525000
H                  4.81049600    2.15046900   -2.43604500
H                  4.80647600   -2.14158900   -2.44648600
H                 -5.62436300    2.13974700    1.20388400
H                 -5.61599300   -2.16986200    1.17733700
H                 -3.04130400    0.22120600   -3.26777700
H                 -1.42418100    0.01296100   -3.52200600
H                  3.12483400   -0.00340900    2.56897100
H                  4.45530500   -0.13263300    1.60554300
H                  5.70364000    0.00574500   -3.31143500
H                 -6.68968600   -0.02354300    1.80732500
H                 -0.86801500    4.96581400   -1.88700000
H                 -3.04168600    5.66145700   -0.84999300
H                 -4.50554300    3.93471200    0.18075000
H                 -0.21865800    2.52943200   -1.79168500
H                 -0.84651600   -4.94567600   -1.93845200
H                 -3.00967200   -5.66086100   -0.90039100
H                 -4.48144900   -3.95050500    0.14189700
H                 -0.20429400   -2.50959300   -1.83239200
H                  3.26510700    5.68865100    0.11866700
H                  1.27390100    5.24111600    1.57373600
H                  1.26400400   -5.22215300    1.58908600
H                  3.24880700   -5.68258400    0.12883000
H                  4.25303800    3.82512400   -1.18713000
H                  0.31725100    2.91967500    1.61139700
H                  4.23730600   -3.83113100   -1.19477900
H                  0.31423700   -2.89044700    1.61389800
C                  1.34500500   -0.02113800    3.93392900
O                  2.55470800   -0.10149500    4.30544000
O                  0.95783800    0.06202400    2.66014500
C                  0.18076200   -0.01704300    4.94018800
H                 -0.42883800    0.87853200    4.79033200
H                 -0.45814300   -0.89278500    4.76762400
H                  0.57662500   -0.04236500    5.96623100

用ts方法优化的命令行#p opt=(ts,calcfc,noeigentest) freq=noraman ub3lyp/3-21g scf(maxcycle=300)
尝试了固定原子，用iop(1/8=7)缩小步长的方法，都不行。不知道是不是哪里出了错。
回复此楼

» 收录本帖的淘帖专辑推荐

过渡态计算 催化与分子模拟 量化计算有意义的讨论

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2012-06-02 09:52:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

scottfreedom

金虫 (正式写手)

【答案】应助回帖

引用回帖:
4楼: Originally posted by 浅秋0321 at 2012-06-04 08:43:36
十分感谢您！请问“在TS前，使用Redundant Coordinate进行冗余坐标设定，将你所要找的相关过渡态的关键原子的坐标固定住，再Opt”，我将所猜测的ts初始结构固定了参与过渡态的原子，也尝试了下固定参与过渡态的键长 ...

这种情况我也遇见过，原因一般都是固定的原子的键长，键角的问题，在Gaussian里面，有时候，可能就是那么一点点差距，得到的结果就会有很大的不同，所以，用Redundant Cordinate优化找过渡态，是一件很耗费脑力和精力的事情，如果实在是很烦，也可以尝试scan，不过这个方法我也只是听说过，具体的，还是要等高手来解决，建议楼主可以在坛子里搜搜过渡态的相关帖子，说不定会有答案的啦！慢慢来啦！不急的！
赞一下 回复此楼
Terse
5楼2012-06-04 20:56:22
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

浅秋0321

新虫 (初入文坛)
怎么没有人回复呢？
赞一下 回复此楼
2楼2012-06-02 14:32:36
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

scottfreedom

金虫 (正式写手)

【答案】应助回帖

★ ★ ★
感谢参与，应助指数 +1
gmy1990: 金币+3 2012-06-03 09:21:26
找过渡态，方法一般用有这么两种，第一，就是楼主你所用的TS，但是，看了楼主你的说明，建议最好在TS前，使用Redundant Coordinate进行冗余坐标设定，将你所要找的相关过渡态的关键原子的坐标固定住，再Opt，做一次优化，再Ts，这样靠谱点，然后别忘了，Ts后，做一个IRC路径追踪，有时候看虚频不一定能看出来。第二，用QST2或者QST3，这两种方法比较偷懒，不过，前提是，你所提供的初始物和产物要进行过Opt，如果你用了QST3，还要给一个你猜测的构型，当然，这个猜测的不用你来优化，个人浅酌，仅供参考。
赞一下(2人) 回复此楼

» 本帖已获得的红花（最新10朵）

浅秋0321
Terse
3楼2012-06-02 20:28:02
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

浅秋0321

新虫 (初入文坛)
送鲜花一朵
引用回帖:
3楼: Originally posted by scottfreedom at 2012-06-02 20:28:02
找过渡态，方法一般用有这么两种，第一，就是楼主你所用的TS，但是，看了楼主你的说明，建议最好在TS前，使用Redundant Coordinate进行冗余坐标设定，将你所要找的相关过渡态的关键原子的坐标固定住，再Opt，做一次 ...

十分感谢您！请问“在TS前，使用Redundant Coordinate进行冗余坐标设定，将你所要找的相关过渡态的关键原子的坐标固定住，再Opt”，我将所猜测的ts初始结构固定了参与过渡态的原子，也尝试了下固定参与过渡态的键长，但是不知道是不是我的初始有问题，优化了好多步（用Guassview查看是450步），结构已经跑掉了，然后是l9999错误。应该怎么办呢？

用qst3算的时候，看很多经验贴上说需要确定反应物和产物一定是在反应坐标的两侧，这一点让我有点困惑，不知道对于我这个体系，应该怎么处理呢？

再次感谢您的帮助哈！
赞一下 回复此楼
4楼2012-06-04 08:43:36
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭