HPC助力超导材料发现

最新推荐文章于 2024-09-14 19:00:09 发布

集群计算技术

最新推荐文章于 2024-09-14 19:00:09 发布

阅读量590

点赞数 22

文章标签： linux 运维服务器

本文链接：https://blog.csdn.net/u013206206/article/details/139946623

版权

据新华社报道，2023年度国家最高科学技术奖6月24日在京揭晓，李德仁院士、薛其坤院士获得中国科技界崇高荣誉。

李德仁是著名的摄影测量与遥感学家，一直致力于提升我国测绘遥感对地观测水平。他攻克卫星遥感全球高精度定位及测图核心技术，解决了遥感卫星影像高精度处理的系列难题，带领团队研发全自动高精度航空与地面测量系统，为我国高精度高分辨率对地观测体系建设作出了杰出贡献。

薛其坤是凝聚态物理领域著名科学家，取得多项引领性的重要科学突破。他率领团队首次实验观测到量子反常霍尔效应，在国际上产生重大学术影响；在异质结体系中发现界面增强的高温超导电性，开启了国际高温超导领域的全新研究方向。

薛其坤的主要研究领域之一高温超导是物理界的“流量密码”，去年七月，韩国的研究团队声称合成了全国首个室温常压超导体，引起了全球的关注。除此之外本次获奖也有其他关于超导的机理研究工作。今天就来了解一下这个物理界的“顶流”-超导。

什么是超导，高温超导？

电阻起源于载流子(电子或者空穴)在材料运动过程中收到各种各样的阻尼。按照电阻率从大到小可以把材料分为绝缘体、半导体、和导体。但有一种特殊的导体，再到达某种特定温度(超导临界温度)以下的时候，将进入超导态，此时电阻为0，同时所有外磁场磁力线被排出，导致体内磁感应强度为零，此时出现零电阻态和完全抗磁性。

超导体图源：罗切特大学

1911由荷兰物理学家海克·卡梅林·昂内斯（Heike Kamerlingh Onnes）首次发现。在实验室内，他成功制备了纯度极高的汞，并将其冷却至接近绝对零度（-273.15℃）；当他对汞的电阻进行测量时，惊奇地发现，在温度降至4.2K（-268.95℃）以下时，汞的电阻竟突然消失，呈现出几乎完美的导电性能。再到1986年，瑞士物理学家卡尔·亚历山大·米勒（Karl Alexander Müller）和他的德国合作者约翰尼斯·格奥尔·贝德诺尔茨（Johannes Georg Bednorz）宣布，他们发现了一类铜氧化物超导体，其临界温度可以高达30K（-243.15℃），这一高温超导现象的发现打破了人们对超导只能存在于极低温的认知。随后科学家们也陆续发现了多个临界温度高于液氮温度的室温超导。并且在之后的很长一段时间内，科学见们都致力于研究寻找具有较高的超导临界温度的高温超导体。

超导的理论原理

目前较为成熟的常规超导微观理论(BCS理论)：金属中自旋和动量相反的电子可以配对成所谓的“库伯对(Cooper pairs)”,库伯对在晶格中可以无损耗的的运动，形成超导电流。在宏观上展示出零电阻。

电子间的直接相互作用是相互排斥的库伦力。如果仅仅存在库伦力直接作用的话，电子不能形成配对。但电子间还存在以晶格振动（声子）为媒介的间接相互作用：电声子交互作用。电子间的这种相互作用在满足一定条件时，可以是相互吸引的，正是这种吸引作用导致了“库珀对”的产生。因此计算电声相互作用是研究超导性质的重要一步。

李政道先生提议的有关BCS超导机理漫画：单翅蜜蜂代表单个电子，题曰：“单行苦奔遇阻力，，双结生翅成超导”，下面为蜂窝状的C60系列超导体。

从超导到EPW

EPW方法能够计算电子-声子相互作用方法，能够准确地描述材料中电子和晶格振动之间地相互作用。通过计算电子结构和声子谱，可以计算电子-声子相互作用的强度，EPW方法将电子-声子相互作用表示为矩阵元素的形式，并使用这些矩阵元素计算超导临界温度等性质。为超导材料的设计和优化提供了重要的理论依据。

以上提出的理论及计算方法均是针对常规超导来说，而对于高温超导体，其理论物理机制尚不明确，更无法通过DFT求解。人们虽然可以通过计算电声子相互作用，预言常规超导体的超导相变，但尚未发展出计算高温超导相变的公认方法。可见通过DFT理论验证高温超导目前还尚未有理论支撑。但对于常规超导来说，计算模拟是预测和发现新的超导的重要手段。

本着认真科普的原则，附上QE中EPW安装方法，以及使用EPW计算超导的基本教程。

首先安装QE,具体可参考公众号或之前发布的文章：

https://mp.weixin.qq.com/s?__biz=MzkxMjY5ODA1OQ==&mid=2247484636&idx=1&sn=bf4567cd170b23ca8606f46865154b9d&chksm=c109bb78f67e326ee43307b4efd2ce34b404c75060e3a0c5e5b8df4859fcd7361775f60ef03c&token=2001984167&lang=zh_CN#rd

EPW模块安装

没有外网的环境下，需要单独下载wannier90-3.1.0：

https://codeload.github.com/wannier-developers/wannier90/tar.gz/v3.1.0

在EPW同级目录install下有个plugins_makefile文件，打开我们发现wannier90部分，文件夹是命名为W90, 而解压wannier90-3.1.0.tar.gz得到的文件夹是wannier90-3.1.0,需要下载wannier90-3.1.0放在EPW同级目录，解压然后修改文件名为W90

mv wannier90-3.1.0/ W90(编译安装wannier需要是这个文件名)
cd EPW
make

wannier90可以顺利安装，但在编译epw时找不到wannier90-3.1.0文件

此时将W90重新命名回wannier90-3.1.0

mv W90 wannier90-3.1.0（但编译安装epw需要这个文件名）
make

以MgB2为例计算超导性能

计算声子

cd /install/EPW/examples/mgb2/
cd phonons/

自洽计算：

mpirun -np 20 pw.x < scf.in >scf.out

计算声子频率和特征模式：

mpirun -np 20 ph.x < ph.in > ph.out

这将为我们提供包含动态矩阵和 28 个不可约q点的自洽势变化的文件。现在我们需要将.dyn和.dvscf文件以及_ph0/diam.phsave文件夹复制到文件夹中save/。这些文件包含 EPW 需要的来自 Quantum Espresso 的所有信息。

EPW 希望根据列表中不可约q点的数量对这些文件进行标记。为了正确标记文件，您可以使用 python 脚本pp.py。如果您决定以这种方式继续，只需发出语句：

python3 /install/EPW/bin/pp.py(生成.dyn,.dvscf等save文件)

如果遇到：ImportError:No module named builtins（报错信息）等类似错误

安装python 的future模块

https://pypi.org/project/future/#files 下载解压

tar -zxvf future-0.18.2.tar.gz
cd future-0.18.2
python setup.py install

运行EPW

cd /public/software/apps/qe-6.8/EPW/examples/mgb2/epw
mpirun -np 20 pw.x < scf.in > scf.out 
mpirun -np 20 pw.x < nscf.in > nscf.out 
mpirun -np 20 epw.x  < epw.in > epw.out

生成MgB2相关文档，包含超导间隙，临界温度等数据。

详细可参考：

https://docs.epw-code.org/doc/MgB2.html

新闻来源：新京报

内容引用：

https://www.kepu.net.cn/blog/luohuiqian/201903/t20190327_475843.html

https://mp.weixin.qq.com/s/MUNcBUkanw2bA4YdPP90aw