Conformal Villarceau Rotors

We consider Villarceau circles as the orbits of specific composite rotors in 3D conformal geometric algebra that generate knots on nested tori. We compute the conformal parametrization of these circular orbits by giving an equivalent, position-dependent simple rotor that generates the same parametric track for a given point. This allows compact derivation of the quantitative symmetry properties of the Villarceau circles. We briefly derive their role in the Hopf fibration and as stereographic images of isoclinic rotations on a 3-sphere of the 4D Clifford torus. We use the CGA description to generate 3D images of our results, by means of GAviewer. This paper was motivated by the hope that the compact coordinate-free CGA representations can aid in the analysis of Villarceau circles (and torus knots) as occurring in the Maxwell and Dirac equations.


Villarceau Circles
In two recent papers, by Consa [3] and Hestenes [8], models of the electron are proposed that involve it moving rapidly internally around a torus. The parametrization of that motion in both papers is done by the two uniformly increasing angles evolving from the usual Euclidean image of generating a torus as a circle moving around a circle. In [9], Hestenes derives this form explicitly as a solution of the Dirac equation when a Lorentz boost is factored out.
A second physical situation in which orbits on a torus occur is in the knotted solutions to Maxwell's equation in vacuum, in which the Poynting vector of the electromagnetic field lies on a torus knot [2]. Such a knot is a   [13] geometric object, not a path, and it no longer requires a parametrization to describe its time dependence-the field lines themselves evolve with the speed of light. Any convenient parametrization works, though the innate conformal symmetry of the Maxwell equation may prefer a conformal parametrization.
In both applications, the simplest orbits or field lines are the (1, 1) torus knots. In this paper, we analyze their remarkable properties (they are circles!) by means of our representation from [5], where we encountered (m, n) torus knots as the orbits of certain conformal rotors. The natural parametrization of that solution was based on an orthogonal (commuting) decomposition, which involves one Euclidean rotation around the torus axis, and a 'conformal rotation' in the orthogonal direction. The resulting motion differs intriguingly from the Consa/Hestenes parametrization, and the canonical orbits of the two parametrizations are not identical.
In mathematics, the path that the (1, 1) knot on the torus traces is also known as a Villarceau circle. 1 Villarceau circles are usually introduced as two intersecting circles that are the cross-section of a torus by a well-chosen plane cutting it (Fig. 1). Picking one such circle and rotating it around the torus axis, the resulting family of circles can be used to rule the torus. By nesting tori smartly, the collection of all such circles then form a Hopf fibration of 3D space. In the theory of Penrose twistors, the Villarceau circles occur as the Robinson congruence, the family of light rays representing a non-null twistor; this is yet another link to interesting physics (see e.g. [1] for a geometric algebra treatment).
In the present paper, we prefer to consider the Villarceau circle as the (1, 1) torus knot rather than as a planar cut, and we generate it by means of a composite rotor in conformal geometric algebra in Sect. 2. We then characterize any Villarceau circle compactly as the orbit of a location-dependent simple rotor (after studying those in detail in Sect. 3) passing through a given point, with a varying speed determined by this rotor (Sect. 4). That compact description permits the derivation of a number of geometric properties (including the planar cut construction), many involving a specific imaginary sphere (Sect. 5), which is a natural completion of the conformal basis. The Hopf fibration then follows simply from letting the rotor act on arbitrary points, and by means of spherical inversion we can relate this directly to isoclinic rotations on the Clifford torus in 4D (in Sect. 6). The conformal nature of the construction allows a straightforward extension to Villarceau circles on arbitrary Dupin cyclides, in Sect. 7.

Rotor Generation
In [5] we showed how 3D conformal motions can be generated as the orbit of a rotor that is the exponential of the sum of two commuting 2-blades of CGA. These 2-blades geometrically represent 'orthogonal' point pairs (with perpendicularity properties described that paper). In the case of the torus, one point pair is the dual of a line, the other an imaginary point pair on that line. Both point pairs have a negative square, so that each generates periodic orbits.
As our prototype for this situation we can take B − = e 1 e 2 and B + = e 3 e + , in the usual orthonormal basis {e 1 , e 2 , e 3 , e + , e − } of the R 4,1 Minkowski space for the description of 3D conformal geometric algebra. Boldly denoted are the basis vectors for the Euclidean space (which can be interpreted as orthogonal planes at the origin), and e + and e − are the unit real sphere and the unit imaginary sphere, respectively. So e 2 − = −1, all other orthogonal basis vectors square to +1. For each of the 'real' (positively squared) basis elements, algebraic orthogonality is equivalent to geometric perpendicularity of the corresponding geometrical primitives at their intersection. For the 'imaginary' sphere, algebraic orthogonality with the real planes and the real sphere geometrically implies that the corresponding intersection should be an equator of the imaginary sphere (see [4] Chapter 14).
We can construct 2-blades from these basis vectors, and exponentiate them, as a canonical way of making rotors. Three different types of rotor result, depending on the sign of the square of the 2-blade. We illustrate the orbits of a point moving by each of them in Fig. 2. When the square is positive, the orbit has a source and a sink; when negative, it is periodic; when zero (for a blade containing a null vector like o = e − + e + or ∞ = e − − e + ), merely a sink. We will analyze these cases and their parametrizations in more detail in Sect. 3.
Since intersection of the dual plane e 3 with the dual unit sphere e + . This 'core' is not the Euclidean center of the circles, but acts as their conformal center. The simple rotors R − and R + can be composed into the product rotor where the addition of the exponent is permitted since the 2-blades commute. When moving in the product of the rotors, the resulting orbit resides on a torus, see [5]. When φ and ψ are taken in a ratio of ψ = ±(m/n) φ, the orbit is a periodic (m, n) torus knot. These knots are inscribed on the torus in a conformal manner, cutting all B + -circles of longitude (or B − -circles of latitude) by equal angles. This is illustrated in Fig. 3a (from [5]), which shows the R rotor orbit in red, and also indicates the R − -orbit (in blue) and the R + -orbit (in green). The knotted red orbit can be seen as locally composed of small steps in each of those circular orbits.
To obtain a Villarceau circle through a point x, we choose ψ = ±φ, generating a (1, 1) torus knot depicted in Fig. 3b. Let us define B as the sum of the unit 2-blades: This B is not a 2-blade (its square B 2 = B 2 − + B 2 + + 2 B − ∧ B + is not a scalar), so it cannot be factored by the outer product. Now form a rotor by exponentiation of B, apply it to a point x: and investigate its orbit-which should be a Villarceau circle. But this is not obvious! The orbit of a simple conformal rotor based on a 2-blade rotor is always a circle, as we illustrated in Fig. 2 (with a line, the orbit of a translation rotor, considered as a circle passing through infinity). The orbit of a rotor based on a general 2-vector, such as B above, is not. To demonstrate that in this exceptional case we still obtain a circle, we show that its orbit x ∧ B is identical to that of a 2-blade (which will depend on x)-even though B is not a 2-blade itself. Before we can do so, we first need to expose the simple rotors and their parametrizations in more detail.

Orbits of Simple Conformal Rotors
To analyze the orbit of the non-simple rotor above, let us briefly derive the properties and parametrization of simple rotors (i.e., exponentials of 2-blades) in conformal geometric algebra. Our main focus is on the exponential of an 'imaginary' 2-blade (i.e. one that squares to a negative scalar). But while we are at it, we can include real and null 2-blades with very little additional effort, even though their characteristics differ considerably (as Fig. 2 showed).

Imaginary Point Pair: Periodic Motions
Consider an imaginary point pair P (a 2-blade from CGA), normalized so that P 2 = −1. For invertible P , we can split x into its projection and rejection relative to P , and write using commutation and anti-commutation properties of those split terms with the P -rotor. For an imaginary point pair, the exponential involves trigonometric functions: exp(φP ) = cos(φ) + P sin(φ). Using those, x(φ) can be expanded to: The vectors (x · P ), (x · P )/P and (x ∧ P )/P form an orthogonal basis for the orbit, consisting of two real spheres (x · P ) and (x · P )/P , and an imaginary sphere (x ∧ P )/P (intersected by the real spheres in its equator). So using CGA, we can describe a moving point on a circle as the weighted sum of three specific orthogonal spheres. The parametrization by φ of the points along the circle is not-equidistant, but equiangular in a conformal manner. Since all the terms in Eq. (3) transform covariantly under conformal transformations, we can design a conformal mapping from a canonical situation that is particularly intuitive to parametrize, and then use that parametrization for the general case. Let us explain this by a specific example, illustrated in Fig. 4a, b.
The Euclidean rotation with rotation rotor determined by the 2-blade P = e 1 e 2 has a natural equi-angular parametrization relative to the center, the origin indicated by the conformal point vector o. The transformed points x(φ) = exp(−φP/2) x exp(φP/2) along the orbit x ∧ P then lie on rays o ∧ x(φ) ∧ ∞, and the angular parameter makes the points traverse this orbit at a uniform speed. Now transforming this situation by inversion in σ = e 2 + e + , the 2blade becomes e 1 e + , and its rotor represents a 'conformal rotation'. The This covariant relationship is illustrated in Fig. 4b, which shows how the Euclidean situation of Fig. 4a is conformally transformed by inversion in σ. Our figure is effectively a re-rendering of Picture 27 in Chapter 6 of the enjoyable book by Needham [10] in Visual Complex Analysis (though we do x Figure 4. a A Euclidean rotation, b a conformal rotation, obtained from a by an inversion. Equiangular rays from the conformal center are indicated in red. In b, (x · P )/P is an imaginary sphere, depicted as dotted (color figure online) not use complex numbers, and therefore are not limited to 2D). Figure 4b also illustrates how the orthogonal basis of Fig. 4a gets transformed in an orthogonal basis for the transformed circle. The imaginary sphere e − remains, the planes through x and the plane perpendicular to that become orthogonal spheres (still cutting the imaginary sphere e − in great circles).
This particular example is the typical situation: we can always design a conformal transformation to transform the Euclidean axis 2-blade to a desired conformal 2-blade, and the Euclidean equiangular parametrization then transforms to the equiangular parametrization from the transformed center. To transform e 1 e 2 to the product of two arbitrary orthogonal unitsquare spheres σ 1 σ 2 , one may use the method from [6] to construct the corresponding transformation versor as (σ 1 − e 1 ) e 2 + σ 2 (σ 1 − e 1 ).

A Real Point Pair: Exponential Motion from Source to Sink
The derivation of Eq. (2) also holds for a real point pair, where P 2 = 1. The expansion now involves hyperbolic functions; substituting exp(γP ) = cosh(γ) + P sinh(γ) we obtain Again the moving point is a weighted sum of elements (x · P ), (x · P )/P and (x ∧ P )/P , which form an orthogonal basis of three orthogonal spheres. See Fig. 5b. For the standard situation, we can actually use the same inversion as for the rotation, but now considering the radial motions. Figure 5a shows an exponential scaling relative to the flat center o ∧ ∞ = e + ∧ e − . The 2-blade for the radial rotor exp(−γ o ∧ ∞/2) transforms any point on a radial ray to a scaled location which is e −γ of the previous distance of the center (simultaneously affecting its weight by a factor e γ ). This rotor therefore provides Figure 5. a A Euclidean exponential scaling: P is a flat point, and the real and imaginary spheres x · P and (x · P )/P coincide. b A conformal scaling, obtained from a by an inversion. Circles of equal logarithmic increase are indicated in red (color figure online) a natural scaling parameter in the radial direction. The exponential allows the full range of γ ∈ (−∞, ∞) without leading to awkwardness at the (unreachable) origin, and therefore this parametrization by γ is more natural than the customary linear radial parameter of polar coordinates. Applying the inversion in the sphere σ = e 2 + e + , the 2-blade e + ∧ e − becomes the real point pair e 2 ∧ e − , and the radial motions become circular motions on circles emanating from the transformed origin. Thus Fig. 5a is transformed into Fig. 5b. Yet the parameter for 'natural' motion along the orbits of the transformed rotor exp(−γ e 2 ∧e − ) is still γ, to produce the transformed point −σ x(γ)/σ. That moving point departs from the source, and approaches the sink, in a hyperbolic pacing along the circular path, by Eq. (4).
An inversion in the origin sphere e + converts this situation to the 2-blade α o e 1 , which characterizes a transversion (in Fig. 6b), and the equidistant x (x · P )P located at P x · P P (b) (a) Figure 6. a A Euclidean translation, b a conformal transversion, obtained from a by an inversion planes e 1 + α∞ now become the spheres e 1 + 2αo passing through o with center e 1 /(2α), which cut the new orbit x ∧ o ∧ e 1 at the α-parametrized point x (α). Uniform translational speed towards infinity in the Euclidean situation naturally becomes a convergence to the sink of the transversion, by circular motions whose angle relative to the Euclidean center of the circular orbit can be shown to vary as the arctangent of the reciprocal of α.

Simplifying the Villarceau Rotor
To resume our tale after Eq. (1), we look for a simple rotor (i.e. the exponential of a 2-blade) of which the orbit through x is the same as that of exp(−Bφ/2), with B = B − + B + where B − and B + are commuting normalized 2-blades. We reiterate that this is not equivalent to rewriting the rotor as the exponential of a single 2-blade independent of its starting point; that cannot be done.
We will need a result from the Taylor expansion of a rotor (not necessarily simple) acting on a vector, in terms of its bivector A (from [4], Chapter 7): For simple rotors based on a 2-blade A, this series can be expressed in terms of trigonometric functions (when A 2 < 0), hyperbolic functions (when A 2 > 0) or linear functions (when A 2 = 0), as we showed in the previous section. The expansion then simplifies to the closed expression Eqs. (3), (4) or (5). However, in our case the rotor is not simple (since our bivector B is not a 2-blade), so we need to proceed more carefully.
We can use Eq. (6) for a small bivector A = B to understand the local orbit. The tangent vector to the orbit at x is: while the local osculating circle is: In our case, the bivector B equals B = e 1 e 2 + e 3 e + . For that particular rotor, the tangent circle is the actual full orbit. To show this, we compute the next term in the Taylor series of Eq. (6): After three terms, we are back at the first (with a minus sign); therefore the rest of the terms follows the same pattern. When substituting these results in the Taylor series Eq. (6), we can collate terms to express the series exactly on a basis of three elements, as: This is very similar to the general form of the orbit of an imaginary unit 2-blade as given in Eq. (3): x(φ) = (x ∧ P )/P + (x · P )/P cos φ + (x · P ) sin φ.
We observe that in the general form we can establish the original 2-blade of the rotor from the ratio of the factors of sine and cosine. Therefore in our case, the unit 2-blade generating the motion along the Villarceau circle is: since due to e − · (x · B) = 0 and x · (x · B) = 0, the geometric product can be replaced by an outer product. The factor x + (x · e − )e − is a real sphere, with center at 2 x/(1 − x 2 ) and with radius (1 + x 2 )/(1 − x 2 ) (where x is the Euclidean location of the conformal point x).
With B V (x) as its local 2-blade, the orbit at C(x) passing through x is therefore the circle.
Because of the final rewriting, it would seem that it should be simpler to use is a real point pair, with a source and a sink, and the parametric (time) dependency it generates along its orbit is therefore completely different (as we derived in Sect. 3). So if the parametrization is relevant, we should use B V rather than B V ; but for properties of the orbit as a whole (such as in the EM field lines of [2]), the simpler B V would be sufficient.
We conclude that the motion performed by x under the rotor with bivector B is equivalent to a pure conformal rotation induced by the imaginary point pair B V (x), and that it generates the Villarceau circle. This point pair B V (x) depends on x itself, so the motion by the rotor exp(B/2) is not globally equivalent to a single conformal rotation.

Geometric Properties
Since the Villarceau circles are orbits of a conformal rotor, they cut the parallels (x ∧ B − ) and meridians (x ∧ B + ) of the torus at a constant angle. This property is illustrated for the (1, 1) knot in Fig. 8. As a consequence of being conformal, all constructions and statements transform covariantly under conformal transformations. Not only can we translate, rotate and scale them from the standard situation, but using spherical inversion we can produce fully analogous properties of Villarceau circles on Dupin cyclides, as we will see in Sect. 7. For now we prefer to continue to use the torus vocabulary, as our token example of these conformally equivalent extensions. For the same reason, we use the unit radius sphere e + in B + = e 3 e + ; other scales can be produced by the usual conformal scaling rotors. The covariance of the conformal rotor formalism makes those extensions unnecessary to prove-they are automatic.  The other Villarceau circle through a point x is generated by reversing one of the rotors, say using R + R − , the exponential of the bivector B + − B − rather than B + + B − . The situation in Fig. 1 is rather special: there apparently exist points in which the two Villarceau circles passing through them are coplanar. We can easily establish where those points are, as follows. We demand that the two Villarceau circles x ∧ e − ∧ x · (B + + B − ) and x ∧ e − ∧ x · (B + − B − ) are coplanar (we can now use their simplified form with e − only, since their parametrization is immaterial to this property). This implies: (remember that x is the Euclidean part of the conformal point x). A trivial non-zero solution is x = CGA(α e 3 ) for any α, when the circles both degenerate to a line. But more interesting is the set of points such that x · e + = 0. Thus in our standard setup, any point x on the real unit sphere e + produces a planar pair of Villarceau circles. We then have to show that all Villarceau circles can be generated from a point on this unit sphere. And indeed they can: the meet of a Villarceau circle with the unit sphere at the origin is the point pair e + · x ∧ (x · B) ∧ e − , of which the square is ( . This is non-negative, so the point pair is real (with the only null degeneracy on the core circle of the set of Villarceau circles, where x 2 = 1 and x · e 3 = 0). Thus each Villarceau circle intersects the unit sphere, and each can therefore be generated from such an intersection point. In view of the pairing of Villarceau circles, starting them from a such point on the unit sphere e + is perhaps the most natural way to generate the circles in our standardized setup. Of course every circle intersects the unit sphere in two points, so in fact only half of the unit sphere is required; though we might use the other half to index Villarceau circles of opposite orientation.
We can obtain some more insight in the geometry of a Villarceau circle by substituting the series expansion terms of Eq. (9). Thus when used in the outer product, the second order term simplifies to This implies that the Villarceau circle can also be defined as 'tangent to the local tangent vector T (x) and intersecting the unit sphere at the origin in two antipodal points' (since it contains the imaginary unit sphere e − as a factor, see [4], Chapter 15 on the 'plunge' operation). This fact is visible in Fig. 7.
Actually, the circle C(x) is invariant under reflection in the imaginary sphere e − . To show this, first note that B changes sign under reflection in Then, rather surprisingly, we find that the sphere x · B (a sphere passing through x and orthogonal to the local tangent vector T (x)) is invariant under the reflection: since (x · B) · e − = 0. Then we can establish that C(x) is invariant modulo a sign In more detail, the inversion in e − exchanges the point x with a point x on the other side of the conformal center of the Villarceau circle (i.e. one of the locations where the real circle (B V ) * intersects the plane C(x) ∧∞), along an orthogonal circle passing through that conformal center and x. The resulting points x and x are symmetric relative to the origin: x ∧ e + ∧ x is a straight line. Inversion in an imaginary sphere like e − is easy to visualize: it is equivalent to an inversion in a real sphere of the same radius and center, followed by a point reflection in that center. For the origin-centered unit spheres of our standard situation, this follows from e − = (e − ∧ e + ) e + . When read as a compound operator for sandwiching, the factor e − ∧ e + = o ∧ ∞ is the flat point at the origin producing the point reflection. It is not a rotor, yet it can be used in the sandwiching product, see [4] Chapter 16.

L. Dorst
Adv. Appl. Clifford Algebras Figure 9. (Graphics and text from wikipedia page: Rotations in 4-dimensional space.) A 4D Clifford torus stereographically projected into 3D looks like a torus, and a double rotation can be seen as a helical path on that torus. For a rotation whose two rotation angles form a rational number, the paths will eventually reconnect, while for an irrational ratio they will not. An isoclinic rotation will form a Villarceau circle on the torus, while a simple rotation will form a circle parallel or perpendicular to the central axis

Space Filling Villarceau Circles
In the field generated by the rotor B, every point x picks out a Villarceau circle on a torus as its orbit. The torus is determined by the circles x ∧ B − and x ∧ B + , and the Villarceau circle is x ∧ B V (x), with B V (x) defined as in Eq. (11). Clearly the points of 3D space can be grouped in equivalence classes of the circles of points that are on the same Villarceau circle generated by B V (x); this class may be indexed by a point on the unit sphere e + (in a double cover, or in two opposite orientations). All Villarceau circles on the same torus are a larger equivalence class, rotationally invariant under the rotor R − , see Fig. 9, and each indexed by a point on a shared parallel of the unit sphere. Moving the index point x to other classes generates Villarceau circles on a different torus, involving the conformal rotation R + . All in all, a space-filling set of circles appears, together forming a Hopf fibration of 3-dimensional space, see

Isoclinic Rotations in 4D
As we have seen, there are two rotors involved in composing the rotations, in our specific example produced by the commuting 2-blades B − = e 1 e 2 and B + = e 3 e + . Since e 2 + = 1, we could consider the four vectors involved as just In that 4-dimensional space, a general rotation has two planes in which it needs to be specified, and we have chosen them orthogonally. But for the special rotor with 2-blade B = B + + B − , those rotations are coupled and performed at the same rate. This composite motion is called an isoclinic rotation, and normally one would consider it primarily in the 4D Euclidean space spanned by four orthogonal unit basis vectors {e 1 , e 2 , e 3 , e 4 }.
Our algebraically equivalent representation on the basis {e 1 , e 2 , e 3 , e + } also consists of four orthogonal unit vectors (though the full basis contains an additional negatively squared vector e − ), but now all can be drawn in 3D (as the three coordinate planes plus the unit sphere, all at the origin). Doing so apparently shows a toroidal image of the unit 3-sphere in 4D, converting the orbits of 3D points under an isoclinic rotation (which reside on that sphere) into the set of circles of the Hopf fibration.
In fact, our replacing e + by e 4 , which algebraically looks like a simple relabeling, is geometrically equivalent to performing an inversion in a sphere, of the 4D space. And the selection of a 3D subspace of the result makes this into a stereographic projection, from the unit sphere in 4D to a flat 3D space. These known facts take on a simple computational form in CGA, as follows.
To perform the desired stereographic projection from 4D to 3D, consider the dual sphere e + ; it resides in either space. Use σ = e 4 + e + (the real dual sphere in 4D with center at e 4 and radius As a consequence, orbits on the unit sphere e + in 4-space of, say, the 4D Euclidean rotation with 2-blade e 3 e 4 become mapped to circles in 3D, the orbits of a conformal rotation in 3-space with 2-blade e 3 e + . Since this no longer involves e 4 , we can depict the results in a (conformal) 3D space. The inversion in σ also maps the rest of the 4D space (i.e. points not on its unit sphere), to points outside the hyperplane e 4 . We do not depict those when interpreting the mapping as a stereographic projection. Thus only the sphere S 3 (i.e. e + in 4D) is depicted in compactified R 3 (which is the space that constitutes all points of the 4D coordinate plane e 4 , plus the point at infinity).
In an attempt to illustrate this by analogy, in Fig. 11 we show the same principle in one dimension less. We map from 3D to 2D, employing the versor e + + e 3 to swap the 2-sphere e + of 3D with the plane e 3 , which is then interpreted as a compactified 2D space.
To express and derive certain properties of Villarceau circles very compactly, we have seen that the orthonormal basis of our algebra needs to be completed by the imaginary sphere e − . The existence of this 5-sphere basis for 3D space was already explored in a 1934 paper by Raynor [11]. The imaginary sphere e − is not affected by the inversion above, and so can be considered to be part of both the 4D and the 3D spaces involved in our construction.

Villarceau Circles on Dupin Cyclides
Within CGA, it is easy to generalize from the specific unit-sphere and unitplane based generators we have employed in this paper. The conformal covariance of all expressions implies that we may apply any conformal versor V to all elements, and preserve their relationships. The corresponding 2-blades In fact, there is not much more to say about this, since obviously any property that can be formulated in a covariant manner transforms covariantly. So we merely show some Villarceau circles on cyclides in Fig. 12 that are the equivalent (by conformal transformations) of Fig. 3b for the torus.
We just have to be a bit careful transferring the property of Fig. 1: namely that the two Villarceau circles can be considered as a planar cross section of a torus. The original cut plane of the figure will transform to a sphere under the conformal versor, and it is a trivial statement that two 3D circles intersecting in a point pair lie on a common sphere. So the desired cut plane of the cyclide is not the transform of the cut plane of the torus. The proper procedure is to transfer not the plane itself, but the condition of Eq. (12): that for each Villarceau circle generated by B + + B − there exists a point x on it so that the circle generated by B + −B − from that point shares a plane with the first circle. The non-trivial solution for our standard situation was x·e + = 0: we can start any Villarceau circle and its coplanar counterpart from the same point on the unit sphere. After the conformal transformation, that condition becomes x ·(V e + /V ) = 0 for transformed points x . Points on this transformed sphere still produce two Villarceau circles in a coplanar pair. Therefore the cyclides, too, can be cut by a well-chosen plane to produce two intersecting Villarceau circles. Figure 13 illustrates this fact, for an interesting special case in which one of the Villarceau circles has become a line.

Maxwell-Dirac Toroids
As stated in the introduction, Villarceau circles occur in physics. When they are specific field lines related to the electromagnetic field in vacuum, as in [2], they are geometric objects without a need for a parametrization. Such field lines themselves move with the speed of light, and are space filling in Figure 13. An unbounded cyclide (generated from the same imaginary point pair generator as in Fig. 12a by a well-chosen spatial point). Two Villarceau circles are indicated in red-one has degenerated to a line-as well as the cutting plane that can be thought to produce them (color figure online) the manner of a Hopf fibration. The electric field E and the magnetic field B deform under this motion (or similarly, if they are considered in a different inertial frame of reference), though their Poynting field line E × B retains its canonical shape of a Villarceau circle in a Hopf fibration. Our way of describing this geometry compactly by a location-dependent circle may be of some interest in that field. It is well known that Maxwell's equation exhibits conformal symmetries (see e.g. [12]), so one would expect conformal rotor descriptions of the geometry it generates to be efficacious. The conformal Villarceau orbits from this paper are a special case of the more general conformal (m, n) knots treated in [5], which may then apply to the more general fields, also exposed in [2].
As for the toroidal motion of the electron (and the hypothesized internal structure of the photon), [9] derives the Consa parametrization of the torus directly from the Dirac equation. To do so, Hestenes employs a factorization of the spinor wave function that splits off a Lorentz boost, leaving a purely Euclidean rotor of two factors that easily give rise to the Consa toroidal solution, in the Euclidean two-angle parametrization by what we would call simple rotors with 2-blades e 1 e 2 and e 2 e 3 . That sequence of Euclidean rotations makes the Consa parametrization of the torus seem natural, though the parameter is not physical travel time-the travel speed is actually hypothesized to be uniform (equal to the speed of light). We suggest that factorizing the wave function into conformal components (in our parlance, the rotors of e 1 e 2 and e 3 e + ) may also be feasible, and in that context the parametrization from this paper may be useful. This is not just a change of representation: should the electron move in a closed (1,1) knot motion along the torus, then the geometry of the spatial orbit differs between the uniform and the conformal parametrizations. The conformal parametrization would have it moving in a plane (as we showed in this paper), but the uniform angle parametrization of the doubly-periodic motion of the Consa model is non-planar. The different accelerations those two motions impose may lead to physically observable consequences, which would force the natural choice between them. Motion along the general (m, n) knots would similarly be geometrically, and hence physically, different. The conformal description of the Robinson congruence (the guise of Villarceau circles in twistor theory) has been given in terms of conformal rotors in [1], for massless particles. Viewing the Villarceau circles of the congruence as a special case of (m, n) conformal torus knots may perhaps allow extension of that geometrical twistor theory to treat other particles as well.
Open Access. This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4. 0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.