Huygens triviality of the time-independent Schr\"odinger equation. Applications to atomic and high energy physics
Arkady L.Kholodenko (1), Louis H. Kauffman (2) ((1) Clemson University, (2) University of Illinois at Chicago)

TL;DR
This paper demonstrates that the time-independent Schrödinger equation is mathematically equivalent to the wave equation through Huygens triviality, clarifying longstanding confusions about Huygens' principle in quantum mechanics and exploring diverse applications in physics.
Contribution
It establishes the Schrödinger equation's Huygens triviality, linking it to the wave equation, and applies this insight to various areas in atomic, high energy physics, and cosmology.
Findings
Schrödinger equation belongs to Huygens' equivalence class with the wave equation.
Clarification of Huygens' principle role in quantum mechanics.
New applications of geometric methods in physics.
Abstract
Huygens triviality-a concept invented by Jacques Hadamard-describes an equivalence class connecting those 2nd order partial differential equations which are transformable into the wave equation. In this work it is demonstrated, that the Schr\"odinger equation with the time-independent Hamiltonian belongs to such an equivalence class. The wave equation is the equation for which Huygens' principle (HP) holds. The HP was a subject of confusion in both physics and mathematics literature for a long time. Not surprisingly, the role of this principle was obscured from the beginnings of quantum mechanics causing some theoretical and experimental misunderstandings. The purpose of this work is to bring the full clarity into this topic. By doing so, we obtained a large amount of new results related to uses of Lie sphere geometry, of twistors, of Dupin cyclides, of null electromagnetic fields, of…
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Huygens triviality of the time-independent Schrödinger equation.
Applications to atomic and high energy physics
Arkady L.Kholodenko1, Louis H. Kauffman2
1375 H.L.Hunter Laboratories, Clemson University, Clemson,SC 29634-0973,USA
2Department of Mathematics, Statistics and Computer Science,University of
Illinois at Chicago, 851 South Morgan Street, Chicago, IL, 60607-7045
Abstract
Huygens triviality-a concept invented by Jacques Hadamard-describes an equivalence class connecting those 2nd order partial differential equations which are transformable into the wave equation. In this work it is demonstrated, that the Schrödinger equation with the time-independent Hamiltonian belongs to such an equivalence class. The wave equation is the equation for which Huygens’ principle (HP) holds. The HP was a subject of confusion in both physics and mathematics literature for a long time. Not surprisingly, the role of this principle was obscured from the beginnings of quantum mechanics causing some theoretical and experimental misunderstandings. The purpose of this work is to bring the full clarity into this topic. By doing so, we obtained a large amount of new results related to uses of Lie sphere geometry, of twistors, of Dupin cyclides, of null electromagnetic fields, of AdS-CFT correspondence, of Penrose limits, of geometric algebra, etc. in physical problems ranging from the atomic to high energy physics and cosmology.
Keywords :
Huygens principle
Time-independent Schrödinger equation
Lie sphere geometry
Dupin cyclides
1. Motivation and background
It is well documented that Schrödinger’s equation in its known form emerged only in the 4th installment of Schrödinger’s papers on quantum mechanics -all published in 1926.The Huygens principle was introduced in the 2nd installment. In the 1st installment the Hamilton-Jacobi (H-J) equation was used as point of departure. This equation was then used as an input in his variational derivation of the stationary Schrödinger’s equation. Although the variational way of obtaining this equation was subsequently endorsed by Courant and Hilbert [2], neither Schrödinger himself nor the rest of physics community were using this (variational) way for obtaining the stationary Schrödinger’s equation. Instead, in the 2nd installment, being guided by results of De Broglie (his PhD thesis was completed in 1924), Schrödinger presented the following arguments. ”Hamilton’s variational principle can be shown to correspond to the Fermat principle for wave propagation in configuration (q)-space, and the H-J equation expresses Huygens’ principle for this wave propagation.” In the same paper he also writes ”The H-J equation corresponds to Huygens’ principle (in its old simple form, not in the form due to Kirchhoff).” Kirchhoff’s way of dealing with Huygens’s principle is discussed, for example, in the book by Baker and Copson [3]. Thus, for Schrödinger the Huygens principle is synonymous with the H-J equation. Later on, in 1948, Feynman [4] made the following comment about Huygens’ principle. In section 7, page 377, of [4] we read that the equation
[TABLE]
”is easily interpreted physically as the expression of Huygens’ principle.” Here the index represents time ticks, is the wave function, and are known constants, is the classical action between space-time points and . Further down he writes: ” Actually Huygens’ principle is not correct in optics. It is replaced by Kirchhoff’s modification which requires that both the amplitude and its derivative must be known on the adjacent surface. This is a consequence of the fact that the wave equation in optics is second order in the time. The wave equation of quantum mechanics is first order in time; therefore, Huygens’ principle is correct for matter waves, action replacing time.” From these quotations the question emerges: How the Huygens principle (HP) by Schrödinger is related to that by Feynman? The first attempt to provide mathematically satisfactory answer to this question was made by Gutzwiller [5]. His work is incomplete though as he acknowledges himself. Thus, this paper (and those which will follow in the sequel) is aimed at providing missing details. By doing so, it will become obvious why this topic is still of such profound importance. To demonstrate the importance, it is sufficient to recall some facts from Feynman’s lectures on physics [6] as well as from his book on path integrals [7]. They both begin with the discussion of the two-slit experiment. For the light this experiment was set up originally by Young in 1801 [8]. Its explanation involves uses of HP [4]. For electrons it was performed initially in 1961[9]. After this, use of other heavier particles, including C60 fullerens, showed the same pattern as that observed by Young for light [10,11]. In view of these experimental results, it is appropriate to make some comments on Gutzwiller’s paper [5]. On page 54 we find the following statement: ”The wave equation
[TABLE]
implies Huygens’s principle which is the relativistic version of Feynman’s path integral111Perhaps Gutzviller had in mind that iteration of Eq.(1.1) (which is the epitome of Huygens’ principle for Feynman) is leading to the nonrelativistic path integral satisfying standard Schrödinger’s equation. This result, apparently, can be generalized to the relativistic case where Huygens’ principle in its conventional form holds. Such path integral is expected to satisfy the wave Eq.(1.2). This program was left unfulfilled by Gutzwiller [5]., valid in its usual form for Schrödinger’s equation
[TABLE]
but neither Huygens’ principle nor the path -integral applies to the sationary waves which are the solutions of Eq.s (1.4) and (1.5)”. In numeration of this paper these equations are respectively given by
[TABLE]
and
[TABLE]
By making such a claim Gutzwiller contradicts Schrödinger who clearly stated that the H-J equation is mathematically restated Huygens’ principle. Eq.(1.4) is directly obtainable from the H-J equation as shown in Schrödinger’s 1st installment on quantum mechanics [1]. Below, in this paper we shall prove that Schrödinger’s definition of Huygens’ principle is indeed correct in a rigorous mathematical sense. If this is so, then if electrons and photons222And other heavier particles produce the same interference patterns in the double slit experiment, why then formalism developed in optics [8] cannot be applied to electrons and heavier particles? What makes use of Born’s probabilistic interpretation of the wave function in quantum mechanics superior to that used in optics for description of the two- slit experiment ? This issue was carefully investigated by David Bohm in his classical monograph on quantum mechanics [12], pages 97-98. He found that differences do exist but just in few places. The updated comparison was recently made by Sanz and Miret-Artes in their book [13], chapters 4 and 7. From these chapters it follows that all objections made by Bohm in the remaining few places can be removed. The results presented in [13] along with references on which these results are based are incomplete to a some extent. In [14] we eliminated this deficiency so that it should be read alongside with reading of this manuscript. Since nowadays the sophisticated quantum mechanical experiments are mainly done optically [15,16], ref. [14] supplies helpful additional guidelines for understanding of these experiments. In view of this, it makes sense to claim that our understanding of all subtleties of quantum mechanics is contingent upon our understanding of optics where Huygens’ principle (HP) is playing a very prominent role. Although the essence of HP is summarized in Definition 3.1. in section 3, details are essential. They are presented in sections 4,5 and Appendix B. In addition, the HP is linked with the conformal invariance of the Huygens-trivial equivalence class of equations. The notion of Huygens’ triviality was formulated by Hadamard (details are presented in sections 4 and 5 and in Appendix B). The equivalence class is made of all 2nd order partial differential equations (PDE’s) which can be transformed in a prescribed way (described in the text) into standard wave equation obeying the HP. This equation is conformally invariant. Its conformal group is It is the largest symmetry group leaving invariant free Maxwell’s equations or, which is equivalent, the massless Klein-Gordon (or wave) equation []. is the conformal group of Minkowski spacetime []. It contains 4 translations, 6 rotations, 1 dilatation and 4 inversions. Hadamard conjectured that Huygens’ -trivial equivalence class which he defined is the only one possible. Results of Appendix B demonstrate that this is not the case. Subsequent studies revealed that there are conformally invariant PDE operators which cannot be transformed into D’Alembertian using the transformation rules set up by Hadamard. But they do respect the HP nevertheless! Unexpectedly, their form is determined by the nature of cosmic gravitational plane-wave background. This is explained in the Appendix B.
Our work is made of 7 sections and 4 appendices. The role of appendices is not just auxiliary. Each of them serves as a nucleus for some further work. Therefore, they cannot be omitted upon first reading. In section 2, following the original but not widely known ideas by Schrödinger, we discuss the Schrödinger-style derivation of the uncertainty relations based on relativistic arguments. These arguments are such that they allow to restore the stationary Schrödinger equation. In section 3 we take into account that inclusion of time-dependence, that is replacement of the hyperbolic (wave) equation from which the stationary Shrödinger equation was derived, by the parabolic (truly Schrödinger) equation had occurred only in the 4th installment of Schrödinger’s papers on quantum mechanics [1]. Such an inclusion was associated with a lot of difficulties for Schrödinger. The biggest of this difficulties was the apparent departure from ideas by De Broglie on the physical nature of waves of matter. Clearly, the relativistic considerations leading to Schrödinger’s formulation of the uncertainty principle also need to be sacrificed (e.g. read the additional comments in the Appendix A). As result, the role of HP in the formalism of quantum mechanics had become obscured resulting in some erroneous statements. Subsequent theoretical and experimental works cited in section 3 contain erroneous claims requiring modification of the superposition principle of quantum mechanics-one of the pillars (e.g. the double slit experiment) of quantum mechanics. These circumstances caused us to write a very detailed section 4 describing the mathematical aspects of the HP. Additional information is contained in the Appendix B. In section 4 and Appendix C we find new solutions of the Schrödinger equation-Dupin cyclides. These are having both micro and macro (cosmological) importance. The micro importance ultimately originates from the seminal work of Madelung on hydrodynamical formulation of quantum mechanics which, to our knowledge, was never reproduced in its entirety in English. The cosmological significance of these solutions originates from the remarkable work by Roger Penrose discussed in the Appendix B and also briefly in the Appendix D. Thanks to work by Ward (e.g. read Appendix B), induced by work of Penrose, instead of cosmic microwave background we may think about the cosmic gravitational plane-wave background. It is this background which is the real cause for the Schrödinger equation to exist in its known form. Indeed, Dupin cyclides emerge from analysis of Madelung equations of quantum mechanics, on one hand and, of Huygens- nontrivial (as compared to Huygens-trivial) conformally invariant class of PDE’s which should exist in the gravitational plane-wave background, on another. Use of conformal transformations (not to be confused with transformations defined by Hadamard discussed in section 4) relates Hugens-nontrivial and trivial equations to each other.
To strengthen these conclusions, we also developed alternative paths for reaching the same goals. These are coming from the detailed study of group-theoretical properties of the Schrödinger equation initiated some time ago by Niederer []. By analogy with Newton’s mechanics, where Newton’s second law is invariant with respect to Galileo-type transformations, one would expect the same for the non relativistic Schrödinger’s equation. This happens not to be the case, however. This fact somewhat complicates the recovery of Newton’s equations from quantum mechanics via Ehrenfest theorem []. Studies by many authors had firmly established that the group of symmetries of, say, hydrogen atom is []. This happens to be symmetry of all atoms of periodic system of elements as well as of molecules (at least diatomic) []. But this symmetry group is the largest symmetry group which leaves invariant free Maxwell’s equations or, which is equivalent, the massless Klein-Gordon equation [].The de Sitter and anti-de Sitter groups and are subgroups of These subgroups are not isometry groups of the Minkowski spacetime though. Details about these groups and their Lie algebras are summarized in []. These groups are also discussed in section 7 and Appendix D. It is fundamentally important that static Einsteinian spcetimes are conformally Minkowskian. Therefore they are invariant under the action of the conformal group []. The group is a typical representative of the so called dynamical symmetry groups and their spectrum generated algebras [].The relationship between geometry of the underlying spacetimes and quantum mechanics was noticed and developed to a some extent already by Schrödinger. In his book [] he made an attempt to extend his apparatus of quantum mechanics for flat spacetimes to curved spacetimes. He was interested in finding those spacetimes which permit quantum mechanics to exist. Since the group supports electromagnetic waves, that is dynamics of massless photons, the same should be true for all massless particles, e.g. neutrino, graviton, gluon, etc. Penrose’s twistor theory was developed initially for description of dynamics of such types of particles []. However, it was not used in the atomic and molecular physics because atoms and molecules are massive. This obstacle is possible to by pass as explained in section 6. Thus, in view of this, in this work we initiate use of twistors in the atomic physics.
The Huygens principle is linked with the conformal invariance. Details are provided in section 4 and the Appendix B. The conformal invariance is very nontrivially associated with the Lie sphere geometry. Details on this geometry are given in section 7. By design, this geometry transforms circles into circles and spheres into spheres. Since points are spheres of zero radius and hyperplanes are (hyper)spheres of infinite radius the Dupin cyclides are covering all surfaces made of spheres whose centres are moving along some prescribed curve c(t) parametrized by time so that the radii of spheres also change with time. By design, thus made surfaces are invariants of the Lie sphere geometry. The simplest possible examples are: spheres, planes, cones and cylinders. Use of Möbius transformations (this is the subgroup of the Lie sphere transformations group) transforms these simple objects into the whole variety of Dupin cyclides which stay invariant under the Lie sphere transformations. In the same section 4 we discovered an unusual property of the stationary Schrödinger equation which makes this equation to be treatable gauge-theoretically by employing deep mathematical results of Andreas Floer. This result came as a by product of our efforts to adopt the progressive wave solution method by Friedlander [70] to the stationary Schrödinger equation. His method was applied initially to the conformally invariant D’Alembert wave equation resulting in Dupin cyclides as solutions. In section 4 using mathematically formulated Huygens’ principle we demonstrated that the Dupin cyclides are also valid solutions of the stationary Schrödinger equation. The same result was reobtained in the Appendix C by different methods. In doing so we were motivated by the results obtained in sections 5 and 6. Section 5 begins with removing the mass parameter from the relativistic Klein-Gordon (K-G) equation using Hadamard -type transformations. These are making the K-G equation Huygens-trivial. Next, we notice that every spinor component originating from the Dirac equation is obeying the K-G equation. This means that the massive Dirac equation can also be made Huygens-trivial. Next, in this and in section 6, we discuss all relativistically invariant massive equations with integer or half integer spins and conclude that they are also Huygens-trivial. We continue our study of Huygens triviality in section 6 where we demonstrate Huygens triviality of the stationary Schrödinger equation being influenced by the seminal paper by Vladimir Fock published in 1935 [101],[102]. His goal was not to demonstrate Huygens triviality of the Schrödinger equation (say, for the hydrogen atom). Nevertheless, he came very close to this task. Huygens triviality of the stationary Schrödinger equation enabled us to use twistor methods for solving hydrogen atom problem in section 6 while in section 7 following ideas of Sophus Lee and Felix Klein we demonstrate the isomorphism between the twistor methods and that of the Lie sphere geometry. In the same section we discuss some physical applications such as the interrelation between the Dupin cyclides and torus knots. Such unexpected utility of the Lie sphere geometry we pushed further. In section 7 and the Appendix D, we demonstrated its major role in establishing the AdS-CFT correspondence and in establishing the Penrose limits.
**2. Schrödinger-style derivation of uncertainty relations leading to **
the Schrödinger equation
2.1. Shrödinger-style derivation of uncertainty relations. Role of relativistic arguments
According to Feynman [6], the formalism of quantum mechanics is lying on two pillars: a) Heisenberg’s uncertainty principle and b) the two-slit experiments with photons, electrons, etc. Being guided by these ideas, we found it very illuminating to develop these ideas from scratch based on some considerably lesser known writings by Schrödinger. In an obscure publication [26 he sketched a derivation of the uncertainty relations not found in any other textbooks on quantum mechanics. We begin, however, with equally interesting Schrödinger’s remarks on page 50 of the same reference. ” Now the special technique by which classical mechanics dodges the awkward fact of indeterminateness (the fact that equal initial conditions are followed by different consequences) consists in including initial velocity within the initial conditions….the initial velocity is taken as forming part333In mechanics (our observation) of the initial condition at any given moment. Velocity, after all, is defined as a differential quotient with respect to time
[TABLE]
This definition refers to two moments of time and not to the state at one moment… It may be that the mathematical apparatus devised by Newton is inadequately adapted to nature; and the modern claim that the concept of velocity becomes meaningless for a precisely defined position in space points strongly in that direction.” These arguments by Schrödinger were taken into consideration seriously only recently, e.g. in 2002 in the book ”Quantum Calculus” . Obviously, his arguments serve as precursors for arguments leading to the uncertainty relations. Surprisingly, no such relations can be found in ”Quantum Calculus” [] while in the book by Schrödinger [], on page 126, we find the following comments : ” …according to the fundamental equation of the Quantum theory:444Here is energy, is Planck’s constant and is frequency.
[TABLE]
To measure frequency we need a certain time. Let us think of the primitive procedure of counting vibrations within a definite time Then,
[TABLE]
but manifestly with a possible error of because the process of counting necessarily results in giving the whole number, which is subject to an error of This entails a possible error with respect to energy of hence Now, relativistically555Emphasis is ours, the energy is the fourth component of the energy-momentum vector….Therefore the uncertainty relation can be transferred to the other components as well, for example:
[TABLE]
where is the momentum in the -direction.” The same result Schrödinger obtains differently on page 129. There, he writes:
”We need but replace the particle by a wave-group and let the wave -length and the momentum have a relation666This is clear from the dispersion relation for the light: or Here is speed of light. Multiplying both sides of this dispersion relation by we obtain: where is the same as in Eq.(2.5). Multiplying both sides of Eq.(2.5) by and taking into account that we obtain: Therefore, it follows that the De Broglie relation is also valid for light, that is for the massless particle.
[TABLE]
In order to build such a group a certain interval is required. Let be the length of the group, then…the ratio must be allowed to wary by one unit…Thus:
[TABLE]
Multiply by , then is the uncertainty And so,
From these extensive quotations from writings by Schrödinger it follows that:
a) the relationship implies the relationship
b) the relationship is obtained with the assumption that the
motion is periodic;
c) the De Broglie wavy relation (2.5) is consistent with the relation (2.4) obtained
with the help of relativistic arguments.
It happens, that the stationary Schrödinger’s equation can be obtained based on these three observations only. This can be seen from reading of Schrödinger’s Physical Review paper.[28 In it, Schrödinger acknowledges that his theory came as result of elaboration on works by De Broglie. It is helpful to reproduce Schrödinger’s arguments in order to emphasize some essential elements which were overlooked subsequently in physics literature. Specifically, the connection between the nonrelativistic Schrödinger equation and the relativistic mechanics sketched by Schrödinger has not found its well deserved place in physics literature to our knowledge.
2.2. From uncertainty relations to stationary Schrödinger’s equation
Relations a) and c) of previous subsection are consistent with the relativistically invariant scalar wave equation
[TABLE]
where the D’Alembertian is given by
[TABLE]
Here, as before, is the speed of light in the vacuum. The simplest solution of (2.6a) is given by the plane wave []. When it is substituted into Eq.(2.6a), it leads to the dispersion relation
[TABLE]
implying . However, and Therefore, we obtain: This result was mentioned in the previous subsection. Let furthermore and then Since we obtain or in accord with Schrödinger [. The question arises: Are there solutions of (2.6a) other that plane wave(s) or their linear combinations This issue was studied in detail in Evidently, the plane waves are valid solutions in the vacuum. But there are other types, e.g. progressing waves, to be discussed below in section 4.2., etc. If the waves are propagating in some (innhomogeneous) medium with refractive index , the results from optics require us to replace by . In such a case the simple plane wave solution is no longer suitable and should be replaced by 777More on this, please, read in Appendix C.. Clearly, Substitution of this ansatz into Eq.(2.6a) leads to the equation
[TABLE]
However, Therefore, since (in view of Eq.(2.5)) , we obtain : . Using the De Broglie relation and keeping in mind that the total energy of the dynamical system made of a particle moving in the potential is
[TABLE]
we obtain:
[TABLE]
With help of this result, Eq.(2.8) acquires the familiar form of stationary Schrödinger’s Eq.(1.4). In such a form this equation was used by Schödinger in his first three communications on quantum mechanics [] (not to be confused with the his paper published in Physical Review [ Only in the 4th paper he found a way of introducing the time dependence correctly. This topic is discussed in the next section.
**3. Inclusion of time-dependence into Schrödinger’s equation **
**by Schrödinger and its impact on the latest experimental **
**and theoretical works **
Inclusion of time-dependence into Schrödinger’s equation was made only in the 4th installment of Schrödinger’s series of papers on quantum mechanics [ There is a good reason why this was done so late. It will be explained at length in the companion paper888In the meantime, please, read [30].. In this paper we employ physical arguments, beginning with those by Schrödinger. On page 103 of [ we find the following comment :
” Thus, when we designated equation (1) or ()999Our Eq.(1.4). on various occasions as ”the wave equation”, we were really wrong101010Emphasis is ours…”. To fix this problem (that is to make things right) Schrödinger suggest to call Eq.(1.4) as the ”amplitude equation.” He acknowledges, that such an equation is valid only for the conservative (that is time-independent) systems. Any process of measurement involves, however, time-dependent perturbation []. Accordingly, ”we must search for the real* wave equation* ” since ” the amplitude equation is no longer sufficient”. Thus, Eq.(1.4) is no longer real wave equation for Schrödinger. Furthermore, on page 102 he writes: ” Equation (1)111111That is our Eq.(1.4) contains the energy-or frequency-parameter , and is valid, as expressly emphasized in Part II121212That is, read pages 13-40 of [1]., with a definite -value.” Here the italics are Schrödinger’s. This statement by Schrödinger is very confusing, nevertheless, since in Eq.(1.4) the parameter is determined by solving the respective eigenvalue problems. And this is exactly what Schrodinger did repeatedly, starting with his 1st installment on quantum mechanics! Next, he suggests several ways of arriving at the ”real wave equation” which we do not want to reproduce here since these results are incorrect (as Schrödinger acknowledges himself). After several wrong attempts, he finally writes: ”The dependence of on the time which must exist if (131313That is our Eq.(1.4). to hold, can be expressed by
[TABLE]
as well as by
[TABLE]
Obviously, ”as well as” only means that: a) either his statement is erroneous (since equations (3.1a) and (3.1b) are mutually exclusive because Eq.(3.1a) is leading to the parabolic-type partial differential equation (PDE) while Eq.(3.1b) is leading to the hyperbolic PDE, b) he had in mind to use two times for the same equation141414Incidentally, two times interpretation of the time-dependent Schrödinger equation was proposed quite recently [], in 2011., or c) he suggested to replace (in his opinion, definite value) the parameter ( present in the combination in his ”amplitude” equation (1.4) by in order to arrive at his ”true wave equation.” Finally, the option c) was selected by Schrödinger so that the obtained equation had become a ”true wave equation”, known now as Schrödinger’s equation151515Incidentally, selection of the Eq.(3.1a)) leads to the correspondence between the Heisenberg and Schrodinger’s version of quantum mechanics.. Selection of option c) disconnects this equation from the wave equation (2.8) and makes use of the De Broglie particle-wave arguments much more difficult to implement. On purely mathematical level, selection of the hyperbolic type equation (based on selection of Eq.(3.1b)) requires knowledge of the two initial conditions for the Cauchy problem to be set up correctly, while selection of Eq.(3.1a) requires to use just one initial condition ] for setting up the Cauchy problem as explained in Appendix A. This difference was noticed by Feynman. In section 1 we brought the following Feynman’s quotation [4], page 377: ”Actually Huygens’ principle is not correct in optics. It is replaced by Kirchhoff’s modification which requires that both the amplitude and its derivative must be known on the adjacent surface. This is a consequence of the fact that the wave equation in optics is second order in the time. The wave equation of quantum mechanics is first order in time; therefore, Huygens’ principle is correct for matter waves…”161616But not corresct in optics according to statement by Feynman cited in section 1. . Unfortunately, this statement by Feynman is not correct as we shall explain below. To begin our explanation, we start with some quotations from the book by Feynman and Hibbs on path integrals [] as well as from Feynman’s lectures on physics []. In section 1-1 of [], on the 1st page, we find the following statement: ” There is one lucky break, however-electrons behave just like light. The quantum behavior of atomic objects (electrons, protons, neutrons, and so on) is the same for all, they are all ”particle waves”… So, what we learn about the properties of electrons…we will apply also to all ”particles” including photons of light171717Here italic empasis is ours. Historically, the events occured just in reverse order though, as explained in section 1.”. The interpretation of outcomes of the double slit interference experiment is central for understanding of quantum mechanics as explained in section 1 and in [14]. Because the light is being treated as ”particle waves” by Feynman,we shall begin with the light in the context of the double slit experiment in optics. This experiment was described in detail long before quantum mechanics was born [3,8, 33]. It involves use of Huygens’ principle. Mathematically accurate definition of this principle is given in the next section. For now it is sufficient to think about this principle as follows.
Definition 3.1. The Huygens principle. The new wave front at later time is the envelope of secondary waves emanating from each point of the original wavefront
.
Stated in such a form it leaves entirely open the problem: Why the envelope of secondary waves is moving only forward? Later on Fresnel added a superposition principle for the amplitudes of the secondary waves to explain the phenomenon of diffraction. Subsequently, the same problem was looked upon by Kirchhoff, Fraunhofer, Beltrami, Volterra, Hadamard, Arnol’d, and many others. In optics literature the description of diffraction begins with study of the scalar Helmholtz equation
[TABLE]
If we replace in (2.8) by and introduce then, Eq.(2.8) becomes the Helmholtz equation, provided that the potential . In such a form Eq.(3.2) is used in scattering theory of quantum mechanics [. It is appropriate to notice at this point that it is commonly believed that only Schrödinger’s ansatz, Eq.(3.1a), makes the wave function complex and, with this, it makes Born’s interpretation of quantum mechanics possible. This opinion is erroneous though since the simplest plane wave solution of Eq.(3.2) is known to be [34]
[TABLE]
leading to the probability current
[TABLE]
The question arises: Is it possible to find an analog to this current in optics? From the book by Born and Wolf [], chapter 8, it follows that the analog of the current in optics can be obtained as result of application of Green’s theorem, also used in quantum mechanics for derivation of the current, Eq.(3.4). In optics is an amplitude of the wave. Experimentally the intensity is measured. It is expressible through the combination In optics it is possible to replace the vector wave equations with the scalar ones when diffraction is discussed. Eq.(3.2) is the scalar wave equation for . In optics the intensity is measured experimentally while in quantum mechanics the very same quantity is the probability. Very detailed analysis of similarities and differences between optics and quantum mechanics was made in our work, ref.[14]. It was stimulated by the footnote on page 387 of [8]. In it, the reference to the 1959 work by Emil Wolf was made. In this work it is demonstrated ”that both the (time averaged) energy density and the energy flow in unpolarized quasimonochromatic (optical) wave field may always be derived from one component complex and time-harmonic scalar wave function.” Neither Born nor Wolf attempted to develop these observations subsequently. This was done in [14] where it is demonstrated that Born probabilistic interpretation of quantum mechanics and the optical interpretation (especially for the double-slit experiments) can be made mathematically completely indistinguishable. Initially, the discrepancy of interpretations of quantum mechanics was noticed by De Broglie. On Page 127 of his book [] we find the following statement: ”The probability that the presence of the photon will be made known by photographic action in the apparatus is everywhere proportional to the resultant intensity of the wavetrain…..It is almost certain that the same considerations are valid for the diffraction of material particles181818The italics are ours.”. Thus, results of [] provide rigorous mathematical support to De Broglie’s intuition.
From Feynman’s lectures [], it follows that in all double slit experiments the detectors were used (connected with loudspeakers) to produce ”clicks.” The number of clicks was counted per unit time as a function of the position () of the detector. Clearly, the detectors register particle hits. These are not exactly the probability amplitudes but, traditionally treated , they are represented through these amplitudes nevertheless! These amplitudes were experimentally calculated with the help of protocol described in section 1-1 of [ ] and, indeed, when compared against the interference patterns for light, the complete agreement was found, e.g. read page 8 of section 1-1[6] and page 5 of []. Because of this, further steps were made. On pages 5 and 6 of [], without saying it explicitly, Feynman formulates the Huygens-Fresnel principle in accord with the results known from optics [], page 371. Specifically, on pages 5-6 of [] we find the following description of the 2 hole experiment: ” Furthermore, is the sum of two contributions: the amplitude of arrival through hole 1, plus , the amplitude of arrival through hole 2. In other words, there are complex numbers and such that and Here we say only that , for example, may be evaluated as a** solution of a wave equation representing waves spreading from the source to 1, and from 1 to x.191919This is exactly Huygens’ principle in the form of the ”major premise” by Hadamard described in Appendix B. This reflects the wave properties of electrons (or in the case of light, **photons)202020Here the emphasis is ours.”. But the ”motion” of photons is described by Eq.(2.6 a), while electrons-by Schrödinger’s wave equation Eq.(1.3). Hence, the mathematical description is visibly different: Two different equations- the hyperbolic (wave) and the parabolic (Schrödinger) cannot describe the same reality in the same way! The confusion is caused by the fact that the Helmholtz Eq.(3.2) is used both in optics and in quantum mechanics, e.g. read survey on atomic optics [].When time-dependence effects can be ignored, then mathematically these equations are indistinguishable! With this observation in our hands we came all the way back to Schrödinger’s dilemma: which one of the equations: Eq.(3.1a) or (3.1b), should be used? Please, recall that the switch from (3.1b) to (3.1a) was caused by Schrödinger’s desire to develop correct description of the time-dependent quantum phenomena. For the time-independent situations, however, it is completely safe to use the wave Eq.(2.8) producing the stationary Schrödinger Eq.(1.4). In this and the following (time-dependent) paper we shall resolve Schrödinger’s dilemma by treating the time-dependent case in such a way that the time-dependence is rigorously eliminated. In such a case, even for the time-dependent situations to be treated in the next paper we can begin with the wave Eq.(2.6a) and, by doing so, we shall bring into correspondence treatments of photons and the rest of ”particle waves” (in Feynman’s terminology).
As an input and guidance for what follows, we need now to comment on several recently published papers-all reflecting the existing confusion associated with uses of various definitions of Huygens’ principle in quantum mechanics. The confusion is exacerbated by several additional sources not described above. To illustrate our point, we begin with the Helmholtz Eq.(3.2). The Green’s function for this equation is obtained from the equation
[TABLE]
From quantum mechanics [], we obtain (for the outgoing standing wave)
[TABLE]
(for the incoming wave should be replaced by - in the exponent). Now, is not at all the Green’s function of the wave equation. Accordingly, one cannot use the Fresnel-Huygens principle in the form hinted by Feynman and implemented in []. Specifically, what the authors of [37] are calling as the Huygens-Fresnel principle
[TABLE]
is not all mathematically valid expression! Accordingly, the results of [] are not valid. Since these results question the double slit experiment lying at the heart of quantum mechanics, we now provide the detailed explanation of what went wrong with the results of []. First of all, the propagator does not admit the path integral representation and, accordingly, it does not posses the Markovian property essential for all path integral treatments. This could be seen by direct calculation of Eq.(3.7) in which Eq.(3.6) is used. For the record we quote the result from the book by Ito and McKean [], paragraph 7.21, containing a remarkable formula connecting the Brownian motion propagator with the Newton (or Coulombic) potential
[TABLE]
Here is the dimensionality of space, On the left hand side of this equality we find well known Feynman’s propagator for the free particle of unit mass upon transition Accordingly, for this propagator one can safely use Eq.(3.7). Unfortunately, Eq.(3.8) cannot be replaced by that involving the Yukawa-type potential, Eq.(3.6). While the attempts to find an identity analogous to Eq.(3.8) for the Yukawa-type potentials, e.g. Eq.(3.6), are still ongoing, it is obvious that, based on current knowledge of this topic, Eq.(3.7) is invalid. Because of this, the main statement of [] questioning correctness of the superposition principle of quantum mechanics is also incorrect. This is so because of the following. The quantum mechanical superposition principle for two slits A and B in its orthodox form reads: However, the authors of [], following some earlier works, write instead . The third term, supposedly, originates from the entanglement of the Brownian paths originating at the source with two slits, e.g. A and B. That is to say, accounts for the possibility of the Brownian path going through the slits and coming back to the source/the origin. However, according to Huygens’ principle the secondary waves can move only forward! According to Feynman and the cited experiments, both the light and the particle waves produce the same interference patterns. In optics the description of such patterns does involve uses of Huygens’ principle. Therefore, the same should be true for the ”particle waves.” Accordingly, the equality is incorrect. In many papers attempting to use Feynman’s path integrals for description of the double and triple slit experiments, e.g. read [] and references therein, the authors used the Markovian analog of Eq.(3.7) (e.g. see Eq.(2.31) on page 37 of the book by Feynman and Hibbs [). Use of this analog is illegitimate though because the Brownian propagator in Eq.(3.8) when integrated over time time is not the propagator for the wave equation. The parabolic time-dependent Schrödinger equation is surely not the same thing as the hyperbolic wave Eq.(2.6a) since in the first case we need just one initial condition while in the second-two. For the wave equation Huygens’s principle rigorously holds (as explained in the next section and Appendix B). An attempt at rigorous implementation of Huygens’ principle in quantum mechanics formulated in terms of Feynman path integrals was made by Gutzwiller [5] whose sketchy and incomplete results were discussed in the 1st section.
Replacement of the hyperbolic (wave) equation by the parabolic (Schrödinger) could be a source of wrong conclusions about the outcome of three-slit experiments with photons. These were reported in ”Science” in 2010 [] with the purpose of rigorous testing of the superposition principle. These experiments found no deviations from the superposition principle of quantum mechanics. The quality of these experimental studies was further improved in [] with the same outcome thus confirming the validity of both, Huygens’ and superposition principles. Subsequent analytical and computer studies of the three slit configurations [], culminating in the latest PRL [] -all criticized the experimental results of Sinha et al [] as well as those by Söllner et al [41]. This critique is invalid however for reasons already explained. We noticed in section 2 that the De Broglie relation is valid also for photons. But the dispersion relation for photons is while that for the Schrödinger’s equation for the free particle of mass is (e.g. read Appendix A). This difference in dispersion relations is equivalent to the replacement of Eq.(3.1b) by (3.1a) as done by Schrödinger. It is fundamental since only under such conditions Heisenberg and Schrodinger pictures of quantum mechanics are in agreement with each other! The paradoxicality of the existing situation can be understood by reading some earlier descriptions of Huygens’ principle, the descriptions of diffraction, etc. []. They all involved the wave equation and time-dependence. However, this time-dependence happened to be not essential as could be seen, for example, by reading the authoritative book on optics by Arnold Sommerfeld [42. Subsequently, the time-dependence was wiped out from books on optics [8,43]. Once it was wiped out, use of the time-independent optics formalism, e.g. describing the double slit diffraction, makes it possible to replace the optical formulas by the quantum mechanical ones. More details are given in []. Inclusion of time dependence leads to dramatic effects. Details will be provided in future publications. This can be seeing already from the observation that switching from the hyperbolic (wave) Eq. (2.6a) requiring two functions for the initial conditions (the Cauchy data) to the parabolic (Schrödinger) equation requiring only one function is highly nontrivial , e.g. read subsection 6.2.1. below. Being parabolic, the Schrödinger equation possess the Markovian property while the wave equation is not. To by pass this difficulty, the semigroup analysis of operators was invented [] 212121Also, read the subsection 6.2.1 and ref.[14] where the Duffin-Kemmer formalism is briefly discussed leading to analogous results. In section 1 we quoted Gutzviller [5] who said that ”neither Hyugens’ principle nor the path -integral applies to the stationary waves which are the solutions of Eq.s (1.4) and (1.5)”. This statement of Gutzviller happens to be incorrect as we shall demonstrate below. Furthermore, in [5] Gutzwiller was using the semigroup analysis [44] for claiming that the time-dependent Schrödinger equation is obtainable from the relativistic (wave-like, a la Dirac, equation ) in the limit when First of all, this limit is unphysical (see, however, Appendix A) and, second of all, we shall demonstrate that Gutzwiller’s claims regarding applicability of Hyugens’ principle to Eq.s (1.4) and (1.5) is also incorrect. However, his claims about the Markovian property of path integrals are completely consistent with ours. Eq.(3.7) is surely non Markovian. There is no way of rewriting Eq.(3.7) in the path integral form. Thus, the Markovian property is not coinciding with the Huygens’ principle property! The Hugens principle (and property) in its refined form permits propagation of wave packets only forward (as explained in the following section) while Schrodinger’s equation allows spreading of the wave packets both ways [45]. Interestingly enough, the Huygens principle had been put by Schrödinger into center of his developments of wave mechanics. This can be seen from reading of Part II [pages 13-40, of his four foundational installments on quantum mechanics. In view of the variety of opinions in physics literature about the essence of Huygens’ principle (as compared to the variety of opinions about Huygens’ principle prevailing in mathematics literature), in the next section (and in Appendix B) we are presenting some basic facts about Huygens’ principle as it is developed by mathematicians.
4. Basic facts about Huygens’ principle (mathematicians perspective
)
4.1. Wavy (that is partial differential equations ) -type versus contact-geometric-type aspects of Huygens’ principle
Definition 3.1. provides the essence of Huygens’ principle. In plain words Huygens’ principle can be concisely formulated as follows. Consider a point at the wave front at the moment as the source of a new (secondary) wave emanating from this point. The Huygens principle (HP) as formulated by Huygens himself leaves entirely open the problem: Why the envelope of secondary waves is moving only forward? In the previous section we stated that to fix this problem subsequently Fresnel added a superposition principle for amplitudes of the secondary waves. This helped him to explain the phenomenon of diffraction. We had already mentioned that the same problem was looked upon later by Kirchhoff, Beltrami, Volterra, Hadamard, Arnol’d, and many others. Different authors used the HP with different purposes in mind. This caused differences in interpretation and confusion among users. Basically, the split in interpretations originates from the split of opinions about what is light. If we take the side of proponents of the wavy nature of light, then it is instructive to read carefully the fundamental work by Hadamard [ad also [29], chapter 6. If we take the side of proponents of corpuscular nature of light, then it is instructive to read woks by Arnol’d [47-49] in which methods of contact geometry and topology are used. To our knowledge, the detailed connections between these two directions of thought still do not exist. Works by Maslov [50,51] indicate that the connection problem is solved, in principle.
It is appropriate again to bring the quotation by Schrödinger on this topic [ On page 154 he writes:222222This is conclusion of his Nobel Address delivered on December 12th, 1933.” The light ray, or track of the particle, corresponds to the longitudinal continuity of the propagating process (that is to say, in the direction of spreading); the wave front, on the other hand, to the transversal one, that is to say, perpendicular to the direction of spreading. Both continuities are undoubtedly real. The one has been proved by photographing of particle tracks, and the other by interference experiments. As yet we have not been able to bring the two together into uniform scheme. It is only in extreme cases that the transversal -the spherical-continuity or the longitudinal -the ray -continuity shows itself so predominantly that we believe we can avail ourselves either of the wave scheme or the particle scheme.” Here the italics are Schrödinger’s. To complete this chain of quotations we cite Arnold’s ”Lectures on Partial Differential Equations” [ At the beginning of Lecture 1 he writes: ”In this lecture we shall consider a case in which there is a complete theory, namely the case of one first order equation. From the physical point of view this case (displays) the duality that occurs in describing a phenomenon using waves* or* particles. The field (of waves) satisfies a certain first-order partial differential equation, the evolution of the particles is described by ordinary differential equations, and there is a method of reducing the partial differential equation to a system of ordinary differential equations; in that way one can reduce the study of wave propagation to the study of the evolution of particles.232323Words put in curly brackets as well as italic emphasis are ours”
In view of particle-wave duality just described it is appropriate to define* the dual of Huygens’ principle. *In the theory of ordinary differential equations [52,53] (that is for ”rays” in optical terminology) there is so called rectification theorem which, when translated into language of Hamiltonian dynamics, is known as symplectic rectification theorem [54]. In plain words it is the statement which can formulated as follows.
Definition 4.1. The dual Huygens’ principle. Any point on the dynamical trajectory242424That is ray trajectory can serve as the origin of (new) motion.
That is to say, at any point along the trajectory it is possible to find canonical variables which are making the Hamiltonian of the dynamical system to vanish. At all such points along the trajectory the canonical Hamiltonian equations become trivial. When such dynamical system is quantized, the unitary operator describing system evolution becomes trivial too. Since the dual Huygens’ principle applies to any point along the trajectory in phase space, this means that Heisenberg’s uncertainty principle fails. The way out of this difficulty was suggested by Maslov [50] who suggested to use the so called Lagrangian manifolds . These manifolds naturally arise in contact geometry [47-49]. An introduction to this field can be found in [55]. In the context of Hamilton-Jacobi formalism the Lagrangian manifolds naturally originate based on familiar equations of the type Here , is the dimension of the configurational (tangent) space, is the action of the dynamical system. For the prescribed just stated set of equations defines the Lagrangian manifold . Given that the dimension of the symplectic manifold is , the dimension of the Lagrangian manifold is . Since the symplectic manifold is determined locally by the 1-form , the Lagrangian manifold is determined by the requirement: At the same time, the contact 1-form is determined by while the contact plane is determined by the condition The Bohr -Sommerfeld quantization conditions for trajectories on are trivial, that is since by design on where is the Poisson bracket. Existence of the Lagrangian manifolds does not preclude the existence of quantum mechanics, at least in the asymptotic sense [50], even though the Groenvold- van Hove theorem [56],[57] causes severe difficulties in systematic development of quantum mechanics.
Resolution, if any, of the above mentioned difficulties caused us to discuss only the wavy (that is, the PDE) side of Huygens’ principle []. The dual (particle) side requires uses of methods of contact geometry and topology as we just explained. Much more detailed account of this side of the Huygens principle is given in ], Chr.10. Reference [] contains some alternative treatment while very readable but rigorous ref.[60] claims along with Schrödinger [1], part II, that Huygens’ principle and the Hamilton-Jacobi equation are equivalent concepts. More on this contact/symplectic aspects of Huygens’ principle will be presented in our forthcoming publications.
We begin our discussion of the wavy aspects of the HP by summarizing results by Hadamard [ In ref.[ he formulated the following
Hadamard problem. How to describe/to classify all second order hyperbolic partial differential equations satisfying HP?
Hadamard discovered that the HP is working only in spaces in which the number of spatial dimensions is odd, e.g. in 3+1 spacetimes of Lorentzian signature. From here, it follows that the behavior of solutions of wave equations in 2+1 and 3+1 dimensional spaces is quite different. A short time point-like perturbation in R3 will produce sharply defined moving spherical wave so that when it reaches the observer he/she will see (or hear) just a short lasting flash/or splash. At the same time, for waves on the surface of water (this is an example used by Feynman, e.g. read the discussion pertinent to Fig.1-2 in [6] for his preliminary description of the two-slit experiment), that is in R the whole region inside the moving (spreading) circle will be disturbed. They say, that the 3 dimensional event is the result of the action of Huygens’ principle while the 2 dimensional event, say, on the surface of water, is caused by the wave diffusion. The following theorem proven independently by Mathisson[ and Asgeirsson[ is playing the central role in studying the HP -type of problems
Theorem 4.1. If the hyperbolic equation in space-time of Lorentzian signature is satisfying HP, then it is equivalent to Eq.(2.6a).
Thus, is the Huygens operator and a question emerges: How to define the equivalence? Following [ this equivalence can be defined as follows. Let L[ be a Huygens’ operator (that is the operator which is satisfying HP) and let L̃[ be another Huygens operator. Then they are equivalent if:
a) L̃[ can be obtained from L[ by non-singular transformations of the independent variables.
b) L̃[L[ for some positive, smooth function in the causal domain (to be defined momentarily).
c) L̃[ L[ for some positive smooth function in
To define the causal domain we need to define the distance function [ . It is given by
[TABLE]
Each pair is called an event. The Euclidean line segment between two events and can be used to formally define the velocity This definition then allows us to define an open (future) set D of those events that can be reached from with the velocity . Analogously, it is possible to define an open(past) set D of events from which can be reached with the velocity Thus, D if D The boundary of D is called forward (future) (respectively *backward (*past)) characteristic cone C (respectively, C The lightcone is defined by the requirement This is an equation of a 3d sphere for a fixed Based on this information, we are now in the position to define
**Definition 4.2. Hadamard criterion. The operator L[ is of Huygens-type if
[TABLE]
Accordingly, the causal domain is determined by the following requirement [ : event only if C
Definition 4.3. A Huygens operator L[ that arises from via operations a), b) and c) is called trivial Huygens operator.
Hadamard conjecture**. **Every Huygens operator is trivial.
Evidently, the Hadamard conjecture is fully compatible with Theorem 4.1. The question arises: If at any point of (pseudo)Riemannian 3+1 spacetime the metric can be brought into the diagonal form of Lorentzian signature, will the Hadamard conjecture be valid in some open domain of such space? Some studies of this problem were made by Friedlander [ McLenagan [ Goldoni [ Ibragimov [ , Wünsch et al []. To discuss these works further from physical standpoint, we need to relate them to Schrödinger’s work [1]. Incidentally, Schrödinger later on in his life studied the problem of wave propagation in curved spacetimes [24].\More details on this topic are given in the Appendix B.
4.2. Connection of mathematical works on Huygens principle with Schrödinger’s foundational papers on quantum mechanics
4.2.1. General background
Already in section 1 we noticed that for Schrödinger (2nd installment in [1]) HP is synonymous with the H-J equation. This point of view is being shared by such famous mathematicians as Gelfand and Fomin [60], pages 208-217. In our opinion, based on [29,33], Gelfand and Fomin results can be presented in a such a way that their consistency with results of Hadamard and other mathematicians who used the PDE methods for description of the HP will become obvious. Following [33], we shall assume that at time a light signal had been originated at the point (x0,y0,z0). At the later time this signal had penetrated into domain of space enclosed by a surface given analytically in the form
[TABLE]
The surface is called wavefront surface. The velocity of the wavefront can be measured with help of the velocity along the trajectories orthogonal to the wavefront (Appendix A). For points P1 and P2 on such trajectory we obtain:
[TABLE]
or, in the differential form,
[TABLE]
Mathematically, this result is identical to the condition for the contact plane mentioned in the previous subsection. From this observation it follows that it is perfectly reasonable to apply methods of contact geometry mentioned in the previous subsection for description of the HP. At the level of PDE this principle can be further elaborated as follows. Since Eq.(4.3a) is describing a wavefront, that is a two -dimensional surface in three-dimensional space, it makes sence to introduce coordinates on this surface so that
[TABLE]
Since the surface is moving, the above parametrization can be extended as follows
[TABLE]
Evidently,
[TABLE]
and
[TABLE]
To use these equations, we need to make some detour into theory of PDE following [29]. In particular, equations like Eq. (4.3a) are equations of characteristics. These are equations for some surfaces, e.g. These are encountered not only for the first order PDE’s known to standardly trained physics professionals but also to the higher order PDE’s too. Suppose, we are having, say, the 2nd order PDE, e.g like Eq.(2.6a), for which the Cauchy probem (Appendix A) can be set up and is well posed. Then, one can think about extending the Cauchy initial values, e.g. see Eq.s(4.5a),(4.5b), prescribed on to solutions of, say, Eq.(2.6a).
For this purpose it is useful to consider some auxiliary problems first. E.g. let us consider a cone, Eq.(4.1), defined on C Following [29], page 558, we introduce the function By direct calculation we obtain ():
[TABLE]
On the cone we have the condition leading to the 1st order PDE
[TABLE]
Just obtained resullt admits broad generalization. For instance, consider instead (c=1) Then, we again obtain
[TABLE]
Clearly, we can make these results to hold in any number of dimensions. Furthermore, Eq.s (4.6b) and (4.7) are looking the same. Now we demonstrate that: a) there are many other functions than satisfying the same 1st order PDE, Eq.(4.7), b) solutions etc. are also solutions of Eq.(2.6a). The demonstration can be achieved with help of the progressive wave solution method developed by Friedlander [] in 1946. Since his resuts, to our knowledge, are unknown in physics literature, we would like to make a detour and to discuss his results. This will enable us to develop a variety of physical applications.
4.2.2. Progressive wave solution by Friedlander. Introduction
We have to look for solutions of Eq.(2.6a) in the form
[TABLE]
where is an arbitrary well behaving function. The ansatz (4.8a) can be further complicated if, instead, we shall look for solutions in the form
[TABLE]
where
[TABLE]
In this work we shall use only the ansatz given by Eq.(4.8a). Substitution of this result into Eq.(2.6a) results in the following three coupled PDE’s :
[TABLE]
[TABLE]
[TABLE]
Clearly, the simplest case is obtained when . In this case we are left only with Eq.(4.9a). Eq.(4.7) will coincide with Eq.(4.9a) if we choose In view of Eq.(2.8) the obtained results can be extended further. For instance, Eq.(2.6a) can be replaced by
[TABLE]
where is some dimensionless function of spatial variables. In optics it is associated with the index of refraction []. With such a replacement, Eq.(4.9a) is changed into
[TABLE]
Next, we return back to Eq.(4.3c) which we rewrite as
[TABLE]
With a wavefront, Eq.(4.3a), a two parameter ( family of rays orthogonal to the wavefront is associated. A current point on a ray trajectory is described by {x(
Let (x( be some yet arbitrary function on this trajectory. Then,
[TABLE]
Without loss of generality, we can let (x(x( Furthermore, let
[TABLE]
where is the absolute value of the velocity. With such identifications we now replace in Eq.(4.12) by and use the r.h.s. of Eq.s(4.14) in Eq.(4.12) to arrive at
[TABLE]
But Therefore, by squaring we reobtain back Eq.(4.11). In addition, however, using Eq.s (4.14) we obtain
[TABLE]
Accordingly,
[TABLE]
By combining Eq.s (4.5b) and (4.16b) we obtain:
[TABLE]
This result is holding on along with Eq.s(4.5a). Therefore Let in Eq.(4.3a) and rewrite Eq.(4.3a) as
[TABLE]
so that at the characteristic surface is . During the time evolution each point ( of creates its own wavefront according to Eq.(4.3b). Thus, a two-parameter set of wavefronts is obtained. The essense of Huygens’ principle lies in demonstration that the envelope of all these wave fronts is again a wavefront such that for we reobtain The procedure for finding an envelope is reduced to eliminating the parameters ( from three equations
[TABLE]
[TABLE]
[TABLE]
Suppose that
[TABLE]
are calculated from Eq.s(4.18a) and (4.18b), respectively, then is a solution of Eq.(4.11). Indeed,
[TABLE]
in view of Eq.s(4.18a),(4.18b). Analogously, and implying that Eq.(4.11) holds. Thus, the surface (the characteristic) at time is converted into surface described by Eq.(4.3d). Evidently, since it is obeying the same H-J equation (4.11) in view of Eq.(4.20) (and those for y and z components). As we know already (e.g. read section 1), the essence of Huygens’ principle for Schrödinger is the H-J Eq.(4.11). With explanations just made superimposed with those in [60], it should be clear to our readers that this is indeed the case.
4.2.3.
Progressive wave solution method by Friedlander. From mechanics of Bohm
to gauge-theoretic mechanics of Floer
Results of previous subsection were obtained under the assumption that in Eq.(4.8a) is a constant. If it is not a constant, situation becomes much more complicated mathematically []. To decide what to do with these complications physically requires some work. First, we notice that the H-J equation was used in the 1st paper by Schrödinger on quantum mechanics [1], pages 1-12, as an input for obtaining the stationary Schrödinger equation. Although his method of deriving this equation was endorsed by Courant and Hilbert [2], pages 445-450, to our knowledge, it was left without attention in physics literature. For the sake of results we shall develop momentarily, we are going to reproduce the 1st Shrödinger method now. For this purpose, we should notice that: a) the ansatz, Eq.(4.8a), was made without account of dimensionality arguments; b) when the dimensionality arguments are taken into account, the H-J equation, e.g. Eq.(4.11), becomes a simple statement about the classical momentum of the particle252525Following Schrödinger, and for the sake of argument, we are discussing only the one and two-body problems.
[TABLE]
Schrödinger makes the following reversible substitution ln into Eq.(4.21) resulting in his Eq.() (ours Eq.(1.4)). For our readers convenience we rewrite it here again:
[TABLE]
Because of noticed reversibility, this is still the classical H-J equation, even though it has in it. It is exactly the same as Eq.(4.21). Next, using Eq.(4.22), instead of Eq.(4.21), Schrödinger considers the following optimization problem: Find the minimum of the functional
[TABLE]
under the subsidiary condition
[TABLE]
The result of such a minimization is the stationary Schrödinger Eq.(1.4). By design, any solution of Eq.(1.4) should be a minimum of If coming as solution of Eq.(1.4) is such that then it is describing the classical trajectory according to Eq.(4.22) since the substitution ln\psi\is reversible. But, in view of results of previous subsection, this is indeed the case! E.g. solutions of Eq.(4.11) are solutions of Eq.(4.10) if in the ansatz Eq.(4.8a) we replace by and take into account Eq.s (2.8), (2.9). This observation allows us to reinterpret variational results of the 1st Schrödinger paper in terms of the gauge-theoretic formalism developed by Floer. A quick introduction to Floer’s theory is provided in [71] while the detailed account can be found, for instance, in [72]. Even though Floer theory studies a multitude of closed orbits on symplectic manifolds classically, the extremely sophisticated computational methods developed by Floer exactly parallel those used in nonperturbative (instanton) treatments of the Yang-Mills theory. These results establish a connection between the non relativistic Schrödinger equation (with time-independent Hamiltonian) and the gauge-theoretic Yang-Mills-type theory. We shall say more on this subject in subsection 7.4.
If in Eq.(4.8a) is not a constant then, at the very least, we end up with the formalism of quantum mechanics developed by David Bohm [73]. See also Appendix C. It also can be considered as classical because quantum mechanical corrections in Bohmian mechanics are treated by methods of classical mechanics (in fact, of hydrodynamics). Bohm formalism is based on representation given by Eq.(4.8a) with an extra restriction It remains to be investigated how results of Bohm formalism might change if Eq.(4.8b) is used instead. This task is left for further study.
4.2.4. Progressive wave solution method by Friedlander. Dupin cyclides
, the Lie
sphere geometry
and the conformal group SO(4,2)
From what was discussed thus far it follows that both the hyperbolic (wave-like) and the parabolic (Schrödinger-like) equations can be treated with help of the progressive wave solution method. If this is so then, according to Friedlander [70], the most general solution in both cases should be expressible in terms of cyclides of Dupin. For the hyperbolic PDE’s this was demonstrated by Friedlander [65,70]. In this subsection, and also in subsection 7.5.3. and Appendix C, we shall demonstrate that this is also true for the stationary Schrödinger equation. By doing so, we shall put the obtained results in a much broader context. This context will allow us to find a place for (thus far) very exotic Dupin cyclides in atomic, high energy physics and cosmology.
To begin, we notice that 1946 result by Friedlander [70] was reconsidered in 2005 by Sym [74]. Not only he reobtained it in a much shorter and simpler form in ref.s[74,75] but, in addition, he was able to find some (not serious) mistakes in the original work. Both Friedlander and Sym were solving the system of Eq.s(4.9a-c). Sym very cleverly used symmetry to solve these equations. His solution strategy can be summarized as follows. First, solve Eq.(4.9a) by cleverly using the symmetry built into the Dupin cyclides. Second, to take the full advantage of this symmetry, rewrite Eq.s(4.9a-c) in curvilinear (actually geodesic) coordinates reflecting this symmetry. The choice of coordinates in Eq.s(4.9b,c) is determined by the fact that Eq.(4.9a) should admit a natural solution respecting the symmetry of Dupin cyclides. Third, use this solution in Eq.(4.9b) and insure that Eq.(4.9c) is solved in such coordinates as well.
Being armed with these results, we apply them to the Helmholtz Eq.(3.2) which is obtainable from both the D’Alembert and the free particle Schrödinger equations. Therefore, the results of separation of variables in Eq.(4.9c) can be used for solving the Helmholtz equation as well. The next step is made by checking what kinds of potentials in the stationary Schrodinger’s equation allow the full separation of variables with help of these geodesic coordinates. This problem is non trivial. It was studied, for example, in [76,77]. We shall not go into full details, however, since we shall develop much broader vision of the role of Dupin cyclides in such kind of problems which, in addition, will make such calculation much simpler and much more physically appealing.
For this purpose we need to explain: a) What are the Dupin cyclides? b) How their presence/absence affects the results of Schrödinger and Bohm? c) What happens to the Huygens principle if Dupin cyclides are taken into account? The answer to b) is given in part in the already cited [76,77]. But it is also given in Apendices B and C from the entirely different standpoint. The positive answer to c) follows in part from works by Sym[74,75] just mentioned and also follows from results of the Appendix B. Therefore, we need only to provide the answer to a). The answers to b) and c) are provided in just mentioned appendices and in section 7. The answer to a) is provided immediately below.
We begin with ref.[78]. It is associated with the notion of canal surface. Such a surface can be designed as follows. Choose some sphere whose radius is changing in time. The the center of is moving along some curve (the trajectory) parametrized by time so that such a motion is described analytically as
[TABLE]
Definition 4.4. Canal surface. The canal surface is an envelope of a 1-parameter family of spheres centered at the spine curve .
The envelope is being defined as joint solution of two equations
[TABLE]
Here is the scalar product in The canal surface is fully determined by a curve in 4 dimensional Minkowski space R Clearly, every point of R3,1 corresponds to a sphere
Dupin cyclides are made of canal surfaces but not another way around. This follows from the ingenious observation by J.C. Maxwell who constructed some of the first Dupin cyclides following the orininal work by Chales Dupin done in 1803. In 1868 Maxwell noticed that, even though every point belongs to a sphere , there could be two different sets: a) c1 and r b) and c2 and r2 producing the same canal surface
.
**Dedinition 4.5 a). Dupin cyclides are the only canal surfaces **which can be designed in two different ways.
Consider the simplest example-the torus
a) the spine curve the radius function
b) the spine curve the radius function
Using of a) and b) in Eq.s(4.26a,b) and eliminating leads to the following quadric (describing the torus-the canal surface)
[TABLE]
From this elementary example the following general definition follows
**Definition 4.5 b). Dupin cyclides are the surfaces whose lines of curvatures are circles. **Elementary examples of (degenerate) Dupin cyclides include spheres, planes (that is spheres of infinite radius), toruses and cones of revolution.
The above results admit generalization to higher dimensions [79]. More on this will be said in section 7. In the meantime, we would like to notice the following. An affine isometry **f **of Rn is described as follows. Let xRn, then f(x)=Qx+b where Q is orthogonal matrix (that is QQT=1) whose entries belong to R, b R
**Definition 4.6. **It f(c)=**c ** (respectively f( ), then **f **is called symmetry of c (respectively of
To get a feeling of this symmetry it is sufficient to consider the space R Since the canal surfaces are made of circles moving along the spine curve, the symmetry transformations should convert circles into circles. This is possible with help of the Möbius (or conformal) transformations C=R It is clear then, that the transformation f( is a conformal transformation. It is possible to extend these ideas to higher dimensions naively. This will yield the multidimensional Möbius transformations. These are not quite yet the transformations having the Dupin cyclides unchanged. They form a subgroup of the Lie sphere geometry group [80] which by design leaves the Dupin cyclides unchanged. In two dimensions the Möbius group is the isometry group of the Poincare*′* hyperbolic upper half plane. The interrelationship between the hyperbolic and the Minkowski spaces is nicely described in the book by Ratcliffe [81]. From it we find that the Möbius group is also group of isometries of the Minkowski spacetime. In sections 6,7 we shall demonstrate that the group of isometries of the compactified Minkowski spacetime is the conformal group SO(4,2). It is described, for example, in [19], pages 345-348. The significance of this group was noticed already in section 1. In section 7 we shall demonstrate that this group coincides with the Lie sphere group. This observation is of fundamental physical significance to be explained in section 7 and the Appendix B. In the meantime we would like to explain the differences between the Möbius and the Lie sphere groups. Since the Möbius group is the subgroup of the Lie sphere group, the additional symmetry elements originate from the fact that a) spheres could have an extra label- orientation. This can be seen already in dimension two. A circle can have the radius -vector with the base sitting on the circle and directed toward the center of circle. Another option is for the same vector to play the role of the outward normal to the circle. Besides, one should pay attention to the fact that the circles can touch each other and each circle could have an infinite family of circles sitting inside or outside of the given one and touching it at the same point. To these objects one should add points-circles of zero radius and, also -lines, the circles of infinite radius (in two dimensions) or the planes in dimensions higher than two. According to Klein’s Erlangen program, the following definition is the most appropriate.
Definition 4.7. The essence of Lie sphere geometry lies in the study of properties of transformations mapping oriented spheres (including points and planes) to oriented spheres while preserving the oriented contact of sphere pairs.
Corollary 4.8. Dupin cyclides are invariants of the Lie sphere geometry.
Corollary 4.9. The conformal group SO(4,2) was used in atomic [] and in high energy physics [] without any reference to the Lie sphere geometry symmetry group. It was discussed already in section 1. Use of conformal group caused introduction of the unphysical two-times formalism in []. Use of Lie sphere geometry removes this deficiency. In the rest of this paper we demonstrate how the conformal group SO(4,2) is identified with the Lie sphere geometry group. This identification opens the door for new results in conformal dynamics, conformal quantum mechanics, conformal quantum field theory, conformal gravity, etc. []. The AdS-CFT correspondence to be discussed in section 7 and Appendix D is part of this ”conformal program”. It is important to remember that the source-free Maxwellian electrodynamics is invariant with respect to SO(4,2) group as well [] while the results of conformal quantum mechanics recently were utilized with great success for recovery of the Regge-like spectrum of hadrons [].
5. Closer look at the Hadamard conjecture. Recovery of some known
equations for massive/massless particles with or without spins
In this paper we are not going to discuss the developments associated with counterexamples to the Hadamard conjecture. All these counterexamples are discussing the validity of the Huygens’ principle in spacetimes of dimensionality higher than four [. In the light of results of section 1 use of spacetimes of dimensionality higher than four for multielectron atoms and for molecules is also not necessary, apparently. Fortunately, in four dimensional spacetimes no counterexample to the Hadamard conjecture was found. Results of Appendix B should be considered as a cosmologically inspired ramification262626Associated with accounting for Penrose limits of physical spacetimes. of the already known results. Because of this, it is convenient to restate the Theorem 4.1. as
follows
Fundamental Principle 5.1. Quantum mechanical behavior of all elementary particles (massive and massless, of integer and half integer spin) is inseparably linked with the Lorentzian signature of ambient spacetimes.
Conjecture 5.2. Although not immediately obvious, the Fundamental Principle is synonymous with the central role of the Lie sphere geometry acting in conformally flat spacetimes of Lorentzian signature. Its influence on physics is ranging from conformal mechanics to conformal wave mechanics and conformal gravity. Group-theoretical classification of physically sensible spacetimes is contingent upon their ability to sustain quantum mechanics272727Further details are presented in section 7 and Appendix B..
5.1. Hadamard triviality and mass generation. Panoramic view
We continue this section with demonstration of Huygens’ triviality of telegrapher’s and the Klein-Gordon equations. Using Eq.(2.6) and the equivalence condition L̃[L[ with where is some constant, we obtain ():
[TABLE]
Next, let eibx and e Substitute these factors into Eq.(5.1) and apply again Hadamard’s equivalence rules. After a short calculation we arrive at
[TABLE]
If now and , we obtain telegrapher’s equation
[TABLE]
One dimensional version of this equation is discussed at length in the book [ In one and two dimensions this equation admits the path integral treatment. In the meantime, the Klein-Gordon (K-G) equation is obtained now if we make a replacement in Eq.(5.2b) with subsequent replacement of by . After this, we obtain:
[TABLE]
Thus, we just demonstrated that the K-G equation is Huygens-equivalent to the D’Alembert Eq.(2.6a). According to ], page 99, every spinor component of the Dirac equation with nonzero mass is satisfying the K-G equation. This fact establishes the Huygens equivalence between the Dirac and D’Alembert equations. Apparently, the particles with higher spin, e.g. spin-2 gravitons, etc. also belong to the same Huygens equivalence class [].
**Corollary 5.3. Huygens **equivalence between the D’Alembert and all relativistic equations of integer and half integer spin explains why the double slit experiments made with photons and massive particles produce the same fringe patterns.
More details on this is given in ref.[14].
Thanks to the seminal work by Mark Kac [ the propagator for telegrapher’s and Dirac’s equations can be presented in the path integral form at least in 1+1 dimensions []. In view of Eq.s (5.2b) and (5.3) the same is true for the K-G equation as explained in detail in [].
It should be noted that in all these cases the associated path integrals do not involve the Gaussian-type random processes. They are designed with help of the Poissonian- type random processes. Excellent description of these types of processes in conjunction with the telegrapher’s equation is given in []. The noticed connection with random walks is helpful but not crucial for the tasks to be completed in the rest of this paper. Furthermore, the above path integrals can be designed rigorously only in 1+1 and 2+1 dimensions. In [] it was rigorously demonstrated that these results cannot be extended to higher dimensions.
To stay focused, we are not going to discuss any further the connection between the random walks of various kinds and PDE’s. Instead, we notice the following
First. The diffusion/Schrödinger is not the wave-type equation studied by Hadamard. The diffusion/Schrödinger equation describes the dispersive waves discussed in the Appendix A. Therefore, contrary to Feynman’s claims made in [], it apparently does not obey the Huygens principle. However, the results of sections 4 and 6 indicate that this apparent deficiency of the Schrödinger equation can be repaired, and quite rigorously, for as long as the quantum Hamiltonian of this equation is manifestly time-independent. What remains to be proven is hyugens triviality of Schrödinger’s equation with time-independent Hamiltonian.
Second. Initially, Hadamard obtained his results in 3+1 dimensions as explained in previous section. Subsequently, he developed the method of descent allowing use of 3+1 dimensional results as an input for obtaining 2+1 dimensional solutions. Still later, these results were extended by others to 1+1 dimensions, again with use of the method of descent,[], pages 315, 316.
Third. Clearly, by going down from 3+1 to 1+1 dimensions some information is lost. Otherwise the distinction between, say, 2 and 3 dimensional wave propagation is going to disappear. But it is not! Therefore, the attempt to use the method of descent in reverse cannot help us in extending 2+1 dimensional results obtainable with help of, say, the Poissonian statistics to 3+1 dimensions. Other methods, e.g. those using Grassmann variables should be used instead []. Different method, also using Grassmann variables but employing differential-geometric considerations for 3 dimensional paths evolving in 3+1 dimensional spacetime (leading to 3+1 Dirac propagator) was developed in []. In connection with [], the following observation is appropriate.
Forth. The Standard Model of particle physics uses the widely accepted Higgs mechanism responsible for the mass generation. At the same time, the twistor formalism describing all massless particles had been extended recently by accounting for the rigidity of the world-lines of the massless particles []. Such differential-geometric mass generation method is analogous to that proposed in []. The latest results in this direction can be found in [].
**Fifth. **After Eq.(5.3) we stated that the K-G, the Dirac and other basic equations for massive particles can be made Huygens-trivial using Hadamard’s transformation rules. Since establishing of this equivalence is reversible process, this means that it is possible to avoid use of the Higgs mechanism. The two-time formalism developed by Itzhak Bars (summarized in []) leads to the same conclusions. Because the De Broglie relation is valid for massless photons as well as for the massive particles as demonstrated in the footnote 6 and because the stationary Schrödinger equation can be restored from the relativistically obtained uncertainty relations (as demonstrated by Schrodinger (e.g. read section 2)) the notion of mass should be linked with the De Broglie relation. Its actual value is controlled by the rigidity of world lines [97],[98]282828The latest results on this topic are mentioned in the ”Note added in proof”.
In view of these remarks, our first task is to address and to solve the problem formulated in the first remark. In section 3 we provided the following Schrödinger’s comments : ” Thus, when we designated equation (1) or ()292929Our Eq.(4.22). on various occasions as ”the wave equation”, we were really wrong…” Well, actually, Schrödinger was not at all wrong! Eq.(1.4) was used systematically by Schrödinger himself in Parts I-III of [ to solve the Hydrogen atom, the harmonic oscillator, the rigid and non rigid rotators and the Stark effect quantum mechanically. He also developed the perturbation theory based on the exact results which he obtained. All these results are correct! In section 3 we explained what made Schrödinger unhappy with his ”amplitude equation”, that is with our Eq.(1.4). In section 4 using methods of characteristics and progressive waves we demonstrated that Eq.(1.3) is indeed the correct wave equation in the sense of De Broglie [ obeying the Huygens’ principle in the sense of Hadamard. Furthermore, in section 6 we demonstrate that for the time-independen Hamiltonians Schrodinger’s ”real wave equation” can be embedded into the Hadamard scheme of calculations relating a given 2nd order PDE with the wave Eq.(2.6a).
In conclusion, we would like to mention that by using the two-time formalism Itzhak Bars established the equivalence (in the sense of gauge equivalence(duality) defined in his paper) between the D’Alembert Eq.(2.6a), on one hand, and the Schrödinger equations for the hydrogen atom and the harmonic oscillator, on another []. In [99] Bars acknowledged that he was not able yet to provide a complete classification of all gauge-equivalent(dual) quantum mechanical systems originating from the same (gauge-invariant) model. The attempt to do so was made later, in []. However, making choices between the two times in such a formalism still remained mysterious. Below, we shall obtain the same results using entirely different methods enabling us to avoid the two-times formalism altogether.
6. Huygens triviality of the stationary Schrödinger equation
6.1. Some comments about 1935 work by V. Fock on hydrogen atom
On February 8, 1935, Vladimir Fock presented his seminal lecture entitled ”On theory of the hydrogen atom” at the theory seminar of the Leningrad State University. English translation of his talk can be found in [101]. It is based on Fock’s article [102] published in German.
To begin our comments on his paper we rewrite Eq.(1.4) in the following standard form
[TABLE]
Without loss of generality we shall consider only potentials for which the above equation is exactly solvable. At the classical level use of canonical transformations makes all such exactly solvable problems equivalent since they all can be brought into standard action -angle form known for the harmonic oscillator Accordingly, it is sufficient to consider only the hydrogen atom problem for which in known system of units the potential is defined as For this (Coulombic) potential Fock replaces the partial differential Eq.(6.1a) by the equivalent integral equation
[TABLE]
Citing Fock, the rationale for using the integral equation method instead of solving the partial differential equation is caused by the following observations.
” It has long been known that the energy levels of the hydrogen atom are degenerate with respect to the azimutal quantum number l… But any degeneracy of eigenvalues is linked to the transformation group of the relevant equation: e.g. the degeneracy with respect to the magnetic quantum number is attributed to the usual rotational group.” Here Fock refers to the degeneracy of the following type. For a given* l* (associated with the rigid rotator energy l(+1)) one has 2l +1 wavefunctions labeled by the magnetic number He calls such a degeneracy ”accidental” and in his paper he finds the symmetry group causing this accidental degeneracy. He demonstrated that the group causing accidental degeneracy is four-dimensional rotational group .
Instead of copying Fock’s arguments, we shall arrive at the same results much more economically. For this purpose we introduce the notations:
[TABLE]
It can be easily demonstrated that in addition to the Hamiltonian H which is constant of motion, both the angular momentum L and the Laplace-Runge-Lentz vector A are also constants of motion [101]. These objects can be looked upon either at the classical level with subsequent quantization or at the quantum level. The last route was chosen initially by Pauli []. The results in both classical and quantum cases depend upon the value of the energy constant H and . In this work the case will not be considered since it is not related to the tasks we would like to accomplish. Traditional analysis involving bound orbits/states typically begins with the case . In this case it is more convenient to replace the vector A with the rescaled vector D defined as ([104], page 421):
[TABLE]
In the case of , the vector D is defined accordingly as
[TABLE]
In terms of such notations the Poisson brackets commutation relations are readily obtained with the result []
[TABLE]
Here Using this Poisson algebra, the respective Lie algebra is obtained via standard Dirac quantization prescription. As result, the so(4) Lie algebra (for ) and so(3,1) Lie algebra (for ) is obtained. The first one is the Lie algebra of the SO(4) rotation Lie group while the second is the Lie algebra of the SO(3,1) Lorentz group. In his talk Fock did mention the Lorentz group but provided no computational details.
6.2. Huygens’ triviality of the stationary Schrödinger equation in the light
of Fock’s work
6.2.1. General consideration
Just obtained results, when superimposed with results by Fock, are sufficient for group-theoretical proof of Huygens’ triviality of the stationary Schrödinger equation. Nevertheless, it is instructive to arrive at the final destination via extremely informative detour. This detour is possible to perform by using, for example, fundamental work by Gelfand, Milnos and Shapiro []. Alternatively, the same results can be obtained using the Duffin-Kemmer formalism []. The authors of [] discussed carefully the problem of classification of all relativistically invariant equations. That is of all Lorentz-invariant equations. Clearly, the protocol of study of these equations is the same as that for study of rotationally invariant equations. Therefore, we shall study both problems simultaneously. The study begins with the equation of the type
[TABLE]
where is the dimensionality of spacetime, is the mass parameter, and the metric tensor is given by
[TABLE]
Typically, but other dimensionalities could be considered as well. With such defined metric tensor, the 2nd order PDE, Eq.(6.5a), is convenient to rewrite in the form of the system of the 1st order PDE’s. Incidentally, in such a case, one can think about the method of characteristics for such PDE’s and the H-J equations, etc.[]. We shall not touch this topic in this section though since the relevant information was already provided in section 4. Instead, following [], we introduce an auxiliary functions via equation
[TABLE]
A simple calculation produces (for the Lorentzian case, )
[TABLE]
and (for the Euclidean case, )
[TABLE]
At this point it is convenient to introduce the extended wave vector (where means ”transpose”) and to rewrite the system of the above equations into the form
[TABLE]
Here are the matrices. E.g. in the Lorentzian case
[TABLE]
etc. Let be the matrix of either Euclidean rotations or of Lorentz transformations acting on coordinates as Let be some unitary operator such that
[TABLE]
Accordingly, Eq.(6.8) can be rewritten now as
[TABLE]
or, equivalently, as
[TABLE]
The condition of invariance follows immediately from Eq.(6.11b)
[TABLE]
Since this condition is unchanged in the massless case, it is sufficient to consider the massless case only. Clearly, the simplest equation invariant with respect to either the Euclidean rotation or Lorentz group transformations is the four-dimensional Laplacian in Euclidean case or D’Alembertian in Lorentzian case respectively. In his paper [] Fock obtained indeed the four -dimensional Lapacian equation as an equation whose solutions are those of the integral Eq.(6.1b). To obtain solutions in the Lorentzian case requires us only to perform the Wick rotation, that is to make a replacement: This innocently looking operation routinely used in physics literature required 126 pages of proof in mathematics literature []. Thanks to its existence, we are spared from the necessity to repeat needed proofs. The same results could be obtained group-theoretically via development of the unified treatment of SO(4) and S(3,1) Lie groups []. A transparent and motivating example of such an interrelation is given in the pedagogically written article by John Milnor [110]. Nevertheless, below we shall provide yet another derivation. It involves the twistor formalism. Use of twistor formalism enables us to present Fock’s results in a different light. To accomplish this task requires several steps which we would like to describe now. We begin with the following.
6.2.2. Grassmannians and Plücker embedding
First of all, we would like to take a careful look at the Hamiltonian H for the hydrogen atom at the classical level
[TABLE]
In particular, let initially . This condition connects the momenta and coordinates. In particular, it should also hold for In such a case, to maintain the equality, the momenta should become infinite. If initially we had , now we must add a point at infinity to keep the relation, Eq(6.13), unchanged. The addition of leads to the compactification of thus converting it to The group is the group of isometries of Next, instead of rescaling the Runge-Lentz vector A in Eq.(6.3) we can rescale the angular momentum . This is effectively done in Fock’s paper as we would like to explain now. For this purpose we follow [55], pages 361-364.
Consider 2 points x and y on the line in 3d space. Since the coordinates and ** are taken with respect to some fixed origin, **these are vectors. We can construct from these vectors two other vectors: and The new element is coming from the following step. We enlarge the embedding into Euclidean space by increasing its dimensionality, that is by replacing vectors x and y by xp={x0,x} and yp={y0,y}. With such an enlargement the vectors xyp can be looked upon either as vectors in R4 or as points in the projective space P Now we take into account that in 3 dimensions [55]. We would like to apply this correspondence to the vectors xp and yp. Specifically, we consider the exterior product (x0, x1, x2, x3) ,. For the sake of illustration, we can choose x0=y0=1, then Clearly, can be looked upon as components of the angular momenta L. In such a case, following Fock [101,102] we can make an identification yp={y0,y} with Evidently, in such a case we no longer can choose . This is not essential, however. In particular, this could be seen if we notice that are proportional to the components of the Laplace-Runge-Lentz vector. If this is so, then we require that
[TABLE]
However, for the Kepler problem we have:
[TABLE]
At the same time, Therefore, where and are known constants. By the appropriate rescaling this linear combination can be made the same as There is another meaning of Eq.(6.14a) however. It is purely mathematical. And, as such, it is totally independent of its relevance to the Kepler problem. It is associated with the concept of Grassmannian manifold and of Plücker’s embedding of this manifold into the projective space.
To shorten our discussion, without loss of generality we can always replace the real space by the complex Let {e be some vector basis in C Introduction of such a basis requires introduction of the scalar product. This can be done by analogy with the scalar products of quantum mechanics. With help of this basis we introduce the exterior products, e.g. ee . It is convenient to keep track of these products by defining the subsets Let the set of basis vectors {ee be associated with such a subset. There are exactly ways to select such a set from the set {ee and to make the exterior products e out of selected vectors. The Grassmaniann Grp(Cn) is the manifold made out of all -dimensional subspaces of Cn. Suppose that in we changed the basis from {e to { Changes of the basis in Cn leads to changes in the basis for subsets, e.g.
[TABLE]
This relation can be looked upon as the equivalence relation defining a point in the complex projective space CPN. Thus, we just (Plücker) embedded the Grassmannian Grp(Cn) into the complex projective space CPN of dimension , where Since we also have Grp(Cn)=Grn-p(Cn) causing the exterior products e and e to be related to each other in the way known from the Hodge theory of differential forms. Specifically, if we introduce the notations: e and e , this definition is implying that in both cases we use the basis {e so that e . If this is so, we can then construct the products of the type ( ) Construction of such products is subject to the Hodge-type constraints ( ) The sign means (Hodge-type) ”equivalence”. Evidently, it is permissible to have as well the following equivalences ( ), and so on. Just described results lead us to the following
Definition 6.1. Let {ee be the basis for C then we define the set ee with If , then is totally decomposable if
[TABLE]
The homogenous coordinates are called Plücker coordinates on CP
From this definition, it follows at once that
[TABLE]
This is the condition for the Plücker embedding. We want to demonstrate now that Eq.(6.14a) is the condition for the Plücker embedding. For this purpose, consider now the Grassmannian Gr2(C4). For it, we have {e0,…,e3} as the basis for Next, we construct built as follows:
[TABLE]
Since and the condition for the Plücker embedding acquires the following form:
[TABLE]
From here we obtain:
[TABLE]
where the determinants are defined by
[TABLE]
Since, , while and so on, the Plücker embedding condition, Eq.(6.18), is exactly the same as the mechanical condition, Eq.(6.14a).
6.2.3. From Plücker embedding to conformal compactification of Minkowski
space via twistor formalism
The results obtained in previous subsection enable us to complete our proof of Huygens triviality. In the previous subsection we discussed the rationale for compactification of the momentum space . Naturally, we ended up with The 3-sphere is living in . There is some advantage in compactification of as well leading to The argument goes as follows. In the previous subsection we introduced vectors x yp living in R We noticed that they can be looked upon either as vectors (in ) or as points (in ). By complexification we end up with these vectors living either in C4 or CP This allows us to develop the compactification of the Minkowski space (whose isometry group is SO(3,1)) and the Euclidean space (whose isometry group SO(4) ) using the same formalism. Furthermore, using the compactification procedure leads us directly to the formalism of twistors and twistor spaces.
Definition 6.2. The twistor space T is C4 with coordinates Z={Z Z Z ZThe projective twistor space is PT=CP3 with homogenous coordinates {Z Z Z Z
Next, we notice that if for the description of the complex plane we have to use the complex numbers, e.g. then for description of C2 we have to use the quaternions, e.g. . Both complex numbers and quaternions admit matrix presentation []. For instance,
[TABLE]
The scalar product in can be presented either as or as The last presentation can be equivalently rewritten as
[TABLE]
The squared norm = of the complex number can be alternatively represented as
[TABLE]
To extend these results to we notice that The quaternion is encoded by the . In view of the correspondence we introduce 2 complex numbers and and then, by analogy with Eq.(6.19a), we can write
[TABLE]
Accordingly, the squared norm = is given by
[TABLE]
Consider now a spacetime of Euclidean signature in which the vector is given componentwise as
[TABLE]
The square of the Euclidean norm for is given by Here Because of this, we also have Furthermore, we also obtain: Accordingly, Next, we recall that the 3-sphere can be analytically represented as The Hopf map can be constructed now as follows. Begin with and consider the ratio (or Since , it is clear that Alternatively, Consider now two quaternions ZZ and ZZ and consider the quaternionic projective space with homogenous coordinates { If in the case of complex numbers we had the correspondence equivalent to the Hopf map in the quaternionic case we analogously have another Hopf map: To demonstrate that this is indeed the case, we recall (from the Definition 6.2.) that PT=CP At the same time, {Z Z Z Z and, using this fact, we construct a quaternion of the type In complete analogy with the complex numbers, where we have: now we have:
[TABLE]
Recall [] that the projective space CPn can be defined as the quotient
[TABLE]
Therefore, each point in CPn should be identified with the circle on From here, we obtain the familiar Hopf fibration (for ): . This logic fails for for where by analogy we formally should expect to have This quotient is not the Hopf map though. For the Hopf map is given by the quotient From here we are obtaining the already mentioned correspondence:
[TABLE]
In view of these results and, taking into account that we can formally write: , In addition, we have as well: and think about the correspondence as defining some analytic functions of and . In terms of such notations Eq.(6.22) can be rewritten as
[TABLE]
Using the multiplication table for quaternions [], p.39, it is possible to rewrite this result in more suggestive form as
[TABLE]
By comparing Eq.s(6.21) and (6.26) we identify with and with This identification can be extended by identifying and with and and and with and In literature on twistors [113] the spinor language is used. To comply with the standard twistor notations we should make an identification : on the one hand and, on the other. In terms of such notations Eq.(6.26) acquires the following form [114]:
[TABLE]
Here and are the two-component spinors.
**Definition 6.3. **Eq.(6.27) is known in twistor literature as incidence relation. It connects points in R4 (or Euclidean) space with points in PT.
By changing the matrix describing the Euclidean space to that describing the Minkowski space
[TABLE]
leads to the appropriately changed determinant: and to the incidence relation for the Minkowski space .
From previous discussion it is clear that the incidence relation, Eq.(6.27), could be used for both the Euclidean and the Minkowski spaces if we use the complexification: leading to the matrix defined by
[TABLE]
in which and should not be treated as complex conjugates of and (unless otherwise specified). In terms of such notations the incidence relation, Eq.(6.27), can be rewritten as
[TABLE]
Geometrically, this is the system of equations for two hyperplanes. In the language of algebraic/projective geometry they can be conveniently rewritten as
[TABLE]
The question of interest is: Under what conditions these two planes intersect? Why should we be interested in this question? Because we would like to connect just described twistor formalism with that presented in previous subsection. To do so, following [] we consider matrix M of the type
[TABLE]
Its 1st and 2nd rows are made of vectors and respectively. All Plücker coordinates are obtainable now as determinants of columns and of M. In view of Eq.(6.16), the matrix **M **can be used in the following matrix equation
[TABLE]
In the system of linear equations, Eq(6.30a), one can identify Z Z Z Z3 with the vectors , , in Eq.(6.32). Then, the matrix M acquires the following form
[TABLE]
Since the Plücker coordinates for this matrix are: the Plücker embedding condition, Eq.(6.18), acquires the following form
[TABLE]
and, is trivially satisfied.
Next, we would like to explain how just obtained result is connected with the intersection of two hyperplanes, Eq.s(6.30b), in PT=CP3. Following [115], page 83, consider a hyperplane through the origin O in CP It is described by the system of two equations
[TABLE]
The description of this plane remains unchanged if instead of we would use where is an arbitrary parameter. This observation makes the system of Eq.s(6.30a) equivalent to the system of Eq.s(6.35). A symmetrical set of coordinates is obtained by defining the six expressions Use of these expressions in Eq.s(6.35) allows us to eliminate successively resulting in the following system of equations
[TABLE]
Elimination of , and from these equations results in the Plücker relation, Eq.(6.18c), and the reminder equation is the equation for a complex projective line in CP3. So that, indeed, two planes, Eq.s(6.30b), are intersecting in a line. More details can be found in [116], pages 141-144. In view of Eq.s(6.15)-(6.17) six Plücker coordinates , represent a point in CP5 so that Eq.(6.18c) represents a complex quadric in this projective space CP5[116]. Following [] it is useful to rewrite this quadric in terms of new variables
[TABLE]
Use of these variables converts Eq(6.18c) into303030In the context of the Lie sphere geometry this equation is discussed in connection with Eq.(7.15) of section 7.
[TABLE]
Since are complex we can adopt them for physically relevant situations. We begin with
Definition 6.4. The of the origin is defined (according to Eq.(4.1)) as
[TABLE]
The of - the are subject to the constraint
[TABLE]
The Lorentz transformation sends the set of generators of into another set of generators of \Gamma\ preserving Following Penrose [117] we define the group as the group of all conformal maps of the compactified complex plane, that is of to itself. The connected component of the identity in consists of the orientation -preserving conformal maps of given by
[TABLE]
Details are given in []. From this reference we find that the homomorphism of the Lorentz group into is . The quadratic form defined by Eq.(6.39a) is obtainable from the quadric defined in Eq.(6.38) if we require The geometrical and topological meaning of these two extra terms is associated with conformal symmetry typical for the Lie sphere geometry to be described below. In the meantime, at this moment, we restrict ourselves by the concepts which were already in use in mathematical physics literature. For this purpose we recall some results from section 1. In it we mentioned that locally the static Einsteinian space-times are Minkowskian. Let be some 3-manifold so that topologically all static Einsteinian spacetimes are of the form [118], where R represents time. The positivity of mass theorem in general relativity, used in [], superimposed with use of the Yamabe theorem [], allows us to replace by The compactification requirement causes us to replace this manifold by But analytically is Accordingly, Eq.(6.38) does contain information about compactification. We already know that the twistor space is which is the complexification of whose compactification is Thus, we end up with the compactified Minkowski space and the compactified Euclidean space Both are described by the Kleinian quadric [] defined by Eq.(6.38). The sphere is described by
[TABLE]
while the (compactified) Minkowski space by313131Eq.(6.41b) is discussed in terms of the formalism of the Lie sphere geometry in section 7, in connection with Eq.s (7.10) and (7.20).
[TABLE]
**Definition 6.5. **The null cone of the origin of the compactified Minkowski space is described by Eq.(6.41b). The symmetry group leaving this null cone invariant is the conformal group
This is the largest symmetry group of the hydrogen atom.
We have reached these conclusions using arguments entirely different from those developed in the group-theoretic [] and high energy physics[83] literature. Use of the concepts of Lie sphere geometry to be discussed in section 7 provides us with solid theoretical framework for dealing with just discussed problems. At the same time, the obtained results are sufficient for finishing our study of Huygens triviality of the stationary Schrödinger equatioin.
6.3. Harmonic analysis and Huygens’ triviality
In Fock’s paper [102] the stereographic projection was used to relate and . The momentum space was one point ()compactified to 3-sphere living in If N is the North pole of the 3-sphere, N:=(0,0,0,1), we define U
Definition 6.6. The stereographic projection is the map U defined by
[TABLE]
Its inverse is defined by
[TABLE]
The map can be extended to defined as follows. By relating the North pole N for the 3- sphere to the compactification point all is covered. Since just defined extended stereographic projection provides a bijective correspondence between points on and points of the compactified momentum space . Using such stereographic projection Fock transformed the integral equation Eq.(6.1b), defined on into the integral equation on As result, he demonstrated that solutions of the integral equation inside are harmonic functions. That is they are solutions of the 4-dimensional Laplacian. The D’Alembertian, Eq.(2.6b), is obtainable from the Laplacian via formal replacement discussed in subsecton 6.2.1. Although Fock did not discuss the case in his paper, he did mention about the relevance of the Lobachevsky (that is hyperbolic) space for the description of case. At present, it is possible to find in literature solutions describing case. The detailed calculations are given, for example, in []. Although these calculations formally solve the problem, they are not revealing its mathematical essence. This essence was described already in Fock’s paper [] but was left undeveloped. Only recently Fock’s remarks were put into plausible mathematical form. In 2008 in the paper by Frenkel and Libine [122] in Advances in Mathematics entitled ”Quaternionic analysis, representation theory and physics” Fock’s results for both and were rederived. As result, it is sufficient to use only results by Frenkel and Libine for proving Huygens’ triviality of the Schrödinger equation (with time-independent Hamiltonian).
At this point it is logical to present a condensed summary of Fock’s results along with their subsequent improvements. Fock noticed that Eq.(6.1b) when stereographically lifted to the 3-sphere is looking very much the same as the Poisson formula for the harmonic functions inside the circle. It should be noticed, though, that Fock was not referring to the Poisson formula explicitly. The description of this formula employs the standard complex analysis [123]. Its derivation begins with the use of the Cauchy-Riemann equations. Recall, that these are defined as follows. Let and be some analytic functions satisfying the Cauchy-Rieman equations
[TABLE]
Since
[TABLE]
both and are harmonic functions. That is they satisfy the Laplace equation: The crucial theorem leading to the Poisson formula exploits this fact.
Theorem 6.7.* Let be the harmonic function in some simply connected domain of the complex plane C. Then, there is a regular in function such that Function is determined with help of with accuracy up to purely imaginary constant.*
The Poisson formula can then be defined as solution of the Dirichlet problem: Find a function harmonic in and coinciding with some prescribed function at the boundary of . Since any connected domain can be converted into a disc of, say, radius R, the solution of the Dirichlet problem can be represented via the Poisson integral formula as follows
[TABLE]
Let now R in Eq.(6.43), then this formula acquires look very similar to formula Eq.(9.19) of Fock’s paper []. Nevertheless, Fock’s Eq.(9.19) cannot be qualified as the Poisson formula. The Poisson formula is just a modification of the Cauchy integral formula []
[TABLE]
In previous subsections we demonstrated how quaternions replace complex numbers when The question arises: Can the rest of results of complex analysis be rewritten in the quaternionic language in R The development of the theory of quaternionic functions was initiated by Rudolf Fueter. In 1935, the year of publication of Fock’s paper, he produced an exact quaternionic counterpart of the Cauchy integral formula [],[]. However, it took another 73 years before this formula was converted into quaternionic analog of the Poisson formula. This happened only in 2008. Frenkel and Libine [] used the quaternionic version of the Poisson formula for demonstration that it correctly reproduces both the discrete and continuous spectrum of the hydrogen atom. The way of derivation of these spectra is opposite to that used by Fock. Specifically, Fock used the stereographic projection to bring the 3 dimensional integral equation, Eq.(6.1b), into 4-dimensional Poisson -like form. Furthermore, the wavefunction satisfying the Laplace Eq.(9.12) of Fock’s paper is not the same as that obtainable from his 4-dimensional Poisson-like integral equation. For the integral Eq.(6.1b), when converted into the 4-dimensional form, acquires the following look [125], page 83,
[TABLE]
to be compared with Eq.(6.43b). Here is known volume element of the 3-sphere. The X coordinates are X={ These are given by Eq.s(6.42b). The relation between in Eq.(6.44) and in Eq.(6.1b) is given by [], page 83,
[TABLE]
It would appear that the analogous relation for formally would solve the Hyugens triviality problem in view of the transformation rules defined in section 4. This is not the case, however. To explain the existing problem it is useful to briefly discuss the case first, since it is easier. In his paper Fock introduces the stereographic projection by defining coordinates and such that + etc. In addition, he introduces in a rather arbitrary fashion still another set of coordinates
[TABLE]
along with the 4-dimensional Laplace equation
[TABLE]
Since the harmonic function satisfying this equation is not a solution of the integral Eq.(6.44) (surprisingly), clearly Eq.(6.44) cannot be qualified as the Poisson formula adopted to four dimensions. At the same time, if in Eq.(7.46a) is the radius of the 3-sphere (eventually ) then, following Fock, it is possible to represent the harmonic function as
[TABLE]
Following Cordani [], page 90, we would like to rewrite Fock’s results a bit differently. Thus, let hl(X) be a homogenous harmonic polynomial of degree in R That is and Then, if
[TABLE]
leading to
[TABLE]
Let where is the variable unit vector on R In 4 dimensions n=e.g.see Eq.(6.46a). The normal derivative of on the unit sphere is given by Application of Green’s third identity [121], page 333, and [], page 87, permits us then to rewrite the l.h.s.of Eq.(6.44) with help of just defined results in the following form ()
[TABLE]
Comparison with the r.h.s. of Eq.(6.44) yields the discrete portion of the spectrum, for the hydrogen atom. To get the continuum portion of the hydrogen atom spectrum requires us to make some redefinitions of the already obtained results. It is convenient to represent the continuum and discrete results side-by-side. Specifically,
[TABLE]
Here stands for scalar product in space of Euclidean signature while stands for scalar product in the space of Lorentzian signature. In particular, in the last case we have
[TABLE]
so that the result should be looked upon either as an option or as an option Therefore,
Use of stereographic projection is defined by [] by analogy with Eq.s (6.42a),(6.42b) leads to the following identifications:
[TABLE]
Accordingly, for Eq.(6.44) we obtain:
[TABLE]
Here for and for
For the Laplace Eq.(6.46b) is replaced by the D’Alembert Eq.(2.6a). Accordingly, Eq.(6.47a) is now being replaced by [121]
[TABLE]
with N is being a real number. Since the solution of Eq.(6.51) presented in [121] contains many gaps in logic it is not going to be discussed any further in this paper. However one should keep in mind that it is this solution which is being used in solving the integral Eq.(6.50b). Thus, the harmonic function solving the Laplace Eq.(6.46b) for is being replaced by where the function is solution of the D’Alembert Eq.(2.6a) (or Eq.(6.51)) with the boundary conditions on . As in the case of E the spherical harmonic solves the integral Eq.(6.50b) for . However, the solution protocol for is very different. For us it is essential that the function is the solution of D’Alembert equation in the Minkowski spacetime having the hyperboloid of two sheets as boundary. Eq.s (6.45a) and (6.45b) are being replaced now by
[TABLE]
Suppose that solution of Eq of the integral Eq.(6.50b) is found, then the solution of the D’Alembert Eq.(6.51) is given by In the case the stereographic coordinates X={ had been artificially extended to so that if defined on the 3-sphere is the solution of the integral equation, Eq.(6.50), is the solution of the Laplace Eq.(6.46b). Analogous replacement in case produces solution of the D’Alembert equation from solution of the integral Eq.(6.50b). When just described extension of stereographic coordinates is combined with Eq.s (6.52) such hyperbolic analog of solves the problem of Hadamard triviality. This follows in view of the fact that the Schrödinger operator L̃[ Eq is obtainable from the D’Alembert operator L[ via sequence of non-singular transformations of independent variables (stereographic projection and extension). Thus, the Schrödinger operator for the time-independent Schrödinger equation is Hyugens’ trivial. At the physical level of rigor the same result was obtained by I.Bars [99] using two time formalism,
Frenkel and Libine solution of hydrogen atom model resulting in spectrum for bound and scattering states [] improve and considerably simplify just described results since the quaternionic Poisson formula is relating the harmonic functions inside and with solutions of the stationary Schrödinger Eq.(6.1) for both and as required.
7. Physical uses of Lie sphere geometry.
7.1. Bird’s view of the two-times formalism of I.Bars
In [83] the program of two -times formalism was outlined while in [99] detailed calculations illustrating general principles are presented. In short, the main idea is to replace 3+1 Minkowski spacetime by more general (fundamental) 4+2 spacetime having 4 space and 2 time variables so that the signature of this a spacetime is (+,+,+,+,-,-). Evidently, the symmetry of such spacetime is either O(4,2) or SO(4,2). In such an extended spacetime there could possibly exist one fundamental model casting countable number of ”shadows”(projections) into much more familiar Minkowski spacetime. Each ”shadow” is perceived as known distinct particle physics model. Under an umbrella of 2 -times formalism, it appears, that apparently different physical models actually are having the same origin in 4+2 spacetime. Mathematically speaking, the idea is to distribute all particles(and their bound states) of high energy physics into equivalence classes in accordance with the basic (fundamental) models living in 4+2 dimensional spacetime. Since the already made calculations connect massless relativistic particles with the massive ones (even accounting for the spin), the concept of a mass of particle and significance of the Higgs mechanism for mass generation loose their central importance in this formalism. Furthermore, the same formalism connects, for example, the relativistic particles (massive or not) with the non relativistic massive particles, the extended (bound) systems, such as hydrogen atom, with the massless relativistic particles, etc.
In section 6 we demonstrated how the nonrelativistic Schrödinger equation for hydrogen atom is connected with the D’Alembert (in 3+1 dimensions) or the Laplace (in 4 dimensions), that is with equations describing the massless and spinless relativistic particle. In section 5 we demonstrated how relativistic equations for massive particles can be reduced to those whose masses are zero. We were guided by the Hadamard ideas of Huygens triviality casting all Huygens-trivial equations (e.g. read section 4) into the same equivalence class. In another paper [126] Bars connected the hydrogen atom with the massless relativistic free particle and with the massive 3-dimensional harmonic oscillator. Thus, he effectively demonstrated Huygens triviality of the hydrogen-atom and harmonic oscillator. His work develops results of [127] published previously in which Bars (with collaborators) raised a question: Which of 2 times in the 2 time formalism is the familiar time coordinate? Unfortunately for his project, no clear-cut answer to this question was given either in [127] or in the basic reference [83] summarizing all accomplishments of the 2-times formalism. Accordingly, the 2-times formalism did not gain much popularity among the high energy physicists. Nevertheless, in our opinion, upon development, e.g. based on results of this work, there could be a chance to make formalism developed by Bars more viable.
We initiate this process of further development with the discussion of the relationship between the dynamics of hydrogen atom and harmonic oscillator. Being driven by the same problem of regularization of singularities (e.g. read again comments to Eq.(6.13) for the Kepler problem), Kustaanhemo and Stiefel (K-S) using the Hopf mapping (section 6) transformed the 3 dimensional classical dynamics of the Kepler problem into the 4-dimensional problem of classical dynamics of the harmonic oscillator on Their work was extended by many authors both at the classical and quantum levels. At the quantum level readable account of uses of the K-S transformation converting hydrogen atom problem into harmonic oscillator problem is given in []. An independent and much simpler treatment of the same conversion problem was given in [130] by Chen. By utilizing results of [130], Chen demonstrated in [] the isomorphism between the largest symmetry group SO(4,2) of the hydrogen atom and the group SU(2,2) used in physics of twistors [25]. Later on he connected group-theoretically hydrogenic, oscillator and free-particle massless relativistic systems in []. In his derivations Chen used only one-time formalism. His results do not use twistors or compactification of the Minkowski space. Thus, they cannot be immediately linked with discussions involving the Lie sphere geometry-a tool essential for generalization of these results. At the same time, the formalism by Bars, perhaps if further developed, can accommodate the concepts of Lie sphere geometry.
7.2. Lie sphere geometry. Fundamentals
The examples discussed in previous subsection indicate that 2-times formalism invented and developed by Bars in some instances can be entirely replaced by more familiar 1 time formalism. The examples of such a replacement can hardly be generalized though. At the same time, the diversity of results obtained with help of the 2 times formalism is appealing. Thus, in this subsection we would like to suggest a reliable direction enabling us to replace questionable 2-times formalism by mathematically well developed formalism of the Lie sphere and the Möbius geometries. In physics literature the formalism based of utilizing results of Möbius geometry is known as ”geometric algebra”.
We begin with the description of Lie geometry since the Möbius geometry (and, therefore, geometric algebra) is more restrictive. The Lie geometry is the geometry of generalized oriented hyperspheres living in the compactified Euclidean space =The elements of this geometry are:
a) Oriented hyperspheres. These are familiar from section 4 spheres of finite radius with center so that
[TABLE]
The hypersphere divides into 2 parts. If we denote one of these parts as ”positive”, then another part is ”negative”. By introducing such a distinction we are introducing the oriented hyperspheres. In such a case a given is replaced by These results generalize 3 dimensional result, Eq.(4.25), needed for definition of canal surfaces and Dupin cyclides.
In higher dimensions the cyclides of Dupin are being replaced by the Dupin hypersurfaces. Such hypersurfaces can live in Euclidean and hyperbolic spaces [].
**Definition 7.1. **Hypersurfaces with constant principal curvatures are Dupin hypersurfaces.
To distinguish the oriented hyperspheres analytically it is convenient to introduce signed radius which is telling us whether we should consider the inward or outward field of unit normals. According to convention, the positive radii are assigned to hyperspheres with the inward field of unit normals while the negative radii are assigned to hypersheres with the outward field of unit normals. By doing so we just had introduced a bijection between the hypersurfaces of non-vanishing radius and tuples
[TABLE]
Here
b) Oriented hyperplanes. A hyperplane in is characterized by the equation
[TABLE]
with a unit normal and Evidently, the tuples and represent the same hyperplane. As a hypersphere, it divides into two half spaces. By declaring one of the two half spaces to be positive, we are getting the notion of oriented hyperplane. Thus, by analogy with hyperspheres, we obtain a splitting: for any .
c) Points. Points are hyperspheres of zero radius.
d) Infinity. Infinity point is making compactified, that is =
e) Contact elements. A set of all hyperspheres through which are in oriented contact with and with one another thus all sharing normal vector n at x.
Remark 7.2. The connection/replacement with/of 2 times formalism by Bars follows now from the observation that all just described elements are modelled as points, respectively lines, in the projective space P(R with the space of homogenous coordinates R Details follow below.
From now on, for simplicity, we shall only discuss the case: . Following Suris [] we equip the space R4,2 with 6 linearly -independent unit vectors whose scalar product is defined as
[TABLE]
To exemplify the properties of Lie sphere geometry in the most efficient way, it is convenient to make the following redefinitions:
[TABLE]
In terms of just made definitions we redefine the objects of Lie sphere geometry as follows
Oriented sphere with center and signed radius
[TABLE]
Oriented plane defined by with and
[TABLE]
Point
[TABLE]
Infinity:
[TABLE]
Contact element:
[TABLE]
**Remark 7.3. **In the projective space P(R the first four types of elements are represented by the points. These are the equivalence classes of the above vectors with respect to the equivalence relation R
**Remark 7.4. **A contact element in Rn is an isotropic line in P(**R**This line is defined in Eq.(7.11) below.
To be specific, we need to introduce the Lie quadric. This is accomplished as follows. First, we define the set of isotropic vectors in Rn+1,2 via
[TABLE]
then the Lie quadric Qn+1,2 is defined as
Remark 7.5.
For the Lie quadric just defined coincides with that defined by Eq.(6.41b) representing properties of the compactified Minkowski space
It was obtained via twistor formalism in previous section thus hinting at the conection between the twistor and the Lie sphere formalisms. That this is indeed the case will be explained below, in subsection 7.4.
Going back to contact elements, choose now two hyperspheres (or spheres if ) Rn+1,2 such that Then, these spheres are in oriented contact with each other. An elementary proof can be found in [], page 15. Using this fact, we can define a line in projective space . For this purpose let and be some real numbers. Using these numbers we define a linear combination a projective line. Using it, we obtain:
[TABLE]
since both and are isotropic vectors. If and in R4,2 represent two spheres in oriented contact, then the line in trough the corresponding points is isotropic as just demonstrated. It lies entirely on the Lie quadric There are no projective subspaces of higher dimensions completely contained in This is proven in[134], page 17. In just described formalism planes are spheres of infinite radii and points are spheres of zero radii. The Definition 4.7. can be restated now as follows
**Definition 7.6. **The essence of Lie geometry lies in study of projective transformations of (for ) leaving the Lie quadric invariant. The group of such transformations is factor group .
This quotient is just the higher dimensional extension of the earlier discussed quotient describing a homomorphism of embedding of the Lorentz group ) into group introduced in connection with Eq.(6.40). These transformations are preserving the isotropy property described by Eq.(7.10). Furthermore, the (non) vanishing of or of component of a point in is not invariant under a general Lie sphere transformation leading to absence of distinction in this geometry between the oriented spheres, oriented planes and points.
Presented background is sufficient for reading of [135], chapter 15, where many additional facts about the Lie sphere geometry are nicely explained.
7.3. Möbius geometry and geometric algebra. Fundamentals
In this subsection we connect the Lie sphere geometry with the Möbius geometry known in physics literature as ”geometric algebra” [136]. To demonstrate interconnection between Möbius geometry and geometric algebra we have to provide some basics on Möbius geometry first. This is easy to do since Möbius geometry is just part of the Lie sphere geometry. Möbius geometry studies subgroups of Lie sphere geometry preserving subsets of with vanishing Thus, it deals with non-oriented spheres, non oriented planes, points the infinity point compactifying to The elements are in one-to-one correspondence with the points on the projectivized light cone where
[TABLE]
to be compared with Eq.(7.10). Points correspond to points of with non-vanishing while the point corresponds to the only one point of with vanishing component.
**Remark 7.7. In view of the Remark 7.2., **by comparing Eq.s(7.10) and (7.12) it is clear that results of Bars can be redone only with help of the Lie sphere geometry. Nevertheless, Möbius geometry is used in geometric algebra which found its way into physics for some time [136]. This fact deserves some further discussion. In particular, we begin with the following definition
Definition 7.8. The essence of Möbius geometry lies in the study of properties of nonoriented (hyper)spheres invariant with respect to projective transformations mapping points to points, The group of such transformations is isomorphic to . It is the group of Lorentz transformations of Rn+1,1 preserving the time-like direction. Every conformal diffeomorphism of is induced by the restriction of Möbius transformation to . Further details can be found in [].
**Remark 7.9. Although the Dupin cyclides **are described by the Lie sphere geometry in 3 dimensions [], additional studies, e.g. see for example [], demonstrated that in 3 dimensions use of the Möbius geometry and, hence, the geometric algebra is sufficient in the sense described in the next subsection. This algebra describes all conformal motions. In [] it is demonstrated how the Dupin cyclides emerge as sets of orbits of conformal motions.
7.4. Crown achievement of Sophus Lie- discovery of the isomorphism between the Lie sphere and Plücker line geometries.
The essence of 2-times formalism invented by I. Bars and described in his book [] is summarized in Fig.7.7. of this book. In it, the 2-times formalism is placed at the center of this figure while the twistor formalism is presented in the upper left corner. It is depicted as some kind of corollary of the 2-times formalism. We have already explained that the 2-times formalism is nothing else but the Lie sphere geometry formalism. Now, following ideas of Sophus Lie [], we are going to demonstrate that the Lie sphere formalism is isomorphic to the Plücker line formalism. Plücker formalism was discussed in section 6 in connection with twistors. That such an isomorphism is possible, is hinted in the Remark 7.5. Now we provide the details. By describing the Lie sphere-Plücker line correspondence we are establishing the isomorphism between the Penrose twistors and the 2-times formalisms. Since the twistor formalism is associated with the most of the exactly integrable systems originating from all kinds of reductions of the Abelian and non Abelian-Yang Mills gauge fields [], it is hoped, that in view of this isomorphism it is sufficient to study properties of these gauge fields in order to address all problems of high energy physics, including those pertinent to the Standard Model and gravity. Furthermore, the unexpected connection between the gauge-theoretic (Floer) and Schrödinger’s formalisms noticed in subsection 4.2.3. now acquires an independent support.
We develop Plücker’s line geometry by analogy with the Mobius and Lie sphere geometries. For this purpose, we recall first Eq.s(6.16)-(6.17) describing representative elements of the exterior algebra as well as Plücker embedding of this algebra into complex projective space. We are adopting general formalism to four dimensional complex (twistor) space C By treating this space as vector space we have {e0,…,e3} as the basis for C Accordingly, we also can construct a complex six dimensional space of bivectors:
[TABLE]
Notice, that in the case of the Lie sphere geometry we used also six dimensional vector space {Therefore, by analogy with the Lie sphere case, we define the scalar product according to the rules []:
[TABLE]
Notice that thus far we nowhere used the fact that the underlying vector space is complex. Thus, following [] we shall initially treat it as real. Then, based on the multiplication table given by Eq.(7.14), we conclude that the signature of such a space is (3,3). Accordingly, In accord with Eq.(7.10) we introduce the set of isotropic vectors
[TABLE]
The Plücker quadric Q3,3 is now defined by P( It was obtained in section 6, Eq.(6.38). Use of transformations given by Eq.(6.37) brings us back to the well known Plücker embedding result, Eq.(6.18c). Notice that the standard Plücker relation, Eq.(6.18c), does not require use of complex variables. This happens to be essential as we shall explain momentarily. For this we recall that Lie discovered line-spherical geometry correspondence in 1870 following ideas of Poncelet, Plücker and Darboux. All these ideas, including contributions of Sophus Lie, are presented in the book by Klein [144] while the detailed account of the Lie Sphere geometry, Dupin cyclides, etc. is presented in the book by Blaschke []. In particular, in addition to the sketchy presentation of the line-sphere geometry correspondence given in [141], Felix Klein provides accessible and comprehensive discussion on the line-sphere geometry correspondence in his book [144], pages 262-274. Since books by Blaschke and Klein are both in German, for readers convenience, we provide only summary of what was accomplished which we are seamlessly connecting with the previous material. We refer to just mentioned books for further details.
Consider a single sphere in R Analytically it can be described as follows
[TABLE]
In addition to parameters it is necessary to introduce the radius of the sphere via
[TABLE]
These parameters are to be considered as coordinates in the space of spheres. In addition, it is convenient to introduce the homogenous coordinates
[TABLE]
enabling us to insert the sphere into projective space. To do so, we insert Eq.(7.18) into Eq.(7.17) with the result:
[TABLE]
Planes (spheres of infinite radius) and points (spheres of zero radius) are included in this parametrization. For points we have to put in Eq.(7.19) while for planes we require To represent Eq.(7.19) in terms of the Lie sphere quadric Q4,2 (defined after Eq.(7.10)), it is sufficient to represent the combination as Then, Eq.(7.19) is replaced by
[TABLE]
which is the same as previously obtained Eq.(6.41b). Questions remains: a) Is it possible to convert the quadric Q4,2 into Q3,3 by some kind of nonsigular and real transformation? b) If this is not possible, can the same result be achieved by using complex transformations? In other words : Is it possible to prove that the projective groups and are isomorphic? The answer is ”NO” in the domain of real numbers and ”YES” in the domain of complex numbers. Explicitly, consider Plücker’s result, Eq.(6.18c), and perform the following transformation (attributed to Lie)
[TABLE]
Substitution of this result into Eq.(6.18c) brings us back to Eq.(7.19). Although algebraically this substitution looks very simple, geometrically, it is at the heart of the line-spherical geometry correspondence. All missing details of this correspondence are given in the book by Klein [144]. Thus, all results of 2-times formalism developed by Bars should be obtainable from the twistor formalism.
7.5. Some physical implications of the Lie sphere geometry
7.5.1. Lie sphere vs Möbius geometry. Important facts
To study physical implications, we need to extend results describing the line- Lie sphere geometry correspondence presented in Klein’s book [144]. This is needed because the Dupin cyclides were not discussed by Klein in the context of this correspondence even though they were discussed in his book. This is excusable in view of Definitions 4.4 and 4.5. However, in describing the Lie sphere geometry we introduced isotropic vectors via Eq.(7.10) and described the oriented contact between two spheres in Eq.(7.11). A simple calculation done in [], page 15, indicates that the oriented contact condition is equivalent to
[TABLE]
We would like to test this result using the simplest Dupin cyclide described after Definition 4.5. Using these results, we obtain: . At the same time . Thus, we just demonstrated that whenever there is an oriented contact between two spheres, there is a Dupin cyclide associated with such a contact. The contact between two spheres takes place along some curve on the surface of the Dupin cyclide. Definition 4.5 indicates that the oriented contact between spheres is the result of coincidence of envelopes originating from two different canal surfaces. According to Maxwell, only Dupin cyclides are being formed in this way. Evidently, the Lie sphere geometry permits both (or and implying that this geometry permits the existence of both the canal surfaces and the Dupin cyclides. At the same time, the spine curve c(t) could be closed or not. Accordingly, the most primitive canal surfaces are cylinders, cones of revolution and tori. It happens that use of Möbius geometry and, in particular, the inversion transformation, is sufficient for generation of all known Dupin cyclides [] starting with cylinders, cones of revolution and tori. For the case of tori of revolution, a very accessible proof is given in [147], Proposition 20.36. This fact allows trivially reobtain the results by Friedlander [65,70] and Sym [74] starting with cylinders, cones of revolution and tori as solutions of the D’Alembert Eq.(2.6)323232More details are given below..
The Dupin cyclides as much as D’Alembertian, Eq.(2.6), are invariants of the Lie sphere transformations, e.g. read again the discussion around Eq.(7.11), while the cyclides deform under Möbius transformations. Complicated Dupin cyclides are formed from the simplest ones via application of succession of Möbius transformations. Since the line-Lie sphere geometry correspondence does not survive under use of Möbius transformations, it appears, that Möbius transformations are playing an auxiliary role. This conclusion is reinforced by the observation that in Möbius geometry the condition is not the contact condition given Eq.(7.22). In Möbius geometry we have instead [80]:
[TABLE]
which is the condition of two spheres to intersect orthogonally. Notice that the above condition does not require use of the parameter. Surprisingly, these conclusions are wrong! To explain why, following [140], we need to describe the conformal motion first.
7.5.2. Fundamentals of conformal kinematics
As in physics, everything begins with mechanics and mechanics begins with kinematics. The Euclidean kinematics is composed of rotations around some axis followed by translation along the same axis. It is a special kind of conformal kinematics. The kinematics of Möbius geometry is based on the known fact that it is an isometry of hyperbolic space [81]. This means that in such geometry circles are geodesics. They are analogs of straight lines in Euclidean geometry. Using this observation, we cover a plane R2 (or sphere with a network of mutually orthogonal circles (instead of Euclidean plane being covered by a lattice made of mutually orthogonal straight lines). On thus formed Euclidean (or non Euclidean) lattice we can move along geodesics, that is along the segments of straight lines or along the segments of circles. In both cases the result of elementary motion along the segment followed by motion along the orthogonal segment can be interpreted as rotation. In going from two to three dimensions, in Euclidean geometry, the planar square lattice is replaced by three dimensional lattice. In hyperbolic space we have to cover** R3** by the network of mutually orthogonal spheres. This explains at once the significance of Eq.(7.22) in Möbius geometry. By avoiding technicalities clearly explained in [140] we only notice that any three dimensional conformal motion can me made of two commuting elementary motions. These are analogs of Euclidean translation along the axis and rotation around the same axis (screw motion). Consider now the most elementary Dupin cyclide-the torus. For the torus the meridians and parallels are mutually orthogonal. Both are circles. But circles are crossections of spheres! Accordingly, if we embed such a torus into R3 (or, better, in foliated by mutually orthogonal spheres in such a way that meridians and parallels of the torus are the crossections of the associated spheres, then we can think about the conformal trajectories originating on the surface of such a torus. In view of results of previous subsection, any nontrivial cyclide of Dupin is obtainable from torus via operation of Möbius inversion. This means the following:
a) The network of mutually orthogonal circles on the torus will transfer to the network of mutually orthogonal circles on the cyclide.
b) The set of all orbits originating at the torus will form two dimensional surface -Dupin cyclide. If we look only at some orbits, they are developing at the surface of Dupin cyclide.
c) When an initial point x(0) is given, it determines both the orbit and the type of cyclide.
Remark 7.10. The above statement a) can be formulated as a theorem. Originally it was known as conjecture (attributed to Ulrich Pincall): Cyclides of Dupin are the only surfaces in Euclidean space on which two families of orthogonal circles lie. Although we explained why this is so above, the full proof was given by Thomas Ivey [148] in 1995. This proof explains why many researchers, e.g. read [140] or [149], use Möbius geometry for generating Dupin cyclides. Surely, the same results are achievable with the help of Lie sphere geometric methods as explained in [150],[151].
Just obtained results can be used for reobtaining in the most economical (physical) way both -the result by by Friedlander [70] discussed in section 4, and that for knots made of null fields discussed both in electrodynamics [152] and in quantum mechanics [14]. Details are given in the next subsection and Appendices B and C.
7.5.3. Knots and Dupin cyclides in quantum mechanics and electrodynamics
We provided an evidence for existence of Dupin cyclides in quantum mechanics and in electrodynamics in section 4 and, from different perspective, in Appendix C. We would like now to reobtain results by Friedlander using physical arguments. Following [14], we begin with the observation that the set of source-free Maxwell equations in the vacuum can be compactly rewritten as
[TABLE]
Here is the speed of light, is the Riemann-Silberstein vector involving both the electric E and magnetic H fields. By representing this vector as
[TABLE]
the above set of Maxwell’s equations is converted into
[TABLE]
where is + or - and . In plasma physics Eq.(7.26a) is known as ”force-free equation” while in hydrodynamics it is known as ”Beltrami equation”. Eq.(7.26a) was discussed in detail in our book [55] in connection with problems emerging in contact geometry. By applying the curl operator to both sides of Eq.(7.26a) and taking into account Eq.(7.26b) we obtain:
[TABLE]
This vector version of the Helmholtz equation should be compared with its more familiar version, Eq.(3.2), discussed in section 3. This (scalar) version is used both in quantum mechanics and in electrodynamics. The way to relate these two equations to each other is nontrivial [14]. In part, it is based on the useful identity
[TABLE]
In [14] it is shown how the vector, Eq.(7.27), can be restored from the scalar Eq.(3.2) in which By applying operation div to both sides of Eq.(7.26a) and by assuming that we obtain, div( Let be some trajectory on the surface Since by identifying we conclude that the ”velocity” is always tangential to the surface Since the vector field v is assumed to be nowhere vanishing, the surface can only be a torus T2. In the case if is rational number the field lines v on T2 are closed implying that these lines are forming torus knots. Using results of previous subsection and applying Möbius transformation we end up with torus-like knots wound around cyclides of Dupin (that is around the distorted tori). Depending on the type of Möbius transformation, different ”wrapped” Dupin cyclides will be formed. But this is exactly the result of the paper by Friedlander [70]!
**Remark 7.11. **Use of conformal transformation for generation of knots/links was suggested initially by Bateman and Cunningham in 1910. Their results were recently rediscovered and utilized for knot/link generation in [152]. See also [153]. No mention of Dupin cyclides, etc. was made in these references.
**Remark 7.12 **Results of ref.s [152] and [153] demonstrate how the torus knots/links can be generated from null and not null electromagnetic fields (discussed in some detail in in the Appendix B). Being technically permissible, such torus knot generation is not allowed for the non null fields on physical grounds.
7.6. Place of Lie sphere geometry in development of foundations of AdS-CFT correspondence
We conclude our paper with a brief discussion of the foundations of AdS-CFT correspondence. An excellent physical introduction to this subject is given in the recent book by Nastase []. Some mathematical aspects of this correspondence are discussed in our work []. Use of these references allow us to reduce our discussion to the absolute minimum. In physical and mathematical literature on AdS-CFT correspondence no mention of relevance of the Lie sphere geometry exist to our knowledge. Thus, what follows below is the fist attempt at elimination of this deficiency.
We begin with familiar examples, e.g. discussion of the simplest model of hyperbolic space H It is known in literature as the Poincare*′* disc model D2 Geodesics in this model are made of the horocycles. These are circular segments whose both ends lie at the boundary of the disc . The boundary is considered to be ”the spatial infinity”. By the appropriate choice of constants in the Möbius transformation given by the Poincare*′* disc model can be transformed into Poincare*′* half plane model of the hyperbolic space H2. Clearly, use of the inverse conformal transformation converts the half plane model back into the disc model. Let be some point inside such disc model. Successive applications of the Möbius transformation will propel this point to the spatial infinity after infinite number of iterations. For detailed examples, please, read []. This two-dimensional model is generalizable to higher dimensions where it is known as the hyperbolic ball model [81]. Thus, when interested in higher dimensions, we need to replace a combination (H by (H In all described models of hyperbolic geometry the boundary at infinity is playing an important role. This role can be seen already in two dimensional model (H In it, the deformations of lead to the Virasoro algebra -central object of the conformal field theory [], the Teichmüller spaces, etc. Thus, already at the level of two dimensions we are dealing with a kind of AdS-CFT correspondence. By analogy with two dimensions, it is expected that the hyperbolic-like interior- actually, the anti–de Sitter spacetime (AdS) (in higher dimensions)- is affected by (linked with) the conformal field theory (CFT) residing at the boundary. The higher dimensional AdS-CFT correspondence in higher dimensions cannot be rigorously developed based on the hyperbolic ball model though. This is because of the Mostow rigidity theorem. The Mostow rigidity makes deformations of conformal 3-sphere impossible to perform. Nevertheless, the mathematics of this model is serving as a gudeline for more realistic models as demonstrated by Frances []. The bottom line is the following. The hyperbolic space Hn is replaced by the Anti-de Sitter space (AdS)-the Einstein space of constant negative curvature. This space is obtainable from the vacuum Einstein equations with added (negative) cosmological constant [155]. The hyperbolic ball at infinity is now replaced by the space of conformally flat solutions Einn,1 of Einstein’s equations. Being conformally flat these are conformally equivalent to Minkowski spacetime, that is to say their conformal symmetry group is SO(4.2). In the Appendix D we rigorously demonstrate that
[TABLE]
Here is the boundary of the Anti-de Sitter space while is the Lie sphere quadric introduced after Eq.(7.10).
**Acknowledgement ** Authors are deeply grateful to Professor (Dr.) David Delphenich for supplying us with the original of Madelung’s paper, ref. [177], as well as with the English-German translation of this paper. Authors are also very thankful to the unknown referee for his thoughtful remarks which lead to considerable improvement of our presentation.
**Note added in proof. The **interplay beween the differential and the Lie Sphere geometries, the rigidity and quantization (and hence mass generation) in the context of the transversal knots (discussed in our earlier work, Ann.Phys.371 (2016) 77) is developed in great detail in two latest works by E.Musso and L.Nicolodi (Comm.Anal.Geom. 25 (2017) 209 and Nonlin.Anal. 143 (2016) 224). The connection between knots and the Regge-mass spectrum is discussed in our earlier work ( Int. J.Mod. Phys. A 30 (2015), id. 1550189-212).
Appendix A
Comparison between the Cauchy problems for parabolic
and hyperbolic equations
-
Phase velocity. Without loss of generality, we consider partial differential equations (PDE) involving just two variables and . Suppose that solution can be represented in the traveling wave form [45] as follows: Here is the velocity of the profile . More generally, in the case when we get It is called plane wave with wavefront normal to k, speed ( ), and profile .
-
Exponential solutions for simplest PDE’s and group velocity**. **When studying linear partial differential equations, it is convenient to specify the profile function in order to study the complex-valued plane wave solutions of the form
[TABLE]
Using such an ansatz in the diffusion equation leads to
[TABLE]
resulting in the dispersion relation Anticipating generalizations, this result can be rewritten as Therefore Eq.(A.1) can be presented in the form
[TABLE]
This representation allows us to introduce the group velocity cg as
[TABLE]
In the case of simplest Schrödinger’s equation ( the dispersion relation is given by
[TABLE]
Therefore, the group velocity is Upon restoring the usual system of units, we obtain: implying Accordingly, Using the De Broglie relation: we conclude that the group velocity** cg** coincides with the particle velocity in quantum mechanics.
Consider now the wave equation using the ansatz (A.1). We obtain instead:
[TABLE]
resulting in
[TABLE]
From here the phase velocity c is obtained as and the modulus of group velocity cg coincides with the phase velocity.
- Cauchy problems**. ** The Cauchy problem for the 1 dimensional diffusion equation is formulated as follows. For the diffusion equation
[TABLE]
defined on the whole line - whose solution at the time is find a solution for and The solution of this problem is given by
[TABLE]
Here , and
[TABLE]
The last two equations can be combined into
[TABLE]
The analogous problem for the Schrödinger’s equation is obtained now by replacing by in the above result. This fact is helpful for the following reason. Suppose that the initial condition is zero everywhere, except on some small interval . Because the exponential function is never zero, the integral, Eq.(A.9a), is not going to be zero, even for and arbitrary far from (). This can be restated as
Infinite propagation speed [45]. The initial condition affects the solution for all no matter how small is. Thus, heat propagates with infinite speed.
Evidently, the replacement of by in the above results is leading to the Schrödinger equation. But use of the same arguments as in the diffusion case are not changing this conclusion. This fact is one of the sources of Einstein’s spooky action at the distance.
Diffraction of light-main evidence of quantum mechanical behavior- uses time-independent Helmholtz Eq.(3.2) of the main text. This equation emerges in both quantum mechanics and in optics. However, if one is willing to study the Cauchy problem for the wave equation defined on -, it should be formulated as follows. In addition to the initial data one has to supply the initial velocity Thus, while for the well posedness of the Cauchy problem for the diffusion equation one needs just one initial condition, one needs two initial conditions for the wave equation. At the mathematical level no more comments are required. At the physical level, more comments are needed. These are originating from the fact that, say, optical waves are generated by atoms []. The two -level atomic systems - primary sources of photons-require for their description the time-dependent Schödinger equation which, is identical in its structure with the Pauli equation for spin 1/2 particle in varying magnetic field []. For such two-level system two initial conditions are required. In addition, technologically it is important to generate single polarized photons and to use the 50:50 beam splitter to cause such a single photon to be self-entangled []. Thus, even though a dispersion relation Eq.(A.7) originating from the wave Eq.(A.6) can be used for description of a single photon 333333E.g. via the de Broglie relation, in quantum optics polarization of photons is exploited essentially [16] causing use of the Pauli-type Schrödinger equation for their description. In such a case quantum phenomenon of entanglement [] coexist with the classical fact of finite speed propagation of wave signals.
Appendix B
Unimaginable universality of Hadamard premises
B.1. Brief review of Hadamard’s premises
There are 3 ”premises” formulated by Hadamard. E.g. read [162], pages 445-450. These are:
a) Major premise. Suppose we are observing events within the time interval In order to find a state at the moment we need to know a state at the time , then find a state at the time and, using this information, find a state at
Minor premise. Suppose within a short period of time there emerges some light disturbance localized in the vicinity of point then, at time this disturbance will be concentrated in the very thin spherical layer enclosing a sphere of radius centered at ;
c) Corollary. In order to evaluate the action of the initial light perturbation, located at the point O at , it is permissible to replace this perturbation by a set of perturbations emerged at time and distributed on the surface of the sphere centered at O and having radius
According to Hadamard, different writers identified Huygens’ principle either with a) (e.g. Feynman, as described in section 1) or with b) (e.g. read our section 4) or with c). In his early work on hyperbolic equations and Huygens’ principle Hadamard considered only such wavefronts which did not contain any singularities known as caustics. The major premise does not exclude existence of caustics. Therefore, what is considered by mathematicians as ”Huygens’ principle”, corresponds to Hadamard’s minor premise and to corollary. In such a form it is discussed in the main text, e.g. read Definition 4.2. Alternatively, read [68], page 138. In short, Hadamard restricted himself by studying of the wake-free waves. Situations when wakes are present (leading to caustics) is considered in [50]. Much more comprehensive and detailed treatment (but not mathematically well supported) is given in [163].The same but mathematically supported is given in [59]. Huygens triviality is defined in Definition 4.3. Following this definition, the Hadamard conjecture is formulated as follows:
Hadamard conjecture. Every Huygens operator is trivial.
For purposes of quantum mechanics it is almost always sufficient to use the results of Theorem 4.1. In this appendix we would like to explain that in the broader (but physically also very relevant) context Theorem 4.1. sometimes fails. This is so if, following [68], we reformulate Theorem 4.1. as
**Theorem 4.1.1. **If equation of hyperbolic type in conformally flat Minkowski (3+1) spacetime satisfies Huygens principle, then it is equivalent of the wave Eq.(2.6a).
The new element here is mention of conformal flatness. The wave Eq.(2.6a) is surely conformally invariant, e.g. read [82]. Based on this fact, a broader question can be posed
:
Hadamard-like conjecture**.** Is it true that the Huygens principle holds if the underlying equation is conformally invariant?
That is to say:Is the conformal invariance and the Huygens principle are equivalent statements (or interdepenent concepts)?
Stated still a bit differently, the above conjecture can be formulated as follows.
Hadamard-like conjecture. Second version.* Is it always true that the Huygens principle is obeyed whenever it is possible to relate the conformally invariant equation of the 2nd order to the wave Eq.(2.6a)?*
B.2. Conformal invariance versus the Huygens principle
In any given dimension , say, in , all conformal groups are classified in flat and curved spaces [164]. Many important additional details are given in [165] and [166]. Accordingly, all conformally invariant equations can be explicitly written for all Einstein spaces of general relativity. If this is so, will equations with different conformal groups lead to the same Eq.(2.6a)?
It is instructive to provide some details to answer this question. Following [166] we introduce the second order differential (Beltrami) operator as follows
[TABLE]
With help of this operator any 2nd order linear differential equation can be written as
[TABLE]
In the same metric conformally invariant 2nd order equation can be written (in 3+1 dimensions) as
[TABLE]
where is the scalar curvature. At the same time, using [165],[166] it is possible to prove the following:
Theorem B.1. Any spacetime with a given metric g possesses nontrivial group of conformal motions if and only if this spacetime is conformally equivalent to the Lorentzian spacetime.
The most general Lorentzian spacetime V3+1 is described by the metric
[TABLE]
Here is the positively definite matrix.
**Corollary B.2. **Any Lorentzian spacetime with nontrivial conformal group is conformally equivalent to the space for which the Ricci tensor
Corollary B.3. *Every 3+1 Lorentzian spacetime V3+1with nontrivial conformal group admits the 2nd order partial differential equation *
[TABLE]
satisfying Huygens’ principle.
**Remark B.4. ** It is clear that for a special choice of functions and Eq.(B.5) will coincide with the wave Eq.(2.6a). However for other choices these two equations are not always conformally equivalent since they might involve different groups of conformal motions. This will be further explained below.
Remark B.5. We brought to our reader’s attention Eq.(B.5) because of the remarkable paper by Ward [167] allowing us to connect the results of this appendix with those in the main text.
In section 4 we introduced and discussed the progressive wave solutions of the wave equation as well as the Dupin cyclides. In section 7 we reobtained Friedlander’s results [70] in much simpler way by using methods of conformal and contact geometries. Now, using these methods, we want to rederive Ward’s results having additional purposes in mind.
In his paper Ward was not discussing the Huygens principle or its connection with the conformal invariance. Ward’s purpose was to demonstrate that if instead of the wave Eq.(2.6a) one considers the 2nd order wave-like equation with metric taken from the plane-wave gravitational background, one still can apply the progressive wave solution method by Friedlander to this equation and to reobtain all types of Dupin cyclides discussed by Friedlander. Based on results of this appendix and those in sections 4 and 7 of the main text, Ward’s goals can be restated as follows:
Question B.6. If the wave Eq.(2.6a) admits progressing waves solutions resulting in Dupin cyclides, is the same is true for the Eq.(B.5)?
In our discussion of Eq.(B.5) we have not touched the subject of plane gravitational waves while Ward was not looking any further at the origins of Eq.(B.5) which he interpreted in terms of the plane gravitational pp waves background. We would like now to complement Ward’s results with missing details. By doing so we simplify Ward’s results as well and explain why plane gravitational waves are fundamentally relevant for our paper.
B.3. Conformal groups and Lie sphere groups, Rainich, Misner and Wheeler
geometrodynamics of coupled Einstein-Maxwell fields, the AdS-CFT
correspondence and Penrose limit
We begin our discussion with results summarized in [168]. From chapter 3 of this book it follows that the flat Minkowski metric
[TABLE]
can be converted into
[TABLE]
with help of substitutions From chapter 17 of the same book we obtain the metric for the pp-gavitational waves (obtained for the first time by Brinkman in 1925 and interpreted in terms of gravitational waves by Peres in 1959):
[TABLE]
The physical meaning of the perturbational term is explained in detail in [169],[169]. From these references it follows that solutions of Einstein’s field equations leading to the metric of the type given by Eq.(B.7) originate from the exact solution of the coupled Einsten-Maxwell fields. E.g.read pp 324-325 of [168] and ch.r 24 of [169]. Analysis of solutions for these fields made in [170]-[171] indicates that coupled Einstein-Maxwell fields can be of two types: a) non null and b) null. The null electromagnetic fields are such for which H and E We used null electromagnetic fields in [71,96] in connection with formation of torus-like knots while we explained their relevance to the Schrödinger fields in [14] and in section 7.5.3. Although the content of just cited references reproduces many results of Ward’s paper [167], we still have to add several important pieces of information. This is so, because we have not made yet a connection between metrics Eq.(B.7) and (B.4).
Toward this goal, we begin with the observation that the exact solution of the coupled Einstein-Maxwell fields was obtained by Rainich [] in 1925 and rediscovereed by Misner and Wheeler at the end of 1950ies [173]. They reformulated results by Rainich in the context of geometrodynamics. In short, this means that the exact solution of Einstein’s equations for the coupled gravity-Maxwell fields is obtainable in terms of the metric from which it is possible to restore both the gravitational and electromagnetic fields which are non null. The solution for the null fields was obtained much later [170],[171].
We are iterested in reinterpreting these null field results in terms of the Lie sphere geometry discussed in section 7 and Appendix D. To do so, initially we are following the paper by Kuiper [174]. From this paper it follows that the (pseudo) Riemannian space of dimensionality is flat if some region about any point if it can be covered by a metric -preferred coordinate system such that
[TABLE]
Above, , or for all points of
Definition B.7.* A (pseudo) Riemannian space is called conformally flat if some region about any point any point can be mapped conformally into a flat space. A conformally -flat space is obtainable from the flat space via change of metric,that is*
[TABLE]
where is the metric tensor defined in Eq.(B.8) while is some positive function of .
When is constant, then the above transformation is called ”similarity” or ”homothety”. When then the othety is an ”isometry”. Obviously, these definitions are not restricted to the diagonal metric tensor defined by Eq.(B.8).
**Theorem B.8.(Liouville theorem) **For spaces of dimensionality 2 under the conformal transformations described by Eq.(B.9) the (hyper) spheres are carried to (hyper) spheres.
Corollary B.9. From the same reference [174] and from section 7, Eq.(7.10), it follows that the conformal transformations are induced by the projective transformations leaving the Lie quadric Qn+1,2 invariant. From here the connection with the Lie sphere geometry follows. More on Liouville theorem can be found in [175].
Definition B.10. Killing vectors X are defined as solutions of the Killing equations
[TABLE]
where is the Lie derivative in the direction of X. These equations are describing the isometric-type of motion on Accordingly, the conformal Killing vectors X are defined as solutions of the conformal Killing equations
[TABLE]
These are describing the conformal-type of motion on . Examples of solutions of Eq.s(B.11) are demonstrated in [166].
**Definition B.11. The **conformal Lie group Ĉ( of a connected dimensional (pseudo) Riemannian manifold is conformally trivial if by using transformations defined by Eq.(B.9) it is possible to find a conformally-equivalent space in which Eq.(B.11) is replaced by Eq.(B.10).
**Definition B.12. **If use of Eq.(B.9) cannot bring Eq.(B.11) into Eq.(B.10), such conformal group is nontrivial.
In [166] the following remarkable theorem is proved
Theorem B.13.* (pseudo) Riemannian space has nontrivial conformal group only if it is conformally equivalent to the conformally flat space*
Remark B.14. Following [174] for conformally flat spaces about any point there is a region in which (for the parameter group of infinitesimal conformal transformations may exist. For we obtain 15 parameters group of conformal motions. This is SO(4,2) group describing conformal symmetries of hydrogen atom described in section 1. Based on Theorem B.8., this is the Lie sphere geometry group of motions.
Remark B.15. Following [176] the result provides maximal number of conformal generators for conformally flat spacetimes. However further studies demonstrated that Theorem B.13. can be extended. For *(pseudo) Riemannian spaces *which are not conformally flat the number of conformal generators is strictly less than the maximal number. It can be demonstrated [177] that for the metric of the type given by Eq.(B.7) the maximal number is 7.
**Remark B.16. Such **not conformally flat metrics are of no physical interest, however, for variety of reasons. First, we immediately loose the connection with the Lie sphere geometry in view of Theorem B.8. Second, once this connection is lost, we also loose the connection with the AdS-CFT correspondence (section 7 and Appendix D).
Remark B.17. In [166] it was proven that every Lorentzian spacetime with nontrivial conformal group is conformally equivalent to the Ricci flat space, e.g. read the Corollary B.2. In [1 ] it is argued (by invoking the content of Brinkmann’s theorem) that the only Ricci flat but nonflat 4-manifolds (that is 4-manifolds whose scalar curvature is not identically zero) admitting nonhomothetic conformal vector fields are manifolds described in terms of the pp waves metric given by Eq.(B.7). In the same reference it is being argued that for nonconformally flat spacetimes the dimension of the conformal group is at most 7.
Based on the last two remarks, it is physically meaningful to consider the conformally flat pp waves. In such a case we obtain the pp wave metric with circle-preserving 15 components conformal group which is Ricci flat. The condition of Ricci flatness is known to be ([169], page 386):
[TABLE]
The pp metric, Eq.(B.7), with determined from Eq.(B.12) is describing the coupled Einstein-Maxwell null fields. Eq.(24.43) of [169], page 384, explains why the condition, Eq.(B.12), is equivalent to the condition of Ricci flatness. The dependence of can be determined from the condition on the conformal Weil tensor to be zero. Such a condition selects spaces which are conformally flat. Examples of such selection are discussed in [166]. Based on just presented results, it should be clear that, although techically permissible, parts of the results of [152],[153] describing non null electromagnetic knots are without physical justification.
The seminal work of Penrose [], pages 271-275, entitled: ”Any space-time has a plane wave as a limit” explains the universal significance of the metric given by Eq.(B7) which is also known in physics literature as Penrose limit metrics. The significance of Theorem B.1. and appreciation of its role in relating the conformal transformations to Huygens’ principle apparently, totally escaped the attention of physics community. In Appendix D we provide solid evidence that the conformally flat space-times supporting the Lie sphere geometry are indeed the boundaries of space-times in which the AdS-CFT correspondence holds.
Appendix C
Dupin cyclides from Madelung’s hydrodynamical reformulation
of the Schrödinger equation
Using progressive waves method by Friedlander [70] we provided enough evidence in sections 4 and 7 that Dupin cyclides can exist as solutions of the Shrödinger equation. In this appendix we shall approach the same problem from yet another direction which is very illuminating in its own right. In our study we were motivated by two papers by Ferapontov [180],[181]. Ferapontov was interested in establishing the connection between the dynamics of Hamiltonian systems of hydrodynamical type and its realization in terms of differential geometry of hypersurfaces embedded in Minkowski-type spacetimes [180]. Using general results developed in [180], in [181] Ferapontov illustrated general results of [180] using hydrodynamical Hamiltonian systems which do not possess Riemann invariants343434Riemann invariants are to be discussed in a separate publication.. As result, he demonstrated the duality between such Hamiltonian systems and Dupin hypersurfaces. These are the hypersurfaces with constant principal curvatures [133]. In 3 dimensions mathematicians call them as ”Dupin surfaces”.
Almost immediately after Schrödinger published the installment of his foundational papers on quantum mechanics in Annalen der Physik in 1926 [1], on 25th of October of 1926 Madelung submitted his paper entitled ”Quantentheorie in hydrodynamischer form” to the Zeitschrift fur Physik [] where it was published in 1927. Although equations (3) and (4) of his paper contain typos, the key new equations (),() and () are correct. It is worth to reproduce these equations in this appendix. To avoid any ambiguities, we shall use exactly the same symbols as Madelung was using. He begins with the stationary Schrödinger equation
[TABLE]
in which is an energy. Next, he writes and uses this result in the time-dependent version of the Schrödinger equation
[TABLE]
Next, he was looking for a solution in the form where he is saying that, in view of Eq.(C.1), it is sufficient to consider only to be linearly-dependent upon while if we use Eq.(C.2), then it makes sense to consider both and to be time -dependent. Surely, the consistency between equations (C.1) and (C.2) requires to be linearly dependent upon Substituting into Eq.(C.2) and separating the real part from imaginary in resulting equation Madelung had obtained the following two equations353535We reproduce these equations without obvious typos which apparently were overlooked by Madelung when he was proofreading the galleys of his paper. His new key equations (),() and () are correct though.
[TABLE]
[TABLE]
By introducing new notation Eq.(C.4) is converted into the continuity equation (this is Madelung’s Eq.(
[TABLE]
Now, following Madelung, we introduce the velocity . In terms of such defined velocity Eq.(C.3) acquires the following form (this is Madelung’s Eq.(
[TABLE]
This equation can be rewritten in recognizable hydrodynamical form by taking into account that and by applying the gradient operator to Eq.(C.6a). In the end, we obtain: (Madelung’s Eq.()),
[TABLE]
Remark C.1. Eq.(C.6b) along with the continuity Eq.(C.5) describes the irrotational fluid moving under the action of conservative forces. In our book [55], page 67, we noticed (and described in various places of the book) the following chain of correspondences :
classical mechanics thermodynamicselectrodynamics
geometrical optics hydrodynamics magnetohydrodynamics
superconductivity non Abelian gauge Yang-Mills theories.
These correspondences are all derivable from the formalism of contact geometry and topology discussed in our book. Evidently, just described result by Madelung is fully consistent with the gauge-theoretic Floer-type description of Schrödinger’s quantum mechanics outlined in main text, subsection 4.2.3.
Now we recall that, according to Madelung, if we are only interested in the description of the stationary Schrödinger Eq.(C.1) the dependence of can be omitted. In such a case using Eq.(C.3) we obtain:
[TABLE]
and
[TABLE]
Finally, Eq.(C.4) now acquires the following form:
[TABLE]
A quick look at Eq.s (4.9 a-c) (producing Dupin’s cyclides solutions ) and comparing these with just obtained Eq.s( C.7)-(C.9) indicates that these two sets of equations are the same when But, as we know, in general .
Following ingenious ideas by Luneburg [33], which we are about to describe, just noticed difficulty can be resolved. By doing so, we are also going to rederive the results of section 6 via entirely different set of arguments.
We begin with Eq.(C.8) which we formally rewrite as
[TABLE]
Consider now a set of wavefronts: An orthogonal trajectory (a ray) through this wavefront at any point is normal to the wavefront through this point. In analogy with Eq.(4.13), we introduce the parameter along the ray trajectory so that this ray can be described in terms of the set {x(),y(),z()}. This allows us to introduce the equations for these rays as follows:
[TABLE]
Here is parameter describing various choices of parametrization of the ray trajectory. This freedom of reparametrization is the key element in achieving our goal-to bring Eq(C.10) into the form of Eq.(4.9a). Toward this goal, using Eq.(C.11) we consider the following chain of equalities
[TABLE]
Clearly,proceeding analogously, we obtain as well:
[TABLE]
and
[TABLE]
In equations (C.11) let then, in view of Eq.(C.10), we obtain
[TABLE]
That is with such choice of the parameter the yet arbitrary parameter becomes a parameter describing natural parametrization along the ray. We shall qualify it as ”time”. One still can do a better job, though, by noticing that, say,
[TABLE]
Here dEvidently, Eq.s (C.12) become the Newtonian equations of motion
[TABLE]
Here If this is so, then in view of Eq.s(C.10) and (C.11) we easily obtain363636Since now and and we can use Eq.s(C.10) and (C.11)..
[TABLE]
This time, however, we can bring Eq.s(C.13) and (C.18) in correspondence with each other by properly selecting ”time”. Introduce as well Using this definition and selecting we convert Eq.(C10) into
[TABLE]
easily recognizable as Eq.(4.9a) of the main text. Now we multiply Eq.(C.9) by to obtain instead
[TABLE]
Notice also that But etc. Since then In view of Eq.(C.20) and taking again into account that it only remains to demonstrate that Recall Eq.(C.11) and present as Use equations of motion, Eq.s(C.12), in order to write etc. Finally, write , etc. Using this result along with Eq.s(C.10),(C.11) we finally must prove that But and for the time -independent the desired result follows. Alternatively, by comparing Eq.s(C.13) and (C.18) and by noticing that Eq.(C.18) is converted to (C.13) for we obtain: for Thus, we just demonstrated that the Dupin cyclides generating set of equations, Eq.s (4.9), of the main text can be made to coincide with the the set of equations obtained by Madelung [177].
Appendix D
Various models of hyperbolic and Anti- de Sitter spaces
D.1. Hyperboloid model of hyperbolic space** **
Following Danciger [159], consider Rn,1 denoting Rn+1 equipped with Minkowski-type metric tensor :
[TABLE]
Using we define the hyperboloid of two sheets Sheets are determined by the sign of the first(or last) coordinate e.g. read [105]or [159] for more details. Traditionally, the hyperbolic space Hn is defined by the sheet for which Alternatively, it can be also defined as
[TABLE]
Taking a quotient identifies two sheets. The hyperboloid inherits the Riemannian metric of constant curvature from that defined by Eq.(D.1). Isometries of are defined by
[TABLE]
This result defines the Möbius-type transformation in accord with Definition 7.8. The orientation -preserving isometries are isometries lying in the identity component of : Isom For even (and this is our case) we have If we had chosen to think about as positive sheet of the hyperboloid then we should think about as the subgroup of that preserves the positive sheet.
D.2. The projective model of hyperbolic space and Mobius geometry** **
This model is easily obtainable from the hyperboloid model. Instead of Eq.(D.2) we have now
[TABLE]
where the symbol denotes equivalence. That is whenever there is some nonzero such that Thus, every hyperbolic structure is also a projective structure defined on some domain of RPn determined by the equivalence relation. This fact then allows us to define the boundary at infinity as
[TABLE]
We had encountered this equation already as Eq.(7.12). It can be interpreted either as the projectivized light cone or as the condition for two spheres to intersect transversally. has invariant flat conformal structure in the sense of Eq.(D.3). Since a geodesic in is determined by two distinct points on this means that in this formalism every geodesic in is in one-to-one correspondence with the orthogonal intersection of two spheres.
D.3. Hyperboloid model of Anti-de Sitter space** **
Let Rn+1,2 denote Rn+3 equipped with Minkowski-like metric tensor
[TABLE]
Using we define the hyperboloid Instead of Eq.(D2) defining the hyperbolic space we are having now the following definition of the Anti-de Sitter space (AdS)[]
[TABLE]
The hyperboloid inherits a Lorentzian metric of constant curvature from the form defined by Eq.(D.6). Isometries of are defined as
[TABLE]
D.4. The projective model of the Anti-de Sitter space and Lie sphere geometry** **
The projective model of AdS is easily obtainable from the hyperboloid AdS model. Instead of Eq.(D.7) we have now
[TABLE]
The boundary of AdS space is defined now by analogy with Eq.(D.5) that is
[TABLE]
In view of Remark 7.4. this result is easily recognizable as the Lie quadric Q It can be interpreted either as the projectivized light cone or as a condition for two spheres to touch each other. Accordingly, the geodesics in Anti-de Sitter space can be of three types:
a) A space-like geodesics is determined by two distinct points lying on
that is on Q In terms of the Lie sphere geometry such
a geodesics corresponds to two spheres touching each other.
b) A light-like geodesics has its both ends lying at the same point
of the quadric Q This situation is characterizes the projectivized
light cone. In the language of the Lie
sphere geometry such situation corresponds to spheres participating
in the formation of canal surfaces other than Dupin cyclides as explained
in section 7.5.1.
c) A time-like periodic geodesics which does not touch the boundary
The action of preserves in accord with Definition 7.6. and Corollary B.9. Furthermore, according to [161], page 184, Einstein’s space Ein Thus, Einstein’s space is the same as the the space where Lie sphere geometry acts and, using the discovered by Sophus Lie correspondence between the Lie sphere and Plücker line geometries described in section 7.4., it is possible to map such defined Einstein’s space into the space of twistors. The connection with Einstein spaces (and with spaces where Lie sphere geometry acts) and, therefore, with Penrose boundary can be easily explained now. For this, following [161] and using Eq.(7.10) (or Eq.(6.41a)) of the main text, we relate the quadratic form
[TABLE]
i to the nullcone as:
[TABLE]
Since each of is never zero, it is possible to divide both sides of Eq.(D.11b) by the positive number so that this equation can be rewritten as
[TABLE]
It describes the product For we obtain: This is an example of static Einstein spacetime with compactified time axis [22]. In fact, Êinn,1 while EinÊinKühnel and Rademacher had classified all (pseudo) Riemannian Einstein spacetimes having local or global conformal conformal groups [183]. These are surely including the pp wave type Ricci-flat spaces which, based on results of Appendix B, can be identified with Penrose limits.
References
[1] E.Schrödinger, Collected Papers on Wave Mechanics*, *
Chelsea Publ.Co, New York, 1978.
[2] R.Courant and D.Hilbert, Methods of Mathematical Physics,Vol.1,
Interscience Publishers, Inc., New York, 1953.
[3] B.Baker and E.Copson, The Mathematical Theory of Huygens’
Principle, Clarendon Press,Oxford, 1939.
[4] R.Feynman, Rev. Mod.Phys. 20 (1948) 367.
[5] M.Gutzwiller, Huygens’ principle and the path integral, in
Path Summation: Achievements and Goals (Trieste, 1987), pp. 47–73,
L.Schulman editor, World Sci. Publishing, Singapore, 1988.
[6] R.Feynman, R.Leighton and M.Sands, The Feynman Lectures on
Physics, Vol.3, Addison -Wesley Publishing Co., New York, 1989.
[7] R.Feynman and A.Hibbs, Quantum Mechanics and Path Integrals,
McGraw -Hill Co, New York, 1965.
[8] M.Born and E.Wolf, Principles of Optics, Cambridge U.Press.
Cambridge, UK, 1980
[9] C. Jönsson, Zeitschrift für Physik, 161(1961) 454.
[10] M.Arndt, O.Nairz, J.Vos-Andreae, C.Keller, G.Zouw and A.Zellinger,
Nature, 401 (1999) 680.
[11] S.Eibenberger, S. Gerlich, M. Arndt, M. Mayor and J. Tüxen, Phys.
Chem. Chem. Phys. 15 (2013) 14696.
[12] D.Bohm, Quantum Theory, Dover Publications Inc., New York, 1989.
[13] A. Sanz, and S. Miret-Artés, A Trajectory Description of Quantum
Processes I. Fundamentals, Springer-Verlag, Berlin, 2012.
[14] A.Kholodenko, arXiv:1703.04674.
[15] U. Leonhardt, Measuring the Quantum State of Light, Cambridge U.
Press, Cambridge, UK, 1997.
[16] M. Fox, Quantum Optics: An Introduction, Oxford U. Press,
Oxford,UK, 2006.
[17] U.Niederer, Helv. Phys. Acta 45 (1972) 802.
[18] R.Littlejohn, Phys. Rep. 138 (1986) 193.
[19] B.Wybourne, Classical Groups for Physicists, Wiley-Interscience Publ.,
New York, 1974.
[20] M.Kibler, Found.Chem. 9 (2007) 221.
[21] R. Campoamor-Stursberg, J.of Phys.: Conference Series
538 (2014) 012004.
[22] A.Keane and R.Barrett, Class.Quantum Grav. 17 (2000) 201.
[23] A. Bohm, Y. Ne’eman and A.Barut, Dynamical Groups and
Spectrum Generating Algebras, Vol’s 1 & 2, World Scientific,
Singapore, 1988.
[24] E.Schrödinger, Expanding Universes,
Cambridge University Press, Cambrige, 1956.
[25] S.Huggett and K.Tod, An Introduction to Twistor Theory
Cambridge University Press, Cambridge, 1994.
[26] E.Schrödinger, Science and the Human Temperament,
George Allen & Unwin LTD, London, 1935.
[27] V. Kac and P.Cheung, Quantum Calculus, Springer-Verlag, Berlin, 2002.
[28] E.Schrödinger, Phys.Rev.28 (1926) 1049.
[29] R.Courant and D.Hilbert, Methods of Mathematical Physics,Vol.2,
Interscience Publishers, Inc., New York, 1962.
[30] F.Cardin, Elementary Symplectic Topology and Mechanics,
Springer-Verlag, Heidelberg, 2015.
[31] D.Tannor, Introduction to Quantum Mechanics. A Time-Dependent
Perspective, University Science Books,Saustalito, Ca., 2007.
[32] M.Gosson and B.Hiley, Found.Phys.41 (2011) 1415.
[33] R.Luneburg, Mathematical Theory of Optics, U.of California Press,
Berkeley, 1966.
[34] A.Davydov, Quantum Mechanics, Pergamon Press, Oxford, 1976.
[35] L.De Broglie , An Introduction to the Study of Wave Mechanics,
Methuen & Co. LTD., 36 Essex Street, W.C. , 1930.
[36] C.Adams, M.Siegel and J. Mlynek, Phys. Reports 240 (1994) 143.
[37] R.Sawant, J.Samuel,A.Sinha,S.Sinha and U.Sinha,
PRL 113 (2014) 120406.
[38] K.Ito and H. McKean Jr., Diffusion Processes and Their Sample Paths,
Springer-Verlag,Berlin, 1965.
[39] H.De Raedt, K.Michelsen and K.Hess, Phys.Rev.A 85 (2012) 012101.
[40] U.Sinha,C.Couteau, T.Jennewein,R.Laflamme and G.Weihs,
Science 329 (2010) 418.
[41] I. Söllner, B.Gschösser, P. Mai, B. Pressl, Z. Vörös and G. Weihs,
Found. Phys. 42 (2012) 742.
[42] A. Sommerfeld, Optics, Academic Press, NY, 1964.
[43] L.Landau and E.Lifshitz, Classical Theory of Fields,
Butterworth-Heinemann, Oxford, 2000.
[44] A.Pazy, Semigroups of Linear Operators and Applications to
Partial Differential Equations, Springer-Verlag,Berlin, 1982.
[45] L.Evans, Partial Differential Equations,
AMS Publishers, Providence, RI, 2010.
[46] J.Hadamard, Lectures on Cauchy’s Problem in Linear Partial
Differential Equations, Nauka Co, Moscow, 1978
(Russian translation from French original).
[47] V.Arnol’d, Contact Geometry and Wave Propagation,
Enseign. Math. 36 (1990) 215.
[48] V. Arnol’d, Lectures on Partial Differential Equations,
Springer-Verlag, Berlin, 2004.
[49] V.Arnol’d, Mathematical Methods of Classical Mechanics,
Springer-Verlag, Heidelberg, 1978.
[50] V.Maslov and M. Fedoryuk, Semiclassical Approximation in Quantum
Mechanics. D. Reidel Publishing Co., Dordrecht-Boston, Mass., 1981.
[51] M.Fedoryuk, Partial Differential equations V,
Springer-Verlag, Berlin, 1999.
[52] V.Arnol’d , Ordinary Differential Equations,
Springer-Verlag,Berlin, 1992.
[53] C.Chicone, Oerdinary Differential equations with Applications,
Springer-Verlag, Berlin, 1999.
[54] A.Fasano and S.Marmi, Analytical Mechanics,
Oxford U.Press, Oxford, 2006.
[55] A.Kholodenko, Applications of Contact Geometry and Topology in Physics,
World Scientific,Singapore, 2013.
[56] M.Gotay, JMP 40 (1999) 2017.
[57] V.Guilleming and S.Sternberg, Variations on the Theme by Kepler,
AMS Publishers, Providence, RI 1990.
[58] M. Giaquinta and S. Hildebrandt,
Calculus of variations II, Springer-Verlag, Berlin, 1996.
[59] S.Benenti, Hamiltonian Structures and Generating Families,
Springer-Verlag, Berlin, 2011.
[60] I.Gelfand and S.Fomin, Calculus of Variations, Prentice-Hall Inc.,
New Jersey, 1963.
[61] M.Mathisson, Acta Math. 71 (1939) 249.
[62] L. Asgeirsson, Comm. Pure Appl. Math. 9 (1956) 307.
[63] P.Günther, Math.Intelligencer 13 (1991) 56.
[64] P.Günther, Huygens’ Principle and Hyperbolic Equations.
Academic Press, Inc., Boston, MA, 1988.
[65] F. Friedlander, The Wave Equations on a Curved Space-Time,
Cambridge U.Press, Cambridge, 1975.
[66] R.McLenaghan and J.Carminati, Ann.Inst.Henri Poincare 44 (1986) 115.
[67] R.Goldoni, J. Math. Phys. 18 (1977) 2125.
[68] N.Ibragimov, A.Oganesyan, Russian Math.Surveys 46 (1991) 137.
[69] M.Belger, R.Schimming and V.Wunsch, J.for Analysis and its
Applications 16 (1997) 9.
[70] F.Friedlander, Proc.Camb.Phil.Soc. 43(1946) 360.
[71] A.Kholodenko, Anal.Math.Phys. 6 (2016) 163.
[72] M.Audin and M.Damian, Morse Theory and Floer Homology,
Springer-Verlag, Berlin, 2014.
[73] D.Bohm and B.Hiley, Undivided Universe, Routledge,
New York and London, 1993.
[74] A.Sym, J.Nonlinear Math.Phys. 12, Supplement 1 (2005) 648.
[75] R.Pruss and A.Sym, Phys.Lett. A 336 (2005) 459.
[76] A.Sym and A.Szereszewski, SIGMA 7 (2011) 095.
[77] P.Broadbridge, C.Chanu and W.Miller Jr. SIGMA 8 (2012) 089.
[78] J. Alcázar, H. Dahl and G. Muntingh, arXiv:1611.06768.
[79] T.Cecil and P.Ryan, Geometry of Hypersurfaces, Springer-Verlag,
Berlin , 2015.
[80] Y.Suris, in Sophus Lie and Felix Klein. The Erlangen Program
and its Impact in Mathematics and Physics, EMS 23 (2015).
[81] J.Ratcliffe, Foundations of Hyperbolic Manifolds, Springer-Verlag,
Berlin, 1994.
[82] C.Codirla and H.Osborn, Ann.Phys. 260 (1997) 91.
[83] I.Bars and J.Terning, Extra Dimensions in Space and Time,
Springer Science +Business Media, LLC 2010.
[84] P.Arvidsson, R.Marnelius, arXiv:hep-th/0612060.
[85] S.Brodsky, G.de Teramond, H. Dosch and J.Elich,
Phys.Reports 584 (2015) 1.
[86] A.Kolesnik and N.Ratanov, Telegraph Processes and Option Pricing,
Springer-Verlag, Berlin, 2013.
[87] W.Greiner, Relativistic Quantum Mechanics.Wave Equations,
Springer-Verlag, Berlin, 2000.
[88] V. Berestetskii, E. Lifshitz and L.Pitaevskii, Relativistic Quantum
Theory, Pergamon Press, Oxford, UK, 1971.
[89] N.Bogoliubov and D.Shirkov, Introduction to the Theory of Quantized
Fields, John Wiley&Sons, New York,1976.
[90] M. Kac, Rocky Mountain J. of Math. 4 (1974) 497.
[91] B.Gaveau, T. Jacobson, M. Kac and L.Schulman, PRL 53 (1984) 419.
[92] C. Dewitt-Morette and P. Cartier, Functional Integration:
Action and Symmetries, Cambridge U.Press, Cambridge, 2006.
[93] A.Kolesnik and M. Pinsky, J.Stat.Phys.142 (2011) 828.
[94] E.Zauderer, Partal Differential Equations of Applied Mathematics,
Wiley-Interscience, New York, 1989.
[95] B.Gaveau and L.Schulman, Il Nuovo Cimento 11 (1989) 31.
[96] A. Kholodenko, Ann.Phys.201 (1990) 186.
[97] S. Deguchi and T.Suzuki, Phys.Lett.B 731 (2014) 337.
[98] S.Deguchi and S.Okano, Phys. Rev. D93 (2016) 045016.
[99] I. Bars, Phys.Rev.D 58 (1998) 066006.
[100] I.Araya and I.Bars, Phys.Rev.D 89 (2014) 066011.
[101] S. Singer, Lineraity,Symmetry and Prediction in the Hydrogen Atom,
Springer-Verlag, Berlin, 2005.
[102] V.Fock, Z.Phys. 98 (1935) 145.
[103] L.Landau and E. Lifshitz, Quantum Mechanics,
Elsevier Publ.Co, Amsterdam, 1981.
[104] H.Goldstein, Classical Mechanics,
Addison-Wesley Publ.Co., Reading, MA, 1980.
[105] I.Gelfand, R.Milnos and Z.Shapiro, Representations of the Rotation
and Lorentz Groups and Their Applications, Martino Fine Books,
P.O. Box 913, Eastford, CT, 2012. (Translation from the original
1958 Russian edition).
[106] J.Bros and G.Viano, Forum Math.8 (1996) 621.
[107] J.Bros and G.Viano, Forum Math.8 (1996) 659.
[108] J.Bros and G.Viano, Forum Math. 9 (1997) 165.
[109] D.Basu and S.Srinivasan, Czech.J.Phys.B 27 (1997) 635.
[110] J.Milnor, The American Mathematical Monthly 90 (1983), 353.
[111] T.de Laat, Regularization and quantization of the Kepler problem,
PhD Thesis, Department of Mathematics, Radboud University,
Nijmegen, Netheerlands, 2010.
[112] A.Fomenko, Symplectic Geometry, Gordon and Breach Co.,
New York, 1988.
[113] M.Dunajski, Solitons, Instantons and Twistors, Oxford U.Press,
Oxford,UK, 2010.
[114] R.Penrose and W. Rindler, Spinors and Space-Time,
Cambridge U.Press,Cambridge, UK, 1984.
[115] D.Sommerville, An Introduction to the Geometry of N Dimensions
Dower Publications , Inc., New York, 1958.
[116] H.Pottmann and J.Wallner, Computational Line Geometry,
Springer-Verlag, Berlin, 2010.
[117] R.Penrose, Relativistic symmetry groups, in
Group Theory and Nonlinear Problems, A.Barut Editor, pp.1-58,
D.Reidel Publ.Co., Boston, 1974.
[118] M.Kriele, Spacetime, Springer-Verlag, Berlin, 1999.
[119] A.Kholodenko, Int’l J.Mod Phys.A 30 (2015) 1550189.
[120] A.Kholodenko, E.Ballard, Physica A 380 (2007) 115.
[121] M.Bander and C.Itzykson, Rev.Mod.Phys. 38 (1966) 346.
[122] I.Frenkel and M.Libine, Adv.Math. 218 (2008) 1806.
[123] A.Hurwitz and R.Courant, Complex Function Theory (in German)
Interscience Publishers, Inc., New York, 1944.
[124] R. Fueter, Comment. Math. Helv. 8 (1) (1935) 371.
[125] B.Cordani, The Kepler Problem, Birkhäuser, Basel, 2003.
[126] I.Bars, Phys.Rev. D 58 (1998) 066006.
[127] I.Bars, C. Deliduman and O.Andreev, Phys. Rev. D 58 (1998) 066004.
[128] P.Kustaanheimo and E.Stiefel, J.Fur.Reine
und Angevante Math. 218 (1965), 204.
[129] F.Cornish, J.Phys.A 17 (1984) 323.
[130] A.Chen, Phys.Rev.A 22 (1980) 333.
[131] A.Chen, Phys.Rev.A 23 (1981) 1655.
[132] A.Chen,Phys.Rev.A 26 (1982) 669.
[133] G.Thorbergsson, Bull. London Math.Soc. 15 (1983) 493.
[134] E.Huhnen-Venedey, PhD Thesis, department of Mathematics,
Technical University,Berlin, 2007.
[135] E.Musso, G.Jensen and L.Nicolodi, Surfaces in Cassical Geometries,
Springer -Verlag, Berlin, 2016.
[136] C.Doran and A.Lasenby, Geometric Algebra,
Cambridge U.Press, Cambridge, 2004.
[137] A.Bobenko and Y.Suris, Discrete Differential Geometry,
AMS Publishers, Providence, RI, 2008.
[138] A.Bobenko and E.Huhnen-Venedey, Geom. Dedicata 159 (2012) 207.
[139] L.Drouton, L.Fuchs, L.Garnier and R.Langevin,
Adv. in Appl.Clifford Algebra 24 (2014) 515.
[140] L.Dorst, Math.Comput. Sci 10 (2016) 97.
[141] S.Helgason, in The Sophus Lie memorial conference (Oslo 1992),
pages 3-21, Scand.U.Press, Oslo, 1994.
[142] L.Mason and N.Woodhouse, Integrability, Self-Duality,
and Twistor Theory, Clarendon Press, Oxford, 1996.
[143] R.Ward, R.Wells,Jr. , Twistor Geometry and Field Theory,
Cambridge University Press, Cambridge, 1990.
[144] F.Klein, Vorlesungen Über Höhere Geometrie,
Spinger-Verlag, Berlin, 1926.
[145] W. Blaschke, Volresungen Über Differential-Geometrie III,
Springer-Verlag , Berlin, 1929.
[146] M.Schrott and B.Odenhal, J. for Geometry and Graphics, 10 (2006)73.
[147] A.Gray, Modern Differential Geometry of Curves and Surfaces
with Mathematica, Taylor&Francis, Roca Baton, FL 2006.
[148] T. Ivey, AMS Proceedings 123 (1995) 865.
[149] P. Colapino, Articulating Space: Geometric Algebra for Parametric
Design-Symmetry Kinematics and Curvature, PhD Thesis,
U.of California, Santa Barbara, 2016.
[150] A.Bobenko and E.Huhnen -Venedey, Geom.Dedicata 159 (2012) 207.
[151] M.Lavicka and J.Vrsek, J.for Geometry and Graphics 13 (2009)145.
[152] M.Arrayas, D. Bouwmeester and J. Trueba, Phys.Reports 667 (2017), 1.
[153] C.Hoyos, N.Sicar and J.Sonnenschein, J.Phys.A 48 (2015) 255204.
[154] H.Nastase, Introduction to the AdS-CFT Correspondence, Cambridge
U.Press, Cambridge, UK, 2015.
[155] A.Kholodenko, J.Geom.Phys.35 (2000),193.
[156] A.Kholodenko, J.Geom.Phys. 38 (2001) 81.
[157] A.Kholodenko, J.Geom.Phys.43 (2005) 45.
[158] C.Frances, Comm.Math.Helv. 80 (2005) 883.
[159] J.Danciger, Geometric Transitions: From Hyperbolic to AdS Geometry,
PhD Thesis, Department of Mathematics, Stanford University, 2011.
[160] S.Takeuchi, Japanese J.of Appl.Physics 53 (2014) 030101.
[161] T.Barbot, V.Charette, T.Drumm, W.Goldman and K.Melnick, in
Recent Developments in Pseudo-Riemannian Geometry,
D.Alekseevsky and H.Baum Editors, pp 179-229,
ESI Lect. Math. Phys., Eur. Math. Soc., Zürich, 2008.
[162] V. Mazya and T.Shaposhnikova, Jaques Hadamard,
A Universal Mathematician,
AMS Publishers, Providence, RI, 1998.
[163] Y.Kravtsov, Geometric Optics in Engineering Physics,
Alpha Sci.International Ltd., Harrow, UK. 2005.
[164] A.Petrov, Einstein Spaces, Pergamon Press, New York, 1969.
[165] A.Petrov, New Methods in General Theory of Relativity,
Nauka, Moscow, 1966 (in Russian)
[166] N.Ibragimov, Transformation Groups Applied to Mathematical Physics,
D.Reidel Publ.Co., Boston, MA, 1985.
[167] R.Ward, Class.Quant.Gravity 4 (1987) 775.
[168] J.Griffits and J.Podolsky, Exact Space-Times in Einstein’s
General Relativity, Cambidge U.Press, Cambridge UK, 2009.
[169] H.Stephani, D.Kramer,M.MacCallum, C.Hoenselaers and E.Herlt,
Exact Solutions of Einstein’s Field Equations,
Cambridge U.Press, Cambridge,UK, 2003.
[170] R.Geroch, Ann.Phys. 36 (1966) 147.
[171] C.Torre, Class.Quantum Grav. 31 (2014) 045022.
[172] G.Rainich, AMS Transactions 27 (1925) 106.
[173] Ch.Misner and J.Wheeler, Ann.Phys. 2 (1957) 525.
[174] N.Kuiper, Ann.Math. 30 (1949) 916.
[175] W.Kühnel and H-B. Rademacher, J.Math. Pures Appl. 88 (2007) 251.
[176] A.Aminova, Russ.Math.Surveys 50 (1995) 69.
[177] A.Keane and B.Tupper, Class.Quantum Grav.21 (2004) 2037.
[178] W.Kühnel and H-B. Rademacher,Geom.Dedicata 109 (2004) 175.
[179] M.Cahen and M.Flato, Differential Geometry and Relativity,
D.Reidel Publ.Co., Boston, MA.1976.
[180] E.Ferapontov, Soviet Jorn.Math. 55 (1991) 1970.
[181] E.Ferapontov, Differential Geom. and Applications 5 (1995) 121.
[182] E.Madelung, Zeit.f.Phys.40 (1927) 322.
[183] W. Kühnel and H-B. Rademacher, Result.Math. 56 (2009) 421.
