【调剂】北京石油化工学院2024年16个专业接受调剂
查看: 1743  |  回复: 6
【悬赏金币】回答本帖问题,作者gx730将赠送您 3 个金币

gx730

金虫 (正式写手)

[求助] 运行scan,进行势能面的扫面,扫面键长变化很正常,但是扫描键角变化总是不能成功已有1人参与

贴出输入文件
%chk=vac ne scan.chk

# scan ub3lyp/6-31g(d,p)

vac ne scan
0 1
C                 
H    1   B1   
H    1   B2    2    A1
C    1   B3    2    A2    3    D1   
H    4   B4    1    A3    2    D2
O    4   B5    1    A4    2    D3
C    6   B6    4    A5    1    D4
O    7   B7    6    A6    4    D5
C    7   B8    6    A7    4    D6
H    9   B9    7    A8    6    D7
H    9   B10   7    A9    6    D8
H    9   B11   7    A10   6    D9
Ne   1   B12   4    A11   6    D10

    B1    1.0820752
    B2    1.0831100
    B3    1.3284079
    B4    1.0839853
    B5    1.3780400
    B6    1.3735799
    B7    1.2037515
    B8    1.5145735
    B9    1.0945255
    B10   1.0945363
    B11   1.0887872
    B12   1.1734602            100           0.1
    A1    118.4864662
    A2    119.9374459
    A3    122.9867885
    A4    120.3105814
    A5    123.6789028
    A6    117.8036106
    A7    118.4460723
    A8    111.5142463
    A9    111.5116736
    A10   107.4392175
    A11   178.9523219         3         20
    D1   -179.9947461
    D2    0.0026112
    D3    179.9770255
    D4    179.8878673
    D5    179.9727032
    D6   -0.0222308
    D7   -60.3265962
    D8    60.4113397
    D9   -179.9626893
    D10  -179.5775593
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gx730

金虫 (正式写手)

额  我好像发现错误了
是不是角度后面的步长也要是浮点数,要带个小数点
2楼2013-12-24 11:12:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gx730

金虫 (正式写手)

但是奇怪的是 我今天一运行就出错误
上个星期投的任务却运行了四天才出错误
3楼2013-12-24 11:13:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gx730

金虫 (正式写手)

势能面扫面是想拟合L-J参数
这些键长的步长,角度的步长设置是不是要先大概优化出能量最低额点
在这个点附近进行扫面,步长设置小点。
然后在趋于平衡的那段步长设置长点
4楼2013-12-24 12:02:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
gx730: 金币+2, ★★★很有帮助, 谢谢您的帮助 2013-12-24 17:36:55
根据你的输入文件来看,我觉得问题可能有以下几点,说的不对望指正:
1、你键角是走3步每步20°吗?如果是20步每步3°,那就是写反了。
2、每步步长应该是要浮点数的。
3、高斯的键角不能大于180°,二面角没有限制。
4、如果你这样的设置导致一些基本错误,比如同一个原子与周围的同一平面内原子形成的键角之和大于360°,也会马上报错。
一个小窍门是可以先用GV打开看看,如果不能打开,就是结构有问题,GV会报错第几行有问题,然后去看相关的数据。如果能打开,说明初始结构是没有问题的,那肯能是在扫描到一定步数后某个数值,比如键角大于180°,导致的报错。
5楼2013-12-24 17:06:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gx730

金虫 (正式写手)

引用回帖:
5楼: Originally posted by 枪下游魂 at 2013-12-24 17:06:42
根据你的输入文件来看,我觉得问题可能有以下几点,说的不对望指正:
1、你键角是走3步每步20°吗?如果是20步每步3°,那就是写反了。
2、每步步长应该是要浮点数的。
3、高斯的键角不能大于180°,二面角没有限 ...

运行就是这样中断了,output文件也没有显示错误,直接中断了
从result里看scan已经显示出来势能图了,但是就是输出文件里没有数据
运行scan,进行势能面的扫面,扫面键长变化很正常,但是扫描键角变化总是不能成功
高斯出错.jpg

6楼2013-12-24 17:36:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

引用回帖:
6楼: Originally posted by gx730 at 2013-12-24 17:36:15
运行就是这样中断了,output文件也没有显示错误,直接中断了
从result里看scan已经显示出来势能图了,但是就是输出文件里没有数据

高斯出错.jpg
...

回溯整个out文件,看下是不是每步都计算了,只是最后的势能曲线数据统计没有出来?我看不清你这个是在第101步的时候出错的么?如果是的话,可能就是键角的问题了,因为你每100步键角才加一次。还有那个startlink我也不清楚是哪个link,等待高人解释~
7楼2013-12-25 08:37:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gx730 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 研二光催化6月底四篇二区什么水平 5+3 wjtab 2024-04-22 7/350 2024-04-23 13:09 by 木炭萨姆
[教师之家] 发现督导在听课,需要跟TA招呼示意吗?有同学看手机、课件有疏漏算大问题吗? +10 河西夜郎 2024-04-19 10/500 2024-04-23 10:44 by Quakerbird
[教师之家] 国内科研环境差导致国内科研水平差的根源 +8 zju2000 2024-04-18 8/400 2024-04-23 10:40 by Quakerbird
[论文投稿] 编辑返稿让改格式,这个时候能大修内容吗? +6 双倍好运锦鲤 2024-04-17 9/450 2024-04-23 09:49 by fslaser
[考研] 0854-0855调剂 +4 shangannum1 2024-04-21 7/350 2024-04-23 09:46 by alkb
[基金申请] 前同事拿我以前未中的国自然标书申请今年的国自然,现在基金委发来调查函 +26 bjdxyxy 2024-04-19 40/2000 2024-04-23 09:38 by bjdlqz
[教师之家] 刚博士毕业,入职普通一本高校 +16 保持肌肤 2024-04-20 16/800 2024-04-23 08:31 by zxf984
[基金申请] 下雨了 +15 zju2000 2024-04-16 22/1100 2024-04-22 16:03 by newfuzzy1
[考博] 青岛科技大学顾玉超教授课题组招收博士生 +3 yanqingwang 2024-04-19 3/150 2024-04-22 15:21 by 理科生yes
[公派出国] 東京大学 应用化学专业柳田研究室 招收修士/博士生 +4 蝶梦123 2024-04-19 4/200 2024-04-22 15:10 by liukang88j
[考博] 申请24博士 材料/化工/环境 +5 满目_星辰 2024-04-17 6/300 2024-04-22 13:17 by 雄关不惧
[考博] 博士申请 +4 zhangxm0525 2024-04-17 11/550 2024-04-22 12:00 by zhangxm0525
[考博] 重庆大学计算机学院2024博士招生 +3 y_y_y_yyy 2024-04-18 5/250 2024-04-21 08:29 by LYS1200
[公派出国] CSC在职申请求问 +8 Amber0919 2024-04-16 8/400 2024-04-20 18:33 by 371122gyd
[考研] 312求调剂 +3 Leroic000 2024-04-19 3/150 2024-04-20 18:19 by xinaishan
[考研] 本科211 277分求调剂 +3 0854通信工程 2024-04-18 5/250 2024-04-20 17:55 by 0854通信工程
[考博] 24/25读博求博导 +4 宝23 2024-04-16 4/200 2024-04-19 17:18 by wangzhe_bs
[论文投稿] CCS Chemistry投稿求助 30+3 wfqtriumph 2024-04-17 4/200 2024-04-19 08:56 by Bletilla
[论文投稿] 发个文章,结果被某学报请的外审PUA了 +3 ziniu168 2024-04-18 8/400 2024-04-18 21:40 by xli1984
[论文投稿] 投稿求助 5+3 我是洲洲啊 2024-04-17 5/250 2024-04-18 17:13 by topedit
信息提示
请填处理意见