版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗(yàn)證，感謝您的理解與支持！

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

jmforever123

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 198.5
散金: 10
帖子: 44
在線: 16.2小時
蟲號: 3121805
注冊: 2014-04-08
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

[求助] 石墨烯上的水滴的模擬已有1人參與

本人模擬的是石墨烯表面的水分子的運(yùn)動情況，但用vmd顯示后，石墨烯片總是破壞掉了，感覺z方向要使用非周期性邊界條件，但是水分子有長程力作業(yè)，pppm又必須使用周期性邊界條件，本人菜鳥，求大神指教。謝謝。
dimension    3
boundary    p p p
atom_style    full
read_data    h2o.data
pair_style    hybrid lj/cut/tip4p/long 2 1 1 1 0.1546 10.0 12.0 lj/cut 10.0  tersoff
pair_coeff    1 1 lj/cut/tip4p/long 0.0000 0.0000
pair_coeff    1 2 lj/cut/tip4p/long 0.0000 0.0000
pair_coeff    2 2 lj/cut/tip4p/long 0.16275 3.16435
pair_coeff    1 3 lj/cut 0.0000 0.0000
pair_coeff    2 3 lj/cut 0.3876 3.262
pair_coeff    * * tersoff SiC.tersoff NULL NULL C
kspace_style pppm/tip4p 1.0e-4

bond_style    harmonic
bond_coeff    1 450 0.9572

angle_style harmonic
angle_coeff 1 55 104.52

timestep    ${dt}
thermo       $d

graphene.png

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

材料	分子動力學(xué)模擬

» 猜你喜歡

346求調(diào)劑一志愿070303有機(jī)化學(xué) 已經(jīng)有9人回復(fù)
085601 329分調(diào)劑已經(jīng)有5人回復(fù)
289求調(diào)劑已經(jīng)有4人回復(fù)
367求調(diào)劑已經(jīng)有6人回復(fù)
271求調(diào)劑已經(jīng)有4人回復(fù)
085600，321分求調(diào)劑已經(jīng)有3人回復(fù)
化學(xué)工程085602 305分求調(diào)劑已經(jīng)有35人回復(fù)
282求調(diào)劑不挑專業(yè) 求收留已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有7人回復(fù)
334分一志愿武理-080500 材料求調(diào)劑已經(jīng)有6人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

1樓 2016-06-13 16:46:31

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

完全是火

銅蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 67.1
散金: 16
紅花: 4
帖子: 275
在線: 117.6小時
蟲號: 2742436
注冊: 2013-10-21
專業(yè): 多相流熱物理學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
jmforever123: 金幣+5, ★★★很有幫助 2016-06-14 15:19:55

這里有個例子，你參考一下能用不

units    metal
dimension  3
boundary p p f
atom_style full

read_data    mixed-graph-water.dat
neighbor    3.0 bin
neigh_modify  delay 0 every 1 check yes

# set potential function and parameters

### Create Groups ###
group hy type 2    # H_water
group ox type 3    # O_water
set group ox charge -0.8476 #給氫氧原子帶電
set group hy charge 0.4238
region box block INF INF INF INF INF INF units box
group graphite type 1
group boxes region box

### Flexible SPC/E Potential Parameters ###
### Zhang et al., Fluid Phase Equilibria, 262 (2007) 210-216 ###
### 設(shè)置氫氧原子間的作用
pair_style hybrid/overlay rebo lj/cut/coul/long 10.0
pair_coeff * * rebo ./potential/CH.airebo C NULL NULL
pair_coeff 1 2 lj/cut/coul/long 0.0 0.0       # C-H
pair_coeff 1 3 lj/cut/coul/long 0.013 2.71       # C-O
pair_coeff 2 2 lj/cut/coul/long 0.0 0.0       # H-H
pair_coeff 2 3 lj/cut/coul/long 0.0 0.0       # H-O
pair_coeff 3 3 lj/cut/coul/long 0.0067 3.166    # O-O

bond_style  harmonic
bond_coeff  1 7.669 1.0
angle_style harmonic
angle_coeff 1 1.829 109.47
kspace_style  pppm 1.0e-5 #final npt relaxation
kspace_modify slab 3.0

thermo       1
thermo_style custom step pe etotal press temp
thermo_modify norm no
dump          myatom all custom 500 atom.dat id

group slab    type 1    # graphite_slab
group water    type 2 3

timestep 0.0001

velocity    all create 10.0 482748 dist gaussian
fix       freaze_slab slab move linear 0.0 0.0 0.0 #This can be useful for boundary or other atoms, whose movement can influence nearby atoms.
#                            move style  args
fix       zwall water wall/reflect  zhi EDGE
#                      wall/reflect face arg ... keyword value ...
fix       1 all nve
fix       2 all temp/rescale 10 10.0 10.0 0.02 1.0          ##溫度設(shè)置為10度
#    temp/rescale N  Tstart Tstop  window fraction
fix       shake_1 water shake 1e-5 500 0    b       1       a    1 # 算法固定鍵長鍵角化學(xué)鍵的斷裂、生成
dump 1 graphite custom  1  graphite.lammpstrj id type xs ys zs
#                         style tol iter N constraint values  constraint values
dump       dcd_1 all dcd 100 nve_10.dcd
run 10000
write_restart restart.10.*
unfix       2
unfix       shake_1
undump    dcd_1
fix       3 all temp/rescale 10 10.0 100.0 0.02 1.0
fix       shake_2 water shake 1e-5 500 0 b 1 a 1
dump       dcd_2 all dcd 100 nve_10_100.dcd
run 10000
write_restart restart.10_100.*
unfix       3
unfix       shake_2
undump    dcd_2
fix       4 all temp/rescale 10 100.0 200.0 0.02 1.0
fix       shake_3 water shake 1e-5 500 0 b 1 a 1
fix       shake_3 water shake 1e-5 500 0 b 1 a 1
dump       dcd_3 all dcd 100  nve_100_200.dcd
run 10000
write_restart  restart.100_200.*
unfix       4
unfix       shake_3
undump    dcd_3
fix       5 all temp/rescale 10 200.0 298.0 0.02 1.0
fix       shake_4 water shake 1e-5 500 0 b 1 a 1
dump       dcd_4 all dcd 100  nve_200_298.dcd
run 10000
write_restart  restart.200_298.*
unfix       5
unfix       shake_4
undump    dcd_4
fix       6 all temp/rescale 10 298.0 298.0 0.02 1.0
fix       shake_5 water shake 1e-5 500 0 b 1 a 1
dump       dcd_5 all dcd 100  nve_298_298.dcd
run 10000
write_restart restart.298.*

贊一下(3人)

回復(fù)此樓

2樓2016-06-13 19:54:50

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

568286346

木蟲 (小有名氣)

應(yīng)助: 31 (小學(xué)生)
金幣: 789
帖子: 288
在線: 153.3小時
蟲號: 1785423
注冊: 2012-04-29
專業(yè): 理論和計算化學(xué)

我覺得你可以把z方向的盒子邊長擴(kuò)大，相當(dāng)于在z方向構(gòu)建一個汽液界面。這樣就相當(dāng)于在三個方向都使用了周期性邊界條件。

發(fā)自小木蟲IOS客戶端

贊一下(2人)

回復(fù)此樓

堅持才能實(shí)現(xiàn)自己的理想！

5樓2016-06-14 19:15:13

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

jmforever123

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 198.5
散金: 10
帖子: 44
在線: 16.2小時
蟲號: 3121805
注冊: 2014-04-08
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

引用回帖:

2樓: Originally posted by 完全是火 at 2016-06-13 19:54:50
這里有個例子，你參考一下能用不

units    metal
dimension  3
boundary p p f
atom_style full

read_data    mixed-graph-water.dat
neighbor    3.0 bin
neigh_modify  delay 0 every 1 check  ...

我試一下，非常感謝

贊一下

回復(fù)此樓

3樓2016-06-13 20:06:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jmforever123

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 198.5
散金: 10
帖子: 44
在線: 16.2小時
蟲號: 3121805
注冊: 2014-04-08
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

引用回帖:

還是不行，老是出現(xiàn)ERROR：Lost atoms: original 1520 current 1073

贊一下

回復(fù)此樓

4樓2016-06-14 15:17:10

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jmforever123

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 198.5
散金: 10
帖子: 44
在線: 16.2小時
蟲號: 3121805
注冊: 2014-04-08
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

引用回帖:

5樓: Originally posted by 568286346 at 2016-06-14 19:15:13
我覺得你可以把z方向的盒子邊長擴(kuò)大，相當(dāng)于在z方向構(gòu)建一個汽液界面。這樣就相當(dāng)于在三個方向都使用了周期性邊界條件。

非常感謝，我這個的主要原因就是z方向是不對稱的對吧，然后就需要把z方向擴(kuò)大�？墒菙U(kuò)大之后，下邊界是C原子，上邊界沒有原子，是不是也需要把下邊界往下擴(kuò)展，不知道可以嗎

贊一下

回復(fù)此樓

6樓2016-06-14 22:25:30

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

568286346

木蟲 (小有名氣)

應(yīng)助: 31 (小學(xué)生)
金幣: 789
帖子: 288
在線: 153.3小時
蟲號: 1785423
注冊: 2012-04-29
專業(yè): 理論和計算化學(xué)

可以的，因?yàn)楦鶕?jù)周期性邊界條件，上方擴(kuò)充之后，下方的C其實(shí)也是處于汽固界面的。這其實(shí)跟你擴(kuò)充下方是一樣的。

發(fā)自小木蟲IOS客戶端

贊一下

回復(fù)此樓

堅持才能實(shí)現(xiàn)自己的理想！

7樓2016-06-15 00:08:32

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jmforever123

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 198.5
散金: 10
帖子: 44
在線: 16.2小時
蟲號: 3121805
注冊: 2014-04-08
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

引用回帖:

7樓: Originally posted by 568286346 at 2016-06-15 00:08:32
可以的，因?yàn)楦鶕?jù)周期性邊界條件，上方擴(kuò)充之后，下方的C其實(shí)也是處于汽固界面的。這其實(shí)跟你擴(kuò)充下方是一樣的。

謝謝

回復(fù)此樓

8樓2016-06-15 11:09:10

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

Niklaus

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1339.5
帖子: 155
在線: 87小時
蟲號: 2795553
注冊: 2013-11-12
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

引用回帖:

老板，一般怎么找代碼哈？

發(fā)自小木蟲IOS客戶端

贊一下

回復(fù)此樓

高溫合金狗一只，求各路大神指點(diǎn)

9樓2016-06-15 15:02:00

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

完全是火

銅蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 67.1
散金: 16
紅花: 4
帖子: 275
在線: 117.6小時
蟲號: 2742436
注冊: 2013-10-21
專業(yè): 多相流熱物理學(xué)

引用回帖:

4樓: Originally posted by jmforever123 at 2016-06-14 15:17:10
還是不行，老是出現(xiàn)ERROR：Lost atoms: original 1520 current 1073...

不做這方面的幫不了你哈，好好研究吧。

贊一下

回復(fù)此樓

10樓2016-06-16 14:36:50

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 jmforever123 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 289求調(diào)劑 +3	Acesczlo 2026-03-29	4/200	2026-03-31 14:48 by 熱情沙漠
[考研] 271求調(diào)劑 +3	勒布朗@ 2026-03-31	4/200	2026-03-31 14:44 by 勒布朗@
[考研] 0710生物學(xué)求調(diào)劑！ +4	敘述文 2026-03-31	4/200	2026-03-31 14:32 by 記事本2026
[考研] 317求調(diào)劑 +11	蛋黃咸肉粽 2026-03-26	12/600	2026-03-31 13:30 by 記事本2026
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +9	柳溪邊 2026-03-26	9/450	2026-03-31 12:16 by oooqiao
[考研] 269求調(diào)劑 +4	我想讀研11 2026-03-31	4/200	2026-03-31 10:04 by cal0306
[考研] 看看我！301求調(diào)劑 081700 +3	一個小明同學(xué) 2026-03-24	3/150	2026-03-31 09:40 by longlotian
[考研] 各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩 +10	Koxui 2026-03-28	10/500	2026-03-30 19:33 by 源_2020
[考研] 311求調(diào)劑 +7	藍(lán)月亮亮 2026-03-30	7/350	2026-03-30 19:24 by 無際的草原
[考研] 一志愿廈門大學(xué)材料工程專碩354找調(diào)劑�。。� +5	貝唄鋇鋇 2026-03-30	5/250	2026-03-30 18:16 by 無際的草原
[考研] 342求調(diào)劑 +4	加油a李zs 2026-03-26	4/200	2026-03-30 16:39 by 晶體之美
[考研] 311求調(diào)劑 +10	lin0039 2026-03-26	10/500	2026-03-30 10:26 by herarysara
[考研] 337求調(diào)劑 +6	《樹》 2026-03-29	6/300	2026-03-30 10:15 by herarysara
[考研] 一志愿南航 335分 \| 0856 \| GPA 4.07 \| 有科研經(jīng)歷 +8	cccchenso 2026-03-29	8/400	2026-03-29 23:53 by 我是小康
[考研] 343求調(diào)劑085601 +3	要努力學(xué)習(xí)x 2026-03-29	3/150	2026-03-29 18:35 by wxiongid
[考研] 327求調(diào)劑 +6	汲亦昊 2026-03-29	6/300	2026-03-29 13:40 by peike
[考研] 311（085601）求調(diào)劑 +4	liziyeyeye 2026-03-28	4/200	2026-03-28 18:50 by 535743368
[考研] 285求調(diào)劑 +4	AZMK 2026-03-27	7/350	2026-03-27 20:59 by AZMK
[考研] 0856調(diào)劑 +5	求求讓我有書讀?/a> 2026-03-26	6/300	2026-03-27 15:12 by caszguilin
[考研] 324求調(diào)劑 +5	hanamiko 2026-03-26	5/250	2026-03-27 10:33 by wangjy2002

24小時熱門版塊排行榜

[求助] 石墨烯上的水滴的模擬 已有1人參與

» 收錄本帖的淘帖專輯推薦

» 猜你喜歡

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

【答案】應(yīng)助回帖

[求助] 石墨烯上的水滴的模擬已有1人參與