应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 【求助】单个氧原子的能量E出现问题，请高手指教！
191/1返回列表
查看: 3603  |  回复: 18
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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次，请高手给分析一下
回复此楼

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

第一性原理

» 猜你喜欢

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

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

查看全部散金贴
1楼 2010-12-08 09:46:46
已阅   回复此楼   关注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:


为什么一定要勾自旋？

单原子 由于对称性被破坏，必须要考虑电子的自旋极化，才能得到准确的基态能量。
赞一下(2人) 回复此楼
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人) 回复此楼
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点就好了
赞一下(1人) 回复此楼
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的回帖

zzy870720z

荣誉版主 (文坛精英)
★ ★
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   ...

不用管那个空间群吧
赞一下(2人) 回复此楼
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的回帖

zzy870720z

荣誉版主 (文坛精英)
★
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的回帖

zzy870720z

荣誉版主 (文坛精英)
★
sunyang1988(金币+1):谢谢交流 2010-12-08 14:31:21
引用回帖:
Originally posted by liang7035 at 2010-12-08 13:56:42:

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

应该不会，单点能应该会很快的
赞一下(1人) 回复此楼
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:

必须要勾。

为什么一定要勾自旋？
赞一下(1人) 回复此楼
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 的主题更新
191/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭