计算化学公社

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

[综合交流] 简单量子化学问题答疑专帖

  [复制链接 Copy URL]

2

帖子

0

威望

73

eV
积分
75

Level 2 能力者

4996#
发表于 Post on 2026-2-4 10:59:17 | 只看该作者 Only view this author

各位老师好,如果在油品里面加入一些添加剂,和原来的油品比较的话,存在电子密度差的概念吗?我想通过模拟研究油品的绝缘性,想从偶极距和电子密度差两个方面着手,但是不是很确定电子密度差能不能这样做,我目前是是根据http://bbs.keinsci.com/thread-46380-1-1.html和《使用Multiwfn作电子密度差图》(http://sobereva.com/113)学习。还有一个问题是,如果我的这个思路是可行的,我应该研究的是油品和加了添加剂的油品的电子密度差,还是研究油品和加了添加剂的油品和添加剂三者的电子密度差呢?如果这个思路不可行的话,请问各位老师觉得有什么办法可以通过模拟研究油品的绝缘性或者介电损耗吗?

1238

帖子

6

威望

2533

eV
积分
3891

Level 5 (御坂)

傻傻的木瓜

4997#
发表于 Post on 2026-2-4 11:36:40 | 只看该作者 Only view this author
本帖最后由 Uus/pMeC6H4-/キ 于 2026-2-13 20:03 编辑
AprilYao 发表于 2026-2-4 10:59
各位老师好,如果在油品里面加入一些添加剂,和原来的油品比较的话,存在电子密度差的概念吗?我想通过模 ...

如果根本目的是研究“油品的绝缘性或者介电损耗”这么大的话题,请去“化学”版块单独发帖提问,而且不要期待量子化学计算的方法就是答案(参见X-Y problem的概念);有具体的既往文献研究了密度差的话可以举出来源链接来让别人看看计算是否合理,而如果只是没有根据的个人猜想那就不要提出来了(因为很容易带偏话题)。

确定此类宏观凝聚相动态的化工问题不是按连续介质建模、用有限元分析之类的方法处理的么?对少数孤立分子体系做计算的结果,未必与凝聚相的绝缘性有因果性甚至相关性,试图设计分子描述符间接建立联系的话也必须当心适用范围。此外施加外电场,特别是接近击穿电压的高场强的条件,处理起来也不简单;有一些电化学反应的相关讨论在http://bbs.keinsci.com/thread-57304-1-1.html及其中链接可供参考。

编辑:另见http://bbs.keinsci.com/thread-43640-1-1.html
√546=23.36664289109

84

帖子

0

威望

303

eV
积分
387

Level 3 能力者

4998#
发表于 Post on 2026-2-4 12:49:13 | 只看该作者 Only view this author
本帖最后由 NHX 于 2026-2-4 13:01 编辑
sobereva 发表于 2026-2-4 05:45
语义不明,说清楚是自己算单点时候还是用sobMECP搜索MECP的时候
直接给完整gjf文件,并且直接说清楚文字 ...

抱歉sob老师,我没说清楚。是在搜索MECP的时候出现的报错。报错信息是No analytic first derivatives for this method.因为我是自己产生的chk文件,因此在./prepare.sh后,直接进行了./runMECP.sh。相关输入输出文件已上传,请老师指点

Input_Header_A

118 Bytes, 下载次数 Times of downloads: 0

Input_Header_B

118 Bytes, 下载次数 Times of downloads: 0

Input_Tail

3.76 KB, 下载次数 Times of downloads: 0

Job0_A.gjf

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

Job0_A.log

4.3 KB, 下载次数 Times of downloads: 0

Job0_B.gjf

3.88 KB, 下载次数 Times of downloads: 0

Job0_B.log

4.3 KB, 下载次数 Times of downloads: 0

6万

帖子

99

威望

6万

eV
积分
125115

管理员

公社社长

4999#
 楼主 Author| 发表于 Post on 2026-2-5 03:28:49 | 只看该作者 Only view this author
NHX 发表于 2026-2-4 12:49
抱歉sob老师,我没说清楚。是在搜索MECP的时候出现的报错。报错信息是No analytic first derivatives for ...

使用int=dkhso时没有解析梯度所致
尝试用force=enonly要求用有限差分方式算受力。但sobMECP没有测试过这种情况,如果实测不行的话,用赝势基组,或者修改sobMECP代码
北京科音自然科学研究中心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
积分
125115

管理员

公社社长

5000#
 楼主 Author| 发表于 Post on 2026-2-5 03:30:56 | 只看该作者 Only view this author
AprilYao 发表于 2026-2-4 10:59
各位老师好,如果在油品里面加入一些添加剂,和原来的油品比较的话,存在电子密度差的概念吗?我想通过模 ...

这个问题和电子密度差毫无关系
北京科音自然科学研究中心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

84

帖子

0

威望

303

eV
积分
387

Level 3 能力者

5001#
发表于 Post on 2026-2-5 08:37:25 | 只看该作者 Only view this author
本帖最后由 NHX 于 2026-2-5 10:32 编辑
sobereva 发表于 2026-2-5 03:28
使用int=dkhso时没有解析梯度所致
尝试用force=enonly要求用有限差分方式算受力。但sobMECP没有测试过这 ...

感谢sob老师回复,我已尝试,能够正常运行了,非常感谢sob老师!

108

帖子

0

威望

1170

eV
积分
1278

Level 4 (黑子)

5002#
发表于 Post on 2026-2-6 16:22:44 | 只看该作者 Only view this author
Hi everyone,
I want to calculate the energy plot as the phenol ring rotates around the axis perpendicular to the benzene ring. I create 2 dummy atoms to create the rotation axis, but the ModRedundant does not allow dummy atoms. Are there any solutions for this?
Thank you

202602061622268728..png (181.78 KB, 下载次数 Times of downloads: 0)

202602061622268728..png

6万

帖子

99

威望

6万

eV
积分
125115

管理员

公社社长

5003#
 楼主 Author| 发表于 Post on 2026-2-7 03:50:34 | 只看该作者 Only view this author
moiumuio 发表于 2026-2-6 16:22
Hi everyone,
I want to calculate the energy plot as the phenol ring rotates around the axis perpend ...

You have to write a script to generate a batch of coordinate corresponding to every point of the (rigid) scan, and perform single point calculation for them, and finally extract energy from the output files.
北京科音自然科学研究中心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

108

帖子

0

威望

1170

eV
积分
1278

Level 4 (黑子)

5004#
发表于 Post on 2026-2-7 12:36:57 | 只看该作者 Only view this author
sobereva 发表于 2026-2-7 03:50
You have to write a script to generate a batch of coordinate corresponding to every point of the ( ...

Hello teacher,
In this case, if I do a rigid scan instead of a relaxed scan, would the information I have be usable since the first point in the optimized structure?
Thank you.

6万

帖子

99

威望

6万

eV
积分
125115

管理员

公社社长

5005#
 楼主 Author| 发表于 Post on 2026-2-7 16:55:11 | 只看该作者 Only view this author
moiumuio 发表于 2026-2-7 12:36
Hello teacher,
In this case, if I do a rigid scan instead of a relaxed scan, would the informatio ...

I cannot understand your meaning...
北京科音自然科学研究中心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

1

帖子

0

威望

19

eV
积分
20

Level 1 能力者

5006#
发表于 Post on 2026-2-7 17:53:37 | 只看该作者 Only view this author
本帖最后由 布兰柯 于 2026-2-7 22:49 编辑

请问有离子配位时,还能称为计算键解离能吗
键解离能指溶剂分子均裂生成自由基的反应中焓的变化,对于烷氧分子与离子配位时,比如研究[A+](EtOEt)中,亚甲基C-H键能,
(1)此时溶剂的键解离能 还存在吗?是不是不能用包含离子的自由基片段+氢自由基-配合物的焓的变化来代表键解离能了(此时总焓变中可能包含了离子与氧作用的焓变?)
(2)有离子配位时,是否就只能说研究的是断键反应的焓变了,谢谢老师
科音学员

1万

帖子

0

威望

9849

eV
积分
22083

Level 6 (一方通行)

5007#
发表于 Post on 2026-2-7 20:52:36 | 只看该作者 Only view this author
布兰柯 发表于 2026-2-7 17:53
请问有离子配位时,还能称为计算键解离能吗
键解离能指溶剂分子均裂生成自由基的反应中焓的变化,对于烷氧 ...

计算结果仍然可以叫键解离能,只不过是“和A配位的溶剂的键解离能”,不能不提A而直接称其为溶剂的键解离能

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
布兰柯 + 2 谢谢

查看全部评分 View all ratings

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
本团队长期招收研究生,有意者可私信联系

84

帖子

0

威望

303

eV
积分
387

Level 3 能力者

5008#
发表于 Post on 2026-2-7 23:17:07 | 只看该作者 Only view this author
本帖最后由 NHX 于 2026-2-7 23:20 编辑

各位老师好,学生目前在计算过渡金属阳离子与各种气体的反应路径,考虑到可能存在自旋翻转,所以计算势能面时不止计算了一个自旋多重度,最终在给出势能面上各结构的电子能量时想同时考虑SOC的影响,即得到考虑SOC后的最终电子能量。经过几天的学习我发现CASSCF能实现我想要的结果,但是里面的nroots参数目前还不知道具体如何设置,该参数的设置和活性空间有关系吗?我发现计算激发态时也涉及到这个参数,我不知道我这个计算目的涉不涉及激发态,好像仅仅涉及不同自旋多重度势能面的考虑,希望各位老师指点。

108

帖子

0

威望

1170

eV
积分
1278

Level 4 (黑子)

5009#
发表于 Post on 2026-2-8 00:52:47 | 只看该作者 Only view this author
sobereva 发表于 2026-2-7 16:55
I cannot understand your meaning...

Hi teacher,
I just want to ask, should I use rigid PES (only do SP) or relaxed PES (Opt and SP)?
Thank you for your answer.

6万

帖子

99

威望

6万

eV
积分
125115

管理员

公社社长

5010#
 楼主 Author| 发表于 Post on 2026-2-8 05:33:35 | 只看该作者 Only view this author
moiumuio 发表于 2026-2-8 00:52
Hi teacher,
I just want to ask, should I use rigid PES (only do SP) or relaxed PES (Opt and SP)?
...

Performing rigid scan is the easiest.
Relaxed scan may be feasible using the GIC feature of Gaussian. GIC allows users to perform a relaxed scan for a user-defined coordinate (including dihredral), which may involve center of ring(s). But the GIC code is not very robust and may fail.
北京科音自然科学研究中心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-16 15:05 , Processed in 0.181621 second(s), 24 queries , Gzip On.

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