24小时热门版块排行榜    

查看: 601  |  回复: 7
当前主题已经存档。

lqllinqiaoli

至尊木虫 (著名写手)

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

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

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

» 猜你喜欢

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

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

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的回帖

superdirac

木虫 (正式写手)

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

原子左右的团簇,这些团簇应该炸开吗??
  
  我看有人用EAM计算的时候没有炸开,但是我用LJ势算的时候就炸开了!
  不知道是我的方法有问题,还时体系就应该炸开!!
我认为,酒一口一口喝,路一步一步走~步子迈大了,喀~容易扯着蛋
4楼2009-03-08 11:44: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的回帖

superdirac

木虫 (正式写手)

★ ★
yuhuobuku(金币+2):欢迎参加讨论
呵呵呵!
   100个左右的原子在真空中,之间通过这种范德华力相联结。应该怎么样!

  大家想在着空中 放30个水分子组成的团簇。肯定会炸开。只是在水分子内部由于有化学键的存在,不会断开,水分子之间的距离肯定会迅速增大! 
我认为,酒一口一口喝,路一步一步走~步子迈大了,喀~容易扯着蛋
6楼2009-03-08 14:38:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)

★ ★
yuhuobuku(金币+2):欢迎参加讨论
”炸开“和”不炸开“,不在于势函数模型,恐怕取决于参数。如果你LJ势中的吸引能大一些,也许就不会”炸开“了。同时温度大概也是一个比较重要的参数。

一般来说,求算金属体系,最好还是用EAM势之类的势函数。
7楼2009-03-08 16:59:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huixiecn

木虫 (正式写手)

引用回帖:
Originally posted by xuefei06 at 2009-3-8 13:45:
EAM计算没有炸开,而lJ势炸开了,这个问题可能有很多可能,
   1.LJ的势参数的选择是不是有问题,降低排斥,增加吸引势看看如何
   2.体系控制参数温度,步长,截断半径是不是有问题

应该是LJ参数选择的不对

EAM没试过,但是LJ是不会炸开的
8楼2009-04-14 11:29:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lqllinqiaoli 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见