计算化学公社

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

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

   关闭 [复制链接 Copy URL]

6万

帖子

99

威望

6万

eV
积分
125208

管理员

公社社长

3766#
 楼主 Author| 发表于 Post on 2017-6-9 21:29:52 | 只看该作者 Only view this author
冰释之川 发表于 2017-6-9 14:42
关于高斯超极化率β输出文件的解读。我将高斯输出文件放入Multiwfn进行parse之后,发现β各个分量关系如下 ...


由于超极化率有张量对称性,所以程序只记录对称唯一的部分。没必要太关注这个
北京科音自然科学研究中心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

581

帖子

0

威望

1700

eV
积分
2281

Level 5 (御坂)

3767#
发表于 Post on 2017-6-10 08:35:34 | 只看该作者 Only view this author
sobereva 发表于 2017-6-9 21:25
不要引这个。非要引IRC可以引Fukui的IRC原文

谢谢sob老师

113

帖子

0

威望

271

eV
积分
384

Level 3 能力者

3768#
发表于 Post on 2017-6-10 12:20:09 | 只看该作者 Only view this author
sob老师,请问gaussview有没有设置快捷键的功能?

6万

帖子

99

威望

6万

eV
积分
125208

管理员

公社社长

3769#
 楼主 Author| 发表于 Post on 2017-6-10 21:48:06 | 只看该作者 Only view this author
erya123 发表于 2017-6-10 12:20
sob老师,请问gaussview有没有设置快捷键的功能?

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

31

帖子

0

威望

134

eV
积分
165

Level 3 能力者

3770#
发表于 Post on 2017-6-11 11:42:05 | 只看该作者 Only view this author
sob老师您好,我想请教几个Gaussian的问题:
1、gaussian怎么计算体系(eg:金属团簇)在变化的电场下的电荷密度分布,http://sobereva.com/113中   4.1   施加电场后COCl2电子密度的变化的等值面图(在外场作用下,电子密度分布会被极化。这里考察一下COCl2分子在C-O轴线方向上加上0.03a.u.均匀电场...)——这里加的是均匀场。
2、按理解来说,在加外场的情况下,如果外场加的很大则难以计算,是否能用gaussian中线性响应理论加小的电场来计算体系的电荷密度分布?或者说加小的单频外场。
感谢sob老师抽空答疑
你生而有翼,为何竟愿一生匍匐前进,形如虫蚁?-——贾拉尔•阿德丁•鲁米

6万

帖子

99

威望

6万

eV
积分
125208

管理员

公社社长

3771#
 楼主 Author| 发表于 Post on 2017-6-11 12:16:15 | 只看该作者 Only view this author
785524519 发表于 2017-6-11 11:42
sob老师您好,我想请教几个Gaussian的问题:
1、gaussian怎么计算体系(eg:金属团簇)在变化的电场下的电 ...


1 直接用关键词没法加
如果你想考察变化的外场各个时刻体系密度的变化,可以自写个脚本实现,每一步用不同的外场强度/方向并作图,最后合并成比如gif动画

2 Gaussian的field关键词加的外场是不含频的,加的是静态的外场。
可以加大的外场,只不过容易SCF不收敛而已,此时优化结构的话,也够也容易变形或者崩溃(这是实际情况,不是bug或数值不稳定)
北京科音自然科学研究中心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

31

帖子

0

威望

134

eV
积分
165

Level 3 能力者

3772#
发表于 Post on 2017-6-11 13:57:05 | 只看该作者 Only view this author
sobereva 发表于 2017-6-11 12:16
1 直接用关键词没法加
如果你想考察变化的外场各个时刻体系密度的变化,可以自写个脚本实现,每一步用 ...

谢谢sob老师这么及时的回答。但可能是我没有说清楚,实际上主要想问考虑与时间有关的电场不知道Gaussian能不能算电子密度的变化。因为这个涉及到含时薛定谔方程,应该不能按照算静电场下的电子密度变化的办法用每一步不同的电场来算。用一般的大小的电场来算,计算随时间变化的电子密度,因涉及到求解含时薛定谔方程,应该很难计算。如果用线性响应理论,相当于含时薛定谔方程中把外场看成微扰并只取一阶近似,就会容易的多。我看您的帖子中,提到过溶剂的极化可以用线性响应理论计算,不知能不能用到这里。
你生而有翼,为何竟愿一生匍匐前进,形如虫蚁?-——贾拉尔•阿德丁•鲁米

78

帖子

0

威望

814

eV
积分
892

Level 4 (黑子)

3773#
发表于 Post on 2017-6-11 17:43:54 | 只看该作者 Only view this author
请问老师, 用opt ModRedundant 关键词扫描二面角,这些输入行是什么意思啊,谢谢老师了。
Lines in a ModRedundant input section use the following syntax:

[Type] N1 [N2 [N3 [N4]]] [A | F] [[min] max]]
[Type] N1 [N2 [N3 [N4]]] B [[min] max]]
[Type] N1 [N2 [N3 [N4]]] K|R [[min] max]]
[Type] N1 [N2 [N3 [N4]]] D [[min] max]]
[Type] N1 [N2 [N3 [N4]]] H diag-elem [[min] max]]
[Type] N1 [N2 [N3 [N4]]] S nsteps stepsize [[min] max]]
N1, N2, N3 and N4 are atom numbers or wildcards (discussed below). Atom numbering begins at 1, and any dummy atoms are not counted.
自己决定自己的DNA

6万

帖子

99

威望

6万

eV
积分
125208

管理员

公社社长

3774#
 楼主 Author| 发表于 Post on 2017-6-11 20:10:54 | 只看该作者 Only view this author
QuantumHan 发表于 2017-6-11 17:43
请问老师, 用opt ModRedundant 关键词扫描二面角,这些输入行是什么意思啊,谢谢老师了。
Lines in a Mod ...

就是基本格式
[]里的是可选的
简单来说,opt=modredundant,然后例如末尾空一行写17 15 1 2 S 18 10.0,就表示扫描17-15-1-2二面角,18步(共19个点),步长10度。
北京科音自然科学研究中心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
积分
125208

管理员

公社社长

3775#
 楼主 Author| 发表于 Post on 2017-6-11 20:13:47 | 只看该作者 Only view this author
785524519 发表于 2017-6-11 13:57
谢谢sob老师这么及时的回答。但可能是我没有说清楚,实际上主要想问考虑与时间有关的电场不知道Gaussian ...


含时外场情况计算一般用octopus做RT-TDDFT计算,Gaussian派不上用场
北京科音自然科学研究中心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

78

帖子

0

威望

814

eV
积分
892

Level 4 (黑子)

3776#
发表于 Post on 2017-6-11 22:36:10 | 只看该作者 Only view this author
sobereva 发表于 2017-6-11 20:10
就是基本格式
[]里的是可选的
简单来说,opt=modredundant,然后例如末尾空一行写17 15 1 2 S 18 10.0 ...

谢谢Sob老师
自己决定自己的DNA

73

帖子

0

威望

144

eV
积分
217

Level 3 能力者

3777#
发表于 Post on 2017-6-11 23:25:45 | 只看该作者 Only view this author
请问sob老师,计算ecd光谱时,增加激发态数,得到的理论结果只能覆盖到光谱图越低的波长吗?如何也覆盖到光谱图越高的波长?

6万

帖子

99

威望

6万

eV
积分
125208

管理员

公社社长

3778#
 楼主 Author| 发表于 Post on 2017-6-12 00:53:26 | 只看该作者 Only view this author
qcn 发表于 2017-6-11 23:25
请问sob老师,计算ecd光谱时,增加激发态数,得到的理论结果只能覆盖到光谱图越低的波长吗?如何也覆盖到光 ...


高波长对应低激发能,程序正常情况是按照能量由低到高计算的,自然而然就把能算出来的高波长的所有态都算出来了
北京科音自然科学研究中心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

31

帖子

0

威望

134

eV
积分
165

Level 3 能力者

3779#
发表于 Post on 2017-6-12 09:31:46 | 只看该作者 Only view this author
sobereva 发表于 2017-6-11 20:13
含时外场情况计算一般用octopus做RT-TDDFT计算,Gaussian派不上用场

感谢sob老师的耐心指点! 那我再去研究下
你生而有翼,为何竟愿一生匍匐前进,形如虫蚁?-——贾拉尔•阿德丁•鲁米

27

帖子

0

威望

399

eV
积分
426

Level 3 能力者

3780#
发表于 Post on 2017-6-12 10:32:41 | 只看该作者 Only view this author
各位好,请问我在优化一个复合物时,其中一个片段的氢总是转移到另一个片段上去,而导致不收敛,这个问题如何解决?谢谢

本版积分规则 Credits rule

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

GMT+8, 2026-2-27 19:54 , Processed in 0.260674 second(s), 20 queries , Gzip On.

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