&!--百合厨SAM--&的推荐 | LOFTER（乐乎）
LOFTER for ipad —— 记录生活，发现同好
&!--百合厨SAM--& 的推荐
{list photos as photo}
{if photo_index==0}{break}{/if}
品牌${make||'-'}
型号${model||'-'}
焦距${focalLength||'-'}
光圈${apertureValue||'-'}
快门速度${exposureTime||'-'}
ISO${isoSpeedRatings||'-'}
曝光补偿${exposureBiasValue||'-'}
镜头${lens||'-'}
{if defined('posts')&&posts.length>0}
{list posts as post}
{if post_index < 3}
{if post.type == 1 || post.type == 5}
{if !!post.title}${post.title|escape}{/if}
{if !!post.digest}${post.digest}{/if}
{if post.type == 2}
{if post.type == 3}
{if post.type == 4}
{if data.msgRank == 1}{/if}
{if data.askSetting == 1}{/if}
{if defined('posts')&&posts.length>0}
{list posts as post}
{if post.type == 2}
{if post.type == 3}
{if post.type == 4}
this.p={ dwrMethod:'querySharePosts', fpost:'1cbe378b_688c0a0',userId:4662138,blogListLength:19};百合网上的花篮骗子，化名：谢永辉，46，广洲人，在深圳开什么伟发家电城。1-8号那天骗网友说新店开张，要其送花篮祝贺，先后骗得钱财16600元！晚上网友想着白天的所有事件过程，发觉不对，马上查看百合网的安全信息，才发觉自己上当了！1-9号早上还狂得发信息给网友问早安！到目前为至，这个骗子的两个联系电话仍开通着，想来它不急着关掉，可能还在行骗，因此天涯的朋友们赶快行动起来，抓住这个恶棍吧！还社会一个清静之地！　　骗子电话为：（仍通着）　　
（仍通着）　　
（此号暂时关机，不过相信它还是会开机）　　
（此号为花店电话，暂关机）
楼主发言：1次 发图：0张
　　我的朋友也被骗了，手法一模一样，深圳的
　　花店电话　　在网易同城约会里叫相信真爱
　　我一个好朋友也被骗了，而且骗子的名字也叫谢永辉，是在网易同城约会里认识的，善良的她被死骗子骗走了10万几，全部积蓄都让骗子骗光了，朋友自被骗后，一直无法平复心情，懊悔、自责、难过。希望大家对这个人渣进行人肉搜索，把他揪出来，不要让他继续骗下去了。死人骗子谢永辉你这个猪狗不如的人渣，你骗到钱也不会有好下场的，你在做，天在看，这不义之财肯定用回在你身上的，你就天天拿这钱吃药吧，还有你那个父亲、姐、花店人都不得好死，你们就如丧家之狗，天天躲藏在一个出租屋里，见不得光。不是不报，时候未到，相信你们的日子不多了。　　骗子谢永辉电话：，上海号码　　花店电话：。深圳号码　　帐号：农行：1193610
王虎　　建行：0168174
　　真不懂怎么那么容易受骗，如果自己目的单纯怎么会被骗呢。　　
　　我也是前几天在百合网上认识的自称叫袁利文的人渣给骗了，刚开始对我嘘寒问暖，关怀备至，也是跟上面朋友说的一样，家具城在香港开业需要庆典，让我以家人的身份送花篮，这样才有面子，还让他爸爸跟我通话（其实就是他本人，用了变音器），感觉他家里都很在乎你，还把他姐姐都搬出来，说家里人都非常喜欢我，我感动不已，我一天前后给花了6万多，还说我给他家里争光了，到了第二天又说家具店门口需要摆2个石狮子，我买母的，他爸买公的，我这才感觉不对劲，他随后就关机了。我真是欲哭无泪呀，明明是想真心实意的找个伴侣，没想到骗子太多了，希望得到警方的重视。　　也怪自己就是太单纯，太不理智，当时自己脑子就跟短路了一样，也不听支配了。　　骗子袁利文电话： 宅电：3 香港宅电：02 香港花店：60，这些电话估计都是网络的虚拟电话，都处于关机状态。我感觉骗子都是一伙的，希望警方尽快抓住他们，让那些在百合网上征婚的单纯男女少受损失。据我分析：自称自己是深圳，广州，香港，台湾的骗子居多，一定要提防上当受骗！
　　如果谁发现线索，互相通知举报，共同抓住骗子，我的qq号是
　　如果一个男的对你很好，你千万不要激动，这种人会让你先爱上他，然后以公司开业，什么庆典为由，让你送花蓝，这个花蓝一个就要几百几千，所以没有见面前，一定不要钱付出，在乎你的人是不会让你还没有见面就出钱的，请大家小心为好，有一些恨嫁的或者失婚女性，是他们攻击的对象他们演戏很真，所以请大家不要上当了！下面我就将在我身边的真人真事给大家描述出来，给各位恨嫁的女孩子提个醒！　　我一个朋友，暂时喊她叫做A吧！A是一个很重感情的女孩子，但是由于长相不是特别出众，加上即将奔三，一直没有找到男朋友，家里又一直催她，所以她注册了好几个婚恋网，想在网上找寻一段美丽的爱情，但是一直没有找到，所以很是苦恼，由于我人比较好说话，这个A什么都愿意和我分享，所以周伟华这个男人利用婚恋来做得“花篮骗局”，我就描述出来给各位恨嫁的美女们一个经验和教训，所以特意写下来和大家共享，希望能帮各位姐妹们安全交友。文采有限，请各位看官见谅。　　故事是这样开始的：A无意中发现“赶集婚恋网”这个网站免费交友，又是大牌明星打的广告。于是注册了账户把自己的个人资料和照片如实的挂了上去。满怀希望等待着爱情来敲门。结果七夕，果真收到不少信件。其中一个叫做“疼你”的男人他真诚的告白触动了A这颗恨嫁的心。　　一、相识　　我查看了他的注册资料，广州佛山人，未婚， 所在行业： 营销/市场/策划 ，但独白比较真诚，“我希望遇到了生命中的另一半，她是我的一个梦；为了她我愿承受一切。我会努力让这个梦成为现实，不管以后的路有多难走，我都会陪在她身边，和她一起散步，一起白头，一起看夕阳。“感觉照片的挺帅气的，就给他回复了，还告诉了对方手机号码（一般骗子都没耐心陪你慢慢聊的）。当时联系的时候，有告诉对方扣扣，但是对方以忙为理由，直接电话联系。　　二、进展　　对方知道A手机号码之后，当天晚上A看到几个未接电话，因为是外地的就没有接，第二天早上，就收到周伟华的短信，嘘寒问暖。至此每天中午13点左右，他就会发送短信问候A；晚上8-9点又来电话，聊聊家常，A也就很兴奋的和他聊天，看到这样帅气的照片，任谁都会有点眩晕。再加上A的婶婶来A的城市看望她，也听到了这个男的和A聊天，和周伟华通过话，觉得这个男的很懂事，A更加放松了警惕，周伟华对A说 ，不要告诉父母，他们是网上认识的，否则别人说一些乱七八糟的话，影响他们之间的感情（后来想想，这个周伟华还真的有反侦察的手段，他自己都知道一说网恋，可信度都很低）　　几天的短信和电话交往之后，他的细心、体贴征服了A，嘘寒问暖确实能够对女方产生不一样的感觉，原来女孩子真的喜欢听甜言蜜语，哪知道是包着糖衣的毒药，A便情不自禁对他产生了好感，觉得他就是A今生要找的人。相信能够给A带来真爱。感觉自己是个公主。熟悉之后，就给A承诺，说自己开了第二家电器店，是世界品牌格力空调专卖店，说是过几天就开张了，A一听，这个人真有本事。再加上对方还发一张在家里的照片，更加相信不已。　　　　几天后周伟华说筹备自己的新店开张，在开张前一日把他的父母接过来，并把A的照片拿出来给家人看。他说家人看过我的照片后感觉非常好,希望能够尽快带A去见他们，期间周伟华还经常发短信 给A，其中一条是：”认识你我很开心，因为这些年都是我一个人孤独走来的，所以一直都希望能够找一个全心全意爱我的人和我爱的人，当然我希望那个人是你。”A还经常在我面前晒幸福，看了这些话，这让我完全没有被骗的感觉，何况是沉浸在幸福里面的A，果然爱情让爱因斯坦退化为猩猩这句话还是有道理的。　　到了新店开张这天，这个自称周伟华的人一大早就打电话给A.告诉A说早上和客人们喝早茶的时候，他感觉好尴尬，因为他爸爸宣布了两件事，一件是感谢亲朋好友来参加格力店开张，第二件是周伟华已经找到了对象希望不久之后他们能来喝儿子的喜酒。然后他又说，为了让大家知道我的存在，父亲希望我也给他送开张花篮。给A说自己在广州番禺南沙区黄阁镇黄阁大道36-38号开了一家世界品牌格力专卖店，新店开张让A送花篮，而且规定A不能送单数，只能是双数，还说他父亲希望摆在最显眼的位置，A当时和我说，我就感觉不对劲了，面都还没有见，怎么就可以有经济往来了。　　于是A继续听这个叫周伟华的人说他的生意伙伴和亲戚都是在广州这边一家和他们比较熟悉的花店电话订购的。我顾虑到汇款的问题，我说在这边不好付款，要不你帮我先垫着回来之后我把钱补偿给你，算我的。结果他说按照风俗开张的时候金钱只能入不能出。说钱下次一定会带来给我，他经济条件不差，到时候他父母会给我留个大红包，不可能要我出钱的，还说今天早上，他妈妈已经给A买了一个玉手镯，（A当时听后和我说，早上电话的时候，还不到八点钟呢！不可能有玉器店开门啊，疑点）周伟华继续解释今天是个特别的日子，需要讨个好彩头。他只有给我打电话了，还说他父母也去买了几个花篮送给他，还把他们那边的花店的号码（）给A，还说大家要求把A送的放在最前面，再三叮嘱一定要写上“未婚妻某某祝贺”，说他父亲会讲A送的花篮摆在店铺最显眼的位置。　　A只有勉强答应，那头的他显得非常高兴这样的风俗我大概也听说过，但是A长了个心眼直接打周伟华留给A的手机号。于是按照他给A的那家花店的电话打了电话，问对方知不道广州番禺南沙区黄阁镇黄阁大道36-38号叫做周伟华的，那个花店人说，认识。A当时还询问对方花店叫做什么名字，座机号码是多少，对方告知了，是.花店人告诉A今天格力空调开张，很多来这里买花，还问这花篮多少钱一个，对方回答档次不一样，价格就不一样，一百到一千的都有，他还强调送的越贵越有面子，而且是送双数，我寻思就挂了电话。A自己琢磨着有点问题，考虑再三，还是上网查查，（因为平时A比较相信度娘，上网查询：花篮。就出来好多这样的骗局信息，看到那句话：凡是没见面就问你借钱的网友均是100%骗子，切记！ 好在没有轻易听信对方的花言巧语）　　后来A发了信息给这个叫做周伟华的，“对不起，在这个特别的日子里发这样的信息，我们还未见面之前还是不要发生经济往来的好，我相信你会有办法解决目前的难题”，结果这个叫做周伟华的回复：“算了，你怕我骗你就别交往了，简直不可理喻！”A看后，庆幸，好在没有上当！　　狐狸尾巴终于露出了了，各位看官，这就是花篮骗局，事后A回复：“真是太遗憾了，等了这么久原来是一个花篮骗局，下次能够玩点新鲜的吗？不过你的演技真的很不错，不做演员真是可惜了”。　　花篮骗子照片　　　　　　　　　　
　　罪魁祸首还是百合网我已经起诉百合网了　　这是判决书，不过我败诉了，我会上诉的
　　我正在经历，还好幸亏在之前先在百度上查过，婚恋网站的骗局案列，他说明天是开店十周年庆典，今天晚上已经让我跟他父亲通过电话，笑死我了，他父亲是在念纸上已经写好的。傻子都能听出来。期待明天他继续表演。
　　陈飞扬
　　花蓝骗
　　@警惕花篮骗局 6楼
18:33　　如果一个男的对你很好，你千万不要激动，这种人会让你先爱上他，然后以公司开业，什么庆典为由，让你送花蓝，这个花蓝一个就要几百几千，所以没有见面前，一定不要钱付出，在乎你的人是不会让你还没有见面就出钱的，请大家小心为好，有一些恨嫁的或者失婚女性，是他们攻击的对象他们演戏很真，所以请大家不要上当了！下面我就将在我身边的真人真事给大家描述出来，给各位恨嫁的女孩子提个醒！　　我一个朋友，暂时喊她叫做A吧！A是一个很重感情的女孩子，但是由于长相不是特别出众，加上即将奔三，一直没有找到男朋友，家里又一直……　　-----------------------------　　
　　-----------------------------Rongr225:我的朋友也遭遇了花篮骗，不过没骗着就被识破了，骗子名叫?伟雄，广东惠州人吧，反正他自己这么说！ 　　
　　中炎投资公司是个什么公司　　
请遵守言论规则，不得违反国家法律法规南京大学白石小百合账号上面显示 错误! 此帐号已经有人使用,请重新选择。! 那是哪一项填写重复了呀？ - 叫阿莫西中心 - 中国网络使得骄傲马戏中心！
南京大学白石小百合账号上面显示 错误! 此帐号已经有人使用,请重新选择。! 那是哪一项填写重复了呀？
的blog阅读
者: [uid]sgsoldog[/uid]
题: [转载] 第一性原理简介
间: Fri Dec 16 15:00:54 2005
【 以下文字转载自 D_Materials 讨论区 】
【 原文由 Tropicstar 所发表 】
myid520 译
X.Gonze,J.M.Beuken,R.Caracas,F.Detraux,
M.Fuchs,G.-M.Rignanese,L.Sindic,M.Verstraete
uniteP.C.P.M.,Universite catholique de Louvain ,Louvain Cla-Neuve,Lelgium
G.Zerah,F.jollet,M.Torrent,A.Roy
Commissariat a l’energic atomique,Bruyeres-le-chatel,France
Mitsubishi chemical Corp.,Yokohama,Japan
PH.Ghosez,J.-Y.Raty
Departement de physique,U,de Liege,Liege,Belgium
Corning Inc. ,Corning NY,USA
用密度泛函理论（DFT）计算电子结构、总能量和材料的其他性能是一个不断发展的领域
。为了站在知识的最前沿，处理恰当的话，一个DFT软件工程能够很好的得益于广泛的合
作。另外，现代软件工程思想可以使软件开发相当容易。ABINIT工程依赖下列信念：资
源的自由性、可靠性、可移植性和在开发复杂的平面波赝势时所强调的个人文档。
在这里，我们介绍在GNU公共软件协议下发布的ABINIT3.0版，给出其功能和一直使用到
现在的各种软件技术：用perl脚本和cpp指令处理处理Fortran90资源文件的独立集合来
为各种操作系统产生并行或串行目标代码；200多个自动测试来保护现有的功能；一直坚
持的严格的编码规则；在线帮助、自学教程、网页资源等文档。
关键词：密度泛函理论，软件工程，电子结构
pacs:71.15.-m,81.05.Zx,31.15.Ew,89.88.+h
在过去的20年中，以量子力学和电磁理论为基础的材料性能的第一性原理计算已经经历
了巨大的进步，而应用广泛的密度泛函理论正处在这个快速发展的领域的核心。他已经
在各种计算机代码和和一般理论（如时间依赖泛函理论）中应用，或者用作很多复杂形
式体系的基础。下面，我们不加综合地列举这一领域80年代中期以来的一些里程碑，来
看一下它的持续发展，并多花一些笔墨在赝势的使用上。
1985年，Car和Parrinello建议了一个把DFT和分子动力学统一起来的算法[4]，同年，H
ybertsen和Louie独自地把Hedin[5]的基态近似应用于单电子的叠加和移动能的计算。1
987年，固体的动力学和介电性能上的线性响应方法被Baroni,Giannozzi 和Testa[7]实
现。Allan和Teter[8]建议用分离Kleinman- Bylander[9]赝势和Car-Parrinello的技术
进行结合，1990年，Vanderbilt设计了一个超软赝势[10]，三年后，和King-Smith一起
指出电子极化是一个几何相[11]；同时，举行串行运算的计算机上的运用开始出现。赝
势方面更加精细的思想也产生了导致了Blochl的投影放大波方法[13]。为了解决激发态
、时间依赖DFT和Bethe-Sapeter方程，格式化[13-18]最近开始实现。
因此，为了和时代同步，用来进行材料到第一性原理计算的计算机程序必须包括越来越
多的功能，随着发展，这些功能变成基本的功能，而在原有的计算机代码上测试一个新
想法也明显变得无效。另外，由于软件越来越多的功能需要合并，一个人很难在不牺牲
自己研究的情况下来继续这一软件的维护和开发。而DFT形式系统及其应用、应用领域的
改进、普遍化需要10年甚至更长的时间。所以，小组开发，甚至国际合作是必要的，在
这种情况下，各种现代软件工程技术证明时非常有用的。
ABINIT软件工程开始于1997年，在这个基础上，作为一个开放的软件工程，它对功能没
有限定，并且使用多种软件工程技术来方便很多开发组之间的国际合作。ABINIT的主程
序使用平面波和赝势，来进行材料性能的密度泛函计算。它的一些特殊的功能还能处理
反应函数。
现在，软件有下面四个特点：1.开放源代码（可以在GNU公共软件协议下使用），2.自我
测试，3.在不同的操作系统下进行并行、串行运算的可移植性，4.个人文档。另外还提
供给使用这一个自学程序。而为了方便国际合作，还制定了旨在要求各个开发组要有统
一明确代码风格的协议。网站[20]提供官方版本、赝势、各种应用程序、基准结果、电
子邮件列表和参考文献信息。
本文将以3.0版为例进行介绍。以前的版本不是在GNU GPL下开发的，网站也将不提供这
和ABINIT软件包所包含的各种应用程序一样，ABINIT主程序也是用Fortran90编写的，它
包括300多个子程序，大约10万行，其中大约三分之一是注释。软件包包括文档文件、自
动测试脚本和赝势文件。
ABINIT v3.0于2000年12月第一次发布，从那时起，程序的错误被定期的报告并在网上公
布（现在的3.06版）另外，软件开发的主流一直超过3.0版，包括现在的3.1和3.2版，本
文的作者对3.0版作出了主要的贡献，另外还有50多人对本版本也作出了贡献，他们改动
了一个或一部分子程序，或发现了一些程序错误。
现在，一些其他代码，如基于平面波和赝势，和基于投影防大波（超软赝势形式的改进
）的代码可以获得，我们提供的列表[21-31]可能不太详尽，他们中的有些是商业软件，
其他的可以通过和主要开发组合作来获得，或者直接从网上下载。
在未来，发展与下列两种团队的合作是有意义的，一种是在GNU GPL下提供他们代码的团
队（这对开发者和使用者都提供了法律的保护），另一种是免费提供源代码的团队。和
组包括相互检测软件精度和速度，以及共享子程序和函数库。显然，在这方面，最根本
的是要开放源代码。
在本文中，我们首先集中介绍ABINIT3.0版的特色功能（第二部分），其中有些是一些基
本功能，已经在相似的代码中使用的很久一段时间，我们在这里不进行太详细地介绍。
其他的是第一次应用的功能，也是我们将详细解释的。当然，全面地描述还将会另外发
表。另外，我们还将粗略的给出软件包的结构和输入文件的组织形式，然后在第三部分
，我们将叙述世界范围内小组开发的组织形式，而现在软件工程的思想（自我测试、脚
本可移植性、个人文档、使用者自学）使这种组织成为可能。在最后一部分，我们将批
判地讨论现在的成就和我们所遵循的开发模型，同时指出需要改进的地方。
2 特色功能
ABINIT软件包包括主程序及应用程序MERGE、IFC和CUT3D。主程序是用来驱动不同密度泛
函的基础计算的。下面我们将区分ABINIT在基态计算、结构相关计算（最优化和分子动
力学）和反映函数计算方面的特色功能。对大部分功能，给出了许多出版我作参考，包
括粗略的描述。然而，有两种情况：自旋轨道影响的处理和绝热联结波动耗散理论，AB
INIT时用艺术态来处理的，同时给出大量的解释。而应用程序MERGE、IFC和CUT3D将在一
个单独的部分描述
2.1 电子基态计算功能
密度泛函计算中的赝势平面波技术是Payne和其同事来审查的。
ABINIT 3.0版支持下列密度泛函近似：Perdew-Zunger[34]和Teter[35]的局域（自旋）
密度近似（LSDA)，Perdew-Buike和Ernzerhof[36]的广义梯度近似(GGA)，Van Leeuwen
-Baerends[37]的GGA时3，以及Gunnarsson-Lundpuist[38],Wigner[39]和Slater(X-alp
ha)[40]的非自旋极化局域密度近似（LDA）。
在可分离的形式下，ABINIT能够使用的赝势有下列各类型：
--在数字格子上的标准正则赝势，例如Troullier-Martins[41]赝势，但是不仅仅这一个
--Goedecker-Teter-Hutter[35]赝势或Hartwigsen-Goedecker-Hutter[42]赝势，这两种
赝势的解析形式特别简洁。
-- Teter[43]推荐的扩展的正则赝势。
网站上给出了LDA赝势的两个完全或者说准完全集，一个是Troullier-Martins类型的，
另一个是Hartwigsen-Goedecker-Hutter类型的。Fritz-Habor-Institute代码可以用来
产生新的赝势，特别是那些用在GGA中的。Hartwigsen-Goedecker-Hutter有一个自旋轨
道的部分，与其相对应，Martins赝势代码[45]的最新版本也能产生自旋轨道依赖赝势。
根据Louie,Froyen和Cohen或Tetet的建议，非线性相关中的核心相关是可以处理的。最
后，需要研究核心水平化学位移的核心洞赝势也是可以使用的。
在每个k点（倒易空间的波矢，通常在第一布里渊区），波函数用平面波有限集的数字系
数来表示，这个系数由动能的截断值决定。k点的集合可按照Monkhorst-Pack[48]系统或
一般形式自动产生。在倒易空间中，由于规则地分隔在不同直线上的k点可以自动的产生
，电子带结构的结果也可以很容易的形成。
对称性用来减少布里渊区取样所需k点数目，所以，只用不可约去的部分才需要取样。
ABINIT 3.0提供不同的可能性来描述倒易空间的波函数（平面波系数）。经常情况下，
波函数的系数是复数，是旋量空间的标量。利用独立的自旋和反自旋波函数，以及自旋
轨道耦合、旋量波函数，ABINIT 3.0能处理共线磁化（铁磁性和反铁磁性）。能处理自
旋轨道耦合是ABINIT的一个想当高级的功能，我们也将在一个单独的部分中介绍。但是
ABINIT 3.0中还没有处理非线性磁化的能力。对具体的波矢在时间反转（和通过倒格子
矢量的转变）下的不变形，如（0 0 0）或（1/2 0 0）……，波函数用适量的平面波系
数的二分集代表。
在实验有效势中，决定波函数的算法是Teter 、Payne 、Allan的连接梯度算法的改进型
算法。电子密度就是通过这些波函数构造，并且用来产生势的Hatree能和交换相关能部
分。这样，输入实验势就和输出势联系了起来。各种算法允许反复地建立向Kohn-Sham势
(输入势和输出势相同的固定点)收敛的实验势，Kohn-Sham势包括：简单混合[50]，And
erson势[51]和基于势的连接梯度算法[52，53]。这种势的预处理通过典型介电函数或近
似介电函数完成。
当电子密度已经自洽地建立，相应的势可以非自洽地用来简并非占有态（需要带结构简
并，或ACFD形式，见2.5部分）。
态依赖占有数可以用多种方式处理，可以手工建立（每个能带同时建立，或每个自旋、
每个k点或能带单独建立），也可以自动加载。对后一种情况，使用者可以选择半导体填
充或金属填充，当然，这要根据不同的抹去方案，看是高斯抹去（Fu和Ho[55]），厄米
―高斯抹去(Methfessel和Paxton[56])，还是冷抹去（Marzari[57]）。
2.2 有关结构方面计算的能力
ABINIT给出了对称规范的不同选择。假如给出了原胞中的原子位置，则对称操作能被自
动识别，相反，如果给出了对称操作和不能约去的原子，原胞中余下的原子也能自动产
空间点群能够根据国际晶体表[58]给出，在这种情况下，对称操作表可以通过一个数据
库初始化。假如晶格常数和对称操作已经给出，ABINIT就能建立布拉菲格子和点对称群
（ABINIT3.0版中还没有点对称群）。
得益于Hellmann-Feynman理论和密度泛函理论中的应力理论，以和应力在所有情况下都
可以计算，除旋转轨道情况外，总能量在这种情况下也都能计算出来，也就是总能量的
计算适合标量波函数和自选极化情况、金属和绝缘体材料，以及各种交换相关函数。
使用Broyden算法[61]、改进型Broyden算法[62]、当离子的运动方向和力的方向相反时
运算立即停止的Verlet算法[63],力和应力可以用来最优化结构（最小化力和应力，任意
限制它们的分量）；使用Verlet算法[64]或Numerov算法[65]，力和应力还可以用来产生
分子动力学轨道。Verlet算法情况下，可以使用Nose-Hoover 和Langevin热浴[66-70]来
取样正则系综。
另外，代码能自动分析键长键角，以及xyz形式的原子坐标。
2.3 相应函数方面计算的能力
我们首先考虑对原子位移和静态同类电子场的相应。
在Berry相形式下，可以计算出介电极化，这一特色功能对绝缘体、磁体和非磁体可用，
但是，自旋轨道分裂情况出现时还不可用。
线性响应技术（密度泛函微扰理论）[7,71-74]可以用来计算对原子位移和同类电子场的
响应。计算将产生在选择的波矢处的动力学矩阵、Born有效电荷、或者（电子）介电常
数。这些量放在一个数据库中，可以用MERGE和IFC来分析（后面还要专门讲到）。对AB
INIT3.0版来说，线性响应在自旋极化情况、旋量波函数情况和对交换相关函数的GGA近
似情况还不可用。
磁化矩阵和介电矩阵（在频率为0，和对倒格矢）可以通过对整个态公式求和[77，78]计
电子激发也可以在时间依赖密度泛函理论下研究。
2.4 自旋轨道耦合方面计算的能力
自旋轨道耦合在平面赝势代码中不是经常用到，在力、应力、响应函数可以利用的ABIN
IT3.0版，它的执行将详细描述[79]。我们在这里给出一个简略的描述，描述他在计算总
能量时的执行情况。从狄拉克方程开始，能够得到和薛定谔方程相类似的方程，这个方
程和相对论波函数的最大最小分量相混合，具有1/c 的数量级。对每个l+ 和l- 轨道，
赝势都能产生，对赝势进行适量的线性组合，将产生标量相对论部分 和自选轨道部分
。按照参考文献[42,81],重新组合赝势，写出电子--离子势，离子势是标量相对论相(S
R)和自旋轨道相（SO）的和：Vei= + (1)
标准波函数基扩展到旋量波函数基，旋量波函数基的元素由 表示的二分量波函数。一般
矩阵元的自旋独立部分为:
自旋轨道部分为:
在Kleinman-Bylander公式中，自旋独立相的结果就是著名的： =4π（2l+1）fl(|G|)f
l(|G’|)Pl( (4),其中 、fl(G)是倒易空间中Kleinman-Bylander形式的因子，Pl是勒让
得多项试。这个结果，从增加的定理： （5）
到得到自旋轨道项的相似算符，我们首先考虑矢量： 。利用定义L=r p,恒等式 和 ，我
们得到： = -i 。因此，我们把矩阵元素变换成：
最后一项的计算模仿自旋独立项的计算。注意，由于与 叉乘， 项无贡献。最后，我们
其中 是 的一阶微分。
通过这个公式，我们能够求出体系的总能量，同时，通过对原子位移和原胞变形微分，
我们能求出力和应力。
2.5 绝热连接波动耗散理论
DFT中的局域密度近似和广义梯度近似经常产生固体、表面、分子物理和化学性能的有效
精确描述。特别是它们能给出原子结构的实际的计算，以及弹性和震动性质。虽然在分
子相互作用和相关势能面方面，GGA比LDA有改进，这也是GGA的典型部分，但是，我们还
是需要更精确的函数来克服临界时的缺陷，例如，GGA还不能在化学精度下预测化学反应
热力学（分子分裂能、反应热和激活能垒），和LDA一样，不能恰当的包括远距离子系统
间的范德瓦尔斯相互作用。
绝热连接波动耗散理论（adiabatic-connection fluctuation-dissipation theorem A
CFDF）允许产生完整的非局域交换相关函数，这个函数在远距离情况下包括范德瓦尔斯
相互作用或精确的非局域相关，因此，这一理论超越了传统的局域密度近似和广义修正
近似。并明确给出了密度为 的电子系统的交换相关能[82]：
其中 （虚频率）是通过一个标量库伦势 产生电子相互作用的系统的动力学密度响应函
数，它进入一个修改的外部势，以便对物理（ ）基态密度保持不变。对 ，我们处理无
相互作用的Kohn-Sham系统,它的响应函数由Kohn-Sham本征态 和本征值 ，即：
其中求和包括所有占有态（ =1）和非占有态（ =0）。当 时，相互作用和Kohn-Sham响
应函数通过代森型屏幕方程相关，代森型屏幕方程为：
= + （10），
其中 = 库伦和交换相关核函数，它在时间依赖密度泛函理论[83]中建立。就像在参考文
献[84]中提到的一样，交换相关能的相关部分和交换部分可以分开。方程8-10可以归结
为ACFDT形式。原则上，它能产生精确的密度泛函交换相关能，但实际上，通过对未知的
时间依赖交换相关核函数、Kohn-Sham响应函数、和（或）代森方程的解使用具体的近似
，从各种近似的动力学密度响应函数都可以产生近似函数。
在ABINIT3.0版中，初始的Kohn-Sham基态（也就是 ）是在LDA或GGA条件下计算的，而把
响应函数和核函数看作 点，并在定义的平面波代表中处理，例如，对 ，通过 = 来处理
，其中G是倒格矢。通过对方程（9）中的所有态求和，得到Kohn-Sham响应函数。使用不
同的时间依赖交换相关核函数，我们不需要进一步近似就解出代森方程（10），因为线
性方程系统如下：
为得到交换相关能，我们给相关能付值为：
因此，增加了交换能泛函[86-88]。至于对 和u积分，我们重复解方程（9）和（11），
并使用高斯求积分公式。
功能的实现还处于试验阶段，还需要最优化。我们的初始目标是评定ACFDT家族的各种函
数的精度。
2.6 并行计算能力
我们已经实现不同水平的并行计算:把不同的k点、和给定k点中不同的态相关的工作、与
各波函数的系数相关的工作分配在不同的处理器上。当然，还有很大的空间来最优化各
水平的并行计算，所以，我们这一部分的目的不是给出加速多少的详细分析，只是给出
并行计算下的基本思想。
最有效的并行计算（和通讯有关的大量计算）基于k点的分配，它利用MPI库实现，对两
个以上k点的情况（电子结构、确定总能量、响应函数）都适用。它工作在大型并行或S
MP机上，或者网络基通讯的机群上。不幸的是，k点数目和原胞尺寸的比值随着原子数目
的增加而减小。尽管如此，还有一些情况是工作量和k点数目都比较大（例如金属、响应
函数下小单元的性质）。下面我们举几个例子，在有25个处理器的SGI Origin 2000（处
理器是R12k,300MHz）上,可以得到比一般情况快20倍的速度，而在有25个Intel奔腾III
处理器（550MHz），操作系统为Linux的机群上，速度依然是一般运算的15倍以上。
在给定k点内各态相关的工作的扩展中，ABINIT3.0版仅仅实现了在响应函数情况下的扩
展。正如k点的分配，它利用MPI库实现，工作在大型并行或SMP机上，以及网络基通讯的
机群上。他的缩放比例和系统尺寸的大小是线性关系，并且和k点的反缩放比例反平衡：
k点数目和态数目的成绩大约是个常数，对相似的系统用相似的精度处理，它的典型的值
在200和1000多之间。对响应函数计算来说，在处理器间通讯的数量方面，态基并行运算
只比k点的分配大一点点，然而，ABINIT3.0版中，每个处理器所需的内存不因处理器的
增加而减小。
最后，使用OpenMP编译指令，可以在SMP机的处理器上分配和波函数的系数有关的工作。
利用这一技术，FFT算法和对波函数的非局域算符的应用已经可以并行运算，然而，这个
应用是三个并行运算中最新的，效率还不是太高，一个典型的速度是在四个处理器上并
行计算时，速度是在一个处理器上运算的2倍。有意思的是，MPIO和penMP并行计算可以
同时使用。
2.7 应用程序
应有程序MERGE和IFC时用来分析数据库的，这些数据库是ABINIT程序产生的总能量对原
子位移的微分数据库和ABINIT程序产生的同类电子场数据库。数据库还包括动力学矩阵
、波恩有效电荷以及电子介电张量。MERGE时用来进行数据库处理的，而IFC时进行数据
库分析的。
根据参考文献[74]，IFC分析数据的功能如下：
（1） 基于波恩有效电荷张量和电子介电张量，建立原子间力常数，包括其渐进行为。
（2） 通过数据库提供的动力学矩阵的傅立叶插值，在布里渊区中的任意点建立动力学
矩阵，从而建立相应的本征矢和本征值，其中本征值形成声子带结构。
（3） 建立在零波矢时绳子的对称特征。
（4） 在准谐波近似的情况下，通过在布里渊区中对声子自由度积分，用波松-爱因斯坦
占有因子建立热力学性质（如自由能、热容和熵）。
（5） 频率低于电子隙的情况下，建立频率依赖介电张量。
应用程序CUT3D用来分析ABINIT程序产生的实空间三位密度或势文件。特别是在一个与三
维原胞相交的二位平面上顺着某一直线或在任意一个点，它都能添加已知的密度或势到
一个三维任意地卡尔体积中去。另外，它还能使原来非格式化的的密度或势文件格式化
，以便作为图形文件MATLAB[75]和MOLEKEL[76]的输入文件。
3 软件开发的组织
就像在前言中提到的一样，ABINIT软件工程依靠一大批开发者，他们属于不同的国际开
发组。为了协调这种软件工的开发，计算机科学家开发了软件工程技术，其中有些软件
工程技术包括辅助软件，一般在GNU一般公共协议些可以得到，另外一些软件工程技术是
引导开发者工作的一些思想。在这些软件工程技术中，有一部分已经在ABINIT中采用，
我们将在这一部分叙述，另外一些还没有使用，这将在讨论部分中叙述。
3.1 软件共享：GNU一般公共协议
GNU组织提出了重要的自由软件思想：软件的版权属于开发者，但是软件在一个协议下发
布，这个协议保证用户有权利得到源代码、更改源代码、甚至有权利重新发布软件。这
一思想在一个法律文件中进行了详尽的阐述，这一文件可以在GNU网站[19]得到。所谓的
GNU一般公共协议已在现在所能得到的绝大部分自由软件中使用。
我们已经选择在这个协议下发布ABINIT,其软件包中源文件和文档文件的版权属于其作者
和发布协议。通过这种方法，开发者保留他们努力得到的知识，同时允许别人在未来的
版本中改进他们的工作，用户也有权得到源代码，并且可以在调试方面做出自己的贡献
3.2 自动测试
开发者修改工作可能引入错误，这个在代码开发中的问题源在一个大的开发作加入时变
得更为严重：大部分开发者仅仅知道代码中有限的部分，同时开发者小心的水平也是千
差万别。为了避免这一缺点，引入了自我测试。
在一个完整的ABINIT软件包中包括200多个测试，每个测试的测试环境包括一个输入文件
、赝势以及为运行测试批和自动分析他们的结果所需的脚本。对ABINIT的每一个功能都
存在一个或多个测试，这些测试都在软件实现时建立。在个人电脑（CPU为800MHz的奔腾
III）上运行时，每个测试一般情况运行几十秒。
自动测试环境的最终目的是当一切正常时，程序返回一个正值（一个字节），否则返回
适当的失败分析（当然超过一个字节）。在ABINIT3.0版中测试是多水平的：
（1） 第一个脚本叫fldiff(浮点差别)，它写在PERL中，它用来计算结果文件和参考文
件的差别， 是用这种方法，由于不同平台和计算机的果有差别引起的浮点错误在一定范
围内将会被忽略。这个工具要求结果文件和参考文件要十分相似，也就是只存在浮点方
面的差别。测试成功的话，只有几行文字写到概要文件中去，集中了自动测试批的测试
（2） 当有不可忽略的差别时，unix命令diff的结果将会自动得到，它允许对更复杂复
杂差别进行简单测试。
（3） 每次ABINIT运行，都有一个详细的log文件，其中列出了运行的错误、警告、注释
当要是他们的开发生效时，每个对ABINIT工作的贡献者都要求建立一个测试，修改的原
文件、测试用到的输入文件和参考文件都将在包含在下一个版本中。开发者把参考文件
送给负责官方版本的人也是必要的，这样能够确保软件的新功能确实是安全的。
最后，还有一小部分测试不是用来测试软件功能的，而是用来标定软件速度的，他们比
较不同版本的需要大量CPU的程序, 并允许用户比较不同机器的速度。
3.3 可移植性
ABINIT3.0版安装在了一下平台：
（1） 使用Linux操作系统，编译器（波兰特集团公司）或富士通编译器，奔腾、奔腾I
I、奔腾III处理器的个人电脑。
（2） 使用Windows 98或NT操作系统，PGI工作站程序组，英特尔486或奔腾II处理器的
个人电脑。
（3） 基于阿尔法处理器（EV56,EV6或EV678），使用操作系统固件的康柏计算机。
（4） 基于阿尔法处理器（EV56），使用Linux操作系统的康柏计算机。
（5） 基于Power2和Power3+处理器（型号590，3终端，夜鹰级）的IBM RS6000计算机。
（6） 基于R1000处理器的SGI Origin2000计算机。
（7） CRAY T3E 计算机。
（8） 富士通VPP700计算机。
（9） 太阳公司ultraspare II计算机。
（10） NEC SX4和SX5计算机。
（11） 日立SR8000计算机.
（12） 苹果公司的麦金托什机。
对平台（1）到（7）来说，主代码的最优化非常高级，对其它平台来说，主代码的最优
化可能更好。开发者开发了用于大部分平台的二进制代码，这些都能在ABINIT网站上得
由于cpp指令，以及MAKE文件和各种脚本的使用，ABINIT已经可以成功的安装在很多平台
上。ABINIT援文件的唯一集是在编译时匆忙的预处理的，它产生时间依赖代码，结果，
ABINIT可以在UNIX型、Windows型以及Macos型的操作系统下工作。
对ABINIT安装的每一个平台，一个包括编译器名、选项和为预处理而选择的信息的机器
依赖文件必须准备（大约15行）。这个机器依赖文件对ABINIT的版本依赖很小，一旦这
个文件安装在一个平台上，进一步安装在相同的平台上就十分相似了。
相同的软件技术用来维持ABINIT串行和并行版，他们都通过源文件特别集的预处理产生
ABINIT MPI版已经在SMP结点的集群上测试和工作，机群由上面提到的计算机组成（太阳
公司ultraspare II计算机除外），工作环境为UNIX型操作系统。ABINIT的OpenMP版已经
在下列SMP计算机上测试，这些计算机包括：SGI Origin2000计算机、44个Power3+处理
器的IBM RS6000计算机（每个节点四个处理器）、以及康柏/数据设备公司的ES40 EV67
计算机（每个节点四个处理器）。
3.4 自动测试的可移植性
可移植性和自动测试的相互影响带来很有意思的问题。就像在3.2部分提到的一样，结果
的分析包括和参考文件的比较，参考文件时曾旭在一个特定的平台上运行产生的，虽然
fldiff能够把浮点运算结果放在一个能忍受的水平（如果来个浮点数的差小于10e-12,它
就认为这两个数相等），这已经足够消除大部分因平台改变带来的误差，但是在设计测
试输出结果方面，还需要额外的注意是，否则，大部分测试情况的结果都会包含大量的
平台以来部分，这很不适合自动测试。
产生大量的数字差别的第一种情况是由于矩阵的对角化。本征值是简并的，则相应的本
征矢可以在得到的矢量字空间中任意取值，即使是非简并的情况，波矢的位相也是没定
义的。大部分的计算量（）不依赖矢量和位相的特定选择，然而，不是总是这样。确实
，下面是我们所感兴趣的：
（1） 通过动力学矩阵对角化产生的输出结果和原子本征位移的可视化。
（2） 从不同的空间几何，通常是有更高对称性的空间几何产生的波函数来初始化电子
波函数（为了能很好的收敛）。
这两种情况，由于相和矢量相等的线性组合的选择，排除任何合理的简单自动分析，输
出文件的具体部分可能完全不同。
ABINIT3.0版用两个方法解决这一问题，第一（程序fxphas.f）通过最大化他们的平方和
，来固定本征矢的位相，然后选择第一个非零元素的符号为正；第二（程序phfeq3.f和
mkkin.f），通过在每1.0e12部分区一个部分的水平，来稍微改变动力学和动能矩阵的对
移植问题的第二种情况由算法本身的不稳定产生。例如，基于Broden算法的原子间距离
的最优化包括海赛形式的矩阵的时间反演，如果从一个分子的两个原子间开始最优化，
而且在这个区域中，原子间势的二阶微分几乎为零的话，很小的平台依赖差别将会被放
大，以至浮点结果不能满足所需的要求。
在ABINIT3.0版中，这个问题是通过下面方法解决的：首先尽量避免在这些区域中进行测
试，如果确实不能避免，就要让这一算法在一两个循环之后，在结果在不同的机器上还
没发散之前就停下来。
移植问题的第三种情况是由于实数的分类问题产生的。假设在一个大部分数不相等的列
表中，有两个数是相等的，如果这种情况不小心的处理，分类算法将会对他们进行对称
分类，但是这将因平台的不同而改变。在ABINIT3.0版中，遇到的这些问题的物理实现如
（1） 以原子间距离为顺序，产生相应原子的近邻列表（对称性经常在相同的距离处给
出原子），如程序bonds.f 和rsiaf9.f。
（2） 以k点的量级为顺序，产生相应k点的列表（为了有效的填满不可简约的布里渊区
），如程序listk.f。
（3） 为了打印，找到密度的最大最小点，如程序prtrhomxmn.f。
为了纠正这一问题，采取的措施是在一定差别允许范围内，让算法知道这两个数是相等
的，同时保留着两个数字原有列表（分类前的列表）的顺序。
可移植性中的第四个问题和物理意义没有直接的关系，假如一些浮点数必须打印，如果
截断精度小于自动分析工具所允许的值（如使用f8.2的FORTRAN输出格式），由于数字和
轴离的太近，舍入情况在不同的机器上也不同。
要解决这一问题，必须在打印数字之前，使用应用机器独立的截断函数，然而这还没在
ABINIT3.0版中应用。
即使这四个问题都解决了，还有一些测试情况，它的结果文件和参考文件中的浮点数差
别超过所选的容许范围，这些测试情况通常包括响应函数计算，这是由于初始的基态计
算有点缺少可移植性，而这少量可移植性的缺乏又在随后的响应函数计算中被放大，这
样在基态计算中可以忽略的部分由于在反应函数中被放大而不能忽略。得到了缺乏可移
植性的观测结果不等于程序中有错误，然而，有关输出文件中不同部分的意义的知识还
是必需的，以此来验证程序是否出现了反常情况。
尽管还存在缺乏可移植性的部分，但是使用自动测试工具，对测试的完全集进行的概要
的检查也就至多需要2到3分钟。
3.5 个人文档
当人们说起大型科学的代码（如ABINIT）时，个人文档的软件思想就经常被提到。我们
应该注意，为了使开发者不需对软件包有详尽的了解就能维持这些程序，和一个程序有
关的文档应该定位在这个程序，或者通过一个链接定位到这个程序。这一思想在“有文
化的程序设计”中实现[83]。个人文档工具允许开发者不需保持两个独立的文档就能够
利用源代码编写程序文档。用户能用自动的方法掌握代码的整个结构，所有文档都有本
地文档产生。在一些局部，必须使用一些标签来引导自动工具。
自动工具能产生有超级链接的HTML文件，例如每个程序的页眉都显示在屏幕上，包括他
功能的描述、调用他的程序(父程序)的名字以及他所调用的程序（子程序）的名字。所
有程序的目录和搜索工具都能够建立。很多其他功能，例如数学方程的自动格式化也能
够建立。为了系统地实现上述功能，开发者必须为他们的注释使用并坚持使用一个标准
的格式。由于这经常被看作是额外的负担，在这方面所作的工作应尽量少。
在ABINIT3.0版中，各种辅助程序用来实现这一思想。
这些程序中的第一个程序是ROBODoc[90]，它是个文档工具，用来从源文件中提取特别格
式的的注释页眉。他自动地从原文本中为每个程序产生HTML文档文件，还包括连接其它
程序的超级链接。ROBODoc能在很多语言环境下工作，包括Fortran90。ABINIT的开发者
必须为每一个函数建立标准的页眉，包括包括函数或称序所有的信息，这样ROBODoc为变
量、类程、函数等建立索引表。
然而，ROBODoc不提供父子程序的自动列表，虽然如果提供了这个的列表，把它们连接起
来的超级链接将由ROBODoc自动建立。一个叫“父页面”的PERL脚本已经设计，它可以用
来在ROBODoc处理之前自动地在每个程序之前插入父程序。
数学方程是ROBODoc遇的的另外一个问题，这个问题主要是通过用免费软件scr2tex[91]
来解决的，这个软件能从子程序的原文件中产生有很好格式方程的LATEX文件。每个程序
中插入LATEX风格的方程格式是ABINIT格式的一部分。在scr2tex中使用的格式可以同时
在ROBODoc中使用。
在ABINIT3.0版中，源代码是从ROBODoc和scr2tex出发来确定格式的，以便HTML和LATEX
文件可以分别作为在线和离线文档。由于例如TtH的转变，在未来的版本中，将写一些像
HTML文件一样易于获得的LATEX格式的文档。
3.6 用户自学
ABINIT3.0版中另一个有意义的方面涉及到不需要联系开发者或到学校去学习就可以学会
怎么使用他们。为了这个目的，开发者建立了自学教程，教程是HTML格式的，共有5课，
每课2小时（包括在一个CPU为800MHz的奔腾III的个人电脑上计算的时间），用来引导学
其中教程的前两课是为了使学生熟悉是用各种算法和交换相关函数，来计算氢原子的格
式化的总能量、电子结构和最优化原子间距离。第三课通过一个绝缘体的范例―硅的电
子结构的计算，来介绍固态思想。第四课集中介绍金属铝以及其表面能的计算。第五课
介绍绝缘体AlAs的动力学性质和介电性质，包括伽玛点的声子、介电常数、波恩有效电
荷、LO-TO分裂以及整个布里渊区中的声子。
在线手册和帮助文件的有效性，以及自动测试中输入文件的有效性使ABINIT基本使用方
法的介绍更完备，在线手册和帮助文件包括了所有输入变量的详细描述，而自动测试中
的输入文件给出了ABINIT所有功能具体适用的例子。
我们已经尽量简化了ABINIT的使用，其中输入文件的组织依赖于两个文件，第一个文件
是“files”型的软件，用户在这个文件中给出了输入文件和输出文件，以及所使用的赝
势。这个文件很短，一般6到7行。在另外一个输入文件中，用户具体给出了那些默认值
将不考虑的所有变量的值，这个文件由ABINIT解释执行，解释执行前ABINIT检查其中使
用的输入变量的关键字。ABINIT中有150多个输入变量，但是在这个文件中要记住的不超
过20个。元字符允许个输入变量集在一个单独的文件中具体指出，甚至为了基于以前计
算的结果初始化计算，而允许在同一个计算中把各个输入变量链接起来。
3.7合作协议
自由软件、自我测试、个人文档和用户自学的思想的思想对国际合作开发来说重要性是
明显的。组织同步工作也是值得的。
在单独开发组的水平，我们鼓励使用软件CVS[94]来开发ABINIT，然而，所有开发者都必
须能得到CVS库，因为这个软件同步开发工作。另为，开发者之间的协商也是需要的。
SoureForge[95]网站允许在一个很大的范围内交换CVS思想，在不同的组，来自不同的国
家的开发者能够得到相同的库。我们从2000年13月开始试者使用SoureForge库，但是，
对ABINIT规模的软件，通向它的速度太慢了。因此，我们从1997年开始一直坚持的合作
开发协议还一直有效。
我们工作在一个反复的分散―集中模型中：反复的第一步是官方本版在ABINIT网站上发
布，一段时间后，来自不同组的非交迭、持续开发工程列表将建立并公布；每个开发组
在ABINIT已经暂时加上去的部分程序上共作2到3个月，然后把各个开发组的贡献（资源
的改变部分和测试情况）汇集到一块，再用两三个周把各个贡献融合到一块，在新的官
方版本发布前，软件还要至少在5个平台上测试和连接。
最后，我们谈一下ABINIT的代码风格。为了使代码风格相同，我们在“代码规则”文档
中明确地写出了每个开发组都要遵守的规则。这个文件在软件包中可以找到，并且经常
更新。它的建了基于在过去几年来的开发中获得经验（包括在不同平台间移植的经验）
，以及为了避免程序出现各种错误。例如在所有FORTRAN程序中，要求使用“含蓄的不”
“代码规则”包括：需要遵守的变量说明风格、变量名的选择风格、FORTRAN援文件格式
的风格（特别地，格式必须能够用ROBODoc处理）、数据流控制的结构、数组的使用风格
、一些好的代码实例、具体的例外处理方面的代码规则、一些旧的应该避免的代码风格
的实例、BLAS和LAPACK子程序的使用、以及当今反射出来的一些主题。
和ABINIT工程的特定相比较，很多还要在3.0版的基础上去做。我们现在首先描述现在已
经实现的功能，然后描述还缺少的功能，最后叙述现在已经部分使用或还没使用的软件
工程思想。现在，对本软件做出贡献的开发组已不限于作者所列出的。可以通过网站[2
0]，来更详细地了解本软件共输出的状态。
下面是和雄心勃勃的软件工程相关的部分：
（1） Blochl[13]的投影放大波技术和超软赝势[10]的使用。
（2） 包含基态代码的界面ABINIT.
（3） 非线性磁化情况下旋量波函数（包括自旋轨道）的处理。
（4） 现在，在所用情况下，响应函数计算的结果都可以为基态计算所利用。
（5） 手工处理并行运算，以及现在实现的并行运算的最优化，特别是不同水平的并行
运算的同时运用。
下面是我们的开发工作所涉及到的部分：
（1） 实空间群数目的识别。
（2） 密度泛函微扰理论中的2n+1定则[71]。
（3） Kohn-sham精确交换。
（4） Shubnikov(反铁磁)对陈群。
（5） 表面偶极子的处理。
（6） 频率依赖传导。
（7） 和极化Berry相计算有关的高级特色功能[97]。
（8） 晶格Wannier函数的自动产生。
下面是需要做但是还没给出具体做的时间的部分：
（1） 非局域算符的实空间处理[97]。
（2） 态密度在原子轨道分量上的分解。
（3） 对弹性和炼金术微扰得相应。
（4） 几何形状最优化的新方法和转变态的计算。
在软件工程水平下，首先是持续的努力来处理个人文档中的方程。然后得到应该使用的
数据类型（在ABINIT3.0版中很少出现）。FORTRAN90面向对象的特色功能也没在ABINIT
3.0版中使用：ORTRAN90面向对象的特色功能虽然没有C++的发达，但是也应该应用。F9
00面向对象的特色功能方面的局限也必定带来现在的ABINIT基本操作的局限：由于历史
的原因（已存在代码的再利用），我们建立的是F90代码。作为另一个可选的语言C++，
现在很多开放源代码的工程同时使用Python脚本语言和其它语言来计算任务（一般是C语
言，但是为什么不是F90呢）。我们仍然这样认为，我们工程的基本选项应该使它在很多
年中一直是可行的和成功的。
我们感觉ABINIT工程是成功的，从工程的开始算起，我们在相对较短的时间能建立了大
量的功能。合作的组织使很多人都可以对这个软件工程做出贡献。我们实际上一直根据
计算机科学通讯中Linux的开发[98]的例子来开发本软件。
和本文不集中在实际地使用ABINIT来得到有用的计算结果对比，本工程的目的是使用从
头计算的方法为与预测材料性能建立可靠的、艺术态的计算程序。在这一思想的指导下
，网站[20]上列出了70多篇已发表的文献，来显示其广泛的应有范围，以及迄今为止已
使用本代码的了个人和团体。有用输出量的增加是我们成功的基本量度
我们首先感谢ABINTI的顾问委员会，包括本文的作者和不是本文的作者的：N.Hill,K.R
abe,和D.Teter.我们还要感谢对本软件工程、ABINIT3.0版及以后的版本做出贡献者以及
给我们提出建议和热情的支持的人们,包括：A.Alouges,S.Bernard,J.-C.Charlier,T.D
eutsh,S.Erwin,A.Filipetti,R.Godby,S.Goedecker,A.Horsfield,A.Khein,D.D.Klug,C
.Lee,Z.Levine,B.Magne,H.Matthis,J.-P.Minet,V.Olevano,G,IbudamA.Quandt,Na.Sai
,M.Stengel,J.Sullivan,P.Tepesch,P.Thibaudeau,J.Tse,D.Vanderbilt,B.Van Renter
ghem , M.Veithen,A.Zaitsev.We(X.G.,J.-Y.R.,和G.-M.R.),感谢国家科学研究基金（
比利时）在财政和计算资源方面（FRFC工程，NB2.4556.99）的支持.我们还要感谢“po
le d’attraction interuniversitaire”(阶段五，比利时联邦政府)和“Action de r
echerche concerte-interaction electron-vibration dans les nanostructures”的
财政支持。
[1] P.Hohenberg and W.Kohn, Phys.Rev.136,B864(1964)
[2] W.Kohn and L.J.Sham, Phys.Rev.140,A)
为什么,爱她爱得如此深沉,
我已无法割舍的热恋,
连灵魂都已熔化其间,
爱你而为你付出一切,
含着泪,带着微笑!
※ 来源:．南京大学小百合站 http://bbs. [FROM: 172.16.27.182]
※ 来源:．南京大学小百合站 http://bbs. [FROM: 202.119.45.252]
(现有评论0条)
说的太好了，我顶！
Copyright & 2014
Corporation, All Rights Reserved
Processed in 0.0256 second(s), 3 db_queries,
0 rpc_queries