24小时热门版块排行榜    

查看: 4579  |  回复: 12
【悬赏金币】回答本帖问题,作者天波稀客将赠送您 45 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

天波稀客

金虫 (小有名气)

[求助] 求助高斯09计算中的stable与stable=opt的用法

本人最近在研究双自由基问题,高斯手册中提到臭氧分子是典型的双自由基分子(两个成单电子自旋相反),我用"UCCSD(T)/aug-cc-pVTZ guess=mix stable"命令在高斯09中检测了它电荷为0自选多重度为1的波函数稳定性(已在相同计算水平下进行构型优化),结果显示The wavefunction has an internal instability,我的哪里出现问题了呢?然后我又用"UCCSD(T)/aug-cc-pVTZ guess=mix stable"检测过它电荷为0自选多重度为3的波函数,结果也是The wavefunction has an internal instability,用"RCCSD(T)/aug-cc-pVTZ guess=mix stable"检测过它电荷为0自选多重度为1的波函数,结果显示The wavefunction has an RHF -> UHF instability。我现在找不到任何一个stable检测结果稳定的波函数,所以想不通。Sobereva老师曾给我指点用stable=opt命令,于是我用UCCSD(T)/aug-cc-pVTZ Stable=opt对电荷为0自选多重度为1的臭氧进行计算,结果提示出错Error termination in NtrErr: NtrErr Called from FileIO. 请懂这方面知识的老师帮忙指点一下,stable与stable=opt的用法,非常感谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天波稀客

金虫 (小有名气)

送红花一朵
引用回帖:
8楼: Originally posted by paramecium86 at 2017-07-28 21:40:13
因为stable只检测 HF 或者 DFT 波函数的稳定性。 所以即使 写成 RCCSD(T) stable 实际上检测的还是 RHF wavefunction 检查完之后直接就结束任务 不进入post-HF部分 所以不会报错。 但是STABLE=opt的时候 检测完一次 ...

非常感谢您的耐心回复,我的问题都被您点拨的很清楚了,谢谢!
您比我高出好多个level,您就是老师
9楼2017-07-28 22:24:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

paramecium86

版主 (著名写手)

★ ★ ★ ★ ★
zhou2009: 金币+5 2017-07-26 14:48:14
Stable命令就是检查一下 波函数稳定不稳定
stable=opt就是检查一下 稳定就okay 不稳定就opt到稳定。
算自旋极化的单重态 波函数不稳定挺正常的。

不过这个错误一般是硬盘不够了才会出现的(查查输出文件里是不是还有 这句话 FileIO operation on non-existent file.)  如果有的话可以检查一下checkpoint所存放的磁盘是否够大? 毕竟ccsdt这个级别的计算所用的硬盘量不小。

另外除了直接算整体也可以尝试在分子里设置片段 两边分别是三重态的 alpha和beta电子,这个如果你看sob的博客上写过一个讲这个的可以找找。

» 本帖已获得的红花(最新10朵)

2楼2017-07-26 11:56:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天波稀客

金虫 (小有名气)

送红花一朵
引用回帖:
2楼: Originally posted by paramecium86 at 2017-07-26 11:56:47
Stable命令就是检查一下 波函数稳定不稳定
stable=opt就是检查一下 稳定就okay 不稳定就opt到稳定。
算自旋极化的单重态 波函数不稳定挺正常的。

不过这个错误一般是硬盘不够了才会出现的(查查输出文件里是不 ...

首先谢谢您的回答!但是我对高斯手册中提到臭氧分子是典型的双自由基分子(两个成单电子自旋相反),用"UCCSD(T)/aug-cc-pVTZ guess=mix stable"命令在高斯09中检测了它电荷为0自选多重度为1的波函数稳定性,结果显示The wavefunction has an internal instability,还是非常不解。另外想问一下您,臭氧分子是不是双自由基分子呢?双自由基分子都是单重态的吗?如果您了解这方面知识拜托相告,非常感谢!至于那个Error termination in NtrErr: NtrErr Called from FileIO的出错,文件里并没有FileIO operation on non-existent file这句话 ,应该不会是磁盘内存不足问题引起的出错,文件尾部信息见下文:
Error termination in NtrErr:
NtrErr Called from FileIO.
Error: segmentation violation
   rax 0000000000000000, rbx 00000000012d74c0, rcx ffffffffffffffff
   rdx 0000000000013a8c, rsp 00007ffc177bef68, rbp 00007ffc177bef90
   rsi 000000000000000b, rdi 0000000000013a8c, r8  00002b58abdd72c0
   r9  0000000000000000, r10 00007ffc177becf0, r11 0000000000000202
   r12 00007ffc177bf138, r13 00007ffc177bf140, r14 000000000131cb48
   r15 00002b58acab3350
  --- traceback not available


-----Phoenix cluster job statistics-----
Time info:

real        0m4.480s
user        1m6.986s
sys        0m1.249s

Disk space usage:
15M        /scratch/2439851.phoenix.chem.wisc.edu
3楼2017-07-26 22:16:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)

引用回帖:
3楼: Originally posted by 天波稀客 at 2017-07-26 22:16:06
首先谢谢您的回答!但是我对高斯手册中提到臭氧分子是典型的双自由基分子(两个成单电子自旋相反),用"UCCSD(T)/aug-cc-pVTZ guess=mix stable"命令在高斯09中检测了它电荷为0自选多重度为1的波函数稳定性 ...

input和output 上传上来看看 可能能帮助确定一下错误来源。
4楼2017-07-27 09:14:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[高分子] 寻找聚酯反应釜 +3 茕茕恭煮 2024-06-15 6/300 2024-06-18 14:15 by 茕茕恭煮
[硕博家园] 关于硕博连读的一些疑问? +5 Lwenter 2024-06-14 6/300 2024-06-18 12:49 by 博csc分享
[高分子] 烧瓶内合成聚酯 +3 大帝国乐 2024-06-17 7/350 2024-06-18 11:10 by 大帝国乐
[硕博家园] 豫北虫友互识 +4 xuhongli903 2024-06-18 4/200 2024-06-18 10:26 by 飘逸的丁香
[基金申请] 博后基金刷到的BUG,图片来的更直观 +12 carolloo 2024-06-17 13/650 2024-06-18 08:35 by q478652742
[基金申请] F口信息学部拿面上,大概需要什么样的成果 +7 _奋黎_ 2024-06-16 12/600 2024-06-18 07:41 by ld_cauchy
[文学芳草园] 累并快乐着 +15 MYHLD521 2024-06-14 15/750 2024-06-18 01:28 by liyeqik
[论文投稿] 审稿问题:为什么荧光激发波长和紫外吸收波长差的大? 10+5 sdawege 2024-06-14 10/500 2024-06-17 18:54 by HH-探针
[有机交流] 苯酚和炔醇反应,mitsunobu反应 5+3 TONGMEIMEI 2024-06-12 9/450 2024-06-17 15:00 by dschong
[基金申请] 博后基金,博管会会提前知道消息吗? +4 yuyiang 2024-06-13 4/200 2024-06-16 11:40 by yangyuzhong4
[基金申请] 为什么我的博后基金还在流动站审核中?不会是学院给我卡了吧? +14 王凯12 2024-06-13 26/1300 2024-06-15 15:22 by 好人与坏人
[论文投稿] 求机械类四区sci推荐 5+3 迷茫小旷 2024-06-14 4/200 2024-06-15 11:25 by bobvan
[基金申请] 博后基金,以往的结果点不开,怎么回事呢?最后一次机会了,两次都没中前面。 +7 kyukitu 2024-06-14 13/650 2024-06-15 06:46 by 我是王小帅
[论文投稿] 投了一篇4区的SCI,审稿人一个拒稿,一个小修,编辑给了大修。 +9 安稳22123 2024-06-13 10/500 2024-06-14 23:45 by jurkat.1640
[基金申请] 博士后基金需要结题吗? +8 zhouchuck 2024-06-13 8/400 2024-06-14 17:27 by liuyupu132
[考博] 申博找导师 +4 疏影横斜水清浅3 2024-06-13 6/300 2024-06-14 14:31 by zxl_1105
[有机交流] ππ堆积会发生在有机溶剂中吗 5+3 zibuyu0420 2024-06-13 4/200 2024-06-14 14:17 by 小肉干
[基金申请] 75批博后基金 +10 kyukitu 2024-06-13 13/650 2024-06-14 10:31 by kyukitu
[基金申请] 工材E10口函评结束了吗 10+3 我1的飞翔 2024-06-13 5/250 2024-06-14 06:35 by nono2009
[论文投稿] 投稿后发现其他作者的邮箱填错了该怎么办呀 10+4 在飞的猪 2024-06-13 6/300 2024-06-14 04:45 by 小虫子咔咔
信息提示
请填处理意见