 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1540 | 回復(fù): 3 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
cgzhang_gg銅蟲 (著名寫手)
|
[求助]
pair_style coul/cut 的測(cè)試
|
|
|
有誰測(cè)試過pair_style coul/cut 的對(duì)勢(shì)， 我通過pair_writes 導(dǎo)出的勢(shì)，與我fortran按照庫倫作用表達(dá)式算出的勢(shì)相差很大 ，定性上都是錯(cuò)的。也就是說，通過pair_writes 導(dǎo)出的 TI_Ti ，Ti_O, O_O 三個(gè)對(duì)勢(shì)都是正的，而且與電荷量沒有關(guān)系。我懷疑是我in.文件里那個(gè)地方設(shè)置的不對(duì)。下面是我的in.文件的內(nèi)容： 大家?guī)臀铱纯茨睦镉胁粚?duì)的地方，或者咱們討論下，好人有好報(bào)！ 注：我測(cè)試的體系TiO2金紅石結(jié)構(gòu)：一個(gè)原胞里面有兩個(gè)Ti和四個(gè)O.,晶格常數(shù)是4.493A。 設(shè)置第一類原子Ti的電荷為2.192e, e為元電荷，一個(gè)質(zhì)子所具有的電量。相應(yīng)的O的電荷設(shè)置為：-1.098e. 單位是按照metal 來確定的。 # TiO2 coul/cut detect right or wrong units metal atom_style charge boundary p p p lattice custom 4.493 a1 1.0 0.0 0.0 a2 0.0 1.0 0.0 a3 0.0 0.0 0.669& basis 0.0 0.0 0.0 basis 0.5 0.5 0.5 & basis 0.303 0.303 0.0 basis 0.697 0.697 0.0 & basis 0.803 0.197 0.5 basis 0.197 0.803 0.5 region mybox block 0 1 0 1 0 1 create_box 2 mybox create_atoms 2 box basis 1 1 basis 2 1 & basis 3 2 basis 4 2 basis 5 2 basis 6 2 mass 1 47.867 # 設(shè)置第一類原子Ti的質(zhì)量：47.867 mass 2 15.9994 # 設(shè)置第二類原子O的質(zhì)量：15.9994 set type 1 charge 2.196 # 設(shè)置第一類原子Ti 的電荷為2.196 set type 2 charge -1.098 # 設(shè)置第二類原子O的電荷為-1.098 # 1-coul/cut 8.5 pair_style coul/cut 8.5 pair_coeff 1 1 pair_coeff 1 2 pair_coeff 2 2 pair_write 1 1 50000 r 0.4 8.0 table.T_T.c Ti_Ti pair_write 1 2 50000 r 0.4 8.0 table.T_O.c Ti_O pair_write 2 2 50000 r 0.4 8.0 table.O_O.c O_O quit |
回復(fù)此樓» 猜你喜歡
08工科275分求調(diào)劑
已經(jīng)有3人回復(fù)
-
288求調(diào)劑 一志愿哈工大 材料與化工
已經(jīng)有14人回復(fù)
-
安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分，英語已過四六級(jí)，六級(jí)521，求調(diào)劑
已經(jīng)有11人回復(fù)
-
070300化學(xué)求調(diào)劑
已經(jīng)有12人回復(fù)
-
化學(xué)308分調(diào)劑
已經(jīng)有12人回復(fù)
-
359求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿武理材料工程302調(diào)劑環(huán)化或化工
已經(jīng)有8人回復(fù)
-
材料工程085601數(shù)二英一335求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿北京科技大學(xué)材料學(xué)碩328分求調(diào)劑
已經(jīng)有3人回復(fù)
-
0856求調(diào)劑
已經(jīng)有9人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- Lammps fixnve 為什么溫度也一直升高��？
已經(jīng)有7人回復(fù)
- LAMMPS中描述C-Si-O勢(shì)能
已經(jīng)有5人回復(fù)
- 運(yùn)行之后水分子怎么亂了,氫氧鍵跑老遠(yuǎn),上圖,上in data文件
已經(jīng)有12人回復(fù)
- 有關(guān)create_atoms相關(guān)命令
已經(jīng)有10人回復(fù)
- 有沒有用 lj勢(shì)函數(shù) 模擬lammps的朋友？
已經(jīng)有6人回復(fù)
- 請(qǐng)幫忙看一下lammps中是否有這個(gè)勢(shì)函數(shù)？
已經(jīng)有3人回復(fù)
- pair_style lj/cut/coul/cut command的使用
已經(jīng)有5人回復(fù)
- 在ubuntu下編譯了lammps,運(yùn)行例子顯示ERROR: Invalid pair style (../force.cpp:175)
已經(jīng)有9人回復(fù)
- pair style hybrid 的用法，請(qǐng)高手幫忙看下對(duì)不對(duì)？
已經(jīng)有6人回復(fù)
- 求助，lammps對(duì)勢(shì)的應(yīng)用（短程為已知函數(shù)，并且?guī)靷愰L(zhǎng)程作用，ewald求和）
已經(jīng)有6人回復(fù)
- lammps 離子液體，離子中相鄰隔的挺遠(yuǎn)的原子成鍵了，不知什么原因
已經(jīng)有8人回復(fù)
- 關(guān)于lammps中的in文件運(yùn)行時(shí)的錯(cuò)誤問題查找
已經(jīng)有10人回復(fù)
- LAMMPS運(yùn)行結(jié)果出現(xiàn)nan
已經(jīng)有20人回復(fù)
- 【討論】lammps運(yùn)行eam/fs出錯(cuò)
已經(jīng)有38人回復(fù)
- 【求助】LAMMPS中非周期性邊界條件和Born-Mayer-Huggins勢(shì)函數(shù)不會(huì)設(shè)置
已經(jīng)有4人回復(fù)
cgzhang_gg
銅蟲 (著名寫手)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 3.7
- 散金: 790
- 紅花: 5
- 帖子: 1096
- 在線: 479.2小時(shí)
- 蟲號(hào): 1521579
- 注冊(cè): 2011-12-03
- 專業(yè): 極端條件下使用的金屬材料
送紅花一朵 |
經(jīng)過測(cè)試，可以很明確的說，用pair_write 導(dǎo)出庫倫對(duì)勢(shì)的時(shí)候沒有考慮電荷項(xiàng)，只是導(dǎo)出了：C/r 這一項(xiàng)， r 按照所用單位值確定，C是 energy-conversion constant 。如units metal r 的單位為 A ，C的值為14.399645 eV.A, （C/r）的單位為eV 。 但是在程序內(nèi)部計(jì)算時(shí)，是考慮電荷的。 atom_style charge 。。。。。。。。。。。。。 set type 1 charge 2.196 # 設(shè)置第一類原子Ti 的電荷為2.196 set type 2 charge -1.098 # 設(shè)置第二類原子O的電荷為-1.098 。。。。。。。。。。。。。 這樣設(shè)置來模擬coul作用是完全正確的。 |
3樓2013-06-18 09:55:54
cgzhang_gg
銅蟲 (著名寫手)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 3.7
- 散金: 790
- 紅花: 5
- 帖子: 1096
- 在線: 479.2小時(shí)
- 蟲號(hào): 1521579
- 注冊(cè): 2011-12-03
- 專業(yè): 極端條件下使用的金屬材料
|
我再補(bǔ)充下，我把用pair_write 導(dǎo)出來的對(duì)勢(shì)，T_T，T_O,O_O的三個(gè)對(duì)勢(shì)與 1/r 相比較，結(jié)果發(fā)現(xiàn)比值恰為：14.399645，這正好是energy-conversion constant C 見lammps手冊(cè) pair_style coul/cut command 大約872頁， 然后我又測(cè)試了通過read_data 命令來讀取構(gòu)型以及每個(gè)原子的帶電量，然后又通過pair_write 導(dǎo)出對(duì)勢(shì)，結(jié)果還是不對(duì)？我現(xiàn)在懷疑lammps 里面有一個(gè)bug, 就是通過pair_write 命令導(dǎo)出對(duì)勢(shì)的時(shí)候根本沒考慮qi 和qj，至于程序內(nèi)部計(jì)算力，以便更新位置和速度，應(yīng)該考慮了，要不然會(huì)出現(xiàn)很多離奇的錯(cuò)誤。我測(cè)試的in.文件如下： # TiO2 TiO2.pair.table detect right or wrong units metal atom_style charge #atom_style atomic boundary p p p read_data data.TiO2 mass 1 47.867 # Ti mass 2 15.9994 # O # 1-coul/cut 8.5 pair_style coul/cut 8.5 #pair_coeff * * pair_coeff 1 1 pair_coeff 1 2 pair_coeff 2 2 pair_write 1 1 50000 r 0.4 9.0 table.T_T.c Ti_Ti pair_write 1 2 50000 r 0.4 9.0 table.T_O.c Ti_O pair_write 2 2 50000 r 0.4 9.0 table.O_O.c O_O fix 1 all nve run 10 quit data.TiO2 的內(nèi)容如下： #TiO2 6 atoms 0 bonds 0 angles 0 dihedrals 0 impropers 2 atom types 0 bond types 0 angle types 0 dihedral types 0 improper types 0 4.493 xlo xhi 0 4.493 ylo yhi 0 3.00582 zlo zhi Atoms 1 1 2.196 0 0 0 2 1 2.196 2.2465 2.2465 1.50291 3 2 -1.098 1.36138 1.36138 0 4 2 -1.098 3.13162 3.13162 0 5 2 -1.098 3.60788 0.885121 1.50291 6 2 -1.098 0.885121 3.60788 1.50291 在這里為了簡(jiǎn)單，我只取了一個(gè)原胞的原子。即，兩個(gè)Ti，一個(gè)O. Atoms 每一列的意義如下： id type q x y z |
» 本帖已獲得的紅花（最新10朵）
2樓2013-06-17 22:49:16
cgzhang_gg
銅蟲 (著名寫手)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 3.7
- 散金: 790
- 紅花: 5
- 帖子: 1096
- 在線: 479.2小時(shí)
- 蟲號(hào): 1521579
- 注冊(cè): 2011-12-03
- 專業(yè): 極端條件下使用的金屬材料
4樓2013-06-18 09:57:10
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿華東理工大學(xué)，080500學(xué)碩，317分，求調(diào)劑 +9 | s1145 2026-03-31 | 9/450 |
|
|---|---|---|---|---|
|
|
[考研] 材料調(diào)劑 +9 | Eujd1 2026-03-31 | 10/500 |
|
|
|
[考研] 考研調(diào)劑求助 +7 | 13287130938 2026-03-31 | 7/350 |
|
|
[考研] 311求調(diào)劑 +9 | 藍(lán)月亮亮 2026-03-30 | 9/450 |
|
|
|
[考研] 343求調(diào)劑 +8 | 愛羈絆 2026-03-28 | 8/400 |
|
|
|
[考研] 求調(diào)劑 +8 | 11ggg 2026-03-30 | 8/400 |
|
|
|
[考研] 085600 295分求調(diào)劑 +13 | W55j 2026-03-30 | 15/750 |
|
|
|
[考研] 070300化學(xué)279求調(diào)劑 +7 | 哈哈哈^_^ 2026-03-31 | 9/450 |
|
|
|
[有機(jī)交流]
 甲基亞磺磺酸鈉和甲基磺酸酯反應(yīng)機(jī)理
 10+3
|
kaobao456 2026-03-29 | 4/200 |
|
|
|
[考研] 11408總分309，一志愿東南大學(xué)求調(diào)劑，不挑專業(yè) +5 | 天賦帶到THU 2026-03-29 | 6/300 |
|
|
|
[考研] 一志愿華中師范化學(xué)332分求調(diào)劑 +3 | Lyy930824@ 2026-03-29 | 3/150 |
|
|
|
[考研] 332求調(diào)劑 +6 | @MZB382400 2026-03-28 | 6/300 |
|
|
|
[考研] 303求調(diào)劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
|
[考研] 085602 化學(xué)工程專碩 340分求調(diào)劑 +4 | qianbai11 2026-03-29 | 4/200 |
|
|
|
[考研] 0856，材料與化工321分求調(diào)劑 +12 | 大饞小子 2026-03-27 | 13/650 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑，各位老師看看我�。�！ +5 | 祁祺祺 2026-03-25 | 5/250 |
|
|
|
[考研] 281求調(diào)劑 +6 | Koxui 2026-03-24 | 7/350 |
|
|
|
[考研] 303求調(diào)劑 +6 | 藍(lán)山月 2026-03-25 | 6/300 |
|
|
|
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
|
[考研] 296求調(diào)劑 +4 | 汪��？！ 2026-03-25 | 7/350 |
|