 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 2567 | 回復(fù): 1 | |||
| 【懸賞金幣】回答本帖問題，作者周三啊將贈(zèng)送您 10 個(gè)金幣 | |||
周三啊鐵蟲 (小有名氣)
|
[求助]
Mismatched compute in variable formula 已有1人參與
|
||
|
各位大佬好，我在使用ILP勢(shì)計(jì)算兩層原子之間的勢(shì)能，官網(wǎng)上說勢(shì)能計(jì)算可以分為排斥和吸引兩部分，并且可以單獨(dú)訪問，語句如下 This pair style tallies a breakdown of the total interlayer potential energy into sub-categories, which can be accessed via the compute pair command as a vector of values of length 2. The 2 values correspond to the following sub-categories: E_vdW = vdW (attractive) energy E_Rep = Repulsive energy To print these quantities to the log file (with descriptive column headings) the following commands could be included in an input script: compute 0 all pair ilp/graphene/hbn variable Evdw equal c_0[1] variable Erep equal c_0[2] thermo_style custom step temp epair v_Erep v_Evdw 但是我在計(jì)算時(shí)卻出現(xiàn)ERROR: Variable ILP: Mismatched compute in variable formula (../variable.cpp:1580)的錯(cuò)誤。希望大家可以解答一下，謝謝！�。。。�！  我的腳本如下： # Initialization units metal boundary p p p atom_style full processors * * 1 # domain decomposition over x and y # System and atom definition # we use 2 atom types so that inter- and intra-layer # interactions can be specified separately # read lammps data file read_data 3.2_.data mass 1 12.0107 # carbon mass (g/mole) | membrane mass 2 12.0107 # carbon mass (g/mole) | adsorbate # Separate atom groups group membrane type 1 group adsorbate type 2 ######################## Potential defition ######################## pair_style hybrid/overlay rebo ilp/graphene/hbn 16.0 #################################################################### pair_coeff * * rebo CH.airebo NULL C # chemical pair_coeff * * ilp/graphene/hbn BNCH.ILP C C # long range #################################################################### # Neighbor update settings neighbor 2.0 bin neigh_modify every 1 neigh_modify delay 0 neigh_modify check yes # calculate the COM variable adsxcom equal xcm(adsorbate,x) variable adsycom equal xcm(adsorbate,y) variable adszcom equal xcm(adsorbate,z) variable adsvxcom equal vcm(adsorbate,x) variable adsvycom equal vcm(adsorbate,y) variable adsvzcom equal vcm(adsorbate,z) #### Simulation settings #### timestep 0.001 velocity adsorbate create 300.0 12345 fix subf membrane setforce 0.0 0.0 0.0 fix thermostat all nve compute 0 all pair rebo compute 1 all pair ilp/graphene/hbn variable REBO equal c_0 variable ILP equal c_1 variable Evdw equal c_1[1] variable Erep equal c_1[2] ############################ # Output thermo 100 thermo_style custom step etotal pe ke v_REBO v_ILP v_Erep v_Evdw temp v_adsxcom v_adsycom v_adszcom v_adsvxcom v_adsvycom v_adsvzcom thermo_modify line one format float %.10f thermo_modify flush yes norm no lost warn dump 1 all custom 100 traj.lammpstrj id mol type xu yu zu dump_modify 1 format line "%7d %3d %3d %15.10g %15.10g %15.10g" flush yes ###### Run molecular dynamics ###### run 1000 |
回復(fù)此樓» 猜你喜歡
0703化學(xué)一志愿211 總分320求調(diào)劑
已經(jīng)有3人回復(fù)
-
0703求調(diào)劑
已經(jīng)有6人回復(fù)
-
（081700）化學(xué)工程與技術(shù)-298分求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天津大學(xué)材料與化工275求調(diào)劑
已經(jīng)有22人回復(fù)
-
0857 資源與環(huán)境 285分
已經(jīng)有6人回復(fù)
-
材料工程,326分，求調(diào)劑
已經(jīng)有5人回復(fù)
-
301求調(diào)劑
已經(jīng)有5人回復(fù)
-
334求調(diào)劑
已經(jīng)有3人回復(fù)
-
材料與化工085600調(diào)劑求老師收留
已經(jīng)有5人回復(fù)
-
一志愿華科071000生物學(xué) 338求調(diào)劑
已經(jīng)有4人回復(fù)
zhizhuceyici
金蟲 (初入文壇)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 578.8
- 帖子: 16
- 在線: 41.9小時(shí)
- 蟲號(hào): 4078417
- 注冊(cè): 2015-09-17
- 性別: MM
- 專業(yè): 固體力學(xué)
2樓2020-10-19 17:11:11
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 周三啊 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考博] 售SCI一區(qū)T0P文章，我:8.O.5.5.1.O.5.4,科目齊全,可+急 +3 | 7x4iz20zqa 2026-03-11 | 8/400 |
|
|---|---|---|---|---|
|
|
[考研] 工科0856專碩化學(xué)工程269能調(diào)劑嗎 +9 | 我想讀研11 2026-03-10 | 9/450 |
|
|
[考研] 293求調(diào)劑，一志愿陜師大生物學(xué) +3 | ??????.?.??? 2026-03-09 | 3/150 |
|
|
|
[考研] 0856材料與化工353分求調(diào)劑 +11 | NIFFFfff 2026-03-09 | 11/550 |
|
|
|
[考研] 0860求調(diào)劑（272分） +3 | lllllcsjsj 2026-03-05 | 4/200 |
|
|
|
[碩博家園] 木蟲好像不熱鬧了，是不是？ +4 | 偏振片 2026-03-10 | 4/200 |
|
|
|
[考研] 考研一志愿長(zhǎng)安大學(xué)材料與化工309分請(qǐng)求調(diào)劑 +6 | dtdxzxx 2026-03-06 | 8/400 |
|
|
|
[考研] 家人們 調(diào)劑不迷路 看這里 +8 | likeihood 2026-03-09 | 13/650 |
|
|
|
[考研] 294 英二數(shù)二物化 求調(diào)劑 +6 | 米飯團(tuán)不好吃 2026-03-09 | 6/300 |
|
|
|
[考研] 一志愿南大化學(xué)339分求調(diào)劑，四六級(jí)已過，有比賽，有文章 +7 | Gallantzhou 2026-03-07 | 7/350 |
|
|
|
[考研] 307求調(diào)劑 +3 | 辛仁豆腐 2026-03-08 | 5/250 |
|
|
|
[考研] 296求調(diào)劑 +4 | Xinyu Wu311 2026-03-09 | 4/200 |
|
|
|
[考研] 334求調(diào)劑 +8 | Trying] 2026-03-06 | 8/400 |
|
|
|
[考研] 346分材料求調(diào)劑 +5 | snow_反季節(jié)版 2026-03-07 | 5/250 |
|
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)0856材料與化工，前三科206，總分283，求調(diào)劑 +7 | 26考研求調(diào)劑 2026-03-06 | 7/350 |
|
|
|
[考研] 求調(diào)劑推薦 +4 | 微辣不吃 2026-03-06 | 4/200 |
|
|
|
[考研] 哈爾濱理工大學(xué)2026年研究生調(diào)劑，材料科學(xué)與化學(xué)工程學(xué)院研究生調(diào)劑 +3 | xinliu866 2026-03-06 | 3/150 |
|
|
|
[考博] 2026申博自薦 六級(jí)440電催化方向 +4 | 櫻落成影花成雙 2026-03-05 | 4/200 |
|
|
|
[考研] 334求調(diào)劑 +6 | Trying] 2026-03-05 | 8/400 |
|
|
|
[考研] 復(fù)試調(diào)劑 +5 | 呼呼？~+123456 2026-03-05 | 5/250 |
|