24小时热门版块排行榜    

Znn3bq.jpeg
查看: 3358  |  回复: 11
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

WJ-MS

新虫 (小有名气)

[求助] lammps中速度自相关函数的计算

14年lammps里面增加了compute vacf命令可以计算速度自相关函数,这个命令计算的结果要怎么分析呢?它的结果也不能进行用trap()命令进行积分,貌似是什么变量性质不匹配之类的。另外,fix ave/correlate 命令计算的结果可以分析相关函数的衰减性什么的,compute vacf怎么做出速度自相关函数的曲线来呢,是要把每一步的计算结果保存下来么?

我另外尝试了直接写命令计算速度子先关函数,vx是per-atom量,fix ave/correlate命令不能用,采用compute reduce命令转换一下,问题是 compute reduce有四个模式,sum/max/min/ave,我试了一下,这四个模式得到的结果都不同,那么该怎么样把vx,vy, vz 转换为fix ave/correlate能够采用的属性呢?

“compute          1 all reduce max vx vy vz
variable         v1 equal c_1[1]
variable         v2 equal c_1[2]
variable         v3 equal c_1[3]
fix                 VV all ave/correlate $s $p $d v_v1 v_v2 v_v3 type auto file V0Vt1100K.dat ave running”

希望大家帮忙答疑解惑!!!
回复此楼
加油!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

一条笨鱼

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by challenger9 at 2016-06-24 10:21:21
我每次用这个命令总是提示Invalid compute style (modify.cpp:812)
,求助,谢谢!...

你那可能是lammps版本过低了。。去官网下个新的就好了应该。
5楼2016-09-29 15:36:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

JackieXiao

铜虫 (初入文坛)

compute vacf all vacf
fix 1 all ave/time 1 1 1 c_vacf[4] file vacf.dat
试试看
2楼2014-10-18 10:16:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

WJ-MS

新虫 (小有名气)

引用回帖:
2楼: Originally posted by JackieXiao at 2014-10-18 10:16:40
compute vacf all vacf
fix 1 all ave/time 1 1 1 c_vacf file vacf.dat
试试看

恩恩 我试一试
加油!
3楼2014-10-19 15:01:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

challenger9

银虫 (初入文坛)

引用回帖:
2楼: Originally posted by JackieXiao at 2014-10-18 10:16:40
compute vacf all vacf
fix 1 all ave/time 1 1 1 c_vacf file vacf.dat
试试看

我每次用这个命令总是提示Invalid compute style (modify.cpp:812)
,求助,谢谢!
4楼2016-06-24 10:21:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 346分,工科0854求调剂,专硕 +6 moser233 2026-04-12 7/350 2026-04-12 22:11 by fqwang
[考研] 268分085602化学工程调剂 +30 月照花林。 2026-04-09 30/1500 2026-04-12 19:17 by laoshidan
[考研] 电子信息270求调剂 +18 terminal469 2026-04-07 18/900 2026-04-12 16:23 by ajpv风雷
[考研] 296求调剂 +14 汪!?! 2026-04-10 16/800 2026-04-12 10:48 by zhouyuwinner
[考研] 求调剂 +10 璃茉一定上岸 2026-04-10 10/500 2026-04-11 13:31 by 1005715100
[考研] 22408调剂求助 +7 毂12 2026-04-09 9/450 2026-04-11 09:23 by 哦哦123
[考研] 284求调剂 +12 archer.. 2026-04-10 13/650 2026-04-11 08:44 by zhq0425
[考研] 一志愿北理工298英一数二已上岸,感谢各位老师 +14 Reframe 2026-04-10 16/800 2026-04-10 23:07 by caotw2020
[考研] 263能源动力专硕求调剂 +3 加大号饭盒袋 2026-04-10 3/150 2026-04-10 22:23 by 286640313
[考研] 263能源动力专硕求调剂 +4 加大号饭盒袋 2026-04-10 4/200 2026-04-10 20:52 by gong120082
[考研] 0702物理学学硕299求调剂 +6 祁柒连 2026-04-06 6/300 2026-04-10 11:10 by Roomoo
[考研] 调剂申请086000一志愿西北农林科技大学生物与医药320分-本科齐鲁工业大学 +3 美美女士 2026-04-09 3/150 2026-04-10 10:31 by liuhuiying09
[考研] 070300化学 求调剂 +13 73372112 2026-04-08 13/650 2026-04-09 20:22 by maddjdld
[考研] 材料工程调剂 +12 小刘同学吖吖 2026-04-06 13/650 2026-04-09 17:07 by luoyongfeng
[考研] 求调剂材料科学与工程一志愿985初试365分 +5 材化李可 2026-04-08 5/250 2026-04-09 17:00 by Lilly_Li
[考研] 349学科化学045106求调剂,化学类都可以 +8 保好懂懂 2026-04-08 8/400 2026-04-09 14:03 by xulei3024
[考研] 296求调剂 +3 汪!?! 2026-04-08 3/150 2026-04-08 22:00 by zhouyuwinner
[考研] 生物学328分求调剂 +9 闪电kkl 2026-04-08 10/500 2026-04-08 21:42 by liuhuiying09
[考研] 一志愿吉大化学327求调剂 +12 王王白石 2026-04-06 13/650 2026-04-08 16:05 by luoyongfeng
[考研] 372分材料与化工(085600)英二数二求调剂 +4 蓝笺片 2026-04-06 4/200 2026-04-07 12:30 by dongzh2009
信息提示
请填处理意见