24小时热门版块排行榜    

查看: 7578  |  回复: 72
当前主题已经存档。

jghe

金虫 (著名写手)

请问楼主对优化六角和正交体系进行结构优化?
候柱峰那本书看过,有点看不懂。主要是怎么确定a/c的值。望指点!
31楼2009-02-26 20:24:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fansw1129

铁杆木虫 (著名写手)

利用FPLAPW(wien2k)在计算掺杂体系的磁性时候
关于k点的选取时候有项 ADD INVERSION 如何选取,选择yes还是no
不同的选项对计算的结果有多大的影响?
32楼2009-03-02 08:50:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vboyz

新虫 (初入文坛)

刚刚接触vasp,老板就让做磁性,一头雾水中,有这个帖子真好啊……
1.我现在在计算磁性时设置的参数磁性参数只有MAGMOM,在网上很少看到其他参数的介绍,不知道还有哪些比较重要,一般不用默认值的参数。
2.我现在在做一些有机分子材料,吸附Fe原子的磁性,我想知道对于两个反铁磁的Fe原子的MAGMOM应该怎样设置,采用1,-1可以吗,还是需要计算初始磁矩呢?如何计算?
体系中的C,H原子,考虑或不考虑铁磁,又该怎么处理呢?没有磁性是不是就设成0呢?
刚刚上手,好不容易装上了,实在是对这个软件了解甚少。期待高手指教!
33楼2009-03-03 16:35:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

灯塔守望者

铁杆木虫 (正式写手)

★ ★ ★
qasd(金币+3,VIP+0):xiexie~ 4-9 21:29
引用回帖:
Originally posted by jghe at 2009-2-26 20:24:
请问楼主对优化六角和正交体系进行结构优化?
候柱峰那本书看过,有点看不懂。主要是怎么确定a/c的值。望指点!

这需要优化两个参数,a, c/a
一般做的话,选择c/a为实验数值,然后优化a,之后再固定a,优化c/a,一般c/a值离实验值偏差很小。
34楼2009-03-04 10:16:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

灯塔守望者

铁杆木虫 (正式写手)

★ ★ ★
qasd(金币+3,VIP+0):xiexie~ 4-9 21:29
引用回帖:
Originally posted by vboyz at 2009-3-3 16:35:
刚刚接触vasp,老板就让做磁性,一头雾水中,有这个帖子真好啊……
1.我现在在计算磁性时设置的参数磁性参数只有MAGMOM,在网上很少看到其他参数的介绍,不知道还有哪些比较重要,一般不用默认值的参数。
2.我现 ...

1  如果你研究共线磁性的话,这两个数值必须得设定
   ISPIN=2
   MAGMOM =1*3 1*3 16*0 24*0
其中n*M, n为多少个原子的磁矩为M,注意必须按POTCAR中元素顺序来设定。
如果是反铁磁性的话,可以设为3,-3(比如)。如果是非磁性元素,设为0就可以了。

2,可以设为1,-1,不过最好取大点,原则上初始磁矩大小对结果没影响,就担心计算过程中不同的初始磁矩收敛到不同的亚稳态,一般这种情况很少,不过你取两次不同的初始磁矩,比较一下最后的收敛磁矩是否相同。

[ Last edited by 灯塔守望者 on 2009-3-4 at 10:24 ]
35楼2009-03-04 10:21:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vboyz

新虫 (初入文坛)

谢谢灯塔守望者,今天看到帖子后一直在算,但是ISPIN打开后,我的计算就好难收敛,问了一些人,把K点调了一下,encut也试着提高过,但是还是没什么改观。
请问磁性计算中影响收敛的有哪些设置?我现在80几步收敛正常吗?
36楼2009-03-04 21:41:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

灯塔守望者

铁杆木虫 (正式写手)

★ ★ ★
qasd(金币+3,VIP+0):xiexie 3-5 12:26
不收敛的可能性有很多种,你首先检查你的体系建的是否正确,然后改变不用的优化方法,再试试看。
以下是我的INCAR,仅供参考。
   ISTART =      0
   ISPIN=2
   MAGMOM =3 -3 16*0 24*0

Electronic Relaxation
   ISMEAR =   0
   POTIM  =0.5
   SIGMA  =  0.05
   NELMIN =     4
   NELM   =    30
   EDIFF  =   1E-05  stopping-criterion for ELM
Ionic Relaxation 1
   NSW    =   30
   IBRION =     2 (或1,适用于优化结构变化不大的情况)
   EDIFFG =   1E-03   stopping-criterion for ELM
   ISIF   =      2
DOS related values:
   RWIGS = 1.10 1.1 0.7

Global parameter

LREAL  =Auto
ALGO  = Fast
NSIM  =4
LPLANE = .TRUE.
NPAR = 1
# PREC   =     Accurate
  ENCUT=450
  LWAVE= F
37楼2009-03-05 10:30:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

miracle3130

铁杆木虫 (小有名气)

你好,想请教一个问题:
使用VASP计算非共线反铁磁问题,在设置MAGOM时,按照要求每个原子的x,y,z三个方向上都设置了初始磁矩.由于超胞原子数近50个,在设置MAGOM这一参数时,总字符数(包括空格)超过了一行所规定的最多255个字符,运行VASP,错误显示在读MAGOM时,数据只能读取总数据的2/3左右.
       试图寻找过VASP的续行问题解决办法,vasp.4.lib中的drdatab.F有这么几行话:
        !     -- data lists can be continued over several lines (with the only     *
        !        restriction that the total input may not be longer than 32767  *
        !        characters ...). Continuations must be marked by some \ at     *
        !        the end of the previous line (but really at the end of the        *
        !        line!! -- so lines to be continued may not contain comments!)  *
试着利用反斜杠"\"续行,一部分数据后面"空格"加"\"加回车、数据后面直接跟"\"加回车、或者是以上两种情况最后不加回车这些情况我都试过,可是问题始终解决不了,MAGOM的所有数据不能读完整.
       不知是否有人也碰到过这样的类似问题!希望能提供正确的续行格式,用简单的几个数据作为例子,不甚感激!
38楼2009-03-11 16:00:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ztwang

铜虫 (小有名气)

请问vasp在计算磁性的时候,oszicar中得到的磁矩和outcar中得到各原子磁矩之和不一致,那么总磁矩取oszicar的磁矩(跟前人的文献一致),但单个原子的磁矩如何取得更准确!谢谢
39楼2009-04-09 08:51:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chgzhong

木虫 (正式写手)

请问版主, 在Linux中如果一行输入的字符太多出错,要想换行继续,那必须用什么指令?
40楼2009-04-09 19:59:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 灯塔守望者 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见