Hexagonal close-packing structure on a cubic cell

N. G. Almarza
2005 Journal of Chemical Physics  
The solid phases at high pressure of atomistic systems with simple interaction models usually present close-packing structures. These structures are usually described 1 in terms of stacking of close-packed triangular lattice layers with stacking sequences ͕ABABAB͖ for the hexagonal-close-packed ͑hcp͒ structure and ͕ABCABC͖ for the face-centered cubic ͑fcc͒. The relative thermodynamic stability at given conditions of temperature and pressure of these two structures depends on the specific
more » ... the specific details of the pair potentials. The freeenergy difference between both phases uses to be small and it is often assumed that fcc is the stable crystalline phase; this is due to the fact that the fcc structure can be built up on a cubic cell, which makes the management of periodic boundary conditions ͑PBCs͒ easier in simulations, 2 whereas in the simulation of hcp phases other geometries for the simulation cells are usually chosen. 3,4 In practice, many simulation codes are restricted to work with cubic boxes. It is therefore interesting to analyze the possibility of building up a hcp structure compatible with a cubic box. Let us, firstly, remind how the fcc structure can be built in a cubic simulation box of length L. We consider systems with N =4m 3 , with N being the number of sites and m a positive integer number. The particle positions are set to be ͑k x l , k y l , k z l͒, with l = L / ͑2m͒, k x , k y , and k z being integers, that fulfill 0 ഛ k ␣ Ͻ 2m and k x + k y + k z being an even number. Within this realization of the fcc structure we can consider the different layers introduced above to be perpendicular to the ͑1, 1, 1͒ vector. We can then identify the layer of a given position, P͑x , y , z ͉ m͒, by adding up the site coordinates. P͑x,y,z͉m͒ When applying Eq. ͑1͒ to the site positions of the perfect lattice ͑as described above͒ we get integer values of P :0,1,2,.... In principle, one can apply certain displacements to all the atoms of a number of layers ͑keeping each atom in its layer plane͒ to modify the stacking sequence and attain different close-packing structures ͑e.g., ͕ABCABC͖ ↔ ͕ABABAB͖͒. However, PBCs have to be considered to apply such a strategy correctly. Let RЈ and R be the positions of the different images of a given particle: RЈ = R + n x Le x + n y Le y + n z Le z , ͑2͒ with n ␣ being integer numbers and e ␣ a unit vector on the ␣ direction. The layer identifiers P will fulfill This result implies that the system size ͑expressed in terms of m͒ must be an integer multiple of the periodicity of the stacking modification. In the current case m must be a multiple of 6. Therefore, we can use cubic cells with m =6,12,..., with fcc arrangement of sites, as starting points to build up hcp lattices; for each lattice site we can compute its layer number ͓Eq. ͑1͔͒ and perform, if required for such a layer, the translation of the particle to the corresponding new position in the hcp lattice, and eventually project back into the simulation box the new coordinates. In practice, it is possible to build up a simple algorithm to perform these operations in a simple way. For a given system size N =4m 3 ͑with m =6,12,...͒ we firstly define on the cubic cell of length L a simple cubic lattice with ͑6m͒ 3 sites located at positions given by ͑k x ␦l , k y ␦l , k z ␦l͒ ͓with k ␣ =0,1, ... ,6m − 1, and ␦l ϵ L / ͑6m͔͒. A given site will be occupied if it fulfills the following conditions: where all the operations in Eq. ͑4͒ are defined over integer numbers. In conclusion, the algorithm proposed in this work can be useful to perform simple and quick tests on the stability of different close-packing crystal structures of simple solids. Nevertheless, for a rigorous comparison between different phases, one should take into account the possible small departures of the hcp phase from its ideal structure. 3
doi:10.1063/1.1997138 pmid:16108696 fatcat:dn423wtnfzcejlnhwj6h4s44ya