Hexagonal to square lattice conversion in bilayer systems

We report the results of extensive molecular dynamics simulations of the reconstructive hexagonal to square lattice conversion in bilayer colloid systems. Two types of interparticle potential were used to represent the colloid-colloid interactions in the suspension. One potential, due to Marcus and Rice, is designed to describe the interaction of sterically stabilized colloid particles. This potential has a term that represents the attraction between colloid particles when there is incipient overlap between the stabilizing brushes on their surfaces, a (soft repulsion) term that represents the entropy cost associated with interpenetration of the stabilizing brushes, and a term that represents core-core repulsion. The other potential we used is an almost hard core repulsion with continuous derivatives. Our results clearly show that the character of the reconstructive hexagonal to square lattice conversion in bilayer colloid systems is potential dependent. For a system with colloid-colloid interactions of the Marcus-Rice type, the packing of particles in the square array exhibits a large interlayer lattice spacing, with the particles located at the minima of the attractive well. In this case the hexagonal to square lattice transition is first order. For a system with hard core colloid-colloid interactions there are two degenerate stable intermediate phases, linear and zigzag rhombic, that are separated from the square lattice by strong first order transitions, and from the hexagonal lattice by either weak first or second order transitions.