 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 1686 | 回復: 1 | ||
liuchenhan新蟲 (小有名氣)
|
[求助]
關于石墨加壓的lammps模擬
|
|
小弟希望做lammps加壓下石墨晶格常數的變化。 要求：c-axis 方向加1GPa，面內不加壓，自由變化，溫度要設在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、周期性邊界條件加的對不對，不對的話應該怎么改; 2、壓強加的對不對，不對的話應該怎么改；加壓強時需要周期性邊界條件嗎？ 3、腳本還有哪些錯誤的地方，麻煩指出 非常感謝大神能為小弟指點迷津，剛剛入門，還望海涵，拜托啦。 |
回復此樓» 猜你喜歡
北科281學碩材料求調劑
已經有13人回復
-
289求調劑
已經有11人回復
-
網絡空間安全0839招調劑
已經有5人回復
-
材料277求調劑
已經有4人回復
-
274求調劑
已經有7人回復
-
085600 材料與化工 329分求調劑
已經有3人回復
-
302求調劑
已經有3人回復
-
312求調劑
已經有5人回復
-
329求調劑
已經有3人回復
-
【2026考研調劑】制藥工程 284分 求相關專業(yè)調劑名額
已經有8人回復
» 本主題相關商家推薦: (我也要在這里推廣)
» 本主題相關價值貼推薦，對您同樣有幫助:
- MS修改鍵長，MS畫石墨烯，MS文件導入LAMMPS
已經有3人回復
- 石墨烯 lammps
已經有11人回復
liuchenhan
新蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 1536.5
- 散金: 4
- 帖子: 170
- 在線: 108小時
- 蟲號: 1508714
- 注冊: 2011-11-25
- 專業(yè): 傳熱傳質學
2樓2014-02-25 17:38:48
相關版塊跳轉
我要訂閱樓主 liuchenhan 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 334分 一志愿武理-080500 材料求調劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 求調劑，一志愿:南京航空航天大學大學 ，080500材料科學與工程學碩，總分289分 +6 | @taotao 2026-03-19 | 6/300 |
|
|
|
[考研] 289材料與化工（085600）B區(qū)求調劑 +4 | 這么名字咋樣 2026-03-22 | 5/250 |
|
|
|
[考研] 一志愿華東理工大學081700，初試分數271 +5 | kotoko_ik 2026-03-23 | 6/300 |
|
|
|
[考研] 環(huán)境學碩288求調劑 +8 | 皮皮皮123456 2026-03-22 | 8/400 |
|
|
|
[考研] 一志愿北京化工大學 070300 學碩 336分 求調劑 +7 | vv迷 2026-03-22 | 7/350 |
|
|
|
[考研] 327求調劑 +5 | prayer13 2026-03-23 | 5/250 |
|
|
|
[考研] 361求調劑 +3 | Glack 2026-03-22 | 3/150 |
|
|
|
[考研] 298求調劑 +8 | 上岸6666@ 2026-03-20 | 8/400 |
|
|
|
[考研] 070300，一志愿北航320求調劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|
|
|
[考研] 317求調劑 +12 | 申子申申 2026-03-19 | 18/900 |
|
|
|
[考研] 一志愿華中農業(yè)071010，總分320求調劑 +5 | 困困困困坤坤 2026-03-20 | 6/300 |
|
|
|
[考研] 求調劑院校信息 +6 | CX 330 2026-03-21 | 6/300 |
|
|
|
[考研] 269專碩求調劑 +6 | 金恩貝 2026-03-21 | 6/300 |
|
|
|
[考研] 085600材料與化工306 +4 | z1z2z3879 2026-03-21 | 4/200 |
|
|
|
[考研] 初試 317 +7 | 半拉月丙 2026-03-20 | 7/350 |
|
|
|
[考研] 311求調劑 +3 | 勇敢的小吳 2026-03-20 | 3/150 |
|
|
|
[考研] 材料與化工（0856）304求 B區(qū) 調劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
|
[考研] 288求調劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
|
[考研] 材料考研調劑 +3 | xwt。 2026-03-19 | 3/150 |
|