 |
 |
 |
| 24小時(shí)熱門(mén)版塊排行榜 |
- >論壇更新日志 (3376)
- >考研 (709)
- >導(dǎo)師招生 (572)
- >蟲(chóng)友互識(shí) (125)
- >文獻(xiàn)求助 (95)
- >碩博家園 (58)
- >考博 (50)
- >博后之家 (48)
- >招聘信息布告欄 (42)
- >論文投稿 (41)
- >休閑灌水 (31)
- >基金申請(qǐng) (24)
- >找工作 (18)
- >教師之家 (15)
- >文獻(xiàn)求助完結(jié)子版 (14)
- >論文道賀祈福 (13)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 2599 | 回復(fù): 11 | |||
[求助]
困擾很久的難題-結(jié)構(gòu)變形嚴(yán)重-SiO2 已有4人參與
|
|||
|
大家好，這個(gè)問(wèn)題困擾了我很長(zhǎng)時(shí)間，請(qǐng)大家?guī)臀铱纯词悄睦锏膯?wèn)題吧，多謝。 問(wèn)題是這樣的，我想模擬水溶液環(huán)境下納米物質(zhì)在SiO2表面的吸附。但是程序運(yùn)行的時(shí)候總是丟失原子，而且我發(fā)現(xiàn)是SiO2變形較大，因而我首先判斷是SiO2的問(wèn)題。所以我將SiO2單獨(dú)拿出來(lái)模擬。首先懷疑是兩方面原因，一個(gè)是初始模型不合理，因而我試用了MS中的alpha-quartz模型及VMD中的SiO2模型，同時(shí)也將模型的大小做了一些調(diào)整，例如試用了15*15*4,14*14*4及10*10*4大小的模型。另一個(gè)是勢(shì)函數(shù)不合理，因而我試用了BKS，tersoff及comb勢(shì)，其中BKS也采用了table。但是結(jié)果仍然不合理。我將其中一組模擬結(jié)果黏貼在附件里，有兩張圖片，第一張是初始構(gòu)型，第二張是運(yùn)行50000步之后的結(jié)構(gòu)。請(qǐng)大家?guī)兔�，謝謝！圖片直接上傳不上去，給大家造成的麻煩請(qǐng)諒解 |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲(chóng)僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶(hù)自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶(hù)本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : 1.jpg
- 附件 2 : 2.jpg
2015-07-17 16:13:45, 688.83 K
2015-07-17 16:13:53, 688.83 K
» 猜你喜歡
262求調(diào)劑
已經(jīng)有5人回復(fù)
-
0856調(diào)劑
已經(jīng)有6人回復(fù)
-
張芳銘-中國(guó)農(nóng)業(yè)大學(xué)-環(huán)境工程專(zhuān)碩-298
已經(jīng)有6人回復(fù)
-
考研生物與醫(yī)藥調(diào)劑
已經(jīng)有3人回復(fù)
-
274求調(diào)劑
已經(jīng)有5人回復(fù)
-
化學(xué)工程085602 305分求調(diào)劑
已經(jīng)有26人回復(fù)
-
334求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩
已經(jīng)有5人回復(fù)
-
看看我！301求調(diào)劑 081700
已經(jīng)有3人回復(fù)
-
求調(diào)劑 生物學(xué) 377分
已經(jīng)有3人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
oxox6085
專(zhuān)家顧問(wèn) (正式寫(xiě)手)
-
專(zhuān)家經(jīng)驗(yàn): +129 - 模擬EPI: 1
- 應(yīng)助: 69 (初中生)
- 金幣: 8061
- 散金: 222
- 紅花: 39
- 帖子: 473
- 在線: 636.8小時(shí)
- 蟲(chóng)號(hào): 612537
- 注冊(cè): 2008-09-25
- 性別: GG
- 專(zhuān)業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
- 管轄: 分子模擬
6樓2015-07-20 19:02:39
chuntao118
新蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 93 (初中生)
- 金幣: 2590.9
- 散金: 1402
- 紅花: 25
- 帖子: 1611
- 在線: 241.2小時(shí)
- 蟲(chóng)號(hào): 1791567
- 注冊(cè): 2012-05-02
- 專(zhuān)業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
jane8556: 金幣+20, ★★★很有幫助 2015-09-17 09:46:50
感謝參與，應(yīng)助指數(shù) +1
jane8556: 金幣+20, ★★★很有幫助 2015-09-17 09:46:50
| 模擬原子丟失應(yīng)該是邊界條件處理的問(wèn)題吧！個(gè)人認(rèn)為和勢(shì)函數(shù)以及模擬區(qū)域沒(méi)有關(guān)系。因?yàn)樵趧?shì)函數(shù)的選取可以說(shuō)別人給定的，模擬區(qū)域的大小只影響計(jì)算的模擬時(shí)間，怎么都不會(huì)影響到模擬體系的原子數(shù)目的變化的。道理上就講不通，建議：不要意味的用軟件，好好的把道理想明白。目的搞清楚。軟件中的細(xì)節(jié)弄明白。一知半解自然錯(cuò)誤多多。 |
» 本帖已獲得的紅花（最新10朵）
jane85562樓2015-07-18 08:39:15
3樓2015-07-18 09:46:36
聲夢(mèng)奇緣001
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 56 (初中生)
- 金幣: 5597.4
- 紅花: 5
- 帖子: 614
- 在線: 218.7小時(shí)
- 蟲(chóng)號(hào): 2406187
- 注冊(cè): 2013-04-07
- 專(zhuān)業(yè): 化學(xué)工程及工業(yè)化學(xué)
4樓2015-07-18 12:57:12
5樓2015-07-19 20:46:30
oooofish
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 帖子: 1
- 在線: 28分鐘
- 蟲(chóng)號(hào): 3981537
- 注冊(cè): 2015-07-21
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
jane8556: 金幣+10, ★★★很有幫助 2015-09-17 09:59:04
jane8556: 金幣+10, ★★★很有幫助 2015-09-17 09:59:04
 Can't type Chinese in this computer, not on purpose . I also think something must be wrong in the input file. I simulated SiO2 by BKS and Charmm force field, and the structure was ok. And also, can you visualize the structure by using VDW? It's hard to see what happens from the snapshot you pasted. |
7樓2015-07-21 02:50:21
oxox6085
專(zhuān)家顧問(wèn) (正式寫(xiě)手)
-
專(zhuān)家經(jīng)驗(yàn): +129 - 模擬EPI: 1
- 應(yīng)助: 69 (初中生)
- 金幣: 8061
- 散金: 222
- 紅花: 39
- 帖子: 473
- 在線: 636.8小時(shí)
- 蟲(chóng)號(hào): 612537
- 注冊(cè): 2008-09-25
- 性別: GG
- 專(zhuān)業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
- 管轄: 分子模擬
8樓2015-07-21 04:05:15
|
您好，麻煩再幫我看看行嗎？ 我把input和data文件貼在下面： dimension 3 boundary p p p units metal atom_style full read_data data.deposition # Buck potentials pair_style hybrid buck/coul/long 11.0 lj/cut 12.0 kspace_style pppm 1.0e-4 kspace_modify mesh 24 24 30 order 6 pair_coeff 1 1 buck/coul/long 1388.7730 0.3623 175.000 pair_coeff 1 2 buck/coul/long 18003.7572 0.2052 133.5381 pair_coeff 2 2 lj/cut 0.0 0.0 neighbor 2.0 bin neigh_modify every 1 delay 0 check yes velocity all create 100.0 9999999 rot yes dist gaussian fix 1 all nvt temp 100.0 100.0 0.01 drag 2.0 restart 50000 restart.*.zuizhong timestep 0.00001 thermo_style custom step temp press etotal ke pe thermo 100 dump 1 all custom 100 dump.deposition1 type x y z run 100000 ####################################### data file: LAMMPS Atom File(deposition) 6075 atoms 0 bonds 0 angles 0 dihedrals 0 impropers 2 atom types 0 bond types 0 angle types 0 dihedral types 5.108 74.253 xlo xhi 2.943 116.85 ylo yhi 3.897 25.474 zlo zhi -36.825 0.0 0.0 xy xz yz Masses 1 15.9994 2 28.0900 Atoms 1 1 2 2.4 10.153 42.767 6.897 2 1 1 -1.2 8.747 43.254 7.526 3 1 2 2.4 10.316 40.312 10.498 4 1 1 -1.2 10.597 43.762 11.128 5 1 2 2.4 8.108 44.136 8.698 6 1 1 -1.2 9.233 45.11 9.327 7 1 1 -1.2 9.872 41.307 11.67 8 1 1 -1.2 11.236 42.656 8.068 9 1 1 -1.2 11.722 40.799 9.869 10 1 2 2.4 10.153 42.767 12.299 11 1 1 -1.2 8.747 43.254 12.928 12 1 2 2.4 10.316 40.312 15.9 13 1 1 -1.2 10.597 43.762 16.53 14 1 2 2.4 8.108 44.136 14.1 15 1 1 -1.2 9.233 45.11 14.729 16 1 1 -1.2 9.872 41.307 17.072 17 1 1 -1.2 11.236 42.656 13.47 18 1 1 -1.2 11.722 40.799 15.271 19 1 2 2.4 10.153 42.767 17.701 20 1 1 -1.2 8.747 43.254 18.33 21 1 2 2.4 10.316 40.312 21.302 22 1 1 -1.2 10.597 43.762 21.932 23 1 2 2.4 8.108 44.136 19.502 24 1 1 -1.2 9.233 45.11 20.131 25 1 1 -1.2 9.872 41.307 22.474 26 1 1 -1.2 11.236 42.656 18.872 27 1 1 -1.2 11.722 40.799 20.673 28 1 2 2.4 10.153 47.678 6.897 29 1 1 -1.2 8.747 48.165 7.526 30 1 2 2.4 10.316 45.223 10.498 31 1 1 -1.2 10.597 48.672 11.128 32 1 2 2.4 8.108 49.047 8.698 33 1 1 -1.2 9.233 50.021 9.327 34 1 1 -1.2 9.872 46.217 11.67 35 1 1 -1.2 11.236 47.566 8.068 36 1 1 -1.2 11.722 45.71 9.869 37 1 2 2.4 10.153 47.678 12.299 38 1 1 -1.2 8.747 48.165 12.928 39 1 2 2.4 10.316 45.223 15.9 40 1 1 -1.2 10.597 48.672 16.53 41 1 2 2.4 8.108 49.047 14.1 42 1 1 -1.2 9.233 50.021 14.729 43 1 1 -1.2 9.872 46.217 17.072 44 1 1 -1.2 11.236 47.566 13.47 45 1 1 -1.2 11.722 45.71 15.271 46 1 2 2.4 10.153 47.678 17.701 47 1 1 -1.2 8.747 48.165 18.33 48 1 2 2.4 10.316 45.223 21.302 49 1 1 -1.2 10.597 48.672 21.932 50 1 2 2.4 8.108 49.047 19.502 51 1 1 -1.2 9.233 50.021 20.131 52 1 1 -1.2 9.872 46.217 22.474 53 1 1 -1.2 11.236 47.566 18.872 54 1 1 -1.2 11.722 45.71 20.673 55 1 2 2.4 10.153 52.588 6.897 56 1 1 -1.2 8.747 53.075 7.526 57 1 2 2.4 10.316 50.133 10.498 58 1 1 -1.2 10.597 53.582 11.128 59 1 2 2.4 8.108 53.957 8.698 60 1 1 -1.2 9.233 54.931 9.327 61 1 1 -1.2 9.872 51.127 11.67 62 1 1 -1.2 11.236 52.476 8.068 63 1 1 -1.2 11.722 50.62 9.869 64 1 2 2.4 10.153 52.588 12.299 65 1 1 -1.2 8.747 53.075 12.928 66 1 2 2.4 10.316 50.133 15.9 67 1 1 -1.2 10.597 53.582 16.53 68 1 2 2.4 8.108 53.957 14.1 69 1 1 -1.2 9.233 54.931 14.729 70 1 1 -1.2 9.872 51.127 17.072 71 1 1 -1.2 11.236 52.476 13.47 ............................................................................. SiO2的勢(shì)函數(shù)為tersoff時(shí)，將勢(shì)函數(shù)類(lèi)型及參數(shù)變動(dòng)如下： pair_style tersoff pair_coeff * * SiO.tersoff Si O 麻煩幫我看看是否正確，多謝！ |
9樓2015-09-17 09:58:33
snowislands
新蟲(chóng) (初入文壇)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 354
- 帖子: 29
- 在線: 39.7小時(shí)
- 蟲(chóng)號(hào): 3287010
- 注冊(cè): 2014-06-22
- 性別: MM
- 專(zhuān)業(yè): 能源化工
10樓2015-12-28 20:24:07
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 jane8556 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0856調(diào)劑 +6 | 曲聽(tīng)筠 2026-03-30 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-30 | 5/250 |
|
|
[考研]
 物理學(xué)調(diào)劑
+3
|
小羊36 2026-03-30 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | 研s. 2026-03-25 | 3/150 |
|
|
|
[考研] 08工科求調(diào)劑286 +5 | tgs_001 2026-03-28 | 5/250 |
|
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +11 | 多久上課 2026-03-27 | 12/600 |
|
|
|
[考研] 277跪求調(diào)劑 +8 | 1915668 2026-03-27 | 12/600 |
|
|
|
[考研] 一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑 +25 | 稻妻小編 2026-03-26 | 25/1250 |
|
|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +18 | Mr. Z 2026-03-25 | 19/950 |
|
|
|
[考研] 311求調(diào)劑 +10 | lin0039 2026-03-26 | 10/500 |
|
|
|
[考研] 本科新能源科學(xué)與工程，一志愿華理能動(dòng)285求調(diào)劑 +7 | AZMK 2026-03-28 | 11/550 |
|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +6 | 蘆lty 2026-03-25 | 7/350 |
|
|
|
[考研] 283求調(diào)劑 +7 | A child 2026-03-28 | 7/350 |
|
|
|
[考研] 0856，材料與化工321分求調(diào)劑 +12 | 大饞小子 2026-03-27 | 13/650 |
|
|
|
[考研] 266分求材料化工冶金礦業(yè)等專(zhuān)業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
|
[考研] 08開(kāi)頭275求調(diào)劑 +4 | 拉誰(shuí)不重要 2026-03-26 | 4/200 |
|
|
|
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
|
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4 | 小魚(yú)愛(ài)有機(jī) 2026-03-25 | 4/200 |
|