24小时热门版块排行榜    

查看: 604  |  回复: 7
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

lqllinqiaoli

至尊木虫 (著名写手)

[交流] 【求助】流体的分子动力学模拟

为什么对于惰性体系的的流体,动力学模拟时添加的是L-J势,而金属的熔体要用EAM(Embedded-Atom Method)模拟呢?

[ Last edited by lei0736 on 2009-11-24 at 19:58 ]
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

superdirac

木虫 (正式写手)

噢!
  我想问问,用L-J势来描述体系的势能,  用分子动力学的方法,来计算真空中100个

原子左右的团簇,这些团簇应该炸开吗??
  
  我看有人用EAM计算的时候没有炸开,但是我用LJ势算的时候就炸开了!
  不知道是我的方法有问题,还时体系就应该炸开!!
我认为,酒一口一口喝,路一步一步走~步子迈大了,喀~容易扯着蛋
4楼2009-03-08 11:44:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

zyy816

至尊木虫 (著名写手)

★ ★
yuhuobuku(金币+2,VIP+0):欢迎参加讨论 2-22 20:16
这些势能函数都是经过实验数据或者第一性原理拟合出来的,有一定的经验性与适用范围,如果你发现L-J模拟结果能与实验结果符合很好,一样可以用到金属体系的。
上面是一个经验的结果,没有非常强的理论根据必须要这么干
2楼2009-02-22 19:58:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)

★ ★
yuhuobuku(金币+2,VIP+0):欢迎参加讨论 2-22 21:04
LJ势是对势,对势总是预言金属的弹性常数存在Cauchy关系,即C12=C44,这对于过渡金属是不成立的。对它若做些简单修正,如添加一个与体积变化相关的项,可以解决这个困难,但又带来了一些新的困难。嵌入原子法则顺利地解决了这个问题,所以嵌入原子法是当前用于金属,尤其是过渡金属的最流行的方法。但也不是唯一的方法。同时,LJ势及其修正形式也不是说完全不可以用于金属,全看你要模拟什么性质。
3楼2009-02-22 21:03:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuefei06

木虫 (正式写手)

★ ★
jghe(金币+2,VIP+0):感谢你参与讨论! 3-8 16:47
EAM计算没有炸开,而lJ势炸开了,这个问题可能有很多可能,
   1.LJ的势参数的选择是不是有问题,降低排斥,增加吸引势看看如何
   2.体系控制参数温度,步长,截断半径是不是有问题
5楼2009-03-08 13:45:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见