| 2 | 1/1 | 返回列表 |
| 查看: 1685 | 回復(fù): 1 | ||
liuchenhan新蟲 (小有名氣)
|
[求助]
關(guān)于石墨加壓的lammps模擬
|
|
小弟希望做lammps加壓下石墨晶格常數(shù)的變化。 要求:c-axis 方向加1GPa,面內(nèi)不加壓,自由變化,溫度要設(shè)在300K.以下是我的腳本。 # simulation # Initialization variable dt equal 0.0005 variable Tstart equal 1000.0 variable Tstop equal 300.0 variable Tdamp equal 100.0*${dt} variable Pstart equal 10.0 variable Pstop equal 10.0 variable Pdamp equal 1000.0*${dt} units metal dimension 3 newton on boundary p p p atom_style atomic neighbor 0.3 nsq neigh_modify check yes read_data nanonew.xyz # tersoff pair_style hybrid tersoff lj/cut 8.5 pair_coeff * * tersoff Cnew.tersoff C C C C C C C C pair_coeff 1 2 lj/cut 0.0024 3.4 #8.5 3.4 pair_coeff 1 3 lj/cut 0.0024 3.4 pair_coeff 1 4 lj/cut 0.0024 3.4 pair_coeff 1 5 lj/cut 0.0024 3.4 pair_coeff 1 6 lj/cut 0.0024 3.4 pair_coeff 1 7 lj/cut 0.0024 3.4 pair_coeff 1 8 lj/cut 0.0024 3.4 pair_coeff 2 3 lj/cut 0.0024 3.4 pair_coeff 2 4 lj/cut 0.0024 3.4 pair_coeff 2 5 lj/cut 0.0024 3.4 pair_coeff 2 6 lj/cut 0.0024 3.4 pair_coeff 2 7 lj/cut 0.0024 3.4 pair_coeff 2 8 lj/cut 0.0024 3.4 pair_coeff 3 4 lj/cut 0.0024 3.4 pair_coeff 3 5 lj/cut 0.0024 3.4 pair_coeff 3 6 lj/cut 0.0024 3.4 pair_coeff 3 7 lj/cut 0.0024 3.4 pair_coeff 3 8 lj/cut 0.0024 3.4 pair_coeff 4 5 lj/cut 0.0024 3.4 pair_coeff 4 6 lj/cut 0.0024 3.4 pair_coeff 4 7 lj/cut 0.0024 3.4 pair_coeff 4 8 lj/cut 0.0024 3.4 pair_coeff 5 6 lj/cut 0.0024 3.4 pair_coeff 5 7 lj/cut 0.0024 3.4 pair_coeff 5 8 lj/cut 0.0024 3.4 pair_coeff 6 7 lj/cut 0.0024 3.4 pair_coeff 6 8 lj/cut 0.0024 3.4 pair_coeff 7 8 lj/cut 0.0024 3.4 velocity all create ${Tstart} 1000001 mom yes rot yes dist gaussian fix NPT all npt temp ${Tstart} ${Tstop} ${Tdamp} x 10.0 10.0 1000.0*${dt} y 10.0 10.0 1000.0*${dt} z 1000.0 1000.0 1000.0*${dt} variable V equal vol variable LX equal lx variable LY equal ly variable LZ equal lz variable T equal temp variable PE equal pe variable E equal etotal fix VOL all print 100 "$V ${LX} ${LY} ${LZ}" file volume.txt screen no fix TE all print 100 "$T ${PE} $E" file TE.txt screen no run 1000000 問題:1、周期性邊界條件加的對不對,不對的話應(yīng)該怎么改; 2、壓強(qiáng)加的對不對,不對的話應(yīng)該怎么改;加壓強(qiáng)時需要周期性邊界條件嗎? 3、腳本還有哪些錯誤的地方,麻煩指出 非常感謝大神能為小弟指點迷津,剛剛?cè)腴T,還望海涵,拜托啦。 |
新蟲 (小有名氣)
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑323材料與化工 +4 | 1124361 2026-03-24 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 278求調(diào)劑 +4 | 我可以上岸的對?/a> 2026-03-19 | 4/200 |
|
|
[考研] 287求調(diào)劑 +10 | 晨昏線與星海 2026-03-19 | 11/550 |
|
|
[考研] B區(qū)考研調(diào)劑 +4 | yqdszhdap- 2026-03-22 | 5/250 |
|
|
[考研] 08工學(xué)調(diào)劑 +12 | 用戶573181 2026-03-20 | 17/850 |
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
[考研]
|
黃粱一夢千年 2026-03-24 | 3/150 |
|
|
[考研] 341求調(diào)劑(一志愿湖南大學(xué)070300) +5 | 番茄頭--- 2026-03-22 | 6/300 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境,總分308求調(diào)劑 +7 | 墨墨漠 2026-03-23 | 8/400 |
|
|
[考研] 328求調(diào)劑 +4 | LHHL66 2026-03-23 | 4/200 |
|
|
[考研] 材料與化工085600,總分304,本科有兩篇sci參與,求調(diào)劑 +4 | 幸運的醬醬 2026-03-22 | 5/250 |
|
|
[考研] 求調(diào)劑一志愿海大,0703化學(xué)學(xué)碩304分,有大創(chuàng)項目,四級已過 +6 | 幸運哩哩 2026-03-22 | 10/500 |
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 4/200 |
|
|
[考研] 303求調(diào)劑 +5 | 安憶靈 2026-03-22 | 6/300 |
|
|
[考研]
|
Grand777 2026-03-21 | 3/150 |
|
|
[考研] 311求調(diào)劑 +3 | 勇敢的小吳 2026-03-20 | 3/150 |
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | 白QF 2026-03-21 | 3/150 |
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
[考研] 一志愿福大288有機(jī)化學(xué),求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|