Quaternion Quantum Mechanics II: Resolving the Problems of Gravity and Imaginary Numbers

: We present a quaternion representation of quantum mechanics that allows its ontological interpretation. The correspondence between classical and quaternion quantum equations permits one to consider the universe (vacuum) as an ideal elastic solid. Elementary particles would have to be standing or soliton-like waves. Tension induced by the compression and twisting of the elastic medium would increase energy density, and as a result, generate gravity forcing and affect the wave speed. Consequently, gravity could be described by an index of refraction.


Introduction
Quaternion quantum mechanics, QQM, is ontic in the sense that it answers the central questions of the interpretation of quantum mechanics; it is directly related to being (the Cauchy elastic continuum), as well as to the basic categories of being and their relations [1].
Quantum mechanics: where we are.From its beginning, "spooky action" irritated Einstein.Present explanations assume that the collapse of the wave function has no observable consequences and is philosophically permissible.However, the unsolved problem is what happens with the mass and energy of a particle when its wave function collapses.The instantaneous jump is not expected in general relativity, and the "string theory" does not help either.
Schrödinger did not like the "probability" interpretation of the wave function and always considered the wave to be a real wave: "Let me say at the outset, that in this discourse, I am opposing not a few special statements of quantum physics held today (1950), I am opposing as it were the whole of it, I am opposing its basic views that have been shaped 25 years ago, when Max Born put forward his probability interpretation, which was accepted by almost everybody."[2] David Bohm and Basil Hiley developed an interpretation of complex quantum mechanics (complex quantum mechanics may be more precisely called operator quantum mechanics), CQM, which gives a clear and intuitive interpretation of its meaning with no need to assume a fundamental role for the human observer [3,4].This deterministic interpretation is commonly considered as basically equivalent to the Copenhagen orthodox understanding.The importance of the Bohm approach, i.e., the fact that it consistently solves the measurement problem and allows the classical description of macroscopic objects, is frequently ignored.Unfortunately, the predictive equivalence of the two theories was recently wiped out [5].
John Bell [6], despite his great impact on our understanding of CQM by his verification that nonlocal features characterize natural processes, also expressed dissatisfaction with the conceptual status of CQM [7]: "Either the wavefunction, as given by the Schrödinger equation, is not everything, or it is not right." There are widely known remarks by Richard Feynman in 1964 [8]: "It is safe to say that no one understands quantum mechanics" and Murray Gell-Mann in his lecture at the 1976 Nobel Conference [9]: "Niels Bohr brainwashed the whole generation of theorists into thinking that the job (of finding an interpretation of quantum mechanics) was done 50 years ago".
There are several concepts that contradict the "probability" interpretation of CQM and are relevant to our QQM.It is known that certain nonlinear Schrödinger (NLS) equations, in one or more space dimensions, possess space-localized solutions ψ = ψ(t, x), e.g., solitons in the one-dimensional case.From numerous attempts, we have selected a few considering such settings.Bodurov has shown that space-localized solutions happen for a large class of complex nonlinear wave equations and NLS equations [10,11].Białynicki-Birula and Mycielski have found that NLS equations admit closed-form space-localized solutions (gaussons) [12].They have also shown that "...in every electromagnetic field, sufficiently small gaussons move like classical particles".Weng's results reveal that the quaternion space is appropriate to describe the gravitational features [13].The Three Wave Hypothesis by Horodecki, which is based on de Broglie's particle-wave duality, and the assumption of covariant aether [14,15] are also consistent with QQM.Close demonstrated a description of rotational waves in an elastic solid as the spin equivalent [16][17][18].
Regardless of recent progress, it is still safe to say that there has been insignificant advances in understanding of CQM.
Quaternion quantum mechanics today.The first suggestion of quaternion quantum mechanics came from Birkhoff and von Neumann [19].Already in 1936, they mentioned that quaternion quantum mechanics has greater logical consistency than classical (complex, operator) quantum mechanics.
Yang [20] shows that it is not necessary to go beyond the three-number systems, real numbers R, complex numbers C, and quaternions Q for the representation of quantum mechanics.It should be noted that this is consistent with the Hurwitz theorem in which real numbers, complex numbers, quaternions, and octonions O , are the only normed division algebras over real numbers.In simple words, e.g., only R, C, Q and O can be used in the models where energy, is conserved.Finkelstein et al. [21] showed that quaternion calculus exists, and it is always possible to represent pure states of a system of quantum mechanics by rays in a vector space over the quaternions, but not over real and complex numbers.Recently, the global effects in quaternionic quantum field theory [22] were applied to analyze the experimental status of quaternionic quantum mechanics [23].
Adler studied downgraded quaternion Lagrangian [24] and quaternionic group representations [25,26].His idea of the trace dynamics relies on using a variational principle based on a Lagrangian method, constructed as a trace of noncommuting operator variables, making systematic use of cyclic permutation under the trace operation.Nottale used the bi-quaternion concept in suggesting an answer to the question of the origin of complex numbers and Clifford algebra in quantum mechanics [27].Recently, Gantner demonstrated the equivalence of complex and quaternionic quantum mechanics [28].
To summarize, QQM has many new features that make it a much richer theory.It is caused generally by the noncommutativity of quaternion-valued wave functions.Our quaternion Klein-Gordon [29] and Schrödinger equations [1] carry much more physical information than their complex equivalents and make QQM a much richer theory.In this work, we combine the model of the Cauchy elastic continuum with the Planck-Kleinert crystal hypothesis and derive the first and extended second-order differential equations of QQM.
The Planck-Kleinert crystal.Elastic waves play a remarkable function in physics.Thomas Young explained the polarization of light as analogue to shear waves, Navier developed his equations by adding dissipative terms to the Cauchy equation of motion, Maxwell constructed equations of electromagnetism by modelling a lattice of elastic cells, etc. [30].The Cauchy model of elastic solids was already published [31] when Maxwell considered the crystal hypothesis.In A Dynamical Theory of the Electromagnetic Field [32], Maxwell explicitly remarked on the aether: "On our theory, it (energy). . .may be described according to a very probable hypothesis, as the motion and the strain of one and the same medium (elastic aether)" and "... what if these molecules, indestructible as they are, turn out to be not substances themselves, but mere affections of some other substance?"[33], less known, if not entirely forgotten, is the remark on gravity: ". ..assumption, therefore, that gravitation arises from the action of the surrounding medium leads to the conclusion that every part of this medium possesses, when undisturbed, an enormous intrinsic energy.As I am unable to understand in what way a medium can possess such properties, I cannot go any further in this direction in searching for the cause of gravitation." Maxwell's idea of solid aether showing "enormous intrinsic energy" was unimaginable in the 19th century.
We consider aether as the Planck-Kleinert crystal, P-KC [1].The macro-properties of such a crystal are approximated by the Cauchy model of an elastic solid continuum in the quaternion representation.
The original arguments to implement classical mechanics equations in the field of wave mechanics in crystalline, granular aether were given by Kleinert [34,35].Soon after, it was shown that quantum gravity effects, when applied to a non-relativistic particle in a one-dimensional box, imply the quantization of length [36].This result was interpreted as an indication of the fundamental discreteness of space itself.Similarly, corrections to the Klein-Gordon and Dirac equations gave rise to area and volume quantizations, again indicative of the fundamentally grainy nature of space.Such an approach modifies all quantum mechanical Hamiltonians [37] and suggests that space itself is discrete, i.e., that all measurable lengths are quantized in units of a fundamental Planck length.
The building blocks of the Planck-Kleinert crystal are Planck particles, m P , that obey the laws of mass, momentum, and energy conservation.Each particle exerts short-range force at the Planck length, l P .The Kleinert concept linked with the Cauchy model of the elastic continuum was analyzed with the arbitrary assumption of the complex potential field [38].
Recently, the Cauchy theory was rigorously combined with the Helmholtz decomposition and Planck-Kleinert crystal hypothesis.The quaternion representation of the deformation, σ, in the Cauchy displacement field, u, produced the system of second-order wave (Klein-Gordon and Poisson) [39] and Schrödinger equations [39].
Dirac equation.Dirac's equation, on the one hand, is a first-order linear differential equation, and on the other hand, the iterated application of the equation yields the Klein-Gordon wave equation and consequently, the invariance under Lorentz transformation.Because of its success in explaining both the electron spin and the fine structure of atomic energy levels, the utmost importance of Dirac's discovery was evident.Several trials were made to avoid the operator method used by Dirac and to bring his equation into a form that could be interpreted in terms of normal vector analytical concepts.Cornel Lanczos made important progress and derived the first-order differential equation using quaternion algebra.
Lanczos was with Einstein in Berlin, working with the great man to whom, in 1919, he dedicated his dissertation: a quaternionic field theory of classical electrodynamics [40].Only a year after Dirac had discovered his relativistic wave equation for the electron, Lanczos published a series of papers on Dirac's equation [41,42].He showed how to derive Dirac's equation from a more fundamental system.He predicted that spin 1  2 particles should come in pairs, as well as the correct form of the wave equation of massive spin 1 particles that would be rediscovered by Proca in 1936.He foresaw the possibility of a nonlinear theory and the origin of mass exactly of the kind that would be developed almost thirty years later.In 1933 in his new derivation [43], there is a doubling in the number of solutions, from which four in Dirac's theory (two for spin and two for particle/antiparticle) increases to eight, a feature that we can today interpret as isospin.The isospin partner of the proton-the neutron-was discovered in 1932.Nobody ever thought of using Lanczos's doubling to explain the existence of isospin particles.His article, over eighty years later, still contains a number of ideas that remain at the forefront of fundamental theory.Sadly, the quaternions were non-popular and Lanczos articles were ignored by the vast majority of his contemporaries.Lanczos himself abandoned quaternions and never returned to quaternionic field theory for the rest of his life.He briefly referred to his quaternion articles of 1929 only twice [44].Over eighty years later, his papers contain ideas that remain at the frontline of fundamental theory.The whole series of Lanczos's articles is a remarkable discussion of the fundamental problems of matter, fields, and the origin of mass, most of which is still pertinent today.The first problem of physical interpretation is due to the fact that Lanczos's equation is much more general than Dirac's.The trouble with Lanczos's fundamental system (from which Dirac's equation can be derived as a special case) is that it allows for spin 1  2 solutions (such as the electron), as well as for spin 0, 1 and 3 2 .Lanczos, like anybody at the time, was completely unaware of the abundance of the elementary particles.It seems that Lanczos was also not aware of the idea that covariance with respect to spatial reversal also had to be included in order to have full relativistic invariance.
Recently, Silvis applied the quaternion formalism to the Dirac equation by making a translation of the Dirac equation as usually stated in quaternion formalism.In his approach to the Dirac equation, a two-component biquaternion and one-component biquaternion wave equations were considered [45].
In the present paper, we solve the problem from a different point of view.We use formalism, which is well adapted to the problem and is based upon the "quaternions" introduced by Hamilton.We do not try to heuristically find analogies with the classical field equations.Using quaternion algebra, we combine the Planck-Kleinert crystal hypothesis and Cauchy theory of an ideal elastic solid.We construct a Hamiltonian with the use of the Cauchy-Riemann operator, acting on quaternionic valued functions.
Quaternion calculus has never really been adopted in physics.The ideas coming from complex quantum theory remain almost completely unfamiliar to most mathematicians, mainly because of the absence of clear definitions and statements of the concepts involved.This paper attempts to close some of these gaps in communication and starts with the fundamentals of the quaternion quantum theory, with the specification of what the theory is basically about.The quaternion section can be skipped by experienced readers.

The Quaternions
Hamilton created the R 4 analog of complex numbers; his unquestionable motivation was the mechanics of solids and liquids.In Hamilton's own words [46]: "Time is said to have only one dimension, and space to have three dimensions.The mathematical quaternion partakes of both these elements; in technical language it may be said to be 'time plus space', or 'space plus time': and in this sense it has, or at least involves a reference to, four dimensions." The beauty of quaternions was immediately recognized.James Clerk Maxwell stated the following [47]: "The invention of the calculus of quaternions is a step towards the knowledge of quantities related to space which can only be compared for its importance, with the invention of triple coordinates by Descartes.The ideas of this calculus are fitted to be of the greatest in all parts of science." Quaternions can be considered a physical reality; they allow the computation of processes in continua, particularly wave mechanics.The reformulation of basic principles in terms of quaternion algebra allows one to understand classical and quantum mechanics.Our review of the basic definitions and formulas of quaternion numbers and functions is limited to those used in the paper [48].
The algebra of quaternions, Q, owns all laws of algebra with unique properties [49].The essentials here are: (1) the multiplication of quaternions is noncommutative; (2) the quaternionic displacement potential, i.e., the displacement four-potential or qpotential, which is a relativistic vector function from which the displacement field can be derived.It combines both a compression scalar potential (pressure) and a torsion vector potential (twist) into a single quaternion (four-vector); (3) the quaternionic displacement potential is the Lorentz invariant.
In the original Hamilton notation, a quaternion is regarded as the sum of a real (scalar q 0 ) and imaginary (vector q) parts: q = q 0 1 + q = [q 0 , q] ∈ Q.The following algebraical notation is useful: e 0 = 1, e 1 = i, e 2 = j, e 3 = k.Thus, an arbitrary quaternion q, i.e., q ∈ Q : = R ⊗ P, can be written in terms of its basis components, q = (q 0 , q 1 , q 2 , q 3 ) = q 0 1 + q Because it is reserved for the scalar quantity (real), for the first component, the notation l is used.The unit vector 1 behaves like the ordinary unit and can be ignored as a factor, q = q 0 + q 1 i + q 2 j + q 3 k.The remaining unit vectors, i, j, k, are usually called imaginary units.
The component-wise addition and component-wise scalar multiplication are the conventional operations.Multiplication is the fundamental operation that is defined by the multiplication of the unit vectors.The Hamilton product (multiplicative group structure) on the quaternions is defined as follows: - The real quaternion 1 is the identity element; - The real quaternions commute with all other quaternions, that is a • q = q • a, for every quaternion q and every real quaternion a; -The Hamilton product is not commutative, p Thus, the quaternions form an associative algebra over the real numbers; -Every nonzero quaternion has an inverse with respect to the Hamilton product; - The product is first given for the unit vectors, and then extended to all quaternions; The quaternions form division algebra.This means that the non-commutativity of multiplication is the only property that makes quaternions different from a real and complex numbers.The unit vectors obey the following relations: The multiplication is associative but not commutative.Instead of the simple commutative law, p • q = q • p, in quaternion algebra we have the following law: From the multiplication law (3) follows the convenient formula: where p = ∑ 3 i=0 p i e i , q = ∑ 3 i=0 q i e i ∈ R 4 ; p = ∑ 3 i=1 p i e i , q = ∑ 3 i=1 q i e i ∈ P; and • and × mean the scalar and vector, i.e., cross-products in P, respectively: A conjugate quaternion is defined as follows: where the asterisk means the following: one goes over to the "conjugate" of the quaternion, that is to say, one gives the imaginary units the opposite sign.The conjugate means one gives the vector components (the space part), q = q 1 i + q 2 j + q 3 k, the opposite sign: It is easy to see that the quantity q • q * is simply a scalar number, and all spatial components vanish.From Equations ( 2)-( 6), it can be seen that q • q * = q * • q = ∑ 3 i=0 q 2 i , and therefore the Euclidian norm can be denoted as follows: Hence, Q is a normed algebra.The multiplication given by ( 2) and ( 3) is noncommutative.The cross-product of p and q relative to the orientation determined by the ordered basis i, j, and k is as follows: The vector space R 4 with the multiplication (3) is a noncommutative algebra with unity usually denoted by Q, and it is named quaternion algebra.The commutator of two elements, p and q, is defined by the following: and can be looked at as a measure of noncommutativity.The noncommutativity of quaternion multiplication stems from the multiplication of vector quaternions.Two quaternions commute [p, q] = 0 if, and only if, their vector parts are collinear.

Representations of quaternions. The quaternions can be represented as follows:
matrices in such a way that quaternion addition and multiplication correspond to matrix addition and matrix multiplication, e.g., as 2 × 2 complex matrices and 4 × 4 real matrices [49].There is a strong relation between quaternion units and Pauli matrices; -exponent functions that have trigonometrical representation: e q = e q 0 (cos| q| + q/| q| sin| q|); -rotors, the generalization of quaternions that represents a rotation about the origin and introduces the concept of bi-vectors.Only in R 3 does the number of basis bivectors equal the number of basis vectors, and each bivector can be identified as a pseudovector.In physics and mathematics, a pseudovector (or axial vector) is a quantity that is defined as a function of some vectors or other geometric shapes, which resemble a vector and behave like a vector in many situations.Geometrically, the direction of a reflected pseudovector is opposite to its mirror image but with equal magnitude [50].
Functions of a quaternion variable.Like functions of a complex variable, functions of a quaternion variable represent useful physical models.For example, the original electric and magnetic fields described by Maxwell are functions of a quaternion variable [51].
Let Ω ⊂ R 3 be a bounded set.The so-called Q-valued functions may be written as where the functions q 0 (x), q l (x), l = 1, 2, 3 are real-valued.
Similarly, the functions q (t, x), depending on time t, may be considered.Properties such as continuity, differentiability, integrability, and so on, which are ascribed to q, have to be possessed by all the components q 0 (t, x), q l (t, x), l = 1, 2, 3.In this manner, the Banach, Hilbert, and Sobolev spaces of Q-valued functions can be defined [51], e.g., in the Hilbert space over Q, and we introduce the inner product as follows: In a similar way, the Sobolev spaces are defined: The definition of self-adjoint operators acting on these spaces is analogous as in the real and complex cases.Moreover, the theories of analytic functions, distributions, Fourier series, Lebesgue measure, Gelfand triples, Laplace transform, and many others on the vector space of Q-valued functions over Q can be defined in a standard way as in the real and complex cases with analogous properties.
Remark 1.Because it is possible to divide quaternions, they form a division algebra, and the norm makes the quaternions into a normed algebra.Hurwitz's theorem says that there are only four normed division algebras: R, C, Q and the octonions algebra.
Lagrange's four-square theorem in number theory states that every non-negative integer is the sum of four integer squares.This theorem may have applications in different areas of mathematics, e.g., quaternion algebra.
Time is scalar and has only one dimension, and space has three.Quaternion might be conceived as "time plus space", and in this sense, it has reference to four dimensions.

Quaternions and Cauchy Elastic Continuum
The displacement vector u has the following standard definition: X denotes the position vectors of material points at t = 0 and x spatial position at other times t of the point that moved, and was X at t = 0.The velocity and acceleration are defined by the following: The Cauchy theory describes the case in which an infinitesimal line element dX of the reference configuration undergoes extremely small rotations and fractional change in length in deforming to the corresponding line element dx, i.e., when ∂ u i /∂X j << 1.
The exponent function has its trigonometrical representation as follows: where σ is a Q valued function.
We also introduce a deformation four-potential as a single quaternion (four-vector) It is the relativistic function defined by the displacement field, u.As measured in a given frame of reference, and for a given gauge (gauge theory is a type of field theory in which Lagrangian is invariant under local transformations), the first component of the deformation four-potential is the compression scalar potential, and the other three components make up the twist vector potential.Note that while both the scalar and vector potential depend upon the frame, the deformation four-potential is the Lorentz covariant.

Quaternion Representation of the Cauchy Classical Theory of Elasticity
In the following section, the mechanical reactions in the real FCC crystal are assessed by means of Cauchy continuum theory, i.e., we approximate the grainy continuum by field variables.The Cauchy model of an ideal elastic continuum [52] constitutes the consistent base used here due to the following: the macroscopic phenomena are expressed in terms of field variables [53]; -from the beginning, the model was applied to study the elementary waves [54]; -the proof of the uniqueness of solutions [55] and the completeness proof are complete [56].
We follow the Planck-Kleinert crystal hypothesis [38] and consider an ideal FCC structure, in which the Poisson number ν = 0.25, l P equals the Planck length and denotes the dimension of the FCC elementary cell that consists of four interacting Planck particles, showing the Planck mass m P .The density of such continuum equals ρ P = 4m P /l 3 P .We reduce the problem, and the continuum is treated as a closed system occupying the constant volume Ω ⊂ R 3 .The Cauchy theory describes the case when any infinitesimal line element dX of the reference configuration undergoes extremely small rotations and fractional change in length in deforming the corresponding line element dx, i.e., when ∂ u i /∂X j << 1.The following considerations should also be taken into account: 1.
The continuum density, ρ P , is high, and we consider the small deformation limit only, l P ∼ = const.;thus, the density changes are negligible and ρ P = 4m P /l 3 P = const.;2.
The small deformation limit implies the invariant wave's velocities, particularly the constant transverse wave velocity in Equation ( 24): where Y is the Young modulus [57].

3.
We consider here the long evolution times, t t P , where t P is the Planck time; 4.
The quasi-stationary wave exists, which may exhibit the velocity of its mass center, υ [10,11] In such a continuum, the equation of motion relates to local acceleration due to the displacement, u, with the field variables, compression (divu), and twist (rotu): where, for the sake of simplicity, we do not consider the external fields.
From Equation (24), the energy per mass unit in the deformation field follows [57,58] where where Note that the Cauchy-Riemann operator D defined by the Equation (20), acting on the quaternion-valued potential A in an ideal elastic continuum, corresponds physically to the force F in R 3 , Equation (26).We use here the weak formulation of Helmholtz's decomposition theorem.Every deformation can be expressed by the curl-free component, u 0 , and a divergence-free component, u φ , and if u belongs to the C 3 class of functions, then u = u 0 + u φ , where rotu 0 = 0 and divu φ = 0 [59].Upon acting on Equation ( 24) by the divergence and rotation operators, we decompose it and obtain well-known transverse and longitudinal wave equations in the usual form a tt = k ∆a: The Cauchy equation of motion combined with the Helmholtz decomposition theorem in (27) leads to four second-order scalar differential equations, i.e., "quattro cluster", which implies the presence of transverse and longitudinal waves in the Cauchy elastic solid.Note that these equations remain coupled by the relation of the energy density (25); however, the more complex wave phenomena are not apparent in (27).
The Cauchy displacement field in the quaternion deformation representation shows the physical reality, the correspondence with Hamilton time-space continuum, and the complexity of wave phenomena.The Hamilton algebra Q allows the curl-free and divergencefree components that are separated in (27) to be coupled.Upon denoting σ 0 = divu 0 and φ = rotu φ we obtain and the energy density per mass unit (25) takes the following form: The decomposition u = u 0 + u φ in ( 27) and the change in variables results in four Equations (28) and allows the use of Hamilton quaternions.Namely, it implies the existence of the deformation field σ = σ 0 + φ that represents the twist and compression fields as a superposition of real (scalar compression σ 0 ) and imaginary (twist vector φ) field parts at each point where the Helmholtz decomposition implies the following constraint [58]: This representation specifies three displacement components in terms of four potential components; furthermore, the divergence of u 0 is arbitrary.It is common to choose φ with zero divergence: div φ = 0.
Adding Equations in (28), and from (30) we obtain the quaternion form of the motion equation where φ must obey constraint (31).
. u 3 , the overall energy of the deformation field, Formula (29) reaches the quaternion form The energy is conserved, so relation (33) leads to the nonlocal boundary condition for Equation (32) [29].
Remark 2. The Cauchy model combined with the Helmholtz decomposition theorem and quaternion algebra results in second-order differential Equation (32) and constraint (31).It infers the transverse, longitudinal, and complex forms of waves and shows Lorentz invariance.Equations ( 32) and ( 33) satisfy the Euler-Lagrange differential equation, i.e., satisfy the fundamental equation of the calculus of variations.

Quaternion Quantum Mechanics, the Planck-Kleinert Model
In this section, we present already published results, namely the use of the quaternion algebra for combining the Cauchy R 3 model with the Planck-Kleinert crystal hypothesis [29, 35,38,39].We regard quantum space as an analog to the Cauchy elastic solid.The properties of ideal elastic aether are presented in Table 1.The relativistic waves in a Cauchy continuum follow from the postulate of existence of the stable wave, showing the energy due to the motion and the strain of one and the same medium.Upon splitting Equation (32) into the system, the nonlinear form of the wave equation follows [1,29] where k 0 = 1/λ 0 (E).In [29], we postulated that k 2 0 = 8πm/( t P ), where = m P c 2 t P .The above second-order equation fulfills the laws of special relativity and the Lorentz invariance and fixes the fundamental problems of negative energy in the complex formulations [61].The energy computed using Formula (33) is per definition always positive due to the constraint (31).The system (34) is a hyperbolic-elliptic quaternion representation of a the wave and has solutions of the following form: The second equation in (34) is the Poisson-type equation [62], which describes the compression potential as a function of energy density in a case of the particle showing the energy E, [29].When expressed as a function of the local mass density ρ = ρ E /c 2 , where The gravitational constant equals: Remark 3. The low deformation limit allows for the simplified assumption of the constant mass density, ρ P ∼ = const, as well as the constant transverse wave velocity.Consequently, gravity in the simplified form of the Poisson equation follows.By considering the nonlinear dependence of the energy density on deformation and its impact on the wave velocity, c := c (ρ E ), one can obtain a more general form of relation (36), i.e., the relations of general relativity [13].

The Quaternion Schrödinger Equation
We treat the wave as a particle in an arbitrary volume Ω [1].The overall wave energy, E = E 0 + Q, where E 0 and Q are the ground and excess energies, respectively, follows from the energy density, E = Ω ρ E dx, where ρ E = ρ p e. The key step in deriving Schrödinger is the symmetrization of the overall energy Formula (33) i.e., e = 1 2 ˆ.
, which can be written in the equivalent symmetrical form: The overall mass of the particle, m, follows from the overall energy density The terms c 2 σ • σ * and ˆ. u • ˆ. u * oscillate and depend on the time and position.The symmetry normalizes the deformation and mass velocity with respect to the overall particle mass: The quaternionic particle mass density ψ can be called the quaternionic probability because the relation Ω ψ • ψ * dx = 1 in ( 39) is satisfied.Obviously, terms ψ = ρ P /m σ (t, x) and ψ • ψ * , vary in time.
We analyzed the evolution of the wave as in relation ( 38) in the time-invariant potential field [1], e.g., the wave in the field generated by other particles.The overall energy is now a sum of the ground and excess energy Q, We considered the low excess energies, and the impact of Q on the overall particle mass was marginal.Thus, relation (40) becomes Both the E 0 and m are constant; thus, it is enough to minimize the relation The above relation contains two unknowns: ˆ. u = ∂ û/∂t and ψ.By relating the local lattice velocity ˆ.
u to the force, the normalized Cauchy-Riemann derivative of deformation l P D σ, the velocity equals By introducing (43) and normalization (39), the relation (42) becomes the functional The functional Q[ψ], that is, the integral above, was minimized with respect to a quaternion function, such that ψ satisfies the normalization introduced in the relation (39) [1].In simple terms, we looked for a differential equation that has to be satisfied by the ψ function to extremize (here minimize) the energies allowed by (44).Given the functional (44) and the constraint in (41), the conditional extreme was found using the Lagrange coefficients method and the Du Bois Reymond variational lemma [63].The whole procedure is presented in [1].It was found that in such a case, ψ satisfies the time-invariant diffusion equation, i.e., the time-independent Schrödinger equation satisfied by the particle wave in the ground state of the energy E, that has to be satisfied together with the condition (31)

Time-Dependent Schrödinger Equation
By analogy to the complex time-dependent Schrödinger equation i ∂Ψ ∂t = − 2m ∆Ψ + 1 h V(x) Ψ, we proposed the quaternion form [1] i 46) and demonstrated that in the diagonal case, both the quaternion ( 46) and the complex timedependent Schrödinger equations are equivalent in same sense.Moreover, it was shown that by suitable natural substitution, the time-dependent Schrödinger Equation ( 46) implies the quaternion stationary Schrödinger Equation ( 45).Upon multiplying Equation (44) by −(i + j + k), it can also be expressed as follows: where Θ P = (i + j + k) /2m m 2 s −1 denotes the imaginary diffusion coefficient.When the external potential V(x) is negligible, then it can be seen that we generated a quaternion form of the diffusion equation: ∂Ψ/∂t = Θ P ∆Ψ. (48)

Second-Order Wave Systems of Equations
The Cauchy equation of motion is a sum of transverse (vector) and longitudinal (scalar) deformation waves, according to Equation (32).
It looks like Equation ( 49) contains two matchless terms.By postulating the existence of a stable wave, we already draw from (49) the Klein-Gordon Equation (29) and subsequently the formulae relating the density of the wave σ • σ * , with the density of the rate of momentum change G 0 (m)σ • σ * : where the q-potential is given by σ (σ 0 , φ 1 , φ 2 , φ 3 ) = σ 0 + φ, and G 0 (m) ∈ R is a scalar function of the particle mass.System (50) represents a boson particle showing positive energy at rest, m > 0 [29] The real meaning of System ( 50) is the postulate of the scalar coupling, G 0 (m)σ • σ * , between the longitudinal and transverse waves.Coupling is more evident upon expressing system (50) in the equivalent form: Systems ( 50) and ( 51) are identical: five equations and five unknowns, σ (σ 0 , φ 1 , φ 2 , φ 3 ) and m.If mass is known, m is the parameter in the Poisson equation above.In this section, we further develop the coupling concept and present the family of second-order quaternion wave equations.
Coupling coefficient.System (51) can be generalized to the following form: It follows that harmonic oscillator implies the following relation: The term G 0 (m)σ • σ * in ( 53) is corresponds to the density of the rate of the momentum change and can be called propagator.In the following G 0 (m) is referred to as the harmonic oscillator.It is evident that at n = 0, the coupling for boson particle follows Systems (50) and (51).For weaker coupling, n = 2, 3, . .., and the q-potentials σ n equal Upon σ n substitution into System (52), the two wave equations are evident:

The Quaternionic Oscillator
The coupling take place in the crystal elementary cell, i.e., at the Planck scale.The oscillator grants the following: That the accelerations at the Planck scale of all q-potential components are equal, . .
The function G 0 σ * • σ∈ R we call the quaternionic oscillator with the following properties: 1.
The q-potential components show common frequencies of the two harmonic oscillations: of the particle wave f and of the local process f P .The oscillations energies obey the equipartition theorem in in the P-KC unit cell; 3.
The slowest process within the particle wave controls the velocity of deformation propagation.In all the systems (55), particle wave propagation depends on the velocity of the transverse wave cl 4.
Following Cauchy, we neglect the dependence of the velocities of transverse and longitudinal waves on the energy density.This implies the constant Planck frequency: The overall mass of the particle controls the frequency of the particle wave, namely the frequency of the compression and twists: f = f (m), where m might be known or computed; 6.
The amplitudes of the q-potential, i.e., the Euclidian norms σ n = σ n • σ * n (t, x), depend on the particle geometry, e.g., its volume, shape, the velocity of the particle center, etc.These are not discussed in this work; 7.
We also neglect the energy of the external field, E QF , that is generated by the particle itself, i.e., we neglect the energy of the force fields generated by Poisson equations, Equations ( 51) and ( 55); 8.
The duration of the particle wave cycle T = 1/ f exceeds the Planck cycle duration (the Planck time is the least analyzed period of time and can be considered as the time unit in QQM) by many orders of magnitude: f P f .We consider processes at t > T t P , i.e., stable particles only, and do not analyze processes at T > t > t P , e.g., the collapse or interactions between particles.
The original Cauchy motion equation in vector form, in R 3 , and its equivalent quaternion form, in R 4 , show that the cycle duration cannot be affected by space dimensions: R 4 vs.R 3 .Accordingly, we simplify the oscillator problem by analyzing the displacement u in R 1 and applying the results to the q-potential σ in R 4 .
The displacement u (t, x) of the Planck mass m P in the unit cell is a function of two simultaneous harmonic processes, i.e., f and f P , namely the displacement due to the Planck cycle u P (decisive for the propagation velocity) and due to the particle cycle u λ .The same is valid in R 1 : u (t, x) = f (u λ , u P ).The simultaneous displacements can be understood as the displacement during the Planck cycle occurring simultaneously with the displacement due to the particle wave cycle.In both cycles, we assume the harmonic approximation, which implies a simple relation between the characteristic velocities: the magnitude .u P and the average velocity .u P : During the Planck cycle, the average velocity, .
u P = c, and the magnitude of the displacement velocity are related by Equation ( 56): .
During each Planck cycle, the velocity changes four times in the range [−1/2π c, 1/2π c], as shown in Figure 1.Thus, the sum of velocity changes at the Planck distance equals the following: Upon dividing the sum of the changes by the Planck length, we obtain the rescaled Planck frequency as follows: The momentum change during the particle wave cycle follows the same schema, as displayed in Figure 1.The average velocity of the particle wave .u λ = f λ and the magnitude of the particle wave velocity follow relation (56): . u λ = 1/2π f λ.The sum of velocity changes at the wavelength solely due to the particle cycle equals the following: Which, upon dividing by the wavelength λ, results in the rescaled frequency solely due the particle cycle: During each Planck cycle, the velocity changes four times in the range      u P cos α (t), where α (t) = 1 ⁄2π f P t and .u P = 1 ⁄2πc is the magnitude of the displacement rate, according to Equation ( 56).
The Planck and particle cycles are simultaneous, and the average displacement acceleration is a product as follows: ..
By noting that σ 0 = divu 0 = lim ∆x→l P ∆u 0 /∆x, we assume that relation (62) holds at the Planck scale for deformation.Thus, the average acceleration of the scalar part σ 0 of the q-potential σ equals the following: The common frequency postulate allows the relation (63) for the all q-potential components to be extended: σ 0 , φ 1 , φ 2 , φ 3 in R 4 : The average acceleration of the q-potential σ equals the estimated average acceleration of changes of the quaternionic oscillator in the particle wave: where f is an unknown particle frequency that may be postulated or computed.
The particle wave frequency f = f (m 0 ) follows from the R 1 schema in Figure 1.The sum of moments of all the Planck masses forming the particle wave (at the arbitrary time and solely due to the particle wave) equals the momentum of particle m 0 itself.To simplify, we estimate the average moment of the arbitrary single Planck mass m P during the particle cycle T = f −1 .The cycle implies that Planck mass returns to its initial conditions: u P (t) = u P (t + T) and .u P (t) = .u P (t + T).The overall distance in which the moment of the mass m P changes equals 2πl P .Consequently, the average momentum of a Planck mass m P is given as follows: p (m P ) = m P 2π l P T = 2π m P l P f .(66) The momentum of the particle m 0 is due to the particle propagation velocity: Both moments ( 66) and ( 67) must be equal, and the frequency of the particle wave becomes: Combining relations (68) and f P = 1/t P , the total power of the quaternionic oscillator equals: Upon replacing m 0 c 2 = E 0 in (68), the Planck-Einstein relation follows: E 0 = h f where h = 2π .

The First-Order Wave Equation in P-KC
Complex quantum mechanics is based on complex number algebra, the matrices, and the matrix algebra [64].Canonical quantization starts from classical mechanics and assumes that the point particle is described by a "probabilistic wave function".Dirac applied complex combinations of displacements and velocities in the linear problem of secondary quantization [65,66].He replaced second-order Klein-Gordon equation by an array of first-order equations, and as a result, separated the different time scales.Dirac immediately recognized the problem of medium for the transmission of waves: "It is necessary to set up an action principle and to get a Hamiltonian formulation of the equations suitable for quantization purposes, and for this the aether velocity is required" [67].We follow a different path and advance quaternion quantum mechanics using simple heuristic considerations based on the concept of the medium as a solid "aether", i.e., we consider the aether as the Planck-Kleinert crystal.We base it on the following: (1) Quaternion representation of the P-KC dynamics and canonical quantization (canonical quantization in the sense that we develop quantum mechanics from quaternion representation of classical mechanics) that yielded the Klein-Gordon Equation (29): where σ = σ 0 + φ ∈ Q is the deformation q-potential; (2) Postulation of the time-invariant harmonic oscillator at the Planck scale G 0 (m) oper- ating at the Planck frequency f P = 1/t P = const (see Section 5.2); (3) Quaternion representation of the deformations (37) and ( 39) velocities (43) that yielded the Schrödinger equation.
From the second-order "electron wave equation" to the first-order equation.Schema of the secondary quantization: Step 1: Second order equation in R 4 , variable : The elastic compression in R 4 : the quaternionic form of the Riccati equation Step 1. From the q-potential σ to the vectorial momentum ˆ. u Using relations ( 54) and the q-potential of the Schrodinger equation, i.e., σ = σ 0 + φ = √ 3σ 0 + φ [1], we generate the following system: The Cauchy wave equations in R 4 & harmonic oscillator : Thus, System (55) for the deformation potential σ = σ 0 + φ = √ 3σ 0 + φ becomes: where c denotes the transverse wave velocity.The relation (40) for the total energy of the free particle and the relation between the mass velocity and the Cauchy-Riemann derivative, D σ = − m ˆ.
u, hint at the displacement velocity (i.e., the normalized momentum ˆ. u.) as the alternative variable: The particle is stable, and its wave is at a steady state.The local changes of deformation potential σ are only due to the wave propagation within the volume occupied by the particle.We know the propagation velocity c; thus, the time derivative of the potential σ in (72) we express as follows: ∂ σ ∂t The first term on the right-hand side is the propagation velocity c, and the term in the bracket is the Cauchy-Riemann derivative.We have already proved that, taking into account the obligatory restriction of the Cauchy continuum, div φ = 0, at any time t, the spatial distribution of deformation potential obeys the following set [1]: The next step might be called secondary quantization.The ∂x/∂t term is the known wave propagation velocity: ∂x ∂t and relations (76) imply the following The relation between the deformation and kinetic energies in P-KC, relations (37) and ( 40), imply the Introduce [1] Introducing the relations (77), ( 78) and G 0 = 8π mc 2 t P in system (72) results in and finally, Relation ( 39) ψ (t, x) = ρ P /m ˆ. u, implies that by multiplying the particle wave equation in the system (80), by ρ P /m, it will be expressed as a function of probability, 1 c The system (80) requires the σ 0 time dependence.This dependence results from the continuity equation presented in the next section.Only upon neglecting the compression In this section, we combine the equations of the mass continuity, CE, and of the state to obtain the time dependence of the scalar potential σ 0 (compression) in an ideal elastic solid (it can be considered an quaternionic equivalent of the Riccati equation that usually is written as follows: ∂/∂t + .u • ∇ σ + σ • σ + σ p = 0).Obviously, a simplified assumption of the constant Planck density is disobeyed (consequently, all the waves' velocities depend on the displacement).The mass density in an ideal elastic solid in R 3 is affected only by displacements u: The mass density can be defined as follows: where m (Ω) denotes the time-invariant mass contained in the deformed volume Ω (t, x).
In an ideal elastic medium, it depends exclusively on the irrotational deformation (divu = divu 0 = σ 0 ).One can relate the deformation and density in a case of sphere in an ideal elastic continuum.In such a case, the mass in a sphere does not depend on its radius, m(Ω[r(t, x)]) = m = const, and the density is affected only by the radius: The quaternion form of the first-order wave equation presented in this work allows one to obtain an insight into the Dirac equation and therefore spin 1 /2.Spin 1 /2 fermions are the cause of the Pauli exclusion principle, and therefore it is important to understand the physical meaning of spin 1 /2 in the Planck-Kleinert model.In order to visualize this concept, a simple interactive simulation of a periodically twisting and compressing 3D grid illustrating spin 1 /2 in an elastic solid for two particles is presented [68,69].

Conclusions
The presented quaternion representation of quantum mechanics allows its ontological interpretation.In simple words, the correspondence between classical and quaternion quantum equations permits one to consider the universe (vacuum) as an ideal elastic solid.Elementary particles would have to be standing or soliton-like waves.Tension induced by the compression and twisting of the elastic medium would increase energy density and consequently, carry out the following: generate a gravity forcing; -affect the wave speed.Consequently, the gravity could be described by an index of refraction [69].
The present theory was created by combining the Cauchy model of the elastic continuum with the Planck-Kleinert crystal hypothesis.The quaternion-imaginary Lagrangian, the quaternion motion equation, and the quaternionic oscillator allowed the following to be derived: -A Schrödinger equation from the functional integral, which identifies the quaternionimaginary quantum Hamiltonian; - The second-order wave equation system describing both the bosons and the gravity in terms of quaternionic Poisson equation; - The first-order quaternionic wave equation system; - The family of the second-order wave equation systems describing both the particles and the generated quaternionic force fields (four-potential); - The Planck constants, = m P c 2 t P = 1.0545727 × 10 −34 , and gravity constant, G = l 3 P / t 2 P m P = 6.674082 × 10 −11 ; -The quaternionic continuity equation in an ideal elastic solid.
The meaning of the particle mass center (particle ≡ wave) is assigned here to "spacelocalized" and is used in the sense given by the Bodurov definition [70].
Quaternion quantum mechanics has many new features that make it a much richer theory.Its great potential is visible, e.g., in the following: - The comparison of the first-order wave equations in quaternion formulation, Equation ( 81 A simple interactive simulation of a periodically twisting and compressing 3D grid illustrating spin 1 2 in an elastic solid for two particles is presented [68].

.u
= ∂u/∂t.Equation(24) and relation(25) obey the Euler-Lagrange relation are sufficient to describe deformation in the ideal elastic contin- uum.The Helmholtz theorem allows the use of quaternion algebra.The strong formulation of the decomposition theorem introduces the four-potential A:

Figure 1 .
Figure 1.The Planck cycle.The velocity of the displacement in the Planck volume of the ideal elastic continuum visualized as 1 displacement rate, according to Equation(56).Upon dividing the sum of the changes by the Planck length, we obtain the rescaled Planck frequency as follows: during the particle wave cycle follows the same schema, as displayed in Figure1.The average velocity of the particle wave uf   = and the magnitude of the particle wave velocity follow relation(56):12 uf   =.The sum of velocity changes at the wavelength solely due to the particle cycle equals the following:

Figure 1 .
Figure 1.The Planck cycle.The velocity of the displacement in the Planck volume of the ideal elastic continuum visualized as R 1 projection of the circular motion: .u P (t) =.

Table 1 .
The physical constants of the Planck-Kleinert crystal (fcc ideal isotropic crystal).