Next Article in Journal
Kinetic Energy of a Free Quantum Brownian Particle
Next Article in Special Issue
Hierarchy of Relaxation Times and Residual Entropy: A Nonequilibrium Approach
Previous Article in Journal
Stochastic Proximal Gradient Algorithms for Multi-Source Quantitative Photoacoustic Tomography
Previous Article in Special Issue
Glass Transition, Crystallization of Glass-Forming Melts, and Entropy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Attraction Controls the Entropy of Fluctuations in Isosceles Triangular Networks

1
School of Mechanical Engineering, Tel-Aviv University, Tel-Aviv 69978, Israel
2
Sackler Center for Computational Molecular and Materials Science, Tel-Aviv University, Tel-Aviv 69978, Israel
*
Author to whom correspondence should be addressed.
Entropy 2018, 20(2), 122; https://doi.org/10.3390/e20020122
Submission received: 26 January 2018 / Revised: 8 February 2018 / Accepted: 10 February 2018 / Published: 12 February 2018
(This article belongs to the Special Issue Residual Entropy and Nonequilibrium States)

Abstract

:
We study two-dimensional triangular-network models, which have degenerate ground states composed of straight or randomly-zigzagging stripes and thus sub-extensive residual entropy. We show that attraction is responsible for the inversion of the stable phase by changing the entropy of fluctuations around the ground-state configurations. By using a real-space shell-expansion method, we compute the exact expression of the entropy for harmonic interactions, while for repulsive harmonic interactions we obtain the entropy arising from a limited subset of the system by numerical integration. We compare these results with a three-dimensional triangular-network model, which shows the same attraction-mediated selection mechanism of the stable phase, and conclude that this effect is general with respect to the dimensionality of the system.

1. Introduction

Geometrically-constrained systems may show peculiar features compared to their unconstrained counterparts. In particular, geometric constraints can lead to frustration because the system cannot simultaneously minimize all local interaction energies, or free energies. Frustrated systems usually show a degenerate ground state and then they may posses a residual entropy. Frustration is relevant in physical and biological systems that range from water [1] and spin ice [2,3] to magnets [4], magnetic island [5], high transition-temperature superconductors [6], elastic beams [7,8] and colloids [9,10,11,12]. The possibility to control colloidal interparticle interactions and to visualize and manipulate each particle and follow its motion in both space and time makes colloidal suspension a powerful tool to study phenomena in condensed-matter physics, ranging from glass formers [13], to crystals and gels [14].
A prototypical geometrically-confined system is composed of short-range repulsive colloids confined in a slit pore of two plates. Varying density and plate separation, discontinuous phase transitions between layered, buckled, rhombic and adaptive prism crystal structures occur [15,16,17]. In the case of a slit pore with a plates interdistance slightly larger than a colloid diameter, when density approaches the close-packing value ρ c p , colloids, due to their free-volume-dominated free energy, tend to touch opposite walls, giving rise to effective antiferromagnetic interactions [9,18], and to glassy dynamics [19]. Multiple configurations corresponding to the same ρ c p can be obtained by alternating straight stripes of up and down spheres (Figure 1a) or by any set of zigzagging stripes (Figure 1b).
This ground-state degeneracy implies a subextensive residual entropy at ρ c p ( S 0 N , with N the number of particles in the system) [9] so that the residual entropy per particle tends to zero in the thermodynamic limit. At ρ = ρ c p , in the straight- or zigzagging-stripes configuration, a colloidal sphere is surrounded by two colloids touching the same wall (giving rise to frustrated bonds) and by four colloids touching the opposite wall (satisfied bonds), except for the limit case of α = 90 , in which case every sphere is surrounded by four colloids touching the same wall and four colloids touching the opposite wall, where α is defined in Ref. [20]. The terms frustrated and satisfied bonds refer both to the packing consideration of neighboring spheres wanting to touch opposite walls, and also from the obvious connection to an antiferromagnetic spin-1/2 Ising model [18]. In this case, nearest neighbors at opposite spin states satisfy the antiferromagnetic interaction between them, while those at the same spin state represent the frustration of this Ising model on the triangular lattice [21]. The three-dimensional (3D) network of links connecting the centers of neighboring spheres is composed of tilted equilateral triangles, since the 3D distance between the centers of each pair of contacting spheres is equal to the sphere diameter. When the centers of the colloids are projected on the plane parallel to the plates, it reduces to a 2D network composed of isosceles triangles. In this network, the longer link of each triangular plaquette corresponds to a frustrated bond, while the two shorter links correspond to satisfied bonds. Straight and zigzagging stripes are the only configurations corresponding to a 2D network of isosceles triangles which can tile the plane [9].
This geometric mechanism underlying the ground state of the buckled colloidal system composed of straight or zigzagging stripes, has been realised, experimentally also by packing a granular system in a container under the effect of gravity [22,23], and theoretically with spins starting from the Wannier antiferromagnetic Ising model on a triangular lattice [21] by allowing for the lattice to elastically deform [24,25]. Zigzagging stripes patterns have been found also in the Ising model on an anisotropic triangular lattice [26]. A relevant phenomenon involving the selection of one between two competing phases related to the distorsion of equilateral into iscosceles triangles and due to angular frustrations concerns stiff, nematically ordered, polymer molecules such as DNA [27]. The Ising antiferromagnet on a deformable triangular lattice has the same degeneracy of the ground state and subextensive entropy at T = 0 as the colloidal system at ρ = ρ c p (that is S 0 N , equivalent to the N 2 / 3 scaling found in Perovskite Oxynitrides [28]). In the following we refer to T = 0 as temperature being arbitrarily close to zero. Indeed, it has been shown that the third law of thermodynamics, implying the unattainability of absolute zero temperature in a finite number of steps and within a finite time, holds for arbitrary classical or quantum systems or involving infinite-dimensional reservoirs [29].
For temperature slightly larger than zero, the degeneracy is removed by the order-by-disorder effect [30,31,32,33,34,35,36,37] and in the elastic Ising model straight stripes represent the stable phase [20,25], while bent stripes are selected in the colloidal system for ρ < ρ c p when colloids are modeled using hard or soft repulsive potentials [20]. By tuning the attractive vs repulsive components of an asymmetric power-law potential used to model colloids, we found that the stable phase in the colloidal system can be turned from bent to straight stripes for attraction strong enough compared to repulsion [20]. We established a connection between the effect of the attraction on the phase stability and the packing of hard spheres and their entropy. We showed that other parameters of these systems are irrelevant to the phase stability, as for example their dimensionalities. Indeed, the elastic Ising antiferromagnet is defined in 2D, while the colloidal system is 3D or quasi-2D due to the buckling of the monolayer.
In this paper we study a 2D isosceles triangular network which shows the same ground-state degeneracy as the Ising elastic antiferromagnet at T = 0 or the colloidal monolayer at ρ = ρ c p . Increasing temperature above zero, the degeneracy is removed through the order-by-disorder effect and we find the straight-stripes phase to be selected for particles linked with harmonic interactions, while the bent-stripes phase is more stable if only the repulsive component of the harmonic inter-particle potential is considerd (that we call repulsive harmonic potential, described in detail in the following section). This result suggests that the inversion of the stable phase by adding attraction to repulsively-interacting particles in triangular networks is a general mechanism, irrespective of the dimensionality of the system.

2. Isosceles Triangular Network Model

The model we consider is composed of particles in a 2D triangular network linked with springs of two different lengths at rest such that every plaquette or triangle is formed by a longer edge 2 a and two shorter identical edges c, and thus has a height b = c 2 a 2 and a head angle α such that a = b tan ( α / 2 ) (see Figure 1e). We will first consider particles interacting through a harmonic potential and then will consider a repulsive harmonic potential, as defined below. The Hamiltonian of the system for harmonic inter-particle interaction can be written as
H h = m , n l = 1 3 K 2 ( d r l d r 0 ) 2 ,
where K is the spring constant, which is assumed to be identical for all springs, 1 m , n L , N = L 2 is the number of particles or nodes of the network and the index l runs over three of the six neighbors each particle has to avoid double counting of bonds. The positions of the particles are described by the coordinates { x i , y i } . At T = 0 multiple degenerate states, represented by straight or zigzagging stripes, minimize the system free energy. For T > 0 the system cannot jump from one configuration to another, but at T = 0 it is at mechanical equilibrium in every state under consideration. Therefore, the present model is not ergodic by construction. In the conclusion we will discuss entropy calculation in non ergodic systems. In order to study the stability of the system at T > 0 , we will consider small fluctuations about the equilibrium position described by small displacements { u i , v i } of all particles in the straight and bent configurations. d r is the distance between particles i and j, and its square is thus given by
d r 2 = ( d x + d u ) 2 + ( d y + d v ) 2 = d r 0 2 + 2 ( d x d u + d y d v ) + d u 2 + d v 2
where d x = x i x j , d y = y i y j , d u = u i u j , d v = v i v j and d r 0 = ( d x 2 + d y 2 ) 1 / 2 is the length at rest of the spring linking particles i and j, which can take the values c or 2 a (see Figure 1e, Table 1 and Table 2). Since we consider the expansion around mechanical equilibrium, we will drop terms linear in d u and d v and write: d r 2 = d r 0 2 + d u 2 + d v 2 . Expanding to harmonic order the expression of d r that we get from Equation (2), we obtain
d r = d r 0 + d u 2 2 d r 0 1 d x 2 d r 0 2 + d v 2 2 d r 0 1 d y 2 d r 0 2 d x d y d u d v d r 0 3 .
The Hamiltonian of the straight-stripes configuration, with particle positions specified for the unit-cell in Table 1, is
H s h = K 2 m , n d u 1 2 + a 2 c 2 ( d u 2 2 + d u 3 2 ) + b 2 c 2 ( d v 2 2 + d v 3 2 ) 2 a b c 2 ( d u 2 d v 2 + d u 3 d v 3 ) .
Using the relations d u l = u l u 0 and d v l = v l v 0 we get
H s h = K m , n u 0 2 u 0 u 1 + a 2 c 2 ( 2 u 0 2 u 0 u 2 u 0 u 3 ) + b 2 c 2 ( 2 v 0 2 v 0 v 2 v 0 v 3 ) 2 a b c 2 ( u 0 v 2 + u 2 v 0 u 0 v 3 u 3 v 0 ) .
The Hamiltonian of the bent-stripes configuration, with particle positions specified for the unit-cell in Table 2, is
H b h = K 2 t , n d u 10 2 + cos 2 α d u 50 2 + sin 2 α d v 50 2 sin ( 2 α ) d u 50 d v 50 + sin 2 ( α 2 ) ( d u 40 2 + d u 41 2 + d u 31 2 ) + sin 2 ( 3 α 2 ) d u 21 2 + cos 2 ( α 2 ) ( d v 40 2 + d v 41 2 + d v 31 2 ) + cos 2 ( 3 α 2 ) d v 21 2 sin α ( d u 40 d v 40 + d u 41 d v 41 d u 31 d v 31 ) + sin ( 3 α ) d u 21 d v 21 ,
where 1 t L / 2 . Indeed, the unit-cell of the bent stripe configuration includes two particles (see Figure 1d): particle 0, which represents the particles with odd m, and particle 1, which represents the particles with even m. Therefore, we set for particle 0, m = 2 t 1 , and for particle 1, m = 2 t . Using the relations d u l 0 = u l u 0 , d u l 1 = u l u 1 and d v l 0 = v l v 0 , d v l 1 = v l v 1 we get
H b h = K 2 t , n 2 sin 2 α + 3 sin 2 ( α 2 ) + sin 2 ( 3 α 2 ) ( u 0 2 + u 1 2 ) + ( sin α sin 2 α + sin 3 α ) ( u 0 v 5 + u 5 v 0 ) sin α ( u 0 v 4 + u 4 v 0 u 1 v 4 u 4 v 1 + u 1 v 3 + u 3 v 1 ) sin 3 α ( u 1 v 2 + u 2 v 1 ) + 1 cos 2 α + 3 cos 2 ( α 2 ) + cos 2 ( 3 α 2 ) ( v 0 2 + v 1 2 ) 2 sin 2 α v 0 v 5 2 cos 2 ( α 2 ) ( v 0 v 4 + v 1 v 3 + v 1 v 4 ) 2 cos 2 ( 3 α 2 ) v 1 v 2 .
For a harmonic interparticle potential, the Hamiltonian for straight and bent stripes configurations we expanded around mechanical equilibrium takes the quadratic form: H = K m , n A m , n q m q n , where { q } = { u , v } represents small displacements about the equilibrium position of every particle. In the canonical ensemble the difference between the entropy per particle of the straight and bent configurations for such Hamiltonian is [20]: Δ s = ( S s S b ) / N = 1 / ( 2 N ) ln ( A b / A s ) where the subscript s refers to straight and b to bent, and A is the determinant of A. The dimensionless matrix A depends only on the deformation angle α and on the zigzagging-stripe realization. In [20] we used a recursive method to obtain the matrix A in the case of the elastic Ising model for any subset of the network composed by shells of particles (see Figure 1). Here we apply the same method to the 2D spring network model. The number of particles n belonging to the shells up to n s is given by n = 1 + 3 n s ( n s 1 ) . In our shell-expansion calculation, these n particles are free to move, while the other N n particles of the network are frozen in their equilibrium position. Increasing n, Δ s should converge to the exact result (see Figure 2a), which includes the simultaneous fluctuation of all particles in the system.
In Figure 2a we show Δ s for the 2D harmonic network model for a number of shells up to n s = 20 , that is n = 1141 particles free to move. From it we can see that Δ s > 0 for every deformation angle α of the network and for every order n s of the expansion (except for very small deviations at small α and small n s ). Considering only one particle free to move, n s = 1 , gives a qualitative indication on the behavior of Δ s for every orders of approximation. This rapid convergence with increasing n s gives confidence in this expansion method when we will apply it for purely repulsive interactions, for which we are technically much more limited in the number of particles that we may numerically calculate the simultaneous fluctuation of.
Now we consider the same 2D triangular network model, but for repulsive harmonic interactions, that is we consider the following Hamiltonian
H r = m , n l = 1 3 K 2 d r l d r 0 2 θ ( d r 0 d r l ) ,
where θ ( ) is the Heaviside step function. Namely, now each spring applies a restoring force only when compressed ( d r l < d r 0 ) and there is no resistance to stretching ( d r l > d r 0 ). In this case we have to numerically integrate the partition function in order to get the entropy of straight- and bent-stripe configurations. In Figure 2b we show Δ s for repulsive harmonic interactions for n = 1 , 2 , 3 . Numerical calculations for n > 3 are beyond our computational reach (numerical integration for a number of variables larger than 6, in our case suffers from fluctuating results for a limited capacity in the precision). For n = 1 particle free to move, each particle can be equivalently chosen to be free. For n = 2 we consider the particles 0 and 1 (see Figure 1c,d). For n = 3 we consider the particles 0, 1, 2 and 0, 1, 4 for the straight and bent configurations, respectively (see Figure 1c,d). As for the 3D spring network of [20], we find that also in our 2D model Δ s < 0 for repulsive interactions.
For the case n = 1 we can show how the inversion of the stable phase when turning from harmonic to repulsive harmonic potential depends on the spatial configurations of straight and bent stripes, and in particular of the angular distribution of neighboring particles around each particle in the network.
For n = 1 the computation of the canonical partition function for the repulsive harmonic system can be easily reduced to the integration of single variable functions using polar coordinates ( ρ , γ ) (see Figure 3). The Hamiltonian of the free particle 0 can be written as H 0 r = 1 / 2 K ρ 2 i = 1 6 g i ( α , γ ) where ρ 2 g i ( α , γ ) is the contribution to H 0 r coming from the neighbor i of the particle 0, and the function g i depends on the coordinates of the particle i, as specified below. The canonical partition function is thus
Z 0 r = 0 0 2 π exp β K 2 ρ 2 i = 1 6 g i ( α , γ ) ρ d ρ d γ = 1 β K 0 2 π d γ i = 1 6 g i ( α , γ ) = I ( α ) β K
where β = 1 / ( K B T ) is the Boltzmann factor and I ( α ) = 0 2 π [ i = 1 6 g i ( α , γ ) ] 1 d γ = 0 2 π f 1 ( α , γ ) d γ with f ( α , γ ) = i = 1 6 g i ( α , γ ) . In this case we have Δ s = ln ( I s ( α ) / I b ( α ) ) .
The difference between the calculation of the harmonic and the repulsive harmonic partition function is that in the former case the functions g i contributes to the integral I ( α ) for any angle 0 γ 2 π , while in the latter case every g i contributes to I ( α ) for a specific range of angles only. For the repulsive harmonic potential in the straight-stripes configuration we need to consider the azimuthal ranges coming from each one of the neighboring particles:
f s r ( α , γ ) = cos 2 γ θ ( γ ) θ ( π / 2 γ ) + θ ( γ 3 π / 2 ) θ ( 2 π γ ) + sin 2 ( α / 2 + γ ) [ θ ( γ ) θ ( π α / 2 γ ) + θ ( γ 2 π + α / 2 ) θ ( 2 π γ ) ] + sin 2 ( α / 2 γ ) θ ( γ α / 2 ) θ ( π + α / 2 γ ) + cos 2 γ θ ( γ π / 2 ) θ ( 3 π / 2 γ ) + sin 2 ( α / 2 + γ ) θ ( γ π + α / 2 ) θ ( 2 π α / 2 γ ) + sin 2 ( α / 2 γ ) θ ( γ π α / 2 ) θ ( 2 π γ ) + θ ( γ ) θ ( α / 2 γ ) .
Due to the symmetry of the straight-stripes configuration, thanks to which the reflection about the origin of each neighbor transforms it in another neighboring particle (see Figure 1c and Table 1), we can write the function f s r by just taking the contribution of every g i without the condition imposed by the Heaviside step functions and dividing it by 2, that is
f s r ( α , γ ) = cos 2 γ + sin 2 ( α / 2 + γ ) + sin 2 ( α / 2 γ ) = cos 2 γ ( 2 cos α ) + cos 2 ( α / 2 )
For the repulsive harmonic potential in the bent configuration we need to consider the azimuthal ranges coming from each one of the neighboring particles:
f b r ( α , γ ) = cos 2 γ θ ( γ ) θ ( π / 2 γ ) + θ ( γ 3 π / 2 ) θ ( 2 π γ ) + sin 2 ( α / 2 + γ ) [ θ ( γ ) θ ( π α / 2 γ ) + θ ( γ 2 π + α / 2 ) θ ( 2 π γ ) ] + cos 2 ( α + γ ) θ ( γ π / 2 + α ) θ ( 3 π / 2 α γ ) + sin 2 ( 3 α / 2 + γ ) θ ( γ π + 3 α / 2 ) θ ( 2 π 3 α / 2 γ ) + sin 2 ( α / 2 + γ ) [ θ ( γ π + α / 2 ) · θ ( 2 π α / 2 γ ) ] + sin 2 ( α / 2 γ ) θ ( γ π α / 2 ) θ ( 2 π γ ) + θ ( γ ) θ ( α / 2 γ )
In Figure 4 we show 1 / f for hamonic and repulsive harmonic interactions for straight and bent configurations as a function of the lattice deformation angle α and the azimuthal direction in space γ , over which we numerically integrate in order to get the value of the function I for that specific angle α , which in turn sets the entropy via Equation (9). From Figure 4 we can see, particularly for big angles α for which Δ s takes its larger values (see Figure 2b), that repulsion accentuates the contribution to the function I for bent stripes (corresponding in Figure 4 to a wider region composed of brighter colors, i.e., white, yellow and red, for repulsive harmonic over harmonic interaction in the case of bent stripes). More in general, we can say that repulsion accentuates differences in the contribution to the partition function and thus to the free energy between symmetric and asymmetric distribution of neighboring particles.

3. Conclusions

We studied a 2D triangular-network model composed of particles interacting through harmonic or repulsive harmonic springs. At T = 0 the ground state is degenerate and composed of straight or any set of zigzagging-stripes configurations. At T > 0 we found that the stable phase is composed of straight or bent stripes depending on the harmonic or repulsive harmonic nature of particle interaction, respectively. This selection mechanism of the stable phase through the order-by-disorder effect is equivalent to that observed in the colloidal [20] and Ising [25] antiferromagnets irrespective of the dimensionality of the system. This suggests that the phase inversion of isosceles triangular networks is controlled by the attraction component of the interparticle interaction. We suggest that this is due to the fact that repulsive interactions accentuate differences in the contribution to the free energy between symmetric and asymmetric distribution of neighboring particles, as we have shown for n = 1 free particle calculation.
Both the Ising antiferromagnet on a deformable triangular lattice and the 2D isosceles triangular network model at T = 0 from one side, and the colloidal monolayer at ρ = ρ c p from the other side, have the same degeneracy and subextensive entropy S N and thus a vanishing residual entropy per particle in the thermodynamic limit. At T > 0 for the triangular networks and at ρ < ρ c p for the colloidal monolayer this degeneracy is removed, but they can still have a residual entropy per particle for finite system size if the ergodicity is broken. Indeed, even at T > 0 or ρ < ρ c p a system may be trapped in a local minimum of the free-energy landscape and thermal fluctuations are not large enough for a small system to overcome energy barriers. For ergodic systems the time average of observables can be computed by using ensemble averages thanks to the Birkhoff theorem [38]. From the other hand, for non-ergodic systems, the phase space is divided into disjoined sets. In this case, states can be counted either following the kinetic view [39], for which only states visited by the system at the observational time scale are taken into account, or following the Edwards approach [40], for which all possible states are considered regardless of whether they are explored or not by the system. Recently, the Edwards hypothesis has been proved to be valid at the un-jamming point [41]. In thermal ergodic systems at equilibrium, the two sampling methods give the same result. Following the Edwards approach we can conclude that an indication of the presence of residual entropy in a system is given by the ergodicity breaking (which can be checked for generic temperature or density) instead of by the degeneracy of the ground state (defined for T = 0 or ρ = ρ c p only).

Acknowledgments

This research was supported by the Israel Science Foundation Grant No. 968/16.

Author Contributions

Authors contributed equally to research and writing of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pauling, L. The structure and entropy of ice and of other crystals with some randomness of atomic arrangement. J. Am. Chem. Soc. 1935, 57, 2680–2684. [Google Scholar] [CrossRef]
  2. Harris, M.J.; Bramwell, S.T.; McMorrow, D.F.; Zeiske, T.; Godfrey, K.W. Geometrical frustration in the ferromagnetic pyrochlore HO2Ti2O7. Phys. Rev. Lett. 1997, 79. [Google Scholar] [CrossRef]
  3. Libál, A.; Nisoli, C.; Reichhardt, C.J.O.; Reichhardt, C. Inner Phases of Colloidal Hexagonal Spin Ice. Phys. Rev. Lett. 2018, 120, 027204. [Google Scholar] [CrossRef] [PubMed]
  4. Bramwell, S.T.; Gingras, M.J. Spin ice state in frustrated magnetic pyrochlore materials. Science 2001, 294, 1495–1501. [Google Scholar] [CrossRef] [PubMed]
  5. Wang, F.; Nisoli, C.; Freitas, R.S.; Li, J.; McConville, W.; Cooley, B.J.; Lund, M.S.; Samarth, N.; Leighton, C.; Crespi, V.H.; et al. Artificial ’spin ice’ in a geometrically frustrated lattice of nanoscale ferromagnetic islands. Nature 2006, 439, 303–306. [Google Scholar] [CrossRef] [PubMed]
  6. Anderson, P.W. The resonating valence bond state in La2CuO4 and superconductivity. Science 1987, 235, 1196–1198. [Google Scholar] [CrossRef] [PubMed]
  7. Kang, S.H.; Shan, S.; Košmrlj, A.; Noorduin, W.L.; Shian, S.; Weaver, J.C.; Clarke, D.R.; Bertoldi, K. Complex ordered patterns in mechanical instability induced geometrically frustrated triangular cellular structures. Phys. Rev. Lett. 2014, 112, 098701. [Google Scholar] [CrossRef] [PubMed]
  8. Coulais, C.; Teomy, E.; de Reus, K.; Shokef, Y.; van Hecke, M. Combinatorial design of textured mechanical metamaterials. Nature 2016, 535, 529–532. [Google Scholar] [CrossRef] [PubMed]
  9. Han, Y.; Shokef, Y.; Alsayed, A.M.; Yunker, P.; Lubensky, T.C.; Yodh, A.G. Geometric frustration in buckled colloidal monolayers. Nature 2008, 456, 898–903. [Google Scholar] [CrossRef] [PubMed]
  10. Libál, A.; Reichhardt, C.; Reichhardt, C.J.O. Realizing colloidal artificial ice on arrays of optical traps. Phys. Rev. Lett. 2006, 97, 228302. [Google Scholar] [CrossRef] [PubMed]
  11. Ortiz-Ambriz, A.; Tierno, P. Engineering of frustration in colloidal artificial ices realized on microfeatured grooved lattices. Nat. Commun. 2016, 7, 10575. [Google Scholar] [CrossRef] [PubMed]
  12. Tierno, P. Geometric frustration of colloidal dimers on a honeycomb magnetic lattice. Phys. Rev. Lett. 2016, 116. [Google Scholar] [CrossRef] [PubMed]
  13. Gokhale, S.; Sood, A.K.; Ganapathy, R. Deconstructing the glass transition through critical experiments on colloids. Adv. Phys. 2016, 65, 363–452. [Google Scholar] [CrossRef]
  14. Lu, P.J.; Weitz, D.A. Colloidal particles: Crystals, glasses, and gels. Annu. Rev. Condens. Matter Phys. 2013, 4, 217–233. [Google Scholar]
  15. Schmidt, M.; Löwen, H. Freezing between two and three dimensions. Phys. Rev. Lett. 1996, 76, 4552–4555. [Google Scholar] [CrossRef] [PubMed]
  16. Schmidt, M.; Löwen, H. Phase diagram of hard spheres confined between parallel hard plates. Phys. Rev. E 1997, 55, 7228–7241. [Google Scholar] [CrossRef]
  17. Oğuz, E.C.; Marechal, M.; Ramiro-Manzano, F.; Rodriguez, I.; Messina, R.; Meseguer, F.J.; Löwen, H. Packing confined hard spheres denser with adaptive prism phases. Phys. Rev. Lett. 2012, 109, 218301. [Google Scholar] [CrossRef] [PubMed]
  18. Shokef, Y.; Lubensky, T.C. Stripes, zigzags, and slow dynamics in buckled hard spheres. Phys. Rev. Lett. 2009, 102, 048303. [Google Scholar] [CrossRef] [PubMed]
  19. Zhou, D.; Wang, F.; Li, B.; Lou, X.; Han, Y. Glassy spin dynamics in geometrically frustrated buckled colloidal crystals. Phys. Rev. X 2017, 7, 021030. [Google Scholar] [CrossRef]
  20. Leoni, F.; Shokef, Y. Attraction controls the inversion of order by disorder in buckled colloidal monolayers. Phys. Rev. Lett. 2017, 118, 218002. [Google Scholar] [CrossRef] [PubMed]
  21. Wannier, G.H. Antiferromagnetism. The triangular Ising net. Phys. Rev. 1950, 79, 357–364. [Google Scholar] [CrossRef]
  22. Harth, K.; Mauney, A.; Stannarius, R. Frustrated packing of spheres in a flat container under symmetry-breaking bias. Phys. Rev. E 2015, 91. [Google Scholar] [CrossRef] [PubMed]
  23. Lévay, S.; Fischer, D.; Stannarius, R.; Szabó, B.; Börzsönyic, T.; Törökad, J. Frustrated packing in a granular system under geometrical confinement. Soft Matter 2018, 14, 396–404. [Google Scholar] [CrossRef] [PubMed]
  24. Ilyin, V.; Procaccia, I.; Regev, I.; Shokef, Y. Randomness-induced redistribution of vibrational frequencies in amorphous solids. Phys. Rev. B 2009, 80, 174201. [Google Scholar] [CrossRef]
  25. Shokef, Y.; Souslov, A.; Lubensky, T.C. Order by disorder in the antiferromagnetic Ising model on an elastic triangular lattice. Proc. Natl. Acad. Sci. USA 2011, 108, 11804–11809. [Google Scholar] [CrossRef] [PubMed]
  26. Dublenych, Y.I. Ground states of the Ising model on an anisotropic triangular lattice: Stripes and zigzags. J. Phys. Condens. Matter 2013, 25, 406003. [Google Scholar] [CrossRef] [PubMed]
  27. Lorman, V.; Podgornik, R.; Žekš, B. Positional, Reorientational, and Bond Orientational Order in DNA Mesophases. Phys. Rev. Lett. 2001, 87, 218101. [Google Scholar] [CrossRef] [PubMed]
  28. Camp, P.J.; Fuertes, A.; Attfield, J.P. Subextensive entropies and open order in Perovskite Oxynitrides. J. Am. Chem. Soc. 2012, 134, 6762–6766. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Masanes, L.; Oppenheim, J. A general derivation and quantification of the third law of thermodynamics. Nat. Commun. 2017, 8, 14538. [Google Scholar] [CrossRef] [PubMed]
  30. Villain, J.; Bidaux, R.; Carton, J.P.; Conte, R. Order as an effect of disorder. J. Phys. 1980, 41, 1263–1272. [Google Scholar] [CrossRef]
  31. Henley, C.L. Ordering by disorder: Ground-state selection in fcc vector antiferromagnets. J. Appl. Phys. 1987, 61, 3962–3964. [Google Scholar] [CrossRef]
  32. Henley, C.L. Ordering due to disorder in a frustrated vector antiferromagnet. Phys. Rev. Lett. 1989, 62, 2056–2059. [Google Scholar] [CrossRef] [PubMed]
  33. Chubukov, A. Order from disorder in a kagome antiferromagnet. Phys. Rev. Lett. 1992, 69, 832–835. [Google Scholar] [CrossRef] [PubMed]
  34. Reimers, J.N.; Berlinsky, A.J. Order by disorder in the classical Heisenberg kagomé antiferromagnet. Phys. Rev. B 1993, 48, 9539–9554. [Google Scholar] [CrossRef]
  35. Bergman, D.; Alicea, J.; Gull, E.; Trebst, S.; Balents, L. Order-by-disorder and spiral spin-liquid in frustrated diamond-lattice antiferromagnets. Nat. Phys. 2007, 3, 487–491. [Google Scholar] [CrossRef]
  36. Patrykiejew, A.; Staszewski, T. Ordering and order-disorder phase transition in the (1 × 1) monolayer chemisorbed on the (111) face of an fcc crystal. Condens. Matter Phys. 2016, 19, 13001. [Google Scholar] [CrossRef]
  37. Guruciaga, P.; Tarzia, M.; Ferreyra, M.; Grigera, L.C.S.; Borzi, R. Field-tuned order by disorder in frustrated Ising magnets with antiferromagnetic interactions. Phys. Rev. Lett. 2016, 117, 167203. [Google Scholar] [CrossRef] [PubMed]
  38. Birkhoff, G.D. Proof of the Ergodic Theorem. Proc. Natl. Acad. Sci. USA 1931, 17, 656–660. [Google Scholar] [CrossRef] [PubMed]
  39. Mauro, J.C.; Gupta, P.K.; Loucks, R.J. Continuously broken ergodicity. J. Chem. Phys. 2007, 126, 184511. [Google Scholar] [CrossRef] [PubMed]
  40. Edwards, S.F.; Oakeshott, R.B.S. Theory of powders. Physica A 1989, 157, 1080–1090. [Google Scholar] [CrossRef]
  41. Martiniani, S.; Schrenk, K.J.; Ramola, K.; Chakraborty, B.; Frenkel, D. Numerical test of the Edwards conjecture shows that all packings are equally probable at jamming. Nat. Phys. 2017, 13, 848–851. [Google Scholar] [CrossRef]
Figure 1. (a) Straight and (b) maximally zigzagging or bent configuration. Shells, which order is indicated with n s , are denoted by increasingly darker colors for increasing n s . Thicker, gray lines correspond to springs of longer (for α > π / 3 ) or shorter (for α < π / 3 ) length at rest; (c,d) are the unit cell for the straight- and bent-stripes confgirations, respectively. Numbers associated to particles correspond to the particle positions in tables in Section 2; (e) Parameters associated to every plaquette: a, b, c and α .
Figure 1. (a) Straight and (b) maximally zigzagging or bent configuration. Shells, which order is indicated with n s , are denoted by increasingly darker colors for increasing n s . Thicker, gray lines correspond to springs of longer (for α > π / 3 ) or shorter (for α < π / 3 ) length at rest; (c,d) are the unit cell for the straight- and bent-stripes confgirations, respectively. Numbers associated to particles correspond to the particle positions in tables in Section 2; (e) Parameters associated to every plaquette: a, b, c and α .
Entropy 20 00122 g001
Figure 2. (a) Entropy difference per particle Δ s , between straight- and bent-stripe configurations vs. deformation angle α for increasing numbers of shells n s of fluctuating particles; (b) Δ s vs. α for n = 1 , 2 , 3 particles free to move for harmonic (open symbols) and repulsive harmonic (filled symbols) interactions.
Figure 2. (a) Entropy difference per particle Δ s , between straight- and bent-stripe configurations vs. deformation angle α for increasing numbers of shells n s of fluctuating particles; (b) Δ s vs. α for n = 1 , 2 , 3 particles free to move for harmonic (open symbols) and repulsive harmonic (filled symbols) interactions.
Entropy 20 00122 g002
Figure 3. Example of straight (a) and bent (b) configuration for n = 1 particle free to move. The deviation of the central particle from its equilibrium position is described by polar coordinates ( ρ , γ ).
Figure 3. Example of straight (a) and bent (b) configuration for n = 1 particle free to move. The deviation of the central particle from its equilibrium position is described by polar coordinates ( ρ , γ ).
Entropy 20 00122 g003
Figure 4. Color diagram of 1 / f , as defined in the text, as a function of the angles α and γ , for harmonic and repulsive harmonic interactions and for straight- and bent-stripes configurations.
Figure 4. Color diagram of 1 / f , as defined in the text, as a function of the angles α and γ , for harmonic and repulsive harmonic interactions and for straight- and bent-stripes configurations.
Entropy 20 00122 g004
Table 1. Distances between the neighboring particles and the central particle 0 in the unit cell of the straight-stripe configuration. Particle positions are graphically shown in Figure 1c.
Table 1. Distances between the neighboring particles and the central particle 0 in the unit cell of the straight-stripe configuration. Particle positions are graphically shown in Figure 1c.
Particles dx dy dr 0
1, 0 2 a 0 2 a
2, 0abc
3, 0 a bc
4, 0 2 a 0 2 a
5, 0 a b c
6, 0a b c
Table 2. Distances between the neighboring particles and 0 and 1 particles in the unit cell of the maximally zigzagging-stripe configuration. Particle positions are graphically shown in Figure 1d.
Table 2. Distances between the neighboring particles and 0 and 1 particles in the unit cell of the maximally zigzagging-stripe configuration. Particle positions are graphically shown in Figure 1d.
Particles dx dy dr 0
1, 0 2 a 0 2 a
4, 0abc
5, 0 2 a 1 8 a 2 c 2 4 b 1 b 2 c 2 2 a
6, 0 a 1 4 b 2 c 2 b 3 4 b 2 c 2 c
7, 0 a b c
8, 0a b c
2, 1 a 1 4 b 2 c 2 b 3 4 b 2 c 2 c
3, 1abc
4, 1 a bc
8, 1 a b c
9, 1 2 a 1 2 b 2 c 2 4 b 1 b 2 c 2 2 a

Share and Cite

MDPI and ACS Style

Leoni, F.; Shokef, Y. Attraction Controls the Entropy of Fluctuations in Isosceles Triangular Networks. Entropy 2018, 20, 122. https://doi.org/10.3390/e20020122

AMA Style

Leoni F, Shokef Y. Attraction Controls the Entropy of Fluctuations in Isosceles Triangular Networks. Entropy. 2018; 20(2):122. https://doi.org/10.3390/e20020122

Chicago/Turabian Style

Leoni, Fabio, and Yair Shokef. 2018. "Attraction Controls the Entropy of Fluctuations in Isosceles Triangular Networks" Entropy 20, no. 2: 122. https://doi.org/10.3390/e20020122

APA Style

Leoni, F., & Shokef, Y. (2018). Attraction Controls the Entropy of Fluctuations in Isosceles Triangular Networks. Entropy, 20(2), 122. https://doi.org/10.3390/e20020122

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop