| 查看: 7125 | 回复: 36 | ||||||||
[交流]
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 ] |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : ELK_2.2.5_MinGW_Win32.zip
2013-10-30 11:02:28, 1.7 M
» 收录本帖的淘帖专辑推荐
 第一性原理 | 第一原理资料汇编 | castep | 第一性原理 |
| DFT calculation |
» 本帖已获得的红花（最新10朵）
hezhu» 猜你喜欢
复试调剂，一志愿南农083200食品科学与工程
已经有4人回复
-
一志愿太原理工安全工程300分，求调剂
已经有3人回复
-
284求调剂
已经有11人回复
-
食品工程专硕求调剂
已经有3人回复
-
324求调剂
已经有7人回复
-
一志愿武汉理工，总分321，英一数二，求老师收留。
已经有4人回复
-
287求调剂
已经有7人回复
-
325求调剂
已经有5人回复
-
343求调剂
已经有4人回复
-
求调剂推荐 材料 304
已经有4人回复
» 抢金币啦！回帖就可以得到:
-
【强推｜可留校】济南大学博士后招聘（水利工程/铁锰碳环境过程）
+1/290
-
南方医科大学发育生物学教研室夏来新教授课题组审核制博士/硕博连读招生启事
+2/78
-
[b]【材料表征】成分含量丨微观形貌丨分子结构丨材料加工丨预存服务[/b]
+1/72
-
欢迎环境科学与工程、水利工程专业研究生加入北师大珠海校区董越课题组
+2/56
-
浙江省人民医院博士后岗位招聘
+1/41
-
安徽工程大学化学与环境工程学院对外招收2026级调剂硕士研究生
+1/40
-
烟台大学精准材料高等研究院26年材料方向研究生招生
+1/39
-
上海大学 “生物有机电子材料及器件”团队硕士、博士研究生招聘
+1/36
-
中国计量科学研究院长度量子基准实验室2026年调剂招生
+1/35
-
江西水利电力大学接收调剂08工学皆可
+1/14
-
【博士招生】天津理工大学国家杰青王铁课题组招收2026年博士研究生
+1/9
-
南京林业大学化学工程学院柏惺峰课题组诚聘科研助理、水杉英才E类/F类青年教师或博后
+1/7
-
湖北大学材料学招26级硕士研究生若干名
+1/5
-
双非有经验之人求调剂
+1/5
-
上海交大招收26年级秋季入学博士生，CET-6过，生物医用材料方向（有高分子合成基础）
+1/4
-
宁波东方理工大学余鹏课题组2026年上海交大联培博士生、博士后、科研助理招聘
+1/3
-
化学工程与技术 国家双一流学科 211高校 国家重点实验室 博士研究生
+1/2
-
2026年北京石油化工学院环境学科-水污染控制工程方向课题组欢迎您
+1/2
-
【博士招生】上海交通大学机动学院-招2026年秋“申请-考核制”博士生1名-3D打印方向
+1/2
-
中国人民大学物理学系蔡鹏课题组招收优秀硕士博士研究生
+1/1
2楼2013-10-26 01:12:07
怎样？觉得我有潜质不？ 6楼2013-10-26 10:16:08
|
Windows哪儿有什么更好的人机对话。 对我自己而言，在Windows下编译是种游戏，尤其用不支持POSIX的MinGW，编译PWscf时会有一些代码无法通过编译，这时自己去修改代码达到编译的目的，这不是很好玩么？ 顺便，也给那些想试试除Materials Studio和MedeA外Linux下DFT软件，却始终胆小或者发懒（不管是脑子还是手脚了）的人一个尝试的机会。 话说，双系统的一台机器，用ifort编译的ELK在OpenMP环境下跑，于Server 2012 DC和CentOS 6.4下速度基本没什么区别（稍微有点，可能是因为Windows版的FFT用的是ELK自带的，但是Linux版的是用的MKL的FFTW3）。不过还是说回来，MinGW编译的ELK的效率真心低了那么一点点，不过好处是，很容易把LibXC加上去～ |
29楼2013-10-30 21:55:42
★ ★
franch: 金币+2, 鼓励交流。 2013-11-04 16:41:49
franch: 金币+2, 鼓励交流。 2013-11-04 16:41:49
|
其实比较VASP和Abinit的计算性质差异是比较不科学的，因为两者的赝势从来没有互通过。说白了，除了VASP开发组的人，这个世界上没人知道VASP赝势的生成参数及赝势结构。所以拿着VASP跟Abinit比就不免有点apple比orange的意思了。 至于能比的话，用Abinit和PWscf到是可以比一比，用完全相同参数生成的NCPP，然后完全相同的计算参数。 而ELK和Wien2k，有赝势问题么？MT半径和valance state都可以自己修改的，修改成一样再比就可以了。 [ Last edited by souledge on 2013-11-3 at 18:11 ] |
36楼2013-11-03 18:10:07
打错字了……是2013年6月底…… 3楼2013-10-26 02:21:43
4楼2013-10-26 09:49:28
5楼2013-10-26 09:56:48
7楼2013-10-26 10:54:00
8楼2013-10-26 11:10:08
9楼2013-10-27 08:00:20
10楼2013-10-27 13:39:29
★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
franch: 金币+2, 鼓励交流。 2013-10-28 21:40:38
souledge: 金币+1, 这个ZnO的弛豫不成功的案例可以单独开贴上elk.in讨论嘛~看起来比较有趣。 2013-10-30 10:48:12
小木虫: 金币+0.5, 给个红包，谢谢回帖
franch: 金币+2, 鼓励交流。 2013-10-28 21:40:38
souledge: 金币+1, 这个ZnO的弛豫不成功的案例可以单独开贴上elk.in讨论嘛~看起来比较有趣。 2013-10-30 10:48:12
|
最近一直在用，现在scf收敛速度快了很多，也稳定了很多，我试着用晶格优化功能，优化zno晶胞一直没有成功过，估计对六方的优化还需要改进，立方的似乎不错。另外，我用TB09算zno能带也一直没有办法收敛，不知道楼主有没有好的建议？ [ 发自手机版 http://muchong.com/3g ] |
11楼2013-10-27 14:00:11
12楼2013-10-27 14:13:02
14楼2013-10-28 08:26:10
15楼2013-10-28 09:43:55
16楼2013-10-28 11:59:56
17楼2013-10-28 12:17:41
18楼2013-10-28 18:02:54
19楼2013-10-28 21:02:43
20楼2013-10-29 17:06:20
21楼2013-10-29 17:09:48
22楼2013-10-29 21:17:11
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
souledge: 金币+1, 基本的编译参数，这个显然是要自己下点课外功夫去学习的东西……其实没什么好讲的，左边都是些宏而已，右边是宏的值和编译器的参数~ 2013-10-30 10:50:04
小木虫: 金币+0.5, 给个红包，谢谢回帖
souledge: 金币+1, 基本的编译参数，这个显然是要自己下点课外功夫去学习的东西……其实没什么好讲的，左边都是些宏而已，右边是宏的值和编译器的参数~ 2013-10-30 10:50:04
|
看了楼主的帖子后，很是激动，就决心加入Elk队伍和大家一起探讨了（尽管用过vasp，但对于Elk目前还是小弱啊），欣喜地告诉楼主，此楼下载的时候，版本号已经更新到2.2.5了。但是在make.inc这一步遇到了一个初级问题，该问题半个小时前已解决，并把解决方案回复到“求助elk并行编译！”帖子中，但是希望路过者看到的时候，顺便说说如下make.inc每项的含义，如何？ MAKE = make F90 = gfortran F90_OPTS = -O3 -ffast-math -funroll-loops -fopenmp F77 = gfortran F77_OPTS = -O3 -ffast-math -funroll-loops -fopenmp AR = ar LIB_SYS = LIB_LPK = lapack.a blas.a LIB_FFT = fftlib.a [ Last edited by panzky on 2013-10-30 at 10:18 ] |
23楼2013-10-29 22:11:47
24楼2013-10-30 07:18:22
25楼2013-10-30 07:31:18
26楼2013-10-30 11:16:53
27楼2013-10-30 14:03:32
28楼2013-10-30 14:04:43
30楼2013-10-31 06:18:31
31楼2013-10-31 13:06:10
32楼2013-10-31 13:08:03
33楼2013-11-01 17:25:17
34楼2013-11-03 18:03:38
35楼2013-11-03 18:05:56
37楼2014-02-26 02:54:28
简单回复
sg1840892613楼
2013-10-27 22:51
回复
