计算化学公社

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

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

  [复制链接 Copy URL]

147

帖子

0

威望

523

eV
积分
670

Level 4 (黑子)

2011#
发表于 Post on 2023-7-28 16:09:36 | 只看该作者 Only view this author
想请问下各位老师,最近想做一个金属配合物的离子的casscf计算,最开始用dft验证了三重态是能量最低体系,请问这种情况下活性空间轨道要怎么选择呢?要把α和β轨道分开考虑吗?

1万

帖子

0

威望

9740

eV
积分
21936

Level 6 (一方通行)

2012#
发表于 Post on 2023-7-28 16:13:12 | 只看该作者 Only view this author
sxh3231796 发表于 2023-7-28 09:09
想请问下各位老师,最近想做一个金属配合物的离子的casscf计算,最开始用dft验证了三重态是能量最低体系, ...

CASSCF的活性轨道永远是restricted的,也就是alpha和beta轨道的数目、形状必须严格一致。
活性空间选择的方法和金属在p、d还是f区、金属轨道和周围配体轨道有没有强烈的相互作用等有关,和自旋多重度基本没关系
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
本团队长期招收研究生,有意者可私信联系

147

帖子

0

威望

523

eV
积分
670

Level 4 (黑子)

2013#
发表于 Post on 2023-7-28 16:47:09 | 只看该作者 Only view this author
wzkchem5 发表于 2023-7-28 16:13
CASSCF的活性轨道永远是restricted的,也就是alpha和beta轨道的数目、形状必须严格一致。
活性空间选择 ...

谢谢老师,还想再问一下那还能使用之前计算三重态的chk文件作为初猜吗?需要用guess=mix作为关键词来处理轨道吗?

1万

帖子

0

威望

9740

eV
积分
21936

Level 6 (一方通行)

2014#
发表于 Post on 2023-7-28 17:22:34 | 只看该作者 Only view this author
sxh3231796 发表于 2023-7-28 09:47
谢谢老师,还想再问一下那还能使用之前计算三重态的chk文件作为初猜吗?需要用guess=mix作为关键词来处理 ...

最好用ROHF轨道等作为初猜。如果用UHF轨道,不同程序的处理方式不一样,有的程序可能会报错,有的程序会用alpha轨道,等等,我不知道高斯是哪种情况。
没有任何理由做guess=mix,你做的又不是自旋破缺计算
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
本团队长期招收研究生,有意者可私信联系

147

帖子

0

威望

523

eV
积分
670

Level 4 (黑子)

2015#
发表于 Post on 2023-7-28 18:57:08 | 只看该作者 Only view this author
wzkchem5 发表于 2023-7-28 17:22
最好用ROHF轨道等作为初猜。如果用UHF轨道,不同程序的处理方式不一样,有的程序可能会报错,有的程序会 ...

好的,好的,多谢老师

6万

帖子

99

威望

5万

eV
积分
124694

管理员

公社社长

2016#
 楼主 Author| 发表于 Post on 2023-7-28 22:34:40 | 只看该作者 Only view this author
sxh3231796 发表于 2023-7-28 16:47
谢谢老师,还想再问一下那还能使用之前计算三重态的chk文件作为初猜吗?需要用guess=mix作为关键词来处理 ...

Gaussian会用alpha轨道

除了RO外,对于非限制性波函数,也可以考虑用Multiwfn产生双正交化轨道当做初猜(导出fch格式,然后unfchk转成chk以供guess=read)
用于非限制性开壳层波函数的双正交化方法的原理与应用
http://sobereva.com/448http://bbs.keinsci.com/thread-11377-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

286

帖子

0

威望

2687

eV
积分
2973

Level 5 (御坂)

2017#
发表于 Post on 2023-7-30 15:50:49 | 只看该作者 Only view this author
本帖最后由 吴双 于 2023-7-30 15:56 编辑

想跟各位老师请教一下,在用高斯做柔性扫描的时候,可以冻结两个原子,保证键长不变么。

6万

帖子

99

威望

5万

eV
积分
124694

管理员

公社社长

2018#
 楼主 Author| 发表于 Post on 2023-7-30 19:40:40 | 只看该作者 Only view this author
吴双 发表于 2023-7-30 15:50
想跟各位老师请教一下,在用高斯做柔性扫描的时候,可以冻结两个原子,保证键长不变么。

可以
按下文相同的方式添加冗余内坐标的冻结设置
在Gaussian中做限制性优化的方法
http://sobereva.com/404http://bbs.keinsci.com/thread-9022-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

286

帖子

0

威望

2687

eV
积分
2973

Level 5 (御坂)

2019#
发表于 Post on 2023-7-31 10:34:01 | 只看该作者 Only view this author
sobereva 发表于 2023-7-30 19:40
可以
按下文相同的方式添加冗余内坐标的冻结设置
在Gaussian中做限制性优化的方法

好的 谢谢老师

100

帖子

2

威望

2089

eV
积分
2229

Level 5 (御坂)

2020#
发表于 Post on 2023-7-31 11:54:22 | 只看该作者 Only view this author
老师,请问下图中的盒子中各个距离如何通过VMD来显示呢?

202307311154073062..png (525.45 KB, 下载次数 Times of downloads: 10)

202307311154073062..png

6万

帖子

99

威望

5万

eV
积分
124694

管理员

公社社长

2021#
 楼主 Author| 发表于 Post on 2023-7-31 21:30:05 | 只看该作者 Only view this author
Lym 发表于 2023-7-31 11:54
老师,请问下图中的盒子中各个距离如何通过VMD来显示呢?

此问题和量子化学无关,请发到分子模拟版块
北京科音自然科学研究中心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

3

帖子

0

威望

123

eV
积分
126

Level 2 能力者

2022#
发表于 Post on 2023-8-1 12:52:51 | 只看该作者 Only view this author

IRC中某个点SCF不收敛

在做IRC过程中出现L502错误,SCF不收敛,通过减少stepsize步长也没用,#p IRC=(calcfc,maxpoints=200,LQA) ub3lyp/genecp em=gd3bj guess=read。求老师帮忙看下怎么回事

202308011252021761..png (83.06 KB, 下载次数 Times of downloads: 10)

202308011252021761..png

6万

帖子

99

威望

5万

eV
积分
124694

管理员

公社社长

2023#
 楼主 Author| 发表于 Post on 2023-8-1 13:01:40 | 只看该作者 Only view this author
673010195 发表于 2023-8-1 12:52
在做IRC过程中出现L502错误,SCF不收敛,通过减少stepsize步长也没用,#p IRC=(calcfc,maxpoints=200,LQA)  ...

照常解决
解决SCF不收敛问题的方法
http://sobereva.com/61
减小步长不可能完全没用。至少发个stepsize=4或者更小的输入输出文件出来以证明
北京科音自然科学研究中心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

威望

230

eV
积分
261

Level 3 能力者

2024#
发表于 Post on 2023-8-1 17:20:39 | 只看该作者 Only view this author
想请问各位老师,如果要比较两个配体和Mn离子配位能力的大小可以通过什么描述符比较呢,或者这种配体交换反应能够计算过渡态吗?

6万

帖子

99

威望

5万

eV
积分
124694

管理员

公社社长

2025#
 楼主 Author| 发表于 Post on 2023-8-1 17:24:43 | 只看该作者 Only view this author
white-album 发表于 2023-8-1 17:20
想请问各位老师,如果要比较两个配体和Mn离子配位能力的大小可以通过什么描述符比较呢,或者这种配体交换反 ...

最严格的是配体解离能
也可以用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

本版积分规则 Credits rule

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

GMT+8, 2026-1-25 10:10 , Processed in 0.362223 second(s), 29 queries , Gzip On.

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