24小时热门版块排行榜    

Znn3bq.jpeg
查看: 3352  |  回复: 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的回帖

三岁薇恩

新虫 (正式写手)

引用回帖:
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的回帖
查看全部 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 290调剂生物0860 +23 哇哈哈,。 2026-04-11 25/1250 2026-04-12 10:22 by 猪会飞
[考研] 305求调剂 +7 玛卡巴卡boom 2026-04-11 7/350 2026-04-12 07:35 by zhouxiaoyu
[考研] 药学专硕调剂 +8 ? 一路生?花? 2026-04-10 10/500 2026-04-11 21:21 by zhouxiaoyu
[考研] 366求调剂 +7 不知名的小卅 2026-04-11 7/350 2026-04-11 16:45 by PeggyPeng17
[考研] 求调剂 +3 胃痉挛累了 2026-04-11 5/250 2026-04-11 14:13 by luhong1990
[硕博家园] 新一代电子信息294求调剂 不挑学校 +6 Ytyt11 2026-04-09 7/350 2026-04-11 10:52 by AA小小木虫
[考研] 300分求调剂 (085501机械专硕,本科扬大) +8 xu@841019 2026-04-11 8/400 2026-04-11 10:46 by qingpingzhu
[考研] 085410-273求调剂 +6 X1999 2026-04-10 6/300 2026-04-11 10:32 by Delta2012
[材料工程] 材料调剂推荐 +8 蛋糕x2 2026-04-07 8/400 2026-04-10 23:13 by Ftglcn90
[考研] 336材料与化工085600求调剂 +21 水星记infp 2026-04-05 24/1200 2026-04-10 15:28 by luoyongfeng
[考研] 一志愿中南大学物理学,英一66,求调剂 +4 长烟旖旎 2026-04-08 5/250 2026-04-10 10:31 by 颖果儿
[考研] 求机械专硕297第二批调剂 +5 拾柒12。 2026-04-08 5/250 2026-04-09 16:43 by 允当适度
[考研] 274求调剂 +5 山阿蔓 2026-04-07 5/250 2026-04-09 15:28 by 18828373951
[考研] 349学科化学045106求调剂,化学类都可以 +8 保好懂懂 2026-04-08 8/400 2026-04-09 14:03 by xulei3024
[考研] 材料工程322 +18 哈哈哈吼吼吼哈 2026-04-07 19/950 2026-04-09 10:44 by cymywx
[考研] 一志愿985初试354分生物调剂 +3 031001 2026-04-06 3/150 2026-04-09 00:30 by Evan_Liu
[考研] 277求调剂 +4 考研调剂lxh 2026-04-06 6/300 2026-04-08 10:40 by 逆水乘风
[考研] 材料调剂 +11 一样YWY 2026-04-07 11/550 2026-04-07 15:13 by shdgaomin
[考研] 328求调剂 +4 ghhh88888 2026-04-06 5/250 2026-04-07 14:45 by ghhh88888
[考研] 调剂 +3 李广火 2026-04-05 3/150 2026-04-05 18:57 by 蓝云思雨
信息提示
请填处理意见