24小时热门版块排行榜    

查看: 449  |  回复: 0

digghost

木虫 (著名写手)

[求助] minimize 为什么不能得到其他的晶格结构

下面的in文件本来是fcc结构作优化,得到稳定值4.05,我将fcc改成了bcc,结果优化后结构还是bcc,这是为什么?
然后我将温度升高到4000 K使体系完全熔化,在降温到100 K 结构还是bcc结构的。这是怎么回事?难道优化还和设定的结构有关?

# 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!"
回复此楼

» 猜你喜欢

呵呵了
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子,希望有用哦~

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 digghost 的主题更新
信息提示
请填处理意见