24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 3369  |  回复: 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的回帖

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的回帖

一条笨鱼

新虫 (初入文坛)

引用回帖:
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的回帖

cuso4852cuoh

新虫 (小有名气)

fix 1 all ave/time 1 1 1 c_vacf[4] file vacf.dat   只需要看第四列数吗??
6楼2019-01-07 11:37:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

三岁薇恩

新虫 (正式写手)

引用回帖:
6楼: Originally posted by cuso4852cuoh at 2019-01-07 11:37:38
fix 1 all ave/time 1 1 1 c_vacf file vacf.dat   只需要看第四列数吗??

你好 请问你会计算速度自相关函数了吗 求指点
7楼2019-11-28 17:31:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
本帖内容被屏蔽

8楼2019-12-31 17:46:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
本帖内容被屏蔽

9楼2020-01-03 11:03:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
本帖内容被屏蔽

10楼2020-01-07 11:07:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 WJ-MS 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 26药学专硕105500求调剂 +6 喽哈加油 2026-04-13 6/300 2026-04-14 16:40 by zhouxiaoyu
[考研] 材料专业344求调剂 +17 hualkop 2026-04-10 22/1100 2026-04-14 16:21 by sxdj2
[考研] 食品与营养(0955)271求调剂 +15 升格阿达 2026-04-12 16/800 2026-04-14 13:18 by 浮若_安生
[考研] 085600材料与化工349分求调剂 +16 李木子啊哈哈 2026-04-12 17/850 2026-04-14 09:11 by fenglj492
[考研] 机械还有还有名额吗?太难了 +8 笑笑袁 2026-04-10 8/400 2026-04-14 08:44 by screening
[考研] 求调剂,985材料与化工348分 +9 涵竹刘 2026-04-11 14/700 2026-04-13 22:26 by 涵竹刘
[考研] 297工科,求调剂? +13 河南农业大学-能 2026-04-12 13/650 2026-04-13 14:12 by dingyanbo1
[考研] 2本,初试303,0860求调剂 +6 floriea 2026-04-12 8/400 2026-04-12 18:13 by zhouxiaoyu
[考研] 0831生医工第一轮调剂失败求助 +12 小熊睿睿_s 2026-04-11 16/800 2026-04-12 16:28 by 钰璞
[考研] 296求调剂 +14 汪!?! 2026-04-10 16/800 2026-04-12 10:48 by zhouyuwinner
[考研] 267求调剂 +8 再忙也要吃饭啊 2026-04-09 8/400 2026-04-11 21:42 by cfdbai
[考研] 农学0904 312求调剂 +6 Say Never 2026-04-10 6/300 2026-04-11 10:33 by wwj2530616
[考研] 085410-273求调剂 +6 X1999 2026-04-10 6/300 2026-04-11 10:32 by Delta2012
[考研] 机械专硕270求调剂,接受跨专业 +12 老师看看我吧aba 2026-04-09 14/700 2026-04-11 10:21 by laoshidan
[考研] 085402通信工程调剂,有4项学科竞赛国奖(电赛国二),硕士研究生调剂自荐信。 +5 m永o不v言o弃m 2026-04-09 5/250 2026-04-11 09:33 by zhq0425
[考研] 一志愿中科大070300化学,314分求调剂 +12 wakeluofu 2026-04-09 12/600 2026-04-10 09:57 by liuhuiying09
[考研] 材料专硕初试分332一志愿西北工业大学, +12 故人?? 2026-04-09 12/600 2026-04-09 18:34 by Ccclqqq
[论文投稿] 求助文献原文 10+3 18500821399 2026-04-08 3/150 2026-04-09 16:56 by 北京莱茵润色
[考研] 286求调剂 +19 Faune 2026-04-08 20/1000 2026-04-09 08:36 by 哦哦123
[考研] 生物学学硕,初试351分,求调剂 +4 …~、王…~ 2026-04-08 5/250 2026-04-08 21:49 by limeifeng
信息提示
请填处理意见