版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

duba1023

金虫 (小有名气)

应助: 1 (幼儿园)
金币: 939.7
帖子: 139
在线: 35.6小时
虫号: 1732152
注册: 2012-04-02
专业: 凝聚态物性 II ：电子结构

[求助] vasp中计算磁性的相关问题已有1人参与

新人求助vasp的问题，在vasp里加磁矩怎么优化呢？我在INCAR里加入了
ISPIN= 2
MAGMOM= 12*2.08 4*0
优化完以后，OSZICAR中的mag值变为0了。
后来我看了些文献，尝试着加入了
GGA = 91
VOSKOWN = 1
mag=22左右
想问一下
1、INCAR里计算磁矩这样写可以吗？还需要加些什么吗？
2、OSZICAR中的mag值需要再除以原子数吗？还是这个就是每个原子的磁矩了？如果要除以每个原子的话，以本题为例，是应该除以12还是16呢？
求高手赐教，不胜感激。

我的INCAR设置：
SYSTEM = Fe-B
PREC = Accurate
ENCUT = 600.0
EDIFF = 1e-5
IBRION = 2
ISTART=0
ICHRAG=2
ISIF = 3
NSW = 100
GGA = 91
VOSKOWN = 1
ISMEAR =1 ; SIGMA = 0.20
POTIM = 0.100
LCHARG = FALSE
LWAVE = FALSE
PSTRESS = 0.0001
EDIFFG = 1e-4
ISPIN= 2
MAGMOM= 12*2.08 4*0

[ Last edited by duba1023 on 2013-12-30 at 14:24 ]

回复此楼

» 猜你喜欢

博士申请都是内定的吗？已经有7人回复
读博已经有5人回复
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？已经有5人回复
博士读完未来一定会好吗已经有29人回复
投稿精细化工已经有4人回复
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！已经有4人回复
导师想让我从独立一作变成了共一第一已经有9人回复
心脉受损已经有5人回复
Springer期刊投稿求助已经有4人回复
小论文投稿已经有3人回复

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

VASP计算中total-force为什么等于零已经有7人回复
VASP 中NELECT 的设置已经有8人回复
关于计算磁性交换常数J的问题已经有20人回复
请问vasp里面有Kubo-Greenwood公式的计算么已经有7人回复
VASP计算磁矩问题已经有8人回复
vasp计算磁性时，MAGMOM的设置已经有9人回复
求CASTEP 和 VASP 计算磁性问题已经有4人回复
问题：vasp计算突然中断，cpu利用率显示为0，这是什么情况，没有断电已经有3人回复
vasp磁性计算问题已经有8人回复
vasp计算磁性的问题已经有5人回复
大家交流一下用VASP计算非共线磁性的方法和结果分析已经有12人回复
VASP 计算功函数已经有10人回复
vasp 那些参数设置可以影响到磁性计算的磁矩已经有12人回复
vasp计算中，如何选择正确的K点已经有7人回复
非共线磁性计算已经有10人回复
vasp计算在进行静态计算时，出错。已经有3人回复
vasp计算磁矩的困惑已经有10人回复
VASP中计算dos的相关设置已经有6人回复
【求助】如何用VASP 计算体系的表面能已经有11人回复
【求助】vasp磁矩计算已经有8人回复

走七步，退七步，三思而后行。

1楼2013-12-30 13:52:21

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

xj20120828

木虫 (正式写手)

应助: 3 (幼儿园)
金币: 2320.6
散金: 476
红花: 1
帖子: 511
在线: 195.4小时
虫号: 2205360
注册: 2012-12-25
性别: MM
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
franch: 金币+2, 谢谢回帖交流 2013-12-30 22:30:37
duba1023: 金币+5, ★★★★★最佳答案, ^_^ 谢谢了 2013-12-31 14:42:36

查看磁矩的大小应该在OUTCAR中最后有一个magnetization (x)，命令为tail -300 OUTCAR ，可以看到磁矩的大小，总磁矩是0也有可能是反铁磁耦合。在INCAR中设置ISPIN和MAGMOM应该就可以了

赞一下

回复此楼

2楼2013-12-30 14:38:36

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

duba1023

金虫 (小有名气)

应助: 1 (幼儿园)
金币: 939.7
帖子: 139
在线: 35.6小时
虫号: 1732152
注册: 2012-04-02
专业: 凝聚态物性 II ：电子结构

引用回帖:

2楼: Originally posted by xj20120828 at 2013-12-30 14:38:36
查看磁矩的大小应该在OUTCAR中最后有一个magnetization (x)，命令为tail -300 OUTCAR ，可以看到磁矩的大小，总磁矩是0也有可能是反铁磁耦合。在INCAR中设置ISPIN和MAGMOM应该就可以了

我查了一下
eigenvalue-minimisations  :  6144
total energy-change (2. order) : 0.1423758E-05  (-0.7814050E-06)
number of electron  108.0000007 magnetization 23.1207282
augmentation part 50.1475618 magnetization 19.2734856
那我的体系里有12Fe和4B，我想知道Fe的磁矩，这个值是算错了吗？还是需要除以原子数的？

赞一下

回复此楼

走七步，退七步，三思而后行。

3楼2013-12-30 14:55:48

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

xj20120828

木虫 (正式写手)

应助: 3 (幼儿园)
金币: 2320.6
散金: 476
红花: 1
帖子: 511
在线: 195.4小时
虫号: 2205360
注册: 2012-12-25
性别: MM
专业: 凝聚态物性 II ：电子结构

引用回帖:

3楼: Originally posted by duba1023 at 2013-12-30 14:55:48
我查了一下
eigenvalue-minimisations  :  6144
total energy-change (2. order) : 0.1423758E-05  (-0.7814050E-06)
number of electron  108.0000007 magnetization 23.1207282
augmentation part 5 ...

在OUTCAR那里可以看到每个原子的磁矩

赞一下

回复此楼

4楼2013-12-30 16:00:53

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

duba1023

金虫 (小有名气)

应助: 1 (幼儿园)
金币: 939.7
帖子: 139
在线: 35.6小时
虫号: 1732152
注册: 2012-04-02
专业: 凝聚态物性 II ：电子结构

引用回帖:

4楼: Originally posted by xj20120828 at 2013-12-30 16:00:53
在OUTCAR那里可以看到每个原子的磁矩...

可是我确实没有看到每个原子的磁矩的显示，这是我前面用你说的命令提取出来的结果，里面没有看到每个原子的磁矩的信息。
5    0.6712    1.00000
   6    0.7358    1.00000
   7    1.0734    1.00000
   8    1.1229    1.00000
   9    2.1683    1.00000
   10    2.2124    1.00000
   11    2.7201    1.00000
   12    2.7549    1.00000
   13    3.4031    1.00000
   14    3.4939    1.00000
   15    3.5761    1.00000
   16    3.6723    1.00000
   17    4.1695    1.00000
   18    4.1800    1.00000
   19    4.2754    1.00000
   20    4.3058    1.00000
   21    4.5121    1.00000
   22    4.5637    1.00000
   23    4.6608    1.00000
   24    4.6939    1.00000
   25    4.7913    1.00000
   26    4.8451    1.00000
   27    4.9410    1.00000
   28    4.9616    1.00000
   29    5.1615    1.00000
   30    5.2005    1.00000
   31    5.2550    1.00000
   32    5.2996    1.00000
   33    5.4237    1.00000
   34    5.4628    1.00000
   35    5.5091    1.00000
   36    5.5867    1.00000
   37    6.1108    1.00014
   38    6.1318    1.00025
   39    6.3804    1.02410
   40    6.4176    1.03208
   41    6.4986    1.02371
   42    6.5127    1.01409
   43    6.6473    0.69919
   44    6.6638    0.63386
   45    7.3085    -0.00007
   46    7.3519    -0.00002
   47    7.4324    0.00000
   48    7.4722    0.00000
   49    7.5124    0.00000
   50    7.5327    0.00000
   51    7.5569    0.00000
   52    7.6152    0.00000
   53    7.7780    0.00000
   54    7.8257    0.00000
   55    7.8720    0.00000
   56    7.8843    0.00000
   57    7.9430    0.00000
   58    8.0020    0.00000
   59    8.0563    0.00000
   60    8.0952    0.00000
   61    8.1593    0.00000
   62    8.2137    0.00000
   63    8.2343    0.00000
   64    8.2831    0.00000
   65    8.4723    0.00000
   66    8.5259    0.00000
   67    8.5423    0.00000
   68    8.5668    0.00000
   69    9.4280    0.00000
   70    9.5185    0.00000
   71    9.7887    0.00000
   72    9.8605    0.00000
   73    10.8057    0.00000
   74    10.9211    0.00000
   75    11.4078    0.00000
   76    11.5144    0.00000
   77    12.7080    0.00000
   78    12.7637    0.00000
   79    13.2055    0.00000
   80    13.3343    0.00000

--------------------------------------------------------------------------------------------------------

soft charge-density along one line, spin component          1
      0       1       2       3       4       5       6       7       8       9
total charge-density along one line

soft charge-density along one line, spin component          2
      0       1       2       3       4       5       6       7       8       9
total charge-density along one line

pseudopotential strength for first ion, spin component:          1
-5.283 0.005 0.000 0.000 0.000  -0.726  -0.001 0.000
  0.005  -5.290 0.000 0.000 0.000  -0.001  -0.724 0.000
  0.000 0.000  -5.371  -0.011 0.019 0.000 0.000  -0.700
  0.000 0.000  -0.011  -5.371 0.040 0.000 0.000 0.003
  0.000 0.000 0.019 0.040  -5.409 0.000 0.000  -0.006
-0.726  -0.001 0.000 0.000 0.000 0.951 0.000 0.000
-0.001  -0.724 0.000 0.000 0.000 0.000 0.950 0.000
  0.000 0.000  -0.700 0.003  -0.006 0.000 0.000 0.944
  0.000 0.000 0.003  -0.700  -0.011 0.000 0.000  -0.001
  0.000 0.000  -0.006  -0.011  -0.689 0.000 0.000 0.002
  0.000 0.000 0.000  -0.004 0.008 0.000 0.000 0.000
  0.000 0.000 0.000 0.001  -0.001 0.000 0.000 0.000
  0.003  -0.003 0.000 0.000 0.000  -0.001 0.001 0.000
  0.000 0.000  -0.005 0.001 0.003 0.000 0.000 0.001
  0.000 0.000 0.001 0.002  -0.004 0.000 0.000 0.000
-0.001 0.001 0.000 0.000 0.000 0.000 0.000 0.000
  0.000 0.000 0.003  -0.001  -0.002 0.000 0.000  -0.001
  0.000 0.000  -0.002  -0.003 0.003 0.000 0.000 0.000
pseudopotential strength for first ion, spin component:          2
-3.796 0.043 0.000 0.000 0.000  -1.176  -0.012 0.000
  0.043  -3.774 0.000 0.000 0.000  -0.012  -1.181 0.000
  0.000 0.000  -3.772 0.044 0.044 0.000 0.000  -1.181
  0.000 0.000 0.044  -3.815 0.011 0.000 0.000  -0.012
  0.000 0.000 0.044 0.011  -3.674 0.000 0.000  -0.013
-1.176  -0.012 0.000 0.000 0.000 1.086 0.003 0.000
-0.012  -1.181 0.000 0.000 0.000 0.003 1.087 0.000
  0.000 0.000  -1.181  -0.012  -0.013 0.000 0.000 1.087
  0.000 0.000  -0.012  -1.169  -0.003 0.000 0.000 0.003
  0.000 0.000  -0.013  -0.003  -1.209 0.000 0.000 0.004
  0.000 0.000 0.010  -0.012 0.010 0.000 0.000  -0.002
  0.000 0.000  -0.001 0.001  -0.001 0.000 0.000 0.000
-0.003  -0.002 0.000 0.000 0.000 0.001 0.000 0.000
  0.000 0.000  -0.005 0.000 0.001 0.000 0.000 0.001
  0.000 0.000 0.004 0.000  -0.004 0.000 0.000  -0.001
  0.000 0.001 0.000 0.000 0.000 0.000 0.000 0.000
  0.000 0.000 0.003  -0.001  -0.001 0.000 0.000  -0.001
  0.000 0.000  -0.002  -0.002 0.003 0.000 0.000 0.001
total augmentation occupancy for first ion, spin component:          1
  1.469  -0.024 0.004 0.000  -0.002 0.102  -0.035  -0.001 0.003 0.004  -0.002 0.000 0.015 0.003 0.000  -0.001
-0.024 1.425 0.003 0.002 0.003  -0.035 0.091 0.000  -0.003  -0.003 0.002 0.000  -0.013 0.001 0.000 0.001
  0.004 0.003 1.364  -0.033  -0.001  -0.001 0.000 0.129  -0.053  -0.026  -0.038  -0.013 0.000  -0.038 0.020 0.000
  0.000 0.002  -0.033 1.373 0.013 0.003  -0.003  -0.053 0.151  -0.036  -0.057 0.027 0.001 0.012 0.037 0.000
-0.002 0.003  -0.001 0.013 1.282 0.004  -0.002  -0.026  -0.036 0.045  -0.011  -0.013  -0.001 0.017  -0.035 0.000
  0.102  -0.035  -0.001 0.003 0.004 0.173  -0.010 0.000  -0.003  -0.003  -0.002 0.000 0.001  -0.001 0.000 0.001
-0.035 0.091 0.000  -0.003  -0.002  -0.010 0.174 0.001 0.002 0.002 0.000 0.000  -0.011 0.000 0.000 0.001
-0.001 0.000 0.129  -0.053  -0.026 0.000 0.001 0.236  -0.001  -0.008 0.053  -0.037 0.000  -0.032 0.026 0.000
  0.003  -0.003  -0.053 0.151  -0.036  -0.003 0.002  -0.001 0.226  -0.038  -0.111 0.046 0.001 0.011 0.026 0.000
  0.004  -0.003  -0.026  -0.036 0.045  -0.003 0.002  -0.008  -0.038 0.204 0.059  -0.037 0.000 0.024  -0.034 0.000
-0.002 0.002  -0.038  -0.057  -0.011  -0.002 0.000 0.053  -0.111 0.059 1.243  -0.371 0.000  -0.014 0.000 0.000
  0.000 0.000  -0.013 0.027  -0.013 0.000 0.000  -0.037 0.046  -0.037  -0.371 0.147 0.000 0.005 0.014 0.000
  0.015  -0.013 0.000 0.001  -0.001 0.001  -0.011 0.000 0.001 0.000 0.000 0.000 0.209 0.000 0.002  -0.014
  0.003 0.001  -0.038 0.012 0.017  -0.001 0.000  -0.032 0.011 0.024  -0.014 0.005 0.000 0.238  -0.030 0.000
  0.000 0.000 0.020 0.037  -0.035 0.000 0.000 0.026 0.026  -0.034 0.000 0.014 0.002  -0.030 0.229 0.000
-0.001 0.001 0.000 0.000 0.000 0.001 0.001 0.000 0.000 0.000 0.000 0.000  -0.014 0.000 0.000 0.001
  0.000 0.000 0.002 0.000  -0.001 0.000 0.000 0.003  -0.002  -0.002 0.002  -0.001 0.000  -0.020 0.003 0.000
  0.000 0.000  -0.001  -0.001 0.002 0.000 0.000  -0.004  -0.003 0.004 0.002  -0.002 0.000 0.003  -0.018 0.000
total augmentation occupancy for first ion, spin component:          2
  0.344 0.051  -0.001  -0.008  -0.004  -0.240  -0.012 0.000 0.004 0.003 0.002 0.000  -0.020 0.000 0.000 0.000
  0.051 0.350 0.000 0.002 0.005  -0.012  -0.226  -0.002  -0.001  -0.003  -0.004 0.000 0.001  -0.001 0.001 0.000
-0.001 0.000 0.441 0.050 0.024 0.000  -0.002  -0.236  -0.013  -0.006 0.014 0.001 0.001  -0.001 0.008 0.000
-0.008 0.002 0.050 0.401  -0.022 0.004  -0.001  -0.013  -0.227 0.002 0.007  -0.004  -0.002 0.003 0.001 0.000
-0.004 0.005 0.024  -0.022 0.570 0.003  -0.003  -0.005 0.004  -0.259 0.009 0.010 0.002  -0.008 0.003 0.000
-0.240  -0.012 0.000 0.004 0.003 0.029 0.005  -0.001  -0.003  -0.003  -0.001 0.000 0.005  -0.001 0.000 0.000
-0.012  -0.226  -0.002  -0.001  -0.003 0.005 0.026 0.001 0.002 0.002 0.001 0.000 0.001 0.000 0.000 0.000
  0.000  -0.002  -0.236  -0.013  -0.005  -0.001 0.001 0.012 0.015 0.007 0.007 0.003 0.000 0.009  -0.004 0.000
  0.004  -0.001  -0.013  -0.227 0.004  -0.003 0.002 0.015 0.005 0.011 0.000  -0.009 0.000 0.001  -0.003 0.000
  0.003  -0.003  -0.006 0.002  -0.259  -0.003 0.002 0.007 0.011 0.027 0.001 0.003 0.000 0.004 0.004 0.000
  0.002  -0.004 0.014 0.007 0.009  -0.001 0.001 0.007 0.000 0.001 0.026 0.083 0.001  -0.001  -0.006 0.000
  0.000 0.000 0.001  -0.004 0.010 0.000 0.000 0.003  -0.009 0.003 0.083  -0.051 0.000  -0.002  -0.001 0.000
-0.020 0.001 0.001  -0.002 0.002 0.005 0.001 0.000 0.000 0.000 0.001 0.000  -0.010 0.000 0.001 0.002
  0.000  -0.001  -0.001 0.003  -0.008  -0.001 0.000 0.009 0.001 0.004  -0.001  -0.002 0.000  -0.003  -0.001 0.000
  0.000 0.001 0.008 0.001 0.003 0.000 0.000  -0.004  -0.003 0.004  -0.006  -0.001 0.001  -0.001  -0.009 0.000
  0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.000 0.000 0.000
  0.000 0.000  -0.001 0.000 0.001 0.000 0.000  -0.001 0.000 0.000 0.000 0.000 0.000 0.003 0.000 0.000
  0.000 0.000 0.001 0.001  -0.002 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.003 0.000

------------------------ aborting loop because EDIFF is reached ----------------------------------------

--------------------------------------------------------------------------------------------------------

CHARGE:  VPU time 2.54: CPU time 2.55
FORLOC:  VPU time 0.01: CPU time 0.01
FORNL :  VPU time 14.97: CPU time 15.03
STRESS:  VPU time 15.40: CPU time 15.48
FORCOR:  VPU time 1.02: CPU time 1.02
FORHAR:  VPU time 0.04: CPU time 0.04
MIXING:  VPU time 0.00: CPU time 0.00
d Force = 0.3013047E-04[-0.508E-05, 0.653E-04]  d Energy = 0.2434815E-04 0.578E-05
d Force =-0.2739393E-01[-0.280E-01,-0.268E-01]  d Ewald  =-0.3150170E+00 0.288E+00

--------------------------------------------------------------------------------------------------------

  FREE ENERGIE OF THE ION-ELECTRON SYSTEM (eV)
  ---------------------------------------------------
  free  energy TOTEN  =    -156.109246 eV

  energy  without entropy=    -156.112698  energy(sigma->0) =    -156.110397
  enthalpy is  TOTEN =    -156.109237 eV P V=       0.000010

  FORCE on cell =-STRESS in cart. coord.  units (eV/reduce length):
  Direction X       Y       Z       XY    YZ    ZX
  --------------------------------------------------------------------------------------
  Alpha Z 757.16 757.16 757.16
  Ewald -2555.28  -2487.04  -2494.29    0.00    0.00    0.00
  Hartree 370.78 406.57 416.60    0.00    0.00    0.00
  E(xc)    -518.47 -518.67 -518.74    0.00    0.00    0.00
  Local    774.67 667.39 662.90    0.00    0.00    0.00
  n-local -309.12 -308.80 -308.93    -0.30    -0.25    -0.21
  augment 463.47 465.35 465.55    0.00    0.00    0.00
  Kinetic 1016.81 1017.99 1019.72    1.10    0.05    -0.81
  -------------------------------------------------------------------------------------
  Total       0.02    -0.04    -0.02    0.00    0.00    0.00
  in kB       0.23    -0.39    -0.26    0.00    0.00    0.00
  external pressure =    -0.14 kB  Pullay stress =       0.00 kB

VOLUME and BASIS-vectors are now :
-----------------------------------------------------------------------------
  energy-cutoff  :    600.00
  volume of cell :    154.41
   direct lattice vectors                reciprocal lattice vectors
   5.343016377  0.000000000  0.000000000    0.187160197  0.000000000  0.000000000
   0.000000000  6.606270236  0.000000000    0.000000000  0.151371343  0.000000000
   0.000000000  0.000000000  4.374636295    0.000000000  0.000000000  0.228590432

  length of vectors
   5.343016377  6.606270236  4.374636295    0.187160197  0.151371343  0.228590432

FORCES acting on ions
electron-ion (+dipol)          ewald-force                   non-local-force                convergence-correction
-----------------------------------------------------------------------------------------------
-.117E+01 0.137E-11 0.181E+01 0.150E+01 0.324E-13 -.203E+01 -.317E+00 0.000E+00 0.226E+00 -.444E-02 -.112E-11 0.177E-02
-.117E+01 -.161E-12 -.181E+01 0.150E+01 -.170E-14 0.203E+01 -.317E+00 0.000E+00 -.226E+00 -.444E-02 -.836E-12 -.177E-02
0.117E+01 0.605E-13 -.181E+01 -.150E+01 -.314E-13 0.203E+01 0.317E+00 0.000E+00 -.226E+00 0.444E-02 0.114E-11 -.177E-02
-.926E+01 -.164E+01 -.229E+01 0.965E+01 0.205E+01 0.211E+01 -.402E+00 -.412E+00 0.179E+00 -.919E-03 0.163E-02 0.179E-02
-.926E+01 0.164E+01 0.229E+01 0.965E+01 -.205E+01 -.211E+01 -.402E+00 0.412E+00 -.179E+00 -.919E-03 -.163E-02 -.179E-02
-.926E+01 0.164E+01 -.229E+01 0.965E+01 -.205E+01 0.211E+01 -.402E+00 0.412E+00 0.179E+00 -.919E-03 -.163E-02 0.179E-02
0.926E+01 0.164E+01 -.229E+01 -.965E+01 -.205E+01 0.211E+01 0.402E+00 0.412E+00 0.179E+00 0.919E-03 -.163E-02 0.179E-02
0.926E+01 0.164E+01 0.229E+01 -.965E+01 -.205E+01 -.211E+01 0.402E+00 0.412E+00 -.179E+00 0.919E-03 -.163E-02 -.179E-02
0.926E+01 -.164E+01 -.229E+01 -.965E+01 0.205E+01 0.211E+01 0.402E+00 -.412E+00 0.179E+00 0.919E-03 0.163E-02 0.179E-02
-.926E+01 -.164E+01 0.229E+01 0.965E+01 0.205E+01 -.211E+01 -.402E+00 -.412E+00 -.179E+00 -.919E-03 0.163E-02 -.179E-02
0.926E+01 -.164E+01 0.229E+01 -.965E+01 0.205E+01 -.211E+01 0.402E+00 -.412E+00 -.179E+00 0.919E-03 0.163E-02 -.179E-02
0.117E+01 0.545E-12 0.181E+01 -.150E+01 0.287E-13 -.203E+01 0.317E+00 0.000E+00 0.226E+00 0.444E-02 0.765E-12 0.177E-02
-.381E+00 -.650E-12 -.153E+01 0.461E+00 -.116E-14 0.165E+01 -.732E-01 0.867E-18 -.107E+00 -.663E-03 0.421E-12 0.172E-02
0.381E+00 -.133E-12 -.153E+01 -.461E+00 0.152E-13 0.165E+01 0.732E-01 0.867E-18 -.107E+00 0.663E-03 -.271E-12 0.172E-02
-.381E+00 -.683E-13 0.153E+01 0.461E+00 -.795E-15 -.165E+01 -.732E-01 -.867E-18 0.107E+00 -.663E-03 0.294E-12 -.172E-02
0.381E+00 0.269E-14 0.153E+01 -.461E+00 0.220E-13 -.165E+01 0.732E-01 -.867E-18 0.107E+00 0.663E-03 -.301E-12 -.172E-02
-----------------------------------------------------------------------------------------------
0.153E-11 0.752E-11 -.148E-10 -.810E-14 0.710E-14 0.311E-14 0.694E-16 0.167E-15 0.000E+00 0.180E-12 -.750E-12 0.288E-12

POSITION                                     TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
   2.55934    4.95470    1.65105       0.005587    0.000000    0.001976
   5.23085    4.95470    0.53627       0.005587    0.000000    -0.001976
   2.78368    1.65157    2.72359       -0.005587    0.000000    -0.001976
   3.61296    2.93668    0.65825       -0.008632    -0.002171    -0.001395
   0.94145    0.36645    1.52907       -0.008632    0.002171    0.001395
   3.61296    0.36645    0.65825       -0.008632    0.002171    -0.001395
   4.40157    3.66959    2.84556       0.008632    0.002171    -0.001395
   1.73006    3.66959    3.71639       0.008632    0.002171    0.001395
   4.40157    6.23982    2.84556       0.008632    -0.002171    -0.001395
   0.94145    2.93668    1.52907       -0.008632    -0.002171    0.001395
   1.73006    6.23982    3.71639       0.008632    -0.002171    0.001395
   0.11217    1.65157    3.83837       -0.005587    0.000000    0.001976
   4.71540    1.65157    1.86191       0.005763    0.000000    0.005001
   3.29913    4.95470    4.04923       -0.005763    0.000000    0.005001
   2.04389    1.65157    0.32541       0.005763    0.000000    -0.005001
   0.62762    4.95470    2.51272       -0.005763    0.000000    -0.005001
-----------------------------------------------------------------------------------
total drift:                               0.000000    0.000000    0.000000

--------------------------------------------------------------------------------------------------------

--------------------------------------------------------------------------------------------------------


reached required accuracy - stopping structural energy minimisation
   LOOP+:  VPU time  236.26: CPU time  237.49


General timing and accounting informations for this job:
========================================================

               Total CPU time used (sec):    8389.433
                        User time (sec):    8352.846
                        System time (sec):    36.586
                     Elapsed time (sec):    8434.046

               Maximum memory used (kb):    137532.
               Average memory used (kb):          0.

                        Minor page faults:       37403
                        Major page faults:          0
               Voluntary context switches:          27

赞一下

回复此楼

走七步，退七步，三思而后行。

5楼2013-12-30 16:06:15

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

xj20120828

木虫 (正式写手)

应助: 3 (幼儿园)
金币: 2320.6
散金: 476
红花: 1
帖子: 511
在线: 195.4小时
虫号: 2205360
注册: 2012-12-25
性别: MM
专业: 凝聚态物性 II ：电子结构

引用回帖:

5楼: Originally posted by duba1023 at 2013-12-30 16:06:15
可是我确实没有看到每个原子的磁矩的显示，这是我前面用你说的命令提取出来的结果，里面没有看到每个原子的磁矩的信息。
5    0.6712    1.00000
   6    0.7358    1.00000
   7    1 ...

magnetization (x)

# of ion    s    p    d    tot
----------------------------------------
  1       0.000 0.000 0.000 0.000
  2       0.001 0.000 0.000 0.001
  3       0.001 0.000 0.000 0.001
  4       0.000 0.000 0.000 0.000
  5       0.002 0.000 0.000 0.002
  6       0.003 0.000 0.000 0.003
  7    -0.010 0.000 0.000  -0.010
  8    -0.003 0.000 0.000  -0.003
  9       0.010 0.000 0.000 0.011
10       0.004 0.000 0.000 0.004
11       0.001 0.001 0.000 0.003
12       0.002 0.002 0.000 0.004
13       0.001  -0.002 0.000  -0.001
14       0.001 0.000 0.000 0.002
15    -0.006  -0.014 0.000  -0.020
16       0.002  -0.002 0.000 0.000
17       0.002 0.011 0.000 0.013
18       0.003 0.014 0.000 0.016
19    -0.008  -0.022 0.000  -0.031
20    -0.002  -0.002 0.000  -0.003
21       0.001 0.007 0.309 0.318
22       0.038 0.028 3.419 3.485
------------------------------------------------
tot       0.044 0.022 3.728 3.794

赞一下

回复此楼

6楼2013-12-30 16:51:04

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

duba1023

金虫 (小有名气)

应助: 1 (幼儿园)
金币: 939.7
帖子: 139
在线: 35.6小时
虫号: 1732152
注册: 2012-04-02
专业: 凝聚态物性 II ：电子结构

引用回帖:

6楼: Originally posted by xj20120828 at 2013-12-30 16:51:04
magnetization (x)

# of ion    s    p    d    tot
----------------------------------------
  1       0.000 0.000 0.000 0.000
  2       0.001 0.000 0.000 0.001
  3    ...

那我能再问一下吗？就对于每个Fe而言，从这里面看每个Fe的磁矩是？ tot的值指的是？

赞一下

回复此楼

走七步，退七步，三思而后行。

7楼2013-12-31 08:45:47

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