计算化学公社

 找回密码 Forget password
 注册 Register

Multiwfn现已可以调用cubegen使静电势分析耗时有飞跃式的下降!

查看数: 15086 | 评论数: 19 | 收藏 Add to favorites 6
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2018-8-18 00:41

正文摘要:

注1:2020-Jul-4及以后版本的Multiwfn的自身的静电势代码计算速度较之前版本有脱胎换骨般提升,见《Multiwfn的计算静电势的内部代码速度得到了极大的提升!》(http://bbs.keinsci.com/thread-18268-1-1.html)、《M ...

回复 Reply

sobereva 发表于 Post on 2020-6-27 01:12:37
winnerwill 发表于 2020-6-26 05:39
一直奇怪我怎么改GAUSS_PDEF,cubegen的调用核数始终与nthreads保持一致的。原来Multiwfn调用cubegen是直接 ...

这个对实际应用不会造成太大不便。
让用户去设GAUSS_PDEF的话,给初学者去说明怎么设环境变量很费劲,初学者还很可能忽略掉
winnerwill 发表于 Post on 2020-6-26 05:39:36
一直奇怪我怎么改GAUSS_PDEF,cubegen的调用核数始终与nthreads保持一致的。原来Multiwfn调用cubegen是直接用nthreads这一变量来控制的(ncubegenthreads=nthreads)。感觉是不是要由GAUSS_PDEF控制会好点,即(ncubegenthreads=GAUSS_PDEF)。这样,当cubegen需要减少并行核数来增加平均内存的时候,Multiwfn前面的其他步骤仍可以有更多的核数来运行。
wuzhiyi 发表于 Post on 2019-4-26 03:51:13
谢谢sob
sobereva 发表于 Post on 2019-4-26 02:42:59
wuzhiyi 发表于 2019-4-26 00:37
一开始以为是弥散函数的问题,现在发现是cp的问题,嫁了cp之后,Total SCF Density就没了,而cubegen依赖To ...

Multiwfn从fch里载入的是轨道信息,不直接使用其中的密度矩阵,所以fch里没有密度矩阵时也可以计算静电势

counterpoise任务产生的chk绝对不要用于波函数分析目的,否则结果可能没意义

CHELPG、RESP等都是基于实空间函数(静电势)计算的,完全不怕弥散函数,而且带了弥散函数结果会更好些

关于做分析能不能带弥散函数这点详见此文Q36
Multiwfn FAQ
http://sobereva.com/452
wuzhiyi 发表于 Post on 2019-4-26 00:37:47
一开始以为是弥散函数的问题,现在发现是cp的问题,嫁了cp之后,Total SCF Density就没了,而cubegen依赖Total SCF Density。
Multiwfn居然可以在没有Total SCF Density的时候计算静电势。。。
wuzhiyi 发表于 Post on 2019-4-25 23:13:11
本帖最后由 wuzhiyi 于 2019-4-25 23:44 编辑

想问一下我对一个东西算单点# M062X/gen EmpiricalDispersion=GD3 Counterpoise=3
然后想用Multiwfn借用cubegen来计算原子电荷
却出现错误
Running: /system/software/linux-x86_64/gaussian/g16/A.03__pgi_16.5__atlas/g16/cubegen    4 potential=SCF "HIP_E127_Y158H.fch" ESPresult.cub -5 h < cubegenpt.txt > nouseout
sh: line 1:  3382 Segmentation fault      (core dumped) /system/software/linux-x86_64/gaussian/g16/A.03__pgi_16.5__atlas/g16/cubegen 4 potential=SCF "HIP_E127_Y158H.fch" ESPresult.cub -5 h < cubegenpt.txt > nouseout
在nouseout文件里显示
Will use up to     4 threads.
Top CubeGen          Thu Apr 25 15:57:34 2019
Aft UFChkP           Thu Apr 25 15:58:09 2019      139.98 cpu,       34.99 elapsed
Could not find " ","Total SCF Density" on formatted checkpoint file.
Error termination via Lnk1e at Thu Apr 25 15:58:10 2019.
想问一下这如何解决? 谢谢
因为体系中有负电荷所以算结合能的时候加上了弥散,我知道如果算Mulliken的话,不能加弥散,但不太清楚算MK,CHELPG或者RESP的时候,加弥散有没有坏处,因为单点已经算好了,所以想能不能循环利用一下chk文件。
sobereva 发表于 Post on 2019-4-18 17:17:56
15whyhq 发表于 2019-4-18 16:52
cubegenpath路径是服务器装g16的路径吗?

不是路径,而是要指向g16那个可执行文件名
15whyhq 发表于 Post on 2019-4-18 16:52:52
sobereva 发表于 2018-9-14 05:20
没版本的事,Multiwfn调用不了cubegen纯粹是cubegenpath路径没写对而已

cubegenpath路径是服务器装g16的路径吗?
sobereva 发表于 Post on 2018-9-14 05:20:03
tianpingpei 发表于 2018-9-13 14:05
sob老师,我下载最新的Multiwfn,计算高斯得到fchk文件,显示the cubegen can not be found,我的高斯是09w ...

没版本的事,Multiwfn调用不了cubegen纯粹是cubegenpath路径没写对而已
sobereva 发表于 Post on 2018-9-14 05:17:19
tianpingpei 发表于 2018-9-13 11:46
sob老师,请问一下,高斯处理chk文件的时候,必须等高斯计算完毕吗?不能设定参数后,高斯计算几分钟,然后 ...

显然不行
思路清奇
tianpingpei 发表于 Post on 2018-9-13 14:05:44
sob老师,我下载最新的Multiwfn,计算高斯得到fchk文件,显示the cubegen can not be found,我的高斯是09w版本的,然后我尝试把settings_ini文件中的您的高斯位置在d盘,我改成我的高斯安装位置c盘,还是显示thecubegen can not be found,是什么原因无法调用cubegen,是高斯版本的原因还是设置原因?
tianpingpei 发表于 Post on 2018-9-13 11:46:39
sob老师,请问一下,高斯处理chk文件的时候,必须等高斯计算完毕吗?不能设定参数后,高斯计算几分钟,然后强行终止,利用产生的chk文件转化为fch文件后在Multiwfn中计算吗?
chunlinxxx 发表于 Post on 2018-9-13 10:13:39
太感谢了,真的是快了
sobereva 发表于 Post on 2018-8-30 02:48:15
今日又更新了Multiwfn,同步地更新了本帖。settings.ini里增加了cubegendenstype参数,可以设定cubegen利用.fch里哪个级别的密度矩阵来计算静电势。而且还支持了借用cubegen计算TrEsp原子跃迁电荷,用此电荷可以快速而且校准地计算激子耦合能。

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

GMT+8, 2024-11-27 11:17 , Processed in 3.949099 second(s), 27 queries , Gzip On.

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