| 查看: 2204 | 回复: 2 | ||
奈落之夜木虫 (正式写手)
|
[求助]
关于lammps 质心速度vcm（）如何输出问题
|
| 新手刚学，要输出质心速度 想用vcm这个函数 根据说明书 我先定义了变量 variable v_m equal vcm(motion,v_m) motion是一个group 然后输出 dump v_data motion cfg 10 dump.snap.*.cfg $v_m 可能很多语法错误 manual里面vcm（group，dim）前面应该是定义的组，后面的自变量不知道是什么定了，也不知道该怎么输出 |
回复此楼» 猜你喜欢
本科南方医科大学 一志愿985 药学学硕284分 求调剂
已经有5人回复
-
295分求调剂
已经有12人回复
-
调剂求收留
已经有25人回复
-
272分材料子求调剂
已经有49人回复
-
211本科材料化工求调剂
已经有16人回复
-
279学硕食品专业求调剂院校
已经有12人回复
-
一志愿双非085400电子信息344 求调剂，对材料和化学方向也感兴趣
已经有11人回复
-
279求调剂
已经有8人回复
-
277求调剂
已经有21人回复
-
290求调剂
已经有15人回复
bluewhale
铁杆木虫 (正式写手)
- 应助: 58 (初中生)
- 金币: 8438.3
- 散金: 50
- 红花: 10
- 帖子: 890
- 在线: 362.2小时
- 虫号: 278289
- 注册: 2006-09-10
- 专业: 凝聚态物性I:结构、力学和
【答案】应助回帖
★ ★ ★ ★
感谢参与，应助指数 +1
月只蓝: 金币+2, 感谢指导！ 2013-11-01 08:35:56
奈落之夜: 金币+2, ★有帮助 2013-11-01 09:10:33
感谢参与，应助指数 +1
月只蓝: 金币+2, 感谢指导！ 2013-11-01 08:35:56
奈落之夜: 金币+2, ★有帮助 2013-11-01 09:10:33
|
参照下面的例子： heated Brownian particle # typical input script used to investigate the Brownian motion of a heated nanoparticle. The diffusion coefficient is computed as the plateau value of the integral of the velocity autocorrelation function. Details in: # L. Joly, S. Merabia and J.-L. Barrat, EPL 94, 50007 (2011) # http://iopscience.iop.org/0295-5075/94/5/50007/ # the data file is generated by the perl script 'data.pl'. # fix ave/sph is similar to fix ave/spatial but using spherical bins (enclosed sources might be broken with latest lammps versions). # to apply a rescaling thermostat only to liquid atoms lying beyond 15 sigma from the center of the particle (see article for details), fix recenter is used to constrain the center-of-mass of the particle to remain at the origin, with keyword 'shift all' so that it does not alter the dynamics of the system. atom_style bond read_data data.brown group liq type 1 group sol type 2 region far sphere 0.0 0.0 0.0 15.9 side out units box region int1 sphere 0.0 0.0 0.0 14.4 side out units box region int2 sphere 0.0 0.0 0.0 15.9 units box region int intersect 2 int1 int2 units box region ins sphere 0.0 0.0 0.0 14.4 units box compute Tsol sol temp/com compute Tfar liq temp/region far compute Tint liq temp/region int group far region far group int region int group ins region ins group insliq intersect ins liq neigh_modify delay 0 every 1 check yes pair_style lj/cut 2.5 pair_modify shift yes bond_style fene special_bonds fene mass * 1.0 pair_coeff * * 1.0 1.0 bond_coeff * 30.0 1.5 1.0 1.0 variable ts equal step variable vx equal vcm(sol,x)-vcm(liq,x) variable vy equal vcm(sol,y)-vcm(liq,y) variable vz equal vcm(sol,z)-vcm(liq,z) fix 3 all print 5 "${ts} ${vx} ${vy} ${vz}" file print.brown screen no title "# ts vx vy vz" fix 4vx all ave/correlate 5 2401 1200000 v_vx type full ave running start 200000 file cvx.brown fix 4vy all ave/correlate 5 2401 1200000 v_vy type full ave running start 200000 file cvy.brown fix 4vz all ave/correlate 5 2401 1200000 v_vz type full ave running start 200000 file cvz.brown thermo_style custom step cpu c_Tsol c_Tint c_Tfar press vol thermo 1000 compute KEatom all ke/atom variable Tatom atom (2.0/3.0)*c_KEatom fix 5 all ave/sph 100 2000 200000 x origin 0.1 density/number v_Tatom units box file sph.brown fix 1 liq nph aniso 0.015 0.015 5.0 dilate partial fix 1far liq temp/rescale 1 0.75 0.75 0.0 1.0 fix_modify 1far temp Tfar fix 1int liq temp/rescale 1 0.75 0.75 0.0 1.0 fix_modify 1int temp Tint fix 1sol sol nvt temp 3.5 3.5 0.5 fix_modify 1sol temp Tsol variable vol equal vol fix avevol all ave/time 10 8000 160000 v_vol variable l equal 0.5*((f_avevol)^(1./3.)) fix 2 sol recenter 0.0 0.0 0.0 shift all units box velocity sol create 0.75 36582 dist gaussian mom yes rot yes velocity int create 0.75 36582 dist gaussian mom yes rot yes velocity far create 0.75 36582 dist gaussian mom yes rot yes velocity insliq create 0.75 36582 dist gaussian mom yes rot yes run 160000 unfix 1 fix 1 liq nve fix dxyz liq deform 1 x final -$l $l y final -$l $l z final -$l $l units box unfix 2 fix 2 sol recenter 0.0 0.0 0.0 shift all units box run 10000 unfix dxyz run 30000 #variable F equal 10.0 #variable Fs equal $F/count(sol) #fix 9s sol addforce ${Fs} 0.0 0.0 #variable Fl equal -$F/count(liq) #fix 9l liq addforce ${Fl} 0.0 0.0 #unfix 2 #fix 2 sol recenter 0.0 0.0 0.0 shift all units box run 12000000 |
2楼2013-10-31 20:45:20
奈落之夜
木虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 1890.6
- 红花: 1
- 帖子: 393
- 在线: 73.9小时
- 虫号: 2550025
- 注册: 2013-07-17
- 专业: 凝聚态物性I:结构、力学和
3楼2013-11-01 09:11:14
