24小时热门版块排行榜    

查看: 2248  |  回复: 5

lsc84

铁杆木虫 (职业作家)

[求助] ORCA计算得到gbw文件后,怎么作自旋密度图?

大家好!
     我用ORCA计算得到了gbw文件,现在需要用molekel4.3作自旋密度图。如果使用orca_plot MyGBWFile.gbw -i命令只能作分子轨道图,没有其他的备选项了。如果在计算的输入文件中输入
%plots Format Cube
       MO("abc.cube",4,0) ;
end
          可以正常计算,并可以得到一个分子轨道的Gaussian_cube文件。
      相应地在输入文件中输入
%plots Format Gaussian_Cube
       SpinDens(“Cu2-SpinDens”) ;
end
          出现的结果是
Done building standard basis
LEAVING ORCA
        小弟尝试了很多次,也不能正常计算。
    请高手帮忙看看是怎么回事?
    非常感谢!
    如果采用简便输入orca_plot gbw-file plot-inputfile,plot-inputfile文件该怎么编写,我按说明书上的%plots关键字输入也不行。期待高手的解答!非常感谢!
    我想要作的自旋密度图如图所示:
回复此楼

» 本帖@通知

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsc84

铁杆木虫 (职业作家)

自旋密度图如下
2楼2011-10-01 10:02:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsc84

铁杆木虫 (职业作家)

图片直接上传

3楼2011-10-01 10:03:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

★ ★ ★ ★ ★ ★ ★ ★
ben_ladeng(金币+4): 谢谢应助,今后注意点应助回帖哦。应楼主要求奖励5金币。1个税收 2011-10-02 00:32:32
ben_ladeng(金币+4): 节日应助,奖励加倍。专家考核存档 2011-10-02 00:33:11
原则上,可以先产生molden格式的文件,然后用gabedit打开,显示自旋密度。

但是,orca的molden文件有些问题,基组部分的收缩因子都乘上了的高斯指数,对这个bug作者坚决不改,不知道目的是什么。个人认为还是应该改为正常的基组。从https://bse.pnl.gov/bse/portal 找到gaussian格式的基组,复制到molden文件,并做适当的排版。
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
4楼2011-10-01 22:53:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsc84

铁杆木虫 (职业作家)

★ ★ ★ ★ ★
ben_ladeng(金币-5): 扣除5个金币奖励。您只需要发放15个金币。剩余的5个可以请斑竹退还。谢谢理解 2011-10-02 00:35:08
beefly ,您好!
     我采用这种方法可以做出自旋密度图,首先在命令提示符下输入orca_2mkl gbw file name -mkl,先创建或生成一个.mkl文件,再用molekel4.3打开,先计算spin density,然后再更改cutoff的数值,这样也可以得到自旋密度图。
     但是,这种方法得到的自旋密度图,在前期处理时不能控制格点的数目,做出来的图很粗糙。如果计算时采用自定义基组或从Frank Neese的文献支持信息里得到的镧系元素基组计算时,采用这种方法作自旋密度图时,出现错误为:can't read basis-set.从而使计算spin density那一步不能进行,得不到自旋密度图。
     我想能否在input file里输入%plots 控制格点数目,使得到的图好一点。
     非常感谢您的关注!同时也感谢大家关注此帖!
5楼2011-10-02 00:04:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

【答案】应助回帖

★ ★ ★ ★ ★
lsc84(金币+8): 3Q 2011-10-02 13:14:35
小红豆(金币+5): 谢谢 2011-10-02 14:32:48
lsc84(金币+4): 3q 2011-10-05 09:27:08
还有一个办法

1. 下一个molden2aim最新版本(2.0),
http://muchong.com/bbs/viewthread.php?tid=1683682

2. orca生成molden文件。在文件里加入一行:
[Program] Orca
然后用molden2aim转化为wfn格式

3. 用multiwfn打开生成的wfn文件,显示自旋密度
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
6楼2011-10-02 10:31:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lsc84 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:50 by 5BDX0d0WFp7t
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 05:50 by 5BDX0d0WFp7t
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 3n8v2C8RimXI 2026-09-13 3/150 2026-09-15 05:49 by 5BDX0d0WFp7t
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-15 05:34 by 5BDX0d0WFp7t
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 05:13 by 5BDX0d0WFp7t
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-15 05:00 by mibUvS8DDCwf
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 00:12 by vZfe6xYu34yj
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-14 23:57 by vZfe6xYu34yj
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-14 23:13 by vZfe6xYu34yj
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-14 22:45 by vZfe6xYu34yj
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-14 22:33 by vZfe6xYu34yj
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 22:25 by vZfe6xYu34yj
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:36 by vZfe6xYu34yj
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:25 by vZfe6xYu34yj
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +4 两三岁ss 2026-09-14 4/200 2026-09-14 20:29 by 班主任的同桌
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-14 17:48 by Aj1rhIDL5ixY
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
信息提示
请填处理意见