应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Turbomole/ORCA » ORCA计算得到gbw文件后，怎么作自旋密度图？
61/1返回列表
查看: 1701  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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关键字输入也不行。期待高手的解答！非常感谢！
    我想要作的自旋密度图如图所示：
回复此楼

» 本帖@通知

@beefly @boylc789 @lihb734 @tanghao86 @yjcmwgk

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-10-01 09:59:01
已阅   回复此楼   关注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文件，并做适当的排版。
赞一下(2人) 回复此楼
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 控制格点数目，使得到的图好一点。
     非常感谢您的关注！同时也感谢大家关注此帖！
赞一下(1人) 回复此楼
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文件，显示自旋密度
赞一下(1人) 回复此楼
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
6楼2011-10-02 10:31:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lsc84 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭