Physical Review E 68, 060401(R) (2003)

Structure and Melting of Two-Species Charged Clusters in a Parabolic Trap

J.A. Drocco1, C.J. Olson Reichhardt2, C. Reichhardt2, and B. Jankó1

1 Department of Physics, University of Notre Dame, Notre Dame, Indiana 46556
2 Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545

(Received 12 September 2003; published 24 December 2003)

We consider a system of charged particles interacting with an unscreened Coulomb repulsion in a two-dimensional parabolic confining trap. The static charge on a portion of the particles is twice as large as the charge on the remaining particles. The particles separate into a shell structure with those of greater charge situated farther from the center of the trap. As we vary the ratio of the number of particles of the two species, we find that for certain configurations, the symmetry of the arrangement of the inner cluster of singly-charged particles matches the symmetry of the outer ring of doubly-charged particles. These matching configurations have a higher melting temperature and a higher thermal threshold for intershell rotation between the species than the nonmatching configurations.
Clusters of repulsive particles in confined traps have attracted considerable recent attention due to their applicability to a wide variety of systems. For example, two-dimensional (2D) clusters can represent electrons in quantum dots [1] or on the surface of liquid helium [2], vortices in superfluids [3], colloidal particles in circular traps [4], confined ferromagnetic particles [5], and charged dust particles in plasma traps [6]. The 2D charged clusters also resemble the problem of charge distribution studied by Thomson in his "plum-pudding" model of the atom [7].
When confined to a parabolic trap, charged particles form a structure of concentric rings, with the inner particles forming a distorted triangular lattice resembling a defected Wigner crystal, and the outer rings taking on a more circular shape that conforms to the radial symmetry of the trap [8,9]. Among the possible charge configurations are several "magic" arrangements in which the number of particles is such that the shells form with only a few symmetrically distributed dislocations, and so have a reduced total energy compared to what is predicted based on a semi-empirical approximation [10]. Bubeck et al. [4] have observed that certain colloidal clusters confined by a circular hard-wall trap exhibit two-stage melting. As the temperature is increased in their experiment, intershell rotation between the outer two shells first occurs, destroying the angular order of the system. Angular order is then restored, until at higher temperatures complete melting occurs.
Several explanations of this two-stage melting phenomenon have been proposed [8,9,11,12], all of which focus on the intershell rotation which occurs prior to the exchange of particles between shells. Most plausible among these is the theory that the rotation is due to an incommensuration between the shapes of the potentials created by the adjacent shells. For this intershell rotation to occur, the inner shell configuration must be sufficiently stable to have a melting temperature higher than the threshold for intershell rotation.
In our simulation, we extend the confined charge system to include particles with two distinct values of charge. We find that the two species separate into shells, with those of greater charge located farther from the center of the trap. Since the particles have a long range repulsive interaction, the lowest energy configuration places the stronger charges as far apart as possible. If two-stage melting and intershell rotation is caused by an incommensuration between the potentials formed by the particles, then we should be able to predict the occurrence of intershell rotation based on the ratio of particles in the outer shell of singly-charged particles to the number of doubly-charged particles.
We consider a system of Ns + Nd charged particles interacting via an unscreened 1/r Coulomb repulsion, where Ns is the number of single-charge particles with qs=1 and Nd is the number of doubly-charged particles with qd=2. The particles are free to move in two-dimensions but are confined by a parabolic trap centered at the origin and increasing radially as r2. The dimensionless Hamiltonian [9] for this system is:
H= Ns+Nd

i=1 
Ns+Nd

j=i+1 
qiqj

|

r
 

i 

r
 

j 
|
+A( Ns+Nd

i=1 
|

r
 

i 
|2)
(1)
where qi (ri) is the charge (position) of particle i and we fix the strength of the parabolic trap to A=10. Charged colloidal particles in a confinement potential, such as those of Wei et al. [13], bear the closest resemblance to the Coulomb interaction used in this study.
Using a molecular dynamics (MD) simulation method, we initialize the system at high temperature, simulated by random Langevin kicks, and then slowly anneal it to a T=0 ground state configuration. We checked the accuracy of this method by reproducing the ground state configurations for single-species particle clusters found by Kong et al.[14], albeit with different confinement strength. Once we have obtained the ground states of the two-species clusters, we slowly increase the temperature and observe the melting of the system.
Fig.  1 shows several examples of the ground-state configurations we obtained for Nd=5, 6, and 7, respectively. The structure of the inner core of singly-charged particles is altered both by the presence and by the number of doubly-charged particles that are present. For Nd=5, configurations with Ns=6, 11, and 31 have five-fold symmetry matching the number of doubly-charged particles in the outer shell. Ns=8 approximates a six-fold symmetric arrangement by substituting a singly-charged particle to fill a potential well in the outermost shell. For Nd=6, configurations with Ns=7 and Ns=19 have six-fold symmetry. Configurations with Ns=12 and Ns=27 have three-fold symmetry and configuration Ns=30 has two-fold symmetry, which also are commensurate with the arrangement of particles in the outer shell. For Nd=7, Ns=8 and 22 are the only seven-fold symmetric inner particle configurations. Ns=23 and 24 are nearly symmetric, with seven-fold symmetry in the outer two shells of the inner particle configuration.
In general we find highly ordered structures when the smaller particles form a commensurate structure with the outer particles. All the configurations are shown at [15]. The general expressions Ns=kNd and Ns=kNd+1, with k a small integer, predict some configurations that have rotational symmetry through an angle 2π/Nd, such as Nd=5, Ns=6, but do not hold for all values of k due to the fact that the outer particles may distort from a uniformly spaced arrangement in order to better accommodate the inner particles. For example, at Nd=5, Ns=8, [Fig. 1(b)] a small charge moves to the outer ring. This can be simply understood in terms of angular stiffness of the rings. The radius of the outer ring increases as Ns increases, but since Nd is fixed, the large charges move further apart, decreasing the angular stiffness of the outer ring and making it susceptible to angular compressions that permit the intrusion of a small charge. A similar but opposite process occurs for the inner ring of small charges. As Ns increases, the number of small charges around the ring increases until the angular stiffness of the ring is so large that it is not possible to compress it enough to add another small charge, and the next charge must move out into a new ring.
Although in this work we only present results for qd=2qs, we have considered other charge ratios, and find that the more highly charged particles always move to the outside of the trap, and the same general commensurate-incommensurate phases occur. We have also performed simulations with larger Nd > 7; however, we do not treat this case here since the outer particles begin to form multiple rings which have their own commensurate-incommensurate transitions. In addition, we have run a set of simulations for Ns < 15 with the same trap but with a 1/r3 interparticle interaction potential, which corresponds to the interaction between magnetic colloids. We found that, in general, the ground states and also the dynamical properties were qualitatively identical to those in the case of Coulomb repulsion for small Ns (Ns < 15).
Fig1.png
Figure 1: Ground-state colloid configurations for: (a-d) Nd=5, (e-h) Nd=6, and (i-l) Nd=7, with Ns = (a) 6, (b) 8, (c) 11, (d) 31, (e) 7, (f) 12, (g) 19, (h) 30, (i) 8, (j) 19, (k) 22, (l) 24. Small dots are the singly charged particles and large dots are the doubly charged particles.
Fig2.png
Figure 2: Plot of σs vs. temperature for Nd=6, Ns=7, averaged over five realizations.
Fig3.png
Figure 3: Particle trajectories (lines) for Nd=6, Ns=7 at temperatures (a) T=3.5 with no intershell exchange and (b) T=3.6 with intershell exchange.
In order to confirm whether the highly ordered commensurate phases are more stable we consider the melting of the two-species system. We determine the temperature Te of the first exchange of particles between shells by measuring intershell exchange of the singly charged particles,
σs= 1

Ns
Ns

i=1 
|

r
 

i 
(t)−

r
 

i 
(0)|
(2)
This gives the mean radial distance of the inner singly-charged particles from their initial T=0 positions. The slow increase in σs at low temperatures corresponds to the small radial distortions of the configuration introduced by the random Langevin kicks, which increase in size as the temperature increases. The overall configuration remains fairly stiff at these low temperatures. σs shows a marked increase at Te when the inner particles begin to jump between shells, as illustrated in Fig. 2 for a system with Ns=7 and Nd=6. This configuration is highly stable, as shown in Fig. 1(e), with one central particle surrounded by an inner hexagonal shell of six singly-charged particles and an outer hexagonal shell of six doubly-charged particles. At T=3.55, σs jumps when one of the first-shell particles exchanges with the central particle, as shown in the trajectories of Fig. 3.
A similar measure, σd, tracks the doubly-charged particles; however, as they never formed more than one shell for the parameters considered here, there were no intershell exchanges and σd only increased linearly with temperature. We observe exchanges of a singly-charged particle with a doubly-charged particle only at very high temperatures T >> Te, and these exchanges occur only for highly asymmetric and disordered configurations.
To track the onset temperature Tr of intershell rotation between the two species when it occurs, we use a second measure, ∆θ:
θ=| 1

Ns
Ns

i=1 
i(t)−θi(0))− 1

Nd
Nd

j=1 
j(t)−θj(0))|  ,
(3)
which gives the difference between the mean angular displacements of the two species from their initial configurations. ∆θ increases when the shells slip past each other, but is insensitive to coherent rotation of the two species. ∆θ becomes meaningless if the particles do not maintain the same orientation with respect to the other particles of the same species, so it can detect intershell rotation only when this occurs before the onset of intershell exchange. ∆θ also detects relative slip between shells of the same particle species, which does not interest us here. Such same-species shell slips were generally limited to the erratic rotation of two particles at the center of the configuration, and produced a sufficiently small change in ∆θ to be distinguished easily from a genuine rotation relative to the other species. For configurations with a single central particle located roughly at the origin, the angular displacement of that particle was excluded from ∆θ in order to reduce noise.
Fig4.png
Figure 4: Plot of ∆θ vs. temperature for Nd=7, Ns=7, averaged over five realizations.
Fig5.png
Figure 5: Particle trajectories over a short time interval for Nd=7, Ns=7 at (a) T=0.6 with no intershell rotation and (b) T=0.7 with intershell rotation.
An example of intershell rotation measured by ∆θ is shown in Fig.  4, where we plot ∆θ vs temperature for a configuration of Ns=7 and Nd=7. In this case, the configurations of each of the two species are highly stable independently, with the singly-charged particles forming a hexagonal ring around a central particle. However, there is an incoherence between the hexagonal ring of singly-charged particles and the seven-particle ring of doubly charged particles which surrounds it. Since the onset of intershell rotation at Tr is more gradual than the intershell exchange transition, we define Tr to occur when ∆θ exceeds a threshold value of ∆θ=[(π)/(Nd)], corresponding to an angular displacement of half of the angular distance between neighboring particles in the outer shell of doubly-charged particles. In Fig. 4, ∆θ=[(π)/(Nd)] falls at T ≈ 0.65, corresponding to a slipping between the outer shell and the inner cluster of particles. The trajectory plot, Fig.  5, also indicates this slipping, as the trails of the inner shell particles become larger than those of the outer shell particles at this temperature.
Fig6.png
Figure 6: Melting temperature Tm vs Ns for (a) Nd=5, (b) Nd=6, and (c) Nd=7.
The matching of the inner particle symmetry with the outer shell symmetry produces an elevated melting threshold for commensurate configurations. The melting temperature Tm is taken to be the lower of Te or Tr. In Fig. 6(b) we plot Tm for a range of Ns at fixed Nd=6. Figs. 6(a) and (c) show Tm for Nd=5 and Nd=7, respectively. For each configuration, we averaged σs and ∆θ over five realizations to reduce error. In Fig. 6(b), Tm for Nd=6 has a peak at Ns=19, another perfect triangular arrangement. We also find an elevated melting temperature for Nd=12, which has threefold symmetry. As seen in Fig. 6(a), configuration Nd=5 and Ns=11 also exhibits a higher than average melting temperature. This is expected since the configuration Ns=11 has five-fold symmetry, commensurate with the number of outer particles. Fig. 6(c) shows high Tm values for Nd=7 at Ns=23 and Ns=24, which is expected due to the matching seven-fold symmetry of the outer three rings. The inner part of the configuration Nd=7, Ns=19 forms the same highly stable perfectly hexagonal arrangement of Nd=6, Ns=19, and thus has a high threshold for intershell exchange between the singly-charged inner particles. However, the incommensuration with the 7 particles on the outside yields a low threshold for intershell rotation.
There is an obvious shift in the peaks of Fig. 6(a-c) as the number of outer-shell doubly-charged particles changes. Most notably, the highest melting temperature for each Nd occurs for Ns=Nd+1, when the singly-charged particles form a central particle surrounded by a shell that is commensurate with the outer particles. As expected, the highest melting temperature occurs for Nd=6 and Ns=7, as this arrangement matches the triangular lattice of the Wigner crystal. However, melting temperatures nearly as high occur for Nd=5 and Ns=5 and 6. We believe that the smaller size of the crystal in this case exaggerates the effect of the commensuration between rings, despite the dissimilarity with the bulk lattice configuration.
In conclusion we have investigated the structure and melting of two species of charged particles in a parabolic trap. The more highly charged particles form an outer ring. Highly ordered clusters occur when the structure of the central cluster of singly charged particles is commensurate with the outer ring of doubly charged particles. We observe variations in the melting temperatures of the two species clusters, with elevated melting thresholds for intershell rotation and intershell particle exchange when a commensuration occurs between the symmetry of the inner cluster and the outer ring. It may be possible to observe similar variations in melting for other confining potentials, including hard wall potentials, and this would be an interesting direction for further research.
This work was supported by the U.S. Department of Energy under Contract No. W-7405-ENG-36 and by the U.S. DoE, Office of Science, under Contract No. W-31-109-ENG-38. BJ and JD were supported by NSF-NIRT award DMR02-10519 and the Alfred P. Sloan Foundation.

References

[1]
M.A. Reed and W.P. Kirk, Nanostructure Physics and Fabrication (Academic Press: Boston, 1989).
[2]
P. Leiderer, W. Ebner, and V.B. Shilkin, Surf. Sci. 113, 105 (1987).
[3]
Y. Kondo, J.S. Korhonen, M. Krusius, V.V. Dmitriev, E.V. Thuneberg, and G.E. Volovik, Phys. Rev. Lett. 68, 3331 (1992).
[4]
R. Bubeck, C. Bechinger, S. Neser, and P. Leiderer, Phys. Rev. Lett. 82, 3364 (1999).
[5]
M. Golosovsky, Y. Saado, and D. Davidov, Phys. Rev. E 65, 061405 (2002).
[6]
W-T. Juan, Z-H. Huang, J-W. Hsu, Y-J. Lai, and Lin I, Phys. Rev. E 58, R6947 (1998).
[7]
J.J. Thomson, Philos. Mag. 7, 237 (1904); B. Partoens and F.M. Peeters, J. Phys.: Condens. Matter 9, 5383 (1997).
[8]
V.M. Bedanov and F.M. Peeters, Phys. Rev. B 49, 2667 (1994).
[9]
V.A. Schweigert and F.M. Peeters, Phys. Rev. B 51, 7700 (1995).
[10]
A.A. Koulakov and B.I. Shklovskii, Phys. Rev. B 57, 2352 (1998).
[11]
I.V. Schweigert, V.A. Schweigert, and F.M. Peeters, Phys. Rev. Lett. 84, 4381 (2000).
[12]
Y-J. Lai and Lin I, Phys. Rev. E 64, 015601 (2001).
[13]
Q.-H. Wei, C. Bechinger, D. Rudhardt, and P. Leiderer, Phys. Rev. Lett. 81, 2606 (1998).
[14]
M. Kong, B. Partoens, and F.M. Peeters, Phys. Rev. E 65, 046602 (2002).
[15]
All configurations can be seen at http://www.t12.lanl.gov/home/olson/Colloid2S.html.



File translated from TEX by TTHgold, version 4.00.

Back to Home