Critical Casimir Force between Inhomogeneous Boundaries

Critical Casimir Force between Inhomogeneous Boundaries


To study the critical Casimir force between chemically structured boundaries immersed in a binary mixture at its demixing transition, we consider a strip of Ising spins subject to alternating fixed spin boundary conditions. The system exhibits a boundary induced phase transition as function of the relative amount of up and down boundary spins. This transition is associated with a sign change of the asymptotic force and a diverging correlation length that sets the scale for the crossover between different universal force amplitudes. Using conformal field theory and a mapping to Majorana fermions, we obtain the universal scaling function of this crossover, and the force at short distances.

11.25.Hf, 05.40.-a, 68.35.Rh

Fluctuation-induced forces are generic to all situations where fluctuations of a medium or field are confined by boundaries. Examples include QED Casimir forces Casimir (1948); Bordag et al. (2009), van der Walls forces Parsegian (2005), and thermal Casimir forces in soft matter which are most pronounced near a critical point where correlation lengths are large de Gennes and Fisher (1978); Krech (1994). The interaction is then referred to as critical Casimir force (CCF). Analogies and differences between these variants of the common underlying effect have been reviewed in Ref. Gambassi (2009).

Experimentally, CCFs can be observed indirectly in wetting films of critical fluids Nightingale and Indekeu (1985), as has been demonstrated close to the superfluid transition of He Garcia and Chan (1999) and binary liquid mixtures Garcia and Chan (2002). More recently, the CCF between colloidal particles and a planar substrate has been measured directly in a critical binary liquid mixture Hertlein et al. (2008); Soyka et al. (2008). Motivated by the possibility that the lipid mixtures composing biological membranes are poised at criticality Veatch et al. (2007); Baumgart et al. (2007), it has been also proposed that inhomogeneities on such membranes are subject to a CCF Machta et al. (2012) which provides an example of a 2D realisation.

The amplitude of the CCF is in general a universal scaling function that is determined by the universality classes of the fluctuating medium Diehl (1986). It depends on macroscopic properties such as the surface distance, shape and boundary conditions of the surfaces but is independent of microscopic details of the system Krech (1994). Controlling the sign of fluctuation forces (attractive or repulsive) is important to a myriad of applications in design and manipulation of micron scale devices. While for QED Casimir forces a generalized Earnshaw’s theorem rules out the possibility of stable levitation (and consequently force reversals) in most cases Rahi et al. (2010), the sign of the CCF depends on the boundary conditions at the confinement. For classical binary mixtures, surfaces have a preference for one of the two components, corresponding to fixed spin boundary conditions ( or ) in the corresponding Ising universality class. Depending on whether the conditions are like ( or ) or unlike ( or ) on two surfaces, the CCF between them is attractive or repulsive. So-called ordinary or free spin boundary conditions are difficult to realize experimentally but can emerge due to renormalization of inhomogeneous conditions as we shall show below Toldin et al. (2013).

Motivated by their potential relevance to nano-scale devices, fluctuation forces in the presence of geometrically or chemically structured surfaces have been at the focus recently. Sign changes of CCFs due to wedge like surface structures have been reported very recently Bimonte et al. (2015). Competing boundary conditions can give rise to interesting crossover effects with respect to strength and even sign of the forces. Here we consider such a situation for the Ising universality class in 2D. At criticality, this system can be described by conformal field theory (CFT) Friedan et al. (1984); Cardy (1989a), and CCFs are related to the central charge of the CFT Cardy (1986); Kleban and Vassileva (1991); Kleban and Peschel (1996), and scaling dimensions of boundary operators Cardy (1989b).

Figure 1: Ising strip of width with alternating fixed spin boundary conditions on one side, with a typical spin configuation indicated by the shading.

In this Letter, we show that boundary conditions which alternate periodically between two spin states (see Fig. 1) give rise to a novel phase transition. Associated with that is a diverging correlation length that sets the scale for a sign change of the CCF on one side of the transition. We obtain the critical exponents and exact expressions for the universal scaling function of the force in the critical region. Consider the Ising model on an infinitely long strip of width , and assume that the system is at its critical temperature so that it is conformally invariant. For homogenous, fixed spin boundary conditions , on the two boundaries, the critical Casimir energy per unit strip length, , is determined by CFT. Since is the only finite length scale, the energy obeys a simple power law. The amplitude is determined only by the central charge of the Ising model and the scaling dimension of the so-called boundary condition changing (BCC) operator from to (see below for details on the BCC operator) Cardy (1989b),


where we measure here and the following energies in units of . For like boundary conditions or one has and hence an attractive force . For unlike boundary conditions one gets and hence a repulsive force. There is one more conformally invariant boundary condition that corresponds to free (f) spins or ordinary boundary conditions. When combined with fixed boundary conditions, the corresponding BCC operator has the scaling dimension which implies a repulsive interaction in Eq. (1). In the following we consider a strip with homogeneous spins on one boundary and alternating regions of and spins of length and , respectively, on the other boundary, see Fig. 1.

If the temperature is slightly different from , the system is in the critical region, where the free energy density can be decomposed into non-singular () and singular () contributions,


that depend on the reduced temperature , the width , and a scaling variable that is specific to the alternating boundary conditions in Fig. 1. While the non-singular part is an analytic function of and , the singular part is not. For homogeneous boundary conditions, is the only relevant scaling variable, and in the critical region the singular part of the free energy density is given by a universal scaling function that depends only on Diehl (1986); Krech (1994) where is the bulk correlation length with amplitude and exponent for the Ising model. As we shall see below, the same renormalization-group (RG) concepts apply to a novel, boundary induced critical region that we identify for inhomogeneous boundary conditions around . To focus on that region, we assume in the following that the system is at its bulk critical point, . For large the singular part of the free energy density can be expressed in terms of a universal scaling function of the new correlation length ,


Below we shall determine and the exponent .

BCC operators have been introduced in CFT to study systems with discontinuous boundary conditions Cardy (1989b). When inserted on a boundary, these local operators interpolate between the different boundary conditions on either side of the insertion point. They are highest weight states of weight and all such states may be realized by an appropriate pair of boundary conditions. For the critical Ising model, the BCC operator that takes the boundary condition from spin to spin corresponds to the chiral part of the energy operator . This can be understood easily in the representation of the Ising model in terms of a free Majorana fermion field out of which the energy operator is composed, di Francesco et al. (1997): The Jordan-Wigner transformation shows that the fermion creation and annihilation operators flip locally the spin orientation.

Now the BCC operators permit us to relate the partition function of the strip with alternating boundary conditions to a correlator for the field at positions where the boundary conditions change. On the upper complex plane, one has which yields (after a conformal map) for the partition function of the strip the Pfaffian,


with , where we used the Wick theorem for fermions, are the positions of the BCC operators on the upper edge of the strip, and is the partition function of the homogenous system with . Due to the symmetry under translations by , the matrix is of block Toeplitz form, , with


where .

The free energy density can be expressed in the thermodynamic limit as


The Szegö-Widom (SW) theorem for block Toeplitz matrices states that the determinant can be expressed in terms of the matrix valued Fourier series as Widom (1974)


where acts now on a matrix. It turns out that this formula can be only applied for the case . The reason for that is a subtle difference between the Toeplitz matrix and the corresponding circulant matrix that describes periodic boundary conditions along the strip. While for the spectra of and become equivalent for , for there exists a pair of eigenvalues of that tend to zero exponentially for , yielding an extra contribution that is determined by the decay of the Fourier integral


and has to be subtracted from the r.h.s. of Eq. (7) for . Here denotes the -element of the matrix . In the following we apply Eqs. (7) and (8) to compute the critical Casimir force in various scaling limits.

When , the function defined below Eq. (5) can be replaced by which yields the exact determinant


For , the SW theorem then yields . For , this is also the correct result as it follows from subtracting the correction which follows from Eq. (8) and as . It follows that the critical Casimir force for is


It has an analytic amplitude that varies continously with . This result is identical to an addition of the amplitudes from Eq. (1) for unlike and like boundary conditions, weighted by and , according to their occurrence. Hence, additivity holds at short distances. This has been observed also for a 3D Ising model in the special case of boundaries with alternating stripes of equal width Toldin et al. (2013).

Next, we consider the case . Using the Abel-Plana summation formula, it can be shown that in this limit the elements of the matrix




with and


where is a hypergeometric function. For , the SW theorem then yields the free energy density


where we have subtracted a -independent contribution that does not change the force, and defined . In the evaluation of the integral, the correlation length defined above Eq. (3) becomes important. The integrand is exponentially localized around over a small range . Also, it can be shown that has the scaling property for any constant . Hence, in the critical region of small , or , the proper scaling is obtained by setting (up to a numerical coefficient), showing that the exponent . In the integral, can be replaced by and one obtains after a simple integration the result for the universal scaling function of Eq. (3) when or ,


where is a polylogarithm function. Outside the critical region , one has so that the force is fully dominated by the boundary regions with like spins. On the contrary, for , and hence , the frustration between almost equal amounts of fixed and spins on the boundaries leads to a renormalization to effectively free boundary conditions with . For , the correction yields an extra contribution determined by


so that the scaling function for is . Since , the scaling function is continuous around . For , however, so that the system asymptotically realizes homogenous unlike boundary conditions with .

Figure 2: Schematic overview of critical Casimir force amplitudes as function of the strip width and the ratio . For the solid curves represent the diverging correlation length . The horizontal dashed line indicates the cut along which the force amplitude is plotted in Fig. 3. Along the red curve the sign of the force changes whereas the blue curve indicates only a change between two universal (repulsive) limits.

Our findings can be summarized by the scheme of Fig. 2. It shows the different scaling regimes and the corresponding asymptotic amplitudes of the Casimir force. At short distance the amplitude varies continuously across the critical point at , with a sign change at . For there exist three distinct regions: around appears a region where where the force is repulsive and approaches for asymptotic the universal amplitude for fixed-free spin boundary conditions. For , the force changes sign from attractive to repulsive when approaches , corresponding to a stable point. For , the force is always repulsive but the amplitude crosses over from to under an increase of beyond .

The dependence of the force on at fixed (see dashed horizontal line in Fig. 2) is determined by with a universal scaling function of the scaling variable that is defined on both sides of the critical point by . This function is shown in Fig. 3 where we used the results for of Eqs. (15), (16). In the critical region , one has the expansions


whereas for outside the critical region, ,


We see that is not analytic around and hence constitutes the singular part of the free energy density, see Eq. (2). This resembles the singular nature of scaling functions describing the bulk transition at .

Figure 3: Universal scaling function for the critical Casimir force as function of the scaling variable .

Our results show the existence of a novel phase transition for the critical Casimir force in the 2D Ising model that is induced by inhomogeneous boundary conditions with a varying ratio of up and down spins. We obtained exact expressions for the universal scaling function of the force. Due to the observed renormalization of boundary conditions, in binary mixtures, ordinary (free spin) boundary conditions can be realized experimentally and “switched” on and off by varying the distance , or an inhomogeneous surface field. The crossover between different universal amplitudes leads to a stable equilibrium point at for . The emergence of the novel phase transition at is related to the relevance of a surface magnetic field at a surface with free spin boundary conditions. This can be seen from the decay of the spin correlations along a single surface, with for free boundary conditions Cardy (1984); Diehl (1986). Since the surface field contributes an energy , the scaling dimension which is identical to our findings above. It is interesting to explore these concepts in general spatial dimensions for Ising and XY models, and tri-critical points which have an even richer spectrum of possible boundary conditions.

We thank M. Kardar for many fruitful discussions.


  1. H. B. G. Casimir, Proc. K. Ned. Akad. Wet. 51, 793 (1948).
  2. M. Bordag, G. L. Klimchitskaya, U. Mohideen, and V. M. Mostepanenko, Advances in the Casimir Effect (Oxford University Press, 2009).
  3. V. A. Parsegian, Van der Waals Forces (Cambridge University Press, 2005).
  4. P.-G. de Gennes and M. E. Fisher, C. R. Acad. Sci. Ser. B 287, 207 (1978).
  5. M. Krech, The Casimir effect in Critical systems (World Scientific, 1994).
  6. A. Gambassi, J. Phys.: Conf. Ser. 161, 012037 (2009).
  7. M. P. Nightingale and J. O. Indekeu, Phys. Rev. B 32, 3364 (1985).
  8. R. Garcia and M. H. W. Chan, Phys. Rev. Lett. 83, 1187 (1999).
  9. R. Garcia and M. H. W. Chan, Phys. Rev. Lett. 88, 086101 (2002).
  10. C. Hertlein, L. Helden, A. Gambassi, S. Dietrich, and C. Bechinger, Nature 451, 172 (2008).
  11. F. Soyka, O. Zvyagolskaya, C. Hertlein, L. Helden, and C. Bechinger, Phys. Rev. Lett. 101, 208301 (2008).
  12. S. L. Veatch, O. Soubias, S. L. Keller, and K. Gawrisch, PNAS 104, 17650 (2007).
  13. T. Baumgart, A. T. Hammond, P. Sengupta, S. T. Hess, D. A. Holowka, B. A. Baird, and W. W. Webb, PNAS 104, 3165 (2007).
  14. B. B. Machta, S. L. Veatch, and J. Sethna, Phys. Rev. Lett. 109, 138101 (2012).
  15. H.-W. Diehl, Phase Transitions and Critical Phenomena (Academic Press, London, 1986), vol. 10, p. 75.
  16. S. J. Rahi, M. Kardar, and T. Emig, Phys. Rev. Lett. 105, 070404 (2010).
  17. F. P. Toldin, M. Tröndle, and S. Dietrich, Phys. Rev. E 88, 052110 (2013).
  18. G. Bimonte, T. Emig, and M. Kardar, Phys. Lett. B 743, 138 (2015).
  19. D. Friedan, Z. Qui, and S. Shenker, Phys. Rev. Lett. 52, 1575 (1984).
  20. J. L. Cardy, in Fields, Strings, and Critical Phenomena, edited by E. Brézin and J. Zinn-Justin (Elsevier, New York, 1989a).
  21. J. Cardy, Nucl. Phys. B 275, 200 (1986).
  22. P. Kleban and I. Vassileva, J. Phys. A: Math. Gen. 24, 3407 (1991).
  23. P. Kleban and I. Peschel, Z. Phys. B 101, 447 (1996).
  24. J. L. Cardy, Nucl. Phys. B 324, 581 (1989b).
  25. P. di Francesco, P. Mathieu, and D. Sénéchal, Conformal Field Theory (Springer, 1997).
  26. H. Widom, Adv. Math. 13, 284 (1974).
  27. J. L. Cardy, Nucl. Phys. B 240, 514 (1984).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description