24小时热门版块排行榜    

查看: 3636  |  回复: 18

liang7035

金虫 (正式写手)


[交流] 【求助】单个氧原子的能量E出现问题,请高手指教!

要计算单个氧原子的能量E,有高手提出建议:”在你原来的计算条件下,建一个10*10*10Angstrom的盒子,里面放一个原子,计算energy,就出来了“ 我建好晶胞了,计算总是失败,出现的情况如下:


*Warning* max. SCF cycles performed but system has not reached the groundstate.
  
Current total energy, E     =  -429.2069682116     eV
Current free energy (E-TS)  =  -429.3123849790     eV
(energies not corrected for finite basis set)

NB est. 0K energy (E-0.5TS)      =  -429.2596765953     eV

  ****************************************************************************
  Warning: electronic minimisation did not converge when finding ground state.
  ****************************************************************************
  Writing checkpoint file...

Writing model to 3D_Atomistic.check
Error check_elec_ground_state : electronic_minimisation of initial cell failed.
Current trace stack:
check_elec_ground_state
castep
Error check_elec_ground_state : electronic_minimisation of initial cell failed.
Current trace stack:
check_elec_ground_state
castep
MPI Application rank 1 exited before MPI_Finalize() with status 1

除了是循环次数不够,还有什么原因啊,我设置是1000次,请高手给分析一下
回复此楼

» 收录本帖的淘帖专辑推荐

第一性原理

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

fanchen021

铁杆木虫 (正式写手)


★ ★
zzy870720z(金币+2):谢谢指教 2010-12-09 10:34:58
引用回帖:
Originally posted by ym23 at 2010-12-08 16:47:00:


为什么一定要勾自旋?

单原子 由于对称性被破坏,必须要考虑电子的自旋极化,才能得到准确的基态能量。
18楼2010-12-09 09:37:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

fanchen021

铁杆木虫 (正式写手)


★ ★
liang7035(金币+1):谢谢参与
zzy870720z(金币+1): 鼓励交流 2011-03-12 16:25:09
Spin polarized 勾了没?
2楼2010-12-08 10:02:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liang7035

金虫 (正式写手)


引用回帖:
Originally posted by fanchen021 at 2010-12-08 10:02:45:
Spin polarized 勾了没?

没有啊?高手
3楼2010-12-08 10:04:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanchen021

铁杆木虫 (正式写手)



sunyang1988(金币+1):谢谢指点 2010-12-08 10:11:25
引用回帖:
Originally posted by liang7035 at 2010-12-08 10:04:57:

没有啊?高手

必须要勾。
4楼2010-12-08 10:06:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liang7035

金虫 (正式写手)


引用回帖:
Originally posted by fanchen021 at 2010-12-08 10:06:32:

必须要勾。

谢谢高手,麻烦你在给看看其他设置是否有问题谢谢


5楼2010-12-08 10:11:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanchen021

铁杆木虫 (正式写手)



zzy870720z(金币+1):谢谢指教,这个很重要,呵呵 2010-12-08 12:02:59
算单原子能量,K点取Gamma点就好了
6楼2010-12-08 10:22:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liang7035

金虫 (正式写手)


好的,谢谢
7楼2010-12-08 10:46:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liang7035

金虫 (正式写手)


引用回帖:
Originally posted by fanchen021 at 2010-12-08 10:22:24:
算单原子能量,K点取Gamma点就好了

高手还得问你一下,计算单个氧原子,那么空间群按O1算呗?

*data for   ICSD #92775
Coll Code   92775
Rec  Date   2003/04/01
Chem Name   Oxygen
Structured  O
Sum         O1
ANX         N
D(calc)     .59
Title       Modulated structure of the pseudohexagonal In Fe1-x-4d Tix+3d O3+x/2
            (x=0.61) composite crystal
Author(s)   Michiue, Y.;Onoda, M.;Watanabe, M.;Brown, F.;Kimizuka, N.
Reference   Acta Crystallographica B (39,1983-)
            (2001), 57, 458-465
Unit Cell   5.835(3) 2.568(6) 12.082(7) 90 90 90
Vol         181.04
Z           4
Space Group C c m mSG Number   63
Cryst Sys   orthorhombic
Pearson     oS4
Wyckoff     c
R Value     .03
Red Cell    C  2.568 3.187 12.082 90 89.999 113.754 90.52
Trans Red   0.000 -1.000 0.000 / 0.500 0.500 0.000 / 0.000 0.000 1.000
Comments    Modulated structure first part structure: collection no.
            92774
Atom  #   OX   SITE      x           y           z           SOF      H
O    1  +0    4 c   -0.171(2)   0.          0.25           1.         0   
*end for    ICSD #92775
8楼2010-12-08 11:19:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
liang7035(金币+1):谢谢参与
sunyang1988(金币+1): 谢谢交流 2011-03-12 18:43:51
引用回帖:
Originally posted by liang7035 at 2010-12-08 11:19:59:

高手还得问你一下,计算单个氧原子,那么空间群按O1算呗?

*data for   ICSD #92775
Coll Code   92775
Rec  Date   2003/04/01
Chem Name   ...

不用管那个空间群吧
9楼2010-12-08 12:03:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liang7035

金虫 (正式写手)


引用回帖:
Originally posted by zzy870720z at 2010-12-08 12:03:50:


不用管那个空间群吧

直接是P1吗?
10楼2010-12-08 13:03:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988(金币+1):谢谢指点 2010-12-08 14:31:14
引用回帖:
Originally posted by liang7035 at 2010-12-08 13:03:40:

直接是P1吗?

可以,没多大关系
11楼2010-12-08 13:27:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liang7035

金虫 (正式写手)


引用回帖:
Originally posted by zzy870720z at 2010-12-08 13:27:10:


可以,没多大关系

谢谢,在问一下,你看我估计这个计算应该很快,但是我一上午没出结果!
12楼2010-12-08 13:56:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988(金币+1):谢谢交流 2010-12-08 14:31:21
引用回帖:
Originally posted by liang7035 at 2010-12-08 13:56:42:

谢谢,在问一下,你看我估计这个计算应该很快,但是我一上午没出结果!

应该不会,单点能应该会很快的
13楼2010-12-08 14:28:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liang7035

金虫 (正式写手)


那你说原因可能出现在哪呢?
我把参数设置都发在上面了啊!拜托给看看谢谢
14楼2010-12-08 15:28:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liang7035

金虫 (正式写手)


引用回帖:
Originally posted by zzy870720z at 2010-12-08 14:28:16:


应该不会,单点能应该会很快的

那你说原因可能出现在哪呢?
我把参数设置都发在上面了啊!拜托给看看谢谢
本文来自: 小木虫论坛 http://muchong.com/bbs/viewthread.php?tid=2677342&pid=706545&page=2#pid706545
15楼2010-12-08 15:29:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ym23

铁杆木虫 (著名写手)



liang7035(金币+1):谢谢参与
引用回帖:
Originally posted by fanchen021 at 2010-12-08 10:06:32:

必须要勾。

为什么一定要勾自旋?
16楼2010-12-08 16:47:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ym23

铁杆木虫 (著名写手)


引用回帖:
Originally posted by fanchen021 at 2010-12-08 10:22:24:
算单原子能量,K点取Gamma点就好了

为什么直接取Gamma点就够了?
17楼2010-12-08 16:47:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hiluokun

禁虫 (小有名气)


liang7035(金币+1):谢谢参与
本帖内容被屏蔽

19楼2011-03-12 15:54:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 liang7035 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
[考研] 求调剂 +4 Mqqqqqq 2026-03-19 4/200 2026-03-20 14:15 by 星空星月
[考研] 材料学硕297已过四六级求调剂推荐 +6 adaie 2026-03-19 6/300 2026-03-20 14:12 by 星空星月
[考研] 一志愿中南化学337求调剂 +4 niko- 2026-03-19 5/250 2026-03-20 13:49 by 促天成
[考研] 317求调剂 +4 申子申申 2026-03-19 8/400 2026-03-20 11:20 by 申子申申
[考研] 0703化学调剂 +10 妮妮ninicgb 2026-03-15 14/700 2026-03-19 22:59 by 学员8dgXkO
[考研] 288求调剂 +15 于海海海海 2026-03-19 15/750 2026-03-19 22:41 by 学员8dgXkO
[考博] 申博26年 +3 八6八68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +3 葵梓卫队 2026-03-18 5/250 2026-03-19 19:35 by 给你你注意休息
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 6/300 2026-03-19 17:18 by fei626-918
[考博] 东华理工大学化材专业26届硕士博士申请 +8 zlingli 2026-03-13 8/400 2026-03-19 16:32 by 轻松不少随
[考研] 一志愿北京化工大学0703化学318分,有科研经历,求调剂 +3 一瓶苯甲酸 2026-03-14 3/150 2026-03-19 15:17 by 尽舜尧1
[考研] 求调剂,一志愿:南京航空航天大学大学 ,080500材料科学与工程学硕,总分289分 +3 @taotao 2026-03-19 3/150 2026-03-19 14:07 by peike
[考研] 0703化学 305求调剂 +4 FY_yy 2026-03-14 4/200 2026-03-19 05:54 by anny19840123
[考研] 304求调剂 +6 司空. 2026-03-18 6/300 2026-03-18 23:03 by 星空星月
[考研] 354求调剂 +4 Tyoumou 2026-03-18 7/350 2026-03-18 21:45 by Tyoumou
[考研] 一志愿武理材料305分求调剂 +5 想上岸的鲤鱼 2026-03-18 6/300 2026-03-18 17:53 by 无际的草原
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[考研] 0703化学调剂 290分有科研经历,论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
信息提示
请填处理意见