【其他】征集本版成功求助贴 - 第一性原理 - 其他 - 第6页小木虫论坛-学术科研互动平台

51楼2009-06-04 22:07:26

mozhui

应助: (幼儿园)
在线:
虫号: 596549

★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
wuli8(金币+5,VIP+0):感谢，参与 6-5 20:51

http://muchong.com/bbs/viewthread.php?tid=1341516&fpage=2

Vasp&MedeA

【求助】DOS相加的fortran代码？

split_dos，大家常用。不知大家有没有编写过DOS相加的代码，能否上传？

表面上的原子较多，平均DOS靠ORIGIN太烦锁，自己大致能看fortran，写基本不会，望经验者指点，谢谢

3楼：有个sum—dos的，上网能下到
4楼：谢谢LS提醒。tools里面就有，没注意到

52楼2009-06-05 12:48:38

mozhui

应助: (幼儿园)
在线:
虫号: 596549

★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
wuli8(金币+5,VIP+0):感谢，参与 6-5 21:13

http://muchong.com/bbs/viewthread.php?tid=1355620&fpage=1
其他
【求助】XCRYSDEN怎么调背景颜色

即background颜色，默认的背景颜色是黑的，
我想调成白色的，找了半天没找到，晶胞，平面，原子等的颜色倒是可以调
请高手指教怎么调整背景颜色。

1楼：哈哈哈.左上角有个按钮可以调的，换各种颜色背景。

53楼2009-06-05 20:59:12

mozhui

应助: (幼儿园)
在线:
虫号: 596549

★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
wuli8(金币+5,VIP+0):感谢，参与 6-5 21:14

http://muchong.com/bbs/viewthread.php?tid=1364716&fpage=1

Vasp&MedeA

【讨论】VASP中POSCAR文件

VASP中的POSCAR文件具体应该怎么写呢？初学者请教大家
是不是不同的晶体结构写的不一样呢？
第二行写成1.00可以吗？

3楼：
of course. you can write 1.0 in the second line.
you could write lattice vector with ang in the follow three lines.
it is very easy to do it.
ex:
POSCAR
case
1.0
ax ay az
bx by bz
cx cy cz
the number of first atom "space" the number of second atom ...
Direct
coordinate of every atom
over
diamond
1.000000000000000
   0.0000000000000000 1.7626797900540120 1.7626797900540120
   1.7626797900540120 0.0000000000000000 1.7626797900540120
   1.7626797900540120 1.7626797900540120 0.0000000000000000
2
Direct
  0.0000000000000000  0.0000000000000000  0.0000000000000000
  0.2500000000000000  0.2500000000000000  0.2500000000000000

  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
primitive cell

4楼：帮助链接：http://muchong.com/bbs/viewthread.php?tid=646969

54楼2009-06-05 21:03:57

mozhui

应助: (幼儿园)
在线:
虫号: 596549

★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
wuli8(金币+5,VIP+0):感谢，参与。希望您继续参与，谢谢！ 6-5 21:14

http://muchong.com/bbs/viewthread.php?tid=1343599&fpage=2

Vasp&MedeA

【求助】vasp如何后台运行

各位朋友，请问vasp如何后台运行呢？比如gaussian，后加个&就可以了，这个vasp呢？我不太懂了，我总不能总开着电脑啊~~呵呵，谢谢帮我这个小忙~~

2楼：一样的，&
3楼：nohup vasp &
这样退出后还能继续算
4楼：其实还可以用shell的script控制，以下是一个例子：
c shell环境，使用MPI并行，8台机器的计算。

如果这个文件名为
run-job.script，在运行是直接 nohup ./run-job.script &即可

#!/bin/csh
#设定vASP路径,此例是说明VASP在/research/bin/vasp
setenv VASP /research/bin/vasp
#
#设定当前工作的路径
setenv WORK `pwd`
#
#
#使用MPI进行并行计算
cd $WORK
wipe -v HOSTS
lamboot -v HOSTS
echo -n "start time  " > time ; date >> time
mpirun -c 8 n0-7 -ssi rpi tcp -O -w  $VASP >& vasp.out
echo -n "end time  " >> time ; date >> time
cat time >> vasp.out
#
wipe -v HOSTS
#
#
exit

楼主在帖子里的其他问题：请问你说的这个事什么意思呢？输入nohup vasp &后，是后台运行还是中途退出啊？你说的继续运算是怎么回事？谢谢！nohup vasp &
这个还是在前台算啊
不过只是可以中途退出

10楼：在后面加个输出文件，如：
nohup vasp &out
你的输出内容就会记录在out里，而不会影响你其他的工作。这样就能随意进行别的任务操作了

55楼2009-06-05 21:11:36

mozhui

应助: (幼儿园)
在线:
虫号: 596549

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
wuli8(金币+5,VIP+0):感谢，参与 6-5 21:21
wuchenwf(金币+5,VIP+0):谢谢 6-5 21:22

http://muchong.com/bbs/viewthread.php?tid=1209018&fpage=8

Vasp&MedeA

【求助】关于居里温度。

我在用VASP计算材料的磁性。请问如何利用VASP得到体系的居里温度？计算居里温度的意义是什么？

谢谢！

2楼回答了居里点的定义：我对VASP一知半解，不会居里温度计算。
就回答你后一半吧

铁磁体和顺磁体之间改变的温度，
磁制冷知道吗，关键就是调整居里温度，
永磁体材料也很关注居里温度高低。

3楼：VASP计算实际上并不能对温度进行设置,属于零温计算.至于居里温度,反映了两种不同的相的转变温度,比如由顺磁相到铁磁相,可以查阅铁磁学之类的资料,这方面讲解较为详细.据我所知,利用VASP进行计算得到相关数值,居里温度是可以根据一些理论(如mean-field approximation)可进行计算估计居里温度的.这方面一些文献曾有过这样的报道.

赞一下(21人)

56楼2009-06-05 21:18:19

mozhui

应助: (幼儿园)
在线:
虫号: 596549

http://muchong.com/bbs/viewthread.php?tid=1192656&fpage=9

Vasp&MedeA

求助：功函数

用VASP怎么算功函数？

2楼：
http://new.quantumchemistry.net/ ... le.asp?ArticleID=88
自己看一下吧。

在VASP的INCAR文件中，指定LVTOT=.TRUE. ,并且在计算时超晶胞大小选取有足够的真空区域，以保证势能函数收敛到真空能级

。指定LVTOT=.TRUE.后，VASP的运行结果会多出一个叫LOCPOT的文件，里面包含的信息是计算中的electrostatic potential。

在真空区域的electrostatic potential即为真空能级。真空能级减去费米能级即为该体系的功函数。

3楼：左右真空能级减去费米能级，然后与实验数值对照。。。

赞一下

57楼2009-06-05 21:29:49

mozhui

应助: (幼儿园)
在线:
虫号: 596549

★ ★ ★
wuli8(金币+3,VIP+0):感谢参与，(*^__^*) …… 6-7 08:37

http://muchong.com/bbs/viewthread.php?tid=1168936&fpage=9

Vasp&MedeA

【求助】大家有没见过这个warning

计算开始不久报错：WARNING: random wavefunctions but no delay for mixing, default for NELMDL，然后刚开始写电子步就停止了，每次提交都是这样，请教下这时什么问题

我看这个值有正负之分，能否解释下有什么区别，才结束vasp，不怎么了解，呵呵

可能的原因
3楼：设置一下 NELMDL即可
5楼：增大NELMDL的值
7楼：NELMDL是指定每个离子步中非自洽电子步数
负值是在第一个离子步有效
正值是每个离子步都有效
8楼：这个不是很要紧，我的任务出现这个提示WARNING，照算不误，最后收敛。

你的计算中断估计有其他原因

楼主自己的结论：多谢楼上各位的回复，问题解决了，应该是空间不够造成的，这个WARNING应该没什么影响

赞一下(20人)

58楼2009-06-05 21:35:52

mozhui

应助: (幼儿园)
在线:
虫号: 596549

★ ★ ★
wuli8(金币+3,VIP+0):感谢参与，(*^__^*) …… 6-7 08:37

http://muchong.com/bbs/viewthread.php?tid=1058778&fpage=12

Vasp&MedeA

vasp计算态密度问题

想问一下高手vasp算态密度的时候如何让DOSCAR的数值密集一些，就是DOSCAR数值多一些！有什么参数吗，我增加K点好像没有任何效果。

2楼：将SIGMA的值设置小一点。
3楼：要使态密度结果精确，接近实际，也就是尖峰比较多，则将 sigma 设小即可。。。
要是DOSCAR中的数值密集（默认是301行），则将 NEDOS 设大即可。。。

赞一下(20人)

59楼2009-06-05 21:40:24

mozhui

应助: (幼儿园)
在线:
虫号: 596549

★ ★
wuli8(金币+2,VIP+0):感谢参与，(*^__^*) …… 6-7 08:38

http://muchong.com/bbs/viewthread.php?tid=1052373&fpage=13

Vasp&MedeA

【求助】关于ENCUT的设置问题

在算band时根据文献我用了ENCUT=400,的设置，可是做完静态计算后，发现runlog中出现下面的提示
One real space projector is optimized for                            |
|    E = 449.69, eV  but you are using a  cutoff of                |
|    ENMAX= 400.00 eV  ( QCUT=    5.42 a.u.)                         |
|    This makes no sense reoptimize the projector                      |
|    with the a.u. value given above

这个要怎么改呢，我需要把ENCUT改成449.69吗

2楼：ENCUT按说需要做收敛测试才能给出可靠值，如果不做测试，至少给出enmax的1.5倍比较保险
6楼：千万不要把encut搞那么高，，浪费机时啊，可以调整其他参数就没有那个建议了，，lreal
你试试看

赞一下(20人)

60楼2009-06-05 21:45:15