计算化学公社

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

[综合交流] 简单量子化学问题答疑专帖(旧版,已关帖)

   关闭 [复制链接 Copy URL]

17

帖子

0

威望

389

eV
积分
406

Level 3 能力者

8446#
发表于 Post on 2019-2-26 21:47:28 | 只看该作者 Only view this author
对于阳离子参与的弱相互作用(非氢键),原子个数大概有40个,前两周期元素非金属体系,计算相互作用能,MP2/aug-cc-pVTZ理论水平是否合适?我想用DFT方法,也不用这么多的弥散函数,老大想用这个理论水平计算,特来请教一下大家,谢谢

17

帖子

0

威望

795

eV
积分
812

Level 4 (黑子)

8447#
发表于 Post on 2019-2-26 22:11:11 | 只看该作者 Only view this author
sobereva 发表于 2019-2-26 19:04
这是激发态与基态间的密度差,和电场没有丝毫关系,你说的那是电场对基态偶极矩的影响

sob老师,那我如何才能得到deltamu呢?拜托啦,最近计算迫切需要deltamu的值。

6万

帖子

99

威望

6万

eV
积分
125156

管理员

公社社长

8448#
 楼主 Author| 发表于 Post on 2019-2-26 22:58:20 | 只看该作者 Only view this author
heheda 发表于 2019-2-26 22:11
sob老师,那我如何才能得到deltamu呢?拜托啦,最近计算迫切需要deltamu的值。

基态的偶极矩直接算单点就有
激发态的偶极矩做电子激发计算加上density就有
有不清楚的直接论坛首页google框搜索
北京科音自然科学研究中心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

6万

帖子

99

威望

6万

eV
积分
125156

管理员

公社社长

8449#
 楼主 Author| 发表于 Post on 2019-2-26 22:59:19 | 只看该作者 Only view this author
qingqing1223 发表于 2019-2-26 21:47
对于阳离子参与的弱相互作用(非氢键),原子个数大概有40个,前两周期元素非金属体系,计算相互作用能,MP ...

MP2算弱相互作用在现在来看很垃圾(氢键除外),在精度上被计算耗时基本相同的双杂化泛函完爆,用MP2计算氢键以外的弱相互作用是严重不关心理论方法进展的表现。就连比MP2便宜得多的B3LYP-D3(BJ)、M06-2X都比它强得多得多,仔细看下文的计算级别推荐,别再试图用过时的方法计算
乱谈DFT-D
http://sobereva.com/83
北京科音自然科学研究中心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

17

帖子

0

威望

795

eV
积分
812

Level 4 (黑子)

8450#
发表于 Post on 2019-2-26 23:10:57 | 只看该作者 Only view this author
sobereva 发表于 2019-2-26 22:58
基态的偶极矩直接算单点就有
激发态的偶极矩做电子激发计算加上density就有
有不清楚的直接论坛首页goo ...

额,sob老师,我问的是激发态与基态的偶极矩之差的求解方法,不是密度差吧。

6万

帖子

99

威望

6万

eV
积分
125156

管理员

公社社长

8451#
 楼主 Author| 发表于 Post on 2019-2-27 00:12:24 | 只看该作者 Only view this author
heheda 发表于 2019-2-26 23:10
额,sob老师,我问的是激发态与基态的偶极矩之差的求解方法,不是密度差吧。



基态的偶极矩直接算单点就有
激发态的偶极矩做电子激发计算加上density就有

两个数的减法会不会?

不知道density是啥把下文看100遍
Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
http://sobereva.com/314http://bbs.keinsci.com/thread-2413-1-1.html

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

17

帖子

0

威望

795

eV
积分
812

Level 4 (黑子)

8452#
发表于 Post on 2019-2-27 07:00:30 | 只看该作者 Only view this author

sob老师,不好意思,我刚问完就想明白了,但是不知道怎么删除。谢谢您啦,那么晚还在帮助我们解答各种问题。

17

帖子

0

威望

389

eV
积分
406

Level 3 能力者

8453#
发表于 Post on 2019-2-27 11:45:20 | 只看该作者 Only view this author
sobereva 发表于 2019-2-26 22:59
MP2算弱相互作用在现在来看很垃圾(氢键除外),在精度上被计算耗时基本相同的双杂化泛函完爆,用MP2计算 ...

谢谢社长,对于像我这种带正电荷的体系,计算相互作用能的时候用PWRB95方法,基组不加弥散函数合适吗?

6万

帖子

99

威望

6万

eV
积分
125156

管理员

公社社长

8454#
 楼主 Author| 发表于 Post on 2019-2-27 14:55:52 | 只看该作者 Only view this author
qingqing1223 发表于 2019-2-27 11:45
谢谢社长,对于像我这种带正电荷的体系,计算相互作用能的时候用PWRB95方法,基组不加弥散函数合适吗?

PWRB95是垃圾泛函,没有任何使用价值

如果是弱相互作用,计算能量的时候弥散必须加
北京科音自然科学研究中心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

188

帖子

0

威望

922

eV
积分
1110

Level 4 (黑子)

8455#
发表于 Post on 2019-2-27 19:37:13 | 只看该作者 Only view this author

收敛

Item               Value     Threshold  Converged?
Maximum Force            0.000006     0.000450     YES
RMS     Force            0.000002     0.000300     YES
Maximum Displacement     0.027817     0.001800     NO
RMS     Displacement     0.009253     0.001200     NO
opt=verytight 优化收敛得到四个YES,但是频率还是两个YES,这样还要继续优化吗?还是可以接受了

6万

帖子

99

威望

6万

eV
积分
125156

管理员

公社社长

8456#
 楼主 Author| 发表于 Post on 2019-2-27 20:53:38 | 只看该作者 Only view this author
liuyjhx 发表于 2019-2-27 19:37
Item               Value     Threshold  Converged?
Maximum Force            0.000006  ...

求助帖必须在帖子标题明确体现出“提问”、“求助”及具体内容(见http://bbs.keinsci.com/thread-9348-1-1.html),不要就写个“收敛”。过于简单、常见的问题直接在答疑专贴里问,不要单独发帖(恰当使用论坛首页google框)

反复阅读
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278http://bbs.keinsci.com/thread-633-1-1.html
北京科音自然科学研究中心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

188

帖子

0

威望

922

eV
积分
1110

Level 4 (黑子)

8457#
发表于 Post on 2019-2-28 19:16:49 | 只看该作者 Only view this author

低基组改为高基组怎么读取低基组的chk?

低基组命令
%chk=TS01.chk
%nprocshared=16
%mem=20GB
#p opt=(calcfc,ts,noeigentest) hf/3-21g freq

现在想改为高基组找,网上查是这样说的:“当用低的基组优化完直接读chk,ts+rcfc进一步用大基组优化”

现在不知道应该怎么写命令,尤其不知道直接读chk应该怎么写?

6万

帖子

99

威望

6万

eV
积分
125156

管理员

公社社长

8458#
 楼主 Author| 发表于 Post on 2019-2-28 19:34:07 | 只看该作者 Only view this author
liuyjhx 发表于 2019-2-28 19:16
低基组命令
%chk=TS01.chk
%nprocshared=16

用hf/3-21g计算根本毫无意义,就算作为“低水平”也不行,起码理论方法得用DFT

guess=read geom=check opt=(readfc,ts,noeigentest)  此时不要输入文件里写坐标
北京科音自然科学研究中心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

188

帖子

0

威望

922

eV
积分
1110

Level 4 (黑子)

8459#
发表于 Post on 2019-2-28 21:07:26 | 只看该作者 Only view this author
sobereva 发表于 2019-2-28 19:34
用hf/3-21g计算根本毫无意义,就算作为“低水平”也不行,起码理论方法得用DFT

guess=read geom=check ...

好的,谢谢

50

帖子

0

威望

472

eV
积分
522

Level 4 (黑子)

8460#
发表于 Post on 2019-2-28 22:52:26 | 只看该作者 Only view this author
sob老师您好,在键解离能的计算中用的原理及公式如图所示。
请问公式是在0k下推导的,然而gaussian中输出的焓却是298.15k下的,那岂不是就不能带入此公式了吗?

QQ截图20190228224414.png (53.86 KB, 下载次数 Times of downloads: 18)

QQ截图20190228224414.png

本版积分规则 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.286664 second(s), 22 queries , Gzip On.

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