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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】键长键角等固定后的优化..悬赏30...【已经解决～感谢～】
51/1返回列表
查看: 4280  |  回复: 22
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

亍寞

金虫 (小有名气)

[交流] 【求助】键长键角等固定后的优化..悬赏30...【已经解决～感谢～】

新手问题...........
现在正在找一个过渡态，没有固定键长键角之前，构型偏离较大最后就不能收敛，或者分子散了，若构型偏离较小就会优化直接到稳定态
现在准备固定某个键长后继续优化，用gaussview建立模型之后，用自带的redundant，然后手动再修改某几个数值，输入文件中有：
D 4 1 2 3 -120.000000 F
这样的描述，计算过程中也没有出错，但是最后键角仍旧偏离较大，键长也是如此

想请教一下这个如果固定键长之后优化，gaussian是在这个输入结构附近寻找收敛点么？若在这个点附近找不到最小点是不是就会破坏这个限制？加上了键长键角限制之后和没有加键长键角有什么区别呢？


或者....悲情的就是固定方法不对，输入文件不对？

求大神详解...或者能麻烦提供详解的链接地址么？....

[ Last edited by 亍寞 on 2011-3-12 at 21:17 ]
回复此楼

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

量化计算有意义的讨论 Gaussian

» 猜你喜欢

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

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

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

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

亍寞

金虫 (小有名气)
引用回帖:
Originally posted by yongma2008 at 2011-03-12 19:23:36:
反映坐标找到了就按照反映坐标，一步一步的固定键长来优化，优化后会在反映坐标过程中找到能量最大的点，然后在用过渡态优化。一般能找到。优化过程不很确定的话，建议先不要算频率，频率回头用优化好的构型的检查 ...

谢谢回复，这个原理我明白，是输入文件的问题，麻烦您看下，这是我刚才用乙烷分子做了尝试，想固定两个碳原子之间的距离之后优化，这是输入文件

%mem=3200MB
%nprocshared=2
#p opt(modredundant) freq b3lyp/6-31g scf=tight

Title Card Required

0 1
C
H                  1              B1
H                  1              B2    2              A1
H                  1              B3    2              A2    3              D1
C                  1              B4    2              A3    3              D2
H                  5              B5    1              A4    2              D3
H                  5              B6    1              A5    2              D4
H                  5              B7    1              A6    2              D5

   B1             1.07000000
   B2             1.07000000
   B3             1.07000000
   B4             2.00000000
   B5             1.07000000
   B6             1.07000000
   B7             1.07000000
   A1           109.47120255
   A2           109.47123134
   A3           109.47120255
   A4           109.47120255
   A5           109.47120255
   A6           109.47123134
   D1           119.99999260
   D2          -120.00001480
   D3           179.88890000
   D4           -60.11108520
   D5            59.88890740
constant:
B 4 5 2.00000000

（4 5 原子就是两个碳原子）



计算过程中没有出错，但是在最后看输出文件的时候发现键长又回去了，不知道应该怎么输入，像上面那个输入文件，应该怎么改可以改对呢？
赞一下 回复此楼
9楼2011-03-12 19:32:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 23 个回答

qlihytc

新虫 (小有名气)
亍寞(金币+10): 谢谢回复，还有些许问题还请您不吝赐教！ 2011-03-12 15:04:29
我建议这样寻找：
找到一个主要的反应坐标，键长或者键角，固定这个反应坐标，然后popt，通过反应坐标的逐渐改变，得到一条反应的PES，再在最高点进行激发态优化
赞一下 回复此楼
2楼2011-03-12 14:31:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

亍寞

金虫 (小有名气)
引用回帖:
Originally posted by qlihytc at 2011-03-12 14:31:06:
我建议这样寻找：
找到一个主要的反应坐标，键长或者键角，固定这个反应坐标，然后popt，通过反应坐标的逐渐改变，得到一条反应的PES，再在最高点进行激发态优化

谢谢，现在您所说的主要的反应坐标已经找到了，是键长，但是就是再固定的时候除了点问题，固定之后计算，但是计算结果中这个希望固定的键长你却恢复了稳定态的键长，因为以前没有固定过键长进行计算，所以也不知道是不是固定方法有问题
这是部分输入文件：（z-matrix）
%mem=3500MB

%nprocshared=2

#p opt=(modredundant,maxcyc=200) freq b3lyp/genecp scf=tight
键说明部分。。。略
B1 4.900000

基组说明部分。。。略

请问这个是一定要用popt命令么？那么popt（参数）中的参数又应该怎么设置比较好？您所说的反应坐标的逐渐改变是用scan命令还是人为的固定不同键长?

十分感谢您的回复！
赞一下 回复此楼
3楼2011-03-12 15:04:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qlihytc

新虫 (小有名气)
不是SCAN，这样，给你篇文献看看吧
羟基负离子与一氟甲烷双分子亲核取代反应的量子化学研究
赞一下 回复此楼
4楼2011-03-12 16:55:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 23 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭