24小时热门版块排行榜

返回列表

dgwfly

新虫 (正式写手)

应助: 0 (幼儿园)
金币: 683.9
散金: 258
红花: 9
帖子: 543
在线: 202小时
虫号: 3454332
注册: 2014-10-04
性别: MM
专业: 数论

[求助] 势函数的编写问题已有1人参与

老板让用fortran做分子动力学的模拟。之前有Lennard-Jones势能函数简单的NVE和NVT系宗的程序。但现在不知道怎么用S-W势函数或者其它描述固体的势函数来编写程序？主要是气固两相的模拟。求助相关Fortran代码的例子

回复此楼

» 收录本帖的淘帖专辑推荐

编程

» 猜你喜欢

上海工程技术大学张培磊教授团队招收博士生已经有3人回复
上海工程技术大学【激光智能制造】课题组招收硕士已经有5人回复
求助院士们，这个如何合成呀已经有4人回复
临港实验室与上科大联培博士招生1名已经有9人回复
想换工作。大多数高校都是评职称时认可5年内在原单位取得的成果吗？已经有7人回复
需要合成515-64-0，50g，能接单的留言已经有4人回复
自荐读博已经有4人回复
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？已经有6人回复
带资进组求博导收留已经有10人回复
最近几年招的学生写论文不引自己组发的文章已经有11人回复

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

» 本主题相关价值贴推荐，对您同样有帮助:

有关编程语言的选择已经有4人回复
势函数问题已经有8人回复
使用Multiwfn绘制原子轨道图形、研究原子壳层结构及相对论效应的影响已经有17人回复
【资源】著名物理学家吴大猷教授的理论物理丛书下载【珍藏版】已经有475人回复
回答一些关于Multiwfn的疑问以及未来Multiwfn的发展打算已经有37人回复
【sobereva个人文集】制作动画分析电子结构特征已经有29人回复
【分享】量化软件 - 输运相关 (zz) 已经有13人回复

1楼 2016-03-04 09:49:56

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

匿名

用户注销 (小有名气)

应助: 1 (幼儿园)
金币: 1383.8
散金: 210
红花: 4
帖子: 283
在线: 67小时
虫号: 0
注册: 2013-05-18
性别: GG
专业: 社会心理学

本帖仅楼主可见

2楼2016-03-04 12:11:05

已阅申请模拟EPI 回复此楼编辑查看我的主页

lsloneil

专家顾问 (正式写手)

专家经验: +49
模拟EPI: 2
应助: 148 (高中生)
金币: 8640.9
散金: 507
红花: 30
帖子: 505
在线: 598.4小时
虫号: 556471
注册: 2008-05-11
性别: GG
专业: 理论和计算化学
管辖: 分子模拟

【答案】应助回帖

感谢参与，应助指数 +1

Stillinger-Webber是三体势，意味着相比起两体势，要多一重循环，比方说你计算LJ势的module里，是
do i=1, N
do j=i, N
      计算i与j之间相互作用
end do
end do
对于SW，循环就变成
do i=1, N
do j=i, N
      do k=j,N
      计算i,j,k的三体势能
      end do
   end do
end do

SW因为是三体势，计算起来比二体势会慢不少。如果你用的程序不是优化得非常好，为了追求速度的话建议用现成的开源软件包，像LAMMPS， DL_POLY里面都有现成的SW势可以用。如果你想参考代码，可以下载LAMMPS(C++)或者DL_POLY (Fortran)，看它们的代码是怎么做的。