| 查看: 768 | 回复: 0 | ||
| 【悬赏金币】回答本帖问题，作者骄傲的九零后将赠送您 50 个金币 | ||
骄傲的九零后银虫 (正式写手)
|
[求助]
求助 lammps上DPD模拟计算
|
|
|
各位虫友好，最近在lammps上，用DPD计算模拟聚合物相分离，总是出现ERROR: Fix mvv/dpd requires atom attribute vest (src/DPD-MESO/fix_mvv_dpd.cpp:66)这句话。这里的vest是什么意思啊，找了不少资料也没解决，在这里卡了好久了，哪位虫友遇到过啊，下面是我的in文件。 units lj dimension 3 boundary p p p timestep 0.005 ############## Atom and Bond Styles ############## atom_style full read_data Equil-2-5-80.data ############## Atom and Bond Styles ############## bond_style harmonic bond_coeff * 10.0 1.0 bond_coeff 4 100.0 1.0 bond_coeff 1 60.0 1.0 angle_style hybrid harmonic cosine angle_coeff 4 harmonic 1000 120 angle_coeff 1 cosine 1 angle_coeff 2 cosine 1 angle_coeff 3 cosine 1 angle_coeff 5 cosine 1 angle_coeff 6 cosine 1 angle_coeff 7 cosine 1 angle_coeff 8 cosine 1 angle_coeff 9 cosine 1 angle_coeff 10 cosine 1 ############## Groups and Fixes ############## #Set neighbor list radius to 10 neighbor 3 bin neigh_modify one 40000 neigh_modify page 400000 neigh_modify every 1 delay 0 check yes comm_modify vel yes # Fix beads of type 3 as rigid bodies group yellow type 3 4 fix 2 yellow rigid/small molecule langevin 1 1 0.5 456447 ############## Tem compute ############## group my type 1 2 5 region boundary block -24 28 -24 28 -28 24 compute mine my temp/region boundary ############## Output Settings ############## thermo 1000 thermo_style custom step temp time ebond eangle edihed pe ke etotal lx ly lz density ############## Minimization ################ velocity all create 1 34387 #minimize 1e-6 1e-6 1000 1000 #write_data min.data ############## Soft Equilibration ############## dump 1 all custom 1000 pe.lammpstrj id mass mol type xu yu zu dump_modify 1 flush yes ############## Pair Styles ############## # Soft potential for initial relaxation pair_style soft 2.5 pair_coeff * * 10 # Initial NvT equilibration fix 3 all nvt temp 1 1 0.5 run 100000 unfix 3 write_data SEquil-2-5-80.data ############## DPD Equilibration ################ pair_style hybrid dpd 1 1 34387 lj/cut 2.5 #pair_style dpd T cutoff seed #pair_style style args #lj/cut args = cutoff pair_coeff 1 1 dpd 25 4.5 pair_coeff 1 2 dpd 60 4.5 pair_coeff 1 3 dpd 60 4.5 pair_coeff 1 4 dpd 60 4.5 pair_coeff 1 5 dpd 60 4.5 pair_coeff 2 2 dpd 25 4.5 pair_coeff 2 3 dpd 80 4.5 pair_coeff 2 4 dpd 80 4.5 pair_coeff 2 5 dpd 25 4.5 pair_coeff 3 3 dpd 25 4.5 pair_coeff 3 4 dpd 25 4.5 pair_coeff 3 5 dpd 80 4.5 pair_coeff 4 4 dpd 25 4.5 pair_coeff 4 5 dpd 80 4.5 pair_coeff 5 5 dpd 25 4.5 pair_coeff 3 3 lj/cut 2 1.2 pair_coeff 4 4 lj/cut 2 1.2 fix 4 all langevin 1 1 0.5 12345 ####################################################################Deform to the density of 3 g/cm^3########### fix 11 all mvv/dpd 0.65 run 1000000 write_data deformed.data run 1000000 write_data 1deformed.data run 100000 unfix 11 write_data 2deformed.data |
回复此楼» 猜你喜欢
杂志要求提供official language-editing certificate
已经有4人回复
-
深圳大学2026年秋博士招生-物理学-活性胶体方向-高永祥课题组
已经有17人回复
-
物理学I论文润色/翻译怎么收费?
已经有174人回复
-
论物质与能量的统一模型及物理现象解释
已经有0人回复
-
基于基元I统一理论的数学相关应用推导
已经有0人回复
-
基元I统一理论：宇宙本质、层级演化与修炼文明的本源规律
已经有1人回复
-
基元I理论下三大核心空间现象精准推导与细节解析
已经有0人回复
-
基于基元 I 统一理论的反重力理论推导
已经有0人回复
-
基于基元I统一理论的量子力学本源推导
已经有1人回复
-
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，AI润色，全免费
已经有15人回复
-
【EI|Scopus 双检索】第六届智能机器人系统国际会议（ISoIRS 2026）
已经有0人回复