24小时热门版块排行榜    

查看: 1470  |  回复: 6

yiniangeqiu

新虫 (初入文坛)

[求助] 【求助】MD模拟-固体液体间L-J potential 参数-20个金币

【求助】MD模拟: 两不溶液体界面处放置一个或一些nanorod后,看看液体界面变化。
全部采用L-J  potential。 设它们为液体1, 液体2, 固体nanorod 3。
我现在面临的问题:
1。 如何选取3-3; 1,2-3的potential的参数:epsilon,sigma,cutoff? 因为我要保证固体particle不会跑出自己的边界;液体原子不会进入到固体原子中。

2。 液体液体界面,目前看起来可以搞定了。我用的是1-1, 2-2 参数都可以取相同, epsilson, sigma=1.0; 1.0。但是1-2 epsilon, sigma=0.1 1.0.

求大神指导上面固体液体间L-J potential 参数阿!
解决问题的,20个金币答谢!其他5个金币答谢。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaoyixiong

荣誉版主 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
chaizhm: 金币+2, 谢谢~ 2012-12-02 10:04:49
yiniangeqiu: 金币+5, ★★★很有帮助 2012-12-02 12:01:00
“要保证固体particle不会跑出自己的边界”固体颗粒之间应该还成键吧,固体颗粒内部原子之间成键 ,就不会跑出自己的边界。
其他参数的选择,你应该多找找相关的文献
2楼2012-12-02 09:36:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chaizhm

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
yiniangeqiu: 金币+5, ★★★很有帮助 2012-12-02 12:01:19
zh1987hs: 金币+2, 谢谢 2012-12-05 08:04:39
我感觉楼上说的成键有点不好判断,因为lj势本来就没考虑化学键。

只要势参数选择合适的话,应该两个就不会太多扩散,只要是表层原子之间有相互作用
壁立万仞,无欲则刚
3楼2012-12-02 10:08:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yiniangeqiu

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by jiaoyixiong at 2012-12-02 09:36:15
“要保证固体particle不会跑出自己的边界”固体颗粒之间应该还成键吧,固体颗粒内部原子之间成键 ,就不会跑出自己的边界。
其他参数的选择,你应该多找找相关的文献

这只是个课程作业, 不是自己科研用。我可以做得simple一点,不必要真实的数据,但得合理。
如果加bond固定的话,怎么加啊?LAMMPS,软件bond,能只对一种原子间(固体原子间)加么?
4楼2012-12-02 11:55:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yiniangeqiu

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by chaizhm at 2012-12-02 10:08:09
我感觉楼上说的成键有点不好判断,因为lj势本来就没考虑化学键。

只要势参数选择合适的话,应该两个就不会太多扩散,只要是表层原子之间有相互作用

目前我先用了一种固体,一种液体(box 各一半)来实验:固体液体边界是明显存在,但是,一种原子扩散把另外一种原子挤压在一个很小区域了。也不是我想要的固体液体的情况。  请问,您有没有这lj参数方面的经验呢。
5楼2012-12-02 11:58:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaoyixiong

荣誉版主 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
yiniangeqiu: 金币+5 2012-12-02 23:46:19
zh1987hs: 金币+2, 谢谢 2012-12-05 08:04:53
引用回帖:
3楼: Originally posted by chaizhm at 2012-12-02 10:08:09
我感觉楼上说的成键有点不好判断,因为lj势本来就没考虑化学键。

只要势参数选择合适的话,应该两个就不会太多扩散,只要是表层原子之间有相互作用

正是因为LJ势没有考虑化学键,所以才需要把一个纳米颗粒内部原子之间加上化学键啊,单考虑LJ势,那么在模拟的过程中,所有的原子之间都是散开的啊
6楼2012-12-02 20:18:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

polypro

木虫 (正式写手)

【答案】应助回帖

★ ★ ★
月只蓝: 金币+3, 感谢指导! 2013-11-29 08:54:52
单纯LJ势也是可能形成液固界面的,要求固相粒子体积分数大于越53%(玻璃态)或更高(结晶态),且固相内粒子之间的距离处于LJ的吸引势阱中。而对于液相,体积分数要小于临界值。
泉涸,鱼相与处于陆,相呴以湿,相濡以沫,不如相忘于江湖。
7楼2013-11-29 08:45:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yiniangeqiu 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 FChKX8meu1ds 2026-11-09 5/250 2026-11-11 04:29 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 04:10 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.5.5.1.O.54,科目齐全,可+急 +4 zSbvVX9JFl9J 2026-11-09 5/250 2026-11-11 04:01 by Z7qfCsdIiD1T
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 zSbvVX9JFl9J 2026-11-09 6/300 2026-11-11 03:56 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一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 cqv1zKnzd7hr 2026-11-09 5/250 2026-11-11 01:24 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 cqv1zKnzd7hr 2026-11-09 6/300 2026-11-11 01:21 by Z7qfCsdIiD1T
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 FChKX8meu1ds 2026-11-09 4/200 2026-11-10 22:39 by Z7qfCsdIiD1T
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 FChKX8meu1ds 2026-11-09 4/200 2026-11-10 20:51 by Z7qfCsdIiD1T
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 zSbvVX9JFl9J 2026-11-09 4/200 2026-11-10 01:29 by FChKX8meu1ds
[考博] 27博士申请 +3 白白白ly 2026-11-09 3/150 2026-11-09 21:39 by szq0819
[文学芳草园] 厉娜 +3 张丹枫1991 2026-11-09 5/250 2026-11-09 16:53 by 张丹枫1991
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 5sJFDKp2kwuI 2026-11-06 16/800 2026-11-09 11:47 by GwrQVVR6Oenx
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 08:00 by IhnEQauGmT3K
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售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 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
信息提示
请填处理意见