Rolling Hypocycloids

The hypocycloid with n cusps is the curve traced out by a point on a circle rolling inside a circle whose radius is n times larger.

The hypocycloid with 2 cusps is sort of strange:

It’s just a line segment! It’s called the Tusi couple.

The hypocycloid with 3 cusps is called the deltoid, because it’s shaped like the Greek letter Δ:

The hypocycloid with 4 cusps is called the astroid, because it reminds people of a star:

I got interested in hypocycloids while writing some articles here:

Rolling circles and balls.

My goal was to explain a paper John Huerta and I wrote about the truly amazing things that happen when you roll a ball on a ball 3 times as big. But I got into some interesting digressions.

While pondering hypocycloids in the ensuing discussion, the mathematician, physicist, programmer and science fiction author Greg Egan created this intriguing movie:

It’s a deltoid rolling inside an astroid. It fits in a perfectly snug way, with all three corners touching the astroid at all times!

Why does it work? It’s related to some things that physicists like: SU(3) and SU(4).

SU(n) is the group of n × n unitary matrices with determinant 1. Physicists like Pauli and Heisenberg got interested in SU(2) when they realized it describes the rotational symmetries of an electron. You need it to understand the spin of an electron. Later, Gell-Mann got interested in SU(3) because he thought it described the symmetries of quarks. He won a Nobel prize for predicting a new particle based on this theory, which was then discovered.

We now know that SU(3) does describe symmetries of quarks, but not in the way Gell-Mann thought. It turns out that quarks come in 3 ‘colors’—not real colors, but jokingly called red, green and blue. Similarly, electrons come in 2 different spin states, called up and down. Matrices in SU(3) can change the color of a quark, just as matrices in SU(2) can switch an electron’s spin from up to down, or some mix of up and down.

SU(4) would be important in physics if quarks came in 4 colors. In fact there’s a theory saying they do, with electrons and neutrinos being examples of quarks in their 4th color state! It’s called the Pati–Salam theory, and you can see an explanation here:

• John Baez and John Huerta, The algebra of grand unified theories.

It’s a lot of fun, because it unifies leptons (particles like electrons and neutrinos) and quarks. There’s even a chance that it’s true. But it’s not very popular these days, because it has some problems. It predicts that protons decay, which we haven’t seen happen yet.

Anyway, the math of SU(3) and SU(4) is perfectly well understood regardless of their applications to physics. And here’s the cool part:

If you take a matrix in SU(3) and add up its diagonal entries, you can get any number in the complex plane that lies inside a deltoid. If you take a matrix in SU(4) and add up its diagonal entries, you can get any number in the complex plane that lies inside an astroid. And using how SU(3) fits inside SU(4), you can show that a deltoid moves snugly inside an astroid!

The deltoid looks like it’s rolling, hence the title of this article. But Egan pointed out that it’s not truly ‘roll’ in the sense of mechanics—it slides a bit as it rolls.

For the details, see this article on my other blog:

Deltoid rolling on an astroid.

But the really cool part—the new thing I want to show you today—is that this pattern continues!

For example, an astroid moves snugly inside a 5-pointed shape, thanks to how SU(4) sits inside SU(5). Here’s a movie of that, again made by Greg Egan:

In general, a hypocycloid with n cusps moves snugly inside a hypocycloid with n + 1 cusps. But this implies that you can have a hypocycloid with 2 cusps moves inside one with 3 cusps moving inside one with 4 cusps, etcetera! I’ve been wanting to see this for a while, but yesterday Egan created some movies showing this:

Depending on the details, you can get different patterns:

Egan explains:

In the first movie, every n-cusped hypocycloid moves inside the enclosing (n+1)-cusped hypocycloid at the same angular velocity and in the same direction, relative to the enclosing one … which is itself rotating, except for the very largest one. So the cumulative rotational velocity is highest for the line, lower for the deltoid, lower still for the astroid, in equal steps until you hit zero for the outermost hypocycloid.

In the second movie, the magnitude of the relative angular velocity is always the same between each hypocycloid and the one that encloses it, but the direction alternates as you go up each level.

Here’s one where he went up to a hypocycloid with 10 cusps:

Note how sometimes all the cusps meet!

Yonatan Zunger is the chief architect of Google+, but before that he used to work on string theory. Reading about this stuff on Google+, he had a number of interesting ideas, like this:

I wonder, can a deltoid roll in a 5-hypocycloid? I haven’t worked through the math of just how the SU(n) → SU(n+1) embedding guarantees a fit here.

(And also, are there corresponding pictures we could draw to illustrate the embeddings of other Lie groups? This could be a lovely way to illustrate a wide range of relationships if it could be done more generally.)

Egan figured out that yes, a deltoid can move nicely in a hypocycloid with 5 cusps:

He writes:

To make this, I took the border of a deltoid in “standard position” (as per the trace theorem) to be:

\mathrm{deltoid}(t) = 2 \exp(i t) + \exp(-2 i t)

Suppose I apply the linear function:

f(s, z) = \exp(-2 i s / 5) z + 2 \exp(3 i s / 5)

to the whole deltoid. Then at the point s=t on the deltoid, we have:

f(s, \mathrm{deltoid}(s)) = 4 \exp(3 i s / 5) + \exp(-4 (3 i s / 5))

which is a point on a 5-cusped hypocycloid in “standard position”. Of course with this choice of parameters you need to take s from 0 through to 10 \pi/3, not 2 \pi, to complete a cycle.

Using the same trick, you can get a hypocyloid with n cusps to move inside one with m cusps whenever n ≤ m. As for the more general question of how different Lie groups give different ‘generalized hypocycloids’, and how fitting one Lie group in another lets you roll one generalized hypocycloid in another, the current state of the art is here:

• N. Kaiser, Mean eigenvalues for simple, simply connected, compact Lie groups.

But there is still more left to do!

On a somewhat related note, check out this fractal obtained by rolling a circle inside a circle 4 times as big that’s rolling in a circle 4 times as big… and so on: astroidae ad infinitum!

69 Responses to Rolling Hypocycloids

  1. jessemckeown says:

    Clockmakers and watchmakers used to make gear teeth using these patterns (or segments of their epicycloidal counterparts). Nowadays involute gears are prefered — they’re more tollerant of varying axle separation — but it’s easier to machine flat pinion leaves by hand, which works because a diameter line is part of a 2-cusp cycloid.
    With more horological remarks

  2. Vít Tuček says:

    So here are few obvious questions that I don’t have time to think about right now. :-(

    Is there actually an S^1 family of embeddings \mathrm{SU}(n) \to \mathrm{SU}(n+1)?

    And what about the \mathrm{Spin}(4n+2)? Can we at least find similar relation between \mathrm{tr}(\mathrm{Spin}(4n+6)) and \mathrm{tr}(\mathrm{Spin}(4n+2))?

    Regarding Zunger’s suggestion on trying out to illustrate embeddings of other Lie groups – the obvious suggestion is to use another symmetric polynomial in eigenvalues, say e_2 or e_{n-1} and look at \{ \bigl(e_1(\lambda_1,\ldots,\lambda_n), e_2(\lambda_1,\ldots,\lambda_n)\bigr) \,|\, \lambda_1, \ldots, \lambda_n \in \mathrm{Spec}(g) \text{ for some } g \in G \}

    • John Baez says:

      Vit wrote:

      Is there actually an S^1 family of embeddings \mathrm{SU}(n) \to \mathrm{SU}(n+1)?

      Usually when we talk about embedding one group in another, we mean doing so via a homomorphism. There are lots of such embeddings \mathrm{SU}(n) \to \mathrm{SU}(n+1); pick the obvious one and get more by conjugating it with elements of \mathrm{SU}(n+1). I believe (but have not proved) that this trick gives all the embeddings.

      However, since the trace is invariant under conjugation, it’s not these embeddings that explain how the set \mathrm{tr}(\mathrm{SU}(n)) rolls around in the set \mathrm{tr}(\mathrm{SU}(n+1)), as you see here:

      <div align = "center"

      The actual explanation can be on my Visual Insight blog. The eigenvalues of a matrix in \mathrm{SU}(n) are of the form

      e^{i\phi_1}, \; \dots , \; e^{i\phi_{n-1}}, \; e^{-i(\phi_1+\phi_2+...+\phi_{n-1})}

      where the angles \phi_i are arbitrary. On the other hand, the eigenvalues of any element of \mathrm{SU}(n+1) can be written as

      e^{i\theta}e^{i\phi_1}, \; \dots , \; e^{i\theta}e^{i\phi_{n-1}}, \;e^{i\theta}e^{-i(\phi_1+\phi_2+...+\phi_{n-1})}, e^{-in\theta}

      where the angles \phi_i and \theta are arbitrary. Thus we have

      \displaystyle{ \mathrm{tr}(\mathrm{SU}(n+1)) = \bigcup_{0 \le \theta \le \pi} e^{i\theta} \mathrm{tr}(\mathrm{SU}(n)) + e^{-in\theta} }

      As \theta ranges from 0 to 2\pi, this gives a filled-in hypocycloid with n cusps rolling snugly inside one with n+1 cusps!

      • Layra Idarani says:

        Perhaps a presentation that is a bit easier to understand is the maximal torii; since we’re only looking at compact groups, every element is conjugate to something in a maximal torus, and the trace is invariant under conjugation, so the image of the trace over the maximal torus is the image of the trace over SU(n). SU(n) has an (n-1)-dimensional maximal torus, and SU(n+1) has an n-dimensional maximal torus, and we can view the n-dimensional torus as the product of the (n-1)-dimensional torus and S^1. This has the helpful property of easily abandoning the group structure, so we don’t worry about the embeddings being homomorphisms or not.

        • John Baez says:

          Yes, this is a good way to think about it. Concretely one can think of the maximal torus of SU(n) as the set

          \{  (e^{i\phi_1}, \dots, e^{i \phi_{n-1}} , e^{-i(\phi_1 + \cdots + \phi_{n-1})}) \colon \phi_i \in \mathbb{R} \} \subseteq \mathbb{C}^n

          and the trace map as the linear functional

          \mathrm{tr} : \mathbb{C}^n \to \mathbb{C}

          that just sums all the components of a vector. We could study this sort of thing for other compact Lie groups, either \mathrm{Spin}(4n+2) using the trace in the complex spinor representation, or more general Lie groups using other linear functionals! I haven’t had time to try it. But it seems this approach should give us various sets

          \{ \sum_{i = 1}^n  e^{i \phi_i} \} \subseteq \mathbb{C}

          where the angles \phi_i are constrained by linear relations, probably with integer coefficients. So the fun would be to understand the sets that arise.

        • I think this is how it goes for the Spin(4n+2) case:

          Viewing the maximal torus of Spin(4n+2) as diagonal matrices, we get the maximal torus of Spin(4n+6) by taking the Kronecker product with the diagonal matrix

          D(a, b) = \mathrm{diag}(e^{ia + ib}, e^{ia-ib}, e^{-ia+ib}, e^{-ia-ib})

          So for A in Spin(4n+2) and a, b \in \mathbb{R},

          \mathrm{tr}(A\otimes D(a, b)) = \mathrm{tr}(A)(e^{ia + ib} + e^{ia-ib} + e^{-ia+ib} + e^{-ia-ib})

          Unfortunately, e^{ia+ib} + e^{ia-ib}+e^{-ia+ib}+e^{-ia-ib} is always real, wandering between 4 and -4, so we get squashing but no rotation or anything of the sort. In other words, the Spin(6) = SU(4) case gets us an astroid, and the embedding of that into Spin(10) gets us an astroid that grows and shrinks, but further Spins don’t get us any new behavior.

    • Greg Egan says:

      I wonder if another useful way of thinking about this is as a group homomorphism T: \mathrm{U}(1) \times \mathrm{SU}(n) \to \mathrm{SU}(n+1), where the image of T is the subgroup of (n,1)-block-diagonal matrices within \mathrm{SU}(n+1).

      Explicitly, T(e^{i \theta}, g) is the matrix whose n \times n block is e^{i \theta} g and whose 1\times 1 block is e^{-i n\theta}. You can then obtain all the group-homorphic embeddings of \mathrm{U}(1) \times \mathrm{SU}(n) in \mathrm{SU}(n+1) by conjugation, but they all look the same when we take the trace.

    • Greg Egan says:

      More generally, there’s a group homomorphism T_k:\mathrm{U}(1) \times \mathrm{SU}(n) \to \mathrm{SU}(n+k), where T_k(e^{i \theta}, g) is an (n,k)-block-diagonal matrix with the n\times n block given by e^{i k \theta} g and the k\times k block given by e^{-i n \theta} I_k. (Putting n=3, k=2, and with a different parameter of s=-(n+k)\theta, that’s how I got the deltoid rolling in the 5-hypocycloid.)

    • Greg Egan says:

      Ah, and we can also have a group homomorphism V_{k,n}:\mathrm{U}(1) \times \mathrm{SU}(k) \times \mathrm{SU}(n) \to \mathrm{SU}(n+k), where V_{k,n}(e^{i \theta}, h, g) is an (n,k)-block-diagonal matrix whose n \times n block is e^{i k \theta} g and whose k \times k block is e^{-i n \theta} h.

      If you chose to fix different values of h\in \mathrm{SU}(k), you would get different ways of embedding an n-hypocycloid moving inside an (n+k)-hypocycloid, depending on \mathrm{tr}(h).

      For example, this movie for n=3, k=2 shows h \in \{I_2, -I_2,i \sigma_x\}:

      • John Baez says:

        Cool! And I think your homomorphism V_{3,2} is the famous homomorphism

        \phi : \mathrm{U}(1) \times \mathrm{SU}(3) \times \mathrm{SU}(2) \to \mathrm{SU}(5)

        which is used to embed the Standard Model into \mathrm{SU}(5) grand unified theory! I guess your movie actually illustrates that case!

        If you don’t remember this stuff, you can see it discussed here. Scroll down until you hit the big equation defining \phi.

        \phi isn’t 1-1, but it comes close: it’s 6-1.

        Since this homomorphism \phi explains why quarks have the weird charges they do—if you believe in the \mathrm{SU}(5) GUT—maybe some movies like the ones you’re making could illustrate some ideas from particle physics!

      • Greg Egan says:

        That explains my sense of déjà vu when I was writing this!

      • John Baez says:

        Here’s the best I can do so far at giving a physical interpretation of this movie:

        In lattice gauge theory we chop spacetime into hypercubes and label each edge with an element g of our gauge group. This describes the gauge field. The action in lattice gauge theory is a sum over all the square faces of our hypercubes. For each square we go around and multiply the group elements g_1, g_2, g_3, g_4 labelling its edges. We take the product

        g = g_1 g_2 g_3 g_4

        This is the holonomy of the gauge field around the square. We then compute the Wilson loop function

        \mathrm{tr}(g)

        in some representation, and then subtract this from its largest possible value (the dimension of that representation). That’s what we sum over squares! It’s zero when our gauge field is flat, that is, when the holonomy around each square is 1.

        So now look at the case of the \mathrm{SU}(5) grand unified theory. There’s a homomorphism

        \begin{array}{rcl} \phi : \mathrm{U}(1) \times \mathrm{SU}(3) \times \mathrm{SU}(2) &\to& \mathrm{SU}(5) \\ (e^{i\theta}, g, h) &\mapsto & \mathrm{diag}( e^{-2i\theta}g, e^{3i \theta}h)  \end{array}

        where the three factors correspond to the hypercharge force, the strong force and the weak force. (The hypercharge force and weak force combined are another way of thinking of the electroweak force.)

        If we let the \mathrm{SU}(5) gauge field be arbitrary, the value of the Wilson loop function

        \mathrm{tr}(g)

        can be anything inside the 5-cusped hypocycloid here:

        If we set the hypercharge and weak force gauge fields equal to zero, the Wilson loop function must lie inside some particular deltoid sitting in the 5-cusped hypocycloid. This is still true if the hypercharge force is set equal to any constant value—but then we get a different deltoid.

        I hope that as we turn up the hypercharge force, the deltoid rolls around precisely as shown in this movie! I haven’t actually done the calculation: I’m posting a bit hastily here. To check, we should study the set of allowed values of

        \mathrm{tr}(e^{-2i\theta}g) + 2 e^{3 i \theta}

        where g is anything in \mathrm{SU}(3), while \theta \in \mathbb{R} takes some fixed value. We’ll get a deltoid for any choice of \theta: I think that’s easy to see. I’m hoping it rolls nicely as we vary \theta.

      • Greg Egan says:

        I’m hoping it rolls nicely as we vary \theta.

        Yes, that’s exactly in accord with the formula I used for that movie of the rolling deltoid (apart from an irrelevant factor in the phase I used, which was chosen to describe the point of contact between the curves rather than anything more physical).

      • John Baez says:

        Great, this makes for a fun interpretation of the movie!

        It would be nice if there were some use for this interpretation. In an attempt to blunder toward this, we could ask:

        • What’s special about an element of \mathrm{SU}(n) whose trace lies at a cusp of a hypocycloid? I think I know the answer—it’s related to how the center of the group \mathrm{SU}(n) is \mathbb{Z}_n. So this one counts as a puzzle rather than a research problem.

        • What’s special about the values of \theta that make a cusp of the deltoid land snugly in the cusp of the 5-cusped hypocycloid? In other words, what’s special about \theta such that for all elements

        (g, h) \in \mathrm{SU}(3) \times \mathrm{SU}(2)

        the element

        \phi(e^{i\theta}, g, h) \in \mathrm{SU}(5)

        has a trace lying in a deltoid with cusp equal to one of the cusp of the 5-cusped hypocycloid

        \{ \mathrm{tr}(x) \colon x \in \mathrm{SU}(5) \}  ?

        • Greg Egan says:

          For the first puzzle: if an element of \mathrm{SU}(n) has n equal eigenvalues (which must be nth roots of unity), its trace will achieve the maximum possible absolute value of n, and lie on the disk that circumscribes the hypocycloid, which only the cusps intersect. And the only elements with equal eigenvalues like that are nth roots of unity times the identity, i.e. elements of the form e^{2 \pi k i/n} I_n for k=0,1,...n-1.

          But I don’t see how the second puzzle can have a solution. I don’t think a choice of \theta alone can guarantee anything about \mathrm{tr}(\phi(e^{i \theta},g,h)). In the movie I made with three choices for h:

          one choice has \mathrm{tr}(h)=0, and the cusps of the corresponding deltoid never touch the cusps of the 5-hypocycloid, for any value of \theta. And even for the two different choices of h with \mathrm{tr}(h)=1, the cusps of the deltoid coincide with the cusps of the 5-hypocycloid at different values of \theta.

        • John Baez says:

          Okay, let me fix my second question! I should ask: what sets of choices of (\theta, g, h) give values of the trace \mathrm{tr}(\phi(\theta, g, h)) that range over a deltoid having a cusp that touches a cusp of the 5-epicycloid?

        • Greg Egan says:

          Oops, I should have said \mathrm{tr}(h)=\pm 2 in my previous comment, not \mathrm{tr}(h)=1. That is, the movie has h \in \{I_2, -I_2, i \sigma_x\}, and hence \mathrm{tr}(h) \in \{2,-2,0\}.

          I need to get that straight to answer the revised second question! The deltoid has cusps where g=e^{2 \pi k i/3} I_3, and the 5-hypocycloid has cusps where \phi(\theta,g,h)=e^{2 \pi m i/5} I_5. So we need \theta=\pi(k/3-m/5) in order that the first three diagonal elements are correct, and then we need h=e^{\pi (m-k) i} I_2.

          This implies h \in \{I_2, -I_2\} and \mathrm{tr}(h) \in \{2, -2\}.

          So I think the answer to the question is that you get a deltoid with one cusp coinciding with one of the 5-hypocycloid cusps by letting g range over all of \mathrm{SU}(3), while fixing \theta=\pi(k/3-m/5) and h=e^{\pi (m-k) i} I_2 for any m, k \in \mathbb{Z}.

      • Vít Tuček says:

        Bravo!

        I think one can push the visualisations even a bit further, since \mathrm{SU}(2) is isomorphic to \mathrm{Spin}(3) and hence any path of rotations could give us a nicely moving n-hypocycloid inside of (n+2)-hypocycloid via

        \mathrm{SU}(2)\times \mathrm{SU}(n) \to \mathrm{SU}(n+2)

        What would Dirac’s belt trick correspond to? :) One can also add Hopf fibration S^1 \to S^3 \to S^2 to better visualize \mathrm{SU}(2) \simeq S^3.

  3. […] Rolling Hypocycloids, par John Baez. […]

  4. Another cool post, John. BTW, an n-cusped object rolling inside an n+1-cusped object is ubiquitous in the automotive world; it’s what’s called a trochoidal oil pump.

    • John Baez says:

      Interesting!

      That led me to the Wikipedia articles on the gerotor:

      A gerotor is a positive displacement pump. The name gerotor is derived from “Generated Rotor”. A gerotor unit consists of an inner and outer rotor. The inner rotor has N teeth, and the outer rotor has N+1 teeth. The inner rotor is located off-center and both rotors rotate. The geometry of the two rotors partitions the volume between them into N different dynamically-changing volumes. During the assembly’s rotation cycle, each of these volumes changes continuously, so any given volume first increases, and then decreases. An increase creates a vacuum. This vacuum creates suction, and hence, this part of the cycle is where the intake is located. As a volume decreases compression occurs. During this compression period, fluids can be pumped, or compressed (if they are gaseous fluids).

      Gerotor pumps are generally designed using a trochoidal inner rotor and an outer rotor formed by a circle with intersecting circular arcs.

  5. arch1 says:

    Thanks Greg & John for creating and sharing these fascinating, beautiful patterns.

  6. Ayman says:

    In a way, something similar is used to drill square holes with drill bits shaped like Reuleaux triangles. Though not a hypocycloid, the Reuleaux triangle does maintain three points of contact with the outer square at all times and this could probably be extended to higher Reuleaux polygons.

    • John Baez says:

      Yes, there’s some resemblance between the deltoid inside the astroid:

      and the Reuleaux triangle, which is a famous curve of constant width:

      However, I don’t know any precise relation, except for the fact that they both have a small curve moving in a big one while keeping three points of contact, and that the small curve is not truly rolling: it’s slipping!

      A while ago Greg noticed that the deltoid is slipping, as are all the higher hypocycloids in this and the other pictures here:

      I’ve fixed this blog article to reflect that.

  7. jzal says:

    Can you guys create a rolling tetrahedral hypocycloid? Does it exist? Or how about a 3D hypocycloid with an internal 2D hypocycloid that looks like a line in a perpendicular view?

    • westy31 says:

      I can think of how to trace a Platonic polyhedron, with cusps at its vertices, and hypocycloids as edges, using a small ball inside a larger ball. But the small ball would only be allowed to trace out certain trajectories. I can’t see yet if you can trace out faces along with the edges.

      Gerard

  8. David Hartwell says:

    Fascinating discussion and movies! Of course the hypocycloid is a special case of the hypotrochoid. Another special case is the ellipse – produced by rolling a circle inside a circle. I have wondered if the inverse is possible, that a circle could be traced by rolling an ellipse inside an ellipse, but the math was beyond me to explore this question.

  9. I’m seeing generalized Wankel engines in there somewhere, too

  10. […] Evelyn Lamb describes John Baez’s Rolling hypocycloids: […]

  11. […] from which I transferred: Rolling Hypocycloids The hypocycloid with n cusps is the curve traced out by a point on a circle rolling inside a […]

  12. Ali Khezeli says:

    Did anybody notice this? Rolling a deltoid inside an astroid is just like rolling a circle of radius 3 inside a circle of radius 4. The contact point moves just like rolling a circle of radius 1 inside a circle of radius 4:

    Same is true for rolling a deltoid inside a hypocycloid with 5 cusps:

    • Ali Khezeli says:

      Similarly, any hypocycloid can roll nicely on a cycloid. The contact point traces the cycloid in an standard way. Do you think it comes from a group embedding?

    • John Baez says:

      These are very nice, Ali! I think these two are implicit in work Egan has already done, but they’re still nice to see:

      The ones involving a cycloid seem more surprising to me!

    • Greg Egan says:

      These are great!

      I never thought about a hypocycloid rolling on a cycloid. In retrospect it’s a really simple idea: if you have a circle of radius 1 and a circle of radius n rolling on a line, with the x-coordinates of their centres staying the same, the small circle must be truly rolling inside the big one, because the linear velocities of points fixed to the perimeter of either circle are both zero at the point of contact.

      A point fixed to the small circle as it rolls must trace out a cycloid in the plane, but in coordinates locked to the larger circle, the same point must trace out a hypocycloid … which in turn is rotating and advancing with the large circle.

      • Ali Khezeli says:

        A similar thing must be true for epicycloids, hypocycloids and the cycloid rolling on each other.

        • westy31 says:

          Yes, I was thinking the same. Here is a picture:


          You and Greg should be able to make these pictures by just changing a couple of minus signs in the hypocycloid case.

          Gerard

        • westy31 says:

          My picture is actually an epicycloid in an epicycloid, so not quite the same. I can think of cool variations, I’ll make some soon.

          Gerard

        • westy31 says:

          More pictures:




          (Click to enlarge.)

        • John Baez says:

          Thanks, these are nice!

          (I turned your links into pictures visible here. When I do this, I usually make it so people can click on them and go to the original pictures, which are often larger.)

          What is the precise mathematical description of what’s going on here?


          Are the shapes rolling inside the hypocycloid a bunch of epicycloids?

      • westy31 says:

        To explain the maths behind the pictures:
        If you look at the generating formula’s for the epicycloid and the hypocycloid, you can notice that they are in fact the same, if you allow negative numbers for the parameters. In the Wikipedia formula, fill in a negative ‘r’ to go from ‘epi’ to ‘hypo’.
        If you take the pictures of Greg and Ali and and fill in a negative r, some interesting features continue to work. An n-cycloid fits in exactly in an (n+1) cycloid. (To make cycloids roll on the *outside*, fill in negative ‘R’ )

        If you want to experiment with this, you can download an Excel sheet:
        http://westy31.home.xs4all.nl/ScrapBook/Cycloid.xls

        I might update this sheet, tidiness is unfortunately not one of my better skills,,,
        Gerard

    • Greg Egan says:

      Dorost! Here’s a deltoid rolling on a cardioid:



      • Ali Khezeli says:

        زیباست!
        And similarly, for any two curves rolling on a third curve, (instead of three circles or lines), we get two moving loci which are always tangent.
        A calculus puzzle: When do the two shapes really “roll” on each other?

      • Greg Egan says:

        A calculus puzzle: When do the two shapes really “roll” on each other?

        If you mean, for what choice of two curves rolling on a third (like the small and large moving circles rolling on the fixed small circle in the movie above) will the moving locus be truly rolling on the fixed locus (as if the blue deltoid were truly rolling on the red cardioid) throughout the entire motion … I think the answer is that no proper choice of curves can give that result.

        But if you mean, when during the motion can a curve like the blue deltoid really be rolling on the red cardioid, it does happen for an instant in the movie above: each cusp of the deltoid has a velocity of zero at the instant it touches the cardioid.

        More generally, if the fixed curve is given by (F_x(s), F_y(s)) and the second curve rolling on it is given by (H_x(s), H_y(s)) at time zero (like the large circle in the movie above), and these functions are both parameterised by arc length s, my calculations give a formula for the magnitude of the velocity of a point on the moving locus, at the point of contact, which contains a factor of:

        H_x'(t) F_y''(t) - F_x'(t) H_y''(t)

        Whenever that factor is zero, true rolling will occur. But I don’t believe it’s possible to choose (F_x(s), F_y(s)) and (H_x(s), H_y(s)) to make this term identically zero throughout the motion, while keeping the two curves distinct. This is probably easier to see by making use of the fact that the magnitude of the tangent vector to either curve is 1, because of the choice of parameterisation, and so the expression above can be rewritten as:

        \pm \sqrt{1-H_y'(t)^2} F_y''(t) \pm \sqrt{1-F_y'(t)^2} H_y''(t)

        Having said all that … it would be much more interesting if I’m wrong and there is a counterexample that demonstrates true rolling throughout the motion!

        • Ali Khezeli says:

          I arrived to almost the same thing. Let the fixed curve be \gamma (s), the second curve (like the big circle above) be \beta (s) (at time zero) and the first curve be \alpha (s) (at time zero). The rolling happens in two cases: at the cusps and when \beta and \gamma have the same curvature at the point of contact. This shows that the shapes cannot roll in a time interval, except when \beta and \gamma have some identical arcs which match together (not an interesting case).
          In your calculations, if we substitute F''_y(t) by \kappa(t) F'_x(t), we arrive to this conclusion.
          I used another method which arrived to the physical interpretation of rolling. Let A(t) be the moving point and B(t) be the common contact point. Let \theta(t) be the angular speed of rotation of the whole curve \alpha at time t and \varphi(t) be the same for \beta. The magnitude of A'(t) is just \theta'(t)\left|A(t)-B(t)\right|. Similarly, the magnitude of the velocity of the moving point on the moving locus is (\theta'(t)-\varphi'(t))\left|A(t)-B(t)\right|. So either we should have A(t)=B(t) or \varphi'(t)=0. The former case corresponds to cusps and the latter case happens when the curvature of \beta and \gamma are equal at time t.

        • Greg Egan says:

          That’s a nice, clear account! It’s much better to rewrite the second derivative in terms of the curvature.

          Here’s the result of two equal-sized circles rolling on an ellipse whose minimum radius of curvature is equal to that of the circles. There are 4 points in each cycle when the moving locus is truly rolling on the fixed locus: three involving cusps, and one that is due solely to the matching radii of curvature.

          http://www.gregegan.net/images/circlesOnEllipse.gif

        • Ali Khezeli says:

          As a result of my above comment, the length of an arc of a n-hypocycloid is (1-\frac 1n) time the length of an arc of a cycloid.

      • jzal says:

        I think this is really pretty. So this is n+1 cusps rolling on n cusps.

        It is interesting that the deltoid has to roll around three times to get back to where it was. Can the Tusi couple line roll around the cardioid or vice versa?

    • Greg Egan says:

      And here’s a chorus line of cardioids rolling on a cycloid:



    • Greg Egan says:

      And here’s all three kinds of curve combined: a cardioid and a deltoid rolling together on a cycloid:


  13. voldrani says:

    These animations are really awesome and fun to look at! My two-year old couldn’t get enough. Would you guys be willing to release these under a Wikimedia Commons-compatible license, and letting me upload them there? CC0 (public domain) would be best, but CC-BY (attribution) or CC-BY-SA (attribution, share-alike) would also work.

  14. David Tanzer says:

    What kind of software tools do people use here for generating these nice animations?

    • Greg Egan says:

      I use Mathematica, which is a (fairly expensive) general purpose program for performing symbolic mathematics.

      I’d be surprised if there aren’t lots of much cheaper programs that can do these kinds of simple animations, but having bought Mathematica for other reasons, I haven’t gone looking for them myself.

    • Ali Khezeli says:

      I used the open source software Geogebra. The deltoid can be drawn using the Curve command:
      Curve[2cos(t) + cos(-2t), 2sin(t) + sin(-2t), t, 0, 2 pi]
      You can rotate and animate it nicely. However, for exporting the animations I had to use screen capture and convert it to gif with other software. It also exports to applets.

    • westy31 says:

      I make pictures with Excel. This is mainly because I use it also for my job as an engineer. You can export single pictures to GIF, then you can combine them into an animated GIF using for example the free program ‘Unfreeze’.

      Gerard

  15. ron davison says:

    Perhaps the cycloid should be called Heartoid?
    :)
    Add a bearing and a clutch to one of the circles to allow sliding.
    Engage clutch at points of zero sliding.

    removes sliding-friction-heat-failure problem

    All power can be transferred in a duty cycle, coasting while sliding.
    In reality this is only how it transfers power anyway in these type designs. If you do surface friction/stiction leads to galling of surface when large forces are attempted to be transferred during sliding at oblique angles.
    very cool images. thanks
    Ron Davison

  16. John Baez says:

    Over on G+, Malin Christersson writes:

    Rolling on the outside

    (follow up from the posts http://goo.gl/iYYd6b and http://goo.gl/Z7eJL7)

    If you can roll epicycloids/hypocycloids inside each other, then you can also roll them on the outside.

    I was a bit puzzled to see that the cusps of a (n+1)-hypocycloid, rolling on the outside of a n-hypocycloid, seemed to trace out a n-epicycloid. But then I realised that it is probably obvious that they do. If you roll a n-epicycloid on the inside of a (n+1)-epicycloid, the cusps of the outer epicycloid lie on the inner epicycloid. So if you instead roll epicycloids on the outside, the cusps of the (n+1)-epicycloid will glide along the inner n-epicycloid. Since the cusps of (n+1)-epicycloids are the same points as the cusps of (n+1)-hypocycloids, it all makes sense.

    The GeoGebra worksheet on GeoGebraTube:
    http://www.geogebratube.org/student/m104829

    Make sure to check out her images at http://goo.gl/iYYd6b and http://goo.gl/Z7eJL7, too! They’re versions of images here, but they’re very well-done—made using Geogrebra.

    • Todd Trimble says:

      These are just exquisite. There’s something so mesmerizing and meditative about these images (I’m listening to ambient music as I observe their unhurried motion). Bravo, Malin!

  17. westy31 says:

    I have a new picture that has an interesting feature; It is an n-hypocycloid moving in an (n-1) epicycloid. It has the property that all cusps (both inner and outer) touch the other curve:



    Gerard

You can use HTML in your comments. You can also use LaTeX, like this: $latex E = m c^2 $. The word 'latex' comes right after the first dollar sign, with a space after it.

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s

Follow

Get every new post delivered to your Inbox.

Join 2,797 other followers