| 查看: 6852 | 回复: 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» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有8人回复
-
存款400万可以在学校里躺平吗
已经有28人回复
-
最失望的一年
已经有11人回复
-
求推荐英文EI期刊
已经有5人回复
-
请教限项目规定
已经有4人回复
-
国自然申请面上模板最新2026版出了吗？
已经有20人回复
-
26申博
已经有3人回复
-
基金委咋了？2026年的指南还没有出来？
已经有10人回复
-
基金申报
已经有6人回复
-
疑惑？
已经有5人回复
» 抢金币啦！回帖就可以得到:
-
南佛罗里达大学化学系刘文奇课题组 2026 Fall 招收有机/超分子方向博士生
+2/954
-
南京大学能源与资源学院-景旭东教授 (英国皇家工程院院士) 团队博士后招聘
+1/473
-
专任教师招聘
+1/177
-
上海交通大学变革分子学中心申涛课题组2026秋季入学推荐-考核制博士招生（有机）
+1/80
-
哈尔滨工业大学王东博课题组/中科院上海微系统所梁丽娟课题组招收2026年博士生1名
+1/78
-
哈尔滨工业大学王东博课题组/中科院上海微系统所梁丽娟课题组招收2026年博士生1名
+1/70
-
211双一流北京工业大学招计算机、AI、自动控制、电子信息等方向博士生（长期有效）
+1/40
-
国重点实验室双一流A类长江学者团队招2026年全日制博士1-2名
+2/34
-
中国地质大学（武汉）分析地球化学团队招收博士生1名、硕士生3名
+2/32
-
华中农大生物催化组招26年博士：代谢工程、酵母工程、药学、有机、酶工程、计算模拟
+1/30
-
[长期合作招募] 同济大学肖倩老师团队诚邀港澳学者学术交流
+1/29
-
沈阳征女友
+1/21
-
专科起点双非硕士，申博求建议
+1/21
-
2026申博自荐 本硕双一流学科 纳米药物递送方向 一篇一区TOP 两个国家奖学金
+1/15
-
电子科技大学「基础与前沿研究院」文明健老师课题组招收博士
+1/12
-
北京工业大学高靓教授课题组2026级博士研究生招生
+1/10
-
2026 年南方医科大学基础医学院李琳课题组招收“申请-审核” 博士研究生
+1/5
-
山东大学（青岛校区）招博士后（COF\MOF\催化\电池）
+1/3
-
Plant communications难度如何？
+1/3
-
澳门科技大学诚招2026年秋季全奖博士研究生（纳米/水凝胶生物材料方向）
+1/2
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
回复
