计算化学公社

 找回密码 Forget password
 注册 Register
Views: 68980|回复 Reply: 100
打印 Print 上一主题 Last thread 下一主题 Next thread

[新手求助] ONIOM怎么对L层指定原子电荷

[复制链接 Copy URL]

586

帖子

0

威望

4191

eV
积分
4777

Level 6 (一方通行)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 15whyhq 于 2020-1-11 14:21 编辑

我想问:这是我之前用uff算低层,用DFT算H层,L层没有指定电荷的输入输出文件,我现在要指定L层的原子电荷,我应该怎么做


c-2m-1-OH-TS1-oniom-ts-2.gjf

12.37 KB, 下载次数 Times of downloads: 62

c-2m-1-OH-TS1-oniom-ts-2.log

6.53 MB, 下载次数 Times of downloads: 28

8

帖子

0

威望

222

eV
积分
230

Level 3 能力者

102#
发表于 Post on 2017-1-5 18:21:10 | 只看该作者 Only view this author
这个帖子不光体现了sob老师的学识渊博,更体现了他的大师风范,有问必答,有耐心。

1294

帖子

0

威望

6236

eV
积分
7530

Level 6 (一方通行)

101#
发表于 Post on 2016-12-23 10:45:14 | 只看该作者 Only view this author
呼呼呼,终于看完楼主的这个帖子了。我是在群里咨询了ONIOM的问题,sob老师给了回复指点,不太理解,然后查资料,看到了这个帖子。
看完了,有帮助,sob老师回复很仔细,关键是很有耐心。但是,看的我都没有耐心了,忍不住吐槽一下。
1:师傅领到门,巧妙在个人。遇到问题,问别人很好,但是不能这样无休止的问下去,老师是给你指点,解决问题还是需要自己来解决的。
2:楼主的帖子,除了楼主和sob,还有一个回复的,就是希望楼主能有这些问的时间,来查查资料,学习一下。但是,楼主没有听取,还在问。
3:最后不得不说,sob老师真的很有耐心,真的真的。我都看得不耐烦了,哪有这样问的?感觉心里憋火,所以,我才来吐槽一下哈。

看完这个帖子收获很多,不仅是知识,还有一些其他的。个人见解。

评分 Rate

参与人数
Participants 3
eV +7 收起 理由
Reason
rainseason + 3
bing0 + 2
sobereva + 2

查看全部评分 View all ratings

https://www.x-mol.com/groups/fan_jianzhong

strive for greatness

586

帖子

0

威望

4191

eV
积分
4777

Level 6 (一方通行)

100#
 楼主 Author| 发表于 Post on 2016-11-25 11:17:10 | 只看该作者 Only view this author
固定高层优化低层的时候优化了十几次,一直出现这个错误,怎么解决?谢谢
Error imposing constraints                                                         
Error termination via Lnk1e in /vol6/home/hn_zjl/software/Gaussion09/g09B01/l103.exe

2-10-for-th-hf-4-mm-10.zip

493.96 KB, 下载次数 Times of downloads: 2

5万

帖子

99

威望

5万

eV
积分
112496

管理员

公社社长

99#
发表于 Post on 2016-11-16 19:40:23 | 只看该作者 Only view this author
15whyhq 发表于 2016-11-16 19:29
我用oniom方法进行优化构型以及寻找过渡态的关键词如图所示,帮忙看下我的关键词有没有什么问题,谢谢!


优化极小点的opt里写noeigentest完全多余
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

586

帖子

0

威望

4191

eV
积分
4777

Level 6 (一方通行)

98#
 楼主 Author| 发表于 Post on 2016-11-16 19:29:19 | 只看该作者 Only view this author
我用oniom方法进行优化构型以及寻找过渡态的关键词如图所示,帮忙看下我的关键词有没有什么问题,谢谢!

opt.png (2.67 KB, 下载次数 Times of downloads: 70)

opt.png

ts.png (3.23 KB, 下载次数 Times of downloads: 87)

ts.png

586

帖子

0

威望

4191

eV
积分
4777

Level 6 (一方通行)

97#
 楼主 Author| 发表于 Post on 2016-11-1 21:25:27 | 只看该作者 Only view this author
老师,我想问下,ONIOM计算时低层使用amber力场时,指定原子电荷时手动给相应的原子指定原子电荷,我不知道这原子电荷可不可以使用MD模拟的时用resp电荷模型中的HF/6-31G计算得到的原子电荷?
是不是不同的方法计算得到的原子电荷不一样?我用ONIOM计算高层的方法是b3pw91/6-31g(d,p)与计算原子电荷的方法不一样,所以我想知道怎么解决低层用amber力场时怎么得到可以指定的原子电荷?我用oniom计算的初始构型是通过MD模拟得到的。
谢谢!

5万

帖子

99

威望

5万

eV
积分
112496

管理员

公社社长

96#
发表于 Post on 2016-10-31 11:40:25 | 只看该作者 Only view this author
15whyhq 发表于 2016-10-31 09:57
geom=connectivity   在高斯说明书里面我只看到这个关键词是用来表达原子之间的连接关系。 我不明白低层用 ...

用力场计算必须知道原子之间是怎么相连的,否则没法确定原子类型
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

586

帖子

0

威望

4191

eV
积分
4777

Level 6 (一方通行)

95#
 楼主 Author| 发表于 Post on 2016-10-31 09:57:55 | 只看该作者 Only view this author
geom=connectivity   在高斯说明书里面我只看到这个关键词是用来表达原子之间的连接关系。 我不明白低层用力场时为什么要定义连接关系?
老师可以帮忙解释下吗?谢谢!

5万

帖子

99

威望

5万

eV
积分
112496

管理员

公社社长

94#
发表于 Post on 2016-10-30 08:45:52 | 只看该作者 Only view this author
15whyhq 发表于 2016-10-29 20:54
我这个输入文件哪有问题?为什么在GaussView里打不开?

诸如H-H_ _0.095897根本不规矩
元素、原子类型、电荷之间用横线分割。不能有莫名其妙的空格。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

5万

帖子

99

威望

5万

eV
积分
112496

管理员

公社社长

93#
发表于 Post on 2016-10-30 08:42:40 | 只看该作者 Only view this author
15whyhq 发表于 2016-10-25 17:18
当有化学反应发生时,QM 区电荷分布会发生变化,在处理静电相互作用中,用固定原子点电荷描述是不准确的。
...

就是指使用的原子电荷不随反应过程的构型/构象的改变而改变,一套电荷从头用到尾。然而原子电荷必定依赖于化学环境,所以固定点电荷有时会有问题,甚至可能严重。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

586

帖子

0

威望

4191

eV
积分
4777

Level 6 (一方通行)

92#
 楼主 Author| 发表于 Post on 2016-10-29 20:54:10 | 只看该作者 Only view this author
我这个输入文件哪有问题?为什么在GaussView里打不开?

2-16-for-n-4-amber.gjf

42.5 KB, 下载次数 Times of downloads: 4

586

帖子

0

威望

4191

eV
积分
4777

Level 6 (一方通行)

91#
 楼主 Author| 发表于 Post on 2016-10-25 17:18:37 | 只看该作者 Only view this author
当有化学反应发生时,QM 区电荷分布会发生变化,在处理静电相互作用中,用固定原子点电荷描述是不准确的。
我想知道固定原子点电荷是什么意思?怎么做的?
谢谢!

586

帖子

0

威望

4191

eV
积分
4777

Level 6 (一方通行)

90#
 楼主 Author| 发表于 Post on 2016-10-19 09:51:07 | 只看该作者 Only view this author
还有个问题低层用Dreiding方法时,所有的输入文件是不是跟低层用UFF方法时的一样?低层用Dreiding方法得指定原子电荷,是不是也是用QEq电荷?
谢谢!

586

帖子

0

威望

4191

eV
积分
4777

Level 6 (一方通行)

89#
 楼主 Author| 发表于 Post on 2016-10-19 09:09:19 | 只看该作者 Only view this author
为什么低层用UFF时没有使用QEq与使用QEq优化完后构型相差那么多,用QEq后环氧化物与离子液体相距那么远,环氧化物直接从离子液体中间跑到边上去了还隔那么远。两个作业都正常结束。谢谢!

po-7-uff-1.zip

1.44 MB, 下载次数 Times of downloads: 3

po-7-uff-QEq-1.zip

1.57 MB, 下载次数 Times of downloads: 1

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-27 06:43 , Processed in 0.203552 second(s), 26 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list