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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » lammps如何实现耗散粒子动力学控温
81/1返回列表
查看: 2473  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者TW九九将赠送您 5 个金币

TW九九

铜虫 (初入文坛)

[求助] lammps如何实现耗散粒子动力学控温

想问下lammps如何采用耗散粒子动力学控温方式，nvt、npt这种好像都用的Nose-Hoover控温，我看了几个关于dpd的命令感觉都不是。

lammps如何实现耗散粒子动力学控温
微信图片_20201111160502.png
回复此楼

» 猜你喜欢
1楼 2020-11-11 16:09:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luyaobao

木虫 (著名写手)
用pair_style dpd/tstat

发自小木虫Android客户端
赞一下(1人) 回复此楼
2楼2020-11-11 18:47:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

TW九九

铜虫 (初入文坛)
好的，这几天来我这看到了这个命令，但还是一直不确定，因为之前控温都是用的fix命令，而且复现的文献中说采用的LJ势，是不是可以用 pair_style hybrid  lj/cut  ……  dpd/tstat ……命令
赞一下 回复此楼
3楼2020-11-15 10:51:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luyaobao

木虫 (著名写手)
引用回帖:
3楼: Originally posted by TW九九 at 2020-11-15 10:51:34
好的，这几天来我这看到了这个命令，但还是一直不确定，因为之前控温都是用的fix命令，而且复现的文献中说采用的LJ势，是不是可以用 pair_style hybrid  lj/cut  ……  dpd/tstat ……命令

是

发自小木虫Android客户端
4楼2020-11-15 14:37:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

TW九九

铜虫 (初入文坛)
引用回帖:
4楼: Originally posted by luyaobao at 2020-11-15 07:37:59
是
...

谢谢，那对应于pair_style dpd/tstat的pair_coeff 参数是什么呢，我看手册说分别是A (force units)、γ(force/velocity units)、cutoff (distance units)，其中的γ是阻尼系数吗，这个数应该怎么得到呢，我做的是氩原子的。
赞一下 回复此楼
5楼2020-11-16 19:10:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luyaobao

木虫 (著名写手)
引用回帖:
5楼: Originally posted by TW九九 at 2020-11-16 19:10:28
谢谢，那对应于pair_style dpd/tstat的pair_coeff 参数是什么呢，我看手册说分别是A (force units)、γ(force/velocity units)、cutoff (distance units)，其中的γ是阻尼系数吗，这个数应该怎么得到呢，我做的是氩 ...

温控没有系数A，γ不同液体粘度不同，看你的体系是什么了

发自小木虫Android客户端
赞一下 回复此楼
6楼2020-11-16 19:37:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

TW九九

铜虫 (初入文坛)
引用回帖:
6楼: Originally posted by luyaobao at 2020-11-16 12:37:35
温控没有系数A，γ不同液体粘度不同，看你的体系是什么了
...

units                  metal
dimension                     3
boundary                     f f p

atom_style                   atomic
neighbor                     2.0 bin
neigh_modify                 every 1  delay 5 check yes

region                             box block 0 169 0 168 0 400
create_box                   2 box

region                       wall_left block 0 4 0 168 0 400
region                       wall_right block 164 169 0 168 0 400
region                 wall_down block 4 164 0 4 0 400
region                 flow block 4 164 4 44 0 400

lattice                             fcc  5.738
create_atoms                 1 region flow

lattice                             fcc  3.483
create_atoms                 1 region wall_left
create_atoms                 1 region wall_right
create_atoms                 1 region wall_down

mass                               1 39.948
mass                               2 39.948

pair_style                   hybrid  lj/cut 8.5  dpd/tstat 86 86 8.5 34387
pair_coeff                   1 1 lj/cut 0.0103125 3.4
pair_coeff                   1 2 lj/cut 0.00515625 3.4
pair_coeff                   2 2 lj/cut 0.00515625 3.4
pair_coeff                   1 1 dpd/tstat 0.0672 8.5
comm_modify            vel yes

group                         1 region wall_left
group                         2 region wall_right
group                         3 region wall_down
group                         wall union 1 2 3
group                  flow region flow
group                         Ar subtract all wall
set                           group wall type 2

velocity               flow create 86 482748 mom yes rot yes
fix                           1 all nvt temp 86 86 0.02
compute                1 Ar chunk/atom bin/2d x lower 2 y lower 2 units box

velocity               wall set 0.0 0.0 0.0
fix                           2 wall setforce 0.0 0.0 0.0
fix                    3 flow ave/chunk 100 5 500 1 density/mass file c_0.5.profile

timestep                     0.005
thermo                             2000
thermo_style           custom step temp pe press density etotal

dump                   1 all custom 1000 c_0.5.dump id type x y z
run                    310000


想问您以下为什么我这个in文件输出的结果从始至终原子都不动，我文件中也设置了velocity命令了，您知道吗
赞一下 回复此楼
7楼2020-11-20 20:34:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

TW九九

铜虫 (初入文坛)
引用回帖:
7楼: Originally posted by TW九九 at 2020-11-20 13:34:05
units                  metal
dimension                     3
boundary                     f f p

atom_style                   atomic
neighbor                     2.0 bin
neigh_modify                 every 1  delay 5 check yes
...

写错了，我的in文件没有fix nvt那一条命令，然后原子始终不动，后来我试着加上nvt，又动的太快了，几千步就原子丢失了
赞一下 回复此楼
8楼2020-11-22 10:40:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 TW九九 的主题更新
81/1返回列表
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿中南大学化学学硕0703总分337求调剂 +7 niko- 2026-03-22 7/350 2026-03-25 20:14 by qingfeng258
[考研] 材料与化工328分调剂 +5 。，。，。，。i 2026-03-23 5/250 2026-03-25 18:40 by xcjcqu
[考研] 085600材料与化工调剂 +9 A-哆啦Z梦 2026-03-23 15/750 2026-03-25 11:18 by Ainin_
[考研] 278求调剂 +4 我可以上岸的对� 2026-03-19 4/200 2026-03-25 11:01 by userper
[考研] 289材料与化工（085600）B区求调剂 +4 这么名字咋样 2026-03-22 5/250 2026-03-25 08:20 by mx.yue
[考研] 食品专硕 一志愿双一流 328 +3 xiaom99 2026-03-21 4/200 2026-03-24 21:20 by lailaisimei
[考研] 求调剂 +5 林之夕 2026-03-24 5/250 2026-03-24 17:16 by dick_runner
[考研] 一志愿华东理工大学081700，初试分数271 +5 kotoko_ik 2026-03-23 6/300 2026-03-24 10:29 by 学术搬砖er
[考研] 一志愿河北工业大学0817化工278分求调剂 +7 jhybd 2026-03-23 12/600 2026-03-24 09:03 by jhybd
[考研] 284求调剂 +3 yanzhixue111 2026-03-23 6/300 2026-03-23 22:58 by pswait
[考研] 276求调剂 +3 YNRYG 2026-03-21 4/200 2026-03-23 08:31 by 醉在风里
[考研] 317求调剂 +12 申子申申 2026-03-19 18/900 2026-03-22 22:23 by luoyongfeng
[考研] 293求调剂 +3 涛涛Wjt 2026-03-22 5/250 2026-03-22 22:21 by jiangpengfei
[考研] 319求调剂 +4 小力气珂珂 2026-03-20 4/200 2026-03-22 15:53 by ColorlessPI
[考研] 285求调剂 +6 ytter 2026-03-22 6/300 2026-03-22 12:09 by 星空星月
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[考研] 材料学硕318求调剂 +5 February_Feb 2026-03-19 5/250 2026-03-19 23:51 by 23Postgrad
信息提示
关闭
请填处理意见
关闭