24小时热门版块排行榜    

Znn3bq.jpeg
查看: 8175  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 灯塔守望者 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 化学专业申博 +3 赵子羊 2026-05-23 4/200 2026-05-24 18:10 by 工大学长
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 1rx34o113h 2026-05-23 3/150 2026-05-24 17:41 by 0i3mu4vkjz
[基金申请] 评审有感 +16 popular289 2026-05-18 27/1350 2026-05-24 17:34 by hhs666
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 hvkbtfonbv 2026-05-23 4/200 2026-05-24 17:21 by 75ui6h7z2t
[博后之家] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 hvkbtfonbv 2026-05-23 3/150 2026-05-24 17:10 by 75ui6h7z2t
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 a2tycdlnq1 2026-05-23 5/250 2026-05-24 16:21 by hhx1yx9evi
[论文投稿] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 a2tycdlnq1 2026-05-23 4/200 2026-05-24 16:16 by hhx1yx9evi
[基金申请] 河北省自然科学基金 +6 Peterchao 2026-05-18 9/450 2026-05-24 16:02 by 130067131
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 pmo95bazuy 2026-05-23 8/400 2026-05-24 15:56 by 1uy1ht2y9r
[基金申请] 西安交大新媒学院副院长用撤稿论文结题 +3 bjvtcliu 2026-05-24 5/250 2026-05-24 10:16 by kudofaye
[教师之家] 某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路! +4 zju2000 2026-05-21 5/250 2026-05-24 09:35 by songwz
[考博] 26/27申博自荐 10+4 ZXW0202 2026-05-22 9/450 2026-05-24 08:47 by bjvtcliu
[考博] 博士申请 +3 焦晓明 2026-05-21 3/150 2026-05-23 11:26 by mlc840311
[论文投稿] 投稿求助,期刊 +4 希冀,有书读 2026-05-20 8/400 2026-05-22 10:16 by 希冀,有书读
[文学芳草园] 献血感触 +7 呀呀好傻 2026-05-19 13/650 2026-05-21 20:15 by 呀呀好傻
[基金申请] 国自然评分 +4 无名者登山 2026-05-20 5/250 2026-05-21 16:35 by swuq
[基金申请] 国自然上会要求 +7 无名者登山 2026-05-18 11/550 2026-05-21 15:50 by draco1987
[基金申请] 提交了我也来说说感想 +9 fummck 2026-05-20 10/500 2026-05-21 14:17 by draco1987
[有机交流] 反应很差,大量原料没有反应 5+3 Mr.Zot 2026-05-19 8/400 2026-05-20 22:19 by Equinoxhua
[考博] 如果工作了想读博,可以边工作边读全日制嘛? 30+3 铁达火车 2026-05-18 5/250 2026-05-20 09:33 by tfang
信息提示
请填处理意见