 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 4238 | 回復(fù): 14 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
[交流]
【求助】lammps中多元化合物的晶格創(chuàng)建問題 已有9人參與
|
|||
|
lattice custom 3.104 a1 0.8660254037844386 -0.5 0 a2 0 1.0 0 a3 0 0 1.602448454 & basis 0.3333333333 0.6666666667 0.3826000000 & basis 0.6666666667 0.3333333333 0.8826000000 & basis 0.3333333333 0.6666666667 0.0000000000 & basis 0.6666666667 0.3333333333 0.5000000000 region box block 0 3 0 3 0 3 create_box 2 box create_atoms 1 box basis 1 1 basis 2 1 create_atoms 2 box basis 3 2 basis 4 2 我想創(chuàng)建一個(gè)二元化合物，建立了一個(gè)晶格，晶格中有四個(gè)basis，我想在第1個(gè)和第2個(gè)的basis上創(chuàng)建A原子，在3、4上創(chuàng)建B原子，我上邊的代碼卻始終不成功，The Fuc*ing Manual 上寫的create_atoms 中basis的參數(shù)是： basis values = M itype M = which basis atom itype = atom type (1-N) to assign to this basis atom 我不知道我是哪里錯(cuò)了，新人，望大家指教，謝謝! |
回復(fù)此樓» 猜你喜歡
309求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷
已經(jīng)有8人回復(fù)
-
求調(diào)劑
已經(jīng)有5人回復(fù)
-
求化學(xué)調(diào)劑
已經(jīng)有5人回復(fù)
-
085600，材料與化工321分求調(diào)劑
已經(jīng)有10人回復(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ù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 在Lammps中安裝xmovie工具時(shí)，Makefile文件怎么修改X庫路徑
已經(jīng)有3人回復(fù)
- lammps的時(shí)間步長問題
已經(jīng)有16人回復(fù)
- lammps中固定鍵長
已經(jīng)有4人回復(fù)
- 【求助】lammps安裝過程中，在安裝fftw的過程中遇到問題
已經(jīng)有8人回復(fù)
- 【討論】lammps中利用velocity命令無法控制住速度問題
已經(jīng)有11人回復(fù)
- 【求助】lammps如何測量模擬以后材料的晶格參數(shù)？
已經(jīng)有7人回復(fù)
- 【討論】lammps中關(guān)于fix ave/spatial命令的輸出問題
已經(jīng)有4人回復(fù)
- 【求助】MS建模導(dǎo)入lammps問題求救
已經(jīng)有14人回復(fù)
- 【求助】Lammps 多機(jī)并行的問題
已經(jīng)有11人回復(fù)
- 【求助】LAMMPS中非周期性邊界條件和Born-Mayer-Huggins勢函數(shù)不會(huì)設(shè)置
已經(jīng)有4人回復(fù)
老虎大王
木蟲 (著名寫手)
- 模擬EPI: 17
- 應(yīng)助: 26 (小學(xué)生)
- 貴賓: 0.17
- 金幣: 4774.1
- 散金: 8
- 紅花: 42
- 帖子: 1361
- 在線: 215.2小時(shí)
- 蟲號(hào): 659094
- 注冊(cè): 2008-11-21
- 專業(yè): 金屬結(jié)構(gòu)材料
4樓2010-05-09 21:13:29
老虎大王
木蟲 (著名寫手)
- 模擬EPI: 17
- 應(yīng)助: 26 (小學(xué)生)
- 貴賓: 0.17
- 金幣: 4774.1
- 散金: 8
- 紅花: 42
- 帖子: 1361
- 在線: 215.2小時(shí)
- 蟲號(hào): 659094
- 注冊(cè): 2008-11-21
- 專業(yè): 金屬結(jié)構(gòu)材料
2樓2010-05-08 21:29:48
|
首先謝謝上邊的回復(fù)，我的模擬是要對(duì)晶格的尺寸做變化的，其實(shí)在我正式模擬的in文件中晶格參數(shù)里面是有變量的，所以我的原子坐標(biāo)沒有采用read_data讀入，上邊貼出的代碼只是我做測試用的，經(jīng)過測試后發(fā)現(xiàn)是create_atoms命令使用不正確，下邊是我測試的代碼，運(yùn)行時(shí)沒有問題的，我的box當(dāng)中原則上只有162個(gè)basis，可是由運(yùn)行的結(jié)果來看兩個(gè)create_atoms都創(chuàng)建了162個(gè)atoms，共324個(gè)atoms，第一個(gè)creat_atoms我是想讓他在在第一種和第二種類型的basis上創(chuàng)建Ga原子上，第一種類型和第二種類型的basis加起來才162/2個(gè)，可是他卻創(chuàng)建了162個(gè)，所以我想請(qǐng)問大家create_atoms帶basis參數(shù)的命令應(yīng)該是怎么用，謝謝！ units metal boundary p p p atom_style atomic lattice custom 3.104 a1 0.8660254037844386 -0.5 0 a2 0 1.0 0 a3 0 0 1.602448454 basis 0.3333333333 0.6666666667 0.3826000000 basis 0.66666666 67 0.3333333333 0.8826000000 basis 0.3333333333 0.6666666667 0.0000000000 basis 0.6666666667 0.3333333333 0.5000000000 Lattice spacing in x,y,z = 2.68814 4.656 4.974 region box block 0 3 0 3 0 3 create_box 2 box Created orthogonal box = (0 0 0) to (8.06443 13.968 14.922) 1 by 1 by 1 processor grid create_atoms 1 box basis 1 1 basis 2 1 Created 162 atoms create_atoms 2 box basis 3 2 basis 4 2 Created 162 atoms pair_style tersoff pair_coeff * * GaN.tersoff Ga N mass 1 69.723 mass 2 14.0067 neighbor 1.0 bin neigh_modify every 1 delay 5 check yes dump 1 all cfg 1 a*.cfg id type xs ys zs #dump_modify 1 element Al N run 0 Memory usage per processor = 1.831 Mbytes Step Temp E_pair E_mol TotEng Press 0 0 nan 0 nan nan Loop time of 0 on 1 procs for 0 steps with 324 atoms Pair time (%) = 0 (0) Neigh time (%) = 0 (0) Comm time (%) = 0 (0) Outpt time (%) = 0 (0) Other time (%) = 0 (0) Nlocal: 324 ave 324 max 324 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 1188 ave 1188 max 1188 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 17364 ave 17364 max 17364 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 17364 Ave neighs/atom = 53.5926 Neighbor list builds = 0 Dangerous builds = 0 |
3樓2010-05-09 12:47:27
老虎大王
木蟲 (著名寫手)
- 模擬EPI: 17
- 應(yīng)助: 26 (小學(xué)生)
- 貴賓: 0.17
- 金幣: 4774.1
- 散金: 8
- 紅花: 42
- 帖子: 1361
- 在線: 215.2小時(shí)
- 蟲號(hào): 659094
- 注冊(cè): 2008-11-21
- 專業(yè): 金屬結(jié)構(gòu)材料
★ ★ ★ ★ ★ ★ ★ ★ ★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
ghcacj(金幣+8):謝謝 2010-05-10 08:23:47
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
ghcacj(金幣+8):謝謝 2010-05-10 08:23:47
|
好的。 你搞的是GaN，建立的是一個(gè)六方的晶格，。原子位置的設(shè)置是正確的。晶格尺寸不是實(shí)驗(yàn)值，因?yàn)槟阋?br /> 研究晶格尺寸的變化，是的吧？ 我對(duì)你這個(gè)script有這樣幾個(gè)評(píng)論： 一、 首先 create_atoms 1 box basis 1 1 basis 2 1 這樣的用法是錯(cuò)誤的。 請(qǐng)看Manual中create_atoms的說明： For the box style, the create_atoms command fills the entire simulation box with atoms on the lattice. 我測試了 create_atoms 1 box basis 1 1 basis 2 1 create_atoms 1 box 這兩個(gè)語句，結(jié)果都是一樣的。說明后面的basis好像不起作用。我不知道這是程序的bug還是這個(gè)地方確實(shí)不能用box。 反正，經(jīng)測試，create_atoms 1 box basis 1 1 basis 2 1 這一句中，不能用box，而要用region。 二、你改了上述錯(cuò)誤之后（即把這兩句中的box 改成region box，其中后面的box是region的ID），仍然 會(huì)發(fā)現(xiàn)有問題。因?yàn)樯傻脑觽�(gè)數(shù)不對(duì)。這是為什么呢？ 你看你的script，其中有： region box block 0 3 0 3 0 3 create_box 2 box 這兩句。 系統(tǒng)執(zhí)行了以后，返回的信息是：Created orthogonal box = (0 0 0) to (8.06443 13.968 14.922)， 這說明你生成的是一個(gè)方盒子。實(shí)際上，你要的盒子是個(gè)什么盒子呢？你是六角晶系，當(dāng)然應(yīng)是一個(gè)平行 六面體（你可以照著你自己寫的lattice語句，在紙上畫畫）。 所以，你region box block 0 3 0 3 0 3這一句當(dāng)中的style用bolck也是不對(duì)的，正確的設(shè)置應(yīng)該用 prism,至于如何用，請(qǐng)你讀Manual。 以上意見供參考。若仍有問題，可以繼續(xù)討論。 |
5樓2010-05-09 21:33:59
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑 +8 | mikasa的圍巾 2026-03-28 | 8/400 |
|
|---|---|---|---|---|
|
|
[考研] 375求調(diào)劑 +4 | 雨夏整夜 2026-03-29 | 4/200 |
|
|
|
[考研] 332求調(diào)劑 +8 | 蕉蕉123 2026-03-28 | 8/400 |
|
|
|
[考研] 332求調(diào)劑 +4 | @MZB382400 2026-03-28 | 4/200 |
|
|
|
[考研] 343求調(diào)劑 +5 | 愛羈絆 2026-03-28 | 5/250 |
|
|
|
[考研] 本科新能源科學(xué)與工程，一志愿華理能動(dòng)285求調(diào)劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8 | 曼111 2026-03-24 | 9/450 |
|
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
|
[考研] 315分求調(diào)劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|
|
|
[考研] 322求調(diào)劑 +4 | 我真的很想學(xué)習(xí) 2026-03-23 | 4/200 |
|
|
|
[考研] 調(diào)劑推薦 +5 | 清酒714 2026-03-26 | 6/300 |
|
|
|
[考研]
 304材料求調(diào)劑
+4
|
鐘llll 2026-03-26 | 4/200 |
|
|
|
[考研] 調(diào)劑求收留 +7 | 果然有我 2026-03-26 | 7/350 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|
|
|
[考研] 332求調(diào)劑 +6 | 032500 2026-03-25 | 6/300 |
|
|
|
[考研] 各位老師您好：本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|
|
|
[考研] 277分求調(diào)劑，跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|