应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 【求助】流体的分子动力学模拟
81/1返回列表
查看: 794  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

lqllinqiaoli

至尊木虫 (著名写手)

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

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

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

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2009-02-22 19:35:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyy816

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

老虎大王

木虫 (著名写手)
★ ★
yuhuobuku(金币+2,VIP+0):欢迎参加讨论 2-22 21:04
LJ势是对势，对势总是预言金属的弹性常数存在Cauchy关系，即C12=C44，这对于过渡金属是不成立的。对它若做些简单修正，如添加一个与体积变化相关的项，可以解决这个困难，但又带来了一些新的困难。嵌入原子法则顺利地解决了这个问题，所以嵌入原子法是当前用于金属，尤其是过渡金属的最流行的方法。但也不是唯一的方法。同时，LJ势及其修正形式也不是说完全不可以用于金属，全看你要模拟什么性质。
赞一下(6人) 回复此楼
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.体系控制参数温度，步长，截断半径是不是有问题
赞一下(10人) 回复此楼
5楼2009-03-08 13:45:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

superdirac

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

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

老虎大王

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

一般来说，求算金属体系，最好还是用EAM势之类的势函数。
赞一下(1人) 回复此楼
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 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 cqv1zKnzd7hr 2026-11-09 6/300 2026-11-11 04:12 by Z7qfCsdIiD1T
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 04:07 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 zSbvVX9JFl9J 2026-11-09 6/300 2026-11-11 03:59 by Z7qfCsdIiD1T
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 FChKX8meu1ds 2026-11-09 6/300 2026-11-11 01:54 by Z7qfCsdIiD1T
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 01:13 by Z7qfCsdIiD1T
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 01:10 by Z7qfCsdIiD1T
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 FChKX8meu1ds 2026-11-09 4/200 2026-11-10 20:51 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 cqv1zKnzd7hr 2026-11-09 5/250 2026-11-10 20:12 by Z7qfCsdIiD1T
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 cqv1zKnzd7hr 2026-11-09 4/200 2026-11-10 20:07 by Z7qfCsdIiD1T
[考博] 27博士申请 +3 白白白ly 2026-11-09 3/150 2026-11-09 21:39 by szq0819
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 15/750 2026-11-09 13:19 by Equinoxhua
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +6 5sJFDKp2kwuI 2026-11-06 14/700 2026-11-09 11:52 by GwrQVVR6Oenx
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +5 5sJFDKp2kwuI 2026-11-06 16/800 2026-11-09 11:47 by GwrQVVR6Oenx
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
信息提示
关闭
请填处理意见
关闭