科研中量子场论的计算都使用哪来的程序?

科研中量子场论的计算都使用哪来的程序?,第1张

生成费曼图:FeynArt或者Qgraf

单圈计算振幅以及张量约化:FeynCalc或者FormCalc

用符号计算系统自己计算振幅:FORM和mathematica

双圈图以上的约化:Reduze和FIRE

相空间积分:用cuba的蒙卡包(我没用过)

生成费曼规则:FeynRules

自动化树图截面:CalcHep

自动化生成树图或者单圈QCD最终截面:MadGraph,还有一堆类似的我没用过的。另外MadGraph还整合了做Parton Shower和对撞机模拟的包。可以直接得到特定标准模型的散射截面的程序包:MCFM,(还有一堆类似的)画图以及一些统计分析:ROOT管理PDF(部分子分布函数,QCD要用):LHAPDF

最简单的方法:首先用FeynRules生成费曼规则,然后再导入CalcHep或者MadGrap,给定具体过程和一些Cut,直接一步步按回车就能得到树图结果。时间上,熟练的话大概就需要一两个小时,如果从头开始学的话大概要三五天吧。比较麻烦的方法:用FeynRules生成费曼规则,再用FeynArt生成费曼图,用FeynCalc把振幅化简好,然后再自己写蒙卡程序作相空间积纤好租分。好处是可以算到一圈,虽说新物理的圈图一般不算。熟练的话,把自己之前的程序包改一改就能比较快的搞定,从毁兆头开始学可能要搞半个月。最麻烦的方法:用FeynRules生成费曼规则,FeynArt生成费曼图,用Form或者Mathematica手动化简振幅,用REDUZE或者FIRE约化积分,并想办法计算主积分,然后再自己写相空间袜姿积分的蒙卡程序。好处是可以算到两圈,虽说新物理没意义了,标准模型的两圈基本上是需要以这样的方法计算的。我已经搞了半年多了,还没出最后结果。

科学探索无止境。然而,当科学家们向太空宇宙深处探索得越深入时,似乎越是感到迷茫,心中的疑问反而变得更多了,对宇宙的认识不是越来越了解熟悉,而是感觉越来越陌生了。正如求知若渴的人学习知识一样——掌握的知识越多越深入的人,反而越觉得自己的知识匮乏不够使用。

量子力学

许多科学家们都开始依靠量子力学来了解,意识到底是什么了。在量子力学里面,有两个科学家们非常注意的实验,其中一个就是薛定谔的猫,而另一个就是延迟实验。前者相对于来说是思维上的实验,而后者则是一个真正意义上的实验,而且科学家们做了许多次。这两个实验,表面上来看都是意识扮演着非常重要敏汪的角色,甚至连观察一方的意识,也会被反过来改变,但是事实情况真的是这个样子吗?

如果我们的意识,真的有决定物质的能力,那么到了后来,他可能就会面临着一个极其严重的问题,那就是可能我们所生活做的宇宙,并不是我们之前所考虑的这一种状态。又或者来说做旅,我们仅仅是一些比较高层次的人,所编写出来的程序而已,这一种想法,真的让人们感觉到非常的桥胡仔心惊。

不要小瞧一个夸克和一个电子里的能量

如果,我们生存的地球,是被某种高智慧生命体早已通过“门槛系数”锁死,那么就等于是锁死了我们的科技文明,也就等于阻断了我们从宏观世界向微观世界探索的通道。这也是为什么我们始终无法寻找到大统一理论的原因,始终无法架起从宏观世界通向微观世界桥梁的原因。

研究者分析道:“门槛系数”越大的微观世界,所拥有的能量与相态大小是成反比例的。也就是说,微观世界的能量,有可能会是宏观世界里能量的几亿兆倍,这就是原子d能被研究出来的原因。所以,不要小瞧一个夸克和一个电子里的能量,有可能远胜于核d。或许,这也是锁死我们科技文明的根本原因——想想核d的威力吧。

当科学家需要的时候,牛顿力学被修正为相对论力学,诸如此类。宗教则相反,一个断言是百分之百正确的。当然,佛教是一个独特的宗教,哲学成分很大。即使如此,量子力学也不能拿来支持某些说法,例如世界是虚幻的,灵魂是存在的。让科学的归科学,宗教的归宗教,才是宗教和科学相处的最好方式。

专门解答量罩好子力学的软件有:vasp、CASTEP、MS、elk、ASW、abinit、Quantum-ESPRESSO、flapw

vasp由于优化算法比较好,计算速度较快,计算时问题出现较少,需要控制的参数也不是很多。而且他还有比较独特的paw势,在dft平面波的软件中,日渐趋于主流软件。其功能也在逐渐完善。发展潜力很大。

abinit计算软件,我感觉其功能还是很强大的(可以说其它第一性原理软件能计算的性质,它基本都没问题,而且它的gw和dfpt独特功能),计算速度也不是很慢。只是用起来太麻烦,控制参数繁多,入门很慢。

至于再具体的,很多时候就得具体的问题,具体分析了。

MS中包括Visualizer、CASTEP、Dmol3、VAMP、 Discover、 Amorphous Cell、Compass等多个建模和计算软件,可进行晶体、非晶电子结构的量子力学计算,也可进行分子的量子力学计算;可进行材料的分子动力学计算;可银坦进行x-ray衍射计算;能够处理稀土元素,功能强大,就是贵。有Linux和Windows版本,便于学习。 VASP 具有很好的赝势,与CASTEP相似,使用平面波基组。 Wien2k是全电子计算的量子化学软件,处理磁性材料较好。abinit、 Siesta是免费软件,提供原代码。处理重金属不准,缺乏相应的赝势。用于计算晶体的电子结构。 Gaussian主要用于分子、离子的计算,可处理激发态,精度物搏铅高,耗资源。

flapw中强烈推荐elk,主要优点:代码清晰,容易了解计算原理。后处理极为简单,像画能带图时,会自动给出高对称点,使用自带的elkband可以很容易得到能带图而不用复杂的后处理过程。同时可以处理的性质比较多。

缺点:文档太少。只有一个输入说明,不适合新手。不过官方论坛提问作者一般都会回答。再就是为方便后处理,输出文件比较多,使用前最好看下给的相关例子。再就是目前不支持mpi(可以用openmp并行),对声子不能用dfpt。

缀加球面波方法-ASW

这个软件是在量子化学网上看到的,之前一直都在使用abinit,但是苦于我要计算的体系所含元素的赝势不全,所以就尝试使用ASW。但是目前至少在小木虫上发现使用此软件的人极少。ASW程序的执行文件需要通过邮件向volker Eyert申请。

我总结ASW的特点:计算速度快;输入文件只有一个而且相对简单;磁性计算比较全面:包括无磁、铁磁、反铁磁。当然个人认为它最大的一个优点就是作者编写了很多计算和后期作图的脚本,使用很方便,特别是处理分波态密度时相当轻松。另外目前也发现了一些问题:个人感觉ASW对计算体系的结构尤其是对称性方面有很严格的限制,并不像VASP或是ABINIT那样相对宽松,还有就是它的优化功能不是很全面,可能是我使用的还不是很熟练地原因。

abinit、MS.的功能非常强大。第一性原理能计算的性能方法,基本都可以计算。但相对与VASP来说,精度方面可能需要加强。VASP计算可以结合其它的程序计算更多的性能。如结合phonopy算声子普。结合ATAT计算激发态的一些性能。

现在常用的第一性原理计算软件中最容易上手的就是MS,虽然它有很多缺点,比如说赝势不好,精度不高,源代码不开放等等,但是对于一般的科研工作而言就够用了,毕竟我们做计算的目的是寻求规律,解释现象,探求本质,而不是一味追求高精度。精度再高也是理想状态,也无法实现复杂实验条件的模拟。至于文章中图谱的效果、好看与否,更大程度上在于个人对于数据的理解程度、后续分析及数据处理,而不是软件本身了。从这个角度讲,ms则是一款比较实用的软件,把时间和精力用于软件的开发和学习还不如用来加深理论功底和数据分析!一点拙见而已。

使用Quantum-ESPRESSO中,与Abinit一样,都是开源的多功能第一性原理计算包,同样支持GW的计算,并且带有SISSA自主研发的TDDFPT,虽然现在发布的还是有很多功能限制和缺陷。另外,声子计算方面,比Abinit要简便,可以直接像给出k点一样给出q网格。事实上,DFPT方法的发明者就是Quantum-ESPRESSO的作者。另外,QE与其他一些软件包都有接口,可以协调工作,比如万尼尔方程。

计算速度也比较快,计算参数设置灵活但是不复杂,并且邮件列表里人不少,编译时也非常简单,基本不用自己设置什么参数。最大的问题就是赝势库过分不完整,很多时候只能靠着转换别的软件的赝势,或者自己生成赝势,这对不了解理论或者赝势的新手非常致命。另外不同的功能分散在不同的可执行文件中,刚开始入门时可能容易犯晕~

flapw中的wien2k也是很不错的。优点主要在于:

1、有图形界面,上手相对容易,输入和后续处理都有比较好的脚本处理,都很方便。

2、手册对各个参数介绍很全,还有mailinglist可以查询和讨论。

3、对很多物理性质直接模拟,比如光学性质,谱,声子谱等。

4、软件价格很便宜,好像是$400,可以在大型服务器上并行,处理上百个原子是没有问题的,当然计算量相对赝势程序要大。

5、常用的各种交换关联势都已经集成,GW方法已经集成只是还没有释放。

Quantum-ESPRESSO, 个人认为对初学者最大的障碍是没有好的manual. 如果以前没有用过其他的第一性原理软件对于参数的设置就比较难以理解. 不过, 确实如souledge 所说邮件列表非常活跃, 问了几次问题都能有人热心的解答. 实在不行了跟软件的作者联系, 会学到很多东西.

个人比较看好 Quantum-ESPRESSO, 作为开源软件最大的优势就是能吸纳最新的研究成果, 并且自己有更大的主动权.


欢迎分享,转载请注明来源:内存溢出

原文地址: https://outofmemory.cn/yw/12313853.html

(0)
打赏 微信扫一扫 微信扫一扫 支付宝扫一扫 支付宝扫一扫
上一篇 2023-05-24
下一篇 2023-05-24

发表评论

登录后才能评论

评论列表(0条)

保存