关于本博
我爱搜集网是一个网站,但更准确的说是一个博客。本博博主 Michael Duan.
博主简介
2006年-2010年:华中科技大学材料科学与工程学院,工学学士;
2010年-2015年:华中科技大学材料科学与工程学院,工学博士,研究方向为计算材料学。
现任教于武汉某高校。
关于博客
本博客主要是博主博士期间写下的一些科研经验,总结了自己在科研中遇到问题的解决方案,以备再用或他用。
自知个人水平相当有限[绝非谦虚],而且写博客时也未详尽考究,所以其中难免错漏之处。不过非常欢迎你能够指出,以便我进行更正。
已有 49 条评论
博主 能给我发下morse的势函数吗?我需要Ti-c,Al-c的作用力,我邮箱244294726@qq.com
抱歉,我也没有。请自行网络搜索。
博主您好,在哪里可以看到您的有关MEAM势的文章呢,麻烦您能不能发我邮箱呢513861106@qq.com
中国知网搜索:晶格反演修正型嵌入原子势
谢谢您!我是一个刚入门的研一学生,现在想自己编程实现渗碳体的力学研究。MEAM势中提到提到屏蔽函数如果给定Cmin和Cmax,则S值是定值。看文献有直接给出了铁的S值,但是目前没有找到碳的S值,我是要继续找还是编程序去计算?
还有您的文章中式2.46中有求和公式,其是从1开始求和的,那在哪结束呢?加到无穷大?
如果你使用2NN MEAM,Fe的Cmin和Cmax应该是有的。至于您提到我论文中的求和,是通过截断半径控制的。
博主您好,我看了您的晶格反演法拟合MEAM势Fe的势函数,我对这方面很有兴趣,我个人在尝试自己编程实现MD的计算,不知道您是否方便能将您的一部分拟合代码给我参考一下。
我拟合主要用的PSO算法,网站有提供下载。
http://www.52souji.net/matlab-code-of-particle-swarm-optimization.html
感谢您的回复,我正在尝试自己编写MEAM的程序,我看了您的毕业论文,其中有个 “参考结构下中每个原子的背景电子密度”由于我是初学者,基础薄弱,不能很好地理解这个概念,不知道您知道这个是什么意思么,我看有些文献好像直接把这个当作1。又或者是第一近邻原子数。
您好,我也有同样的疑问,不知道您是怎么处理rho0的,其计算公式和rhoi有什么不同?
在meam理论中,有一个“参考结构(reference structure)”的概念,是为了帮助计算对势(pair potential)。参考结构中每个原子都是在晶格格点上,所以每个原子的rho0都是一样的,其计算方法与rhoi类似。“在式(2.27)中,rhoi/rho0,即任意体系中某一原子的背景电子密度与参考结构中某一原子的背景电子密度的比值,大致可以反映出任意体系中该原子附近的原子结构与参考结构的偏离情况。”
感谢您的回复,对我很有帮助,在我自己编程计算X.W.Zhou的EAM势之后就想尝试MEAM势的计算,发现MEAM势的计算更为复杂。另外我利用X.W.Zhou的EAM势计算其受力与lammps使用相同势函数计算出来的力不能完全对应,我认为是lammps为了简化计算而列出Ndr,dr,Nrho,drho的分段数目不足,而我的是完全计算的力,不知道我的理解是否有误?
这个我不确定。以我的了解,X.W.Zhou的很多势开发都是基于lammps的,应该不会出现不兼容的问题。我建议你再细致的看看是不是哪里出错了。
好的,我再仔细看看lammps的manual和我自己的程序,虽然我坚信自己的程序是没问题的但是还是希望能够与lammps上的力对应上。再次感谢您的热心回复!
我不太明白你的疑问。
就是MEAM势能中的潜入能计算时,您的论文的21页公式(2.27)的rho0的物理意义我不是特别清楚,不知道您知道是什么意思么?
楼主,你好。请问你用lammps建过孪晶模型吗?我尝试了很久,还是没有做出来,希望楼主解惑。
最好用materials studio建模。即便如此,也不太容易建。
呵呵,博主不错,有心了!钱不多,相比于你的付出,真是寥寥!
谢谢支持!