计算化学公社

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

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

   关闭 [复制链接 Copy URL]

259

帖子

0

威望

851

eV
积分
1110

Level 4 (黑子)

6856#
发表于 Post on 2018-6-3 00:28:58 | 只看该作者 Only view this author
本帖最后由 lufenghxu 于 2018-6-3 00:35 编辑
sobereva 发表于 2018-6-2 20:43
跑什么曲线?
如果是弛豫过程的能量极小路径,高斯没法跑

老师,找到的交叉点是不是就是两个不同态反应路径的交叉点?找到后再怎么办(分析)?

6万

帖子

99

威望

6万

eV
积分
125190

管理员

公社社长

6857#
 楼主 Author| 发表于 Post on 2018-6-3 01:17:45 | 只看该作者 Only view this author
lufenghxu 发表于 2018-6-3 00:28
老师,找到的交叉点是不是就是两个不同态反应路径的交叉点?找到后再怎么办(分析)?


不是
多看相关研究文献的讨论
北京科音自然科学研究中心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

198

帖子

0

威望

450

eV
积分
648

Level 4 (黑子)

6858#
发表于 Post on 2018-6-3 23:04:40 | 只看该作者 Only view this author
我用TD-DFT算gradient,SCF done 之后 density matrix 一直不收敛,试了好多种方法都不行,然后我把收敛标准降低(scf=conver=7),就收敛了,可是这样算出的gradient是不是就不准了?目标就是算原子所受的力。

95

帖子

0

威望

237

eV
积分
332

Level 3 能力者

6859#
发表于 Post on 2018-6-3 23:26:58 | 只看该作者 Only view this author
老师,下图描述的是不同重态的交叉点。下图中的纵坐标(能量)可以从log文件中提取,但是不同重态的反应路径的横坐标是如何提取,如何统一的?

11111.jpg (39.82 KB, 下载次数 Times of downloads: 49)

11111.jpg

6万

帖子

99

威望

6万

eV
积分
125190

管理员

公社社长

6860#
 楼主 Author| 发表于 Post on 2018-6-4 04:46:32 | 只看该作者 Only view this author
算一算 发表于 2018-6-3 23:04
我用TD-DFT算gradient,SCF done 之后 density matrix 一直不收敛,试了好多种方法都不行,然后我把收敛标 ...


降到7还是可以用的。Gaussian默认的收敛限很严,对密度矩阵变化要求稍微有点过头,同样叫做“tightSCF”,而Gaussian的就比ORCA的严很多。
北京科音自然科学研究中心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
积分
125190

管理员

公社社长

6861#
 楼主 Author| 发表于 Post on 2018-6-4 05:02:01 | 只看该作者 Only view this author
lily19830130 发表于 2018-6-3 23:26
老师,下图描述的是不同重态的交叉点。下图中的纵坐标(能量)可以从log文件中提取,但是不同重态的反应路 ...


图中体系很简单,直接对R(N-O)进行扫描,就能得到两个态的势能曲线
而对于复杂情况,需要对优化/搜索MECP过程的各个点依次计算其某个几何变量的数值,然后放到图上
北京科音自然科学研究中心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

198

帖子

0

威望

450

eV
积分
648

Level 4 (黑子)

6862#
发表于 Post on 2018-6-4 08:08:37 | 只看该作者 Only view this author
sobereva 发表于 2018-6-4 04:46
降到7还是可以用的。Gaussian默认的收敛限很严,对密度矩阵变化要求稍微有点过头,同样叫做“tightSCF ...

谢谢回答。我试了几个相似的体系,用8和7两种收敛限分别算,算出来的gradient确实差别不大,但是别的东西有时差别很大,比如分子轨道

11

帖子

0

威望

336

eV
积分
347

Level 3 能力者

6863#
发表于 Post on 2018-6-4 10:22:41 | 只看该作者 Only view this author
老师,我现在算的体系有160个原子比较大,先用6-31G*对CHONP,优化了结构,由于体系中有磷氧键,考虑增大ONP用的基组做优化,是否可以只对ONP使用6-311G*基组或者6-31+G*?体系中还含有Zn,试了非金属原子都用6-311G*,基函数有2000多,实在算的慢。

44

帖子

0

威望

193

eV
积分
237

Level 3 能力者

6864#
发表于 Post on 2018-6-4 10:30:29 | 只看该作者 Only view this author
老师,在Multiwfn中绘制填色图时,怎么将坐标轴上字体变大呢?

6万

帖子

99

威望

6万

eV
积分
125190

管理员

公社社长

6865#
 楼主 Author| 发表于 Post on 2018-6-4 10:35:47 | 只看该作者 Only view this author
新手囧囧 发表于 2018-6-4 10:30
老师,在Multiwfn中绘制填色图时,怎么将坐标轴上字体变大呢?


没法变
非要更大的,就导出数据到sigmaplot里重新绘图,届时随意调整。或者自行修改源代码。或者自行ps
Multiwfn的问题请到Multiwfn板块咨询,本帖只用于回答纯量化问题。
北京科音自然科学研究中心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
积分
125190

管理员

公社社长

6866#
 楼主 Author| 发表于 Post on 2018-6-4 10:36:35 | 只看该作者 Only view this author
算一算 发表于 2018-6-4 08:08
谢谢回答。我试了几个相似的体系,用8和7两种收敛限分别算,算出来的gradient确实差别不大,但是别的东西 ...

至少占据的分子轨道的系数应该相差很小
北京科音自然科学研究中心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
积分
125190

管理员

公社社长

6867#
 楼主 Author| 发表于 Post on 2018-6-4 10:41:59 | 只看该作者 Only view this author
yu947937636 发表于 2018-6-4 10:22
老师,我现在算的体系有160个原子比较大,先用6-31G*对CHONP,优化了结构,由于体系中有磷氧键,考虑增大ON ...

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

39

帖子

0

威望

235

eV
积分
274

Level 3 能力者

6868#
发表于 Post on 2018-6-4 10:56:24 | 只看该作者 Only view this author
sobereva 发表于 2018-6-2 21:13
用论坛首页google框搜索“量子产率”

谢谢老师,我之前搜索以后好像是说只靠Gaussian是没法模拟的
鸟语虫声,总是传心之诀;花英草色,无非见道之文。学者要天机清澈,胸次玲珑,触物皆有会心处。

11

帖子

0

威望

336

eV
积分
347

Level 3 能力者

6869#
发表于 Post on 2018-6-4 11:11:24 | 只看该作者 Only view this author

sob老师,看了您的博文写到6-311G*基组有时会比6-31+G*得到的结果更好,在网上也看到一篇博文计算的磷氧键用三重分裂基组和多极化基组算的键长更接近实验值,又查了相关体系的文献对ONP用的是加弥散的6-31+G*基组。比较纠结该用哪个?

20

帖子

0

威望

343

eV
积分
363

Level 3 能力者

6870#
发表于 Post on 2018-6-4 11:17:56 | 只看该作者 Only view this author
sob老师,我的分子是带有ONO基团的烷氧自由基,根据tddft计算激发态得到第四态符合实验波长,对应C-O成键轨道跃迁到O原子的非键轨道,在用casscf优化基态可以收敛,但激发态总是不能收敛,想麻烦您帮我看看结构

m_trans_aO_td.fchk

2.78 MB, 下载次数 Times of downloads: 0

m_trans_aO_td.log

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

s0.fchk

1.1 MB, 下载次数 Times of downloads: 0

s0.log

1.87 MB, 下载次数 Times of downloads: 0

s4.fchk

1.06 MB, 下载次数 Times of downloads: 0

s4.log

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

本版积分规则 Credits rule

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

GMT+8, 2026-2-26 15:48 , Processed in 0.664195 second(s), 22 queries , Gzip On.

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