| 查看: 8668 | 回复: 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» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有19人回复
-
国社科又开始会评了，不知道这次命运如何
已经有9人回复
-
能否申诉？
已经有5人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
-
面上意见出来了
已经有18人回复
-
找导师
已经有7人回复
-
29号明天会评吗
已经有3人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有11人回复
-
我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？
已经有11人回复
» 抢金币啦！回帖就可以得到:
-
面上中了，散金币。
+5/1545
-
云南大学丁军桥团队2027年博士招生
+1/947
-
八年奋战，面上中了，同时一个读了8年的博士（中途接手）答辩通过了，都是涉险过关
+2/382
-
东南大学能源与环境学院-林日琛教授课题组-27级博士招生
+1/96
-
东南大学能源与环境学院-林日琛教授课题组-27级博士招生
+1/96
-
东南大学能源与环境学院-林日琛教授课题组-27级博士招生
+1/95
-
东南大学能源与环境学院-林日琛教授课题组-27级博士招生
+1/92
-
原子层沉积（ALD）磁控溅射PECVD等微纳代工服务：18817872921
+1/92
-
很认可这句话
+1/84
-
今日少年明日老。
+1/81
-
深圳大学刘翼振课题组诚聘博士后（CRISPR分子诊断）
+1/43
-
2024级硕士求推荐博导-计算模拟背景
+3/35
-
北京师范大学（珠海）化学传感与成像课题组2026年研究生（统招/推免）招生
+1/13
-
中国科学院广州生物医药与健康研究院课题组招收博士研究生1名（2027年9月入学）
+1/7
-
香港大学智能制造实验室招收27年春秋季入学全奖博士和硕士研究生
+1/4
-
2026年工程管理与智慧城市国际会议 (EMSC 2026)
+1/4
-
香港大学先进制造实验室招收全奖博士和硕士研究生
+1/3
-
项目助理招聘启事
+1/2
-
专利审查意见答复，是AI辅助生成中很难的一个功能...
+1/1
-
复旦大学博士招生：Micro-LED光互连芯片设计（2026年11月截止）
+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
回复
