24小时热门版块排行榜    

查看: 3630  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 298-一志愿中国农业大学-求调剂 +9 手机用户 2026-03-17 9/450 2026-03-20 14:24 by 无懈可击111
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 12/600 2026-03-20 12:25 by lidianxing
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 14/700 2026-03-20 10:46 by EBSD
[考研] 材料080500调剂求收留 +6 一颗meteor 2026-03-13 6/300 2026-03-20 10:41 by EBSD
[考研] 081700化工学硕调剂 +3 【1】 2026-03-16 3/150 2026-03-19 23:40 by edmund7
[考研] 0703化学调剂 +10 妮妮ninicgb 2026-03-15 14/700 2026-03-19 22:59 by 学员8dgXkO
[考博] 申博26年 +3 八6八68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[考研] 286求调剂 +6 lemonzzn 2026-03-16 10/500 2026-03-19 14:31 by lemonzzn
[考研] 0703化学调剂 +5 pupcoco 2026-03-17 8/400 2026-03-19 13:58 by houyaoxu
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[考研] 一志愿天大材料与化工(085600)总分338 +5 蔡大美女 2026-03-13 5/250 2026-03-19 10:44 by 是小刘呀~
[教师之家] 焦虑 +9 水冰月月野兔 2026-03-13 13/650 2026-03-19 09:50 by otani
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[考研] 301求调剂 +4 A_JiXing 2026-03-16 4/200 2026-03-17 17:32 by ruiyingmiao
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
信息提示
请填处理意见