Maximum Entropy (Most Likely) Double Helical and Double Logarithmic Spiral Trajectories in Space-Time.

The ubiquity of double helical and logarithmic spirals in nature is well observed, but no explanation is ever offered for their prevalence. DNA and the Milky Way galaxy are examples of such structures, whose geometric entropy we study using an information-theoretic (Shannon entropy) complex-vector analysis to calculate, respectively, the Gibbs free energy difference between B-DNA and P-DNA, and the galactic virial mass. Both of these analytic calculations (without any free parameters) are consistent with observation to within the experimental uncertainties. We define conjugate hyperbolic space and entropic momentum co-ordinates to describe these spiral structures in Minkowski space-time, enabling a consistent and holographic Hamiltonian-Lagrangian system that is completely isomorphic and complementary to that of conventional kinematics. Such double spirals therefore obey a maximum-entropy path-integral variational calculus ("the principle of least exertion", entirely comparable to the principle of least action), thereby making them the most likely geometry (also with maximal structural stability) to be adopted by any such system in space-time. These simple analytical calculations are quantitative examples of the application of the Second Law of Thermodynamics as expressed in geometric entropy terms. They are underpinned by a comprehensive entropic action ("exertion") principle based upon Boltzmann's constant as the quantum of exertion.

Introduction

In 1823 Olbers formulated the problem of the dark night sky as a paradox of cosmological geometry[1], but its overriding significance to us as living beings is its functioning as an entropy engine[2]. Landauer’s seminal work[3] (following Shannon[4] and Brillouin[5] teaches us that information has calculable entropy and obeys physical laws, while the introduction by Jaynes[6] of maximum entropy (MaxEnt) as the basis of the rules of thermodynamics (for example, the determination of the partition function) is now recognised as far-reaching. The associated variational approach to entropy production first described by Onsager[7] also provides critical insights into issues of thermodynamic reciprocity and symmetry in systems far from equilibrium.

Today, the entropic treatment of information is standard in the analysis of the efficiency of communications networks in the presence of noise[8], and it has become clear that information and its transfer are associated with discontinuities[9], implying non-adiabatic (entropy changing) conditions. Indeed, Brillouin considered information (negative entropy, or negentropy[5]) to be anti-correlated with entropy, and Bennett[10] showed elegantly how information erasure has an entropy cost: note that perfect information copying is excluded by the “no-cloning theorem”[11]. Applying Landauer’s Principle[3] to a computation involves the transfer of information and therefore also results in a rise in entropy[12].

Parker & Walker[9] have also shown that holomorphic functions (which allow analytical continuation, with the sum of the Cauchy residues being zero) by themselves cannot transmit information. Instead, network stability[13] considerations require a filter function’s denominator be a Hurwitz polynomial[14], resulting in the necessity for meromorphic functions (that is, complex functions analytical everywhere except at isolated poles, with a non-zero sum of Cauchy residues) being used to represent information transfer.

We will show that certain geometrical structures with simple analytical representations – the double helix and the double logarithmic spiral – can be treated formally as holomorphic; and further, we calculate their geometric entropy with Lagrangian methods (based on a calculus of spatial gradients) showing that the appropriate Euler-Lagrange equations are satisfied, that is, they are maximum entropy structures. Then, to verify the formalism developed, we will calculate certain observable quantities conforming to the Hamiltonian and Lagrangian equations of state and show consistency with real observations.

Holomorphic Info-Entropy

The simplest meromorphic function is functionally equivalent to an isolated singularity, that we place in a Minkowski space-time, described by basis vectors (γ, μ ∈ {0, 1, 2, 3}) which obey a Clifford algebra that formally distinguishes the special behaviour of the time axis γ0, being characterized by a real time axis and imaginary space axes (see Penrose, ch.18[15]; we follow Penrose’s choice of metric). An information vector h can be defined in Minkowski 4-space, and can be shown to be obtained from the sum of the temporal residues h associated with each spatial basis vector γ, n ∈ {1, 2, 3}, (see Appendix A Eq. A.4 in Supplementary Information) given by:

Note that we use Einstein’s summation convention in Eq. 1a (and subsequently, where stated) using tensor index notation where the lower index indicates the row and the upper the column. The bivectors σ ≡ γ γ0 also represent unit vectors along the co-ordinate axes of the 3-dimensional space, forming a quaternion sub-algebra isomorphic to the Pauli spin vectors with the associated pseudoscalar I = σ1σ2σ3, where I2 = −1. Mathematically, this has transformed our starting Euclidean geometry into what will turn out to be a much more useful hyperbolic geometry. Penrose (§2.4)[15] emphasises that such a logarithmic representation is characteristic of hyperbolic geometry, and we see here its intimate relationship with entropic quantities.

We choose to define the entropy s as the Hodge-dual *h of the information since this definition can be shown to have the correct properties; note that Penrose (§19.2[15]) points out that Maxwell’s equations are self-dual in the orthogonal complement sense of the Hodge-dual operation, with σ = *σ = Iσ:

Thus we amplify Brillouin’s assertion of the close relation of information with entropy by treating entropy mathematically as an orthogonal complement of information.

We choose entropic structures exhibiting a transverse helical geometry, that is, s3 = h3 = 0, with a “trajectory” axis (plane waves travelling) in the γ3 direction. Then, given that s and h are conjugate (that is, the orthogonal complements of each other), the entropy eigenvector can be written as (for the right-handed chirality; see Appendix A Eq. A.6b in Supplementary Information)and its (conjugate) information term similarly written as (Eqs A.7b and A.9c)

Note that Eqs. 1 treat the generalised singularity of an isolated pole, whereas Eq. 2 constrain this singularity into a geometry isomorphic with the double-helix implied by Maxwell’s equations.

Courant & Hilbert[16] point out that the Maxwell equations are a hyperbolic version of the Cauchy-Riemann equations, and Salingaros points out that the vacuum electromagnetic (EM) field is holomorphic[17]. To form a holo-morphic info-entropy function we combine together the expressions in Eqs. 2 for information and entropy in the same way (and for the same reason) that is done in the Riemann-Silberstein[18,19] complex-vector (holomorphic) description of the EM field:where and are the 1-vector electric and magnetic fields; is a bivector (see Penrose[15] §19.2), hence the need for γ0. The equivalent complex-vector for the bivector info-entropy case is:so that we have, from Eq. 2 (see Appendix A, Eq. A.10b in Supplementary Information):

Note that the argument of the logarithm is now dimensionless, as is conventional. Note also that meromorphic functions are only piecewise holomorphic, so they can transmit information.

Just as Maxwell’s equations have a complementary (dual, in a strong sense) helical structure of the electric and magnetic fields, we continue to choose a similar double-helical structure to the info-entropic geometry, such that the loci of the x1 and x2 co-ordinates of the info-entropic trajectory are related to each other by a pair of coupled differential equations:where the coupling parameter is given by κ0 ≡ 2π/λ0 with λ0 being the helical pitch along the γ3-axis (that is, the x3 direction) and the prime indicating the differential with respect to x3 (the trajectory axis) x′ ≡ dx/dx3 as usual.

In the entropic domain the x3 co-ordinate plays a role analogous to that normally played by time t in conventional kinematics: to amplify this point, note that x0 ≡ ct and x3 are also commensurate conjugates in the Pauli algebra (see Eq. A.6a in Appendix A in Supplementary Information). Considering only the functional part of the complex-vector, Eq. 4 allows us to write the ‘local’ geometric entropy for a double-helical structure as (Eq. A.12):which is functionally equivalent to Boltzmann’s equation for entropy; where the quantity W ≡ x′/κ0 x therefore represents the number of states available for the nth plane wave.

We now consider the case of the double helix in more detail, and in particular as exhibited by the structure of DNA (which is naturally right-handed). Without loss of generality, we define the locus in space l1 of the first information-bearing helix of DNA with its axis aligned to the γ3 direction:where R0, κ0 and x3 represent respectively the radius, pitch, and axial co-ordinate of the helix. The second helix l2, with its complementary base-pairing and anti-parallel (C2 space group) symmetry contains the same entropic information content as l1, but π/2 phase-shifted and propagating in the opposite (i.e. negative) γ3 direction:

These expressions are mathematically equivalent to those for the electric and magnetic fields of an EM wave, with l1 and l2 being complementary. Equivalent to Eqs. 3, we now express the double-helix as the complex-vector ∑ = l1 + il2 to describe a single holomorphic trajectory in Euclidean coordinates with spatial basis vectors γ (n ∈ {1, 2}):

We therefore see in Eq. 8 the functionals represented by x1 = R0 exp(iκ0x3) and x2 = −iR0 exp(iκ0x3), from Eq. (5), where the phase and sign difference between x1 and x2 are typical for a pair of coupled mode equations, and which together form a holomorphic function (see Appendix B Eq. B.1 in Supplementary Information).

Hyperbolic Geometry & Entropic Momentum

We now exploit Penrose’s assertion (§2.7 p.48[15]) that there is a “hyperbolic overall geometry of the spatial universe … the space of velocities … is certainly a three-dimensional hyperbolic geometry” (his italics; this assertion is underpinned by extensive observations of the cosmic microwave background). So we define for our helix the “hyperbolic position” vectors q in the simplest possible way that involves the logarithm characteristic of the hyperbolic geometry (see Eqs. 1), where the logarithm is kept dimensionless by the normalising (Euclidean) metric R (see Appendix B, Eq. B.2 in Supplementary Information):

For small geometry (x ≪ R) and for x having its origin at R such that x tends to R + x (that is, where x is localised in the vicinity of R) the hyperbolic geometry is approximately Euclidean, q ≈ x, and also independent of the metric R. For the double helix geometry we take R = R0 for n ∈ {1, 2}.

The conjugate quantity for position q is the momentum p, so that moving towards a Lagrangian formalism, we therefore also define the “entropic momentum” p vectors in terms of an “entropic mass” m and the velocity q′, where as before q′ ≡ dq/dx3. Note that q′ is dimensionless, so that either q′ or its inverse 1/q′ can be used as a “velocity” (this ambiguity is a feature of hyperbolic velocities). It turns out that the inverse definition is more fruitful (see Appendix B, Eq. B.6):where the entropic mass m is defined as:and the Boltzmann constant k is introduced on dimensional grounds as the entropic analogue to Planck’s constant in kinematics. We use the subscript ‘S’ as a reminder that a quantity is entropic. Clearly iκ0k is a geometric quantity intrinsically based upon the pitch of the double helix. Simple calculus on Eq. 9a allows us to create the useful auxiliary identity q′ = R · x′/x, again highlighting the intimate relationship between Eqs 6 and 9; we will show elsewhere[20] how Liouville’s theorem allows the conjugate variables p and q to be used to calculate the entropy of the geometry.

We will use Eqs. 9 as the basis for a set of Hamiltonian and Lagrangian equations. We consider first the entropic equivalent to kinetic energy, i.e. ‘kinetic entropy’ (KE) T, based upon the conventional definition of kinetic energy (Appendix B, Eq. B.8b in Supplementary Information):where the additional negative sign accounts for the inverse velocity. For the three spatial directions, we therefore have:

We also define an entropic potential field V(q) as a function of hyperbolic position q (the ‘potential entropy’). However, for the present case of a double helix, Eq. 8 clearly represents a pair of plane waves travelling in space; which is analogous to the kinematic “free-particle” situation, such that there is therefore no associated entropic potential field, V = 0. The entropic Hamiltonian H(q(x3), p(x3), x3) is defined as usual as H = T + V, and (as shown in Appendix B, see Supplementary Information) is also a conserved quantity in hyperbolic space.

Using the canonical Legendre transformation, the entropic Lagrangian is given by (Eq. B.14):

such that the required canonical equations of state are obeyed: ∂L/∂x3 = −∂H/∂x3, as well as and (see Appendix B, Eqs B.15 & B.16).

Double Helix Geometry: Photons & DNA

Exertion

In analogy to the action integral (with units of J·s) we now define the exertion X (units of J/K) as the integration of the entropic Lagrangian L along the spiralling double-helical trajectory:where we note the Pythagoras relationship due to the helical geometry.

For the double-helix plane-waves description of Eq. 8, the associated entropic Lagrangian L has no entropic potential term (that is, VS = 0) since such a system is equivalent to that of a free particle. Appendix C (Supplementary Information) provides the proof that the entropic Lagrangian functional (see Eq. C.3):as employed in Eq. 12 satisfies the Euler-Lagrange equationsdemonstrating that the exertion X is at an extremum (or at least stationary) at any point along the length of the double helix since (see Appendix C, Eq. C.22). Also, Appendix D (Eq. D.3b; both Appendices are in Supplementary Information) shows that the entropic Lagrangian for a double helix can be given by ; that is, in this case L is indeed a constant (invariant with x3). Note also that the exertion X is scaled by the quantum of entropy, Boltzmann’s constant, just as the Lagrangian itself is.

Entropy

Having defined the exertion integral, Eq. 12, we can also now see that the equivalent space-trajectory integral of the entropic Hamiltonian H (see Eq. 11) yields a quantity directly proportional to the entropy:Whereas Eq. 6 describes a ‘local’ entropy s, the integrated quantity S can be considered as the ‘global’ or the overall system entropy. Eq. 14 indicates that the overall entropy S depends not only on the centroidal trajectory of the double helix axis as described by x3, but principally upon the spiralling path described by l with its radial dependency such that the entropy is a function of the full spatial extent (in all spatial dimensions) of the double helix structure. For convenience, we offset the entropic Hamiltonian H by the constant term m ln(κ0 R0) (see Appendix D in Supplementary Information, text prior to Eq. D.2a) which is an invariant for a double helical geometry – any Hamiltonian can be offset by a fixed (constant) amount to enable more convenient manipulation – such that the entropic Hamiltonian for a double helix can therefore be given as H = πκ0k; that is, each KE component (n = 1, 2) of the double helix contributes ½πκ0k. We can also exploit the Fourier (periodic) nature of S along the double helix as characterized by the parameter iκ0 to write the Fourier differential operator as:

Since the Lagrangian and Hamiltonian are inversely related (through the Legendre transformation) and the exertion integral X (Eq. 12) is at an extremum (Eq. 13b), δX = 0, then the closely connected Hamiltonian trajectory integral Eq. 14 (that is, the entropy S) must also be at an extremum, δS = 0. Given that the double helix of DNA represents a highly stable structure we infer from the Second Law that the entropy S is at a maximum; ergo the exertion X is at a minimum and the double helix topology represents a MaxEnt (most likely) trajectory in space. In summary, the overall entropy S of the double helix is given by (see Appendix D Eq. D.4):

It is clear that the entropy S is proportional to the length L of the double helix. However, in the case of a photon its proper length is actually zero relativistically, since it travels at the speed of light: L = 0, therefore S = 0.

B-DNA and P-DNA

In an extraordinary mechanical experiment, Bryant et al.[21] made a controlled transformation of B-DNA to P-DNA, where the latter is an artificial form called after Linus Pauling and discussed at length by Allemand et al.[22] (see Fig. 1). Essentially, Bryant et al. held the B-DNA molecule (of length 4.681 μm) straight in tension (45 pN), and twisted it (4800 turns, with a torque of 34 pN·nm) until it had entirely transformed into the P-DNA form (with an extension of 2.8 μm). Thus, the mechanical energy expended to turn this B-DNA molecule into a P-DNA one is 1151 aJ (126 aJ from the extension and 1025 aJ from the torque). To calculate the conformational energy changes with standard methods is computationally heavy: a recent molecular dynamics calculation by Liebl & Zacharias[23] to determine free energies actually mimicked Bryant et al.’s experiment.

Figure 1

Two forms of DNA, with dimensions. Modified from Fig. 5 of Allemand et al. Proc. Natl. Acad. Sci. USA 95, 14152–14157 (1998)[22]. Copyright (1998) National Academy of Sciences USA.

But determining the change in structural entropy (in this context equivalent to the Gibbs free energy change) is now straightforward for these holomorphic structures. Using Eq. 16, and κ0 ≡ 2π/λ0, where R0 = {1.0, 0.6} nm respectively for the B- and P- forms; λ0 = {3.32, 1.28} nm, and L = {4681, 7286} nm, we obtain from the geometric entropy (at 23 °C) the energies of the two forms of {244, 1428}aJ, yielding a change of 1184 aJ. (Note that the Type A standard uncertainty just from Bryant et al.’s torque measurement is about 70 aJ.)

It is not entirely clear which values to assign to R0, especially for the case of P-DNA, with plausible values for the latter ranging between 0.4 and 0.8 nm. In any case, it is clear that a very simple calculation using the apparatus of geometrical thermodynamics is capable of a result entirely consistent with experiment, where this result is not available without heavy computation using standard methods in physical chemistry.

To explain the stability of fullerene molecules a similar comparison can be made between this simple geometrical thermodynamics and the heavy computation required by the standard physical chemistry methods (which now have a very extensive literature)[24].

The Double-Armed Logarithmic Spiral

Figure 2 shows NGC 1566, an intermediate spiral galaxy 40 million light-years away in the constellation of Dorado (southern hemisphere) and the second brightest Seyfert galaxy known. The Milky Way is known to have a similar geometry (but of course we have no comparable image of it) and some parameters of our galaxy, including a double-armed logarithmic spiral, are overlaid on the Figure. We will show that such a double-armed spiral is holomorphic, just as is the double-helix of the photon or of DNA.

Figure 2

A spiral galaxy with an overlaid double-armed logarithmic spiral. Parameters given (see text) are of the Milky Way (MW): and MBH are the (black body) temperature and mass of the super-massive black hole at the galactic centre; is the cosmic microwave background temperature; R, L, are the MW radius, thickness and “temperature”; MMW is the MW observed virial mass; M is the mass implied by the MW entropy and . All masses in units of solar mass M◉. Underlying Hubble image of NGC 1566, taken 2nd June 2014 by NASA Goddard Space Flight Center https://www.flickr.com/photos/gsfc/14172908657/; licenced under CC BY 2.0 (https://creativecommons.org/licenses/by/2.0/).

However, we can immediately comment on the parameters of Fig. 2, which are largely determined by the mass M of the central galactic black hole whose entropy S (in SI units)

is given by the Bekenstein-Hawking equation: see for example Eq. 2.2 in the review of Bousso[25], or explicitly by Penrose (§27.10 p.716[15]) for a stationary black hole. As usual, ħ is the reduced Planck constant, k is Boltzmann’s constant, and G is the gravitational constant. But it is well known that the galactic entropy is dominated by the entropy of its central super-massive black hole (S ≈ S: see for example the discussion in Penrose 2010[2] §2.6 p.127), which has an equivalent black body temperature of (Bousso[25] Eq. 2.8):

M is given by Gillessen et al.[26] as 4.3 ± 0.4 million solar masses M◉, where this 10% uncertainty is entirely due to the uncertainty in the galactic position of the Sun: the measurement actually has a precision better than 2% (the mass of the Sun is known very accurately, to about 10−4: M◉ = 1.989 × 1030 kg). Applying this temperature to S to obtain the energy (given by the product of entropy and temperature expressed as a mass through E = mc2) we naturally recover M.

In a standard model of the Milky Way[27] (a barred galaxy) the stellar disc is modelled as distinct “thin” and “thick” discs, with the “cut-off bulge radius” of the “thick disc” (or “bulge”) given as 1.9 kpc, and the total (virial) mass within a radius of 60 kpc being 4.0 ± 0.7 × 1011 M◉. Rix & Bovy[28] explain that there is no well defined distinction between the “thin” and “thick” discs, but the characteristic “scale height” of the “thin disc” can be given approximately as 0.3 kpc. Patsis et al.[29] show that the “bar” (bulge) of the galaxy can be described in orbital dynamics terms, and Saito et al.[30] map the bulge (the “bar”) from observational data. We will assume an approximate galactic half-thickness L/2 = 1 kpc. The parsec is defined as the distance from the Sun of a star observed to have one second of arc annual stellar parallax, and is therefore relative to the diameter of Earth’s orbit (1 kpc = 3.09 × 1019 m).

Another study has concentrated on the galactic mass[31]: giving the observed stellar galactic mass and the virial mass respectively as 6.4 ± 0.6 × 1010 M◉, and 1.26 ± 0.24 × 1012 M◉. “Virial” mass includes so-called “dark matter” and is derived from the observations of stellar proper motions in large scale star surveys using the Virial Theorem (Clausius, 1870).

Consider a doubled-armed logarithmic spiral, as frequently observed for galaxies (see Fig. 2 which is a plan view projecting the 3-D object onto a plane). Here, the holomorphic functionals describing the x1 and x2 locus co-ordinates are now (in contrast to Eq. 8):

such that the instantaneous radius is for n = 1, 2 (see Eq. B.24 in Appendix B), where r is the Schwarzschild radius (the event horizon) of the central black hole. The logarithmic radial parameter Λ is given by the requirement that the galactic radius R and the Schwarzschild radius r are related logarithmically by the half-thickness L/2 (see Eq. D.15a in Appendix D, Supplementary Information): r = R exp(−ΛL/2), or

where for the Milky Way Λ = 26 kpc−1. The radius r of the central galactic black hole is determined by the black hole mass M to be 1.270 × 107 km (about 18 times the solar radius; see Eq. D.14). The coupling coefficient κ is assumed to vary similarly to the radius, that is (Eq. B.24d), where κ is the pitch at the black hole event horizon. The associated hyperbolic co-ordinates (using Eq. 9a, see Eq. B.33) are:

All quantities clearly revert to their respective double-helical quantities when the logarithmic spiral parameter Λ = 0. We find that a logarithmic spiral is associated with an entropic potential field V ≠ 0 causing a hyperbolic acceleration; indeed, as the entropic analogy to Newton’s second law of kinematics (F = mẍ), we solve the Euler-Lagrange equations (defined in hyperbolic space q) , where the final term in the equation (the entropic potential gradient) is therefore equivalent to the entropic force F. The associated entropic acceleration is given by , the minus sign being due to the inverse velocity nature of q′. The proof that the double-armed logarithmic spiral satisfies the Euler-Lagrange equations in hyperbolic space q (that is, obeys the principle of least exertion) is given in Appendix C (Eq. C.47, see Supplementary Information).

In Euclidean (x) space, we find that the entropic potential field V for the logarithmic double spiral is expressed as (see Eq. B.42 in Appendix B, Supplementary Information; K0 and K3 are dimensionless):

It is indeed interesting to note the existence of an inverse-square law (in Euclidean space) for the γ1 and γ2 directions at the heart of this entropic potential field; the entropic force varies as

that is, F ∝ x−2, with F also being proportional to the entropic mass m assumed located at the centre of the system and to be the cause of the entropic potential field. We emphasise, however, that although Eqs 22 and 23 express the entropic field in a more intuitive Euclidean form, the entropic Hamiltonian and Lagrangian equations are only correctly applied in hyperbolic space.

The general shape of the Milky Way is closely determined by these holomorphic logarithmic spirals. In particular, it can be shown that the ratio of the extremal radius R to the full-thickness L is (see Appendix D, Eq. D.17, Supplementary Information):

The radius R is rather poorly defined observationally, and the estimate L/2 = 1 kpc given above implies (from Eq. 24) R = 79 kpc, which is within the range usually given: therefore Eq. 24 has some observational support. The Milky Way in reality has a complex structure involving multiple spiral arms, a central “bulge”, and oscillating star densities reported recently[32] as persisting to much larger distances than R. None of this is considered in our zeroth order model. The present treatment should also be regarded as a static approximation neglecting the dynamic mechanisms of galactic formation and evolution. Figure 2 shows only the plan view of the model: the galactic cross-section is here modelled as a disc of essentially uniform thickness L, dimpled at its centre (that is, ignoring the “bulge” altogether).

The spiral coordinate x3 projected onto the plane in Fig. 2 is associated with an azimuthal angle θ = κx3 where the appropriate wavelength scale λ for the galaxy is given by the galactic wavenumber κ = 2π/λ. This can be calculated from the galactic structural entropy S, well approximated by (Appendix D, Eq. D.13b, Supplementary Information):where A ≡ 2π R L closely approximates the area of an ellipsoid of radii L/2 and R. Eq. 25 for the logarithmic spiral is exactly equivalent to Eq. 16 for the double-helix. Thus we get an expression for the galactic wavelength λ (see Eq. D.16 in Appendix D)where l = 1.616 × 10−35 m is the Planck length.

Eq. 25 is startling. The galactic entropy (which is almost exactly the central supermassive black hole entropy) is given “holographically” (see Appendix D, Eq. D.5 and subsequent discussion) by the surface area of the galaxy, just as the black hole entropy is determined by the surface area of the event horizon. We now postulate that just as for the black hole, a “temperature” can be defined at this holographic surface. This temperature must lie between the central supermassive black hole temperature (15 fK) and the cosmic microwave background temperature of 2.73 K.

To obtain a reasonable estimate of we note that the power radiated from a spherical black body surface of radius R and temperature is 4πσ (R2)2, where σ is the Stefan-Boltzmann constant. We therefore highlight here the appearance of the composite quantity R2 that appears as a consequence of the Stefan-Boltzmann law. With both R and increasing exponentially with distance from the galactic centre, it is clear the resulting large temperature gradient along x3 implies a large energy flow towards and into the black hole: the galaxy is not in thermal equilibrium! However, to at least maintain some thermal ‘stability’ along the γ3 axis, we might assume that 2 varies similarly to R, so that the black body power inwardly radiated from each spherical surface along the γ3 axis maintains some continuity. Relying on the isomorphism between the double-helix and the logarithmic double-spiral in hyperbolic space we therefore consider 2 to vary with exp(−Λx3) just as R and the galactic wavelength λ do (see Eq. 19 and Eq. B.24 in Appendix B):

Then, at x3 = −L/2 with Λ = 26 kpc−1, we have  =  = 6.3 nK for L = 2 kpc giving MMW = 0.94 × 1012 M◉ (and R = 79 kpc) consistent with observation. To obtain the central observed value for the virial mass of the Milky Way of MMW = 1.26 × 1012 M◉, we need L = 3.6 kpc (giving R = 142 kpc and  = 8.4 nK).

To summarise: we have shown that the structure of the galaxy for which we have detailed experimental observations (that is, the Milky Way) is consistent with a holomorphic representation in geometric algebra. In particular, we have shown that the galactic shape, aspect ratio, and structural stability (which are all highly constrained by the algebra) are consistent with observation; and we have also shown that the total galactic mass is also consistent with observation. Note that this is a simplified (“zeroth order”) analytical approximation to reality: for example, the black hole angular momentum is neglected, as are the dynamics driving the galactic evolution. Also, we have not started to consider the perturbation problem implied by deviations of the star population from the ideal logarithmic spiral; although we would anticipate that the principle of least exertion causes an entropic force to be exerted so as to maintain the MaxEnt galactic structure.

Notwithstanding the approximations, these results are very surprising, because they underline the dominant effect that the central super-massive black hole has on the galactic structure. In fact, this treatment gives the proper weight to the effect of the black hole entropy, which is certainly not hidden away behind the event horizon.

Isomorphism between Mechanics and Entropy

Table 1 shows the multiple isomorphisms that exist between kinematic and entropic quantities revealed by our treatment. There has been significant recent interest in comparable methods. Baez & Pollard[33] argue for an “analogy” between thermodynamics and quantum mechanics, giving rise to a quantity they call “quantropy” (quantum entropy, which they call “mysterious”). They also give a Table of “analogies” between statistical and quantum dynamics comparable to our Table of isomorphisms. We believe that our results confirm and extend this approach. Velazquez[34] has also tabulated some consequences of the complementarity of the Planck and Boltzmann constants. Dixit et al.[35] have reviewed the use of “Maximum-Caliber” to characterise trajectories (“world-lines”) in non-equilibrium thermodynamics (where “caliber” is a term introduced by Jaynes[36] to characterise the evolution in space-time of the ensemble of trajectories of microstates; it is proportional to our “Exertion”).

Table 1

Isomorphism between kinematic and entropic quantities.

Quantity	Kinematic (Conventional)	Entropic Equivalent	See near Eqs
Physical constant	Planck, ħ [Js]	Boltzmann, kB [JK−1]	9c, 12
Space-time co-ordinates, q	x1, x2, x3, x0 (≡ct) (Euclidean, Minkowski)	q_n=R_nln(x_nR_n) (hyperbolic, Minkowski)	1, 9a, 21
Differential operator	∇ = ∂/∂xn	∇q = ∂/∂qn	11
Wavelength, λ & wavenumber κ	λ(=2πc/ω) & k = 2π/λ	λ = helical pitch; κ = 2π/λ	5
Time-like axis & associated Fourier differential	t,ddt≡iω	x_3,ddx_3≡iκ	11, 15
Momentum, p	p=mx˙=2πℏλ	p=m_Sv^−1=m_Sq′=k_BR	9b, 10
Velocity, v	x˙_n=dx_ndtn=1,2,3	v&v^−1≡q′_n=dq_ndx_3=R_nx′_nx_n (i.e. also has inverse velocity, v−1, characteristics)	9b, 10
Mass, m	m(≡kℏc)[kg]	m_S≡iκk_B[JK^−1m^−1]	9c, 15
Kinetic term T=∫pdv	T=∫pdx˙= ½ mx˙^2	T_S=−∫pdq′=−m_Slnq′= ½ m_SInq′^2	10
Potential term, V	V=mx¨x≡mgx	V_S(q)=m_S(−q′′/q′^2)q≡m_SΓq	11, 22
Hamiltonian, H	H = T + V	H_S=T_S+V_S=∑n=13−m_Slnq′_n+V_S(q_n)	11
Lagrangian, L	L=∑n=13p_nx˙_n−H=T−V	L_S=∑n=13p_nq′_n−H_S=3m_S−H_S	11, 13a
Newton’s 2nd Law: [Euler-Lagrange formulation]	ddt∂L∂x˙_n−∂L∂x_n=0	ddx_3∂L_S∂q′_n−∂L_S∂q_n=0	13b, 21
F = ma	F = −∇V = mx¨	F_S=−∇_qV_S=−m_Sq′′/q′^2=m_SΓ	21, 23
Acceleration	x¨ = g[ms^−2]	−q′′/q′^2=Γ[m^−1]	21, 23
Action/Exertion Integral	= ∫Ldt[ℏ]	X=χ∫L_Sdx_3[k_B]	12
Variational Principle	δ = δ∫Ldt=0, Least Time/Action	δX=0,δ∫L_Sdx_3=0, Least Exertion	12, 13b
Entropy	S=∫dET	S=χ∫H_Sdx_3	14, 16, 25
Maximum Entropy	δ∫Hdt=0, Stationary Phase (Group Velocity)	δS=0,δ∫H_Sdx_3=0	15

Considering Table 1, we have already observed that the hyperbolic Minkowski space (generated through the normalising Euclidean metric, R) is the entropic analogue to the Euclidean Minkowski space of kinematics, with consequent del operators; that Boltzmann’s constant is the entropic quantum analogue to Planck’s quantum of action (also pointed out by Córdoba et al.[37]) with consequently analogous definitions for momentum; and that the helical pitch (or wavelength) implies the space-like entropic analogue of time in kinematics. As we have seen, this latter also implies holographic properties of the treatment (that is, properties of an area being fully equivalent to parameters of a volume).

Both mass and its entropic equivalent m have natural units of inverse length, but m is imaginary as a consequence of the holomorphism (Eq. 9c). The parameter Λ describing a logarithmic spiral contributes to the entropic (hyperbolic) acceleration Γ as a consequence of an entropic force, in analogy to Newton’s 2nd Law; and the double-helix can be seen as a special case (Λ = 0) of the double-armed logarithmic spiral. The Hamiltonian and Lagrangian formulations then follow equivalently for both energy and entropy, with the Exertion integral equivalent to the classical Action integral, both obeying the principle of stationary “action”.

Summary

Formal mathematics establishes tautologies which are frequently very surprising, and we have used well-established formal methods in a properly quantitative treatment of entropy, revealing that measurable (and measured) quantities from the molecular to the galactic scale can be readily calculated in a simple analytical treatment. We have considered systems of high symmetry which are amenable to our simplified analytical approach, but we expect the method to be readily generalisable to more complex systems.

The computational demands of conformational chemistry are very severe; perhaps this approach will stimulate algorithmic advances to speed the calculations for static problems, or even to address dynamic geometrical problems (like protein folding) in new ways?

We have used a “toy” model of the Milky Way, which ignores the central “bulge” and multiple arms, but a more realistic model already available would simply take a linear combination of a spherical central feature[24] and multiple double-spiral arms. The difficulty here is not in the modelling but in the choice of realistic observational data for the model parameters.

Appendices A-D