计算化学公社

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

[综合交流] LJ参数ε和σ的量子化学计算

[复制链接 Copy URL]

13

帖子

0

威望

92

eV
积分
105

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
请问分子模拟力场中的非键相互作用L-J函数的εσ值是否可以通过量化计算进行拟合?如果可以应该怎样拟合,或者说通过量化计算什么物理量进行拟合?请教计算思路,或者相关文献,谢谢!

6万

帖子

99

威望

6万

eV
积分
125181

管理员

公社社长

2#
发表于 Post on 2017-10-9 22:13:38 | 只看该作者 Only view this author
看具体体系。
比如就是简单的Ar-Ar的双原子LJ势,直接用比如CCSD(T)/aug-cc-pVQZ做个势能面扫描,然后用Origin根据LJ公式拟合那两个参数就完了。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

13

帖子

0

威望

92

eV
积分
105

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2017-10-10 09:54:29 | 只看该作者 Only view this author
sobereva 发表于 2017-10-9 22:13
看具体体系。
比如就是简单的Ar-Ar的双原子LJ势,直接用比如CCSD(T)/aug-cc-pVQZ做个势能面扫描,然后用Or ...

老师,如果想计算制冷剂分子(如R134a,CH2FCF3)中原子的ε和σ值,是否可以用量化的方法计算呢?

1171

帖子

7

威望

6862

eV
积分
8173

Level 6 (一方通行)

4#
发表于 Post on 2017-10-10 11:21:03 | 只看该作者 Only view this author
nxh9681 发表于 2017-10-10 09:54
老师,如果想计算制冷剂分子(如R134a,CH2FCF3)中原子的ε和σ值,是否可以用量化的方法计算呢?

一般的做法都是固定力场(比如cgenff,gaff)自己生成的vdw参数,通过量化计算拟合静电参数的

6万

帖子

99

威望

6万

eV
积分
125181

管理员

公社社长

5#
发表于 Post on 2017-10-10 13:57:04 | 只看该作者 Only view this author
nxh9681 发表于 2017-10-10 09:54
老师,如果想计算制冷剂分子(如R134a,CH2FCF3)中原子的ε和σ值,是否可以用量化的方法计算呢?

这一般不是简单量化方法就能算的,而且算也不好算,因为同时牵扯多种原子对,还牵扯静电作用。建议直接用普适型力场里现成的原子LJ参数
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

13

帖子

0

威望

92

eV
积分
105

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2017-10-10 16:21:14 | 只看该作者 Only view this author
sobereva 发表于 2017-10-10 13:57
这一般不是简单量化方法就能算的,而且算也不好算,因为同时牵扯多种原子对,还牵扯静电作用。建议直接用 ...

好的,明白了,多谢老师

13

帖子

0

威望

92

eV
积分
105

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2017-10-10 16:21:46 | 只看该作者 Only view this author
fhh2626 发表于 2017-10-10 11:21
一般的做法都是固定力场(比如cgenff,gaff)自己生成的vdw参数,通过量化计算拟合静电参数的

嗯嗯,多谢指导!

5

帖子

0

威望

15

eV
积分
20

Level 1 能力者

8#
发表于 Post on 2020-12-10 14:29:24 | 只看该作者 Only view this author
sobereva 发表于 2017-10-10 13:57
这一般不是简单量化方法就能算的,而且算也不好算,因为同时牵扯多种原子对,还牵扯静电作用。建议直接用 ...

请问这个普适型力场的原子LJ参数去哪里查

1万

帖子

0

威望

9892

eV
积分
22140

Level 6 (一方通行)

9#
发表于 Post on 2020-12-10 17:05:26 | 只看该作者 Only view this author
weaey 发表于 2020-12-10 14:29
请问这个普适型力场的原子LJ参数去哪里查

比如UFF力场的话,可以查UFF力场的原始文献
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

6万

帖子

99

威望

6万

eV
积分
125181

管理员

公社社长

10#
发表于 Post on 2020-12-10 19:18:42 | 只看该作者 Only view this author
weaey 发表于 2020-12-10 14:29
请问这个普适型力场的原子LJ参数去哪里查

UFF原文J. Am. Chem. SOC. 1992, 114, 10024-10039的表1里就是。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

3

帖子

0

威望

83

eV
积分
86

Level 2 能力者

11#
发表于 Post on 2023-6-28 15:17:06 | 只看该作者 Only view this author
sobereva 发表于 2020-12-10 19:18
UFF原文J. Am. Chem. SOC. 1992, 114, 10024-10039的表1里就是。

请问sob老师,粗粒化基团的σ和ε怎么得到

6万

帖子

99

威望

6万

eV
积分
125181

管理员

公社社长

12#
发表于 Post on 2023-6-29 00:04:29 | 只看该作者 Only view this author
Ecust-zy 发表于 2023-6-28 15:17
请问sob老师,粗粒化基团的σ和ε怎么得到

不同力场有不同做法。比如模拟凝聚相性质,通过让误差最小化来调节σ和ε。看你要用的粗粒化力场或者类似目的的粗粒化力场原文里的说明
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

4

帖子

0

威望

17

eV
积分
21

Level 1 能力者

13#
发表于 Post on 2023-10-17 12:10:03 | 只看该作者 Only view this author
sobereva 发表于 2023-6-29 00:04
不同力场有不同做法。比如模拟凝聚相性质,通过让误差最小化来调节σ和ε。看你要用的粗粒化力场或者类似 ...

老师,我是电气方向准备做一个气体C4F7N的comsol仿真,需要这两个参数,需要怎么计算啊,气体在高压下的流注放电模拟

1万

帖子

0

威望

9892

eV
积分
22140

Level 6 (一方通行)

14#
发表于 Post on 2023-10-17 18:12:04 | 只看该作者 Only view this author
吞吞12138 发表于 2023-10-17 05:10
老师,我是电气方向准备做一个气体C4F7N的comsol仿真,需要这两个参数,需要怎么计算啊,气体在高压下的 ...

你要用σ和ε描述C4F7N分子之间的总相互作用,还是只描述其中的范德华作用部分?是否打算描述Keesom项和Debye项,还是只描述London色散项和交换排斥项?
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

4

帖子

0

威望

17

eV
积分
21

Level 1 能力者

15#
发表于 Post on 2023-10-17 21:50:19 | 只看该作者 Only view this author
wzkchem5 发表于 2023-10-17 18:12
你要用σ和ε描述C4F7N分子之间的总相互作用,还是只描述其中的范德华作用部分?是否打算描述Keesom项和D ...

老师请问这两个参数对于某个分子来说不是常数吗,仿真中只是因为材料库中没有这种气体,需要输入这两个参数,来自定义这个气体的基本物理量。例如软件中材料库中自带的氩气Ar这个两个参数分别是σ:0.365纳米和σ/Kb:178k。这两个参数不是常数吗对某一确定分子。其他的气体大部分都能从材料库中直接选中,这个比较冷门。

本版积分规则 Credits rule

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

GMT+8, 2026-2-26 04:26 , Processed in 0.182391 second(s), 26 queries , Gzip On.

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