 |
 |
 |
| 24小時(shí)熱門版塊排行榜 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 7261 | 回復(fù): 11 | ||
wangsong1016銅蟲 (初入文壇)
|
[求助]
fix deposit沉積原子被反彈
|
|
|
本人目前利用Lammps軟件進(jìn)行薄膜生長(zhǎng)模擬。目前參考文獻(xiàn)J. Phys. D: Appl. Phys. 42 (2009) 065107進(jìn)行重復(fù)性工作。當(dāng)然這個(gè)文獻(xiàn)中采用的是DL_Ploy軟件，我利用Lammps應(yīng)該也能重復(fù)該工作。利用Lammps做薄膜生長(zhǎng)很重要的一個(gè)命令式fix deposit命令，能夠以一定的時(shí)間步長(zhǎng)為周期，向體系中周期性的加入一個(gè)原子。位置，角度，速度可控。計(jì)算采用的勢(shì)能及參數(shù)來源于上面的參考文獻(xiàn)，勢(shì)能采用Buckingham+Coulumb勢(shì)，氧離子間的相互作用采用純庫(kù)倫勢(shì)。該文獻(xiàn)研究的是MgO表面上生長(zhǎng)MgAlO合金，我重復(fù)該工作的時(shí)候出現(xiàn)不合理的情況，請(qǐng)大家?guī)兔⒖家幌隆?br />
1.最基本的方面入手，我模擬MgO基底上生長(zhǎng)MgO薄膜。從文件中讀取MS建模得到的MgO坐標(biāo)。勢(shì)參量如上文獻(xiàn)設(shè)置。共四層原子，在第三層位置利用Langevin熱浴，維持一定溫度，上面兩層原子nve系綜自由弛豫。最底層原子固定不動(dòng)。（物理上可能不盡合理，作為初步測(cè)試應(yīng)該足夠） 2.在大概5Ang之外的位置設(shè)定為沉積區(qū)域，原子Mg和O在該區(qū)域產(chǎn)生并沉積下來。 3.方便起見，我目前都只在windows下利用單核進(jìn)行測(cè)試，模型初步建立起來之后再轉(zhuǎn)移到Linux系統(tǒng)下并行計(jì)算。 模擬結(jié)果發(fā)現(xiàn)利用deposit產(chǎn)生的這個(gè)原子行為很古怪，原子在向表面入射的過程中，會(huì)被表面反彈回來，不能吸附在基底上。換一種沉積產(chǎn)生原子的方式：這個(gè)原子不是利用fix deposit命令產(chǎn)生，而是在讀取MgO原子坐標(biāo)的時(shí)候，在MgO基底之外再多讀取一個(gè)Mg或者O的坐標(biāo)，那么這個(gè)原子就能夠吸附到基底表面上。對(duì)這樣的結(jié)果感到很困惑，希望大家給予幫助。 輸入腳本如下，另外坐標(biāo)數(shù)據(jù)見附件： #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 模擬所需時(shí)間極短，希望大家能夠運(yùn)行測(cè)試，多多幫忙，謝謝！ |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : MgO_structure.in
2014-07-10 09:30:34, 24.91 K
» 猜你喜歡
求化學(xué)調(diào)劑
已經(jīng)有5人回復(fù)
-
085600，材料與化工321分求調(diào)劑
已經(jīng)有10人回復(fù)
-
求調(diào)劑
已經(jīng)有4人回復(fù)
-
0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分
已經(jīng)有5人回復(fù)
-
一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑
已經(jīng)有3人回復(fù)
-
085600，專業(yè)課化工原理，320分求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料與化工272求調(diào)劑
已經(jīng)有16人回復(fù)
-
290求調(diào)劑
已經(jīng)有3人回復(fù)
-
327求調(diào)劑
已經(jīng)有3人回復(fù)
-
085600 286分 材料求調(diào)劑
已經(jīng)有11人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 原子層沉積制備TiO2三維納米結(jié)構(gòu)及其光催化性能研究
已經(jīng)有15人回復(fù)
- 求助有關(guān)原子層沉積用鉑源
已經(jīng)有2人回復(fù)
- 原子層沉積復(fù)合薄膜厚退火的作用？誰知道�。�
已經(jīng)有10人回復(fù)
- 泵油可以反復(fù)過濾后重復(fù)使用嗎，原子層沉積系統(tǒng)的真空抽不了，已排除漏氣可能
已經(jīng)有7人回復(fù)
- MS模擬中固定原子Fix cartesian Position和Fix fractional position 不同點(diǎn)
已經(jīng)有7人回復(fù)
- ACCOUTS OF CHEMICAL RESEARCH:原子層沉積在制備和提高負(fù)載型金屬催化劑方面的應(yīng)用
已經(jīng)有63人回復(fù)
- 原子總數(shù)怎么看
已經(jīng)有3人回復(fù)
- 關(guān)于Atomic Layer Deposition（原子層沉積技術(shù)）的一篇綜述
已經(jīng)有322人回復(fù)
- 關(guān)于原子吸收分光光度法檢測(cè)沉積物和海水的檢出限問題
已經(jīng)有7人回復(fù)
- 原子層沉積技術(shù)-ALD
已經(jīng)有7人回復(fù)
wangsong1016
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 169
- 散金: 4
- 帖子: 45
- 在線: 26.3小時(shí)
- 蟲號(hào): 408613
- 注冊(cè): 2007-06-21
- 性別: GG
- 專業(yè): 物理
|
無論利用set type 1 charge 1這種命令來設(shè)置原子電荷，還是利用read_data命令從文件中讀取電荷，都只能對(duì)體系中當(dāng)前存在的原子進(jìn)行電荷賦值，而在這兩個(gè)命令之后，再利用fix deposit命令產(chǎn)生的原子還是不具有電荷。因此每一次fix deposit命令之后，run一次，產(chǎn)生了一個(gè)原子之后立刻用set命令進(jìn)行電荷賦值，這樣原子就具有了電荷。原子速度不大的情況下，run一步原子不會(huì)跑多遠(yuǎn)，這一步有沒有電荷所導(dǎo)致的影響也應(yīng)該不大。 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
zhangjiekm
銀蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 1446.9
- 紅花: 2
- 帖子: 98
- 在線: 251小時(shí)
- 蟲號(hào): 497558
- 注冊(cè): 2008-02-02
- 專業(yè): 無機(jī)材料化學(xué)
7樓2017-06-17 10:32:48
wangsong1016
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 169
- 散金: 4
- 帖子: 45
- 在線: 26.3小時(shí)
- 蟲號(hào): 408613
- 注冊(cè): 2007-06-21
- 性別: GG
- 專業(yè): 物理
2樓2014-07-10 10:56:19
wangsong1016
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 169
- 散金: 4
- 帖子: 45
- 在線: 26.3小時(shí)
- 蟲號(hào): 408613
- 注冊(cè): 2007-06-21
- 性別: GG
- 專業(yè): 物理
4樓2014-07-11 10:26:05
huoxingpopo
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1928
- 帖子: 18
- 在線: 13.4小時(shí)
- 蟲號(hào): 4039980
- 注冊(cè): 2015-08-27
- 專業(yè): 光學(xué)
5樓2016-07-06 17:27:31
feng0603
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2.9
- 散金: 12
- 帖子: 108
- 在線: 61.4小時(shí)
- 蟲號(hào): 2973885
- 注冊(cè): 2014-02-17
- 專業(yè): 微/納機(jī)械系統(tǒng)
6樓2016-10-22 13:45:41
shawn_zhang
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 3.5
- 帖子: 3
- 在線: 57分鐘
- 蟲號(hào): 7389167
- 注冊(cè): 2017-10-29
- 專業(yè): 微/納機(jī)械系統(tǒng)
8樓2018-05-22 10:30:07
19861226ya
木蟲 (著名寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2019
- 紅花: 6
- 帖子: 1351
- 在線: 318小時(shí)
- 蟲號(hào): 842817
- 注冊(cè): 2009-09-07
- 性別: GG
- 專業(yè): 動(dòng)力學(xué)與控制
9樓2018-06-05 17:46:16
10樓2018-10-18 17:35:31
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 wangsong1016 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600，材料與化工321分求調(diào)劑 +10 | 大饞小子 2026-03-28 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 296求調(diào)劑 +6 | 彼岸t 2026-03-29 | 6/300 |
|
|
|
[考研] 化學(xué)工程085602 305分求調(diào)劑 +22 | RichLi_ 2026-03-25 | 22/1100 |
|
|
|
[考研] 一志愿211，335分，0856，求調(diào)劑院校和導(dǎo)師 +5 | 傾____蕭 2026-03-27 | 6/300 |
|
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +7 | 李李不服輸 2026-03-26 | 7/350 |
|
|
|
[考研] 298求調(diào)劑 +3 | 種圣賜 2026-03-29 | 3/150 |
|
|
|
[考研] 求調(diào)劑，一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +7 | @taotao 2026-03-29 | 7/350 |
|
|
|
[考研] 0703化學(xué) +11 | 妮妮ninicgb 2026-03-27 | 11/550 |
|
|
|
[考研] 329求調(diào)劑 +7 | 星野? 2026-03-26 | 7/350 |
|
|
|
[考研] 322求調(diào)劑 +7 | 宋明欣 2026-03-27 | 7/350 |
|
|
|
[考研] 347求調(diào)劑 +3 | 山頂見α 2026-03-25 | 3/150 |
|
|
|
[考研] 085701環(huán)境工程，267求調(diào)劑 +16 | minht 2026-03-26 | 16/800 |
|
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫(kù) 2026-03-26 | 5/250 |
|
|
|
[考研] 材料求調(diào)劑 +5 | .m.. 2026-03-25 | 5/250 |
|
|
|
[考研] 0703化學(xué)338求調(diào)劑！ +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
|
[考研] 327求調(diào)劑 +7 | prayer13 2026-03-23 | 7/350 |
|
|
|
[考研] 318求調(diào)劑 +3 | plum李子 2026-03-23 | 3/150 |
|
|
|
[考研] 一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|