现在模拟碳管中金属元素的生长情况,要计算金属元素的msd,可是一直弄不通,请大家帮助看看不胜感激!输入文件中将原子分成两类:
group C type 1
group Au type 2
在计算Au的msd时,用了compute mymsd Au msd com yes
输出msd用
fix 1 Au ave/time 1 100 50000 c_mymsd[1] c_mymsd[2] c_mymsd[3] c_mymsd[4] file msd.tmp mode vector
可是报错:ERROR: Fix ave/time compute does not calculate an array;
后来将ave/time那句中的vector换成了scalar,倒是可以运行起来,但是结果感觉完全不对!