24小时热门版块排行榜    

Znn3bq.jpeg
北京石油化工学院2026年研究生招生接收调剂公告
查看: 2766  |  回复: 11

ludeng

新虫 (小有名气)

[求助] 固定一个键长的同时 扫描二面角?

如图所示,我想在C1-C2键长不变的情况下,扫描C5-C1-C2-N3二面角,主要目的是看 C1-C5和C2-N3的何种取向下,整体结构的能量最低。
请问 怎样才能固定一个键长的同时 扫描二面角?
回复此楼

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

过渡态计算

» 猜你喜欢

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

可惜不是你
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

070241078

金虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
gmy1990: 金币+2 2012-04-10 15:36:37
可以先优化出来一个构型,然后以该构型为基础进行势能面扫描,采用Z-矩阵定义变量。扫描时的关键字opt=modredundant,在输入文件的变量定义部分指定要扫描的坐标,如D 5 1 2 3 1.00(初始值) S 10(步数) 10.00(步长)
2楼2012-04-10 10:51:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

ludeng

新虫 (小有名气)

引用回帖:
2楼: Originally posted by 070241078 at 2012-04-10 10:51:22:
可以先优化出来一个构型,然后以该构型为基础进行势能面扫描,采用Z-矩阵定义变量。扫描时的关键字opt=modredundant,在输入文件的变量定义部分指定要扫描的坐标,如D 5 1 2 3 1.00(初始值) S 10(步数) 10.00(步 ...

那我C1-C2的键长该如何固定?我想C1-C2的键长在扫描过程中保持不变
可惜不是你
3楼2012-04-10 12:43:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

飘逸110

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
gmy1990: 金币+2 2012-04-10 15:36:33
建议采用刚性扫描势能面!柔性的势能扫描健长会改变的!070241078给你的就是柔性的扫描,每改变一次二面角度,相应的构型通过健长,键角的改变优化后得到一个最小的能量,不能满足你的要求的!
4楼2012-04-10 15:08:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ludeng

新虫 (小有名气)

引用回帖:
4楼: Originally posted by 飘逸110 at 2012-04-10 15:08:13:
建议采用刚性扫描势能面!柔性的势能扫描健长会改变的!070241078给你的就是柔性的扫描,每改变一次二面角度,相应的构型通过健长,键角的改变优化后得到一个最小的能量,不能满足你的要求的!

请问 该怎么设置呢? O(∩_∩)O谢谢~
#  b3lyp/3-21g scf=(maxcyc=400,conver=5) scan
0 1
C
C                  1                            2.26578598
C                  2              B2    1              A1
C                  3              B3    2              A2    1              D1
...
...
   D51         -138.49152090
   D52          106.36177546
   D53           58.89876470
   D54          177.81778366
   D55          -63.38028521

   D33 36 1 2 8         74.04638379  24 -15

C1-C2的键长我在上面用光了常量,然后在最后面写出了二面角的变化信息,但跑出来的结果貌似很不对...
可惜不是你
5楼2012-04-10 16:31:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ludeng

新虫 (小有名气)

请问 改怎么设置呢?
#  b3lyp/3-21g scf=(maxcyc=400,conver=5) scan
0 1
C
C                  1                            2.26578598
C                  2              B2    1              A1
C                  3              B3    2              A2    1              D1
...
...
   D51         -138.49152090
   D52          106.36177546
   D53           58.89876470
   D54          177.81778366
   D55          -63.38028521

   D33 36 1 2 8         74.04638379  24 -15

C1-C2的键长我在上面用光了常量,然后在最后面写出了二面角的变化信息,但跑出来的结果貌似很不对...


可惜不是你
6楼2012-04-10 16:32:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

飘逸110

木虫 (正式写手)

【答案】应助回帖

最笨的方法就是,改变一个角度计算一个能量!弄一堆输入文件,在linux计算,将计算的能量提取出来作图!个人见解。。。还望高手赐教!
7楼2012-04-10 18:51:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ldzh0531

木虫 (小有名气)

我用SCAN时也遇到了一个问题,就是扫描二面角的时候出来的结果只是一个原子动,我想要的结果是单键的旋转。请问怎么改正。
8楼2012-04-10 20:51:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

飘逸110

木虫 (正式写手)

【答案】应助回帖

引用回帖:
8楼: Originally posted by ldzh0531 at 2012-04-10 20:51:47:
我用SCAN时也遇到了一个问题,就是扫描二面角的时候出来的结果只是一个原子动,我想要的结果是单键的旋转。请问怎么改正。

呵呵。。
我就是这样出错了!没办法就人为的改角度算能量咯!
同道中人!
9楼2012-04-10 21:32:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ludeng

新虫 (小有名气)

引用回帖:
8楼: Originally posted by ldzh0531 at 2012-04-10 20:51:47:
我用SCAN时也遇到了一个问题,就是扫描二面角的时候出来的结果只是一个原子动,我想要的结果是单键的旋转。请问怎么改正。

# opt=modredundant ....
...
...
...

D 6 5 4 3 30.0(初始角度) s(scan简写) 5(步数) 10.0(步长)

比如你要扫描二面角6-5-4-3,在z-矩阵中将要扫描的二面角那一行的键长、键角、二面角都用数值代替,然后将下面变量信息中的数值删掉,就可以了
可惜不是你
10楼2012-04-10 21:37:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ludeng 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 363求调剂 +8 zh096 2026-04-04 8/400 2026-04-07 14:53 by 1104338198
[考研] 材料工程专业日语生求调剂 +7 111623 2026-04-07 7/350 2026-04-07 14:33 by 诗与自由
[考研] 085600材料与化工专硕329 求调剂 +16 额cc 2026-04-06 17/850 2026-04-07 12:45 by rl1980
[考研] 求调剂 +7 wwwwabcde 2026-04-07 7/350 2026-04-07 12:32 by dongzh2009
[考研] 0703化学调剂325分 +12 15771691647 2026-04-04 13/650 2026-04-06 12:00 by lijunpoly
[考研] 315求调剂 +14 欣喜777 2026-04-04 15/750 2026-04-06 06:25 by houyaoxu
[考研] 322求调剂 +3 嗯哼哼恒 2026-04-05 3/150 2026-04-05 19:52 by nepu_uu
[考研] 材料专硕(0856) 339分求调剂 +10 哈哈哈鹅哈哈哈 2026-04-04 10/500 2026-04-05 18:51 by 蓝云思雨
[考研] 282求调剂 +7 aaa车辆 2026-04-02 11/550 2026-04-05 17:24 by yulian1987
[考研] 材料专硕 调剂 +14 CXN123456 2026-04-03 14/700 2026-04-05 17:18 by Hdyxbekcb
[考研] 材料化工306分找合适调剂 +14 沧海轻舟e 2026-04-04 14/700 2026-04-05 09:53 by 朱云虎202
[考研] 调剂 +9 19945159693 2026-04-03 10/500 2026-04-04 20:16 by dongzh2009
[考研] 309分085801求调剂 +11 MY_angel 2026-03-31 11/550 2026-04-04 19:11 by 蓝云思雨
[考研] 材料科学与工程考研 +10 拯救皮特托先生 2026-04-02 10/500 2026-04-03 23:57 by userper
[考研] 311求调剂 +20 zchqwer 2026-04-01 22/1100 2026-04-03 22:09 by lglzsd
[考研] 考研求调剂 +3 木心想继续深造 2026-04-03 3/150 2026-04-03 21:56 by 啵啵啵0119
[考研] 求调剂 +3 心想事成可 2026-04-03 3/150 2026-04-03 11:22 by wangjy2002
[考研] 08工科求调剂290分 +5 1314捧花 2026-04-02 8/400 2026-04-02 13:16 by 乔哒哒哒
[考研] 349求调剂 +6 吃的不少 2026-04-01 6/300 2026-04-01 17:55 by JYD2011
[考研] 复试调剂 +7 双马尾痞老板2 2026-03-31 7/350 2026-03-31 19:49 by Dyhoer
信息提示
请填处理意见