24小时热门版块排行榜    

查看: 4584  |  回复: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 豫北虫友互识 +10 xuhongli903 2024-06-18 11/550 2024-06-20 00:25 by xuhongli903
[海外博后] 墨尔本大学博后offer要不要接 +3 kyxblmm 2024-06-18 3/150 2024-06-19 22:39 by blake1111
[基金申请] F口信息学部拿面上,大概需要什么样的成果 +7 _奋黎_ 2024-06-16 18/900 2024-06-19 20:02 by _奋黎_
[基金申请] 青年基金会评专家到底是怎么会评的呀?主审专家是不是一般不会改动系统按函评给的顺序 5+3 他山攻玉之石 2024-06-18 14/700 2024-06-19 19:43 by liliwang215
[考博] 关于读博感觉自己很抓马 +8 小九月 2024-06-19 11/550 2024-06-19 19:22 by 鱼翔浅底1
[基金申请] 75批博后基金 +11 kyukitu 2024-06-13 14/700 2024-06-19 18:44 by Siuhong Yiu
[基金申请] 基金得中 +4 woaini0218 2024-06-18 4/200 2024-06-19 17:27 by mengzl
[基金申请] Nature 11日发文,中国著名学者们称造假迫不得已 +8 babu2015 2024-06-14 8/400 2024-06-19 15:25 by 风今25
[论文投稿] 求机械类四区sci推荐 5+4 迷茫小旷 2024-06-14 5/250 2024-06-19 14:08 by tangjie12345
[找工作] 杭电、天津科技、青农和宁波工程学院如何选? +12 味道很好啊 2024-06-13 24/1200 2024-06-18 23:09 by 药材
[教师之家] 请问事业编制和年薪制冲突吗? +12 ZHONGWU_U 2024-06-14 12/600 2024-06-18 19:31 by fangyl2005
[考博] 2025考博 +8 自强不息a?a 2024-06-15 13/650 2024-06-18 18:12 by 投必得科研顾问
[有机交流] 跑板能跑开,过柱过不纯怎么办 +4 小胡在努力 2024-06-18 6/300 2024-06-18 15:20 by long8811
[基金申请] 太卷了 +13 laoyuefubio 2024-06-17 26/1300 2024-06-18 10:56 by shuigubio
[有机交流] 求教最后一步问题出在哪 20+3 苏州浩恒 2024-06-15 4/200 2024-06-18 08:33 by like_2005
[论文投稿] 投稿被一个审稿人恶意评审了怎么样? +5 1chen 2024-06-14 7/350 2024-06-15 23:15 by xy66xy
[基金申请] 关于博后基金的bug问题 +6 lxr1991 2024-06-14 9/450 2024-06-15 21:17 by since—2010
[论文投稿] 投稿时忘记修改一作 +7 gll123456 2024-06-13 11/550 2024-06-15 11:49 by gll123456
[基金申请] E12面上申请 +4 汉风之遗 2024-06-13 4/200 2024-06-14 15:28 by 天外飞去来
[基金申请] 博士后面上项目状态还是专家评审吗 10+9 Thatcheremu 2024-06-13 55/2750 2024-06-13 21:23 by 乌合麒麟
信息提示
请填处理意见