24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

souledge

专家顾问 (著名写手)

[交流] ELK！这是卧薪尝胆，十月磨砺的节奏啊！（已上传Windows版）

ELK，作为免费的FP-LAPW代码，作为TD-DFT鼻祖级大牛实验室的作品，一直以来更新非常快，并且代码结构非常清晰，功能庞大，是难得的DFT佳品。
但是不知道是否是FP-LAPW的理论上比较难以实现突破，DFPT和Cell Relax一直都不存在（其实很久以前有人做过DFPT在FP-LAPW中的应用，但是并没有推广开来；至于Cell Relax，就从来没见过有代码说可以了，甚至于连个计算stress的LAPW代码都没听说），这点可以参考FP-LAPW最著名的商业代码Wien2k。
而ELK则在2012年5月21日发布了1.4.22后，持续一年多没有任何更新，险些以为ELK就要退出历史舞台了。
没想到，2013年6月底，ELK给我们带来了一个大的惊喜：density functional perturbation theory (DFPT) phonons now available; only for spin-unpolarised calculations at the moment; spin-polarisation will be added for the next release; use with task=205; still experimental and generates incorrect results for certain symmetries - we're still working on this; see the example 'Al-phonon-DFPT'，同时，其实还有一个：hybrid functionals should now work with libxc (experimental)。要知道，以前的ELK虽然也可以挂LibXC，但是却只能用简单的LDA和GGA，无法使用Hybrid xc。这时的ELK还是个beta版。
没过几天，在2013年7月中旬，ELK更新到了2.1系列的最终版本2.1.25，去掉了beta字样。
对于一个文章中从来离不开声子的人，看到这里感觉已经非常开心了，因为有了DFPT，内存总算不是制约计算的障碍了。本以为ELK的惊喜会持续一段时间，没想到，ELK又很快升级了，这次的升级更是高大上了。
2013年9月20日，ELK发布了2.2.1版，该版除了除去之前2.1版本的很多bug之外，最大的成就，就是引入了Cell Relax：added full structural optimisation; setting 'latvopt=1' performs unconstrained optimisation of the lattice vectors; 'latvopt=2' performs iso-volumetric optimisation; see the example 'MgB2-lattice-opt'。从此，做晶格的弛豫（也就是所谓的优化），不再只是CASTEP、VASP、Abinit、PWscf这些Plane-Wave Pseudo-Potential代码的专利，也不用再为那么多赝势究竟选择哪个头疼，直接在ELK中进行优化即可了。或者说，如果对一些优化的结果有疑问，完全可以以ELK的结果作为一个参考。要知道，FP-LAPW的准确性可是非常有名的，不管是VASP还是Abinit和PWscf，出点PAW的赝势都是要跟FP-LAPW对比的。只是收费的VASP可以与Wien2k比，而免费的Abinit和PWscf只能跟ELK比。
PWPP方法的代码现在大行其道，占据了DFT文章的绝大部分，这一方面是由于号称Accurate DFT的LAPW速度相对较慢，也是由于LAPW代码的功能太少。现在计算机硬件的性能不断提升，ELK又很争气的加了很多PWPP中才有的功能，相信以后只用ELK就做全套文章的想法也很快可以实现。
弟兄们，还等什么，来学习ELK吧！下载地址：http://sourceforge.net/projects/elk/files/?source=navbar

上传了最新的ELK 2.2.5，已经使用MinGW自带的gfortran和gcc编译成支持32-bit Windows的可执行文件（需要把dll文件放在同文件夹或者C:\Windows下）。支持OpenMP。由于MinGW本身的一些问题，运行效率貌似没有直接用Intel Fortran for Windows编译出来的高。而且使用的BLAS和LAPACK库也都是自带的，不如MKL的给力。先凑合着测试一下吧，将来会把更好的编译版发上来的。[ Last edited by souledge on 2013-10-30 at 11:06 ]

回复此楼