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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 用ONIOM关键词算吸收，请问可以看分子轨道吗？
51/1返回列表
查看: 4529  |  回复: 21
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者youxx189将赠送您 10 个金币
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

youxx189

铜虫 (小有名气)

[求助] 用ONIOM关键词算吸收，请问可以看分子轨道吗？ 已有1人参与

用ONIOM关键词算吸收，请问可以看分子轨道吗？得到的chk文件转化成fchk文件会造成数据丢失，在Gaussview里无法打开。


输入文件：

%chk=packing1_guest_abs_ONIOM.chk
%nprocshared=20
%mem=12000MB
#p oniom(PBE1PBE/genecp td=(singlets,nstates=20,direct):uff)=embedcharge  test

absorption

0 1 0 1 0 1
C        26.8564499443,11.2974579201,22.535374973        L
C        26.4312706227,12.5897273583,22.1656391631       L
C        25.079756231,12.9545921491,22.2751903325        L
C        24.1694372817,11.9921553489,22.7275839214       L
C        24.5970630014,10.7317662379,23.1036871816       L
C        25.9390511173,10.3569276904,23.0333609336       L
H        27.8993300059,11.0327314821,22.4387729637       L
中间略到好多原子
C        33.1316875537,25.2458454493,20.5017776507       L
H        34.6438528841,23.5215656048,22.9626666632       L
H        34.5537458777,23.6020072384,20.5486274853       L
C        33.0912226367,25.4397359482,19.0253450849       L
O        33.7582638392,24.7213374764,18.3023322806       L
O        32.3461046073,26.4006128366,18.4698156726       L

C H N O 0
6-31G**
****
Ir 0
LANL2DZ
****

Ir 0
LANL2DZ
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2020-06-21 14:58:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)
引用回帖:
4楼: Originally posted by youxx189 at 2020-06-22 00:40:46
谢谢老师     我加了环境变量还是报错
$ module unload gauss/g09d01
$ module load gauss/g16a03
$ export GAUSS_MEMDEF=4000000000
$ formchk packing1ONIOM_dimer1-2.chk
Read checkpoint file packing1O ...

这种错误一般要么是 1. formchk命令调用的版本不是当前计算出chk的高斯使用的版本。必须完全一致才可以
                                2. chk文件本身有问题。 这个就再算一次单点能就可以了。td 的chk文件本身含有的波函数信息就是基态的，你重新计算一次单点能就可以 不用重新算一次TD这样省时间。
                               3. 高斯当前的版本有bug，换个别的版本的高斯。
赞一下(1人) 回复此楼
5楼2020-06-22 14:02:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 22 个回答

wzhsun

新虫 (小有名气)
个人觉得是你内存设的太小了

发自小木虫IOS客户端
赞一下 回复此楼
2楼2020-06-22 11:38:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)
如果你的chk文件本身很大 比如说上G大小了 那转换一般会内存不够用导致出错。可以在环境变量里加上 export GAUSS_MEMDEF=4000000000  把内存调成可用4 GB应该就能解决问题。 TD计算chk文件中的分子轨道的信息就是基态分子轨道的信息 并没有什么特殊之处

发自小木虫IOS客户端
赞一下 回复此楼
3楼2020-06-22 12:03:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

youxx189

铜虫 (小有名气)
引用回帖:
3楼: Originally posted by paramecium86 at 2020-06-22 12:03:37
如果你的chk文件本身很大 比如说上G大小了 那转换一般会内存不够用导致出错。可以在环境变量里加上 export GAUSS_MEMDEF=4000000000  把内存调成可用4 GB应该就能解决问题。 TD计算chk文件中的分子轨道的信息就是基 ...

谢谢老师     我加了环境变量还是报错
[youxx@mu01 guest-abs]$ module unload gauss/g09d01
[youxx@mu01 guest-abs]$ module load gauss/g16a03
[youxx@mu01 guest-abs]$ export GAUSS_MEMDEF=4000000000
[youxx@mu01 guest-abs]$ formchk packing1ONIOM_dimer1-2.chk
Read checkpoint file packing1ONIOM_dimer1-2.chk
Write formatted file packing1ONIOM_dimer1-2.fchk
Missing data for FChkPn.
Error termination via Lnk1e at Mon Jun 22 13:36:29 2020.
Error: segmentation violation
   rax 0000000000000000, rbx 0000000002b377e0, rcx ffffffffffffffff
   rdx 000000000001ef4c, rsp 00007ffe7319ed28, rbp 00007ffe7319ed80
   rsi 000000000000000b, rdi 000000000001ef4c, r8  00002b4517b60ac0
   r9  00002b4517b60ac0, r10 00007ffe7319e120, r11 0000000000000206
   r12 000000000000000c, r13 0000000000000001, r14 0000000000000001
   r15 0000000000000000
  /lib64/libpthread.so.0(+0xf6d0) [0x2b4517d366d0]
  /lib64/libc.so.6(kill+0x7) [0x2b451827b547]
  formchk() [0x497ec9]
  formchk() [0x46f560]
  formchk() [0x46b8af]
  formchk() [0x453755]
  formchk() [0x4516a8]
  formchk() [0x450af4]
  /lib64/libc.so.6(__libc_start_main+0xf5) [0x2b4518267445]
  formchk() [0x4509e9]
Aborted
赞一下(1人) 回复此楼
4楼2020-06-22 13:40:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 22 个回答
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭