计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: diaok
打印 Print 上一主题 Last thread 下一主题 Next thread

[综合交流] 求教高分子电荷拟合

[复制链接 Copy URL]

384

帖子

0

威望

3364

eV
积分
3748

Level 5 (御坂)

16#
发表于 Post on 2019-2-26 05:20:21 | 只看该作者 Only view this author
本帖最后由 yihanxu 于 2019-2-26 02:20 编辑
sobereva 发表于 2019-2-23 08:51
用Multiwfn算RESP电荷多好,下文就有算聚合物单体片段的原子电荷的例子,比用MMFF94质量好得多,计算极其 ...

目的是算高分子的RESP电荷:用gview画了一个高分子三个重复单元的结构(两端用氢原子封闭了),请问老师:之后可以直接优化、计算电荷吗?还是要多考虑一些构象?
附件是正在优化中的模型,可以请老师帮忙看一看吗:可以直接用中间的一个重复单元算高分子的RESP电荷、两端的重复单元(被氢封闭的)分别算高分子链两端的RESP电荷吗?
谢谢老师。

PANI_RESP.rar

318 KB, 下载次数 Times of downloads: 3

正在优化中的模型

repeating_unit.gjf

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

重复单元的结构

live each day as if it was your last

6万

帖子

99

威望

6万

eV
积分
125156

管理员

公社社长

17#
发表于 Post on 2019-2-26 16:46:51 | 只看该作者 Only view this author
yihanxu 发表于 2019-2-26 05:20
目的是算高分子的RESP电荷:用gview画了一个高分子三个重复单元的结构(两端用氢原子封闭了),请问老师 ...

1 可以。多考虑几种构象有益,但不是必须的
2 可以

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
yihanxu + 5 谢谢

查看全部评分 View all ratings

北京科音自然科学研究中心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

384

帖子

0

威望

3364

eV
积分
3748

Level 5 (御坂)

18#
发表于 Post on 2019-2-27 13:56:38 | 只看该作者 Only view this author
本帖最后由 yihanxu 于 2019-2-27 00:06 编辑
sobereva 发表于 2019-2-26 02:46
1 可以。多考虑几种构象有益,但不是必须的
2 可以

老师好,我分别尝试了两种约束方法计算了RESP电荷:
1. 左侧链末端电荷为0,右侧链末端电荷为0,重复单元电荷为0,得到RMSE:    0.002266   RRMSE:    0.193826 。但是左侧链末端带了电荷3E-6 。2. 链两端整体电荷为0,重复单元电荷为0,得到RMSE:    0.001685   RRMSE:    0.144163 。但是两个被约束的整体最后都带电荷-1E-6 。
请问这两种方法中是否有一种的结果可以用呢?谢谢。

chg_2.txt

3.73 KB, 下载次数 Times of downloads: 1

方法二,Multiwfn窗口拷下来的计算结果

chg_1.txt

2.83 KB, 下载次数 Times of downloads: 1

方法一,Multiwfn窗口拷下来的计算结果

Worksheet.xlsx

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

电荷加和

PANI_RESP_1.chg

3.59 KB, 下载次数 Times of downloads: 1

方法一的计算结果

PANI_RESP_2.chg

3.59 KB, 下载次数 Times of downloads: 1

方法二的计算结果

live each day as if it was your last

6万

帖子

99

威望

6万

eV
积分
125156

管理员

公社社长

19#
发表于 Post on 2019-2-27 16:02:56 | 只看该作者 Only view this author
yihanxu 发表于 2019-2-27 13:56
老师好,我分别尝试了两种约束方法计算了RESP电荷:
1. 左侧链末端电荷为0,右侧链末端电荷为0,重复单 ...

用1,比较有物理意义
1E-6程度的数值可以认为是数值误差,可忽略不计

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
yihanxu + 5 谢谢

查看全部评分 View all ratings

北京科音自然科学研究中心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

43

帖子

0

威望

207

eV
积分
250

Level 3 能力者

20#
发表于 Post on 2023-11-1 15:20:10 | 只看该作者 Only view this author
sobereva 发表于 2016-1-11 20:01
就用AM1-BCC吧,质量不错,不逊于RESP,算得很快。重复单元肯定是需要人工调整一下,这一点的任意性免不了
...

老师,对于普通线性聚合物分子,使用sobtop自动分配MMFF94电荷做分子动力学模拟,这样做的话在目前计算化学领域中发文章会受到审稿人质疑吗?(我看您在sobtop中说RESP10分,MMFF94 7分,但是还是不确定发文章受不受影响) 自己有点懒...体系又多不想一个个算,想问问老师意见。

6万

帖子

99

威望

6万

eV
积分
125156

管理员

公社社长

21#
发表于 Post on 2023-11-1 15:44:35 | 只看该作者 Only view this author
liu840397562 发表于 2023-11-1 15:20
老师,对于普通线性聚合物分子,使用sobtop自动分配MMFF94电荷做分子动力学模拟,这样做的话在目前计算化 ...

有一定概率会被审稿人吐槽
北京科音自然科学研究中心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

43

帖子

0

威望

207

eV
积分
250

Level 3 能力者

22#
发表于 Post on 2023-11-1 17:01:52 | 只看该作者 Only view this author
sobereva 发表于 2023-11-1 15:44
有一定概率会被审稿人吐槽

那还是老实算一下吧 谢谢sob老师

43

帖子

0

威望

207

eV
积分
250

Level 3 能力者

23#
发表于 Post on 2023-11-7 20:53:08 | 只看该作者 Only view this author
sobereva 发表于 2016-1-11 20:01
就用AM1-BCC吧,质量不错,不逊于RESP,算得很快。重复单元肯定是需要人工调整一下,这一点的任意性免不了
...

sob老师,我看了您在物理化学学报关于各种电荷的比较,产生一些问题想咨询您。我看到其中MMFF94电荷在偶极矩计算方面,其与HF/6-31G*的期望值的MUE还不错,但是与实际实验值的MUE就较大了。
这会不会导致使用MMFF94电荷的聚合物体系,在使用gromacs的gmx current计算静态介电常数时候与实验值偏差较大呢?

6万

帖子

99

威望

6万

eV
积分
125156

管理员

公社社长

24#
发表于 Post on 2023-11-8 10:58:34 | 只看该作者 Only view this author
liu840397562 发表于 2023-11-7 20:53
sob老师,我看了您在物理化学学报关于各种电荷的比较,产生一些问题想咨询您。我看到其中MMFF94电荷在偶 ...

MMFF94只能算是将就用的程度
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2026-2-24 20:30 , Processed in 0.215356 second(s), 23 queries , Gzip On.

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