New Research In
Physical Sciences
Social Sciences
Featured Portals
Articles by Topic
Biological Sciences
Featured Portals
Articles by Topic
 Agricultural Sciences
 Anthropology
 Applied Biological Sciences
 Biochemistry
 Biophysics and Computational Biology
 Cell Biology
 Developmental Biology
 Ecology
 Environmental Sciences
 Evolution
 Genetics
 Immunology and Inflammation
 Medical Sciences
 Microbiology
 Neuroscience
 Pharmacology
 Physiology
 Plant Biology
 Population Biology
 Psychological and Cognitive Sciences
 Sustainability Science
 Systems Biology
Symmetry breaking in smectics and surface models of their singularities

Edited by Noel A. Clark, University of Colorado, Boulder, CO, and approved July 20, 2009 (received for review May 13, 2009)
Abstract
The homotopy theory of topological defects in ordered media fails to completely characterize systems with broken translational symmetry. We argue that the problem can be understood in terms of the lack of rotational Goldstone modes in such systems and provide an alternate approach that correctly accounts for the interaction between translations and rotations. Dislocations are associated, as usual, with branch points in a phase field, whereas disclinations arise as critical points and singularities in the phase field. We introduce a threedimensional model for twodimensional smectics that clarifies the topology of disclinations and geometrically captures known results without the need to add compatibility conditions. Our work suggests natural generalizations of the twodimensional smectic theory to higher dimensions and to crystals.
Many of the features and properties of ordered media, such as crystalline solids, magnets, and liquid crystals, are controlled by their defects: those points or lines of discontinuity at which the order changes abruptly or is illdefined. For example, defects have played prominent roles in understanding the plastic deformation of solids (1), the Kosterlitz–Thouless transition of the twodimensional XY model (2), the dislocationunbinding melting of solids (3), and the formation of colloidal crystals in nematic liquid crystals (4). One of the unifying features of the defects in all of these systems is that they are topological, in the sense that many of their properties depend only on the symmetry of the ordered medium. Based on this concept a general framework has been developed for the identification and classification of the topological defects in any broken symmetry medium in terms of group elements (or more generally conjugacy classes) of homotopy groups of a space of degeneracy, representing the ground state manifold of the system (5–11). This framework has provided an understanding of the ways in which defects can combine or split, and it has highlighted the path dependence of such processes when the fundamental group of the ground state manifold is nonAbelian (12).
It is known, however, that the standard topological theory of defects is incomplete in systems with spontaneously broken translation and rotation invariance, such as crystals (9, 11, 13). More specifically, a description of such systems that treats disclinations and dislocations (translational defects) and their interactions in a mathematically coherent fashion is missing. A topological classification of treating dislocations alone is available via the use of Burgers vectors and the like (14, 15). The standard way (16, 17) to add disclinations in these systems is to adjoin extra components representing the orientation of the system to the order parameter. For instance, the order parameter space of the twodimensional smectic becomes the Klein bottle. But in this formalism, the neat correspondence between the fundamental group of the order parameter space and disclinationtype defects breaks down, as we will demonstrate. For instance, it is impossible to have a disclination with index greater than +1 in a twodimensional smectic system (9, 18), yet there is no indication of this in the standard formalism. This breakdown is related to the number of Goldstone modes in these systems and the restrictions they impose on paths in the ground state manifold.
The goal of this article is thus to reconsider the topological theory of defects for systems with broken translational symmetry and argue that “less is more,” at least for components of the order parameter. We will show that it is mathematically consistent to keep the order parameter as a phase field of the sort whose topological defects are wellknown to generate the Burgers vector for dislocations and to look at disclinations in this picture as defects in the derivatives of the phase field. Modeling disclinations as critical points was considered by Trebin (11); we complete this approach and present a method that admits dislocations as well by noting that phase fields in layered systems are not maps to ℝ but rather to periodic, possibly nonorientable, spaces. We also suggest connections to singularity theory and Morse theories that can shed further light on our formulation.
We focus on liquid crystalline systems because they have lowdimensional order parameter spaces and use as our prototypical examples the twodimensional smectic and directed line systems. Smectic liquid crystals consist of rodshaped molecules that spontaneously form both directional (nematic) order and a onedimensional density wave, commonly described as a layered system; the spacing between the layers is approximately the rod length, a. In two dimensions, we can picture the layers as a set of nearly equally spaced curves lying in the plane. The ground states are characterized by both equal spacing between these curves and vanishing curvature. The free energy, accordingly, has two terms; the first is a compression energy that sets the spacing, and the second is a bending energy. Because they are periodic, smectics can be expressed as a density wave ρ(x) with wavelength 2π/q_{0}. Nearly uniform plane waves with wavevector q_{0} = q_{0}q̂_{0} can be written as ρ(x) = ρ_{0} + ρ_{1}Re e^{i(q0·x−q0u(x))}, where u is the Eulerian displacement function and φ = q̂_{0}·x − u is the phase field. Drawing the layers at places of maximum density, we see that the layers are those level sets where φ is an integer multiple of a. The energy may be written in terms of the phase field or the displacement; hence, we will look at the geometry of the smectic in terms of these functions. We note that smectics have an extra symmetry that the φ(x) description does not have, namely that φ → −φ is a symmetry that leads to their lack of orientability, which we will discuss below; for now it suffices to say that a configuration of layers arising from φ(x) carries an orientation. Hence, we call this the directed line system; indeed, it is topologically equivalent to wave systems [see, for instance, Nye and Berry's work on wave dislocations (19)]. Familiar examples of true twodimensional smectic configurations may be found on your fingerprints. Many other realizations of these systems have been studied experimentally (20, 21); even electronic systems can realize quantum smectic states (22).
Additionally, we introduce a geometrically intuitive model for looking at the topology of smectic defects that allows us to demonstrate our results. Because the phase field is the basic ingredient to study a smectic, it is natural to view the smectic as a surface placed over the plane. The singularities of this surface will then correspond to different defects. Our model motivates a highly geometric expression for the nonlinear smectic free energy functional and allows us to consider compression and curvature energies around defects.
We briefly review the standard topological theory of defects via examples, namely disclinations in nematics and dislocations in smectics. We recapitulate the limitations of the standard formalism when extended to disclinations in the directed line system and the smectic. Next, we present our resolution of these problems in terms of our construct and discuss a geometric formulation of the free energy for twodimensional smectics and use this to find compressionfree defects.
Finally, we conclude with applications and a discussion of avenues for extension.
Defects, Fundamental Groups, and Goldstone Modes
Systems with broken symmetries are described by nonvanishing order parameters. When the symmetry is continuous there is a degeneracy of ground states continuously connected by the underlying symmetry. Each of these equivalent but different ground states is represented by a different value of the order parameter. The order parameter field is a local measure of the system's order; in other words, each configuration of a medium defines a map from physical space (the plane, in this article) to the set of order parameters. The goal of the topological theory of defects is to classify the defect structures in physical media by analyzing the properties of these maps that are invariant under continuous deformation, or homotopy.
As an example, recall the familiar XY model of twodimensional spins on the plane. In the broken symmetry state the spins align along a common direction, making an angle θ with the +x̂ direction. The ground state corresponds to a constant value of θ, and the ground state manifold is thus the unit circle, S^{1}. Locally, the orientation at a point P in any texture corresponds to some point in the ground state manifold, and we can think of the ordered state as a map from points in the sample to directions given by points in S^{1}. Upon traversing a closed loop in the sample the orientation must return to its initial value so that θ can only change by an integer multiple of 2π, known as the winding number. A nonzero winding indicates that the loop encircles a topological defect, or vortex, whose strength is given by the winding number. Thus, defect states of the twodimensional XY model are characterized by a single integer.
These notions are formalized in the homotopy theory of topological defects (9–11). There we start with two groups, G and H, which are the symmetry groups of the system in a disordered and ordered state, respectively, so H ⊂ G. The ground state manifold is the quotient space G/H, and the topological defects are classified by the conjugacy classes of the fundamental group, π_{1}(G/H). In the XY model the disordered state has full rotational symmetry, so G = SO(2), which the ordered state breaks completely, so H is the trivial group. The order parameter manifold is G/H ≅ S^{1}, the unit circle, and the fundamental group is π_{1}(S^{1}) = ℤ, reproducing our intuitive result. Twodimensional nematics are almost identical, except that, because nematics have no heads or tails, we must identify θ and θ + π, so that now H = ℤ_{2}. We then have G/H = ℝP^{1} and π_{1}(ℝP^{1}) = (1/2)ℤ, giving the familiar Frank indexing of nematic disclinations (23). In this way, the topological theory provides a general framework for the classification of defects in any broken symmetry ordered medium. Among the key insights homotopy theory affords is that products of loops in the fundamental group yield information on how defects combine or split and, particularly, the path dependence of these processes when the fundamental group is nonAbelian (12).
The same approach may be applied, through a “naïve generalization” (9), to systems with broken translational symmetry, such as directed lines in the plane and twodimensional smectics. Beginning with the directed line model, the ground state consists of parallel layers, perpendicular to x̂ and with equal spacing a and a welldefined orientation as in Fig. 1A. The symmetries of the ground state consist of a translation of the layers by a in the x̂ direction, arbitrary translations in the ŷ direction, and rotations by 2π of the system. Thus, G is the twodimensional Euclidean group consisting of two translations and a rotation, H = ℝ × ℤ and G/H = S^{1} × S^{1}, the twotorus, where the first factor is translations modulo a and the second is rotations modulo 2π. Coordinates for the ground state manifold are naturally provided by the phase φ (acted upon by translations) and the direction θ of the layer normal (acted upon by rotations). If we consider homotopy classes of loops with nonzero winding only in the φ direction, we recover dislocations. In particular, the usual construction for the Burgers vector via counting layers is equivalent to looking at the winding number of φ, because we draw a new layer for every multiple of a in φ. The fundamental group of the torus is ℤ × ℤ, so in this generalization, defects in the full system should be characterized by pairs of integers (m,n), m corresponding to dislocation charge and n corresponding to disclination charge.
In the case of the smectic, the ground state is again a set of parallel lines perpendicular to x̂, but without orientation. Therefore, the symmetries include translation by a in the x̂ direction, and rotations by π owing to the nematiclike symmetry. A rotation by π followed by a translation by δ is equivalent to translation by −δ, which leads to a twist in the order parameter space. Again with φ,θ being the coordinates of the ground state manifold, we must now identify φ with −φ when θ → θ + π. The order parameter space is no longer the twotorus but is the Klein bottle. To construct the fundamental group of the Klein bottle we start with the free group of two elements S, representing a shift by a, and F, a rotation or flip by π. Note, however, that FS^{−1}F^{−1} = S, and it follows that a presentation of the fundamental group is 〈S,FFSF^{−1}S = e〉, where e is the group identity. This group is the semidirect product ℤ × ℤ_{2} (17), so one might classify loops by two numbers (m,n), where m counts the number of dislocations and n counts the number of disclinations. Algebraically, one can show that if n is odd, (m,n) is conjugate to (m+2,n), i.e., S^{−1}(FS^{m})S = FS^{m+2}. This can be interpreted graphically and elegantly as the generation of a dislocation pair in the presence of a +1/2 disclination, as in figure 3 of ref. 24.
Although the algebraic structure we have presented is the natural generalization of the standard study of topological defects (16, 25), it fails to correctly characterize the defects of either system. Mermin recognized this in ref. 9 and pointed out the difficulty of using the standard approach when translational symmetry is broken. Part of the issue is that the homotopy theory predicts entire classes of defects that are not present in the physical system. Poénaru (18) proved that directed line and smectic models cannot have disclinations with index greater than +1, and yet they are predicted by the homotopy theory. More generally, the problem lies in the allowed changes of the order parameter as one moves along a path in the sample. For a system with liquidlike order, such as the nematic, any variation in the order parameter can be accommodated, i.e., any path in the ground state manifold can be mapped onto the specified path in the sample. However, the same is not true for translationally ordered media, like the smectic. A continuous symmetry implies that two degenerate ground states can be connected via an arbitrarily lowenergy Goldstone excitation. This implies that the dimensionality of the allowed excitations corresponds to the number of Goldstone modes. Although the smectic breaks both rotational and translational symmetry, it has only one Goldstone mode, the Eulerian displacement field u(x), a signal of something different. Not any path in the ground state manifold can be realized, but only those that correspond to this mode. The reduction in the number of Goldstone modes and consequent restriction of the realizable paths in systems of broken translational symmetry arises because the rotational and translational degrees of freedom are coupled into a composite object (3, 13, 26) in a manner akin to the Higgs mechanism for gauge fields.*
Only when the number of Goldstone modes is equal to the dimensionality of the ground state manifold do we have sufficient freedom to be able to match any path in the ground state manifold to any path in the sample, as is required for the application of the standard homotopy theory. Our main result is to present a method by which the topological defects of translationally ordered media can be faithfully captured, without the need to remove entire classes of nonrealizable defects. The coupling of rotational and translational broken symmetries tells us that we should no longer consider the coordinates φ,θ of a point in the ground state manifold as freely variable. Recall that they are both given by the order parameter φ(x), the former as the value of the phase field at a point and the latter as the direction of its gradient. To identify the position in the ground state manifold of a configuration at any point P, therefore, we must consider the Taylor series for the configuration^{†} φ(P + x) ≈ φ(P) + ∇φ_{x=P}·x + ⋯. The ground states of the directed line and smectic systems are linear functions φ_{φ0,k}(x) = φ_{0} + k·x, where k is a unit vector. We assign the φ component of the local ground state to be φ_{0}, as usual, and use the directionality of ∇φ_{x=P} to set the second component, k. Note how this generalizes the usual procedure. If the ground state has structure up to order n, strictly speaking, we should keep data at every point of the configuration also up to order n.^{‡} There is an immediate implication regarding the smoothness of the phase field as well. When φ has a nontrivial winding about a point, we know it must be discontinuous there. Suppose, however, that a loop has nontrivial winding in ∇φ/∇φ. In this case, the unit normal must be discontinuous, but this only implies that ∇φ is discontinuous or zero. It follows that at a disclination φ is either singular or critical.
Of course, because any configuration can still be identified locally with a point in the ground state manifold, the defects that do occur will still correspond to nontrivial loops in G/H, and we can continue to label them by the winding numbers of these loops. However, the possible defects that can arise should be determined from the homotopy of φ(x) and ∇φ(x) together, and not from the structure of the ground state manifold. With this construction in hand, we move on to our surface model.
Surface Model
In this section we provide the first steps toward a theory where the objects being homotoped are not simply paths in the ground state manifold, but rather are smooth maps from the physical space to the space of the translational phase field. Thus, the theory of defects for these systems should be given in the language of singularity theory (27) and flavors of Morse theory (28), which may combine those on circlevalued maps (ref. 29, in particular see page 37) and on singular spaces (30). Rather than attempt to lay out precise mathematical conditions and proceed into the realm of rarefied mathematics, we discuss our ideas in an intuitive way, exploiting the following fact about smectics in two dimensions. Because the smectic is represented by level sets of a phase field on ℝ^{2}, we may visualize the phase field as a graph of a twosurface. From this perspective it becomes possible to sidestep the difficulties that we have discussed and to understand, for instance, Poénaru's rigorous result afresh. We begin with the simpler case of directed lines.
Starting with a ground state for the lines, φ = φ_{0} + k·x is the equation for a plane and defines a line for fixed values of φ. Taking level cuts at φ = na, n ∈ ℤ we find lines in the xyplane, with uniform spacing, a. Note, however, that the φ → φ + a symmetry implied by the phasefield construction of the layers means that instead of one plane, we have an infinite stack of equivalent planes each one shifted along the φ direction by a multiple of a as shown in Fig. 1A. Because of this, we can instead take a cut through all of the planes at φ = 0, for instance, and equivalently generate the ground state, viewing φ as an element of S^{1}. As we shall see, this multiplane description is natural from the point of view of defects; they connect the different sheets together to form the required topology. For instance, when taking a Burgers circuit around a +2 dislocation, we know that φ changes by 2a, which means that the sheets must be connected together with the topology of a helicoid so that dislocations may be thought of as branch points, as shown in Fig. 1B. Multiple dislocations are given by surfaces connected by multiple helicoids, e.g., Fig. 1C. This approach has been usefully used to visualize dislocations in waves (31),^{§} and we extend it to incorporate disclinations. Disclinations, as we have discussed, are either zeroes or singularities of ∇φ, and so in this height model, the disclinations correspond to critical points or cusps in the surface (11). For instance, a +1 disclination is represented by a “mountain” or “basin” in the height function, and the disclination sits at the zenith or nadir. Concretely, an equally spaced, +1 disclination is represented by a regular cone φ = x and the level sets at φ = na, giving circular layers with equal spacing, a. Similarly, negative charge disclinations correspond to saddles for creatures with various numbers of appendages.
The notion of topological equivalence in this model is intuitive. First, note that deformations of these surfaces (with the restriction that horizontally or vertically tangent surfaces are singular) correspond to the class of allowed paths in the φ,∇φ space that we discussed in the previous section.^{¶} Hence, results from Morse theory and singularity theory on the topology of these surfaces and their defects up to topological equivalence classes may be regarded as global consequences of the local restrictions developed earlier. As an example, consider Poénaru's result on the lack of disclinations of charge higher than +1. In our approach, Poénaru's result becomes transparent: two +1 disclinations cannot join because two mountains cannot come together on a surface without a mountain pass between them. A mountain pass has a critical point with the geometry of a saddle, and the contours of constant height draw out a −1 disclination in the twodimensional line system. Similar arguments hold for two basins and even for a basin and a mountain.
Smectics can also have 1/2 disclinations because they are nonorientable. Thus, in addition to the infinite set of planes normal to (−k,1) in the line system, the smectic ground state also contains the infinite set of planes normal to (k,1), as shown in Fig. 2A, because they generate the same set of unoriented layers. These two sets of planes intersect whenever φ is a halfinteger multiple of a, which reflects the fixed points of the two symmetry operations φ → φ + a and φ → −φ. At these heights it is possible to cross continuously from one set of planes to the other, thereby changing the local orientation of the layers. Defects again connect the different planes together, although now this can be done in two different ways; either by connecting planes with the same orientation, as in the directed line system, or planes with the opposite orientation. The latter leads to a global unorientability of the smectic and occurs when there are odd halfinteger index disclinations. The two prototypical halfinteger disclinations are depicted in Fig. 2B (+1/2) and C (−1/2). The +1/2 disclination consists of a pair of oppositely oriented planes terminated by a halfcone connecting them, whereas the −1/2 disclination is a threeway junction of pairs of oppositely oriented planes smoothly joined together. These examples provide a visual demonstration of the fact that odd halfinteger index disclinations must have their singular points at φ = 0 or φ = a/2. Because these singular points terminate layers and connect the two sets of parallel planes, it follows that the layers must arise as projections from the fixed point values of φ—in other words, the presence of charge 1/2 disclinations forces the defect cores to sit at φ = na/2, n ∈ ℤ, or the heights of the level sets and halfway between them. Properly, for smectics φ is an element of the orbifold S^{1}/ℤ_{2}, just as it was an element of S^{1} for directed lines.
The way in which these singularities can join depends on the geometry of the height representation. For instance, consider the dislocation and the pincement. We are used to the idea that a disclination dipole creates a dislocation (32). Indeed in Fig. 3A we depict a ±1/2 disclination pair that forms a charge +1 dislocation. However, consider the pincement shown in Fig. 3B. In this configuration the two singular points both lie at the same height and hence they can be cancelled. The smectic dislocation cannot be cancelled because the singular layers do not line up. Although this follows from the algebraic structure of paths on the Klein bottle, [i.e., (S^{1/2}FS^{−1/2})F^{−1} = S, but FF^{−1} = e], the height construction provides a geometric way to see the difference between these two complexions. We are also able to construct a dislocation directly. Because a Burgers circuit around a dislocation changes φ by a multiple of a, it follows that the height function must have the topology of a helicoid. As shown in Fig. 3C and D, the composite dislocation has the structure of two intertwined halfhelicoids. Note that the charge of a dislocation depends on both the sign of the helicoid (left or righthanded) and the tilt of the layers. A righthanded helicoid moved smoothly around a +1/2 disclination remains righthanded, but because the tilt of the layers has reversed, the charge of the dislocation changes sign. This is the move that generates the ambiguity modulo 2 for the total dislocation charge in the presence of a 1/2 index disclination (24). In order for dislocation/antidislocation pairs to form, a neck connecting adjacent surfaces must be formed. As long as this is smaller than the core size of the defects, it can be created without upsetting the overall topology. As the inchoate dislocation pair separates, these necks expand and, in a manner similar to continuous deformations of Riemann's minimal surface, break up into two oppositely handed helicoids (33), as depicted in Fig. 1C. Note that this type of process, which necessarily includes local surgery for combining defects, is implied in the theory of topological defects as well (9). In smectic systems, however, we must perform local surgery whenever we pull a defect through a layer—this follows from the above result that the defect cores must remain at a constant height.
Energetics from a Geometric Point of View
The height function approach provides a purely geometric way to formulate the nonlinear energetics of a twodimensional smectic. We start with the compression energy. In terms of the phase field the constraint of equal spacing is embodied in the condition that ∇φ = 1 (34). Thus any nonlinear strain u_{zz} must vanish when ∇φ = 1. We are free, however, to choose any rotationally invariant function with this property; theories based on different nonlinear strains will differ from each other in anharmonic powers of the strain, and only the coefficient of the harmonic term is definitionindependent. The height function defines a surface with a unit normal N = (−∂_{x}φ,−∂_{y}φ,1)/
Goldstone's theorem requires that the bending energy, on the other hand, reduce to ∇^{2}φ to harmonic order in φ. Possible nonlinear generalizations include the curvature of the smectic lines κ = ∇·(∇φ/∇φ). Note however, that the mean curvature of the surface defined by φ is H = (1/2)∂_{i}N^{i}, which reduces to the required expression to linear order in φ. We thus propose the following geometric free energy for the twodimensional smectic: where B is the bulk modulus and K_{1} is bending elastic constant. The free energy is defined as an integral over the height surface that we have introduced and is the Willmore functional in a field (35). The singularities associated with disclinations and dislocations specify the boundary conditions and topology for the associated variational problem.
Independent of the energy functional, we know that any equally spaced, twodimensional smectic has a corresponding height surface with constant N·ẑ. It follows that the Gauss map of such a surface sweeps out a latitude on the unit sphere, and so the Gaussian curvature K of the height surface must vanish. When K = 0, the surface must be locally isometric to the plane; the surface is either a cone, half of a cone attached to planes, planes, or the tangent developable surface generated by a cylindrical helix (36). Respectively, these four structures have disclination charges of +1, +1/2, 0, and +1. Thus, we provide a direct demonstration that it is only possible to have equally spaced layers for these three defect charges.
Discussion
We have outlined an approach to the study of topological defects in systems with broken translational invariance; topological equivalence requires more than just homotopic paths. One should consider homotopy classes of smooth(er) maps. Whereas the usual homotopy theory of defects would use continuous homotopy in the φ and θ components of the order parameter, we have argued that it is necessary to look at homotopy of φ(x) and ∇φ(x) together. Our method focuses on singularities and critical points in a phase field viewed as a height function over ℝ^{2}. Although the fundamental group of the ground state manifold can be constructed, it is known that when the loops involve both rotations and translations there are homotopy classes of loops implied that cannot be realized from complexions of the physical system. This arises from a mismatch of the dimensionality of the manifold of ground states with the number of Goldstone modes. To remedy this, we have constructed a local map from the configurations to the ground state manifold by employing Taylor series data at each point. We have shown that the homotopy theory of defects only works for singularities in φ (dislocations) and not consistently for singularities in ∇φ. Although the critical points in our theory behave somewhat like disclinations, they are not; as we demonstrated, two +1 disclinations cannot be brought together.
In the case of directed lines, circlevalued maps of the sort we have described above are equivalent (ref. 29, in particular see p. 37) to closed forms (37). For unoriented lines, that line of formalism may best be extended by considering quadratic differentials instead (38). Combining these results with ours may be a way of extending our work to the topology of such patterns on surfaces of nontrivial topology (rather than just the plane). The layer normal of smectics on curved surfaces with planar topology can develop caustics (34); extending our approach here to those systems is complicated by this additional holonomy and how it interacts with nonorientability near 1/2 disclinations.
Although we have focused here on twodimensional smectics and directed lines, our approach easily generalizes to threedimensional smectics that can be viewed as a height function over ℝ^{3} into S^{1}/ℤ_{2}. From this perspective the topological equivalence of edge and screw dislocations becomes particularly vivid; within the fourdimensional space spanned by (x,y,z,φ) the former are helicoids in x,y,φ for a fixed z, whereas the latter are helicoids in x,y,z at constant φ. Likewise, the theory of defects in solids can be formulated as maps from ℝ^{3} to T^{3}/𝕏, where 𝕏 is one of the 230 threedimensional space groups. The dislocations will continue to correspond to cycles on T^{3}, whereas the disclinations will be controlled by the critical points of these maps. Whether the study of singularities can lead to bounds on the total free energy of smectic complexions is an open question, as is the general set of rules for the combination of dislocations and disclinations.
Acknowledgments
We thank Jim Halverson, Elisabetta Matsumoto, Carl Modes, and Ari Turner for stimulating discussions and Koenraad Schalm for bringing ref. 26 to our attention. This work was supported in part by National Science Foundation Grants DMR0547230 and DMR0520020 and by gifts from L. J. Bernstein and H. H. Coburn.
Footnotes
 ^{1}To whom correspondence should be addressed. Email: kamien{at}physics.upenn.edu

Author contributions: B.G.C., G.P.A., and R.D.K. designed research, performed research, and wrote the paper.

The authors declare no conflict of interest.

This article is a PNAS Direct Submission.

↵* The astute reader might wonder about the restriction of paths for a spontaneously broken gauge theory where the Goldstone mode is also absent. There are no restrictions: In the gauged case, the winding is not a nontrivial loop in a ground state manifold connecting equivalent but distinct ground states, but rather a winding on a gauge orbit connecting different representations of the same physical state.

↵† Here and throughout ∇ acts on the twodimensional xyplane.

↵‡ More precisely, we are considering the njet of the configuration at P. Strictly speaking, for 1jets, we should keep track of the magnitude ∇φ_{x=P} as well. But because ∇φ is an element of ℝ, which is contractible, the topological behavior of this component would be trivial. An argument of Kléman, Michel, and Toulouse (8) (suitably reinterpreted) shows that jets higher than first order are not needed for a similar reason.

↵§ Note that, although the branch point of a helicoid is quite apparent when viewed in three dimensions, in any given slice of it, the branch point is not apparent because the layers all remain smooth; see, particularly, images in ref. 19.

¶ There may be other types of deformations that we choose to exclude due to energetic or other physical concerns. For instance, it might be convenient to consider deformations up to those that preserve the number of layers, which would prevent us from pushing the surfaces in the vertical direction.

Freely available online through the PNAS open access option.
References
 ↵
 Taylor GI
 ↵
 ↵
 ↵
 Muševič I,
 Scaron;karabot M,
 Tkalec U,
 Ravnik M,
 Zcaron;umer S
 ↵
 ↵
 Volovik GE,
 Mineev VP
 ↵
 Fichera G
 Rogula D
 ↵
 ↵
 ↵
 ↵
 ↵
 ↵
 Nagel L,
 Stein D
 Sethna JP,
 Huang M
 ↵
 Chaikin PM,
 Lubensky TC
 ↵
 Sethna JP
 ↵
 ↵
 ↵
 ↵
 Nye JF,
 Berry MV
 ↵
 Harrison C,
 et al.
 ↵
 ↵
 ↵
 ↵
 ↵
 Kléman M,
 Lavrentovich OD
 ↵
 ↵
 Golubitsky M,
 Guillemin V
 ↵
 Milnor J
 ↵
 Farber M
 ↵
 Goresky M,
 MacPherson R
 ↵
 Balian R,
 Kléman M,
 Poirier JP
 Berry MV
 ↵
 ↵
 Colding TH,
 Minicozzi WP, II
 ↵
 Santangelo CD,
 Vitelli V,
 Kamien RD,
 Nelson DR
 ↵
 Willmore T
 ↵
 Nistor AI
 ↵
 ↵
 Strebel K
Citation Manager Formats
Sign up for Article Alerts
Article Classifications
 Physical Sciences
 Physics