| 查看: 1023 | 回复: 11 | |||
| 本帖产生 1 个 QC强帖 ，点击这里进行查看 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
renren123123金虫 (正式写手)
|
[求助]
分子优化求助
|
||
|
# opt freq b3lyp/6-31++g(2d,2p) nosymm geom=connectivity A3 1 1 C -0.97218474 -0.04841131 0.00000164 C 0.74654581 0.02078778 0.00004819 N 0.00908182 1.08732390 0.00003223 N -0.08521263 -1.08127384 -0.00003048 N -2.25808564 0.02767855 -0.00010128 H -2.82238882 -0.81322348 0.00033844 H -2.71249356 0.93088470 -0.00019287 1 3 2.0 4 1.0 5 1.0 2 3 1.0 4 2.0 3 4 5 6 1.0 7 1.0 6 7 这样优化始终优化不出我想要的如opt-model.jpg图结构，不自恰，麻烦哪位高手用高斯03或09帮我找到正确优化的原因，如果有时间，还请优化出来贴上结果，谢谢！  opt-model.JPG |
回复此楼» 猜你喜欢
书籍求助：汽车市场营销理论与实务（电子版）——章小平
已经有0人回复
-
东方理-中科大联合博士生招聘
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有166人回复
-
钴酸锂半电池小倍率容量上不去
已经有1人回复
-
吉林大学材料物理本科生求问调剂信息
已经有23人回复
-
2026第二届光电子与半导体器件前沿技术研讨会——光电子赋能·半导体创芯！
已经有1人回复
-
光电子赋能·半导体创芯！
已经有1人回复
-
光电子赋能·半导体创芯！
已经有0人回复
-
QE计算电声耦合的时候报错Error in routine lambda (100)wrong or too many modes
已经有2人回复
-
分子动力学模拟合作需求
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
huilaoshu999
金虫 (著名写手)
- 应助: 18 (小学生)
- 金币: 941.9
- 散金: 5
- 红花: 2
- 帖子: 1734
- 在线: 304.8小时
- 虫号: 228200
- 注册: 2006-03-24
- 专业: 有机分子功能材料化学
【答案】应助回帖
|
建议： 1、删除命令中的geom=connectivity，同时将笛卡尔坐标下的部分全部删除 # opt freq b3lyp/6-31++g(2d,2p) nosymm geom=connectivity A3 1 1 C -0.97218474 -0.04841131 0.00000164 C 0.74654581 0.02078778 0.00004819 N 0.00908182 1.08732390 0.00003223 N -0.08521263 -1.08127384 -0.00003048 N -2.25808564 0.02767855 -0.00010128 H -2.82238882 -0.81322348 0.00033844 H -2.71249356 0.93088470 -0.00019287 1 3 2.0 4 1.0 5 1.0 2 3 1.0 4 2.0 3 4 5 6 1.0 7 1.0 6 7 2、你的计算体系中一个碳原子是六电子的，自旋多重度是1还是3？是否要用Ub3lyp/6-31++g(2d,2p)？ 3、用大基组计算最好先用小机组优化，可能省时，而且效果好些 |
12楼2012-12-03 16:37:39
lxhlxh052c
木虫 (著名写手)
小木虫之无虫
- 应助: 21 (小学生)
- 贵宾: 0.025
- 金币: 1153.8
- 散金: 2681
- 红花: 26
- 沙发: 10
- 帖子: 1914
- 在线: 617.9小时
- 虫号: 759932
- 注册: 2009-04-29
- 性别: GG
- 专业: 化学动力学
2楼2012-11-27 10:47:07
gauss98
禁虫 (正式写手)
- QC强帖: 1
- 应助: 112 (高中生)
- 金币: 1123.5
- 红花: 17
- 帖子: 317
- 在线: 1060.8小时
- 虫号: 2067376
- 注册: 2012-10-17
- 性别: GG
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
renren123123: 回帖置顶 2012-11-27 23:12:43
cenwanglai: 金币+5, QC强帖+1, 感谢热心回帖~ 2012-11-28 10:21:31
感谢参与，应助指数 +1
renren123123: 回帖置顶 2012-11-27 23:12:43
cenwanglai: 金币+5, QC强帖+1, 感谢热心回帖~ 2012-11-28 10:21:31
|
你原来的结构两个碳原子离得太近了。拉远一点形成一个接近正方形就可以了 我给你调整的初始坐标输入文件： % nproc=8 # opt freq b3lyp/6-31++g(2d,2p) nosymm geom=connectivity A3 1 1 C -1.13450913 -0.01127887 -0.00000220 C 0.95631432 0.02923363 0.00005361 N -0.03836219 1.04174624 0.00003058 N -0.08521300 -1.08127400 -0.00003000 N -2.42041013 0.06481113 -0.00010520 H -2.98471313 -0.77609087 0.00033380 H -2.87481813 0.96801713 -0.00019720 1 3 1.0 4 1.5 5 2.0 2 3 2.0 4 1.5 3 4 5 6 1.0 7 1.0 6 7 优化结束后的结果算频率的坐标： C,0,-1.0994478143,0.0776653719,-0.1739757193 C,0,0.7433080732,-0.0902050705,0.5030218806 N,0,-0.0435391597,0.9003683894,-0.3201383759 N,0,-0.0614862126,-0.8627745918,-0.2995524251 N,0,-2.3634015779,0.0714470718,0.0499234862 H,0,-2.8801374252,-0.8041990431,0.0808993623 H,0,-2.8770072735,0.9428622623,0.1599051813 注意结果是个局部能量极小点，能量很高，比线性构型能量高很多！  help1.JPG |
3楼2012-11-27 11:02:53
renren123123
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 924.3
- 散金: 100
- 帖子: 599
- 在线: 81.5小时
- 虫号: 303284
- 注册: 2006-12-03
- 专业: physics
4楼2012-11-27 23:10:08