 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 1210 | 回復(fù): 1 | |||
Lionzaka新蟲 (初入文壇)
|
[交流]
小白提問 關(guān)于單位轉(zhuǎn)換
|
|
剛剛接觸LAMMPS，本身物理基礎(chǔ)也不是很好，不太會設(shè)置一些參數(shù)。本人主要使用LAMMPS研究聚電解質(zhì)溶液。 選用LJ相互作用勢，將LJ表達(dá)式中ε作為能量單位；將庫倫能作為靜電相互作用，庫倫強度系數(shù)Γ=lb/l0，lb是Bjerrum長度，l0是長度單位。Γ取值在0.1-7.0之間。在使用郎之萬動力學(xué)進(jìn)行模擬時，需提供溫度，想知道如何計算溫度。 如果按照能量單位來說，ε=kbT，那這里T豈不是1/kb? 以及在Bjerrum長度中有真空介電常數(shù)，這個要如何轉(zhuǎn)換為單位為lj下的數(shù)值？ 下面是in文件。 # BD simulation of polyelectrolyte polymer in counterion ion echo screen units lj # 單位 dimension 3 # BOX維度 #newton on # 默認(rèn)是on，體系特別大的時候可以off boundary p p p # 周期性邊界 atom_style full # 原子類型 重點需要搞清楚 bond_style harmonic # 設(shè)置鍵、角相互作用 angle_style harmonic read_data data.file # 讀取data文件 pair_style lj/cut/coul/long 1.12246 8.0 # 設(shè)置對相互作用 pair_coeff 1 1 1.0 1.0 pair_coeff 1 2 1.0 0.8 0.898 # 0.8*1.12246 pair_coeff 2 2 1.0 0.6 0.673 # 0.6*1.12246 pair_modify shift yes #保證能量在截斷點時為0 #pair_style soft 1.12246 #pair_coeff * * 10.0 bond_coeff 1 5000.0 1.0 #鍵、角的系數(shù) angle_coeff 1 50.0 180.0 kspace_style pppm 1.0e-4 kspace_modify order 6 ########################################### # 能量最小化 minimize 1.0e-4 1.0e-6 10000 100000 # 寫fix語句，實現(xiàn)分子動力學(xué)模擬-郎之萬動力學(xué) neighbor 0.3 bin #neigh_modify every 2 delay 10 check yes page 100000 #delay 10 代表隔10步 velocity all create 1.2 4928459 rot yes dist gaussian fix 10 all nve fix 20 all langevin 1.2 1.2 1000.0 699483 #zero yes thermo 1000 thermo_style custom step temp pe etotal press vol # 弛豫過程 dump 1 all atom 1000 dump1.lammpstrj #收集數(shù)據(jù),隔1000步收集一個 run 10000 # data collection undump 1 dump 2 all custom 1000 dump2.lammpstrj id type mol x y z vx vy vz fx fy fz run 100000 |
回復(fù)此樓» 猜你喜歡
070300調(diào)劑一志愿南京大學(xué)化學(xué)專業(yè)，初試273分，有論文有專利有競賽，求調(diào)劑B區(qū)211
已經(jīng)有8人回復(fù)
-
322分 085600求調(diào)劑，有互聯(lián)網(wǎng)+國金及主持省級大創(chuàng)經(jīng)歷
已經(jīng)有4人回復(fù)
-
320材料一志愿西工大專碩英二數(shù)二 有兩年光伏工作經(jīng)驗
已經(jīng)有3人回復(fù)
-
一志愿鄭州大學(xué)，學(xué)碩，物理化學(xué)， 333求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿中石油，化學(xué)工程，301求調(diào)劑
已經(jīng)有3人回復(fù)
-
0703化學(xué)求調(diào)劑，總分320分，一志愿中南民族大學(xué)
已經(jīng)有3人回復(fù)
-
320材料與化工，求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿武漢理工大學(xué)-085602-總分296分-求調(diào)劑
已經(jīng)有5人回復(fù)
-
286
已經(jīng)有3人回復(fù)
-
070300化學(xué) 280 一志愿太原理工 求調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
Lionzaka
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 20
- 帖子: 2
- 在線: 7.7小時
- 蟲號: 15121290
- 注冊: 2019-05-17
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
2樓2023-03-02 17:33:38
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 Lionzaka 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703 東華大學(xué) 理學(xué)化學(xué)均過 a 線，總分 281 求調(diào)劑 +3 | 香香 Lu 2026-03-03 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 290求調(diào)劑 +3 | Bananaiy 2026-03-04 | 3/150 |
|
|
|
[考研] 一志愿東北大學(xué)080500，英一數(shù)二343分求調(diào)劑 +3 | Seele_v5 2026-03-04 | 3/150 |
|
|
|
[考研] 304分材料專碩求調(diào)劑 +8 | qiuzhigril 2026-03-03 | 10/500 |
|
|
|
[論文投稿]
 主編終審卡者不動
 100+3
|
Stray2021 2026-03-03 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | G212 2026-03-03 | 5/250 |
|
|
|
[考研] 求調(diào)劑院校 +6 | 云朵452 2026-03-02 | 11/550 |
|
|
|
[考研] 材料工程269求調(diào)劑 +7 | 白刺玫 2026-03-02 | 7/350 |
|
|
|
[考研] 接收調(diào)劑 +15 | 津萌津萌 2026-03-02 | 23/1150 |
|
|
|
[考研] 環(huán)境工程學(xué)碩288求助調(diào)劑 +7 | 多吃億口芝士 2026-03-02 | 7/350 |
|
|
|
[考研] 291求調(diào)劑 +4 | Afy123456 2026-03-03 | 7/350 |
|
|
|
[考研] 276求調(diào)劑 +8 | 路lyh123 2026-02-28 | 10/500 |
|
|
|
[基金申請] 沒有青基直接申請面上，感覺自己瘋了 +5 | kevin63t 2026-03-02 | 6/300 |
|
|
|
[考研] 0805總分292，求調(diào)劑 +12 | 幻想之殤 2026-03-01 | 12/600 |
|
|
|
[考研] 一志愿中科大能動297求調(diào)劑，本科川大 +3 | 邵11 2026-03-03 | 3/150 |
|
|
|
[考研] 0854復(fù)試調(diào)劑 276 +5 | wmm9 2026-03-01 | 7/350 |
|
|
|
[考研] 0856材料求調(diào)劑 +12 | hyf hyf hyf 2026-02-28 | 13/650 |
|
|
|
[考研] 0856化工專碩求調(diào)劑 +15 | 董boxing 2026-03-01 | 15/750 |
|
|
|
[考研] 材料類求調(diào)劑 +11 | wana_kiko 2026-02-28 | 14/700 |
|
|
|
[論文投稿]
Optics letters投稿被拒求助
30+3
|
luckyry 2026-02-26 | 4/200 |
|