24小时热门版块排行榜    

查看: 3307  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 北科281学硕材料求调剂 +6 tcxiaoxx 2026-03-20 6/300 2026-03-22 20:23 by edmund7
[考研] 08工科 320总分 求调剂 +11 梨花珞晚风 2026-03-17 11/550 2026-03-22 17:42 by luoyongfeng
[考研] 070300化学求调剂 +3 苑豆豆 2026-03-20 3/150 2026-03-22 14:48 by ColorlessPI
[考研] 一志愿北京化工大学070300 学硕336求调剂 +5 vv迷 2026-03-21 8/400 2026-03-22 14:20 by ColorlessPI
[考研] 303求调剂 +5 安忆灵 2026-03-22 6/300 2026-03-22 12:46 by 素颜倾城1988
[考研] 求调剂 +7 Auroracx 2026-03-22 7/350 2026-03-22 12:38 by 素颜倾城1988
[考研] 初试 317 +7 半拉月丙 2026-03-20 7/350 2026-03-21 22:26 by peike
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[考研] 化学求调剂 +4 临泽境llllll 2026-03-17 5/250 2026-03-21 02:23 by JourneyLucky
[考研] 一志愿中国石油大学(华东) 本科齐鲁工业大学 +3 石能伟 2026-03-17 3/150 2026-03-21 02:22 by JourneyLucky
[考研] 332求调剂 +4 ydfyh 2026-03-17 4/200 2026-03-21 02:20 by JourneyLucky
[考研] 材料 336 求调剂 +3 An@. 2026-03-18 4/200 2026-03-21 01:39 by JourneyLucky
[考研] 271材料工程求调剂 +8 .6lL 2026-03-18 8/400 2026-03-21 00:58 by JourneyLucky
[考研] 085600材料与化工 +8 安全上岸! 2026-03-16 8/400 2026-03-20 22:13 by luoyongfeng
[考研] 材料学硕318求调剂 +5 February_Feb 2026-03-19 5/250 2026-03-19 23:51 by 23Postgrad
[考研] 081700化工学硕调剂 +3 【1】 2026-03-16 3/150 2026-03-19 23:40 by edmund7
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[考研] 本科郑州大学物理学院,一志愿华科070200学硕,346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 0854可跨调剂,一作一项核心论文五项专利,省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
信息提示
请填处理意见