版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

anviacc

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 3
帖子: 1
在線: 49分鐘
蟲號: 7415044
注冊: 2017-11-02
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

[求助] 求助用lammps模擬銅納米熔化過程，不同溫度下的rdf不變？已有1人參與

初學(xué)lammps，求銅納米熔化過程不同溫度下的rdf函數(shù)？不知哪里出了錯誤,求大神指點

clear
units metal
dimension 3
boundary p p p
atom_style atomic
atom_modify map array

# ---------- Create Atoms ---------------------

lattice fcc 3.61
region box block 0 200 0 200 0 200 units lattice
create_box 1 box

lattice fcc 3.61
region my1 sphere 100.0 100.0 100.0 2.8
create_atoms 1 region my1

# ---------- Define Interatomic Potential ---------------------

pair_style eam/alloy
pair_coeff * * Cu_mishin1.eam.alloy Cu

# ---------- Run Minimization ---------------------
reset_timestep 0
velocity all create 700.0 4928459
#fix 1 all npt temp 900 900 1 iso 0 0 1000

neighbor 1.0 bin
neigh_modify delay 5 check yes

compute myRDF all rdf 50
fix 1 all ave/time 200 1 200 c_myRDF
file tmp.rdf mode vector

thermo 1000
thermo_style custom step temp pe lx ly lz
dump my all movie 10 my.avi type type
dump my1 all xyz 1000 file.xtc
run 50000

print "Total energy (eV) = ${teng};"
print "Number of atoms = ${natoms};"
print "Lattice constant (Angstoms) = ${length};"
print "Cohesive energy (eV) = ${ecoh};"

print "All done!"

回復(fù)此樓

» 猜你喜歡

081700 調(diào)劑 267分已經(jīng)有8人回復(fù)
請教下大家 2026年國家基金申請是雙盲審嗎？已經(jīng)有5人回復(fù)
276求調(diào)劑。有半年電池和半年高分子實習(xí)經(jīng)歷已經(jīng)有10人回復(fù)
一志愿南航材料專317分求調(diào)劑已經(jīng)有4人回復(fù)
求調(diào)劑已經(jīng)有4人回復(fù)
生物學(xué)學(xué)碩求調(diào)劑已經(jīng)有5人回復(fù)
284求調(diào)劑已經(jīng)有10人回復(fù)
一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑已經(jīng)有4人回復(fù)
07化學(xué)280分求調(diào)劑已經(jīng)有4人回復(fù)
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑已經(jīng)有12人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

1樓 2018-03-24 18:49:43

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

brucefan

專家顧問 (著名寫手)

專家經(jīng)驗: +218
模擬EPI: 4
應(yīng)助: 422 (碩士)
金幣: 2294.8
散金: 1720
紅花: 88
帖子: 1423
在線: 670.2小時
蟲號: 3048965
注冊: 2014-03-14
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
管轄: 分子模擬

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

你為什么把fix npt 注釋掉？沒它怎么積分？

[ 發(fā)自手機版 http://www.gaoyang168.com/3g ]

贊一下

回復(fù)此樓

2樓2018-03-24 20:37:30

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 anviacc 的主題更新