 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 913 | 回復: 1 | |||
| 【懸賞金幣】回答本帖問題，作者驕傲的九零后將贈送您 10 個金幣 | |||
驕傲的九零后銀蟲 (正式寫手)
|
[求助]
求助 lammps上DPD模擬計算
|
||
|
各位蟲友好，最近在lammps上，用DPD計算模擬聚合物相分離，總是出現(xiàn)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 |
回復此樓» 猜你喜歡
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復
-
物理學I論文潤色/翻譯怎么收費?
已經(jīng)有70人回復
-
幫我的英語口語老師找學生
已經(jīng)有3人回復
-
計算三階力常數(shù)用thirdorder擴胞時參數(shù)如何選擇？臨近原子數(shù)怎么取值
已經(jīng)有5人回復
-
2026年循環(huán)經(jīng)濟功能材料國際會議（ICFMCE 2026）
已經(jīng)有0人回復
-
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）
已經(jīng)有0人回復
-
哈爾濱理工大學物理系招收物理學考研調(diào)劑
已經(jīng)有0人回復
驕傲的九零后
銀蟲 (正式寫手)
- 應助: 2 (幼兒園)
- 金幣: 640.6
- 散金: 65
- 帖子: 319
- 在線: 104.4小時
- 蟲號: 1691240
- 注冊: 2012-03-14
- 性別: GG
- 專業(yè): 有機化工
2樓2024-12-23 11:57:46
相關版塊跳轉(zhuǎn)
我要訂閱樓主 驕傲的九零后 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
 調(diào)劑的同學，走過路過，不要錯過
+5
|
likeihood 2026-03-06 | 9/450 |
|
|---|---|---|---|---|
|
|
[考研] 材料專碩調(diào)劑 +5 | 慕辰123 2026-03-05 | 8/400 |
|
|
|
[考博] 售瀆SCI文章，我:8⊙ 55 1⊙ 54。備注【⊙=0】 +4 | ipys00zm58 2026-03-05 | 9/450 |
|
|
|
[考博] 售瀆SCI文章，我:8⊙ 55 1⊙ 54。備注【⊙=0】 +4 | ipys00zm58 2026-03-06 | 6/300 |
|
|
[考研] 275求調(diào)劑 +4 | 大爆炸難民 2026-03-06 | 5/250 |
|
|
|
[考研] 085600，一志愿鄭州大學，280分求調(diào)劑 +7 | Wuqi725 2026-03-05 | 7/350 |
|
|
|
[考研] 復試調(diào)劑 +5 | 呼呼？~+123456 2026-03-05 | 5/250 |
|
|
|
[考研] 紡織、生物、化學、材料等專業(yè) +3 | Eember. 2026-03-05 | 7/350 |
|
|
|
[考研] 302材料工程求調(diào)劑 +7 | Doleres 2026-03-01 | 8/400 |
|
|
|
[考研] 085701環(huán)境工程 求調(diào)劑 +7 | xiiiia 2026-03-04 | 7/350 |
|
|
|
[考研] 材料085600 303求調(diào)劑 +7 | 1bygone 2026-03-04 | 7/350 |
|
|
|
[考研] 0703化學求調(diào)劑，總分320分，一志愿中南民族大學 +3 | 19286712440 2026-03-04 | 3/150 |
|
|
|
[考研] 材料工程269求調(diào)劑 +7 | 白刺玫 2026-03-02 | 7/350 |
|
|
|
[考研] 291求調(diào)劑 +4 | Afy123456 2026-03-03 | 7/350 |
|
|
|
[考研] 298求調(diào)劑一志愿中海洋 +3 | lour. 2026-03-03 | 3/150 |
|
|
|
[考研] 0703化學306調(diào)劑 +4 | 26要上岸 2026-03-03 | 4/200 |
|
|
|
[考研] 321求調(diào)劑一志愿東北林業(yè)大學材料與化工英二數(shù)二 +5 | 蟲蟲蟲蟲蟲7 2026-03-01 | 9/450 |
|
|
|
[考研] 261求調(diào)劑 +3 | 陸lh 2026-03-01 | 3/150 |
|
|
|
[考研] 284求調(diào)劑 +10 | 天下熯 2026-02-28 | 11/550 |
|
|
|
[考研] 299求調(diào)劑 +3 | Y墨明棋妙Y 2026-02-28 | 5/250 |
|