第一性原理研究MgZn2相的电子结构及磁性质

运用密度泛函理论的第一性原理计算分析了MgZn₂相的电子结构及相关磁性质。能带结构和态密度分析表明Zn₄s和Zn₄p轨道、Mg₃s和Mg₃p轨道分别发生sp态杂化,然后杂化态之间相互作用而形成Zn-Mg键;Mulliken布居分布计算显示:Zn₁-Mg(Zn₁是处于晶格边缘的Zn原子)和Zn₂-Mg(Zn₂是处于晶格内部的Zn原子)电子云重叠布居数接近0,电子密度分析显示Zn-Mg之间电子密度分布具有明显的定域性。结合上述结果与Zn、Mg原子的电负性差异,确定Zn-Mg键为极性共价键。分态密度(PDOS)分析显示,Zn₁-Mg键和Zn₂-Mg键的差异主要表现在Zn₂4s轨道在-10~-6eV区域对成键的贡献度高于Zn₁4s轨道,而Zn₁4s轨道在2~5eV区域对成键的贡献度高于Zn₂4s轨道。进一步对MgZn₂的积分自旋态密度和磁矩计算表明:MgZn₂磁性质表现为顺磁性,其磁性主要来源于Zn₁-Mg键中的2个自旋相同的未配对电子;MgZn₂的顺磁性特性将使Al-Zn-Mg-Cu(7×××系)高强铝合金产生磁致塑性效应。     

第一性原理研究MgZn2相的电子结构及磁性质

运用密度泛函理论的第一性原理计算分析了MgZn₂相的电子结构及相关磁性质。能带结构和态密度分析表明Zn4s和Zn4p轨道、Mg3s和Mg3p轨道分别发生sp态杂化,然后杂化态之间相互作用而形成Zn-Mg键;Mulliken布居分布计算显示:Zn₁-Mg(Zn₁是处于晶格边缘的Zn原子)和Zn₂-Mg(Zn₂是处于晶格内部的Zn原子)电子云重叠布居数接近0,电子密度分析显示Zn-Mg之间电子密度分布具有明显的定域性。结合上述结果与Zn、Mg原子的电负性差异,确定Zn-Mg键为极性共价键。分态密度(PDOS)分析显示,Zn₁-Mg键和Zn₂-Mg键的差异主要表现在Zn₂4s轨道在-10~-6 eV区域对成键的贡献度高于Zn₁4s轨道,而Zn₁4s轨道在2~5 eV区域对成键的贡献度高于Zn₂4s轨道。进一步对MgZn₂的积分自旋态密度和磁矩计算表明:MgZn₂磁性质表现为顺磁性,其磁性主要来源于Zn₁-Mg键中的2个自旋相同的未配对电子;MgZn₂的顺磁性特性将使Al-Zn-Mg-Cu(7xxx系)高强铝合金产生磁致塑性效应。     

ZnSe掺Cu与Zn空位缺陷的稳定性、电子结构与光学性质

采用基于密度泛函理论框架下的第一性原理平面波超软雁势方法, 对ZnSe闪锌矿结构本体、掺入p型杂质Cu(Zn_0.875Cu_0.125Se)及Zn空位(Zn_0.875Se)超晶胞进行结构优化处理. 计算并详细分析了缺陷体系的形成能和三种体系下ZnSe材料的态密度、能带结构、集居数、介电和吸收光谱. 结果表明: 在Zn空位与Cu掺杂ZnSe体系中, 由于空位及杂质能级的引入, 禁带宽度有所减小, 吸收光谱产生红移; 单空位缺陷结构不易形成, Zn_0.875Se结构不稳定, Cu掺杂ZnSe结构相对更稳定.     

超临界水中醋酸锌水解反应的分子动力学模拟

氧化锌(ZnO)是一种重要的化工原料, 超临界水热合成法制备纳米ZnO的第一步是锌盐与碱或水发生水解反应生成Zn(OH)₂, 后者接着脱水生成ZnO. 以Zn(CH₃COO)₂为原料, 直接和超临界水(SCW)反应能够制备纳米级的ZnO颗粒, 但对反应机理的探讨较少. 本研究利用分子动力学模拟超临界条件下Zn(CH₃COO)₂水解反应过程中的结构和能量变化, 发现Zn(CH₃COO)₂在SCW中容易聚集成无定形的团簇, 1个Zn²⁺平均和5个CH₃COO^-和1个H₂O配位, 形成6配位的八面体结构. 处于Zn(CH₃COO)₂团簇和SCW界面的Zn²⁺能够和更多的H₂O配位. 水解反应后系统的势能降低, 同时伴随Zn(CH₃COO)₂团簇结构的改变. 反应产物OH^-分布在Zn(CH₃COO)₂团簇内部, 富集Zn²⁺, 而CH₃COOH则分布在SCW中. 本文的工作为超临界水热合成的反应过程提供了基本的理论依据.     

Al 掺杂对Mg2Ni 合金的电子结构及贮氢性能的影响

采用基于密度泛函理论(DFT)的平面波赝势(PW-PP)方法, 计算Mg_2-xAl_xNiH₄ (x=0, 0.125, 0.25)合金的晶胞体积、电子态密度、键序、电荷布居、生成焓, 分析原子间成键和结构的稳定性, 研究Al 部分替代Mg 对Mg₂Ni 合金及其氢化物的结构和储氢性能的影响. 结果表明: 随着Al 含量的增加, Mg₂Ni 合金晶胞体积减小, 不利于氢原子进入合金中, 导致合金的储氢容量降低. 在Mg_2-xAl_xNiH₄ (x=0, 0.125, 0.25)中, Mg-H和Al-H相互作用远小于Ni-H的相互作用, 随着Al 含量的增加, 氢化物生成焓减少以及Ni-H的相互作用减弱, 氢化物的结构稳定性降低, Al 部分替代Mg能有效改善Mg₂Ni 合金释氢动力学性能.     

空气中合成固溶体荧光粉Ba2(Zn, Mg)Si2O7:Ce3+,Eu2+,Eu3+及其发光特性

采用高温固相法在空气中合成了Ba_1.97-yZn_1-xMg_xSi₂O₇:0.03Eu,yCe³⁺系列荧光粉。分别采用X-射线衍射和荧光光谱对所合成荧光粉的物相和发光性质进行了表征。在紫外光330~360 nm激发下,固溶体荧光粉Ba_1.97-yZn_1-xMg_xSi₂O₇:0.03Eu的发射光谱在350~725 nm范围内呈现多谱峰发射,360和500 nm处有强的宽带发射属于Eu²⁺离子的4f⁶5d¹-4f⁷跃迁,590~725 nm红光区窄带谱源于Eu³⁺的⁵D₀-⁷F_J (J=1,2,3,4)跃迁,这表明,在空气气氛中,部分Eu³⁺在Ba_1.97-yZn_1-xMg_xSi₂O₇基质中被还原成了Eu²⁺;当x=0.1时,荧光粉Ba_1.97Zn_0.9Mg_0.1Si₂O₇:0.03Eu的绿色发光最强,表明Eu³⁺被还原成Eu²⁺离子的程度最大。当共掺入Ce³⁺离子后,形成Ba_1.97-yZn_0.9Mg_0.1Si₂O₇:0.03Eu,yCe³⁺荧光粉体系,其发光随着Ce³⁺离子浓度的增大由蓝绿区经白光区到达橙红区;发现名义组成为Ba_1.96Zn_0.9Mg_0.1Si₂O₇:0.01Ce³⁺,0.03Eu的荧光粉的色坐标为(0.323,0.311),接近理想白光,是一种有潜在应用价值的白光荧光粉。讨论了稀土离子在Ba₂Zn_0.9Mg_0.1Si₂O₇基质中的能量传递与发光机理。     

Ti和Zr的取代对MgNi型储氢合金电极电化学性能的影响

使用机械合金化法成功地合成了镁基储氢合金MgNi, Mg_0.9Ti_0.1Ni,和Mg_0.9Ti_0.06Zr_0.04Ni，并对其结构和电化学性能进行了研究。X射线衍射（XRD）表明，这一系列合金的主相为非晶态，透射电镜（TEM）表明，Ti和Zr取代的合金的颗粒直径约为2~4 μm。Ti和Zr的取代提高了合金电极的循环寿命。50周充放电循环后， Mg_0.9Ti_0.06Zr_0.04Ni合金电极的放电容量高于MgNi合金电极91.74%，高于Mg_0.9Ti_0.1Ni合金电极37.96%。电极容量衰减的主要原因是在合金电极表面形成Mg的腐蚀产物Mg(OH)₂。动电位扫描结果显示，Ti和Zr的添加提高了合金电极在碱液中的抗腐蚀性能。交流阻抗(EIS)测试表明，适量Ti和Zr的添加可以明显提高合金电极的电催化活性。     

三维锌配位聚合物的合成、晶体结构及与DNA的作用

在水热条件下,由联苯-2,4,4',6-四甲酸(H₄bptc),4,4'-联吡啶(bipy),合成了3种锌配位聚合物[Zn(bptc)_0.5]_n (1),[Zn₂(bptc)(H₂O)₃]_n·nH₂O (2),[Zn₂(bptc)(H₂O)(bipy)_1.5]_n·nH₂O (3),用元素分析、红外光谱等方法对配合物的组成进行了表征,并通过单晶X-射线衍射方法测定了配合物的晶体结构.结果表明:配合物1具有双核结构,八元环金属簇Zn₂(COO)₂²⁺自组装成具有(6,6)-连接拓扑结构;配合物2具有(4,5,6)-连接拓扑结构;配合物3在辅助配体的构筑下形成三维网络结构.用溴化乙锭荧光探针法测试了配合物对EB-DNA复合体系的荧光猝灭效应,实验结果显示配合物均能使EB-DNA复合体系的荧光发生不同程度的猝灭,由此推测配合物均与DNA发生了不同程度的插入作用,引入具有刚性平面辅助配体之后的配合物3,其作用力又强于不加辅助配体的配合物1和2.     

含能配合物[Zn(DAT)6](ClO4)2(DAT＝1,5-二氨基四唑)的合成、晶体结构及性质

采用直接法合成了新型高氮含能配合物[Zn(DAT)₆](ClO₄)₂(DAT＝1,5-二氨基四唑), 并用元素分析、傅立叶变换红外光谱对其结构进行了表征. 利用缓慢蒸发溶剂法培养出其单晶, 采用X射线单晶衍射仪测定其晶体结构, 结果表明该晶体属于三方晶系, 空间群, a＝b＝1.18398(9) nm, c＝0.65700(10) nm, γ＝120°, V＝0.79760(15) nm³, Z＝1. 在目标配合物的最小不对称单元中有1个Zn^2＋, 6个DAT分子和2个. 来自6个DAT分子的6个N原子分别与中心Zn^2＋配位, 形成一个六配位、非中心对称的畸变八面体结构. 用差示扫描量热分析、热重-微分热重分析结合红外光谱研究了标题化合物的热分解机理以及分解反应动力学参数. 测定了标题配合物的感度性能, 结果表明标题配合物具有一定的摩擦感度.     

Mg0.9Ti0.1Ni1－xCox (x＝0.05, 0.1, 0.15, 0.2)合金的合成及电化学性能研究

采用机械合金化法合成了Mg_0.9Ti_0.1Ni_1－xCo_x (x＝0.05, 0.1, 0.15, 0.2)系列四元合金, 并对该系列合金的结构和电化学性能等方面进行了研究. 球磨100 h的该系列合金, XRD结果表明, X射线衍射峰均呈现宽化趋势, 基本呈非晶态. 充放电结果表明, 该系列合金具有较好的活化性能, 它们的循环稳定性明显好于MgNi合金, 其中Mg_0.9Ti_0.1Ni_0.8Co_0.2最大放电容量最高, 为427.5 mAh•g^－1. 在充放电循环过程中, Mg在合金表面形成了Mg(OH)₂是合金电极衰减的主要原因. 腐蚀曲线的测试结果表明, Co的添加可以提高合金电极在碱液中的抗腐蚀能力, 从而提高了电极的循环稳定性.     

Zn‐rich Zincides of the Ternary System Ca–Mg–Zn: Synthesis,Crystal structure,Chemical Bonding

In the course of a study on the role of magnesium in polar zincides of the heavier alkaline‐earth elements, three intermetallic phases of the ternary system Ca–Mg–Zn were synthesized from melts of the elements and their structures were determined by means of single‐crystal X‐ray data. Starting from the binary zincide CaZn₁₁, the phase width of the BaCd₁₁‐type structure reaches up to the fully ordered stoichiometric compound CaMgZn₁₀ [tI48, space group I4₁/amd, a = 1082.66(6), c = 688.95(5) pm, Z = 4, R₁ = 0.0239]. The new compound CaMgZn₅ (oP28, space group Pnma, a = 867.48(3), b = 530.37(5), c = 1104.45(9) pm, Z = 4, R₁ = 0.0385) crystallizes in the CeCu₆‐type structure, exhibits no Mg/Zn phase width and has no binary border equivalent in the system Ca–Mg–Zn. Similar to the situation in CaMgZn₁₀, one M position of the aristotype has a slightly larger coordination sphere (CN = 14) and is accordingly occupied by the larger Mg atoms. The third phase, Ca_2+xMg_6–x–yZn_15+y (hP92, space group P6₃/mmc, a = 1476.00(5), c = 881.01(4) pm, Z = 4, R₁ = 0.0399 for Ca_2.67Mg_5.18Zn_15.15) forms the hexagonal Sm₃Mg₁₃Zn₃₀‐type structure also known as μ‐MgZnRE or S phase. A small phase width (x = 0–0.67; y = 0–0.58) is due to the slightly variable Ca or Zn content of the two Mg positions. The structure is described as an intergrowth of the hexagonal MgZn₂ Laves phase and the CaZn₂ structure (KHg₂‐type). All compounds exhibit strong Zn–Zn and polar Mg–Zn covalent bonds, which are visible in the calculated electron density maps. Their structures are thus herein described using the full space tilings of [Zn₄] and [MgZn₃] tetrahedra, which are fused to polyanions consisting of tetrahedra stars, icosahedra segments etc. and the large (CN = 18–22) Ca cation coordination polyhedra. Pseudo bandgaps apparent in the tDOS are compatible with the narrow v.e./M ranges observed for other isotypic members of the three structure types.     

Alloy formation processes at electrochemical intercalation of lithium into intermetallic compounds of magnesium with zinc

A comparative study of alloy formation processes that occur during the electrochemical intercalation of lithium from lithium chloride solutions in dimethylformamide into intermetallic compounds of magnesium with zinc (MgZn₂, Mg₂Zn₃) and the corresponding individual metals is studied by chronopotentiometric and voltammetric methods. Lithium-containing phases are formed in all samples studied; moreover, for MgZn₂ and Mg₂Zn₃ electrodes, the phases formed are preferentially in the Li-Zn system. The largest number of lithium-containing phases is formed in zinc. It is shown that the electrochemical behavior of intermetallic electrodes is associated with their nature, where a single alloy component plays the key role, namely, zinc for MgZn₂ and magnesium for Mg₂Zn₃. The cathodic intercalation of lithium into MgZn₂ is characterized by anomalously low polarizability as compared with the other electrodes. The lithium extraction coefficient K _ex ^Li increases from the first to the tenth cycle for all electrode studied. The highest K _ex ^Li are typical of Zn and the lowest are typical of Mg₂Zn₃.     

Melting and homogeneous crystallization of a Lennard-Jones system

Melting and homogeneous crystallization in a Lennard-Jones system of 10,976 atoms in a model box with periodic boundary conditions were investigated by the molecular dynamics method in an NVE ensemble. Crystal melting occurs by arbitrary generation and growth of local defects transformed into regions of a disordered phase. These defects gradually span the entire space of the sample, absorbing the residual islands of crystal. Homogeneous crystallization of a liquid starts with generation of crystal nuclei which grow into defective crystals. The resulting crystal varies in structure between different realizations of the model. Face-centered cubic (fcc) structures prevail. A hexagonal close packing (hcp) structure is present on the boundaries of fcc regions and arises from disordering in alternation of atomic planes. Multiple twinning of the fcc structure is observed, and aggregates with fivefold symmetry have been found.     

Modeling of half-Heusler crystals with the chalcopyrite structure: Li<Subscript>2</Subscript>MgZn<Emphasis Type="Italic">X</Emphasis><Subscript>2</Subscript> (X = N,P, As,Sb)

A model of Li₂MgZnX ₂ half-Heusler compounds with the chalcopyrite structure is considered. The electronic structure is studied from first principles, showing that Li₂MgZnX ₂ are direct-gap crystals, except for pseudo-direct-gap Li₂MgZnP₂, with a band gap of 2.7 eV, 2.2 eV, 3.3 eV, and 2.5 eV for X = N, P, As, and Sb, respectively. The band structure and chemical bonding in the model crystals are found to be similar to those in LiMgX and LiZnX half-Heusler crystals. Total electron density and deformation electron density distributions are obtained. It is found that Mg–X and Zn–X ionic-covalent bonds are stronger than Li–X ionic bonds in Li₂MgZnX ₂ crystals, which allows Li atoms to move in the space between MgX ₄ and ZnX ₄ cation tetrahedra.     

Homogeneous crystallization of the Lennard-Jones liquid. Structural analysis based on Delaunay simplices method

The crystallization process of a simple liquid upon slow cooling has been modeled by the Monte-Carlo method. The model contains 10,000 Lennard-Jones atoms in the model box with periodic boundary conditions. The model structure is investigated at different stages of crystallization using Delaunay simplices. The simplex belonging to one or another particular crystal structure was determined by the shape of the given simplex taking into account the shape of its neighboring simplices. Simplices typical of the fcc and hcp crystal structures, as well as of polytetrahedral aggregates, not typical of crystals, were studied. The analysis has shown that the “precursors” of a hcp structure are strongly dominating over the “precursors” of a fcc structure in liquid phase before the beginning of crystallization. When crystallization starts, small embryos of the fcc structure are observed; the simplices peculiar to hcp are present at that in great amount, but they are distributed over the sample more uniformly. As crystallization proceeds, the portion of the fcc phase grows faster than hcp. However, no unified crystal appears in our case of slow cooling of the model. A complex polycrystalline structure containing crystalline regions with multiple twinning, pentagonal prisms and elements of icosahedral structures arises instead.     

The crystal structure of Mg51Zn20

The crystal structure of Mg₅₁Zn₂₀, a phase designated conventionally as “Mg₇Zn₃,” has been determined by the single-crystal X-ray diffraction method. It was solved by the examination of a Patterson synthesis, and refined by the ordinary Fourier and least-squares method; the R value obtained was 4.8% for 1167 observed reflections. The crystal is orthorhombic, space group Immm, with a = 14.083(3), b = 14.486(3), c = 14.025(3) Å, and Z = 2. There are 18 independent atomic sites, Zn1Zn6, Mg1Mg10, A, and B, and the last two sites are statistically occupied by Zn and Mg atoms with the occupancies; 0.46(2)Zn7 + 0.52(2)Mg11 and 0.24(2)Zn8 + 0.74(2)Mg12, for A and B, respectively. The structure of the crystal is described as an arrangement of icosahedral coordination polyhedra, to which all the atomic sites but Zn3 site belong. In this arrangement the Zn atoms other than the Zn3 and Zn8(B) center the icosahedral coordination polyhedra with coordination number 12. The Zn3, Zn8 atoms, and all the Mg atoms except Mg11(A) are located at the centers of various coordination polyhedra with the coordination numbers from 11 to 15. The distances between neighboring atoms are 2.71–3.07, 2.82–3.65, and 2.60–3.20 Å for ZnZn, MgMg, and ZnMg, respectively.     

Geometric and topological analysis of icosahedral structures of samson Mg2Zn11 (cP39) Phases, K6Na15Tl18H (cP40), and Tm3In7Co9 (cP46): Nanocluster precursors, self-assembly mechanism, and superstructure ordering

Geometric and topological analysis and 3D reconstruction of self-assembly of icosahedral structures of Samson Mg₂Zn₁₁ clusters (space group Pm[`3]Pm\bar 3, cP39, 10 compounds) and the K<sub>6</sub>Na<sub>15</sub>Tl<sub>18</sub>H and Tm<sub>3</sub>In<sub>7</sub>Co<sub>9.29</sub> structures were performed by computer methods (the TOPOS program package). The complete decomposition of the 3D graph of the crystal structures into cluster substructures showed the existence of the crystal-forming nanocluster precursor A comprising 45 atoms (A-45). The S-6 cluster spacers were identified in Mg<sub>2</sub>Zn<sub>11</sub>, and the S-7 cluster spacers were found in K<sub>6</sub>Na<sub>15</sub>Tl<sub>18</sub>H. In Tm<sub>3</sub>In<sub>7</sub>Co<sub>9.29</sub>, the S-6 and S-7 cluster spacers with the centers statistically occupying the same position were determined. The A-45, S-6 (octahedron), and S-7 (centered octahedron) clusters have symmetry [`3]m\bar 3m. The A-45 nanocluster contains an inner Zn(Zn)₁₂ template icosahedron and an external quasi-spherical shell composed of 32 atoms (deltahedron D32). A-45 is equivalent to the Bergman cluster used as the approximant of the local structure of quasicrystals. For deltahedron D32, the existence of a hierarchical structure was identified as a result of self-assembly involving two types of cyclic clusters: K-7 with an atom in the center of the sixth ring and three-atom cyclic clusters K-3. The atoms of the K-3 and K-7 clusters occupy all possible positions over the 12 vertices and 20 faces of an icosahedron and thereby form an edge net of bonds made of triangles. For the K₆(Na₁₄MTl₁₈) structures (M = Mg, An, Cd, Hg), the cluster nature of superstructure ordering of three chemically different atoms (14Na, M, and 18Tl) over 33 positions of the Zn atoms in the unit cell of the basis Mg₂Zn₁₁(Mg₆Zn₃₃) structure was considered.     

Effect of agar hydrogel and magnesium ions on the biomimetic mineralization of calcium carbonate

Here, agar hydrogel was selected as diffusion medium and template to control the biomimetic mineralization of calcium carbonate (CaCO₃). Due to three dimensional network structures and abundant functional groups (such as, hydroxyl groups), Ca²⁺ ions were uniformly distributed in the network and electrostatically attracted. The diffusion speed and range of CO₃^2? ions were mediated by the concentration of hydrogel medium. Under the synergistic effect of Mg²⁺ ions, the crystal CaCO₃ was induced by gas phase diffusion method in the hydrogel system. The results showed that the concentrations of Mg²⁺ ions and agar hydrogel had no obvious effect on the calcite phase of CaCO₃, but the morphologies and sizes changed with concentrations of medium and Mg²⁺ ions. Attribute to template effect, the crystallization behavior and growth rate of CaCO₃ crystals were regulated. Since Mg²⁺ ions were easily adsorbed on the surfaces of unit cell, the unique structure of CaCO₃ was precisely controlled. This study provides a useful reference and inspiration for the understandings of the contributions of ion supply rate in bio-mineralization and hydrogel medium in biomimetic mineralization.     

Phase structure and morphology of zinc-iron alloy electrodeposits

Zinc-iron alloy electrodeposits are providing higher corrosion resistance to steel components and also having better mechanical properties when compared to zinc deposits. This is due to the unique phase structure of the alloy formed. This study elucidates the phase structure of the electrodeposited alloy, based on the deposition kinetics and morphological characteristics. Deposition of iron was hindered by charge-transfer process, at low current densities. But zinc deposition was prevailed through diffusion control, only at high current densities. The probability of substitution of iron in hcp lattice along c-axis is more, than a-axis. This is because the linear density along c-axis is lower than a-axis. Intermetallic compounds of variable compositions were identified. Compounds such as FeZn₁₆, FeZn₁₃, Fe₅Zn₇₈ have dominantly “η” phase structure and FeZn₆, Fe₅Zn₂₂, Fe₂Zn₇, Fe₅Zn₂₉, Fe₃Zn₁₃ have “Γ” phase structure.     

Polymorph selection in the crystallization of hard-core Yukawa system

Colloid-colloid interactions in charge-stabilized dispersions can to some extent be represented by the hard-core Yukawa model. The crystallization process and polymorph selection of hard-core Yukawa model are studied by means of smart Monte Carlo simulations in the region of face-centered-cubic (fcc) phase. The contact value of hard-core Yukawa potential and the volume fraction of the colloids are fixed, while the Debye screening length can be varied. In the early stage of the crystallization, the precursors with relatively ordered liquid structure have been observed. Although the crystal structure of thermodynamically stable phase is fcc, the system crystallizes into a mixture of fcc and hexagonal close-packed (hcp) structures under small Debye screening length since the colloidal particles act as effective hard spheres. In the intermediate range of Debye screening length, the system crystallizes into a mixture of fcc, hcp, and body-centered-cubic (bcc). The existence of metastable hcp and bcc structures can be interpreted as a manifestation of the Ostwald’s step rule. Until the Debye screening length is large enough, the crystal structure obtained is almost a complete fcc suggesting the system eventually reaches to a thermodynamically stable state.