| 1 | 1/1 | 返回列表 |
| 查看: 1031 | 回復(fù): 0 | ||
| 【懸賞金幣】回答本帖問題,作者袁遠(yuǎn)淵將贈送您 10 個金幣 | ||
袁遠(yuǎn)淵新蟲 (初入文壇)
|
[求助]
想向各位大佬請教一下怎么正確使用網(wǎng)站提供的勢函數(shù)文件。
|
|
|
各位大佬,我導(dǎo)師想讓我做Co-Pt沉積的分子模擬,我在interatomic potentials repository上找到了Co-Pt的MEAM勢函數(shù)文件,通過將數(shù)據(jù)拷貝到清空內(nèi)容的原本勢函數(shù)文件里,但是程序一直無法運(yùn)行,每次到pair_coeff就會出現(xiàn)error,顯示如下:incorrect args for pair coefficient。 然后我對查了錯誤手冊,說是初始模型問題描述很短。然后我對in文件做了一些改動,比如將meam/c改為meam,不行。將網(wǎng)站的liburary.meam 中關(guān)于兩個元素的數(shù)據(jù)復(fù)制進(jìn)原本的liburary.meam也不行,將pari style與pair_coeff 移動到分層建模的命令前面,也不行。 模擬卡在這里已經(jīng)很久了,我實(shí)在不知道怎么辦了,希望各位大佬幫幫孩子吧,球球了,研二了我現(xiàn)在只想畢業(yè)。 這是我的in文件 ##Co deposited on Co #mpirun -np 16 ./lmp_g++_mpich<Co.txt #定義數(shù)值 variable lat_a equal 2.51 variable lat_c equal 4.10 variable lat_x equal 20*${lat_a} variable lat_y equal 20*${lat_a} variable lat_z equal 50*${lat_c} variable lat_z1 equal 1*${lat_c} variable lat_z2 equal 2*${lat_c} variable lat_z3 equal 3*${lat_c} variable lat_z4 equal 49*${lat_c} #定義晶格 units meta ldimension 3 atom_style atomic boundary p p f neighbor 0.5 bin neigh_modify every 1 delay 0 check yes #創(chuàng)建格子 region bbox block 0 ${lat_x} 0 ${lat_y} 0 ${lat_z} units box create_box 4 bbox #創(chuàng)建Co基底 lattice hcp ${lat_a} region lbox block 0 ${lat_x} 0 ${lat_y} 0 ${lat_z3} units box create_atoms 1 region lbox #定義基板功能區(qū) region sput block 2 50 2 50 ${lat_z4} ${lat_z} units box #沉積原子產(chǎn)生區(qū) group sput region sput region stable block INF INF INF INF 0.0 ${lat_z1} units box #固定層 group stable region stable region relax block INF INF INF INF ${lat_z1} ${lat_z2} units box #熱控制層 group relax region relax region free block INF INF INF INF ${lat_z2} ${lat_z3} units box #自由層 group free region free #摩爾質(zhì)量設(shè)置 mass 1 58.93 mass 2 58.93 mass 3 195.1 #力場設(shè)置 pair_style meam/c pair_coeff * * library.meam Pt Co PtCo.meam Co Co Pt #計(jì)算原子溫度 compute KE all ke/atom variable KB equal 8.625e-5 variable TEMP atom c_KE/1.5/${KB} #輸出原子信息 dump 1d all custom 1000 dump.diposite id type x y z v_TEMP #屏幕輸出信息設(shè)置 reset_timestep 0 thermo 5000 thermo_style custom step pe ke etotal temp vol press #固定底部原子,防止在沉積過程中移動 fix s1 stable setforce 0.0 0.0 0.0 #能量最小化 min_style cg minimize 1e-8 1e-8 10000 10000 #溫度初始化 velocity other create 298 39849 mom yes rot yes dist gaussian #nvt系綜下進(jìn)行弛豫 fix 1 all nvt temp 298.00 298.0 10.0 run 20000 #沉積前設(shè)定 unfix 1 fix 1 relax nvt temp 298.00 298.00 10.0 fix 2 free nve #沉積原子類型為2,數(shù)量500個 fix 3 all deposit 500 2 100 95485 region sput near 2.51 vz -0.5 -1.5 #為生成的原子添加溫度 fix 4 sput temp/csvr 1653.00 1653.00 100 #刪除飛出box外的原子 thermo_modify lost ignore flush yes run 100000@smutao@oxox6085 |
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 331求調(diào)劑 +3 | zzZ&zZ 2026-03-03 | 3/150 |
|
|---|---|---|---|---|
|
[基金申請]
|
xhuama 2026-03-02 | 11/550 |
|
|
[考博] 26申博 +5 | north, 2026-02-28 | 5/250 |
|
|
[考研] 085701環(huán)境工程 求調(diào)劑 +5 | xiiiia 2026-03-04 | 5/250 |
|
|
[考研] 一志愿西工大材料學(xué),英一數(shù)二總分321分,求調(diào)劑。 +3 | zz05zz 2026-03-02 | 5/250 |
|
|
[考研] 276求調(diào)劑 +8 | 路lyh123 2026-02-28 | 10/500 |
|
|
[基金申請] 沒有青基直接申請面上,感覺自己瘋了 +5 | kevin63t 2026-03-02 | 6/300 |
|
|
[考研] 0857調(diào)劑 +6 | 一ll半 2026-02-28 | 8/400 |
|
|
[考研] 課題組接收材料類調(diào)劑研究生 +6 | gaoxiaoniuma 2026-02-28 | 9/450 |
|
|
[考研] 化工335求調(diào)劑 +5 | 摸摸貓貓頭 2026-03-02 | 5/250 |
|
|
[考研] 282求調(diào)劑 +4 | 2103240126 2026-03-02 | 7/350 |
|
|
[考研] 0856材料調(diào)劑 +5 | 沿岸有貝殼OUC 2026-03-02 | 5/250 |
|
|
[考研] 0856材料求調(diào)劑 +12 | hyf hyf hyf 2026-02-28 | 13/650 |
|
|
[考研] 302材料工程求調(diào)劑 +5 | Doleres 2026-03-01 | 6/300 |
|
|
[考研] 材料085601調(diào)劑 +5 | 多多子. 2026-03-02 | 5/250 |
|
|
[考研] 材料與化工328求調(diào)劑 +3 | 。,。,。,。i 2026-03-02 | 3/150 |
|
|
[考研] 284求調(diào)劑 +10 | 天下熯 2026-02-28 | 11/550 |
|
|
[考研] 279求調(diào)劑 +3 | dua1 2026-03-01 | 4/200 |
|
|
[考研] 299求調(diào)劑 +3 | Y墨明棋妙Y 2026-02-28 | 5/250 |
|
|
[考研] 328求調(diào)劑 +3 | aaadim 2026-03-01 | 5/250 |
|