Applications of Differential Geometry Linking Topological Bifurcations to Chaotic Flow Fields

: At every point p on a smooth n -manifold M there exist ( n + 1 ) skew-symmetric tensor spaces spanning differential r -forms ω with r = 0,1, · · · , n . Because d ◦ d is always zero where d is the exterior differential, it follows that every exact r -form (i.e., ω = d λ where λ is an ( r − 1 ) -form) is closed (i.e., d ω = 0) but not every closed r -form is exact. This implies the existence of a third type of differential r -form that is closed but not exact. Such forms are called harmonic forms. Every smooth n -manifold has an underlying topological structure. Many different possible topological structures exist. What distinguishes one topological structure from another is the number of holes of various dimensions it possesses. De Rham’s theory of differential forms relates the presence of r -dimensional holes in the underlying topology of a smooth n -manifold M to the presence of harmonic r -form fields on the smooth manifold. A large amount of theory is required to understand de Rham’s theorem. In this paper we summarize the differential geometry that links holes in the underlying topology of a smooth manifold with harmonic fields on the manifold. We explore the application of de Rham’s theory to (i) visual, (ii) mechanical, (iii) electrical and (iv) fluid flow systems. In particular, we consider harmonic flow fields in the intracellular aqueous solution of biological cells and we propose, on mathematical grounds, a possible role of harmonic flow fields in the folding of protein polypeptide chains.


Introduction
A well-known quote of Einstein says that as far as the laws of mathematics refer to reality, they are not certain, and as far as they are certain, they do not refer to reality.While this may be seen as controversial, it certainly is the case that mathematical theories of physical processes require interpretation in terms of experimental observations.Differential geometry is a logical structure of provable propositions, theorems and lemmas whose truth derives from axioms, definitions and the elementary connectives of formal logic.To apply differential geometry to physical systems requires the construction of bridges between events and variables in the physical world (preferably, but not necessarily, measurable) and well-defined notions in the deductively logical structure of differential geometry.These bridges are intuitive hypotheses linking reality to mathematics and their value derives from the quality of the intuition of the practitioner.If the hypothesis is true, then certain laws of mathematics can be applied to the real-world system.But even if it is demonstrated experimentally that the mathematical laws do describe the real-world behavior, this does not prove the hypothesis.Other hypotheses may do equally well.Indeed, it is not possible to prove if-then conjectures, although they can be easily proven false.In applied mathematics, the 'if-then bridges' between reality and mathematics remain forever intuitive.With respect to certainty, it is the mathematics that is primary.
That said and understood, there is value in showing that the same mathematics can apply usefully to a range of real-world systems.
De Rham's cohomology theorem of differential forms is a remarkable piece of mathematics that links together areas of mathematics that were previously thought unrelated.We believe that the mathematical facts provided by de Rham's theorem have great relevance to the nonlinear behavior and the presence of critical points and singularities in the state spaces of many different real-world systems.In this paper, we consider an interpretation of the theorem in terms of experimental observations of the dynamics of fluid flow.In essence, de Rham's cohomology theory relates the existence of holes of various dimensions in the underlying topology of smooth and Riemannian manifolds with a third type of differential form flow field on the manifold, one that is closed but not exact.Such a flow field is called an harmonic field.In other words, our aim is to demonstrate the links between the harmonic fields of mathematics and physical flows of real-world nonlinear systems with singularities and critical points observed in nature.
To illustrate the sort of nonlinear dynamical processes in the world that can be understood through this mathematical linking, consider the bubbles and fluid flow associated with pouring a glass of champagne, or the generation and flow of foam associated with an ocean wave breaking on a beach, or the chaotic turbulent flow of a mountain stream around rocks.Even the metamorphosis of a caterpillar into a moth may be within the scope of the theory.Certain singularities in the nonlinear dynamics of visual, mechanical, electrical, and fluid-flow systems can be related to harmonic fields.In particular, we will examine the possible role of vortices and eddies in harmonic flow fields induced by the presence of organelles and other molecules immersed in the intracellular aqueous solution of biological cells on the kinetics of protein folding.
Our search for this linking was inspired by two seemingly unrelated happenings.Firstly, in our analysis of the geometry of 3D binocular visual space [1][2][3], we came to realize that three different 3D geometries have to be taken into account.There is of course the 3D Euclidean geometry of the outside world.But the geometry of the 3D visual space encoded by sensory receptors is warped (or curved) by the fact that the size of images on the retinas of the eyes vary in inverse proportion to the Euclidean distance between the object in the environment and the nodal point of the eye [4].This intrinsic 3D space corresponds to a curved Riemannian manifold with a metric proportional to the inverse of the Euclidean distance squared [1].This gives rise to a singularity at the origin where the metric goes to infinity.The consciously perceived 3D visual space is different again because many different cognitive mechanisms of depth perception are brought into play.The singularity at the origin is taken care of by the existence of a hole about the egocenter in visual space.Existence of such a hole is confirmed by the fact that it is not possible to see your own head, let alone your own egocenter.This raises the question, what effect does such a hole in visual space have on the perceived flow of visual images?
The second event concerned an unusual cluster of breast cancer diagnoses that occurred in women working together at an Australian television studio [5,6].It was speculated that the cancers might be caused by leakage of electromagnetic radiation from the TV transmitter.Experts were brought in to test this hypothesis.They found no evidence of any problematic radiation levels in the women's workplace.In physical theories of electromagnetism, there are at least two types of electric field.One is given by voltage gradients and the curl of the other is given by changing magnetic fields.These two types of field are independent because the curl of a gradient is always zero.The question arises, is there a mysterious third type of electric field that would have gone undetected by the standard measurement techniques?Differential geometry suggests that there is.According to de Rham's theorem, the curl of an electric field can be zero without it being the gradient of a voltage.The question becomes, does such a third type of field exist in nature and if so, what properties does it have?
Applying differential geometry to physical systems requires a thorough knowledge of the mathematical theory, because the sequence in which it is applied does not necessarily follow the beautiful logical sequence in which the theorems were developed and proven.Application often requires jumping from one point in the mathematical theory to a distant point with a large amount of mathematics linking the two points.Moreover, an understanding of de Rham's theorem is not trivial.It requires familiarity with a large amount of differential geometry.Consequently, in Appendix A we provide detail of all the steps in our proposed mathematical linking between the appearance and disappearance of holes in the underlying topology of smooth and Riemannian manifolds and the appearance and disappearance of harmonic flow fields on the manifold.Meanwhile, in the following section we provide the reader with a more conceptual introduction to the differential geometry used in relating topological bifurcations to chaotic flow fields in nature.

Topological Structure
An important property of topological spaces is that they are not equipped with a metric.They have no measure of size.Size and topology can be thought of as independent properties.A topological space, like the set of all points on a 2-sphere, can be expanded to the size of a planet or contracted to the size of an atom without changing its topological structure.So long as it is not cut, torn, punctured or any parts of it glued together, it can be deformed in any way by any amount without changing its topological structure.For example, the 2-sphere can be stretched into a long thin rope and tied into a knot without changing its topology.
All the different shapes a topological space can be deformed into without cutting, tearing, puncturing or gluing can be mapped in a continuous, one-to-one, onto and continuously invertible fashion to each other.This is a special type of map called a homeomorphism.Any two homeomorphic spaces have all the same topological properties.Indeed, a topological property can be defined as one that is preserved by homeomorphic maps.Transformations that leave the topological structure unchanged are called homotopies, a well-known example being the transformation of a coffee mug with a handle into a donut and vice versa [7].The mathematical mapping ([8], p. 348]) is shown in Figure 1.
AppliedMath 2024, 4, FOR PEER REVIEW 3 gradient of a voltage.The question becomes, does such a third type of field exist in nature and if so, what properties does it have?Applying differential geometry to physical systems requires a thorough knowledge of the mathematical theory, because the sequence in which it is applied does not necessarily follow the beautiful logical sequence in which the theorems were developed and proven.Application often requires jumping from one point in the mathematical theory to a distant point with a large amount of mathematics linking the two points.Moreover, an understanding of de Rham's theorem is not trivial.It requires familiarity with a large amount of differential geometry.Consequently, in Appendix A we provide detail of all the steps in our proposed mathematical linking between the appearance and disappearance of holes in the underlying topology of smooth and Riemannian manifolds and the appearance and disappearance of harmonic flow fields on the manifold.Meanwhile, in the following section we provide the reader with a more conceptual introduction to the differential geometry used in relating topological bifurcations to chaotic flow fields in nature.

Topological Structure
An important property of topological spaces is that they are not equipped with a metric.They have no measure of size.Size and topology can be thought of as independent properties.A topological space, like the set of all points on a 2-sphere, can be expanded to the size of a planet or contracted to the size of an atom without changing its topological structure.So long as it is not cut, torn, punctured or any parts of it glued together, it can be deformed in any way by any amount without changing its topological structure.For example, the 2-sphere can be stretched into a long thin rope and tied into a knot without changing its topology.
All the different shapes a topological space can be deformed into without cutting, tearing, puncturing or gluing can be mapped in a continuous, one-to-one, onto and continuously invertible fashion to each other.This is a special type of map called a homeomorphism.Any two homeomorphic spaces have all the same topological properties.Indeed, a topological property can be defined as one that is preserved by homeomorphic maps.Transformations that leave the topological structure unchanged are called homotopies, a well-known example being the transformation of a coffee mug with a handle into a donut and vice versa [7].The mathematical mapping ( [8], p. 348]) is shown in Figure 1.

Figure 1.
A diagram illustrating homotopy.The 2-torus (surface of a donut) in the topological space Y and the coffee mug with handle in the topological space Y are equivalent topological structures.These are not only homeomorphically related (i.e., continuous, one-to-one, onto and invertible) but there are continuous homotopic maps f 0 , f 1 : X → Y that map the topological space X to both the 2-torus and the mug.The fact that the donut can be continuously transformed into the mug implies the existence of all the intermediate shapes in the transformation.The hole in the donut morphs into the hole between the mug and its handle.This is illustrated in the diagram by the special case in which Y = X × I, f 0 = i 0 , f 1 = i 1 , where the interval I varies from 0 to 1.For each point x ∈ X, the map f 0 (x) = i 0 (x) = (x, 0) corresponds to the 2-torus, the map f 1 (x) = i 1 (x) = (x, 1) corresponds to the mug and all the other points along the interval I correspond to the intermediate shapes.This is represented in the diagram by the homotopy H.More generally, while the topological spaces X and Y can include holes of various dimensions, they must have the same topological structure in order for the homotopy H to exist.In other words, for f 0 , f 1 : X → Y to be homotopic maps, the induced homomorphisms ( f 0 ) * , ( f 1 ) * : H r (X) → H r (Y) for all r ≥ 0 between the homology groups for X and Y must be equal, i.e., ( There is no way, however, that a 2-sphere can be transformed into a 2-torus without puncturing holes into the sphere.The 2-sphere and the 2-torus do not have the same topological structure.It was shown by A.A. Markov in 1958 (see [8]) that there is no algorithm that can determine all the different topological structures a topological space with dimension 4 or greater can have.The best we can say is that such spaces have a large number of possible different topological structures.Basically, what makes topological structures different is the different number of holes of various dimensions that exist in the topological space (e.g., the difference between the topology of a sphere and the topology of a torus is attributable to the hole(s) in the torus).

Topological Bifurcations
Based on Smale [9], Abrahams and Marsden [10] provide a description of bifurcations in the underlying topology of a smooth manifold.A smooth map f : M → N between smooth manifolds M and N is locally trivial at y 0 ∈ N if there is a neighborhood U of y 0 ∈ N such that f −1 (y) is a smooth properly embedded submanifold (i.e., a level set) in M for each y ∈ U in N and there are diffeomorphic maps h| f −1 (y) : f −1 (y) → f −1 (y 0 ) for each y ∈ U (i.e., h projects nearby level sets f −1 (y) in M diffeomorphically onto the level set f −1 (y 0 ) in M).
The bifurcation set of f : M → N is Σ f := y 0 ∈ N f fails to be locally trivial at y 0 .In other words, as y passes through Σ f , f −1 (y) changes its topological structure.Let σ( f ) denote the closed set of critical points x of f in M, i.e., By Sard's theorem, the singular values f (σ( f )) in N have measure zero (i.e., they are isolated points with empty interior and their union is nowhere dense in N).The singular values of f (σ( f )) are the bifurcation values of f in N [10].In other words, the level sets in M of singular bifurcation values in N have a topological structure that differs from the topological structures of neighboring level sets in M. The number of holes in the level set in M of a bifurcation value in N differs from the number of holes in the neighboring level sets.This can be related to isotropy subgroups, the slice theorem, the tube theorem, the (H)-orbit type submanifold, the (H)-isotropy type submanifold and the (H)-fixed point type submanifold described by Ortega and Ratiu [11] that extend the notion of critical points in submanifolds.

De Rham's Theory
We provide here only a brief description.A more expanded summary of differential geometry linking bifurcations with harmonic flow fields is given in Appendix A.
Every exact differential r-form ω with r = 0, 1, • • • , n on a smooth n-dimensional manifold M or an n-dimensional Riemannian manifold (M, g) is closed (i.e., dω = 0) because d • d is always zero.But a closed differential r-form is not necessarily exact.Differential r-forms that are closed but not exact give rise to a third type of differential r-form.There are thus three types of differential r-forms, (i) r-forms ω whose exterior derivative is not zero (i.e., dω ̸ = 0), (ii) closed r-forms that are exact differentials (i.e., ω = dλ where λ is an (r − 1)-form) and (iii) r-forms ω that are closed but not exact.This third type of r-form is called an harmonic r-form and is contained in the r-th harmonic space H r (M) of the manifold with r = 0, 1, • • • , n.It is also an element in the r-th de Rham cohomology group H r dR (M) with r = 0, 1, • • • , n. H r dR (M) is isomorphic with H r (M).Both H r (M) and H r dR (M) are vector subspaces as well as being groups.If ω is a closed r-form (i.e., dω = 0) on a smooth manifold M, and z r is an r-dimensional cycle in the underlying topology X of the manifold, then the integral of ω over z r is called a period of ω (the term period derives from the periods of elliptic integrals).From Stokes' theorem, if the r-cycle z r happens to be a boundary b r of an (r + 1)-chain c r+1 (i.e., b r = ∂c r+1 ), then the period vanishes since (See Appendix A.1 for definitions of chains, cycles and boundaries).
From this we obtain de Rham's first theorem: a closed form is globally exact if and only if all of its periods vanish [12].Knowing which forms are globally exact and which are not has important consequences.For example, (i) a smooth 1-form is conservative (i.e., its line integral is path-independent) if and only if it is exact, (ii) the integral of every globally exact form over a compact oriented smooth manifold without boundary is zero, (iii) the integral of globally closed forms over boundaries is zero and (iv) if the integral of a closed r-form over an r-dimensional submanifold without boundary is not zero, then the r-form is not exact and the submanifold is not the boundary of a compact oriented submanifold with boundary.
A remarkable property of harmonic spaces and de Rham cohomology groups is that they are homotopically and topologically invariant.They are isomorphic to singular cohomology groups H r (M), and these in turn are isomorphic to singular homology groups H r (M) on the underlying topological structure of the smooth manifold M (see Appendix A.12).This is depicted in Figure 2.This is indeed a remarkable property, because it is not immediately obvious that closeness and exactness of differential forms on a smooth manifold should have anything to do with the underlying topological structure of the manifold.But de Rham's theory shows mathematically that they do!The r-th homology group H r (M) is isomorphic to the r-th de Rham cohomology group H r dR (M) with r = 0, 1, • • • , n. H r (M) detects the presence of r-dimensional holes in the underlying topology of the smooth manifold M. The dimension of H r (M) equals the number of r-dimensional holes in the topology, and this is known as the r-th Betti number b r of the manifold.Thus, de Rham's theory of differential forms provides the link we are looking for between holes in the topological structure of a smooth or Riemannian manifold and harmonic flow fields on the manifold.manifold with boundary.
A remarkable property of harmonic spaces and de Rham cohomology groups is that they are homotopically and topologically invariant.They are isomorphic to singular cohomology groups  (), and these in turn are isomorphic to singular homology groups  () on the underlying topological structure of the smooth manifold  (see Appendix A.12).This is depicted in Figure 2. In other words, if the maps F and G are diffeomorphic, they have the same homeomorphic topological structure and so they are homotopy-invariant.This means that homotopy-equivalent manifolds have isomorphic de Rham groups.The labels Ω r−1 (M), Ω r (M) represent the skew-symmetric tensor spaces of degree (r − 1) and degree (r) for differential forms at each point on the smooth manifold M. The labels Ω r−1 (N), Ω r (N), Ω r+1 ( N) represent the skew-symmetric tensor spaces of degree (r − 1), degree (r) and degree (r + 1) for differential forms at each point on the smooth manifold N. The maps h : Ω r+1 ( N) → Ω r (M) for all values of r are called homotopy operators.Examination of the diagram shows that if the r-form ω is closed, i.e., dω = 0, then d(h ω) = G * ω − F * ω.In other words, G * ω and F * ω differ by an exact differential d(h ω).The maps F * , G * : H r dR (N) → H r dR (M) are said to be cohomologous.This is exactly how equal cohomology groups are defined.Therefore, cohomologous groups have isomorphic de Rham cohomology groups.

The Hodge Laplacian Operator
By definition, the Hodge Laplacian operator on a Riemannian manifold is ∆ : where * is the Hodge star operator * : Ω r → Ω n−r and Ω r is the skew-symmetric tensor space of r-forms for r = 0, 1, • • • , n at every point p on an n-dimensional smooth manifold M or a Riemannian manifold (M, g).The Laplacian operator ∆ plays an important role in de Rham's theory of differential forms because it is the only operator that can be used to compute harmonic forms on a Riemannian manifold.The harmonic space H r (M) = ker(∆Ω r ) for r = 0, 1, • • • , n. (Keep in mind that harmonic r-forms for all values of r on a smooth manifold or a Riemannian manifold are associated with the presence of r-dimensional holes in the underlying topology of the manifold).

Harmonic 0-Forms and 1-Forms
Harmonic 0-forms are given by the equation ∆ φ = d * d φ = div grad φ, where φ is a smooth real-valued function on the smooth manifold M or Riemannian manifold (M, g).Potential fields in a fluid are irrotational and incompressible (i.e., divergence-free) harmonic 0-forms.We focus here on computing harmonic 1-forms on Riemannian manifolds because of their importance in the nonlinear dynamics of the physical systems to be discussed below.Let λ be a differential 1-form in Ω 1 = T * p M at each point p ∈ (M, g).First, we compute the Hodge Laplacian of a 1-form λ on a Riemannian manifold (M, g), i.e., ∆λ := dd * λ + d * dλ, and then we compute the harmonic 1-form λ by setting the Laplacian ∆λ = 0.
The Ricci curvature tensor can be described geometrically as follows.Let (M, g) be a Riemannian manifold and p a point belonging to (M, g) (i.e., p ∈ (M, g)).For every unit vector ξ ∈ T p M, Ric(ξ, ξ) is the sum of the sectional curvatures of the 2-planes in Following Lang [13] (pp.423-424), the inner product of the Ricci curvature tensor with the vector ξ (i.e., ⟨Ric(λ), ξ⟩) is defined to be the scalar-valued form such that, with respect to an orthonormal frame ξ 1 , • • • , ξ n , and any unit vector field ξ with ξ 1 = ξ, we have where we denote by ⟨λ, ξ⟩ the value of the 1-form λ on the vector field ξ.As an operator on 1-forms, the Laplacian ∆ : To compute the Hodge Laplacian inner product ⟨∆λ, ξ⟩ = ⟨(dd * + d * d)λ, ξ⟩, we compute (i) ⟨dd * λ, ξ⟩ and (ii) ⟨d * dλ, ξ⟩ separately and then add: , where λ ∨ is the dual of the 1-form field λ and ξ i for i = 1, • • • , n is a g-orthonormal frame of vector fields.These orthonormal vectors have to be computed using the Riemannian metric Gram-Schmidt orthogonalizing algorithm at each point.We see that (ii) To obtain d * dλ, we first note that if ω is an r-form, we can use the proposition where λ is a 1-form and ξ and η are vectors.Using the fact that for any r-form ω where where From Equations ( 3) and (2) with η = ξ i , we obtain [13] (p. 424) Finally, we add Equations ( 1) and ( 4) to obtain Thus, for a 1-form on an n-dimensional Riemannian manifold (M, g), the Hodge Laplacian operator has the form ∆ = (− ∑ n i=1 ∇ 2 ξ i ) − Ric, where Ric is the Ricci curvature operator.The harmonic 1-form is then obtained by setting the Laplacian ∆ λ equal to zero.

Physical Examples
There exists an endless variety of behaviors in nature of nonlinear dynamical systems associated with different singularities, critical points and degeneracies in their differential equations.Most of these behaviors have not yet been described, let alone solved.One family of such behaviors considered in this paper is associated with the existence of holes of various dimensions in the underlying topology of the state space of the system.Such holes create differential forms and vector fields that are closed but not exact and create 'no-go' spaces in the dynamical state of the system.Such holes can dramatically change the dynamical behavior of the system.Isotropy subgroups in the action of a Lie group on a dynamical system, for example, create dynamic states that a vector field (or subset of vector fields) can flow into but not out of; a type of equilibrium point dynamical-state trap.But nonlinear dynamical systems can transition to neighboring vector fields and flow around the trap, creating novel nonlinear behaviors [11].In what follows, we consider the consequences on dynamical behavior of singularities in the differential equations that give rise to differential forms that are closed but not exact in visual, mechanical, electrical and fluid-flow physical systems.

The Human Visual System
Human binocular 3D visual space encoded by visual, proprioceptive and vestibular sensory receptors is warped (or curved) because the size of 2D images on the retina of each eye varies in inverse proportion to the Euclidean distance between the nodal point of the eye and the object in the environment.Geometrically speaking, 3D visual space corresponds to an egocentric negatively curved Riemannian manifold with a 3D ball of points missing from about the origin (i.e., egocenter) [1,2].This intrinsic hole corresponds to the fact that it is not possible to see your own head.What influence does such a hole about the origin in 3D warped visual space have on the perceived flow of visual images?Differential geometry and de Rham's theorem provide a possible answer to this question.(Note that other holes in 3D visual space are possible due to neural pathology.Known as scotomas, these are not considered here.) Given a 2D hole about the origin (i.e., a missing 3D ball of points about the origin enclosed within a 2-cycle), de Rham's theorem implies that closed differential r-forms for r = 0, 1, 2, 3 are not globally exact.The unit volume form Ω = dx 1 ∧ dx 2 ∧ dx 3 is a closed 3-form on the warped 3D visual space, and therefore the hole about the origin prevents this volume 3-form from being globally exact.This is just as well, because if there was no hole about the origin and the closed volume form Ω was globally exact, then all the periods c 3 Ω equal to the volume form integrated over all 3-chains (i.e., volumes) enclosed by 2-cycles would vanish.The underlying topology of the warped 3D visual space would be such that the perceived volume c 3 Ω of any enclosed 3D space would equal zero.In other words, it would not be possible to sense the size (i.e., volume) of enclosed 3D spaces.We know from everyday experience that we can visually sense the size of the room we are in.This is as predicted by de Rham's theorem.

Mechanical Systems
In the absence of all external forces, Hamiltonian mechanics plays out on the 2ndimensional cotangent bundle T * Q = P (or momentum phase space) of an n-dimensional configuration manifold Q.Using position (configuration) q and momentum p as canonical coordinates, a canonical 1-form θ = p i dq i and a canonical closed 2-form ω = dq i ∧ dp i exist at every point (q, p) on T * Q = P (i.e., on the momentum phase space of the Hamiltonian system).The symplectic structure of the 2-form ω = dq i ∧ dp i defines a symplectic geometry on the system.Contraction of the symplectic 2-form field ω by a vector field ξ gives a 1-form field λ = C ξ ω on T * Q = P.However, there is a wonderful discovery that applies to all Hamiltonian mechanical systems (and to all systems whose dynamics can be described by 2nd-order differential equations).There exists a special vector field ξ h for which the 1-form field λ = C ξ h ω is an exact 1-form field equal to dh, where h is a smooth real-valued function on the momentum phase space T * Q = P.The real-valued function h is called the Hamiltonian function, and the vector field ξ h is called the Hamiltonian vector field of the Hamiltonian mechanical system.If the real-valued function h is equal to the total energy of the mechanical system (i.e., equal to the sum of the kinetic energy and potential energy of the system expressed as a function of position and momentum), then the integral flow of the Hamiltonian vector field ξ h in the momentum phase space T * Q = P is equal to the free motion of the mechanical system.Since by definition there is no energy-injecting or energy-dissipating external force, the integral flow of the Hamiltonian vector field ξ h conserves the total energy h of the system.The integral flow is confined to a (2n − 1)dimensional constant energy (i.e., constant h) contour shell in T * Q = P. Different energy levels correspond to flows in different constant energy contour shells in the momentum phase space of the system.
However, if there exists a 1-dimensional hole in the momentum phase space, then, from de Rham's theorem, the 1-form λ = C ξ ω cannot be globally exact, and hence neither the Hamiltonian function h nor the Hamiltonian vector field ξ h can be globally defined!The exact 1-form field λ = dh is replaced by an harmonic 1-form field that is closed but not exact.

Electrodynamical Systems
In electrodynamics on Minkowski 4D spacetime, Maxwell's equations can be written in terms of differential forms as * F µν where F is a skew-symmetric electromagnetic 2-form at each point in 4D Minkowski spacetime, J is a current 4-vector at each point, c is the speed of light in a vacuum and * is the Hodge star operator (all in Gaussian units).The source-free equations are true if and only if there exists a 4-covector field A (i.e., a 1-form field) called the electromagnetic potential field such that F = dA (i.e., F is a closed, globally exact 2-form field on 4D flat Minkowski spacetime).
Such a 1-form field A can be altered by the addition of an exterior differential dψ of a smooth real-valued function ψ on the 4D Minkowski manifold (i.e., A = A + dψ) without changing the closed 2-form electromagnetic field F. F is unchanged because d • d = dd is always zero, and F can be written as Any such transformation of the electromagnetic potential A that has no effect on the electromagnetic field F is called a gauge transformation.Any choice of ψ that satisfies the inhomogeneous wave equation is called a Lorentz gauge; i.e., where the operator □ is the wave operator or the d'Alembertian This shows that an electromagnetic potential A in the Lorentz gauge can carry an electromagnetic wave □ψ propagating at the speed of light c in a vacuum without changing the electromagnetic field F. The electromagnetic field F is the medium in which the wave propagates with no need of an aether.It is important to notice that i.e., the d'Alembertian wave operator □ on 4D flat Minkowski spacetime is equal to the Laplacian ∆ on a flat 4D spacetime.
However, if a 1-dimensional hole exists in the 4D spacetime manifold, then, from de Rham's theorem, the electromagnetic potential A cannot be globally exact, and consequently the Lorentz gauge cannot be defined.The global electromagnetic potential 1-form field A is replaced by an harmonic 1-form field.This leaves the question about the existence of an aether uncertain.

Fluid Flow
The principles of fluid mechanics have long been a standard in the teaching of physics [14][15][16].Consider a slowly flowing stream of water.Water particles close to the boundary adhere to the boundary and so have zero downstream velocity.Water molecules in this boundary layer form hydrogen bonds with each other that increase surface tension and alter the dipole distribution of electric charge.Points away from the boundary towards the center of the stream have higher average downstream flow velocities.This pattern of flow is known as lamina flow.The random motion of water particles due to thermal energy and electric fields causes an interchange of water particles between the various laminae, with the higher-velocity particles speeding up the slower ones and vice versa.This gives rise to shear stresses and transfer of kinetic energy and momentum between the various layers producing a property of the water known as its dynamic viscosity µ.Viscosity varies between different fluids and provides a measure of the ease of flow.Honey, for example, has a higher viscosity than water and flows less easily.Now consider a fast-flowing mountain stream flowing around rocks.Water particles cannot flow through the rocks, so they have to flow around them.The rocks form 2D holes of various sizes in the space of the water (i.e., the 2D surface of each rock is a 2-cycle that forms a boundary about a hole in the space of water flow).If the stream is flowing fast enough (due to an ample supply of kinetic energy provided by gravity), flow around the rocks induces large eddies and vortices in the flow.These large eddies propagate outwards, transferring kinetic energy to smaller eddies and vortices.These in turn generate even smaller eddies.This cascade of kinetic energy from larger to smaller eddies depends on the velocity of flow (i.e., available kinetic energy) and the rate of energy dissipation in the water due to viscosity.The eddies of various sizes appear and disappear in an unpredictable fashion, causing local tangential velocities and pressures to vary from point to point in a chaotic unpredictable fashion.This increases the amount of mixing.Increased local flow velocities reduce local pressure, allowing gases in the water to vaporize and produce bubbles.Formation and popping of bubbles dissipate energy and produce a 'gurgling' sound.The chaotic pattern of flow associated with such a cascade of energy from large to small eddies with minimal energy dissipation is called turbulent flow.Turbulent flow is observed in many fluid flows in nature, such as ocean waves, storm clouds, smoke, bubbles in champagne and more.

Navier-Stokes Equations
The Navier-Stokes equations are nonlinear partial differential equations that describe the motion of viscous fluids in 3D Euclidean space.They were developed over many years during the first half of the 1800s by Navier and Stokes.In essence, they are an expression of Newton's second law of motion; that is, force equals mass times acceleration.The incompressible Navier-Stokes equations can be derived from the convective form of the Cauchy momentum equation and have the form Equation ( 6) is a continuity equation expressing conservation of momentum and applies to incompressible fluids.Solution of the equations gives the flow velocity u as a vector field over the 3D space of the fluid.The vector u gives the direction and magnitude of the velocity of the fluid at every point in the 3D space of the fluid at every moment in time.∂u /∂t is time acceleration.(u•∇)u or u•(∇u) is spatial or convective acceleration (e.g., eddies and vortices).Convective acceleration is nonlinear, describing the change in velocity u with position.It makes the equations difficult to solve and is the predominant cause of turbulence.∇u is the velocity gradient.It is a Jacobian matrix or a second-order tensor.Its symmetrical component 1 /2(∇u + (∇u) T ) describes stretching and shearing, while its skew-symmetric component 1 /2(∇u − (∇u) T ) describes the rate of rotation.ρ is fluid density, while ∇p is the pressure gradient.µ∇ 2 u is the diffusion divergence of shear stress due to viscosity, and µ is the dynamic viscosity of the fluid.The symmetrical shear stress tensor τ equals µ(∇u + (∇u) T ), and the divergence of τ is given by where ∆u = ∇ 2 u is the Laplacian of u in Euclidean space.F represents all the external forces acting on the fluid, including gravity and electromagnetic forces.The Navier-Stokes equations are difficult if not impossible to solve analytically.However, computer simulation studies using simplifying assumptions such as temporal and spatial averaging or band-pass filtering do a good job of matching simulated fluid flows with those measured experimentally [16].

Reynolds Number
The Reynolds number is a dimensionless number defined by Osborne Reynolds (1842-1917) to be Re = inertial f orces viscous f orces = ρuL µ where ρ = density, u = velocity, L = characterisitc length (e.g., size of rocks in a mountain stream) and µ = dynamic viscosity.The Reynolds number plays an important role in interpreting the Navier-Stokes equations and de Rham's theorem because as Re is increased experimentally from a small number to a large number, fluid flow is observed to undergo a transition from lamina flow to turbulent flow.The transition can be divided into stages depending on the size of the Reynolds number.When viscous forces due to shear stresses in the fluid dominate the inertial forces, the Reynolds number is small (e.g., Re for honey ≈ 0.001).When Re is small, fluid flow is dominated by energy dissipation due to the viscous shear stresses, and consequently eddies and vortices are strongly damped.The result is a stable lamina flow.As Re is increased, outward-flowing 3D pressure waves are observed between lamina layers.This increases mixing and increases energy-momentum transfer between layers.For large Reynolds numbers, inertial forces dominate and viscous damping is minimal.Turbulent flow results from the relatively low level of energy dissipation and damping.When injection of kinetic energy exceeds energy dissipation, instabilities in nonlinear convective interactions create a cascade of energy transfer from large to small eddies and vortices.
For very large Reynolds numbers, the diffusion divergence of shear stress due to viscosity µ∇ 2 u in the Navier-Stokes equation effectively can be ignored.Extreme turbulent flow results, and the cascade of energy produces chaotic rotational flows at the small scale of individual molecules.Ignoring µ∇ 2 u is equivalent to setting the Laplacian of the velocity vector field in Euclidean space equal to zero (i.e., ∆u = ∇ 2 u = 0), where ∇ 2 u equals the Laplacian ∆u in Euclidean space where the Ricci curvature is zero.According to de Rham's cohomology theory, such a velocity vector field is a global harmonic flow field associated with the existence of holes in the underlying topology of the 3D Euclidean fluid space.Given an ample supply of kinetic energy, the existence of holes (such as particles, rocks and bubbles) in the space of fluid flow is associated with extreme turbulence, due to the lack of energy dissipation in the fluid and the presence of chaotic rotational kinetic energy at molecular levels.

Biological Cells
Biological processes regulate human body temperature to 37 • C.There is thus an ample supply of thermal kinetic energy for motion of aqueous solution within biological cells.Biological cells are densely packed with organelles and other molecules.Intracellular fluid cannot flow through these organelles, so it has to flow around them.In other words, the 3D space of the intracellular aqueous solution contains a large number of holes of various shapes and sizes.As described above, de Rham's theorem shows that the presence of such holes in the underlying topology of the 3D space of a fluid has a dramatic influence on the vector and 1-form flow fields.Based on mathematical theory, we propose that the thermal motion of intracellular fluid around organelles creates a transfer of kinetic energy with minimal energy dissipation from large to very small eddies and vortices at a molecular level, and that these rotations are appropriately scaled to play a role in protein folding.

Protein Folding
The study of protein folding began some 60 years ago.Much is known, but much is yet to be learned [17][18][19].Proteins are polypeptide chains of amino acids assembled by ribosomes within biological cells.Hundreds of thousands of atoms can be involved in long peptide chains.The single-file order of amino acid side chains in each peptide chain is specified by codons within mRNA molecules.The sequence of amino acids determines the 3D structure that the peptide chain spontaneously folds into.Folding happens within milliseconds and takes place even as the peptide chain is being assembled by a ribosome.The native 3D conformation of the peptide chain is critically important.It determines the function of the protein.Different functions of thousands of different proteins provide the molecular machinery for living systems.If the proteins are not folded correctly, then they do not function properly, causing many disorders such as amyotrophic lateral sclerosis, sickle cell anemia, Alzheimer's disease and mad cow syndrome, to name a few.But what are the mechanisms controlling the folding of proteins and what causes it to go wrong?
Given the extremely large number of possible bonding angles between the atoms and molecules within the peptide chain, the fact that correct folding happens very quickly indicates that it is not a random trial-and-error process.There must be a physical process directing the folding pathway.There is also a paradox.The large number of possible configurations of the unfolded peptide chain folds into a unique 3D structure of the functioning protein.This violates the second law of thermodynamics.The conformational entropy of the protein decreases as its folded structure becomes more ordered!
The paradox can be resolved by considering the intracellular aqueous solution in which the peptide chains are immersed.This is known to play a crucial role in protein folding.If a folded protein is immersed in urea (which acts as a detergent), it unfolds, and it refolds correctly when returned into intracellular fluid.Ion concentrations, pH and temperature of the intracellular fluid have to be carefully regulated to appropriate levels by biological processes.If the regulation fails, the ribosome-assembled peptide chains do not fold correctly.A prevailing thermodynamic theory for protein folding invokes the role for intracellular aqueous solution as follows.It proposes that the lack of polarization of hydrophobic amino acid side chains distributed along the peptide chain cause the water molecules surrounding these side chains to become structured, driven by their electrical dipoles and hydrogen bonding between water molecules.Thermodynamic free-energy gradients cause the peptide chain to fold so as to enclose its hydrophobic side chains within its core, thereby hiding them from the surrounding water molecules.This releases the water molecules from the ordered structure imposed by the hydrophobic side chains, thereby increasing the entropy of the water.The increase in the entropy of the water exceeds the decrease in the configurational entropy associated with protein folding, so spontaneous folding is thermodynamically feasible.
In addition to thermodynamic free-energy gradients within the intracellular fluid, there are many other chemical and electrical potential energy fields surrounding the peptide chain that can influence its folding, including covalent, ionic and hydrogen bonds, voltage gradients and van der Waals forces between atoms in the protein.Different configurations of the peptide chain bring attractive sites along the chain into alignment.These alignments allow the formation of hydrogen bonds, disulphide bonds, salt bridges and others.Associated potential energy fields drive the folding of the protein towards the minimum energy configuration specified by the sequence of amino acids in the peptide chain.There is no guarantee, however, that a particular set of molecular bonds and salt bridges will correspond to the functional minimum-energy configuration.If it is only a local minimum-energy configuration, then work has to be done on the molecule to push its configuration out of the local energy trap and into a deeper energy valley.

The Proposal
It does not seem unreasonable to suggest that the jiggling (Brownian motion) of atoms and molecules in peptide chains caused by motion of the intracellular aqueous solution plays an important role in the kinematics of protein folding.De Rham's theorem shows that holes in the underlying topology of the intracellular fluid are associated with global harmonic flow fields in the 3D space of the fluid.The experimental study of fluid flow together with Reynolds number, the Navier-Stokes equations and de Rham's theorem provide a connection (consistent with experimental observation) among holes, extreme turbulent flow and harmonic flow fields.The term µ∇ 2 u = µ∆u in the Navier-Stokes equation describes the diffusion divergence of shear stress due to dynamic viscosity µ of the fluid.Experimental studies show that large ratios of inertial forces to viscous forces (i.e., large Reynolds number) are associated with extreme turbulent flow, and the term µ∇ 2 u = µ∆u in the Navier-Stokes equation effectively can be ignored.This is equivalent to setting the Laplacian of the velocity field ∆u equal to zero.In differential geometry, this is exactly how harmonic vector and 1-form fields are defined (i.e., in Euclidean space ∆u = ∇ 2 u = 0).In other words, experimental observations of fluid flow imply that harmonic flow fields associated with flow around organelles and other molecules (including protein molecules) correspond to extreme turbulent flow.Existence of holes in the 3D fluid space prevents closed 1-form and vector flow fields from being globally exact and replaces them with 1-form harmonic flow fields.In other words, given an ample supply of kinetic energy, de Rham's theory provides a formal description of intracellular fluid flow responsible for the Brownian motion of molecules in proteins and other structures within biological cells.
As described in Section 3.4, mathematically speaking, harmonic flow fields associated with holes in the underlying topology of the fluid space can be computed by equating the Laplacian of the flow field to zero (i.e., ∆λ := dd * λ + d * dλ = 0).But in Euclidean space, the Laplacian of a vector flow field is given by ∆u = ∇ 2 u, and according to the Navier-Stokes Equations ( 6) and (7), setting ∆u = ∇ 2 u = 0 is equivalent to the absence of viscous shear stress and absence of energy dissipation in the fluid.Extreme turbulent flow results.An ample supply of thermal kinetic energy with negligible energy dissipation in the fluid causes a cascade of kinetic energy transfer from large eddies and vortices to smaller and smaller ones until finally a large amount of kinetic energy is available in the form of chaotic 3D rotations of individual water molecules.Remembering that water molecules form hydrogen bonds to each other in the neighborhood of hydrophobic amino acid side chains and adhere to hydrophilic amino acid side chains, it seems reasonable to propose that this kinetic energy is available to do work (i.e., inject energy) on the atoms and molecules of protein polypeptide chains.
The rotational torques applied to the protein molecules are at the right energy level and scale to operate on the bonding angles between atoms and molecules within the peptide chain.Moreover, they operate simultaneously on all parts of the peptide chain and, when combined with thermodynamic, chemical and electrical potential energy gradients surrounding the peptide chain, can play an important role in protein folding.They not only supply the energy needed to break away from local energy traps but also can accelerate folding towards minimum energy configurations.Like 'annealing' in minimum-variance search algorithms for neural networks, the chaotic harmonic flow fields within intracellular fluids enable the folding to find the unique absolute minimum energy configuration of the protein.Of course, this requires the appropriate biologically regulated level of available energy.If body temperature is too low, the intracellular fluid flow has insufficient energy to break peptide chains away from local energy traps, and if body temperature is too high, it can break bonds and denature proteins.
A nice analogy for the Brownian motion of protein molecules and side chains caused by turbulent flow of intracellular fluid about the densely packed organelles within the cell is given by the chaotic movement of the leaves and branches of a tree generated by the flow of wind through the tree.

Discussion
To recap, de Rham's theorem shows that holes in the underlying topology of a smooth manifold are associated with closed differential form fields on the manifold that are not exact.Known as harmonic form fields, these can be computed on a Riemannian manifold by setting their Hodge Laplacian ∆ := dd * + d * d equal to zero, and on a Euclidean space by setting the Laplacian ∆ = ∇ 2 equal to zero.The question is, do such harmonic form fields exist in nature and, if so, what properties do they possess?
The mathematics tells us that the integral flow fields on a sphere, for example, will be very different from the flow fields on a torus, the difference being attributable to the hole in the torus.Likewise, for the physical systems described in Section 3, the same mathematics predicts with certainty that the behavior of vector and differential form fields on their smooth manifolds should be profoundly altered by the appearance of any hole or holes in the underlying topology of the manifold.Do we have experimental evidence that this is the case?And what are the implications for modeling these systems?

Mechanical Systems
In Section 3.2, we showed mathematically that, according to de Rham's theorem, the presence of holes in the momentum phase space of a mechanical system prevents the Hamiltonian function h and the Hamiltonian vector field ξ h from being globally defined on the phase space.However, despite this being a definitive mathematical prediction, until appropriate experimental data become available, it remains uncertain as to how exactly a physical mechanical system with holes in its momentum phase space might behave.What type of flow field will replace the globally exact 1-form field λ = C ξ h ω in a Hamiltonian system is yet to be discovered.In other words, de Rham's theorem implies that changes in the topology of the state space of a mechanical system require a different mathematical model from that of Hamiltonian mechanics.
This problem applies similarly to Lagrangian mechanics.While Hamiltonian mechanics plays out on the cotangent bundle (T * Q) over the configuration manifold (Q), Lagrangian mechanics plays out on the tangent bundle (TQ) over the configuration manifold (Q).Legendre's transformation FL = TQ → T * Q maps points in TQ to points in T * Q. Lagrangian and Hamiltonian mechanics are often described as being dual to each other because, if the equations are everywhere regular, both derive the same trajectory of motion in the configuration manifold Q.However, duality between Lagrangian and Hamiltonian systems requires the Legendre transformation FL( i → p i to be a smooth one-to-one diffeomorphism, where L(q, .q) is the Lagrangian equal to the difference between the kinetic energy and potential energy, .q i is the i − th component of the velocity vector .
q and p i is the i − th component of the momentum vector p.Duality also requires that both the symplectic Lagrangian second differential form on TQ and the Hamiltonian symplectic second differential form on T * Q be nondegenerate.However, many real-world mechanical systems do not satisfy these requirements.Indeed, it is this fact that underlies the importance of Emmy Ntre's theorem; namely, if a Lagrangian L(q, .q) is degenerate and does not vary along a coordinate q i , then the associated component of momentum p i is conserved.Removing degeneracies and singularities from the manifolds introduces holes, and these, as shown by de Rham's theory, change the behavior of the mechanical system away from that predicted by classical Lagrangian and Hamiltonian mechanics.As far as we know, the modeling of real-world mechanical systems to account for holes in the manifold is yet to come.Hopefully this setting out of the theoretical background can assist in stimulating the collection of the necessary experimental data.

Electrodynamical Systems
In Section 3.3, we showed mathematically that the presence of holes in Minkowski 4D spacetime prevents Maxwell's source-free equations from being true.The closed 2-form electromagnetic field F is prevented from being globally exact, and the electromagnetic 1-form potential field cannot be defined globally.Consequently, global gauge transformations of the potential field, including the Lorentz gauge, cannot exist.Will the 1-form electromagnetic field be replaced with a harmonic 1-form field and, if so, what properties will it possess?Again, appropriate experimental observations of physical electrodynamical systems are required to answer these questions.
Indeed, the implications go to the very nature of electromagnetic radiation.The churning motion of charged particles in the sun creates changing electric currents J and changing charge densities ρ.But, as shown by Maxwell's equations of electromagnetism, J and ρ are exactly the energy sources for magnetic fields B and electric fields E. In the vacuum of space between the sun and the earth, J and ρ are both zero.Nevertheless, using the curl of the curl relationship and the div and curl relationships between E and B in Maxwell's equations, even with both J and ρ set equal to zero, the coupled E and B fields satisfy wave equations ∇ 2 E − 1 /c 2 ∂ 2 E/∂t 2 and ∇ 2 B − 1 /c 2 ∂ 2 B/∂t 2 .In other words, E and B form coupled electromagnetic waves that can travel at the speed of light c through a vacuum despite both J and ρ being zero.
But in physics there remain unresolved questions.By analogy with waves in water and sound waves in air, it can be seen that waves require a medium in which to travel.But what is the medium in which an electromagnetic wave in a vacuum travels?In the late 1800s, it was proposed by Lorentz that something called the aether filled all of space and that electromagnetic waves travel at the speed of light in the aether.However, the 1887 interferometer experiment of Michelson and Morley showed that there was no difference in the speed of light coming from a star when the earth was moving towards the star compared to when the earth was moving away from the star.This challenged the idea that both the light and the earth are moving in space relative to a fixed aether.More recently, the mathematical approach to particle physics known as gauge theory offers an alternative explanation, namely that the electric and magnetic fields themselves provide the medium in which light travels.Yet as described in Section 3.3, de Rham's theory of differential forms shows that gauge theory itself fails to hold in the vicinity of holes in the manifold.And this is just what the presence of stars and planets create in the fabric of spacetime.
Life on earth depends fundamentally on the ability of electromagnetic waves to carry energy through empty space.Maxwell was quick to assert that light is an electromagnetic wave.Via photosynthesis, plants absorb energy from the light that falls upon them, transforming this into chemical energy stored in sugar molecules.This stored chemical energy in turn provides most of the biological energy needed to sustain life on earth.The understanding of electromagnetic fields and the consequences of the occurrence of holes within them is therefore an important pursuit.

The Human Visual System
For both mechanical and electrodynamic systems, the implications of de Rham's theorem are yet to be addressed by experiment.For the human visual system examined in Section 3.1, however, the situation is somewhat different.In this case, we have extra relevant information derived from our own conscious awareness.We know, for example, that there is a hole about the egocenter in our perceived binocular 3D visual space because we are unable to see our own head.We also know that the size of a perceived visual image increases as the Euclidean distance between the egocenter and the object in the environment decreases.This is consistent with experimental data showing that the size of an image on the retina changes in inverse proportion to the Euclidean distance between the nodal point of the eye and the object in the environment [4].This creates a singularity at the egocenter in 3D visual space where the size of the image on the retina becomes infinite.But this singularity is not a problem, because it is contained within the hole about the egocenter.De Rham's theorem shows that such a hole should prevent perceived image flow fields from being globally exact, replacing them with harmonic flow fields.
The question as to whether or not such harmonic flow fields exist in the human visual system and what properties they possess can be answered to some extent from personal experience.As an object in the outside world is brought increasingly close to the egocenter, we can keep it in the binocular visual field by converging the eyes.The perceived size of the image gets bigger and bigger and starts to exceed the size of the perceived visual space.Eventually, we reach a limit to cross-eyed vision.The perceived visual image becomes defocused and disintegrates in a chaotic fashion.It can no longer be perceived as a coherent single image.It becomes physically impossible to bring the object any closer to the egocenter.It has reached the edge of the hole in perceived 3D binocular visual space.Thus, personal experience shows that harmonic visual-image flow fields in 3D visual space disintegrate in an unpredictable chaotic fashion in a region close to the boundary of the hole about the egocenter; in other words, this is a boundary layer phenomenon.
At this point, it is salient to note that experimentation in human visual science tells us that visual space can have multiple mappings.We argue that the fundamental mapping is that revealed by the surprising fact that there are people who, as the result of brain damage, are visually blind but nevertheless can move about in their 3D environment and perform visually guided movements with precision just as if they could see.This phenomenon has been called blind sight [20].It demonstrates that the brain has a system able to form a subconscious egocentric 3D visual space that has a one-to-one mapping with the Euclidean outside world.This subconscious visual system has to take into account the warping of this 3D visual space, a warping that occurs because image detection by the retinas occurs on a 2D surface.More precisely, the warping is due to the size of an image on the 2D retinas varying in proportion to the angle subtended at the nodal point of the eye by the corresponding object in the 3D Euclidean outside world [4].As a result of this warping, changes in the place of the head in the environment and changes in the posture of the body change not only the computed size and position of encoded visual images of objects in the environment and of the body in that environment, but they also change the outline, occlusions, shape, curvatures, relative velocities and relative accelerations of what is seen.
In other words, there exists a conformal mapping between the Euclidean outside world and the warped subconscious 3D visual space.A conformal mapping preserves the infinitesimal shape of objects and angles between coordinates, but sizes of the encoded images are not preserved.Images of objects shrink in size as the relative distance between the egocenter and the object in Euclidean space increases.Nevertheless, despite this warping, the conformal mapping with the outside world is one-to-one, onto and smooth, and the same warping applies to visual images of the environment and visual images of the body or parts of the body in that environment.Thus, this warped, subconscious, egocentric 3D visual space can be used for the planning and execution of visually guided movements.
That said, there is much evidence that the subconscious and conscious perceptions of the visual world are not the same.The latter can be mediated by top-down cognitive processing of visual images, held to be based on Bayesian statistical hypotheses derived from experience over time.This introduces all sorts of illusions into the consciously perceived visual world.A favorite is the concave face-mask illusion.Because of a lifetime of looking at faces that are convex in shape, when a person looks at the inside concave surface of a face mask, they see it consciously as being convex.Moreover, the illusory convex face seems to turn and follow them as they move about.A consciously perceived visual space containing illusions is not what is needed for the planning and control of precision visually guided movements.The subconscious 3D visual space, derived directly from afferent signals and free of statistical cognitive intervention, is the more appropriate for visually guided movement [1][2][3].
Experimental support for this description of the human visual system comes from an experiment first carried out by Taylor in 1941 but repeated many times since (see [3]).In the dark, a person holds their hand in front of their face.A brief flash of bright light then burns (or bleaches) an image of the hand onto the retinas.This gives rise to a perceived afterimage of the hand that can persist for several tens of seconds.If, in the dark, the person moves the hand towards the face or moves the face towards the hand, a remarkable thing happens.The afterimage of the hand appears to decrease in size even though the image bleached onto the retinas remains unchanged.This is the opposite of what usually happens because, as the hand is moved closer to the face in light, it appears to increase in size.Our explanation is that the nervous system can transform proprioceptively sensed movements of the body into anticipated changes in its visual image [3].If the anticipated visual image is subtracted from the actual visual image, the perceived body image does not appear to change its size but only its position in space (so called size-constancy).During the afterimage experiment, the movement of the hand generates the anticipated change in the retinal image of the hand, but the actual retinal image is burnt onto the retinas and so does not change.Subtracting the anticipated size of the image from the actual retinal image generates an illusory increase in the size of the afterimage.
We have constructed a Riemannian geometry description of the calculations required to transform movements of the body sensed proprioceptively into anticipated changes in 3D visual images of the body [1][2][3].This is not a simple calculation, because the relationship is a multivariable nonlinear one.It requires 116 variables to specify the configuration of the body moving in a local 3D environment.For each place and posture of the body there exists a 3D warped visual space, and at each point in that visual space there exists an encoded visual image of the environment.These visual images change from one place and posture to another.Geometrically speaking, each 3D visual space corresponds to an egocentric negatively curved Riemannian manifold.
This relationship between movement and vision is complicated by the fact that the visual space is missing the 3D ball of points from about the egocenter [1,2].This intrinsic hole is what inspired us in the first place to consider how differential geometry and de Rham's theorem in particular might illuminate this hole's influence on the perceived flow of visual images.The chaotic breakdown of visual images near the boundary of the hole was described above.

Fluid Flow
Turning to Section 3.4, our best interpretation to date of the properties of real-world harmonic flow fields comes from the experimental studies of fluid flow.As set out in that section, Navier-Stokes equations describe fluid flow in 3D Euclidean space where the Laplacian of the vector flow field u is given by ∆u = ∇ 2 u.The harmonic vector field associated with the presence of holes (e.g., particles, rocks, bubbles, etc.) in the Euclidean space of a fluid is obtained by setting the Laplacian to zero (i.e., ∆u = ∇ 2 u = 0).When this is substituted into the Navier-Stokes equation, extreme turbulence is obtained consistent with experimental observations of extreme turbulence associated with large Reynolds numbers.Inertial forces ρuL, velocities, pressures, temperatures and kinetic energy all change in a chaotic fashion in time and from point to point in space.As the Reynolds number of a fluid is increased experimentally, the energy-dissipating viscous forces within the fluid become negligible and the behavior of the fluid flow approaches that of a theoretical harmonic flow field.Because of the lack of viscous damping within the fluid, a cascade of kinetic energy transfer occurs from large to small eddies and a chaotic flow field results.Rotational vortices and rotational energy become available at the small scale of individual molecules in the fluid.The presence of such small-scale eddies in extreme turbulent fluid flow created by the presence of holes in the space of the fluid suggests a role in protein folding for extreme turbulence in the intracellular fluid of biological cells due to holes in the fluid created by the densely packed organelles within the fluid.
Inertial forces changing from point to point in the fluid give rise to a symmetrical mass-inertia matrix corresponding to a kinetic energy Riemannian metric tensor [21].This provides an opportunity to improve on the Navier-Stokes equations.A Riemannian kinetic energy metric allows the fluid flow to be described by vector fields on a 3D curved Riemannian manifold rather than in the flat 3D Euclidean space of the Navier-Stokes equations.As shown in Section 3.4, the harmonic 1-form flow field λ on a Riemannian manifold is obtained by setting the Laplacian of the 1-form λ on the Riemannian manifold to zero (i.e., ∆λ = (− ∑ n i=1 ∇ 2 ξ i )λ − Ric (λ) = 0).The Ricci curvature at each point on the Riemannian manifold takes into account nonlinear dynamics attributable to inertial forces changing from point to point in the fluid.These are the nonlinear forces within the fluid that are predominantly responsible for turbulent flow.Thus, describing fluid flow on a Riemannian manifold endowed with an appropriate Riemannian metric will give a more detailed description of fluid flow than given by the Navier-Stokes equations.However, because of the chaotic nature of turbulent flow, stochastic methods are required to obtain estimates of an appropriate kinetic energy Riemannian metric as a function of position in the manifold.A stochastic time series model of the changing Riemannian metric appears necessary.
Deep learning neural networks have been used to model turbulent fluid flow [16].Unlike a "black box" solution, however, differential geometry offers a more general mathematical account of turbulent flow that can provide an understanding of the underlying physical phenomena.

Protein Folding
There exists ample experimental evidence showing that intracellular aqueous solution plays an important role in protein folding.However, the availability of data that can be used to verify or refute working hypotheses about the precise nature of water-protein interactions is limited by the experimental difficulties involved.That said, available results provide no support for the popular hydrophobic theory of protein folding.Instead, the data support water forming hydrogen bonds between points in the peptide chain and then pulling those points towards each other.Work to measure the role of water in protein folding by measuring neutron diffraction patterns is promising, but the technique is not an easy one, requiring stochastic analysis of measured data and use of computer simulations to find the model that best fits the data.Meanwhile, the general conclusion that movement of water does play a role in protein folding seems beyond doubt [19,22,23].
Our proposal regarding the Brownian motion of bonding angles within peptide chains associated with turbulent rotation of water molecules flowing around organelles within the cell is a novel hypothesis.It is based solely on the mathematics of fluid flow (i.e., the Navier-Stokes equations) in conjunction with de Rham's theorem about the influence of holes on flow fields.It awaits experimental testing that is beyond the resources of the present authors.Meanwhile we are convinced of the following.The logical formalism and certainty within the mathematical formulation set out here make it highly likely that the experimentally observed Brownian motion in the intracellular aqueous solution during protein folding is contributing fundamentally to the folding process.We hope that this can provide a guide to future research.

Conclusions
An aim of this paper has been to demonstrate the worth of applying differential geometry to gain understanding of any physical system where holes can occur in a surrounding flow field.According to de Rham's theory, holes alter the character of the surrounding fields, transforming them into harmonic fields that are closed but not exact.As a result, the consequent flow of a fluid becomes chaotic.We have discussed this phenomenon in relation to 3D binocular vision, which breaks down due to a hole in the visual flow field about the egocenter.We now point out that the presence of electromagnetic screening about a source of electromagnetic radiation, such as a television transmitter, creates a hole in the surrounding 3D space of electromagnetic fields.We suggest that it was the creation of chaotic flow fields that underlaid the television studio issue [5,6] that in part led us to examine this theory in the first place.No cause of the cancer cluster came to light in the investigations carried out, but the studio was ultimately abandoned [24].
For mechanical and electrodynamical systems where flow data are apparently lacking, we hope that the ideas presented here will be of interest to those working directly in those fields.In particular, we are optimistic that our proposal about protein folding might stimulate a desire to know more about differential geometry and de Rham using the tutorial information written here.In summary, we trust that our demonstration of the reach of de Rham's cohomology theorem across diverse real-world systems can help to focus future attention in general on this extraordinary piece of mathematics.After all, the notion that the same mathematics can be applied across a range of systems is what underlies a wealth of modeling and simulation, be it with mechanical, electrical, chemical, economic, physiological, sociological or cosmological variables, etc. Differential geometry has found a key place in physics for those on a quest to understand nature.For those elsewhere who have mathematical inclinations along with a foothold in a relevant field, the work involved to understand and apply differential geometry promises ultimately to be very rewarding.
To illustrate how the Mayer-Vietoris theorem can be combined with other properties of topological spaces to compute singular homology groups, consider computing the singular homology groups for all values of r and for all dimensions n of a unit sphere.Consider a unit n-sphere S n with open subset U equal to S n with the south pole removed and open subset V equal to S n with the north pole removed.A part of the Mayer-Vietoris exact sequence is where the connecting homomorphism is All points in U can be contracted to the north pole and all points in V can be contracted to the south pole.This means that both U and V are contractible and therefore simply connected (i.e., contain no holes).It follows that Consequently, when r > 1 and n ≥ 1, the exact sequence can be rewritten as This is an exact 4-term sequence as described above and, consequently, the connecting homomorphism ∂ * : H r (S n ) → H r−1 (U ∩ V) is a bijective map (i.e., an isomorphism).Just as the equator of a unit 2-sphere S 2 is topologically equivalent to a unit circle S 1 , the equatorial sphere (U ∩ V) of a unit n-sphere S n is topologically equivalent to a unit (n − 1)sphere S n−1 .This means H r (S n ) ∼ = H r−1 (U ∩ V) ∼ = H r−1 (S n−1 ) for r > 1 and n ≥ 1 and H r (S n ) ∼ = 0 if 0 < r < n.We already know that H 0 (S n ) ∼ = Z and H r (S n ) ∼ = 0 if r > n, so we have all the singular homology groups for higher-dimensional spheres.
Appendix A.7. Homotopy Invariance A 2-torus T 2 has two one-dimensional holes and one two-dimensional hole.In other words, it has a two-dimensional degree-one singular homology group H 1 (T 2 ) and a onedimensional degree-two singular homology group H 2 (T 2 ).The rubber sheet geometry of the 2-torus allows it to be deformed into many different sizes and shapes.As long as it is not cut, torn, punctured or glued, all the different sizes and shapes of T 2 can be mapped homeomorphically onto each other and the topological structure remains unchanged (i.e., the number of holes and hence the singular homology groups H 1 (T 2 ) and H 2 (T 2 ) remain unchanged).A similar property holds for all topological spaces.The number of holes of various dimensions, and hence their singular homology groups, are homotopically and topologically invariant.Let f 0 , f 1 : X → Y be continuous homotopic maps (i.e., f 0 (X) in Y and f 1 (X) in Y can be morphed into each other without changing their topological structure), then f 0 (X) and f 1 (X) (together with all the intermediate shapes along the homotopy transformation pathway between f 0 (X) and f 1 (X)) have the same singular homology groups.In other words, continuous homotopy maps f 0 , f 1 : X → Y induce singular homology group homomorphisms ( f 0 ) * , ( f 1 ) * : H r (X) → H r (Y) for all values of r.A continuous map F: X → Y between topological spaces X and Y is said to be a homotopy equivalence if there is a continuous pull-back map G: Y → X such that the composition F • G is homotopic to the identity map I Y on Y and the composition G • F is homotopic to the identity map I X on X.If F: X → Y is a homotopy equivalence between topological spaces X and Y, then for each r ≥ 0, the map F * : H r (X) → H r (Y) is an isomorphism.star operator, that is, an isomorphism between Ω r and Ω n−r .The operators d and d * are not inverses of each other, they are adjoints.The direct sum of skew-symmetric tensor spaces forms a graded algebra Ω = Ω 0 • • • Ω n , known as an exterior algebra or Grassmann algebra.Any r-form ω for which dω = 0 is said to be a closed r-form, while any r-form ω equal to dλ where λ is an (r − 1)-form is said to be an exact r-form.Appendix A.11. De Rham Cohomology Groups From Grassmann algebra and the fact that d • d = dd = 0, it follows that every exact differential form is closed.It does not follow, however, that every closed differential form is exact.De Rham's theory of differential forms makes it clear that whether or not a closed differential form is exact depends uniquely on the underlying topological structure of the manifold.Specifically, it depends on the presence or absence of holes of various dimensions in the underlying topology.As illustrated by the amount of theory summarized in this appendix, a significant amount of differential geometry is needed to appreciate this remarkable fact.Because the exterior differential d : Ω r → Ω r+1 is linear, its kernel ker(d r ) = Z r (M) = closed r-forms and its image im(d r−1 ) = B r (M) = exact r-forms are linear vector spaces.Since every exact form is closed, it follows that the exact r-forms are a subspace of the closed r-forms (B r (M) ⊆ Z r (M)).Thus, the r-th de Rham cohomology group can be defined to be the quotient vector space H r dR (M) = Z r (M) /B r (M) (i.e., the space of all closed forms that are not exact).While H r dR (M) is a real vector space, it is also a Lie group under addition.Since Ω r (M) = 0 for r < 0 and r > n, it follows that H r dR (M) equals zero for these values of r.It also follows that H r dR (M) = 0 for 0 ≤ r ≤ n if and only if every globally closed r-form on M is globally exact.
The vector quotient space Z r (M) /B r (M) partitions the space of closed r-forms into equivalence classes.For any closed r-form ω on M that is not exact, we let [ω] denote the equivalence class of ω.If ω and ω ′ belong to the same equivalence class, [ω] then they differ by an exact r-form (i.e., ω − ω ′ = dη), and we say that ω and ω ′ are cohomologous.This is easy to visualize in terms of homology groups, where it simply means that there can exist multiple r-cycles that contain the same r-dimensional hole.Just as a continuous map f : X → Y between topological spaces X and Y induces a pull-back homomorphism f * : H r (Y; G) → H r (X; G) between the r-th singular cohomology group H r (Y; G) on Y and the r-th singular cohomology group H r (X; G) on X, a smooth map F : M → N between smooth manifolds M and N induces a pull-back F * : H r dR (N) → H r dR (M) between the r-th de Rham cohomology groups of N and M, respectively.This corresponds to a contravariant functor between the category of smooth maps between smooth manifolds and the category of de Rham cohomology groups.
We can form an interesting picture of a relation between degree-one de Rham cohomology groups H 1 dR (M) on a connected smooth manifold and the underlying topology X of the manifold as measured by the fundamental group Π 1 (X, q).All the 1-cycles in all the possible 2-chains in the underlying topology of the smooth manifold M that contain a 1-dimensional hole correspond not only to H 1 dR (M) but also to the fundamental group Π 1 (X, q).The fundamental group consists of all the loops associated with the point q in the underlying topology X of M that cannot be contracted to a point.The set of group homomorphisms Hom(Π 1 (X, q), R) from the fundamental group Π 1 (X, q) to the additive group of real numbers R is a vector space under addition and multiplication by constants.We can define a linear map Φ : H 1 dR (M) → Hom(Π 1 (X, q), R) as follows: Given a cohomology class [ω] ∈ H 1 dR (M), define Φ[ω] : Π 1 (X, q) → R by Φ[ω] = γ ω = a real number, where [γ] is any path homotopy class in Π 1 (X, q) and γ is any piecewise smooth curve representing the same loop path class [γ].The map Φ : H 1 dR (M) → Hom(Π 1 (X, q), R) is a well-defined homomorphism between groups.It follows that H 1 dR (M) = 0 when the topol- = 0 (since dω = 0).

Figure 2 .
Figure 2. A diagram illustrating the fact that if F, G : M → N are diffeomorphic maps between smooth manifolds M and N, then the pull-back maps F * , G * : H r dR (N) → H r dR (M) between the de Rham cohomology groups H r dR (N) and H r dR (M) for all values of r are equal, i.e., F * = G * .In other words, if the maps F and G are diffeomorphic, they have the same homeomorphic topological structure and so they are homotopy-invariant.This means that homotopy-equivalent manifolds have isomorphic de Rham groups.The labels Ω r−1 (M), Ω r (M) represent the skew-symmetric tensor spaces of degree (r − 1) and degree (r) for differential forms at each point on the smooth manifold M. The labels Ω r−1 (N), Ω r (N), Ω r+1 ( N) represent the skew-symmetric tensor spaces of degree (r − 1), degree (r) and degree (r + 1) for differential forms at each point on the smooth manifold N. The maps h : Ω r+1 ( N) → Ω r (M) for all values of r are called homotopy operators.Examination of the diagram shows that if the r-form ω is closed, i.e., dω = 0, then d(h ω) = G * ω − F * ω.In other words, G * ω and F * ω differ by an exact differential d(h ω).The maps F * , G * : H r dR (N) → H r dR (M) are said to be cohomologous.This is exactly how equal cohomology groups are defined.Therefore, cohomologous groups have isomorphic de Rham cohomology groups.