24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知(长期有效)
查看: 3249  |  回复: 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 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见