计算化学公社

 找回密码 Forget password
 注册 Register

用NWChem做SODFT在DFT计算中考虑旋轨耦合效应

查看数: 19136 | 评论数: 15 | 收藏 Add to favorites 13
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2017-4-12 23:23

正文摘要:

用NWChem做SODFT在DFT计算中考虑旋轨耦合效应Use NWChem to do SODFT and consider spin-orbit coupling effect in DFT calculations 文/Sobereva @北京科音First release: 2017-Apr-12   Last update: ...

回复 Reply

beefly 发表于 Post on 2024-8-18 12:01:29
本帖最后由 beefly 于 2024-8-18 12:04 编辑
冰释之川 发表于 2017-4-18 13:01
问题解决了,是Linear dependencies的锅:

Unfortunately, the spin-orbit code does not correctly h ...

nwchem判断基组线性依赖的默认阈值是1.0e-5。由于sodft不支持消除冗余积函数的算法,这么高的阈值显然不适用。可以用手册里没有讲的一个命令降低阈值。根据出错信息里的“Smallest S eigenvalue”,假如此值在1.0e-8以上,可以用

set lindep:tol 1.0e-8

但是如果出现更低的特征值,降低阈值的做法就不保险了,可能会带来其他问题,比如scf不收敛。
sobereva 发表于 Post on 2023-12-3 00:53:24
funyboy 发表于 2023-12-2 21:11
老师,您好,第3步高精度的考虑标量相对论的单点,还是在gaussian算的吧?是不是也可以在NWchem计算?但是 ...

可以
看手册就完了
funyboy 发表于 Post on 2023-12-2 21:11:52
老师,您好,第3步高精度的考虑标量相对论的单点,还是在gaussian算的吧?是不是也可以在NWchem计算?但是不知道NWchem中CCSD(T)单点计算的输入文件怎么写?
funyboy 发表于 Post on 2023-9-22 10:51:56
sobereva 发表于 2023-9-21 01:44
这是我当时的输入输出文件,对照着看吧。也没准是版本问题

好的,谢谢sob老师
sobereva 发表于 Post on 2023-9-21 01:44:02
funyboy 发表于 2023-9-20 16:29
老师,您好,早上在群里咨询您的我单独算Tl的时候每次结果不一样,

这是我当时的输入输出文件,对照着看吧。也没准是版本问题

Tl.nw (2.89 KB, 下载次数 Times of downloads: 5)
Tl.out (12.9 KB, 下载次数 Times of downloads: 1)

并行计算重现性问题见下文的讨论。重现性过差可能是程序本身问题,也可能是编译问题,可以咨询开发者
数值误差对计算化学结果重现性的影响
http://sobereva.com/88

funyboy 发表于 Post on 2023-9-20 16:29:32
老师,您好,早上在群里咨询您的我单独算Tl的时候每次结果不一样,


经您提醒关闭并行,能重复了,结果是:Total SO-DFT energy = -172.724904317135,与您算出的结果是Tl=-172.724034492809差别比较大。
1、我把输入文件和计算结果上传,请您帮我分析一下问题出在哪里,是不是我输入文件编写有问题?
2、请问如何解决多核并行结果不能严格重现的问题?


谢谢


NWChem_SODFT.rar (9.73 KB, 下载次数 Times of downloads: 0)
冰释之川 发表于 Post on 2017-4-18 13:01:59
本帖最后由 冰释之川 于 2017-4-18 13:03 编辑
sobereva 发表于 2017-4-15 19:20
核数设少点,内存给大点再试

问题解决了,是Linear dependencies的锅:

Unfortunately, the spin-orbit code does not correctly handle the linear dependencies associated with your case.

The work-around is to add the line
  1. set lindep:n_dep 0
复制代码
just before the task line, i.e.

  1. set lindep:n_dep 0
  2. task sodft
复制代码

看来NWChem 做SODFT没法处理基函数的线性相关问题,会报错。
sobereva 发表于 Post on 2017-4-15 19:20:24
冰释之川 发表于 2017-4-15 14:37
我这边测试了 把其它元素改成def2-SVP,能够正常迭代。奇怪的是就是设置成def2-TZVP,两分钟后就歇菜了。 ...

核数设少点,内存给大点再试
冰释之川 发表于 Post on 2017-4-15 14:37:29
sobereva 发表于 2017-4-15 13:40
由于def2-TZVP耗时太长,我把其它元素改成def2-SVP,把print low去掉,已经能正常开始迭代了
另外建议 ...

我这边测试了 把其它元素改成def2-SVP,能够正常迭代。奇怪的是就是设置成def2-TZVP,两分钟后就歇菜了。我memory total 1000 mb 开的是4个逻辑核。
sobereva 发表于 Post on 2017-4-15 13:40:28
冰释之川 发表于 2017-4-15 12:21
现在测试下来,就该输入文件没法跑SODFT(如果跑DFT能跑),感觉应该是程序问题?不知道Sob老师能协助测 ...


由于def2-TZVP耗时太长,我把其它元素改成def2-SVP,把print low去掉,已经能正常开始迭代了
另外建议用memory关键词把内存设大一些,默认的很小。
冰释之川 发表于 Post on 2017-4-15 12:21:37
sobereva 发表于 2017-4-14 20:49
如果输入文件只提交一个能跑,同时提交两个没必然性导致出错,不过物理核心只有四个,提交两个4核并行的 ...

现在测试下来,就该输入文件没法跑SODFT(如果跑DFT能跑),感觉应该是程序问题?不知道Sob老师能协助测试一下,大概2分钟内报错
sobereva 发表于 Post on 2017-4-14 20:49:53
冰释之川 发表于 2017-4-14 13:12
ga_orthog: hard zero                1093
--------------------------------------------------------- ...

如果输入文件只提交一个能跑,同时提交两个没必然性导致出错,不过物理核心只有四个,提交两个4核并行的任务效率肯定很烂
冰释之川 发表于 Post on 2017-4-14 13:12:46
本帖最后由 冰释之川 于 2017-4-14 13:25 编辑

ga_orthog: hard zero                1093
------------------------------------------------------------------------
------------------------------------------------------------------------
  current input line :
ga_orthog: hard zero                1093
------------------------------------------------------------------------
------------------------------------------------------------------------
  current input line :
ga_orthog: hard zero                1093
------------------------------------------------------------------------
------------------------------------------------------------------------
  current input line :
------------------------------------------------------------------------
ga_orthog: hard zero                1093
------------------------------------------------------------------------
------------------------------------------------------------------------
current input line :
   129: task sodft
------------------------------------------------------------------------
------------------------------------------------------------------------
This error has not yet been assigned to a category
------------------------------------------------------------------------

这个报错是指MPI并行问题吗?因为我们服务器每个机子是四核八线程,难道同时提交两个任务,每个任务个逻辑4核会导致任务break down吗?
Ben_2Te_SODFT.rar (7.55 KB, 下载次数 Times of downloads: 15)
chemhou 发表于 Post on 2017-4-13 12:06:20
赞!

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

GMT+8, 2024-11-27 21:19 , Processed in 0.325719 second(s), 27 queries , Gzip On.

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