Journal of Physics A.: Mathematical and General 36, 5841 (2003)

Frustration and Melting of Colloidal Molecular Crystals

C J Olson Reichhardt1 and C Reichhardt2

1 T-12, Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87545, USA
2CNLS, Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87545, USA

Received 17 October 2002, in final form 6 January 2003
Published 22 May 2003
Online at stacks.iop.org/JPhysA/36/5841
Abstract
Using numerical simulations we show that a variety of novel colloidal crystalline states and multi-step melting phenomena occur on square and triangular two-dimensional periodic substrates. At half-integer fillings different kinds of frustration effects can be realized. A two-step melting transition can occur in which individual colloidal molecules initially rotate, destroying the overall orientational order, followed by the onset of interwell colloidal hopping, in good agreement with recent experiments.
PACS numbers: 82.70.Dd, 73.50.-h
Colloidal particles are an ideal system for studying 2D ordering and melting for different kinds of substrates as the individual particle positions and dynamics can be directly visualized. Crystallization and melting on 2D periodic substrates is also relevant to vortices in superconductors with periodic pinning arrays [1,2,3] and the ordering of molecular monolayers on atomic substrates [4,5,6,7,8]. Colloidal crystallization on 2D periodic substrates has been the subject of considerable recent interest. In several recent experimental studies, 2D substrates for colloids have been created. Optical tweezers attract colloidal particles due to the refraction of the light by the colloid. Arrays of optical tweezers [9] and 2D crossed laser arrays [12] have thus been used to generate ordered substrates. Additionally, substrates with arbitrary ordering have been created by templating a pattern into the confining wall [10,11]. Colloidal ordering and melting on 2D periodic square and triangular substrates has been studied through simulation [13] and experiment [14]. It was shown that a rich variety of novel colloidal crystalline states, referred to as colloidal molecular crystals (CMC's), form at integer filling of the periodic substrates.
In Ref. [13], integer matching states up through a filling of four colloids per minima were illustrated for both square and triangular arrays, and the two-stage melting behavior of the dimer state at a filling of two colloids per minima was demonstrated for the square substrate. Here, we present the structures observed at half-integer filling for the square and triangular substrates, showing that the perfect regularity of the integer filling state is destroyed away from integer filling. We also demonstrate a two-stage melting transition for the trimer state at a filling of three colloids per minima for the square substrate.
We simulate a 2D system of Nc colloids with periodic boundary conditions in the x and y directions, using Langevin dynamics as employed in previous colloidal simulations [15]. The overdamped equation of motion for a colloid i is

dri
dt
 =fi + fs + fT.
Here fi = −∑i ≠ jNc ∇i V(rij) is the interaction force from the other colloids. The colloid-colloid interaction is a Yukawa or screened Coulomb potential, V(rij)=(Q2/|ri−rj|)exp(−κ|ri− rj| ), where Q=1 is the charge of the particles, 1/κ is the screening length, and ri(j) is the position of particle i (j). The system length is measured in units of the lattice constant a0 and we take the screening length 1/κ = a0/2. For the force from the 2D substrate, we consider both square and triangular substrates with strength A, period a0, and Nm minima. For square substrates, fs = Asin(2πx/a0)∧x + Asin(2πy/a0)∧y, and for triangular substrates, fs=∑i=13 A sin(2 πpi / a0) [cos(θi)∧x− sin(θi)∧y], where pi = x cos(θi) − y sin(θi) + a0 / 2, θ1 = π/6, θ2 = π/2, and θ3 = 5π/6. The thermal force fT is a randomly fluctuating force from random kicks. We start the system at a temperature where all the colloids are diffusing rapidly and gradually cool to T=0.0. We do not take into account hydrodynamic effects or possible long-range attractions between colloids.
Fig1.png
Figure 1: The colloid configurations (black dots) at T=0.0 for a square 2D periodic substrate with A=2.5, for a half-integer colloidal density, Nc=1.5 Nm. Every other minima captures a single colloid, while the remaining minima capture two colloids in a dimer state. There is no long-range rotational order of the dimers.
The colloidal positions for a system with a square substrate at integer matching between the colloidal periodicity and the substrate periodicity, Nc = n Nm, are illustrated in Fig. 1 of Ref. [13] for n = 1 to 4. In particular, at Nc = Nm, each colloid is located at the center of the potential minima and a square colloidal crystal forms, while at Nc = 2Nm, each minima captures two colloids which can be regarded as a colloidal dimer with a rotational degree of freedom. Over a range of substrate strengths, the colloidal dimers form a rotationally ordered state, with neighboring dimers perpendicular to one another. The orientational ordering of the dimers is due to the colloidal repulsion, and allows the distance between the colloids to be maximized under the constraint of the square substrate. In Fig. 1, we show the colloidal positions for the square substrate system at a half-integer filling, Nc = 1.5 Nm. Here we find a combination of the single colloids and dimer colloids. The single colloids are located at every other site with intervening dimers, creating a checkerboard ordering. The strong orientational ordering of the dimers observed at Nc = 2Nm is lost in this state. The dimers alternate their orientation between vertical and horizontal in a disordered fashion, forming a pattern of grain boundaries.
Fig2.png
Figure 2: The colloid configurations (black dots) at T=0.0 for a triangular 2D periodic substrate with A=2.5, for a half-integer colloidal density, Nc=1.5 Nm. It is not possible to arrange a state such that minima with only one colloid will be between every other dimer. Therefore there is no dimer ordering.
The ordered colloidal crystalline states that form on a triangular substrate at integer matching from n=1 to 4 are illustrated in Fig. 2 of Ref. [13]. As in the case of the square substrate, at Nc=Nm each minima captures a single colloid, while at Nc=2Nm each minima captures two colloids which form a dimer state. The dimers again have an additional orientational ordering in which the dimers in each row have the same orientation, which is rotated 45° with respect to the adjacent rows. This herringbone ordering resembles patterns observed experimentally for N2 monolayers adsorbed on graphite [16], or H2 molecules on a triangular lattice [17]. In Fig. 2 we illustrate the colloidal positions for the triangular substrate at half-integer filling, Nc=1.5Nm. Unlike the square substrate, where dimers could tile the lattice in an orderly way, filling every other minima, for the triangular substrate the positioning of the dimers is frustrated. It is not possible to arrange a state such that minima with only one colloid will be between every other dimer. As a result, the colloids form a very disordered state in which the dimers show a wide range of orientations, and not merely two as in the case of the square substrate at Nc=1.5Nm, or the triangular substrate at Nc=2Nm.
Fig3.png
Figure 3: The colloid positions (black dots) and trajectories (lines) over fixed time intervals at different temperatures for the trimer state on a square substrate shown in Fig. 1(c) of Ref. [13]. (a) The ordered solid phase at T/Tm0=0.25, where Tm0 is the melting temperature at zero substrate strength. (b) T/Tm0=1.5. Orientational order is destroyed as the trimers rotate within the substrate minima, but the colloids remain trapped inside each minima so the system is still in a solid phase. (c) T/Tm0=4.0. Individual colloidal diffusion occurs throughout the sample and the system is in the liquid phase.
In Fig. 3 we show the two-stage melting of the CMC at Nc=3Nm on the square lattice. As illustrated in Fig. 3(a) at low temperatures T/Tm0 < 0.25 (where Tm0 is the melting temperature at zero substrate strength), both orientational and translational order of the trimers are present and the system is frozen. This is the "ordered solid" phase. In Fig. 3(b), at T/Tm0=1.5, the trimers begin to rotate within the minima; however, diffusion of individual colloids throughout the sample does not occur. The system is still frozen but the trimer orientational order is lost. This is the "partially ordered solid" phase. In Fig. 3(c), for a higher temperature T/Tm0, the system enters a modulated liquid phase. Here the colloids begin to diffuse throughout the system. We note that in the recent experiments of Brunner and Bechinger [14], the same melting phenomena was also found for trimers on a triangular substrate. Modulated liquids have also been observed in simulations of N2 adsorbed on graphite [5].
To summarize, we have shown the rich variety of novel colloidal crystalline and disordered states that can be achieved with square and triangular two-dimensional substrates. The colloidal molecular crystal states that appear at integer filling of the substrate minima persist at half-integer filling, appearing as a mixture of the two neighboring integer filling states. For the case of Nc=1.5Nm, the single and dimer states coexist. An ordered filling of dimers can be arranged on the square substrate at Nc=1.5Nm, but the dimer arrangement is frustrated on the triangular substrate and a strongly disordered state results. The triangular lattice at Nc=1.5Nm is thus a realization of a frustrated system which would provide an interesting topic for further study. For integer filling of Nc=3Nm on a square substrate, we demonstrate the multistage melting of the CMC, where the orientational order of the colloidal molecule states is lost first, followed by the translational order. Since the colloids within a minima can act as a single particle with a rotational degree of freedom, our results also suggest that certain canonical statistical mechanics models, such as Ising, XY, Potts, and frustrated models, may be realized with colloids on two-dimensional periodic substrates. The states predicted here should be observable for colloids interacting with crossed-laser arrays or optical tweezer arrays, dusty plasmas in 2D with periodic potentials, molecules adsorbed on substrates, and vortices in superconductors with periodic substrates.

References

[1]
Baert M, Metlushko V V, Jonckheere R, Moshchalkov V V, and Bruynseraede Y 1995 Phys. Rev. Lett. 74 3269.
[2]
Harada K, Kamumura O, Kasai H, Matsuda T, Tonomura A, and Moshchalkov V V 1996 Science 271 1393.
[3]
Reichhardt C, Olson C J, Scalettar R T, and Zimányi G T 2001 Phys. Rev. B 64 144509.
[4]
Coppersmith S N, Fisher D S, Halperin B I, Lee P A, and Brinkman W F 1982 Phys. Rev. B 25 349.
[5]
Roth M and Etters R D 1991 Phys. Rev. B 44 6581.
[6]
Schick M, Walker J S, and Wortis M 1977 Phys. Rev. B 16 2205.
[7]
Ecke R E, Shu Q-S, Sullivan T S, and Vilches O E 1985 Phys. Rev. B 31 448.
[8]
Wang S C, Senbetu L, and Woo C-W 1979 Phys. Rev. B 20 3641.
[9]
Dufresne E R and Grier D G 1998 Rev. Sci. Instrum. 69 1974; Dufresne E R, Spalding G C, Dearing M T, Sheets S A, and Grier D G 2001 ibid. 72 1810.
[10]
Lin K-H, Crocker J C, Prasad V, Schofield A, Weitz D A, Lubensky T C, and Yodh A G 2000 Phys. Rev. Lett. 85 1770.
[11]
Pertsinidis A and Ling X S (unpublished).
[12]
Bechinger C et al. (unpublished).
[13]
Reichhardt C and Olson C J 2002 Phys. Rev. Lett. 88 248301.
[14]
Brunner M and Bechinger C 2002 Phys. Rev. Lett. 88 248302.
[15]
Schweigert I V, Schweigert V A, and Peeters F M 2000 Phys. Rev. Lett. 84 4381.
[16]
Diehl R D, Toney M F, and Fain S C 1982 Phys. Rev. Lett. 48 177.
[17]
Berlinsky A J and Harris A B 1978 Phys. Rev. Lett. 40 1579.


File translated from TEX by TTHgold, version 4.00.
 Back to Home