 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 7267 | 回復: 11 | ||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
wangsong1016銅蟲 (初入文壇)
|
[求助]
fix deposit沉積原子被反彈
|
|
|
本人目前利用Lammps軟件進行薄膜生長模擬。目前參考文獻J. Phys. D: Appl. Phys. 42 (2009) 065107進行重復性工作。當然這個文獻中采用的是DL_Ploy軟件，我利用Lammps應該也能重復該工作。利用Lammps做薄膜生長很重要的一個命令式fix deposit命令，能夠以一定的時間步長為周期，向體系中周期性的加入一個原子。位置，角度，速度可控。計算采用的勢能及參數來源于上面的參考文獻，勢能采用Buckingham+Coulumb勢，氧離子間的相互作用采用純庫倫勢。該文獻研究的是MgO表面上生長MgAlO合金，我重復該工作的時候出現(xiàn)不合理的情況，請大家?guī)兔⒖家幌隆?br />
1.最基本的方面入手，我模擬MgO基底上生長MgO薄膜。從文件中讀取MS建模得到的MgO坐標。勢參量如上文獻設置。共四層原子，在第三層位置利用Langevin熱浴，維持一定溫度，上面兩層原子nve系綜自由弛豫。最底層原子固定不動。（物理上可能不盡合理，作為初步測試應該足夠） 2.在大概5Ang之外的位置設定為沉積區(qū)域，原子Mg和O在該區(qū)域產生并沉積下來。 3.方便起見，我目前都只在windows下利用單核進行測試，模型初步建立起來之后再轉移到Linux系統(tǒng)下并行計算。 模擬結果發(fā)現(xiàn)利用deposit產生的這個原子行為很古怪，原子在向表面入射的過程中，會被表面反彈回來，不能吸附在基底上。換一種沉積產生原子的方式：這個原子不是利用fix deposit命令產生，而是在讀取MgO原子坐標的時候，在MgO基底之外再多讀取一個Mg或者O的坐標，那么這個原子就能夠吸附到基底表面上。對這樣的結果感到很困惑，希望大家給予幫助。 輸入腳本如下，另外坐標數據見附件： #This script is used to simulate the MgO structure units metal #The time step is picoseconds atom_style charge boundary p p p pair_style buck/long/coul/long cut long 8.0 read_data MgO_structure.in region thermostatregion block 0.00 25.2672 0.00 25.2672 2.00 2.20 group thermostatlayer region thermostatregion region surfaceregion block 0.00 25.2672 0.00 25.2672 2.20 14.90 group surfacelayer region surfaceregion region ParticleSourceRegion block 0.00 25.2672 0.00 25.2672 15.00 15.20 group ParticleSourceGroup region ParticleSourceRegion mass 1 24 #Mg mass 2 16 # O set type 1 charge 2.0 set type 2 charge -2.0 pair_coeff 1 2 1279.69 0.29969 0.0 pair_coeff 2 2 9547.96 0.21916 32.0 pair_coeff 1 1 0.0 1.1 0.0 #this term calculate the coulombic interaction between two Mg2+. Parameter A and C are set to be zero, while rho can't be set to 0.0, so here a tiny value are added. kspace_style ewald 1.0e-4 min_style cg minimize 1.0e-4 0.01 1000 10000 compute temperature surfacelayer temp #unfixedsubstrate temp fix 1 thermostatlayer nve fix 2 thermostatlayer langevin 300 300 0.10 4876 fix 3 surfacelayer nve fix 4 ParticleSourceGroup nve region depositregion block 0.00 25.2672 0.00 25.2672 10.20 10.40 group depositatoms region depositregion fix 5 depositatoms nve fix 6 depositatoms deposit 1 2 500 9375 region depositregion vz -14.0 -14.0 velocity ParticleSourceGroup set 2.0 2.0 -1.0 thermo_style custom step atoms temp epair etotal press thermo 100 thermo_modify temp temperature lost/bond ignore lost warn dump 3 all movie 10 movie.mpg type type & axes yes 0.8 0.02 view 80 -30 dump 2 all xyz 100 MgO.xyz run 2000 模擬所需時間極短，希望大家能夠運行測試，多多幫忙，謝謝！ |
回復此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內容負責。
本內容由用戶自主發(fā)布，如果其內容涉及到知識產權問題，其責任在于用戶本人，如對版權有異議，請聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : MgO_structure.in
2014-07-10 09:30:34, 24.91 K
» 猜你喜歡
070300化學求調劑
已經有9人回復
-
085404 22408 315分
已經有6人回復
-
材料與化工272求調劑
已經有22人回復
-
347求調劑
已經有10人回復
-
354求調劑
已經有4人回復
-
一志愿華東理工大學，080500學碩，317分，求調劑
已經有4人回復
-
一志愿a區(qū)211，085601-307分求調劑
已經有11人回復
-
284求調劑
已經有7人回復
-
求調劑：085600材料與化工，考材科基，總分319
已經有18人回復
-
362求調劑
已經有7人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 原子層沉積制備TiO2三維納米結構及其光催化性能研究
已經有15人回復
- 求助有關原子層沉積用鉑源
已經有2人回復
- 原子層沉積復合薄膜厚退火的作用？誰知道��？
已經有10人回復
- 泵油可以反復過濾后重復使用嗎，原子層沉積系統(tǒng)的真空抽不了，已排除漏氣可能
已經有7人回復
- MS模擬中固定原子Fix cartesian Position和Fix fractional position 不同點
已經有7人回復
- ACCOUTS OF CHEMICAL RESEARCH:原子層沉積在制備和提高負載型金屬催化劑方面的應用
已經有63人回復
- 原子總數怎么看
已經有3人回復
- 關于Atomic Layer Deposition（原子層沉積技術）的一篇綜述
已經有322人回復
- 關于原子吸收分光光度法檢測沉積物和海水的檢出限問題
已經有7人回復
- 原子層沉積技術-ALD
已經有7人回復
shawn_zhang
新蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 3.5
- 帖子: 3
- 在線: 57分鐘
- 蟲號: 7389167
- 注冊: 2017-10-29
- 專業(yè): 微/納機械系統(tǒng)
8樓2018-05-22 10:30:07
wangsong1016
銅蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 169
- 散金: 4
- 帖子: 45
- 在線: 26.3小時
- 蟲號: 408613
- 注冊: 2007-06-21
- 性別: GG
- 專業(yè): 物理
2樓2014-07-10 10:56:19
wangsong1016
銅蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 169
- 散金: 4
- 帖子: 45
- 在線: 26.3小時
- 蟲號: 408613
- 注冊: 2007-06-21
- 性別: GG
- 專業(yè): 物理
|
無論利用set type 1 charge 1這種命令來設置原子電荷，還是利用read_data命令從文件中讀取電荷，都只能對體系中當前存在的原子進行電荷賦值，而在這兩個命令之后，再利用fix deposit命令產生的原子還是不具有電荷。因此每一次fix deposit命令之后，run一次，產生了一個原子之后立刻用set命令進行電荷賦值，這樣原子就具有了電荷。原子速度不大的情況下，run一步原子不會跑多遠，這一步有沒有電荷所導致的影響也應該不大。 region depositregionMg block 0.00 25.2672 0.00 25.2672 25.20 25.40 group depositMgatoms region depositregionMg fix 4 depositMgatoms nve fix 5 depositMgatoms deposit 80 1 10000 9375 region depositregionMg vz -3.0 -3.0 dump 1 depositMgatoms custom 100 MgOCharge.dat id type mass q dump 3 all movie 50 movie.mpg type type & axes yes 0.8 0.02 view 80 -30 dump 2 all xyz 300 MgO.xyz run 1 set type 1 charge 2.0 set type 2 charge -2.0 set type 3 charge 3.0 run 4999 region depositregionO block 0.00 25.2672 0.00 25.2672 25.45 25.60 group depositOatoms region depositregionO fix 6 depositOatoms nve fix 7 depositOatoms deposit 80 2 10000 1837 region depositregionO vz -3.0 -3.0 run 1 set type 1 charge 2.0 set type 2 charge -2.0 set type 3 charge 3.0 run 5000 set type 1 charge 2.0 set type 2 charge -2.0 set type 3 charge 3.0 run 5000 |
3樓2014-07-11 10:25:54
5樓2016-07-06 17:27:31
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 311求調劑一志愿合肥工業(yè)大學 +9 | 秋二十二 2026-03-30 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 070300一志愿211，312分求調劑院校 +8 | 小黃鴨寶 2026-03-30 | 8/400 |
|
|
|
[考研] 南京大學化學調劑 +11 | 景隨風 2026-03-29 | 16/800 |
|
|
|
[考研] 266求調劑 +3 | 哇呼哼呼哼 2026-03-29 | 3/150 |
|
|
|
[考研] 0710生物學求調劑！ +3 | 敘述文 2026-03-31 | 3/150 |
|
|
|
[考研] 一志愿西交288化工專碩求調劑 +4 | 好運好運接接 2026-03-24 | 4/200 |
|
|
|
[考研] 085601一志愿中山大學深圳材料工程330求調劑 +5 | pipiver 2026-03-30 | 5/250 |
|
|
|
[考研] 一志愿鄭大材料工程290求調劑 +12 | Youth_ 2026-03-30 | 12/600 |
|
|
|
[考研] 370求調劑 +3 | 080700調劑 2026-03-30 | 3/150 |
|
|
|
[考研] 286求調劑 +5 | Faune 2026-03-30 | 5/250 |
|
|
|
[考研] 290求調劑 +3 | dfffsar 2026-03-29 | 3/150 |
|
|
|
[考研] 330分求調劑 +5 | qzenlc 2026-03-29 | 5/250 |
|
|
|
[考研] 復試調劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
|
[考研]
 一志愿廈門大學化學學碩307求調劑
+10
|
y7czhao 2026-03-26 | 10/500 |
|
|
|
[考研] 求調劑 +6 | 蘆lty 2026-03-25 | 7/350 |
|
|
|
[考研] 275求調劑 +10 | jjjjjjjjjjl 2026-03-27 | 10/500 |
|
|
|
[考研] 265求調劑 +8 | 小木蟲085600 2026-03-27 | 8/400 |
|
|
|
[考研] 一志愿211院校 344分 東北農業(yè)大學生物學學碩，求調劑 +5 | 丶風雪夜歸人丶 2026-03-26 | 8/400 |
|
|
|
[考研]
 一志愿 西北大學 總分282 英語一62 求調劑
+7
|
18419759900 2026-03-25 | 8/400 |
|
|
|
[考研] 調劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|