| 查看: 619 | 回复: 5 | |||
| 当前主题已经存档。 | |||
tyronezhang铁杆木虫 (正式写手)
|
[交流]
【求助】200金币求vasp算介电
|
||
| 200金币求vasp 4.6.26利用OPTICS算介电的详细教程，包括OPTICS的编译和使用步骤，最好有一个例子，或者有人愿意用vasp5.2代算一下介电常数亦可，5原子小体系很快的，急，有意者回复。 |
回复此楼» 猜你喜欢
一志愿211，化学学硕，310分，本科重点双非，求调剂
已经有13人回复
-
285求调剂
已经有8人回复
-
22408 266求调剂
已经有6人回复
-
264求调剂
已经有8人回复
-
材料工程专业日语生求调剂
已经有9人回复
-
297求调剂
已经有15人回复
-
277、学硕，求调剂 数一104，
已经有12人回复
-
机械专硕273请求调剂
已经有3人回复
-
080100力学316求调剂
已经有3人回复
-
288环境专硕,求调材料方向
已经有39人回复
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zxzj05(金币+1,VIP+0):鼓励讨论！ 9-28 17:09
zxzj05(金币+20,VIP+0): tyronezhang 转给 9-28 17:13
小木虫(金币+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
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zxzj05(金币+1,VIP+0):鼓励讨论！ 9-28 17:10
小木虫(金币+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
★
zxzj05(金币+1,VIP+0):鼓励讨论！ 9-28 17:10
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
5楼2009-09-28 15:02:34
tyronezhang
铁杆木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 5839
- 散金: 1065
- 红花: 6
- 帖子: 935
- 在线: 204.9小时
- 虫号: 293227
- 注册: 2006-11-04
- 专业: 催化化学
6楼2009-09-28 16:12:43
