3-sphere

Last updated
Stereographic projection of the hypersphere's parallels (red), meridians (blue) and hypermeridians (green). Because this projection is conformal, the curves intersect each other orthogonally (in the yellow points) as in 4D. All curves are circles: the curves that intersect <0,0,0,1>  have infinite radius (= straight line). In this picture, the whole 3D space maps the surface of the hypersphere, whereas in the next picture the 3D space contained the shadow of the bulk hypersphere. Hypersphere coord.PNG
Stereographic projection of the hypersphere's parallels (red), meridians (blue) and hypermeridians (green). Because this projection is conformal, the curves intersect each other orthogonally (in the yellow points) as in 4D. All curves are circles: the curves that intersect 0,0,0,1 have infinite radius (= straight line). In this picture, the whole 3D space maps the surface of the hypersphere, whereas in the next picture the 3D space contained the shadow of the bulk hypersphere.
Direct projection of 3-sphere into 3D space and covered with surface grid, showing structure as stack of 3D spheres (2-spheres) Hypersphere.png
Direct projection of 3-sphere into 3D space and covered with surface grid, showing structure as stack of 3D spheres (2-spheres)

In mathematics, a hypersphere or 3-sphere is a 4-dimensional analogue of a sphere, and is the 3-dimensional n-sphere. In 4-dimensional Euclidean space, it is the set of points equidistant from a fixed central point. The interior of a 3-sphere is a 4-ball.

Contents

It is called a 3-sphere because topologically, the surface itself is 3-dimensional, even though it is curved into the 4th dimension. For example, when traveling on a 3-sphere, you can go north and south, east and west, or along a 3rd set of cardinal directions. This means that a 3-sphere is an example of a 3-manifold.

Definition

In coordinates, a 3-sphere with center (C0, C1, C2, C3) and radius r is the set of all points (x0, x1, x2, x3) in real, 4-dimensional space (R4) such that

The 3-sphere centered at the origin with radius 1 is called the unit 3-sphere and is usually denoted S3:

It is often convenient to regard R4 as the space with 2 complex dimensions (C2) or the quaternions (H). The unit 3-sphere is then given by

or

This description as the quaternions of norm one identifies the 3-sphere with the versors in the quaternion division ring. Just as the unit circle is important for planar polar coordinates, so the 3-sphere is important in the polar view of 4-space involved in quaternion multiplication. See polar decomposition of a quaternion for details of this development of the three-sphere. This view of the 3-sphere is the basis for the study of elliptic space as developed by Georges Lemaître. [1]

Properties

Elementary properties

The 3-dimensional surface volume of a 3-sphere of radius r is

while the 4-dimensional hypervolume (the content of the 4-dimensional region, or ball, bounded by the 3-sphere) is

Every non-empty intersection of a 3-sphere with a three-dimensional hyperplane is a 2-sphere (unless the hyperplane is tangent to the 3-sphere, in which case the intersection is a single point). As a 3-sphere moves through a given three-dimensional hyperplane, the intersection starts out as a point, then becomes a growing 2-sphere that reaches its maximal size when the hyperplane cuts right through the "equator" of the 3-sphere. Then the 2-sphere shrinks again down to a single point as the 3-sphere leaves the hyperplane.

In a given three-dimensional hyperplane, a 3-sphere can rotate about an "equatorial plane" (analogous to a 2-sphere rotating about a central axis), in which case it appears to be a 2-sphere whose size is constant.

Topological properties

A 3-sphere is a compact, connected, 3-dimensional manifold without boundary. It is also simply connected. What this means, in the broad sense, is that any loop, or circular path, on the 3-sphere can be continuously shrunk to a point without leaving the 3-sphere. The Poincaré conjecture, proved in 2003 by Grigori Perelman, provides that the 3-sphere is the only three-dimensional manifold (up to homeomorphism) with these properties.

The 3-sphere is homeomorphic to the one-point compactification of R3. In general, any topological space that is homeomorphic to the 3-sphere is called a topological 3-sphere.

The homology groups of the 3-sphere are as follows: H0(S3, Z) and H3(S3, Z) are both infinite cyclic, while Hi(S3, Z) = {} for all other indices i. Any topological space with these homology groups is known as a homology 3-sphere. Initially Poincaré conjectured that all homology 3-spheres are homeomorphic to S3, but then he himself constructed a non-homeomorphic one, now known as the Poincaré homology sphere. Infinitely many homology spheres are now known to exist. For example, a Dehn filling with slope 1/n on any knot in the 3-sphere gives a homology sphere; typically these are not homeomorphic to the 3-sphere.

As to the homotopy groups, we have π1(S3) = π2(S3) = {} and π3(S3) is infinite cyclic. The higher-homotopy groups (k ≥ 4) are all finite abelian but otherwise follow no discernible pattern. For more discussion see homotopy groups of spheres.

Homotopy groups of S3
k012345678910111213141516
πk(S3)000ZZ2Z2Z12Z2Z2Z3Z15Z2Z2Z2Z12Z2Z84Z2Z2Z2Z2Z6

Geometric properties

The 3-sphere is naturally a smooth manifold, in fact, a closed embedded submanifold of R4. The Euclidean metric on R4 induces a metric on the 3-sphere giving it the structure of a Riemannian manifold. As with all spheres, the 3-sphere has constant positive sectional curvature equal to 1/r2 where r is the radius.

Much of the interesting geometry of the 3-sphere stems from the fact that the 3-sphere has a natural Lie group structure given by quaternion multiplication (see the section below on group structure). The only other spheres with such a structure are the 0-sphere and the 1-sphere (see circle group).

Unlike the 2-sphere, the 3-sphere admits nonvanishing vector fields (sections of its tangent bundle). One can even find three linearly independent and nonvanishing vector fields. These may be taken to be any left-invariant vector fields forming a basis for the Lie algebra of the 3-sphere. This implies that the 3-sphere is parallelizable. It follows that the tangent bundle of the 3-sphere is trivial. For a general discussion of the number of linear independent vector fields on a n-sphere, see the article vector fields on spheres.

There is an interesting action of the circle group T on S3 giving the 3-sphere the structure of a principal circle bundle known as the Hopf bundle. If one thinks of S3 as a subset of C2, the action is given by

.

The orbit space of this action is homeomorphic to the two-sphere S2. Since S3 is not homeomorphic to S2 × S1, the Hopf bundle is nontrivial.

Topological construction

There are several well-known constructions of the three-sphere. Here we describe gluing a pair of three-balls and then the one-point compactification.

Gluing

A 3-sphere can be constructed topologically by "gluing" together the boundaries of a pair of 3-balls. The boundary of a 3-ball is a 2-sphere, and these two 2-spheres are to be identified. That is, imagine a pair of 3-balls of the same size, then superpose them so that their 2-spherical boundaries match, and let matching pairs of points on the pair of 2-spheres be identically equivalent to each other. In analogy with the case of the 2-sphere (see below), the gluing surface is called an equatorial sphere.

Note that the interiors of the 3-balls are not glued to each other. One way to think of the fourth dimension is as a continuous real-valued function of the 3-dimensional coordinates of the 3-ball, perhaps considered to be "temperature". We take the "temperature" to be zero along the gluing 2-sphere and let one of the 3-balls be "hot" and let the other 3-ball be "cold". The "hot" 3-ball could be thought of as the "upper hemisphere" and the "cold" 3-ball could be thought of as the "lower hemisphere". The temperature is highest/lowest at the centers of the two 3-balls.

This construction is analogous to a construction of a 2-sphere, performed by gluing the boundaries of a pair of disks. A disk is a 2-ball, and the boundary of a disk is a circle (a 1-sphere). Let a pair of disks be of the same diameter. Superpose them and glue corresponding points on their boundaries. Again one may think of the third dimension as temperature. Likewise, we may inflate the 2-sphere, moving the pair of disks to become the northern and southern hemispheres.

One-point compactification

After removing a single point from the 2-sphere, what remains is homeomorphic to the Euclidean plane. In the same way, removing a single point from the 3-sphere yields three-dimensional space. An extremely useful way to see this is via stereographic projection. We first describe the lower-dimensional version.

Rest the south pole of a unit 2-sphere on the xy-plane in three-space. We map a point P of the sphere (minus the north pole N) to the plane by sending P to the intersection of the line NP with the plane. Stereographic projection of a 3-sphere (again removing the north pole) maps to three-space in the same manner. (Notice that, since stereographic projection is conformal, round spheres are sent to round spheres or to planes.)

A somewhat different way to think of the one-point compactification is via the exponential map. Returning to our picture of the unit two-sphere sitting on the Euclidean plane: Consider a geodesic in the plane, based at the origin, and map this to a geodesic in the two-sphere of the same length, based at the south pole. Under this map all points of the circle of radius π are sent to the north pole. Since the open unit disk is homeomorphic to the Euclidean plane, this is again a one-point compactification.

The exponential map for 3-sphere is similarly constructed; it may also be discussed using the fact that the 3-sphere is the Lie group of unit quaternions.

Coordinate systems on the 3-sphere

The four Euclidean coordinates for S3 are redundant since they are subject to the condition that x02 + x12 + x22 + x32 = 1. As a 3-dimensional manifold one should be able to parameterize S3 by three coordinates, just as one can parameterize the 2-sphere using two coordinates (such as latitude and longitude). Due to the nontrivial topology of S3 it is impossible to find a single set of coordinates that cover the entire space. Just as on the 2-sphere, one must use at least two coordinate charts. Some different choices of coordinates are given below.

Hyperspherical coordinates

It is convenient to have some sort of hyperspherical coordinates on S3 in analogy to the usual spherical coordinates on S2. One such choice — by no means unique — is to use (ψ, θ, φ), where

where ψ and θ run over the range 0 to π, and φ runs over 0 to 2π. Note that, for any fixed value of ψ, θ and φ parameterize a 2-sphere of radius , except for the degenerate cases, when ψ equals 0 or π, in which case they describe a point.

The round metric on the 3-sphere in these coordinates is given by [2]

and the volume form by

These coordinates have an elegant description in terms of quaternions. Any unit quaternion q can be written as a versor:

where τ is a unit imaginary quaternion; that is, a quaternion that satisfies τ2 = −1. This is the quaternionic analogue of Euler's formula. Now the unit imaginary quaternions all lie on the unit 2-sphere in Im H so any such τ can be written:

With τ in this form, the unit quaternion q is given by

where x0,1,2,3 are as above.

When q is used to describe spatial rotations (cf. quaternions and spatial rotations), it describes a rotation about τ through an angle of 2ψ.

Hopf coordinates

The Hopf fibration can be visualized using a stereographic projection of S to R and then compressing R to a ball. This image shows points on S and their corresponding fibers with the same color. Hopf Fibration.png
The Hopf fibration can be visualized using a stereographic projection of S to R and then compressing R to a ball. This image shows points on S and their corresponding fibers with the same color.

For unit radius another choice of hyperspherical coordinates, (η, ξ1, ξ2), makes use of the embedding of S3 in C2. In complex coordinates (z1, z2) ∈ C2 we write

This could also be expressed in R4 as

Here η runs over the range 0 to π/2, and ξ1 and ξ2 can take any values between 0 and 2π. These coordinates are useful in the description of the 3-sphere as the Hopf bundle

A diagram depicting the poloidal (x1) direction, represented by the red arrow, and the toroidal (x2) direction, represented by the blue arrow, although the terms poloidal and toroidal are arbitrary in this flat torus case. Toroidal coord.png
A diagram depicting the poloidal (ξ1) direction, represented by the red arrow, and the toroidal (ξ2) direction, represented by the blue arrow, although the terms poloidal and toroidal are arbitrary in this flat torus case.

For any fixed value of η between 0 and π/2, the coordinates (ξ1, ξ2) parameterize a 2-dimensional torus. Rings of constant ξ1 and ξ2 above form simple orthogonal grids on the tori. See image to right. In the degenerate cases, when η equals 0 or π/2, these coordinates describe a circle.

The round metric on the 3-sphere in these coordinates is given by

and the volume form by

To get the interlocking circles of the Hopf fibration, make a simple substitution in the equations above [3]

In this case η, and ξ1 specify which circle, and ξ2 specifies the position along each circle. One round trip (0 to 2π) of ξ1 or ξ2 equates to a round trip of the torus in the 2 respective directions.

Stereographic coordinates

Another convenient set of coordinates can be obtained via stereographic projection of S3 from a pole onto the corresponding equatorial R3 hyperplane. For example, if we project from the point (−1, 0, 0, 0) we can write a point p in S3 as

where u = (u1, u2, u3) is a vector in R3 and u2 = u12 + u22 + u32. In the second equality above, we have identified p with a unit quaternion and u = u1i + u2j + u3k with a pure quaternion. (Note that the numerator and denominator commute here even though quaternionic multiplication is generally noncommutative). The inverse of this map takes p = (x0, x1, x2, x3) in S3 to

We could just as well have projected from the point (1, 0, 0, 0), in which case the point p is given by

where v = (v1, v2, v3) is another vector in R3. The inverse of this map takes p to

Note that the u coordinates are defined everywhere but (−1, 0, 0, 0) and the v coordinates everywhere but (1, 0, 0, 0). This defines an atlas on S3 consisting of two coordinate charts or "patches", which together cover all of S3. Note that the transition function between these two charts on their overlap is given by

and vice versa.

Group structure

When considered as the set of unit quaternions, S3 inherits an important structure, namely that of quaternionic multiplication. Because the set of unit quaternions is closed under multiplication, S3 takes on the structure of a group. Moreover, since quaternionic multiplication is smooth, S3 can be regarded as a real Lie group. It is a nonabelian, compact Lie group of dimension 3. When thought of as a Lie group, S3 is often denoted Sp(1) or U(1, H).

It turns out that the only spheres that admit a Lie group structure are S1 , thought of as the set of unit complex numbers, and S3, the set of unit quaternions (The degenerate case S0 which consists of the real numbers 1 and −1 is also a Lie group, albeit a 0-dimensional one). One might think that S7, the set of unit octonions, would form a Lie group, but this fails since octonion multiplication is nonassociative. The octonionic structure does give S7 one important property: parallelizability . It turns out that the only spheres that are parallelizable are S1, S3, and S7.

By using a matrix representation of the quaternions, H, one obtains a matrix representation of S3. One convenient choice is given by the Pauli matrices:

This map gives an injective algebra homomorphism from H to the set of 2 × 2 complex matrices. It has the property that the absolute value of a quaternion q is equal to the square root of the determinant of the matrix image of q.

The set of unit quaternions is then given by matrices of the above form with unit determinant. This matrix subgroup is precisely the special unitary group SU(2). Thus, S3 as a Lie group is isomorphic to SU(2).

Using our Hopf coordinates (η, ξ1, ξ2) we can then write any element of SU(2) in the form

Another way to state this result is if we express the matrix representation of an element of SU(2) as a exponential of a linear combination of the Pauli matrices. It is seen that an arbitrary element U ∈ SU(2) can be written as

[4]

The condition that the determinant of U is +1 implies that the coefficients α1 are constrained to lie on a 3-sphere.

In literature

In Edwin Abbott Abbott's Flatland , published in 1884, and in Sphereland , a 1965 sequel to Flatland by Dionys Burger, the 3-sphere is referred to as an oversphere, and a 4-sphere is referred to as a hypersphere.

Writing in the American Journal of Physics , [5] Mark A. Peterson describes three different ways of visualizing 3-spheres and points out language in The Divine Comedy that suggests Dante viewed the Universe in the same way; Carlo Rovelli supports the same idea. [6]

In Art Meets Mathematics in the Fourth Dimension, [7] Stephen L. Lipscomb develops the concept of the hypersphere dimensions as it relates to art, architecture, and mathematics.

See also

Related Research Articles

<span class="mw-page-title-main">Laplace's equation</span> Second-order partial differential equation

In mathematics and physics, Laplace's equation is a second-order partial differential equation named after Pierre-Simon Laplace, who first studied its properties. This is often written as or where is the Laplace operator, is the divergence operator, is the gradient operator, and is a twice-differentiable real-valued function. The Laplace operator therefore maps a scalar function to another scalar function.

<span class="mw-page-title-main">Stereographic projection</span> Particular mapping that projects a sphere onto a plane

In mathematics, a stereographic projection is a perspective projection of the sphere, through a specific point on the sphere, onto a plane perpendicular to the diameter through the point. It is a smooth, bijective function from the entire sphere except the center of projection to the entire plane. It maps circles on the sphere to circles or lines on the plane, and is conformal, meaning that it preserves angles at which curves meet and thus locally approximately preserves shapes. It is neither isometric nor equiareal.

In mechanics and geometry, the 3D rotation group, often denoted SO(3), is the group of all rotations about the origin of three-dimensional Euclidean space under the operation of composition.

<span class="mw-page-title-main">Lorentz group</span> Lie group of Lorentz transformations

In physics and mathematics, the Lorentz group is the group of all Lorentz transformations of Minkowski spacetime, the classical and quantum setting for all (non-gravitational) physical phenomena. The Lorentz group is named for the Dutch physicist Hendrik Lorentz.

<span class="mw-page-title-main">Hopf fibration</span> Fiber bundle of the 3-sphere over the 2-sphere, with 1-spheres as fibers

In differential topology, the Hopf fibration describes a 3-sphere in terms of circles and an ordinary sphere. Discovered by Heinz Hopf in 1931, it is an influential early example of a fiber bundle. Technically, Hopf found a many-to-one continuous function from the 3-sphere onto the 2-sphere such that each distinct point of the 2-sphere is mapped from a distinct great circle of the 3-sphere. Thus the 3-sphere is composed of fibers, where each fiber is a circle — one for each point of the 2-sphere.

<span class="mw-page-title-main">Bloch sphere</span> Geometrical representation of the pure state space of a two-level quantum mechanical system

In quantum mechanics and computing, the Bloch sphere is a geometrical representation of the pure state space of a two-level quantum mechanical system (qubit), named after the physicist Felix Bloch.

In physics, the Hamilton–Jacobi equation, named after William Rowan Hamilton and Carl Gustav Jacob Jacobi, is an alternative formulation of classical mechanics, equivalent to other formulations such as Newton's laws of motion, Lagrangian mechanics and Hamiltonian mechanics.

Spatial rotations in three dimensions can be parametrized using both Euler angles and unit quaternions. This article explains how to convert between the two representations. Actually this simple use of "quaternions" was first presented by Euler some seventy years earlier than Hamilton to solve the problem of magic squares. For this reason the dynamics community commonly refers to quaternions in this application as "Euler parameters".

In mathematics, the group of rotations about a fixed point in four-dimensional Euclidean space is denoted SO(4). The name comes from the fact that it is the special orthogonal group of order 4.

In physics and fluid mechanics, a Blasius boundary layer describes the steady two-dimensional laminar boundary layer that forms on a semi-infinite plate which is held parallel to a constant unidirectional flow. Falkner and Skan later generalized Blasius' solution to wedge flow, i.e. flows in which the plate is not parallel to the flow.

In special functions, a topic in mathematics, spin-weighted spherical harmonics are generalizations of the standard spherical harmonics and—like the usual spherical harmonics—are functions on the sphere. Unlike ordinary spherical harmonics, the spin-weighted harmonics are U(1) gauge fields rather than scalar fields: mathematically, they take values in a complex line bundle. The spin-weighted harmonics are organized by degree l, just like ordinary spherical harmonics, but have an additional spin weights that reflects the additional U(1) symmetry. A special basis of harmonics can be derived from the Laplace spherical harmonics Ylm, and are typically denoted by sYlm, where l and m are the usual parameters familiar from the standard Laplace spherical harmonics. In this special basis, the spin-weighted spherical harmonics appear as actual functions, because the choice of a polar axis fixes the U(1) gauge ambiguity. The spin-weighted spherical harmonics can be obtained from the standard spherical harmonics by application of spin raising and lowering operators. In particular, the spin-weighted spherical harmonics of spin weight s = 0 are simply the standard spherical harmonics:

In geometry, various formalisms exist to express a rotation in three dimensions as a mathematical transformation. In physics, this concept is applied to classical mechanics where rotational kinematics is the science of quantitative description of a purely rotational motion. The orientation of an object at a given instant is described with the same tools, as it is defined as an imaginary rotation from a reference placement in space, rather than an actually observed rotation from a previous placement in space.

In mathematics, prolate spheroidal wave functions are eigenfunctions of the Laplacian in prolate spheroidal coordinates, adapted to boundary conditions on certain ellipsoids of revolution. Related are the oblate spheroidal wave functions.

In mathematics, the spectral theory of ordinary differential equations is the part of spectral theory concerned with the determination of the spectrum and eigenfunction expansion associated with a linear ordinary differential equation. In his dissertation, Hermann Weyl generalized the classical Sturm–Liouville theory on a finite closed interval to second order differential operators with singularities at the endpoints of the interval, possibly semi-infinite or infinite. Unlike the classical case, the spectrum may no longer consist of just a countable set of eigenvalues, but may also contain a continuous part. In this case the eigenfunction expansion involves an integral over the continuous part with respect to a spectral measure, given by the Titchmarsh–Kodaira formula. The theory was put in its final simplified form for singular differential equations of even degree by Kodaira and others, using von Neumann's spectral theorem. It has had important applications in quantum mechanics, operator theory and harmonic analysis on semisimple Lie groups.

In fluid dynamics, the Oseen equations describe the flow of a viscous and incompressible fluid at small Reynolds numbers, as formulated by Carl Wilhelm Oseen in 1910. Oseen flow is an improved description of these flows, as compared to Stokes flow, with the (partial) inclusion of convective acceleration.

<span class="mw-page-title-main">Cnoidal wave</span> Nonlinear and exact periodic wave solution of the Korteweg–de Vries equation

In fluid dynamics, a cnoidal wave is a nonlinear and exact periodic wave solution of the Korteweg–de Vries equation. These solutions are in terms of the Jacobi elliptic function cn, which is why they are coined cnoidal waves. They are used to describe surface gravity waves of fairly long wavelength, as compared to the water depth.

<span class="mw-page-title-main">Symmetry in quantum mechanics</span> Properties underlying modern physics

Symmetries in quantum mechanics describe features of spacetime and particles which are unchanged under some transformation, in the context of quantum mechanics, relativistic quantum mechanics and quantum field theory, and with applications in the mathematical formulation of the standard model and condensed matter physics. In general, symmetry in physics, invariance, and conservation laws, are fundamentally important constraints for formulating physical theories and models. In practice, they are powerful methods for solving problems and predicting what can happen. While conservation laws do not always give the answer to the problem directly, they form the correct constraints and the first steps to solving a multitude of problems. In application, understanding symmetries can also provide insights on the eigenstates that can be expected. For example, the existence of degenerate states can be inferred by the presence of non commuting symmetry operators or that the non degenerate states are also eigenvectors of symmetry operators.

In image analysis, the generalized structure tensor (GST) is an extension of the Cartesian structure tensor to curvilinear coordinates. It is mainly used to detect and to represent the "direction" parameters of curves, just as the Cartesian structure tensor detects and represents the direction in Cartesian coordinates. Curve families generated by pairs of locally orthogonal functions have been the best studied.

A proper reference frame in the theory of relativity is a particular form of accelerated reference frame, that is, a reference frame in which an accelerated observer can be considered as being at rest. It can describe phenomena in curved spacetime, as well as in "flat" Minkowski spacetime in which the spacetime curvature caused by the energy–momentum tensor can be disregarded. Since this article considers only flat spacetime—and uses the definition that special relativity is the theory of flat spacetime while general relativity is a theory of gravitation in terms of curved spacetime—it is consequently concerned with accelerated frames in special relativity.

Moffatt eddies are sequences of eddies that develop in corners bounded by plane walls due to an arbitrary disturbance acting at asymptotically large distances from the corner. Although the source of motion is the arbitrary disturbance at large distances, the eddies develop quite independently and thus solution of these eddies emerges from an eigenvalue problem, a self-similar solution of the second kind.

References

  1. Lemaître, Georges (1948). "Quaternions et espace elliptique". Acta. 12. Pontifical Academy of Sciences: 57–78.
  2. Landau, Lev D.; Lifshitz, Evgeny M. (1988). Classical Theory of Fields. Course of Theoretical Physics. Vol. 2 (7th ed.). Moscow: Nauka. p. 385. ISBN   978-5-02-014420-0.
  3. Banchoff, Thomas. "The Flat Torus in the Three-Sphere".
  4. Schwichtenberg, Jakob (2015). Physics from symmetry. Cham: Springer. ISBN   978-3-319-19201-7. OCLC   910917227.
  5. Peterson, Mark A. (1979). "Dante and the 3-sphere". American Journal of Physics. 47 (12): 1031–1035. Bibcode:1979AmJPh..47.1031P. doi:10.1119/1.11968. Archived from the original on 23 February 2013.
  6. Rovelli, Carlo (9 September 2021). General Relativity: The Essentials. Cambridge: Cambridge University Press. ISBN   978-1-00-901369-7 . Retrieved 13 September 2021.
  7. Lipscomb, Stephen (2014). Art meets mathematics in the fourth dimension (2 ed.). Berlin: Springer. ISBN   978-3-319-06254-9. OCLC   893872366.

Further reading