24小时热门版块排行榜    

CyRhmU.jpeg
查看: 577  |  回复: 5
当前主题已经存档。

tyronezhang

铁杆木虫 (正式写手)

[交流] 【求助】200金币求vasp算介电

200金币求vasp 4.6.26利用OPTICS算介电的详细教程,包括OPTICS的编译和使用步骤,最好有一个例子,或者有人愿意用vasp5.2代算一下介电常数亦可,5原子小体系很快的,急,有意者回复。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

顺水

金虫 (小有名气)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zxzj05(金币+1,VIP+0):鼓励讨论! 9-28 17:09
zxzj05(金币+20,VIP+0): tyronezhang 转给 9-28 17:13
http://muchong.com/bbs/viewthread.php?tid=1475470关于从vasp4.5、4.6计算光学性质(组合帖)
2楼2009-09-28 10:12:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

apzhough

铁虫 (小有名气)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zxzj05(金币+1,VIP+0):鼓励讨论! 9-28 17:10
一、关于VASP中光学optics.F程序的编译(ifort+mkl)

原帖地址http://new.quantumchemistry.net/ ... ion/200512/187.html

从VASP 4.4.4和VASP.4.4.5开始,VASP支持简单的光学性质的计算。其具体实现程序是optics.F。其中一些必需的后处理程序,可以向Jürgen Furthmüller发信要。

   此程序可以使用ifort和mkl编译出来,用pgi编译它说optics.f语法错,不知道为什么。然后用编译成功的程序处理已经得出OPTIC文件的一个目录,满心欢喜,结果给我了这些东西:
   scanning POSCAR, KPOINTS and OPTIC
   reading KPOINTS and OPTIC
   set up transition energies and amplitudes
   DOS and joint DOS
   forrtl: severe (174): SIGSEGV, segmentation fault occurred
   Image PC Routine Line Source
   optics 0804A8DE Unknown Unknown Unknown
   optics 080502A7 Unknown Unknown Unknown
   optics 0804A594 Unknown Unknown Unknown
   libc.so.6 00314E33 Unknown Unknown Unknown
   optics 0804A4A1 Unknown Unknown Unknown
   forrtl: severe (174): SIGSEGV, segmentation fault occurred
   Image PC Routine Line Source
   libc.so.6 00360ECE Unknown Unknown Unknown
   libc.so.6 003614BA Unknown Unknown Unknown
   optics 08086FFD Unknown Unknown Unknown
   optics 0806A899 Unknown Unknown Unknown
   optics 0806C183 Unknown Unknown Unknown
   optics 0806B296 Unknown Unknown Unknown
   optics 0806D634 Unknown Unknown Unknown
   libpthread.so.0 0054A8A0 Unknown Unknown Unknown
   optics 080502A7 Unknown Unknown Unknown
   optics 0804A594 Unknown Unknown Unknown
   libc.so.6 00314E33 Unknown Unknown Unknown
   optics 0804A4A1 Unknown Unknown Unknown
   貌似段错误。

   具体编译方法是先编译optics上一级目录中的文件,得到prec.o constant.o error.o tet.o
   ,并把它们拷贝到optics目录下。
   然后再到optics目录下编译optics
   其中makefile的头几行修改成如下:
   # THIS IS AN EXAMPLE FOR INTEL-LINUX PLATFORMS RUNNING PGF COMPILERS. CUSTOMIZE!
   
   .SUFFIXES: .inc .f .F
   
   MAKE = /usr/bin/make
   CPP_ = /usr/bin/cpp -P -C
   CPP = $(CPP_) \
   -Dkind8 -DNGXhalf -DCACHE_SIZE=1000 -DPGF90 -Davoidalloc <$*.F >$*.f
   FC = pgf90 -Mfree
   DEBUG = -g -O0
   OFLAG = -O2 -tp p6
   DOUBLE = -r8
   SINGLE =
   TIMING =
   #TIMING = -pg
   
   # binding
   BIND =
   LIB = -L../../vasp.4.lib -ldmy -L/usr/local/lib -llapack -lblas
3楼2009-09-28 15:01:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

apzhough

铁虫 (小有名气)


zxzj05(金币+1,VIP+0):鼓励讨论! 9-28 17:10
编译出了optics,问题已经不大了
首先到http://www.freeware.vasp.de/VASP/
去下载最新版的optics和其他的一些东西。注意,其中的一些F程序似乎有语法错误,起码在pgf90编译的情况下提示出错,它们是:
smear.F的第29行,SUBROUTINE中的S写到了第6列,同一行的换行符应该写在第30行第7列。
error.F的第8行还有下面一个类似的地方,&字符越界,应该换到下一行做换行符用。
symlib.F缺少了sym_prec模块,添加进去。
MODULE sym_prec
  USE prec
  REAL(q) :: TINY=1.E-5_q
END MODULE
然后我很顺利地用pgf90加上pgf90自己的blas和lapack编译出了vasp。然后手动编译error.F,prec.F,smear.F,就用pgf90 -c smear.F这样的格式即可,产生相应的.o文件。然后进入optics目录,make clean一下后,修改makefile,用和编译vasp一样的blas和lapack来编译optics,成功。但是此时的vasp.4.6目录下的 vasp是有问题的,执行起来没有任何输出,让我伤透了脑筋。用原先存着的ifort编译的vasp计算后得到的OPTIC二进制文件optics程序执行后介电函数实部很多nan。还是用以前pgf90编译的vasp运行,这才没有了问题。

optics程序需要用paw的势,可以用USPP计算其他的东西,当要计算光学性质的时候改用pgf90编译的vasp和paw势即可。

OPTIC计算出来后,一定要记得cp IBZKPT KPOINTS。

程序规定要将ISMEAR设为采用Tetrahedron k mesh。因此,k点数起码要超过4个。我计算Se supercell用的1*1*2就不行。真正要计算的时候要将K点取大。

目前来看得到的结果和CASTEP的结果类似。
4楼2009-09-28 15:01:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

apzhough

铁虫 (小有名气)


zxzj05(金币+1,VIP+0):鼓励讨论! 9-28 17:10
先用ISMEAR=-4 or -5,INCAR里面设置LOPTICS=.TRUE.即可。注意要用paw的势才行。
然后cp IBZKPT KPOINTS,建立一个叫OPTCTR的文件,里面写一些控制参数,一个很简单的例子是:
OMMIN = 0
OMMAX = 20
NEDOS = 1000
ISYMM = 1
LJDOS = .TRUE.
LDOS  = .FALSE.
LKRAMERS = .TRUE.
GAMMA = 0.2
最后运行那个optics程序即可。
5楼2009-09-28 15:02:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tyronezhang

铁杆木虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
zxzj05(金币-20,VIP+0):鼓励讨论! 9-28 17:12
引用回帖:
Originally posted by 顺水 at 2009-9-28 10:12:
http://muchong.com/bbs/viewthread.php?tid=1475470关于从vasp4.5、4.6计算光学性质(组合帖)

辛苦了!
资料很多,这些资料我也查到过,觉得有点乱。请版主先从我的帐号下转20金币给顺水,谢谢。
我是希望能有详细的步骤(包括编译的),如果能解决问题,再追加金币。
6楼2009-09-28 16:12:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tyronezhang 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见