ISRN Mathematical Physics

Volume 2012 (2012), Article ID 601749, 16 pages

doi:10.5402/2012/601749

Research Article

Center-Vortex Loops with One Self-Intersection

Julian Moosmann and Ralf Hofmann

Institut für Theoretische Physik, Universität Karlsruhe (TH), Kaiserstrasse 12, 76131 Karlsruhe, Germany

Received 7 December 2011; Accepted 1 April 2012

Academic Editors: E. Akhmedov and S. Ketov

Copyright © 2012 Julian Moosmann and Ralf Hofmann. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We investigate the 2D behavior of one-fold self-intersecting, topologically stabilized center-vortex loops in the confining phase of an SU(2) Yang-Mills theory. This coarse-graining is described by curve-shrinking evolution of center-vortex loops immersed in a flat 2D plane driving the renormalization-group flow of an effective “action.” We observe that the system evolves into a highly ordered state at finite noise level, and we speculate that this feature is connected with 2D planar high superconductivity in FeAs systems.

1. Introduction

The idea of a nontrivial ground state being responsible for the emergence of “elementary” particles is a rather old one: already Lord Kelvin proposed that atoms and molecules should be considered knotted lines of vortices representing distortions in a universal medium (or ground state)—the ether [1]. As we know now, the physics of atoms and molecules described in terms of a much more efficient and elegant framework quantum mechanics. The agent responsible for the chemical bond—Lord Kelvin's electron—is considered a spinning point particle in quantum mechanics, and this yields an excellent description of atomic physics, collider physics, and in the bulk of condensed matter physics.

There are, however, theoretical discrepancies with the concept of the electron being a point particle, and there are exceptional experimental situations pointing to the limitations of this concept to describe reality. As for the former, we have the old problem of a diverging classical self-energy not resolved in quantum electrodynamics where the electron mass is introduced as a free parameter whose running with resolution needs an experimental boundary condition. On the other hand, the two-dimensional dynamics of strongly correlated electrons in condensed matter physics signals the relevance of nonlocal effects possibly related to the nontrivial anatomy of the electron becoming relevant in collective phenomena [2–4]. Also, recent high-temperature plasma experiments indicate unexpected explosive behavior not unlikely related to the mechanism for lepton emergence, see [5] and references therein.

Early papers, mostly investigating on a 4D Euclidean lattice the role of center vortices in forming a confining ground state at low temperatures/resolution [6–11], are based on the definition of a dual order parameter for quark confinement by [12]. Recent developments in understanding the confining phase of an SU(2) Yang-Mills theory suggest that Lord Kelvin's ideas may actually be realized in Nature, see also [13–16]. The authors of [13–16] construct a plausible effective low-energy action for the 4D SU(2) Yang-Mills theory with solutions to the associated field equations representing closed confining strings knotted into stable solitons. In the thermodynamic approach of [17] the emergence of magnetic center-vortex loops (CVLs) is related to discontinuous phase changes of a complex order parameter for confinement across the (downward) Hagedorn transition and the fact that no magnetic charges exist where these flux lines could end. Also, it was discussed in [17] how the locations of topologically stabilized self-intersection represent isolated, spinning magnetic charges. (Notice that with respect to the electromagnetic U(1) of the Standard Model there is a dual interpretation of magnetic charges emerging in an SU(2) Yang-Mills theory).

In our previous article [18] we have investigated the sector with self-intersections by considering a resolution-dependent ensemble average. The corresponding weight-functional is defined purely in terms of the planar curves's geometry. The resolution dependence of this geometry, in turn, is determined by a curve-shrinking equation (heat-equation) [19, 20]. The validity of this description of spatial coarse-graining is motivated by considerations relating local curvature with the direction and speed of “motion” of the associated line-segment. The requirement that the partition function over a given ensemble of planar curves is invariant under a change of the resolution then yields the renormalization-group evolution of the weight-functional which is written as the exponential of an “action.” Here the term “action” is slightly misleading since we do not aim at describing the time-evolution of the system by demanding stationarity of the “action” under curve variation. To do the latter, a model, which relates resolution and time (being a macroscopic concept associated with the measuring apparatus), needs to be introduced. We thus regard resolution over time as the more fundamental quantity to describe certain subatomic systems. Our observation is that the effective “action” exhibits a transition towards dilational invariance after a finite, critical decrease of resolution. On average, CVLs with are shrunk to circular points for a resolution less than the critical value which de facto removes them from the spectrum and thus generates an asymptotic mass gap. (CVLs with are massive [17, 21]). Knowing the evolution of the weight-functional, one is in a position to compute the resolution dependence of “observables” as ensemble averages of the associated (nonlocal or local) “operators.” As for the evolution of the initially sharp center-of-mass position, we observe a spread of the variance with decreasing resolution saturating at a finite value. This is similar to the unitary free-particle evolution of a position eigenstate in quantum mechanics.

The purpose of the present paper is to extend the procedure of [18] to the case of . We now have a singled-out point on the curve: the location of the self-intersection where practically the entire mass of the soliton resides [17]. Setting the Yang-Mills scale of the SU(2) theory equal to the electron mass keV, which in turn determines the mass of the intersection point, we interpret this soliton as an electron or a positron [5, 21]. In the presence of a static electric or magnetic background field it is physically possible to lift the two-fold degeneracy w.r.t. the two possible directions of center-flux: the soliton exhibits a two-fold spin degeneracy. Notice that as long as both wings of center flux are of finite size the position of the intersection point can be shifted at almost no cost of energy. In particular, if the inner angle between in- and outgoing center-flux at the intersection is sufficiently small, then a motion of points on the vortex line directed perpendicular to the bisecting line of the angle easily generates a velocity of the intersection point which exceeds the speed of light, see Figure 1. Recall that the path-integral formulation of quantum mechanics admits such superluminal motion in the sense that the according trajectories sizably contribute to transition amplitudes.

Figure 1: Points on the center flux lines moving oppositely on a line perpendicular to the bisecting line of the angle with velocity modulus . For sufficiently small the velocity modulus of the intersection point is superluminal: .

The paper is organized as follows. In Section 2 we discuss the physics associated with the emergence of topologically stabilized CVLs with intersection number , and how their spatial 2D coarse-graining is captured by a curve-shrinking flow. Some mathematical results on the properties of this flow for immersed curves, which are relevant for our subsequent numerical analysis, are briefly discussed. Also, we repeat our discussion in [18] of how the renormalization-group flow of an effective “action” is driven by the curve-shrinking evolution of the members of a given ensemble of curves. In Section 3 we explain our numerical analysis concerning the computation of the effective “action,” the variance of the location of the self-intersection, and the entropy associated with a given ensemble. Finally, in Section 4 we summarize our results and interpret them in view of certain 2D layered, quasimetallic systems exhibiting high- superconductivity.

2. Conceptual Framework

2.1. Self-Intersecting Center-Vortex Loops

The transition from the non-self-intersecting to the self-intersecting CVL sector is by twisting of non-self-interesecting curves. The emergence of a localized (anti)monopole in the process is due to its capture by oppositely directed center fluxes in the intersection core (eye of the storm). By a rotation of the left half-plane in Figure 2(a) by an angle of , see Figure 2(b), each wing of the CVLs forms a closed flux loop by itself thereby introducing equally directed center fluxes at the intersection point. This does not allow for an isolation of a single, spinning (anti)monopole in the core of the intersection and thus is topologically equivalent to the untwisted case Figure 2(a). However, another rotation of the left-most half-plane in Figure 2(c) introduces an intermediate loop which by shrinking is capable of isolating a spinning (anti)monopole due to oppositely directed center fluxes. Notice that in the last stage of such a shrinking process (short distances between the cores of the flux lines), propagating dual gauge modes are available. (On large distances these modes are infinitely massive which is characteristic of the confining phase): There is repulsion due to Biot-Savart which needs to be overcome. This necessitates an investment of energy manifesting itself in terms of the mass of the isolated (anti)monopole (eye of the storm). Alternatively, the emergence of an isolated (anti)monopole is possible by a simple pinching of the untwisted curve, again having to overcome local repulsion in the final stage of this process.

(a)

(b)

(c)

(d)

Figure 2: (Topological) transition from the sector (a), (b), (c) to the sector (d) by twisting and subsequent capture of a magnetic (anti)monopole in the core of the final intersection. Arrows indicate the direction of center flux.

For the analysis performed in the present work we solely regard the situation depicted in Figure 2(d) and thus no longer need to discuss the direction of center flux within a given curve segment. This is not relevant for the process of a spatial coarse-graining microscopically described by the same curve-shrinking flow as applied to sector with [18].

2.2. Euclidean Curve Shrinking Flow

The restriction of evolution of a CVL with self-intersection number to the plane is an essential constraint on generality if we aim at a fundamental description of the effectively quantum mechanical behavior of a charged lepton (electron, muon, tau-lepton). An electron bound inside a hydrogen atom certainly “moves” in 3D and quantum mechanics describing it as a spinning, relativistic point particle is in accurate agreement with experiment. The quantum mechanically incompletely understood condensed-matter physics of strongly correlated 2D electrons, however, is a possibly fertile application field of our restriction of curve evolution to the plane, see Section 4.

Notice that by immersing an SU(2) CVL with finite core size and mass of the dual gauge field into a flat 2D surface at , , a hypothetic observer measuring a positive (negative) curvature of a segment of the vortex line experiences more (less) negative pressure in the intermediate vicinity of this curve segment leading to its motion towards (away from) the observer, see Figure 3.

Figure 3: Highly space-resolved snapshot of a CVL segment. The pressure in the region pointed to by the normal vector is more negative than the pressure thus leading to a motion of the segment along .

The (inward directed) speed of a point in the core of the vortex will be a monotonic function of the curvature at this point. On average, this shrinks the CVL. Alternatively, one may globally consider the limit , , that is, the confining phase of an SU(2) Yang-Mills theory, but now take into account the effects of an environment which locally relaxes this limit (by collisions) and thus also induces curve shrinking. This situation is described by the following equation for a flow in the (dimensionless) parameter :

where is arc length, is a point on the CVL in the plane, and is a string tension effectively expressing the distortions induced by the environment. After a rescaling, , , (2.1) assumes the following form

where is a (dimensionless) curve parameter, the (inward-pointing) Euclidean unit normal, the scalar curvature, defined as

, and denotes the Euclidean scalar product of the vectors and . In the following we resort to a slight abuse of notation by using the same symbol for the functional dependence on or .

We now consider curves with one self-intersection, that is, , in the sense of the stable situation of Figure 2(d). This situation was mathematically analysed in [22]. Since the direction of center flux is inessential for the shrinking process we may actually treat this situation in a way as depicted in Figure 2(b), where the curve is defined to be a smooth immersion into the plane with exactly one double point and a total rotation number zero, . Here the (dimensionless) curve length is given by the smooth integration . Notice that this is topologically distinct from the case Figure 2(d) where one encounters a nonvanishing rotation number which is not smoothly deformable to zero.

In the case a smooth, embedded curve shrinks to a circular point under the flow for [19, 20]. That is, the isoperimetric ratio approaches from above. The curve in situation Figure 2(b) separates the plane into three disjoint areas two of which are finite and denoted by and . We understand by the finite, critical value of where either or or both vanish. This corresponds to a singularity encountered and thus terminates the flow.

Recall that in the case the rate of area change is a constant, . This is no longer true for . However, we have that

Also, for we have in comparison to the case the more relaxed constraint that .

In contrast to the case the isoperimetric ratio for the case is bounded for if and only if . Notice that the case physically is extremely fine-tuned.

2.3. Effective “Action”

We now wish to interpret curve-shrinking as a Wilsonian renormalization-group flow taking place in the  CVL sector in the sense defined in Section 2.2. A partition function, defined as a statistical average (according to a suitably defined weight) over  CVLs, is to be left invariant under a decrease of the resolution determined by the flow parameter . Notice that, physically, is interpreted as a strictly monotonic decreasing (dimensionless) function of a ratio where () are mass scales associated with an actual (initial) resolution applied to the system. The role of can also be played by the finite temperature of a reservoir coupled to the system.

To devise a geometric ansatz for the effective “action” , which is a functional of the curve representable in terms of integrals over local densities in (reparametrization invariance), the following reflection on symmetries is in order. (i) Scaling symmetry : for , implying at fixed , the “action” should be invariant under further finite rescalings (decoupling of the fixed length scales and ), (ii) Euclidean point symmetry of the plane (rotations, translations, and reflections about a given axis): Sufficient but not necessary for this is a representation of in terms of integrals over scalar densities w.r.t. these symmetries. That is, the “action” density should be expressible as a series involving products of Euclidean scalar products of , or constancy. However, an exceptional scalar integral over a nonscalar density can be devised. Consider the area , calculated as

The density in (2.5) is not a scalar under translations.

We now resort to a factorization ansatz as

where in addition to Euclidean point symmetry () is (is not) invariant under . In principle, infinitely many operators can be defined to contribute to . Since the evolution homogenizes the curvature except for a small vicinity of the intersection point higher derivatives of w.r.t. should not be of importance. We expect this to be true also for Euclidean scalar products involving higher derivatives . To yield conformally invariant expressions such integrals need to be multiplied by powers of and/or or the inverse of integrals involving lower derivatives. At this stage, we are not capable of constraining the expansion in derivatives by additional physical or mathematical arguments. To be pragmatic, we simply set equal to the isoperimetric ratio:

We conceive the nonconformal factor in as a formal Taylor expansion in inverse powers of or due to the property of conformal invariance for .

Since we regard the renormalization-group evolution of the effective “action” as induced by the flow of an ensemble of curves, where the evolution of each member is dictated by (2.2), we allow for an explicit dependence of the coefficient of the lowest nontrivial power or . In principle, this sums up the contribution to of certain higher-power operators which do not exhibit an explicit dependence. Hence we make the following ansatz:

The initial value is determined from a physical boundary condition such as the mean length at .

2.4. Geometric Partition Function

Let us now numerically investigate the effective “action” resulting from a partition function w.r.t. a nontrivial ensemble . The latter is defined as the average

over the ensemble . Let us denote by an ensemble consisting of curves where is obtained from by adding a new curve . The effective “action” in (2.6) (when associated with the ensemble we will denote it by ) is determined by the function , compare with (2.8), whose flow follows from the requirement of -independence of :

This is an implicit, first-order ordinary differential equation for which needs to be supplemented with an initial condition . A natural initial condition is to demand that the quantity

coincides with the algebraic mean defined as

From a value for follows. We also have considered a modified factor

While, due to (2.13), the ansatz for the geometric effective "action" in (2.6) thus is profoundly different for such a modification of physical results such as the evolution of the variance of the intersection agree remarkably well, see Section 3.

3. Results of Simulation

3.1. Preparation of Ensembles

Similar as in [18] we normalize all curves to have the same initial area , and, since we are now interested in the position of the intersection where the (anti)monopole, is localized, we have applied a translation to each curves in the ensembles such that the location of the intersections initially coincides with the origin.

Since the critical value of the flow parameter varies from curve to curve, we order the members of the maximal-size ensemble into subensembles such that . The types of ensembles obtained in this way are referred to as -ordered. We also have performed all simulations with ensembles whose members are picked randomly from and have obtained strikingly similar results for ensemble averages of “observables” using and for the evolution to the left of .

The maximal-size ensemble at is depicted in Figure 4 with the universal choice . The curves in Figure 4 are arranged in a -ordered way. We have . In Figure 5 the evolution of an initial curve under curve shrinking is shown from two viewpoints. The flow is started at and stopped at a value of shortly below . In Figure 6 the flow of the intersection points, corresponding to the initial curves depicted in Figure 4, is shown. The search for solutions to the second-order partial differential equation (PDE) (2.2) subject to periodic boundary conditions in the curve parameter, , and for the initial conditions depicted in Figure 4 was performed numerically using the method of lines. That is, the PDE was discretized on a uniform grid in the parameter yielding a semidiscrete problem in terms of a system of ordinary differential equations (ODEs) in which was solved using Mathematica. Figure 5 indicates why this technique is called the numerical method of lines. As one can also see from Figure 5, a set of discrete points on the curve, although remaining equidistant in , may evolve under the flow such that the spatial distances between next-neighbours-points fall below the numerical precision. Numerically, the flow then encounters a singularity (not to confused with the earlier mentioned nonfictitious singularities). To recognize such a situation automatically, (2.4) was exploited: the evolution was stopped as soon as a sizable deviation occured from what (2.4) predicts. The configuration obtained at this point in was fitted in such a way that a new discretization in yielded well-separated points to restart the methods of lines. (2.4) was also used as an indicator for the final singularity at where or or both vanish.

Figure 4: Initial curves contributing to the ensemble . The intersection points coincide with the origin, and all curves have the same area . By definition is -ordered.

Figure 5: Plot of the evolution of an CVL (curve 12 of Figure 4) under (2.2). The thick central line indicates the trajectory of the intersection point which coincides with the origin at .

Figure 6: Flow of the intersection points for the initial curves depicted in Figure 4.

3.2. Renormalization-Group Invariance of Partition Function

For all ensembles , the dependence of the coefficient in (2.8) roughly behaves like a square root , where is the weakly ensemble-dependent minimal resolution. For the modified “action” the coefficient is well approximated by a linear function . Again, denotes a weakly ensemble-dependent minimal resolution. For -ordered ensembles, the results for for the “actions” (2.8) and (2.13) are shown in Figures 7 and 8, respectively. The results for ensembles do not differ sizably from those presented in Figures 7 and 8.

Figure 7: The squares of the coefficients entering the ansatz for effective “action” of (2.6) specializing to (2.8) for -ordered ensembles up to .

Figure 8: The coefficient entering the ansatz for the effective “action” of (2.6) specializing to (2.13) for -ordered ensembles up to .

3.3. Variance of Location of Self-Intersection

The mean intersection over the ensemble is defined as

where is the location of self-intersection (intersection point) of curve at . The scalar statistical deviation of over the ensemble is defined as

where

and similarly for the coordinate . In Figure 9 plots of are shown when evaluated over the ensembles subject to the “action”

and the initial condition . In Figure 10 the according plots of are depicted as obtained with the “action”

and subject to the initial condition . Relaxing the constraint of -ordering () does not entail a qualitative change of the results.

Figure 9: Plots of for the -ordered ensembles with . We have employed the ansatz for the “action” .

Figure 10: Plots of for the -ordered ensembles with . We have employed the ansatz for the “action” .

The results presented in Figures 9 and 10 are unexpected since in the sector the variance of the “center-of-mass” saturates rapidly to finite values. In contrast, for the sector the variance of the location of the self-intersection initially increases, reaches a maximum, and decreases to zero at a finite value of . This is readily confirmed by the evaluation of the entropy, see Section 3.4.

3.4. Evolution of Entropy

Let us now evaluate the flow of entropy defined as

where is given by (2.6). In Figure 11 plots are shown for . when evaluated with the “action” for -ordered ensembles of size . These graphs look very much alike to the ones generated using the “action” . Notice the continuous approach to zero at finite values of . This implies that order emerges spontaneously in the system with decreasing resolution: starting at a finite value of , a particular member of is singled out by its weight approaching unity. Judging from our results for the sector [18], this behavior is highly unexpected. Therefore the nontrivial topology of induces qualitative differences to the coarse-graining process.

Figure 11: Flow of the entropies for -ordered ensembles of size when evaluated with the “action” . The situation does not change qualitatively if the “action” is used.

4. Summary, Interpretation of Results, and Conclusion

In this paper we have investigated the spatial coarse-graining of CVLs, immersed in a flat 2D plane, of an SU(2) Yang-Mills theory being in its confining phase. The focus was on the sector with one topologically stabilized self-intersection (existence of an isolated magnetic charge at its location, ). We have analysed this coarse-graining process in terms of curve ensembles generated by evolving an initial situation under the curve-shrinking flow [19, 20, 22]. The idea here is to suppose that curve shrinking in the parameter represents an exact coarse-graining of a given initial state and to reconstruct the associated ensemble-weight of the statistical approach (exponential of effective “action”) by demanding invariance of the corresponding partition function under the flow in (renormalization-group evolution). Notice that is related to a physical resolution, applied to probing the system, in a strictly monotonic decreasing manner. This resolution may be associated with a local momentum transfer exerted by an observer or a globally defined temperature of an environment. The functional dependence of on these physical parameters depends on the given experimental situation. It is, however, reasonable to assume that finite values of universally correspond to finite values of these physical parameters.

In Sections 3.3 and 3.4 we have obtained the unexpected result that a statistical ensemble of renormalization-group evolved curves spontaneously orders itself in the sense that, starting from finite values of , only a particular member of the ensemble survives the process of 2D spatial coarse-graining. That is, the entropy attributed to the ensemble is practically zero for sufficiently large values of . For the location of self-intersection (charge of an electron) this means that no dissipation of energy, provided by the environment, can be mediated by the monopole situated within the core of the intersection if the resolution falls below a critical, finite value. This result must drastically depend on the two-dimensionality of space and the fact that we consider the sector with , compare with [18].

The recently discovered, unconventional systems do not appear to exhibit an explicit, strong correlation between the electrons contained in their theoretically suggested, 2D-superconducting layers, see [23] and references therein. If the two-dimensional behavior of noninteracting electrons, subject to an environment represented by the parameter , indeed is described by the coarse-graining process investigated in the present work, then the sudden decrease of entropy that we observe at a finite value of should ultimately be connected to this particular kind of high- superconductivity. Here is a monotonically decreasing function of temperature.

Acknowledgment

The authors would like to thank Francesco Giacosa and Markus Schwarz for useful conversations.

References

W. T. Kelvin and P. G. Tait, Treatise on Natural Philosophy, vol. 2, Cambridge University Press, 1867.

J. G. Bednorz and K. A. Müller, “Possible high superconductivity in the Ba-La-Cu-O system,” Zeitschrift für Physik B Condensed Matter, vol. 64, no. 2, pp. 189–193, 1986.

P. W. Anderson, “Present status of the theory of high cuprates,” http://arxiv.org/abs/cond-mat/0510053

P. W. Anderson, “Twenty years of talking past each other: the theory of high ,” Physica C, vol. 460-462, pp. 3–6, 2007.

F. Giacosa, R. Hofmann, and M. Schwarz, “Explosive Z pinch,” Modern Physics Letters A, vol. 21, no. 36, pp. 2709–2715, 2006.

G. Mack and V. B. Petkova, “Comparison of lattice gauge theories with gauge groups and ,” Annals of Physics, vol. 123, no. 2, pp. 442–467, 1979.

G. Mack, “Predictions of a theory of quark confinement,” Physical Review Letters, vol. 45, no. 17, pp. 1378–1381, 1980.

G. Mack and V. B. Petkova, “Sufficient condition for confinement of static quarks by a vortex condensation mechanism,” Annals of Physics, vol. 125, no. 1, pp. 117–134, 1980.

G. Mack and E. Pietarinen, “Monopoles, vortices and confinement,” Nuclear Physics B, vol. 205, no. 2, pp. 141–167, 1982.

H. B. Nielsen and P. Olesen, “A quantum liquid model for the QCD vacuum. Gauge and rotational invariance of domained and quantized homogeneous color fields,” Nuclear Physics B, vol. 160, no. 2, pp. 380–396, 1979.

E. Tomboulis, “'t Hooft loop in lattice gauge theories,” Physical Review D, vol. 23, no. 10, pp. 2371–2383, 1981.

G. 't Hooft, “On the phase transition towards permanent quark confinement,” Nuclear Physics B, vol. 138, no. 1, pp. 1–25, 1978.

L. Faddeev and Antti J. Niemi, “Stable knot-like structures in classical field theory,” Nature, vol. 387, no. 6628, pp. 58–61, 1997.

L. Faddeev and A. J. Niemi, “Partially dual variables in Yang-Mills theory,” Physical Review Letters, vol. 82, no. 8, pp. 1624–1627, 1999.

L. Faddeev and A. J. Niemi, “Aspects of electric and magnetic variables in Yang-Mills theory,” Physics Letters. B, vol. 525, no. 1-2, pp. 195–200, 2002.

L. Faddeev and A. J. Niemi, “Spin-charge separation, conformal covariance and the Yang-Mills theory,” Nuclear Physics B, vol. 776, no. 1-2, pp. 38–65, 2007.

R. Hofmann, “Nonperturbative approach to Yang-Mills thermodynamics,” International Journal of Modern Physics A, vol. 20, no. 18, pp. 4123–4216, 2005, Erratum vol. 21, pp. 6515–6523, 2006

J. Moosmann and R. Hofmann, “Evolving center-vortex loops,” http://arxiv.org/abs/0804.3527

M. Gage and R. S. Hamilton, “The heat equation shrinking convex plane curves,” Journal of Differential Geometry, vol. 23, no. 1, pp. 69–96, 1986.

M. A. Grayson, “The heat equation shrinks embedded plane curves to round points,” Journal of Differential Geometry, vol. 26, no. 2, pp. 285–314, 1987.

R. Hofmann, “Yang-Mills thermodynamics,” http://arxiv.org/abs/0710.0962

M. A. Grayson, “The shape of a figure-eight under the curve shortening flow,” Inventiones Mathematicae, vol. 96, no. 1, pp. 177–180, 1989.

H.-H. Klauss and B. Büchner, “Neuer Goldrausch in der Supraleitung,” Physik Journal, vol. 7, p. 18, 2008.