| 查看: 501 | 回复: 0 | ||
[求助]
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!" |
回复此楼» 猜你喜欢
这样的filecode谁见过
已经有5人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有13人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
fileCode有新解读？
已经有16人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
找到一些相关的精华帖子，希望有用哦~
外行人帮人求碳化氮晶体5种相的晶格结构图（要求能清楚知道轨道、能带结构、晶格）
已经有3人回复
- 使用CASTEP进行结构优化的问题
已经有25人回复
- 这种晶格结构是用什么软件做出的
已经有4人回复
- 求助超晶格体系的能带结构计算
已经有7人回复
- 如何导出高分辨晶格结构图，
已经有9人回复
- 晶格弛豫与结构弛豫
已经有18人回复
- 关于钢中马氏体的晶格结构
已经有5人回复
- 求教hcp结构的晶格矢量和k点的设置
已经有3人回复
- VASP结构优化ISIF=3，为什么晶格常数不变
已经有12人回复
- 怎样建立掺杂晶格的结构模型？
已经有18人回复
- 如何在MS里构造Co的晶格结构并计算其体系能量
已经有5人回复
- 一次型NiAl结构优化，晶格矢量变化非常大，求助~
已经有5人回复
- 求Cu3Si的晶体结构和晶格常数
已经有6人回复
- 求查询金属化合物的晶格参数,晶体结构,空间群,以及原子位置相关的网站或者书籍
已经有15人回复
- 画晶格结构小软件（绿色版）
已经有169人回复
- 各位高人！孪晶有没有可能是超晶格结构？
已经有10人回复
- 【500金币本版活动】“经验力场分子模拟+第一性原理计算问题讨论互助及资源共享”
已经有82人回复
- 【讨论】想学lammps是否应该先学习C++
已经有28人回复
- 【求助】75BiFeO3-25BaTiO3 相结构如何分析与确定？晶格常数如何确定？
已经有6人回复
- 【求助】VASP怎么做晶格参数的结构优化
已经有17人回复
- 【求助】铟单质、氧化铟、氮化铟 具体的晶格结构？不同晶格结构对晶体形貌的影响
已经有16人回复
- 【求助】请教 Al2O3 所有形态的结构 晶格常数
已经有9人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
