 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 910 | 回復(fù): 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 |
回復(fù)此樓» 猜你喜歡
西安電子科學(xué)大學(xué)杭州研究院劉麗香教授招收智能多模態(tài)傳感器和微型儲能器件方向博士
已經(jīng)有13人回復(fù)
-
重慶交通大學(xué)光子學(xué)微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息
已經(jīng)有1人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有131人回復(fù)
-
一志愿鄭大材料學(xué)碩298分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
尋合作：應(yīng)力腐蝕多尺度模擬
已經(jīng)有3人回復(fù)
-
考研交流
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
驕傲的九零后
銀蟲 (正式寫手)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 640.6
- 散金: 65
- 帖子: 319
- 在線: 104.4小時
- 蟲號: 1691240
- 注冊: 2012-03-14
- 性別: GG
- 專業(yè): 有機化工
2樓2024-12-23 11:57:46
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 驕傲的九零后 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 320材料與化工，求調(diào)劑 +5 | 鶴遨予卿 2026-03-04 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿武漢理工大學(xué)-085602-總分296分-求調(diào)劑 +5 | 紫川葡柚 2026-03-04 | 5/250 |
|
|
|
[考研] 學(xué)碩材料275調(diào)劑 +9 | 路三三 2026-03-03 | 9/450 |
|
|
[考研] 085600 材料與化工 298 +14 | 小西笑嘻嘻 2026-03-03 | 14/700 |
|
|
|
[考研] 0703 東華大學(xué) 理學(xué)化學(xué)均過 a 線，總分 281 求調(diào)劑 +3 | 香香 Lu 2026-03-03 | 5/250 |
|
|
|
[考研] 0857求調(diào)劑 +4 | CDDD1013 2026-03-04 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | 泡了個椒 2026-03-04 | 3/150 |
|
|
|
[考研] 304分材料專碩求調(diào)劑 +8 | qiuzhigril 2026-03-03 | 10/500 |
|
|
|
[考研] 材料工程269求調(diào)劑 +7 | 白刺玫 2026-03-02 | 7/350 |
|
|
|
[考研] 291求調(diào)劑 +4 | Afy123456 2026-03-03 | 7/350 |
|
|
|
[考研] 理學(xué)，工學(xué)，農(nóng)學(xué)調(diào)劑，少走彎路，這里歡迎您！ +8 | likeihood 2026-03-02 | 11/550 |
|
|
|
[考研] 267求調(diào)劑 +6 | 釣魚佬as 2026-03-02 | 6/300 |
|
|
|
[考研] 291求調(diào)劑 +3 | MuoLuo1312 2026-03-02 | 6/300 |
|
|
|
[考研] 281電子信息求調(diào)劑 +3 | jhtfeybgj 2026-03-02 | 6/300 |
|
|
|
[考研] 清華大學(xué) 材料與化工 353分求調(diào)劑 +5 | awaystay 2026-03-02 | 6/300 |
|
|
|
[考研] 306分材料調(diào)劑 +5 | chuanzhu川燭 2026-03-01 | 6/300 |
|
|
|
[考研] 材料085601調(diào)劑 +5 | 多多子. 2026-03-02 | 5/250 |
|
|
|
[考研] 一志愿東北大學(xué)材料專碩328，求調(diào)劑 +3 | shs1083 2026-03-02 | 3/150 |
|
|
|
[考研] 295求調(diào)劑 +8 | 19171856320 2026-02-28 | 8/400 |
|
|
|
[考研] 調(diào)劑 +3 | 13853210211 2026-03-02 | 4/200 |
|