计算化学公社

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

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

  [复制链接 Copy URL]

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

676#
发表于 Post on 2023-2-26 18:28:16 | 只看该作者 Only view this author
lurensan 发表于 2023-2-26 06:41
各位老师好,我使用B3LYP-D3-def2svp优化包含Fe和Ir的107个原子的结构,优化两天了还没有收敛,只有两个Y ...

优化刚刚进入平台期,还不算是出现长时间没有进展的情形。个人觉得可以再优化个20~30步,如果还是不收敛就杀掉任务,从最后一帧开始用较小步长继续优化
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

677#
 楼主 Author| 发表于 Post on 2023-2-27 04:05:12 | 只看该作者 Only view this author
lurensan 发表于 2023-2-26 13:41
各位老师好,我使用B3LYP-D3-def2svp优化包含Fe和Ir的107个原子的结构,优化两天了还没有收敛,只有两个Y ...

放大末期部分仔细看
如果有收敛趋势,再跑个几十步
如果震荡明显了,按下文解决
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

678#
 楼主 Author| 发表于 Post on 2023-2-27 04:08:30 | 只看该作者 Only view this author
糖糖糖豆9988 发表于 2023-2-26 09:16
%chk=eTS_C2-H7-CH3.chk
%nprocs=8
#p UCCSD(T,T1diag) cc-pVTZ

就8个核算什么CCSD(T)
48个核都给CCSD(T)用
并且不要同时运行多个CCSD(T)任务免得硬盘I/O争抢导致哪个都算不痛快
内存怎么分配的要明确说清楚,CCSD(T)耗内存很大

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

477

帖子

0

威望

887

eV
积分
1364

Level 4 (黑子)

679#
发表于 Post on 2023-2-27 09:41:18 | 只看该作者 Only view this author
请教各位老师,输入文件如图所示,但是out文件显示没有计算。扫描的时候我都是只换结构,关键词套用原来的,但是这个计算不下去。请问这个是高斯的问题吗?谢谢

202302270939067647..png (25.4 KB, 下载次数 Times of downloads: 23)

out文件

out文件

202302270939416528..png (47.22 KB, 下载次数 Times of downloads: 19)

输入文件

输入文件

40

帖子

0

威望

565

eV
积分
605

Level 4 (黑子)

680#
发表于 Post on 2023-2-27 10:12:20 | 只看该作者 Only view this author
本帖最后由 Lyy1 于 2023-2-27 17:25 编辑

请问各位老师用高斯算自旋耦合的时候,无报错断了之后,如何续算呢?

202302271725064904..png (97.1 KB, 下载次数 Times of downloads: 23)

202302271725064904..png

477

帖子

0

威望

887

eV
积分
1364

Level 4 (黑子)

681#
发表于 Post on 2023-2-27 11:06:49 | 只看该作者 Only view this author
本帖最后由 糖糖糖豆9988 于 2023-2-27 11:07 编辑
sobereva 发表于 2023-2-27 04:08
就8个核算什么CCSD(T)
48个核都给CCSD(T)用
并且不要同时运行多个CCSD(T)任务免得硬盘I/O争抢导致哪个 ...

社长,我现在申请到一个64核的服务器,管理员说物理内存是128G,我用free查看了一下,结果如图。①那么我就用这个服务器,一次计算一个CCSD(T)计算,输入文件写明%nprocs=64,%mem=60GB可以吗?②如果这样设置了,对于我想算的联环丁烷结构氢提取反应(25个原子左右),您觉得一次计算任务大概耗时多久呢?谢谢社长

202302271106158296..png (62.55 KB, 下载次数 Times of downloads: 19)

202302271106158296..png

165

帖子

0

威望

836

eV
积分
1001

Level 4 (黑子)

682#
发表于 Post on 2023-2-27 13:04:37 | 只看该作者 Only view this author
sobereva 发表于 2023-2-26 03:30
你的模型根本没法用
用CP2K算这个比用Gaussian强百倍
一般都需要4层表面原子,边界还得加氢,得考虑表 ...

sob老师好,我将金刚石换成了过渡金属体系,请这样构建是否合理呢?链接:https://pan.baidu.com/s/1i_eD9e-5-DYpDw1BYHfQ5w?pwd=0000
提取码:0000
社会浮躁

573

帖子

0

威望

1652

eV
积分
2225

Level 5 (御坂)

683#
发表于 Post on 2023-2-27 15:36:21 | 只看该作者 Only view this author
各位老师好,哪里可以下载windows版的32位gview呀?

5

帖子

0

威望

153

eV
积分
158

Level 3 能力者

684#
发表于 Post on 2023-2-27 15:46:23 | 只看该作者 Only view this author
老师你好,请问计算这个分子的福井函数时,由于分子本身是中性,双重态,那么N+1和N-1的电荷和重态分别是1 3和-1 1吗?

202302271546164695..png (51.74 KB, 下载次数 Times of downloads: 19)

202302271546164695..png

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

685#
发表于 Post on 2023-2-27 18:26:37 | 只看该作者 Only view this author
1000ninus7 发表于 2023-2-27 08:46
老师你好,请问计算这个分子的福井函数时,由于分子本身是中性,双重态,那么N+1和N-1的电荷和重态分别是1  ...

N指的是电子数,所以N+1的电荷是-1,N-1的电荷是+1。
这个体系的阳离子和阴离子的自旋多重度(注意没有“重态”这个说法,“重态”两字不能单独成词,必须叫自旋多重度)应该都是1,画一下共振式就知道了
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

165

帖子

0

威望

836

eV
积分
1001

Level 4 (黑子)

686#
发表于 Post on 2023-2-27 18:57:37 | 只看该作者 Only view this author
sobereva 发表于 2023-2-26 03:30
你的模型根本没法用
用CP2K算这个比用Gaussian强百倍
一般都需要4层表面原子,边界还得加氢,得考虑表 ...

另外还想请教一下sob老师,请问博文540中的体系为何银是3层呢?边界需要加氢饱和的情况有哪些呢?
社会浮躁

269

帖子

0

威望

2060

eV
积分
2329

Level 5 (御坂)

687#
发表于 Post on 2023-2-27 20:28:31 | 只看该作者 Only view this author
想跟大家请教一下,想看一下Pt原子负载在一个金属氧化物载体上(其中Pt与金属氧化物载体的4个O成键),想确定一下这个结构中Pt-O的成键轨道分别是那些,有什么方法么(感觉从轨道电子云的分布不是太确定)。
努力做计算的小菜鸡

5

帖子

0

威望

153

eV
积分
158

Level 3 能力者

688#
发表于 Post on 2023-2-27 21:06:49 | 只看该作者 Only view this author
wzkchem5 发表于 2023-2-27 18:26
N指的是电子数,所以N+1的电荷是-1,N-1的电荷是+1。
这个体系的阳离子和阴离子的自旋多重度(注意没有 ...

好的,谢谢老师

269

帖子

0

威望

2060

eV
积分
2329

Level 5 (御坂)

689#
发表于 Post on 2023-2-27 21:35:49 | 只看该作者 Only view this author
各位老师好,想跟大家请教一下在使用Multiwfn+VMD快速绘制高质量分子轨道等值面图(http://sobereva.com/447),可以指定只显示某些原子之间的轨道图么?
努力做计算的小菜鸡

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

690#
发表于 Post on 2023-2-27 21:53:41 | 只看该作者 Only view this author
吴双 发表于 2023-2-27 14:35
各位老师好,想跟大家请教一下在使用Multiwfn+VMD快速绘制高质量分子轨道等值面图(http://sobereva.com/44 ...

请先定义什么叫“某些原子之间的轨道图”
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

本版积分规则 Credits rule

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

GMT+8, 2024-11-24 09:48 , Processed in 0.190280 second(s), 24 queries , Gzip On.

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