| 查看: 15749 | 回复: 314 | |||
| 当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖 | |||
gavinliu7390木虫 (著名写手)
叶落鹰飞
|
[交流]
【其他】在线答疑:VASP的问题 已有104人参与
|
||
|
vasp是一个很不错的软件,尤其称道的是他的计算速度和稳定性。 由于wuchenwf版主即将隐退,特开此贴。 本人能力也极其有限,所以在此抛砖引玉,标题虽是答疑,但是我觉得更应该是讨论。大家术业有专攻,VASP的每个方面不可能都掌握。希望和大家也共同进步。 老帖地址(供大家参考)(wuchenwf):http://muchong.com/bbs/viewthread.php?tid=1740019&fpage=1 老帖地址(供大家参考)(九尾鱼):http://muchong.com/bbs/viewthread.php?tid=679096&fpage=1 老帖地址(供大家参考)(wuchenwf):http://muchong.com/bbs/viewthread.php?tid=1083523 |
» 本帖已获得的红花(最新10朵)
» 猜你喜欢
有院领导为了换新车,用横向课题经费买了俩车
已经有9人回复
酰胺脱乙酰基
已经有13人回复
博士延得我,科研能力直往上蹿
已经有8人回复
同年申请2项不同项目,第1个项目里不写第2个项目的信息,可以吗
已经有4人回复
有时候真觉得大城市人没有县城人甚至个体户幸福
已经有10人回复
天津大学招2026.09的博士生,欢迎大家推荐交流(博导是本人)
已经有5人回复
CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生|MSCA经费充足、排名优
已经有6人回复
面上项目申报
已经有3人回复
遇见不省心的家人很难过
已经有22人回复

★
小木虫(金币+0.5):给个红包,谢谢回帖交流
小木虫(金币+0.5):给个红包,谢谢回帖交流
|
最近跑vasp,遇到了一些问题,过来求教高手们,谢谢各位了。 (1)关于Si 体弹性模量的计算: 根据体模量的定义,我通过单体积下充分弛豫后得到一系列的不同晶格常数和对应能量值的关系(E-A),然后通过一个脚本murn.x利用Murnaghan方程来拟合得到最低能量对应的晶格常数大小(alat)以及体模量值(B0);可是发现,虽然晶格常数可以与实验值吻合的很好(误差在1%左右);但体模量却出现较的偏差(误差达10%以上)。后来也通过计算弹性常数来得到体模量,可是准确性没什么改观。。。请问,是不是有什么原因会很大程度的影响体模量的准确性。还是我拟合体模量的方法有问题,是否有其他的方法能得到较精确的模量值? 附上Si的原子优化脚本: #!/bin/sh rm WAVECAR for i in 5.15 5.20 5.25 5.30 5.35 5.40 5.43 5.455.50 5.53 5.55 5.60 5.65 5.70 do cat >POSCAR< Si-fcc $i 0.0 0.5 0.5 0.5 0.0 0.5 0.5 0.5 0.0 2 Direct 0.0 0.0 0.0 0.25 0.25 0.25 ! echo "a = $i angstrom"; time vasp E=$(grep "TOTEN" OUTCAR |tail -1| awk '{printf "%12.6f \n", $5}') V=$(grep "volume" OUTCAR |tail -1| awk '{printf "%12.4f \n", $5}') echo $i $V $E >>EtVo.dat done INCAR: SYSTEM = Si-fcc ENCUT = 450 ISTART = 0; ICHARG = 2 ISMEAR = -5 IBRION = 2; NSW = 60 ISIF = 2 EDIFF = 1E-4; EDIFFG = -1E-3 PREC = Accurate (2)关于体材料的自旋极化设置问题: 计算Co基磁性材料时,在INCAR中设置了参数ISPIN=2;可是静态自洽scf计算、能带结构band计算等得到的OUTCAR中并未记录magnetization信息;只有在态密度dos计算的OUTCAR文件里有s、p、d 以及Total的magnetization?而且在同种Co基磁性体系的不同结构中,有的结构似乎没有自旋极化( Total magnetization=0.0);有的结构却有明显的自旋极化现象(Total magnetization=2.0); 附上其中一个结构的算DOS的INCAR: SYSTEM = Co-work ENCUT = 450 ISTART = 1; ICHARG = 11 ISMEAR = 0; SIGMA = 0.1 GGA = 91; VOSKOWN =1 ISPIN = 2 RWIGS = 2.163 1.512 PREC = Accurate (3)关于LDA+U的设置: 在计算Co基材料时,由于有3d电子的存在,所以想对GGA势进行“+U”计算。可是不知道其中的U值和J值该怎么设置?(LDAU、LDAUTYPE、LDAUL、LDAUU、LDAUJ)。。vasp手册说要参考实验文献;具体要参考什么实验值?论坛上有人说带隙值,那如果是金属,怎么办?还有,如果没有实验文献值,那以什么来作为判据,选择正确的U和J值。 |
148楼2010-12-20 21:51:00
lishufei
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2640.1
- 散金: 250
- 红花: 2
- 帖子: 680
- 在线: 299.5小时
- 虫号: 867947
- 注册: 2009-10-10
- 性别: GG
- 专业: 理论和计算化学
2楼2010-04-13 21:12:08
lishufei
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2640.1
- 散金: 250
- 红花: 2
- 帖子: 680
- 在线: 299.5小时
- 虫号: 867947
- 注册: 2009-10-10
- 性别: GG
- 专业: 理论和计算化学
★ ★
xiaohunhun(金币+2):xiexie 2010-04-13 22:53
gavinliu7390:你的编译过程应该没问题。 2010-04-14 07:41
xiaohunhun(金币+2):xiexie 2010-04-13 22:53
gavinliu7390:你的编译过程应该没问题。 2010-04-14 07:41
|
说明一下我安装的是串行的安装方法如下: vasp4.6 安装全部过程 串行所需文件 l_fc_c_9.1.036.tar l_mkl_p_9.1.023.tar makefile 串行 NCOM_L_MKL__NP5T-TXVGVMK8.lic vasp.4.6.tar.gz vasp.4.lib.tar.gz NCOM_L_CMP_FOR_NC5M-MNS7LCG3.lic 1. 解压l_mkl_p_9.1.023.tar 命令 tar -xvf l_mkl_p_9.1.023.tar 进入l_mkl_p_9.1.023 命令cd l_mkl_p_9.1.023 安装 ./install.sh 以后按照安装说明进行。 2 解压 l_fc_c_9.1.036.tar 命令 tar -xvf l_fc_c_9.1.036.tar 进入l_fc_c_9.1.036 命令 cd l_fc_c_9.1.036.tar 安装 ./install.sh 以后按照安装说明进行。 注意所有安装都在root下 进入安装说明里会有提示。 3 改变量环境 退出后 进入根目录su 口令 新目录要cd 再回车 显示~ 输入 vi .bashrc 回车 进入后改成 # .bashrc # User specific aliases and functions alias rm='rm -i' alias cp='cp -i' alias mv='mv -i' PATH="$PATH:/usr/local/bin" (该部分为mpich2的路径)注意安装并行才用 # Source global definitions if [ -f /etc/bashrc ]; then . /etc/bashrc fi (以下为mkl及fortran安装后需设置的环境变量) export LD_LIBRARY_PATH=/opt/intel/mkl/9.1.023/lib/em64t:/opt/intel/fce/9.1.036/lib . /opt/intel/fce/9.1.036/bin/ifortvars.sh export LD_LIBRARY_PATH=/opt/intel/fce/9.1.036/bin/ifortvars.sh:$LD_LIBRARY_PATH 注意进入后按i才能对内容修改。 改好后按esc键 :wq 这是保存并退出。 保存后在终端执行命令source /root/.bashrc就可以了 执行没错即可。 4 进入放文件的目录 cd 解压vasp.4.6.tar.gz 和 vasp.4.lib.tar.gz tar xzvf vasp.4.6.tar.gz 和tar xzvf vasp.4.lib.tar.gz 进入vasp.4.lib目录 cd vasp.4.lib cp makefile.linux_ifc_P4 makefile 打开makefile文件修改FC=ifort 保存文件 make 编译成功后,得到libdmy.a文件。 5 退出vasp.4.lib目录 进入vasp.4.6目录 把上面的makefile拷贝到vasp.4.6里去。(这里面是修改过的不需要改动) 进行make 运行结束安装成功。如果显示某行错误可在打开makefile文件找到那行,如果前面没有#就加上#,如果有#就去掉#。在进行make。 |
3楼2010-04-13 21:14:24
gavinliu7390
木虫 (著名写手)
叶落鹰飞
- 1ST强帖: 3
- 应助: 51 (初中生)
- 贵宾: 0.481
- 金币: 1564.7
- 散金: 2975
- 红花: 34
- 帖子: 1278
- 在线: 241.4小时
- 虫号: 675997
- 注册: 2008-12-17
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和

4楼2010-04-14 07:39:56













回复此楼
Copent