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

登錄注冊

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

24小時(shí)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

digghost

木蟲 (著名寫手)

應(yīng)助: 343 (大學(xué)生)
金幣: 120.2
散金: 510
紅花: 16
帖子: 1237
在線: 215.5小時(shí)
蟲號: 576827
注冊: 2008-06-21
性別: GG
專業(yè): 金屬材料的力學(xué)行為

[求助] minimize 為什么不能得到其他的晶格結(jié)構(gòu)

下面的in文件本來是fcc結(jié)構(gòu)作優(yōu)化，得到穩(wěn)定值4.05，我將fcc改成了bcc，結(jié)果優(yōu)化后結(jié)構(gòu)還是bcc，這是為什么？
然后我將溫度升高到4000 K使體系完全熔化，在降溫到100 K 結(jié)構(gòu)還是bcc結(jié)構(gòu)的。這是怎么回事？難道優(yōu)化還和設(shè)定的結(jié)構(gòu)有關(guān)？

# Find minimum energy fcc configuration
# Mark Tschopp, 2010

# ---------- Initialize Simulation ---------------------
clear
units metal
dimension 3
boundary p p p
atom_style atomic
atom_modify map array

variable natoms equal "count(all)"
variable teng equal "c_eatoms"
variable length equal "lx"
variable ecoh equal "v_teng/v_natoms"

# ---------- Create Atoms ---------------------
##############lattice bcc 3
region box block 0 1 0 1 0 1 units lattice
create_box 1 box

##########lattice bcc 3 orient x 1 0 0 orient y 0 1 0 orient z 0 0 1
create_atoms 1 box
replicate 1 1 1

# ---------- Define Interatomic Potential ---------------------
pair_style eam/alloy
pair_coeff * * Al99.eam.alloy Al
neighbor 2.0 bin
neigh_modify delay 10 check yes

# ---------- Define Settings ---------------------
compute eng all pe/atom
compute eatoms all reduce sum c_eng

# ---------- Run Minimization ---------------------
reset_timestep 0
fix 1 all box/relax iso 0.0 vmax 0.001

thermo 1000
thermo_style custom step pe lx ly lz press pxx pyy pzz c_eatoms

min_style cg
minimize 1e-25 1e-25 5000 10000

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

print "All done!"

回復(fù)此樓

» 猜你喜歡

本科新能源科學(xué)與工程，一志愿華理能動285求調(diào)劑已經(jīng)有5人回復(fù)
一志愿華北電力大學(xué)能動專碩，293，求調(diào)劑已經(jīng)有5人回復(fù)
調(diào)劑310 已經(jīng)有7人回復(fù)
食品工程專碩一志愿中海洋309求調(diào)劑已經(jīng)有6人回復(fù)
275求調(diào)劑已經(jīng)有14人回復(fù)
283求調(diào)劑已經(jīng)有3人回復(fù)
085602 化工專碩 338分求調(diào)劑已經(jīng)有12人回復(fù)
312，生物學(xué)求調(diào)劑已經(jīng)有3人回復(fù)
化學(xué)0703 調(diào)劑 306分一志愿211 已經(jīng)有4人回復(fù)
311（085601）求調(diào)劑已經(jīng)有3人回復(fù)

呵呵了

1樓 2014-08-27 23:05:28

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

智能機(jī)器人

Robot (super robot)

我們都愛小木蟲

找到一些相關(guān)的精華帖子，希望有用哦~

外行人幫人求碳化氮晶體5種相的晶格結(jié)構(gòu)圖（要求能清楚知道軌道、能帶結(jié)構(gòu)、晶格）已經(jīng)有3人回復(fù)
使用CASTEP進(jìn)行結(jié)構(gòu)優(yōu)化的問題已經(jīng)有25人回復(fù)
這種晶格結(jié)構(gòu)是用什么軟件做出的已經(jīng)有4人回復(fù)
求助超晶格體系的能帶結(jié)構(gòu)計(jì)算已經(jīng)有7人回復(fù)
如何導(dǎo)出高分辨晶格結(jié)構(gòu)圖，已經(jīng)有9人回復(fù)
晶格弛豫與結(jié)構(gòu)弛豫已經(jīng)有18人回復(fù)
關(guān)于鋼中馬氏體的晶格結(jié)構(gòu) 已經(jīng)有5人回復(fù)
求教hcp結(jié)構(gòu)的晶格矢量和k點(diǎn)的設(shè)置已經(jīng)有3人回復(fù)
VASP結(jié)構(gòu)優(yōu)化ISIF=3，為什么晶格常數(shù)不變已經(jīng)有12人回復(fù)
怎樣建立摻雜晶格的結(jié)構(gòu)模型？已經(jīng)有18人回復(fù)
如何在MS里構(gòu)造Co的晶格結(jié)構(gòu)并計(jì)算其體系能量已經(jīng)有5人回復(fù)
一次型NiAl結(jié)構(gòu)優(yōu)化，晶格矢量變化非常大，求助~ 已經(jīng)有5人回復(fù)
求Cu3Si的晶體結(jié)構(gòu)和晶格常數(shù) 已經(jīng)有6人回復(fù)
求查詢金屬化合物的晶格參數(shù),晶體結(jié)構(gòu),空間群,以及原子位置相關(guān)的網(wǎng)站或者書籍已經(jīng)有15人回復(fù)
畫晶格結(jié)構(gòu)小軟件（綠色版）已經(jīng)有169人回復(fù)
各位高人！孿晶有沒有可能是超晶格結(jié)構(gòu)？已經(jīng)有10人回復(fù)
【500金幣本版活動】“經(jīng)驗(yàn)力場分子模擬+第一性原理計(jì)算問題討論互助及資源共享” 已經(jīng)有82人回復(fù)
【討論】想學(xué)lammps是否應(yīng)該先學(xué)習(xí)C++ 已經(jīng)有28人回復(fù)
【求助】75BiFeO3-25BaTiO3 相結(jié)構(gòu)如何分析與確定？晶格常數(shù)如何確定？已經(jīng)有6人回復(fù)
【求助】VASP怎么做晶格參數(shù)的結(jié)構(gòu)優(yōu)化已經(jīng)有17人回復(fù)
【求助】銦單質(zhì)、氧化銦、氮化銦具體的晶格結(jié)構(gòu)？不同晶格結(jié)構(gòu)對晶體形貌的影響已經(jīng)有16人回復(fù)
【求助】請教 Al2O3 所有形態(tài)的結(jié)構(gòu) 晶格常數(shù) 已經(jīng)有9人回復(fù)

點(diǎn)擊這里搜索更多相關(guān)資源

科研從小木蟲開始，人人為我，我為人人

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