| 查看: 12345 | 回复: 33 | |||||||||||
| 【奖励】 本帖被评价19次，作者顺水增加金币 17.5 个 | |||||||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||||||||
[资源]
【转帖】关于从vasp4.5、4.6计算光学性质（组合帖）
|
|||||||||||
|
一、关于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 ----------------------------------------------------------------------- optics.F与vasp的版本无关，我4.5.5到4.6.25都可以编译的过去。因为它link时，需要的 只是vasp.4.lib中的文件，而vasp.4.lib中的文件基本是不升级或更新的。 [ Last edited by wuli8 on 2009-8-12 at 18:00 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 仿真建模与计算 | 第一性原理计算资源总汇 | VASP | Optical |
| 计算-vasp | first principle | vasp光学 | 模拟计算课题的学习 |
| 光学性质计算 |
» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有17人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
基金不中，共勉
已经有5人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有3人回复
-
国自然面上复盘~欢迎讨论
已经有5人回复
-
我不理解！
已经有23人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有10人回复
-
系统查不到
已经有10人回复
-
怎么查啊
已经有5人回复
-
范进中举一文的中心思想
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- VASP中如何查看正在运行程序的源文件
已经有3人回复
- vasp5.2计算光学性质
已经有22人回复
- 介电常数计算方法
已经有13人回复
- vasp5.2光学计算中怎么加入scissor operation?
已经有13人回复
- vasp算出来的文件基本都是0字节的
已经有24人回复
- 同个任务同样的节点数VASP4.6能算，而VASP5.2不能
已经有9人回复
- vasp编译通过，前台算ok，去qsub到后台，找不到libmkl_lapack.so
已经有19人回复
- 【求助】为什么vasp5.2版本用不了ave-mesh
已经有8人回复
- 【求助成功】vasp5.2并行安装不上
已经有18人回复
- 【求助】关于vasp5.2计算色散力
已经有5人回复
- 【求助】vasp5.2光学性质问题
已经有18人回复
- 【求助完毕】IBZKPT,KPOINTS与DOS计算
已经有4人回复
- 【求助】vasp5.2用B3LYP计算
已经有21人回复
- 【求助】PBE0 计算Si、Ge能带结构
已经有9人回复
- 【原创】用pgi10.0编译vasp5.2
已经有12人回复
- 【求助】vasp5.2异常退出的问题
已经有16人回复
|
optics程序incompatible number of kpoints错误 转自：http://new.quantumchemistry.net/ ... .asp?ArticleID=1143 运行vasp的计算光学性质的程序optics的时候，虽然K点是用的ismear=4之后生成的IBZKPT里面的list,也生成了OPTIC文件。但是运行optics程序时有时候会出现以下的错误：incompatible number of kpoints in file KPOINTS and OPTIC。OPTCTR也有。 出现此错误的原因是vasp和optics程序的big-endian/little-endian模式不同。所谓的big/little endian是一种存储模式，代表高位字节存在高位地址或低位地址里面。在optics和 vasp的 makefile 里面都有这样的编译选项 -byteswapio。 所以解决方法就是只需在vasp和optics两个makefile里面指定相同的模式即可。具体方法为在optics的makefile里面加上或去掉 -byteswapio 选项。要么都有-byteswapio，要么都没有，总之要保持一致。 [ Last edited by 顺水 on 2009-8-9 at 13:15 ] |
3楼2009-08-09 13:10:52
|
OPTIC 编译2 编译出了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的结果类似。 |
2楼2009-08-09 13:10:11
4楼2009-08-09 13:11:16
6楼2009-08-09 13:15:57
