24小时热门版块排行榜    

查看: 1705  |  回复: 10

pelargonium

新虫 (初入文坛)

[求助] ELFCAR坐标混乱,与原子结构不匹配 已有1人参与

VASP做结构优化是加了LELF=TRUE,得到ELFCAR用VESTA处理时,发现ELFCAR导出的数据和原子结构不匹配,在有的有原子的地方显示是空位,有的是空位的地方显示有原子,但是CHGCAR导出数据和原子结构是符合的。
截了几个面,发现ELFCAR的010面试CHGCAR的001面,ELF的001面是CHGCAR的100面,ELF的100面是CHGCAR错开的010面,试了几次都这样,不清楚是怎么回事。
010面,d=0.5a,图给大家看一下对比,蓝底的是ELF,绿底是CHG,还有一张是原子结构。
拜托大家告诉我这是怎么回事吧!!!!!!万分感谢。

CHGCAR.jpg



ELFCAR.jpg



STRUCTURE.jpg
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

魔鬼中的天使

金虫 (小有名气)

【答案】应助回帖

★ ★
fzx2008: 金币+2, 谢谢交流! 2013-01-19 20:05:40
ELF反映的电子的局域化程度
可以用于判断成键的情况
没有原子的地方,也有可能形成多中心键

不知道你这两种原子都是什么,以及其他层面都怎么样,不好进一步分析了

» 本帖已获得的红花(最新10朵)

2楼2013-01-19 19:57:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pelargonium

新虫 (初入文坛)

送鲜花一朵
引用回帖:
2楼: Originally posted by 魔鬼中的天使 at 2013-01-19 19:57:06
ELF反映的电子的局域化程度
可以用于判断成键的情况
没有原子的地方,也有可能形成多中心键

不知道你这两种原子都是什么,以及其他层面都怎么样,不好进一步分析了

之前文献的情况是空位处显示蓝色的,而且现在已经基本确定是ELFCAR的坐标轴与原子结构和CHGCAR相错,ELF的a,b,c轴分别为原子结构的b,c,a轴,不知道这种坐标混乱有没有办法可以解决?
3楼2013-01-20 12:28:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fzuvivi

金虫 (小有名气)

★ ★
franch: 金币+2, 谢谢回帖交流,, 2013-01-21 21:42:03
我觉得你的结构有点怪异。原子间的间隙巨大啊~
ELF最好是用弛豫后的scf得到。不知道用弛豫得到的会不会一样。
建议贴出你的INCAR和POSCAR
4楼2013-01-20 16:53:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pelargonium

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by fzuvivi at 2013-01-20 16:53:41
我觉得你的结构有点怪异。原子间的间隙巨大啊~
ELF最好是用弛豫后的scf得到。不知道用弛豫得到的会不会一样。
建议贴出你的INCAR和POSCAR

INCAR

SYSTEM =GSTGO

Start parameter for this run:
ISTART = 0
ICHARG = 2
NWRITE = 2
ISPIN = 1

Electronic Relaxation
ENCUT = 320
IALGO = 38
NELM = 60
EDIFF = 1E-05
LREAL = .TRUE.
NELMIN = 8

Ionic Relaxation
EDIFFG = -0.01   Stopping-Criterion for IOM (HF-FORCE)
NSW = 200
IBRION = 2
ISIF = 3
POTIM = 3

DOS related values:
ISMEAR = 0
SIGMA = 0.1

Write flags
LWAVE = .FALSE.
LCHARG = .TRUE.
LELF = .TRUE.


LPLANE = .TRUE.
NPAR = 16
LSCALU = .FALSE.
NSIM = 4



POSCAR

GST-225-cell                           
   1.00000000000000     
    12.2260291973393400    0.1071386669248578    0.1068990913524132
     0.1062901342114842   12.2971433711523677    0.1266777097852780
     0.1071629105308827    0.1285186331593658   12.4085723481011669
   Ge   Sb   Te
  13  13  32
Direct
  0.2519444564005303  0.2457650359981172  0.2589318176278494
  0.9927053793093025  0.9851977837201132  0.2377197587234996
  0.9847656567988778  0.2441427315650469  0.9859656537539453
  0.7532801304209262  0.2405193933916095  0.2365323037081040
  0.2402338531882204  0.7317140768334295  0.2500469243345699
  0.7585458972719255  0.5026892086543936  0.0140526432258342
  0.7488639557907306  0.7300586050060833  0.2435272347899303
  0.5124517623788239  0.4901326923991314  0.2343940508892826
  0.2464789415580683  0.9966328293509528  0.5020853638061832
  0.9896531867357761  0.9896993819502333  0.7313560950104891
  0.7578957624066532  0.2403120816158282  0.7362630466060839
  0.2467912369537892  0.7402204498240136  0.7331464874811989
  0.4853241725839526  0.7437039613269170  0.4976385054579157
  0.2446537617615977  0.0101704672063675  0.9990239284972979
  0.4954447844712491  0.9886963908846171  0.2457399807276524
  0.4779860929624358  0.2437909840228637  0.9906872509279020
  0.0222978516260149  0.4940013608414566  0.2419840098388959
  0.9901153404764922  0.7410570523987483  0.0016173329053297
  0.4879393101776366  0.7393977408814439  0.9982058770527577
  0.9876706905788493  0.2637763644961745  0.4968714751714727
  0.7549128404365534  0.9950845213535415  0.5020269979281050
  0.4806749429222791  0.0012382276111660  0.7467919525564011
  0.4783393303541885  0.2588258627011283  0.4918010438989920
  0.9989970643581659  0.4851712292950179  0.7271421997035031
  0.7527816405010023  0.7368405125863526  0.7325543658277024
  0.4817467914679128  0.4955063641435581  0.7263492414510502
  0.9995052005851364  0.0105595111571409  0.0075771143015635
  0.0111381844645109  0.2526350184622870  0.2537269248999988
  0.2529633263429800  0.0086620331218877  0.2630359439461608
  0.2449699822131942  0.2583896810417008  0.0129869751112543
  0.5100361831251892  0.0036238242955188  0.9972282484001543
  0.5093765691924282  0.2561539933374071  0.2459015795385913
  0.7611950653331003  0.0074763795033628  0.2523437702375352
  0.7585466502820453  0.2605063310871414  0.0091509414561707
  0.0080808804359646  0.5004152073476789  0.0059334469243291
  0.0002474538275418  0.7568838571405621  0.2578349987132699
  0.2604060566781884  0.5062093367456050  0.2628974700532522
  0.2500801542345424  0.7565511235223725  0.0155052129586275
  0.5129979956374742  0.5018297484045962  0.0087276696241538
  0.5012379857196617  0.7502477556259579  0.2614003317637340
  0.7484634736392360  0.5027604652114496  0.2542926995097022
  0.7502333873297276  0.7559447522281058  0.0141444212019392
  0.0068069837567921  0.9999011504486340  0.5056746910041416
  0.0092744233091103  0.2456527385910226  0.7554944601660344
  0.2472848544007048  0.0162934618380240  0.7543049508788309
  0.2495370184462199  0.2509602457234247  0.4984650384033009
  0.4981482922097248  0.0010717880027854  0.4990170422813112
  0.5094273359802668  0.2501536854750596  0.7447353012090167
  0.7598255044207077  0.0050704018401942  0.7540840647062297
  0.7509258963240159  0.2561924300758582  0.5081148019252496
  0.9784486739513182  0.5053378403855117  0.5002774584947939
  0.0062984808006569  0.7561554278386354  0.7548978855564014
  0.2510887217941662  0.5147773422790127  0.7652279773662356
  0.2564402820479537  0.7617933353737597  0.5063840621064575
  0.5193177584487135  0.5103125109590432  0.4992852729611504
  0.4954255066751551  0.7531821694545787  0.7496872726425334
  0.7540879316347565  0.4976036947879685  0.7672703110791221
  0.7556889528668411  0.7523474486354006  0.5019361186768055



这个是经过优化的晶体结构GeSbTe,有的地方间隙巨大的应该是空位
5楼2013-01-21 08:59:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pelargonium

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by fzuvivi at 2013-01-20 16:53:41
我觉得你的结构有点怪异。原子间的间隙巨大啊~
ELF最好是用弛豫后的scf得到。不知道用弛豫得到的会不会一样。
建议贴出你的INCAR和POSCAR

还有,什么叫弛豫后的SCF?
6楼2013-01-21 09:00:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fzuvivi

金虫 (小有名气)

【答案】应助回帖

引用回帖:
6楼: Originally posted by pelargonium at 2013-01-21 09:00:50
还有,什么叫弛豫后的SCF?...

我的意思是先弛豫,达到精度后,再做一次自洽,在自洽中加入LELF。然后再通过VESTA读取。看了你的POSCAR,首先佩服下你们的服务器作业能力!呵呵~。应该没什么太大的问题。可是你给出来的图形结构有点怪。。是不是对称性的不足。刚用VESTA看了下,总觉得你那个胞的第一层原子间距有点怪异

7楼2013-01-21 22:10:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

魔鬼中的天使

金虫 (小有名气)

【答案】应助回帖

引用回帖:
3楼: Originally posted by pelargonium at 2013-01-20 12:28:52
之前文献的情况是空位处显示蓝色的,而且现在已经基本确定是ELFCAR的坐标轴与原子结构和CHGCAR相错,ELF的a,b,c轴分别为原子结构的b,c,a轴,不知道这种坐标混乱有没有办法可以解决?...

坐标轴还会混乱么?
直接文本打开CHGCAR和ELFCAR,对比头上的abc看看
同一次计算的结果,应该坐标轴一样的啊
8楼2013-01-22 00:25:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

DF11G0005

新虫 (初入文坛)

我的也出现了这个问题,和楼主一样一样的。查看过了ELFCAR和CHGCAR头上的若干行,坐标和原子位置都是一样的。非常郁闷。
9楼2013-07-24 13:32:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pelargonium

新虫 (初入文坛)

我的ELF坐标混乱的问题后来解决了,之前用vasp5算的坐标混乱,后来换了一台机器用vasp5.2计算就再没出现坐标错位的问题,可能是程序的问题?具体的我也不太清楚,你看你要不也换一台机器试一下?
10楼2013-12-01 16:05:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 pelargonium 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] filecode +8 documentary 2026-08-10 8/400 2026-08-10 19:20 by loufangrui
[基金申请] 帮忙看看fileCode +3 wwncly 2026-08-10 5/250 2026-08-10 19:06 by 2000zf36392
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +10 yufeiwaner 2026-08-09 11/550 2026-08-10 18:47 by yufeiwaner
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
信息提示
请填处理意见