 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 1209 | 回復(fù): 1 | |||
Lionzaka新蟲 (初入文壇)
|
[交流]
小白提問 關(guān)于單位轉(zhuǎn)換
|
|
剛剛接觸LAMMPS，本身物理基礎(chǔ)也不是很好，不太會(huì)設(shè)置一些參數(shù)。本人主要使用LAMMPS研究聚電解質(zhì)溶液。 選用LJ相互作用勢(shì)，將LJ表達(dá)式中ε作為能量單位；將庫(kù)倫能作為靜電相互作用，庫(kù)倫強(qiáng)度系數(shù)Γ=lb/l0，lb是Bjerrum長(zhǎng)度，l0是長(zhǎng)度單位。Γ取值在0.1-7.0之間。在使用郎之萬動(dòng)力學(xué)進(jìn)行模擬時(shí)，需提供溫度，想知道如何計(jì)算溫度。 如果按照能量單位來說，ε=kbT，那這里T豈不是1/kb? 以及在Bjerrum長(zhǎng)度中有真空介電常數(shù)，這個(gè)要如何轉(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，體系特別大的時(shí)候可以off boundary p p p # 周期性邊界 atom_style full # 原子類型 重點(diǎn)需要搞清楚 bond_style harmonic # 設(shè)置鍵、角相互作用 angle_style harmonic read_data data.file # 讀取data文件 pair_style lj/cut/coul/long 1.12246 8.0 # 設(shè)置對(duì)相互作用 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 #保證能量在截?cái)帱c(diǎn)時(shí)為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語(yǔ)句，實(shí)現(xiàn)分子動(dòng)力學(xué)模擬-郎之萬動(dòng)力學(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步收集一個(gè) 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ù)此樓» 猜你喜歡
0856材料與化工，270求調(diào)劑
已經(jīng)有15人回復(fù)
-
085600 材料與化工 298
已經(jīng)有12人回復(fù)
-
312求調(diào)劑
已經(jīng)有7人回復(fù)
-
304分材料專碩求調(diào)劑
已經(jīng)有9人回復(fù)
-
320材料與化工，求調(diào)劑
已經(jīng)有5人回復(fù)
-
成果系統(tǒng)訪問量大，請(qǐng)15分鐘后再嘗試。由此給您造成的不便，敬請(qǐng)諒解。
已經(jīng)有12人回復(fù)
-
281求調(diào)劑
已經(jīng)有5人回復(fù)
-
347求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿314求調(diào)劑
已經(jīng)有6人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
Lionzaka
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 20
- 帖子: 2
- 在線: 7.7小時(shí)
- 蟲號(hào): 15121290
- 注冊(cè): 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ā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿314求調(diào)劑 +6 | 202111120625 2026-03-03 | 6/300 |
|
|---|---|---|---|---|
|
[考博] 26超級(jí)電容器申博 +4 | dhdjdjend 2026-02-25 | 4/200 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | G212 2026-03-03 | 5/250 |
|
|
|
[論文投稿]
EST拒稿重投
 5+3
|
15102603076 2026-03-02 | 3/150 |
|
|
|
[考研] 復(fù)試調(diào)劑 +4 | 杰尼龜aaa 2026-03-03 | 4/200 |
|
|
|
[考研] 298求調(diào)劑一志愿中海洋 +3 | lour. 2026-03-03 | 3/150 |
|
|
|
[考研] 289求調(diào)劑 +7 | BrightLL 2026-03-02 | 9/450 |
|
|
|
[考研] 0805總分292，求調(diào)劑 +12 | 幻想之殤 2026-03-01 | 12/600 |
|
|
|
[考研] 環(huán)境調(diào)劑 +5 | 柒槿levana 2026-03-01 | 5/250 |
|
|
|
[考研] 化工270求調(diào)劑 +10 | 什么名字qwq 2026-03-02 | 10/500 |
|
|
|
[考研] 278求調(diào)劑 +3 | 滿天星11_22 2026-03-02 | 3/150 |
|
|
|
[考研] 沒上岸的看過來 +3 | tangxiaotian 2026-03-01 | 3/150 |
|
|
|
[考研] 材料工程求調(diào)劑 +3 | 1431251 2026-03-03 | 3/150 |
|
|
|
[考研] 0857調(diào)劑 +6 | 一ll半 2026-02-28 | 8/400 |
|
|
|
[考研] 課題組接收材料類調(diào)劑研究生 +6 | gaoxiaoniuma 2026-02-28 | 9/450 |
|
|
|
[考研] 307求調(diào)劑 +6 | wyyyqx 2026-03-01 | 6/300 |
|
|
|
[考研] 271求調(diào)劑 +4 | Ricardo1113 2026-03-02 | 4/200 |
|
|
|
[考研] 338求調(diào)劑 +5 | 18162027187 2026-03-02 | 6/300 |
|
|
|
[考研] 【2026 碩士調(diào)劑】課題組 招收調(diào)劑生 +3 | 考研版棒棒 2026-03-02 | 5/250 |
|
|
|
[考研] 317一志愿華南理工電氣工程求調(diào)劑 +6 | Soliloquy_Q 2026-02-28 | 11/550 |
|