text:
# New correlation duality relations for the planar Potts model

## 1 Introduction

The Potts model, a generalization of the Isling model to $`q`$ spin components introduced by Potts  in 1952, has been at the forefront of research interest for almost five decades. While the critical properties of planar models are now largely known  \- , much less is known about their correlation functions.

A key element of planar spin models is the notion of duality relations. Every planar graph, or lattice, $`G`$ has an associated dual planar graph $`G^{(\mathrm{D})}`$, and it is well-known  that the partition functions of the Potts model on $`G`$ and $`G^{(\mathrm{D})}`$ are proportional to each other when their couplings are correctly related. In the case of $`q=2`$, the Ising model, it is also known  that a duality relation exists between certain two-spin correlation functions. It has long been suspected that the correlations of the Potts models on $`G`$ and $`G^{(\mathrm{D})}`$ should be similarly related. Indeed, in a series of recent papers we have obtained correlation dualities for vertices (spins) residing on the boundary of one face  \-  or two faces  of $`G`$, and reformulated the 1–face results in graph-theoretical terms . In this paper we consider the $`N`$-face problem, namely, correlation duality relations for vertices residing on $`N`$ faces of $`G`$ for general $`N`$. As we shall see, a new picture emerges for $`N2`$. Instead of duality relations between individual correlation functions, we now obtain relations between linear combinations of correlation functions.

The organization of this paper is as follows. In Sec. 2 we establish notations and briefly review the Fortuin-Kasteleyn (F-K) representation of the Potts model, and derive the duality relation for the partition function. Using two identities which we establish as lemmas, we introduce in Section 3 a method which allows us to generate correlation duality relations. While our method does not necessarily generate all such duality relations, we show that generally the duality relates linear combinations of, rather than individual, correlation functions. Explicit examples are given in Section 4 and the relation to earlier results is explained. In Section 5 we reformulate our results in graph-theoretical terms as rooted dichromatic, or Tutte, polynomials.

## 2 The setup

### 2.1 The Fortuin-Kasteleyn representation

Consider a finite, planar connected graph $`G`$ with vertex set $`V`$, edge set $`E`$, and face set $`F`$. We write $`|V|`$ for the size of the set $`V`$ and similarly for $`E`$, $`F`$, and other terms. Let $`\sigma `$ denote an assignment of colors, or states, to the vertices of $`G`$, which is a map from $`V`$ to the color set $`𝒞=\{1,\mathrm{},q\}`$. Then the partition function of the Potts model on $`G`$ is

$$Z=\underset{\sigma }{}\underset{i,j}{}e^{K\delta (\sigma _i,\sigma _j)}$$
(1)
where the summation is over the $`|V|^q`$ coloring assignments $`\sigma `$, the product runs over all nearest neighbor spin pairs interacting with the coupling parameter $`K`$, $`\delta `$ is the Kronecker delta function, and $`\sigma _i`$ is the color of the $`i`$-th vertex.

The Fortuin-Kasteleyn (F-K) representation  of the Potts model replaces the configurational sum in (1) by a sum over edge sets. Let $`SE`$ be an edge set on the graph $`G`$. If we remove from $`G`$ all the edges which are not in $`S`$, we obtain a subgraph which contains all the vertices, and which in general has many connected components. This is called a spanning subgraph of $`G`$. We write $`|S|`$ to denote the number of edges in $`S`$, and write $`p(S)`$ to denote the number of connected components in the spanning subgraph defined by $`S`$. In this picture, any isolated vertex in the spanning subgraph counts as a connected component, and so contributes a summand 1 in the counting of $`p(S)`$. Define

$$v=e^K1$$
(2)
so that we have

$$e^{K\delta (\sigma _i,\sigma _j)}=1+v\delta (\sigma _i,\sigma _j)$$
(3)
The substitution of (3) into (1) now gives rise to the F-K representation of the partition function

$$Z(G;q,v)Z=\underset{SE}{}v^{|S|}q^{p(S)}$$
(4)
The number of colors $`q`$ appears as a variable in (4), rather than a summation limit as in (1). For this reason the F-K representation (4) forms the basis of essentially all analyses of the Potts model.

### 2.2 Dual partition functions

As a precursor to discussions in later sections, we give here a graph-theoretical derivation of the well-known duality relation for the partition function.

The standard dual graph of $`G`$ is constructed by placing vertices in the faces of $`G`$, and connecting them with edges which are transverse to the edges of $`G`$. A pair of transverse edges are said to be dual to each other. The resulting dual graph $`G^{(\mathrm{D})}`$ has $`|F|`$ vertices, $`|E|`$ edges and $`|V|`$ faces, and is also planar. For a given edge set $`S`$ on $`G`$, we define its dual edge set $`S^{(\mathrm{D})}`$ on $`G^{(\mathrm{D})}`$ by requiring that an edge of $`G^{(\mathrm{D})}`$ is in $`S^{(\mathrm{D})}`$ if and only if its dual edge is not in $`S`$. Thus,

$$|S|+|S^{(\mathrm{D})}|=|E|$$
(5)
Moreover, if we start from $`|V|`$ isolated vertices of $`G`$, and add one by one the edges of $`S`$, then each edge reduces $`p(S)`$ by 1 unless it forms an independent cycle, in which case $`p(S)`$ is unchanged. Therefore we have the relation

$`p(S)`$ $`=`$ $`|V||S|+c(S)`$ (6)
$`=`$ $`|V||S|+p(S^{(\mathrm{D})})1`$
where $`c(S)`$ is the number of independent cycles in the spanning subgraph of $`S`$, and is related to $`p(S^{(\mathrm{D})})`$, the number of components in $`S^{(\mathrm{D})}`$, through the topological relation $`c(S)=p(S^{(\mathrm{D})})1`$. Using (5) and (6) we have the identity

$`v^{|S|}q^{p(S)}`$ $`=`$ $`v^{|E||S^{(\mathrm{D})}|}q^{|V||S|+p(S^{(\mathrm{D})})1}`$ (7)
$`=`$ $`k(q/v)^{|S^{(\mathrm{D})}|}q^{p(S^{(\mathrm{D})})}`$
where

$$k=q^{1|F|}v^{|E|}$$
(8)
and we have used the Euler relation

$$|V|+|F|=|E|+2$$
(9)
Substituting (7) into (4), we see that the Potts model partition functions on $`G`$ and $`G^{(\mathrm{D})}`$ are proportional, when the couplings $`K`$ and $`K^{}`$ are related by the equation

$$vv^{}=q$$
(10)
where $`v^{}=e^K^{}1`$. The partition function duality is the statement that

$$Z(G;q,v)=kZ(G^{(\mathrm{D})};q,v^{})$$
(11)
where we have made use of the one-to-one correspondence between edge sets $`S`$ and dual edge sets $`S^{(\mathrm{D})}`$. Note that we have

$$kk^{(\mathrm{D})}=q^{1|F|}v^{|E|}q^{1|V|}(v^{})^{|E|}=1$$
(12)
which says that the dual of the dual partition function is the partition function itself.

### 2.3 The correlation function

The Potts correlation function is defined as the probability that a given set of vertices are assigned fixed colors . Specifically, let $`RV`$ be any subset of the vertices which are assigned fixed colors. We shall call these the roots for the correlation functions, and the graph with roots a rooted graph. Let the $`|R|`$ roots lie on the boundaries of $`N`$ distinct faces of $`G`$, which will be called the external faces of the graph. For simplicity, we shall assume that external faces are not adjacent to each other, namely they do not share any vertices or boundary edges.

Let $`c:R𝒞`$ be any assignment of colors on the roots, so that $`c_i`$ is the color assigned to the root $`i`$, and we write

$$\delta _c(\sigma )=\underset{iR}{}\delta (\sigma _i,c_i)$$
(13)
Then the correlation function for the assignment $`c`$ on the root set $`R`$ is defined to be the expectation value of $`\delta _c(\sigma )`$, namely,

$$P_c\delta _c(\sigma )=Z(G;q,v)^1\underset{\sigma }{}\delta _c(\sigma )\underset{i,j}{}e^{K\delta (\sigma _i,\sigma _j)}$$
(14)
where the summation in the right-hand side is a partial partition function. Clearly, the explicit expression of $`P_c`$ depends on the location of the roots as well as the explicit color assignment $`c`$.

Every color assignment $`c`$ defines a partition $`X`$ of the roots, namely a division of the $`|R|`$ roots into disjoint blocks, by the rule that two vertices $`i,j`$ in $`R`$ belong to the same block of $`X`$ if and only if $`c_i=c_j`$. By the symmetry of the Potts interactions, the correlation (14) depends only on this partition defined by $`c`$, not on the specific colors. We can therefore associate each partial partition function to a partition $`X`$. Alternately, for any partition $`X`$ of the roots, we have a partial partition function

$$Z_XZ_X(G;q,v)=\underset{\sigma }{}\delta _c(\sigma )\underset{i,j}{}e^{K\delta (\sigma _i,\sigma _j)}$$
(15)
where $`c`$ is any color assignment that produces the partition $`X`$. Since $`P_c`$ and $`Z_X`$ are proportional, it is sufficient to discuss duality relations for the partial partition function $`Z_X`$.

The F-K representation also extends to the partial partition functions. For any edge set $`S`$, we define a “connected” partition $`\pi (S)`$ of $`R`$ by the rule that two roots $`i,j`$ belong to the same block if and only if they belong to the same connected component of the spanning subgraph defined by $`S`$. Using this notion, we can separate the edge sets into classes

$$𝒮(X)=\{SE:\pi (S)=X\}$$
(16)
labeled by their corresponding connected partitions of $`R`$. Clearly every edge set $`S`$ belongs to exactly one of these classes $`𝒮(X)`$.

Now the role of the factor $`\delta _c(\sigma )`$ in (14) and (15) is to restrict the summations to fixed colors for vertices in $`R`$. The net result is that each cluster containing roots contributes a factor 1, instead of the factor $`q`$ in (4). This leads us to define, for any partition $`X`$, the summation over edge sets $`S`$ with $`\pi (S)=X`$,

$$F_X(G;q,v)=\underset{S𝒮(X)}{}v^{|S|}q^{p(S)|X|},$$
(17)
where $`|X|`$ is the number of blocks in $`X`$. The set of partitions of $`R`$ is partially ordered. We write $`YX`$ to mean that every block in $`Y`$ is wholly contained in a block in $`X`$ or, equivalently, $`Y`$ is a refinement of $`X`$. Then we have the identity

$$Z_X(G;q,v)=\underset{YX}{}F_Y(G;q,v)$$
(18)
### 2.4 Möbius inversion and planar identities

It is well-known  that the partially ordered sum (18) can be inverted, allowing us to write $`F_X`$ as a linear combination of $`Z_Y`$. For $`YX`$ define the Möbius inversion coefficient

$$\mu (Y,X)=(1)^{|Y||X|}\underset{bX}{}(n_b(Y)1)!$$
(19)
where the product on the right side runs over blocks $`b`$ in $`X`$, and $`n_b(Y)`$ is the number of blocks of $`Y`$ that are contained in $`b`$. Then we have for every partition $`X`$,

$$F_X(G;q,v)=\underset{YX}{}\mu (Y,X)Z_Y(G;q,v)$$
(20)
One peculiarity for planar graphs is that $`F_X`$ can be zero for certain partitions $`X`$. This phenomenon was explored in detail in  in the case of a single external face. In this case, one has $`F_X=0`$ if the partition $`X`$ is non-planar (also called a crossing partition in the mathematics literature ). It then follows from (20) that whenever $`X`$ is a non-planar partition, the partial partition functions $`Z_Y`$ must satisfy the identity

$$\underset{YX}{}\mu (Y,X)Z_Y(G;q,v)=0$$
(21)
giving rise to “sum-rule” relations for the partial partition functions . In general, when there are several external faces, identities of this type also hold for any root set $`R`$, although it seems to be a hard problem to decide if a given partition is planar or not.

## 3 Correlation duality relation

### 3.1 Dual rooted graph $`G^{}`$

In Section 2.2 we introduced the dual graph $`G^{(\mathrm{D})}`$, and described the well-known duality for partition functions on $`G`$ and $`G^{(\mathrm{D})}`$. In this section we introduce the dual rooted graph $`G^{}`$ which plays the corresponding role for correlation functions.

Starting from a rooted graph with $`|R|`$ roots on $`N`$ external faces, one can construct a dual rooted graph $`G^{}`$. The procedure of constructing $`G^{}`$ for $`N=1`$ has been described in detail in . The construction of $`G^{}`$ for general $`N`$ repeats the $`N=1`$ process for each external face, a prescription we now briefly describe.

Place $`N`$ extra vertices, one each in the center of each of the $`N`$ external faces. For each external face, let $`f`$ be the extra vertex and connect it to each root on the boundary of this face by an edge. This gives a new graph $`G^{\prime \prime }`$ which has $`N`$ more vertices than $`G`$ and $`|R|`$ additional edges. The dual graph of $`G^{\prime \prime }`$ is also planar, and it has $`N`$ faces containing the $`N`$ extra vertices $`f`$. Now remove all edges of these $`N`$ faces, and the resulting graph is the dual rooted graph $`G^{}`$.

The vertices of $`G^{}`$ residing inside the $`N`$ external faces of $`G`$ are now the dual roots and we denote this set by $`R^{}`$. Clearly, the number of edges of $`G^{}`$ is $`|E|`$, and we have $`|R^{}|=|R|`$ (because external faces are non-adjacent). It is also clear that the total number of vertices of $`G^{}`$ is

$$|V^{}|=|F|+|R|N$$
(22)
We can recover the standard dual graph $`G^{(\mathrm{D})}`$ from $`G^{}`$ by fusing the dual roots $`R^{}`$ inside each external face into a single dual vertex. This gives a one-to-one correspondence between edges on $`G^{(\mathrm{D})}`$ and $`G^{}`$, and hence we define the dual edge set $`S^{}`$ on $`G^{}`$ to be the same edge set as $`S^{(\mathrm{D})}`$ on $`G^{(\mathrm{D})}`$. Therefore in particular we have

$$|S^{}|=|S^{(\mathrm{D})}|$$
(23)
Let $`G^{(\mathrm{F})}`$ denote the graph generated by fusing all roots $`R`$ in $`G`$, and $`(G^{})^{(\mathrm{F})}`$ the graph generated by fusing all roots $`R^{}`$ in $`G^{}`$. Then, we have the identities

$$(G^{})^{(\mathrm{F})}=G^{(\mathrm{D})}\mathrm{and}G^{(\mathrm{F})}=(G^{})^{(\mathrm{D})}$$
(24)
### 3.2 A preliminary relation

Our results will depend crucially on the topology of the connectedness of the partitions $`\pi (S)`$ and $`\pi (S^{})`$. Fix an edge set $`S`$ on $`G`$, and let $`X=\pi (S)`$, $`Y=\pi (S^{})`$. Recall that $`|X|`$ is the number of blocks in $`X`$. Also, let $`n(X)`$ be the number of clusters of external faces on $`G`$ connected by $`S`$ (where we define two external faces to be in the same cluster if a connected component of $`S`$ contains roots on both faces). Similar definitions apply to $`|Y|`$ and $`n(Y)`$.

Now $`S^{}`$ and $`S^{(\mathrm{D})}`$ are the same edge set, and $`G^{(\mathrm{D})}`$ is obtained from $`G^{}`$ by fusing. Hence we have

$`p(S^{})`$ $`=`$ $`|Y|+\mathrm{\Delta }`$
$`p(S^{(\mathrm{D})})`$ $`=`$ $`n(Y)+\mathrm{\Delta }`$ (25)
where $`\mathrm{\Delta }`$ is the number of connected components in the spanning graph generated by $`S^{}`$ that are not connected to any dual roots. It follows that we have

$$|Y|n(Y)=p(S^{})p(S^{(\mathrm{D})})$$
(26)
and similarly

$$|X|n(X)=p(S)p(S_{}^{}{}_{}{}^{(\mathrm{D})})$$
(27)
where $`S_{}^{}{}_{}{}^{(\mathrm{D})}`$ on $`G_{}^{}{}_{}{}^{(\mathrm{D})}`$ is the edge set dual to $`S^{}`$ on $`G^{}`$. Combining (26) and (27) and making use of (6), we obtain the identity

$`|X|+|Y|n(X)n(Y)`$ $`=`$ $`[p(S)p(S^{(\mathrm{D})})]+[p(S^{})p(S_{}^{}{}_{}{}^{(\mathrm{D})})]`$ (28)
$`=`$ $`[|V||S|1]+[|V^{}||S^{}|1]`$
$`=`$ $`|V^{}||F|`$
$`=`$ $`|R|N`$
where use has also been made of (22), the Euler relation (9), and the relation

$$|S|+|S^{}|=|E|$$
(29)
For $`N=1`$, we have $`n(X)=n(Y)=1`$, and (28) reduces to

$$|X|+|Y|=|R|+1$$
(30)
a relation used in . For $`N=2`$, (28) reduces to Eq. (21) of Ref. .111After the replacements of $`n(X)2\chi (S),n(Y)2\chi (S^{})`$, and $`|R|NN2`$.

### 3.3 First identity

Our construction of correlation duality relations is based on the use of two identities. The first identity, which we now state as a lemma, is a generalization of the identity (7). The new ingredient here is that the generalization should permit the consideration of the factor $`v^{|S|}q^{p(S)|X|}`$ appearing in the F-K representation (17).

###### Lemma 1

For a fixed edge set $`S`$ on $`G`$, define $`X=\pi (S)`$ and $`Y=\pi (S^{})`$. Then we have the identity

$$v^{|S|}q^{p(S)[|X|n(X)]/2}=k^{}(v^1q)^{|S^{}|}q^{p(S^{})[|Y|n(Y)]/2}$$
(31)
where

$$k^{}=v^{|E|}q^{1|F|[|R|N]/2}$$
(32)
Remark: Notice that $`v^1q=e^K^{}1`$, so we define $`v^{}=v^1q`$ as before, and we can rewrite (31) as follows to indicate the explicit dependence on $`S`$ and $`S^{}`$,

$$v^{|S|}q^{p(S)[|\pi (S)|n(\pi (S))]/2}=k^{}(v^{})^{|S^{}|}q^{p(S^{})[|\pi (S^{})|n(\pi (S^{}))]/2}$$
(33)
Proof: The proof of the lemma parallels that of the partition function duality given in Section 2.2. First, the relation (6) still holds so that after using (29) the first line of (7) becomes

$$v^{|S|}q^{p(S)}=v^{|E||S^{}|}q^{|V||S|+p(S^{(\mathrm{D})})1}$$
(34)
The lemma (31) now follows after the use of the identities (26) and (28). Notice that we have

$$k^{}(k^{})^{}=v^{|E|}q^{1|F|[|R|N]/2}(v^{})^{|E|}q^{1|V^{}|[|R|N]/2}=1$$
(35)
where $`|V^{}|`$, the number of faces of the dual of $`G^{}`$, is given by (22).

### 3.4 Second identity

Let $`S`$ and $`T`$ be edge sets which define the same partition of $`R`$, i.e., $`\pi (S)=\pi (T)`$, then a key feature of the case $`N=1`$ is that their dual edge sets also define the same partition of $`R^{}`$. Namely,

$$\pi (S)=\pi (T)\mathrm{if},\mathrm{and}\mathrm{only}\mathrm{if}\pi (S^{})=\pi (T^{}),N=1$$
(36)
This feature permits us to derive correlation duality relations for $`N=1`$  \- . But (36) is false for $`N2`$ . However a weaker statement holds which we now state as a lemma.

First, some definitions: A block of $`X`$ is local if it contains roots from one face only, and is non-local otherwise. A partition $`X`$ is local if every block of $`X`$ is local, and is non-local if any block in $`X`$ is non-local. Given a partition $`X`$, we can construct from it a unique local partition $`X_{\mathrm{loc}}`$ as follows: We split up every non-local block into a collection of local blocks, by separating its roots which lie on different faces. The collection of all these local blocks is the new partition $`X_{\mathrm{loc}}`$. See Fig. 1 for an example.

Any local partition is a collection of planar partitions of the roots on each of the external faces. Each of these planar partitions on a single face has a unique dual partition. So for any local partition $`\xi =X_{\mathrm{loc}}`$, we will define $`\xi ^{}`$ to be the collection of these dual partitions on $`G^{}`$, and by extension we will call it the dual partition of $`\xi `$. Notice that $`\xi ^{}`$ is always a local partition on $`G^{}`$.

###### Lemma 2

For any edge set $`S`$ on $`G`$, we have

$$\pi (S^{})_{\mathrm{loc}}\pi (S)_{\mathrm{loc}}^{}$$
(37)
Furthermore, let $`S`$ and $`T`$ be two edge sets on $`G`$, and suppose that $`\pi (S)=\pi (T)`$. Then we have

$$\pi (S^{})_{\mathrm{loc}}=\pi (T^{})_{\mathrm{loc}}$$
(38)
Remark: The result corresponding to (38) also holds on $`G^{}`$, namely, if $`\pi (S^{})=\pi (T^{})`$ then also $`\pi (S)_{\mathrm{loc}}=\pi (T)_{\mathrm{loc}}`$.

Proof: It is sufficient to establish the relations (37) and (38) for each external face separately. So pick one external face, call it $`f`$, and let $`R_f`$ denote the roots which belong to this face, and $`R_f^{}`$ their dual roots. We will construct two different $`N=1`$ rooted graphs from $`(G,R)`$, as follows (recall that a rooted graph is an ordered pair, consisting of a graph together with a set of roots. Mostly we have ignored the distinction between graph and rooted graph, but we will distinguish them throughout this proof).

First, by ignoring the roots on the other faces we obtain the rooted graph $`(G,R_f)`$. As a shorthand we will denote this pair by $`G_f`$. Second, by fusing the roots on the other external faces (but not on $`f`$) we obtain a different rooted graph $`(G_f^{(F)},R_f)`$, which also has roots only on the face $`f`$. (Recall that we fuse roots on a face by merging them into a single vertex). Again as a shorthand we denote this by $`G_f^F(G_f^{(F)},R_f)`$.

Each of these rooted graphs has a dual rooted graph, which we denote by $`G_f^{}`$ and $`(G_f^F)^{}`$ respectively. These are constructed according to the procedure described in Section 3.1 for the $`N=1`$ case, using only the roots on the face $`f`$. Notice that the underlying graph of $`(G_f^F)^{}`$ is identical to the graph $`G^{}`$, because the process of fusing roots on a face in $`G`$ automatically splits the dual vertex in $`G^{(\mathrm{D})}`$ to produce the dual roots. As a convenient shorthand we write $`(G^{})_f`$ to denote the pair $`(G^{},R_f^{})`$, which is obtained from $`G^{}`$ by ignoring the roots on the other faces. It follows then that

$$(G_f^F)^{}=(G^{})_f$$
(39)
Now let $`S`$ be an edge set on $`G`$. Since the underlying graphs of $`G_f`$ and $`G_f^F`$ have the same edges as $`G`$, $`S`$ also defines partitions of the roots $`R_f`$ on $`G_f`$ and $`G_f^F`$. For clarity in this proof, we write these partitions as $`\pi (S;G_f)`$ and $`\pi (S;G_f^F)`$ respectively, where we include the graph in order to distinguish them. Note that $`\pi (S;G_f)`$ is one part of the local partition $`\pi (S)_{\mathrm{loc}}`$, namely the part that contains the roots on $`f`$. Let $`\pi (S)_f`$ denote this part of the partition $`\pi (S)_{\mathrm{loc}}`$, and similarly $`\pi (S^{})_f`$ the part of $`\pi (S^{})_{\mathrm{loc}}`$ which contains the dual roots on $`f`$. So we have

$$\pi (S;G_f)=\pi (S)_f,\pi (S^{};(G^{})_f)=\pi (S^{})_f$$
(40)
To prove the first result (37), note that since $`G_f^F`$ is obtained from $`G_f`$ by fusing vertices we have

$$\pi (S;G_f)\pi (S;G_f^F)$$
(41)
These are both $`N=1`$ partitions, hence they have dual partitions, and these satisfy

$$\pi (S;G_f)^{}\pi (S;G_f^F)^{}=\pi (S^{};(G_f^F)^{})$$
(42)
Using (39) and (40) this gives

$$\pi (S)_f^{}=\pi (S;G_f)^{}\pi (S^{};(G^{})_f)$$
(43)
Combining these for every external face gives

$$\pi (S)_{\mathrm{loc}}^{}\pi (S^{})_{\mathrm{loc}}$$
(44)
To prove the second result (38), note that since $`\pi (S)=\pi (T)`$ it follows also that

$$\pi (S;G_f^F)=\pi (T;G_f^F).$$
(45)
Taking the dual of both sides gives

$$\pi (S^{};(G_f^F)^{})=\pi (T^{};(G_f^F)^{})$$
(46)
and hence

$$\pi (S^{};(G^{})_f)=\pi (T^{};(G^{})_f)$$
(47)
Again combining these for every external face gives

$$\pi (S^{})_{\mathrm{loc}}=\pi (T^{})_{\mathrm{loc}}$$
(48)
Q.E.D.

### 3.5 Form of duality relations

As a convenient shorthand, for partitions $`X`$ of $`R`$ and $`Y`$ of $`R^{}`$ we write

$$Z_X=Z_X(G;q,v),Z_Y^{}=Z_Y(G^{};q,v^{})$$
(49)
and similarly we write

$$F_X=F_X(G;q,v),F_Y^{}=F_Y(G^{};q,v^{})$$
(50)
where $`F_X`$ is defined in (18). We seek duality relations which express each given $`Z_X`$ in terms of a linear combination of the dual $`Z_Y^{}`$. Recall the Möbius inversion (20) relating $`Z_X`$ to $`F_X`$, it is sufficient to obtain duality relations for the functions $`F_X`$. Indeed, for $`N=1`$, for example, the identity (36) ensures that there is a one-to-one correspondence between the (connected) partitions $`X=\pi (S)`$ on $`G`$ and $`Y=\pi (S^{})`$ on $`G^{}`$. Then, using Lemma 1 one obtains the desired duality relation

$$q^{|X|/2}F_X(G;q,v)=k^{}q^{|Y|/2}F_Y(G^{};q,v^{})N=1$$
(51)
where we have used $`n(X)=n(Y)=1`$, and the coefficient $`k^{}`$ was defined in (32). This duplicates a result of .222It can be readily verified that Eq. (51) is the same as Eq. (49) of , after the substitution of $`F_X=D_X`$ and $`F_Y^{}=q^{|E|/2}D_Y^{}`$.

For $`N2`$, however, the mappings $`XY`$ and/or $`YX`$ are not necessarily one-to-one . See Fig. 2 for an example. Generally, for fixed $`\pi (S)=X,\pi (S^{})=Y`$, we define by analogy to (17) the functions

$`F_{X,Y}(G;q,v)`$ $`=`$ $`{\displaystyle \underset{S𝒮(X),S^{}𝒮(Y)}{}}v^{|S|}q^{p(S)|X|}`$
$`F_{Y,X}(G^{};q,v^{})`$ $`=`$ $`{\displaystyle \underset{S𝒮(X),S^{}𝒮(Y)}{}}(v^{})^{|S^{}|}q^{p(S^{})|Y|},`$ (52)
Then, one has the duality relation

$$q^{|X|/2}F_{X,Y}(G;q,v)=k^{}q^{|Y|/2}F_{Y,X}(G^{};q,v^{})$$
(53)
which is a generalization of (51). The partial partition functions are then

$`Z_X(G;q,v)`$ $`=`$ $`{\displaystyle \underset{X^{}X}{}}{\displaystyle \underset{Y}{}}F_{X,Y}(G;q,v)`$
$`Z_Y(G^{};q,v^{})`$ $`=`$ $`{\displaystyle \underset{Y^{}Y}{}}{\displaystyle \underset{X}{}}F_{Y,X}(G^{};q,v^{})`$ (54)
For $`N=1`$, the $`XY`$ mapping is unique, so that (3.5) becomes (18) and the relation (3.5) can be inverted, permitting one to express $`Z_X`$ as a linear combination of $`Z_Y`$’s. For $`N2`$, however, (3.5) cannot be inverted, and so we cannot express the $`F_{X,Y}`$ in terms of the $`Z_X`$. Hence the duality relations (53) do not help us. As a result, we will instead derive duality relations of the form

$$\underset{X}{}\alpha (X)F_X=k^{}\underset{Y}{}\beta (Y)F_Y^{}$$
(55)
The sum on the left runs over partitions of the root set $`R`$, and the sum on the right side runs over partitions on $`R^{}`$. The coefficients $`\alpha `$ and $`\beta `$ depend on $`q,v`$ as well as their arguments $`X`$ and $`Y`$, and are different for each duality relation.

Using the Möbius inversion (20), (55) provides duality relations for the partial partition functions, and we end up with

$$\underset{W}{}a(W)Z_W=k^{}\underset{U}{}b(U)Z_U^{},$$
(56)
where

$`a(W)`$ $`=`$ $`{\displaystyle \underset{XW}{}}\mu (W,X)\alpha (X)`$
$`b(U)`$ $`=`$ $`{\displaystyle \underset{YU}{}}\mu (U,Y)\beta (Y)`$ (57)
Note that the number of independent relations depends on the root set $`R`$. If all roots lie on one face, so that $`N=1`$, then our result produces the same number of relations as the number of planar partitions of the root set. This duplicates the results of . If all roots lie on distinct faces, so that $`|R|=N`$, then our method produces just one relation, namely the original duality (11).

### 3.6 Correlation duality

We use the two lemmas to generate correlation duality relations. The definition (3.5) was not useful because it involved a sum over all edge sets $`S`$ with fixed partition $`X=\pi (S)`$ and fixed dual partition $`Y=\pi (S^{})`$. As we now show, one way to generate duality relations is by summing instead over all edge sets with fixed local partitions $`X_{\mathrm{loc}}`$ and $`Y_{\mathrm{loc}}`$. According to (37) each of these must be a refinement of the dual of the other, so we can restrict attention to such pairs of local partitions.

Accordingly, writing

$$\xi =X_{\mathrm{loc}},\eta =Y_{\mathrm{loc}}$$
(58)
we say that $`(\xi ,\eta )`$ are compatible if

$$\eta \xi ^{}\mathrm{and}\xi \eta ^{}$$
(59)
Let $`(\xi ,\eta )`$ be compatible local partitions, and define a collection of edge sets on $`G`$ as follows:

$$𝒮(\xi ,\eta )=\{SE:\pi (S)_{\mathrm{loc}}=\xi ,\pi (S^{})_{\mathrm{loc}}=\eta \}$$
(60)
For ease of notation, we write $`𝒮(\xi ,\eta )^{}`$ for the collection of edge sets on $`G^{}`$ which are dual to those in $`𝒮(\xi ,\eta )`$.

In general, for a given pair $`(\xi ,\eta )`$, the set $`𝒮(\xi ,\eta )`$ may be empty. (For example, in the case $`N=1`$ this happens unless $`\eta =\xi ^{}`$). However if $`𝒮(\xi ,\eta )`$ is not empty, then (38) implies that it must have the following form (recall the definition (16))

$$𝒮(\xi ,\eta )=𝒮(X_1)\mathrm{}𝒮(X_m)$$
(61)
for some collection of partitions $`\{X_1,\mathrm{},X_m\}`$. The reason is clear: suppose that $`S𝒮(\xi ,\eta )`$ and that $`\pi (S)=X`$. Let $`T`$ be any other edge set with $`\pi (T)=X`$. Then by definition also $`\pi (T)_{\mathrm{loc}}=\xi `$, and by (38), $`\pi (T^{})_{\mathrm{loc}}=\eta `$. Hence also $`T𝒮(\xi ,\eta )`$, hence $`𝒮(X)𝒮(\xi ,\eta )`$.

We postpone for a moment the question of determining which partitions occur on the right side of (61). First we will use this relation to write the sum over edge sets in $`𝒮(\xi ,\eta )`$ as a sum over factors $`F_X`$ defined in (17). Recalling the left side of (33), and using (61) and (17) we get the following identity:

$`{\displaystyle \underset{S𝒮(\xi ,\eta )}{}}v^{|S|}q^{p(S)[|\pi (S)|n(\pi (S))]/2}`$ $`=`$ $`{\displaystyle \underset{j=1}{\overset{m}{}}}{\displaystyle \underset{S𝒮(X_j)}{}}v^{|S|}q^{p(S)[|\pi (S)|n(\pi (S))]/2}`$
$`=`$ $`{\displaystyle \underset{j=1}{\overset{m}{}}}q^{[|X_j|+n(X_j)]/2}F_{X_j}`$
Next we use our first identity. Using (33) we can rewrite the left side of (3.6) as

$$\underset{S𝒮(\xi ,\eta )}{}v^{|S|}q^{p(S)[|\pi (S)|n(\pi (S))]/2}=k^{}\underset{S^{}𝒮(\xi ,\eta )^{}}{}(v^{})^{|S^{}|}q^{p(S^{})[|\pi (S^{})|n(\pi (S^{}))]/2}$$
(63)
Now we repeat for $`G^{}`$ the argument leading to (61), and obtain

$$𝒮(\xi ,\eta )^{}=𝒮(Y_1)\mathrm{}𝒮(Y_l)$$
(64)
for some collection of partitions $`\{Y_1,\mathrm{},Y_l\}`$ of $`R^{}`$. Hence we end up with the analog of (3.6), namely

$$\underset{S^{}𝒮(\xi ,\eta )^{}}{}(v^{})^{|S^{}|}q^{p(S^{})[|\pi (S^{})|n(\pi (S^{}))]/2}=\underset{i=1}{\overset{l}{}}q^{[|Y_i|+n(Y_i)]/2}F_{Y_i}^{}$$
(65)
Putting together (3.6), (63) and (65) we get

$$\underset{j=1}{\overset{m}{}}q^{[|X_j|+n(X_j)]/2}F_{X_j}=k^{}\underset{i=1}{\overset{l}{}}q^{[|Y_i|+n(Y_i)]/2}F_{Y_i}^{}$$
(66)
This is our new duality relation. As stated earlier, it relates a linear combination of $`F_X`$ to a linear combination of $`F_Y^{}`$.

Now we turn to the question of which partitions can occur on the left-hand and right-hand sides of (66). These are determined by the compatible pair $`(\xi ,\eta )`$ via the relations (61) and (64). In order to decide which partitions occur in (61) and (64), it seems to be necessary to work on a case by case basis, by starting with $`(\xi ,\eta )`$ and examining each planar partition $`X`$ to see if it produces this pair. So partly to hide our ignorance, and partly to tidy up the notation, we define

$$\theta (X;\xi ,\eta )=\{\begin{array}{cc}1\hfill & \text{if }\text{ }SE\text{ such that }\pi (S)=X\text{}\pi (S)_{\mathrm{loc}}=\xi \text{}\pi (S^{})_{\mathrm{loc}}=\eta \hfill \\ 0\hfill & \text{otherwise}\hfill \end{array}$$
(67)
Similarly if $`Y`$ is a partition of $`R^{}`$, then define $`\theta ^{}(Y;\xi ,\eta )=1`$ if there is an edge set $`S^{}`$ with $`\pi (S^{})=Y`$ and also $`\pi (S)_{\mathrm{loc}}=\xi `$, $`\pi (S^{})_{\mathrm{loc}}=\eta `$, and otherwise $`\theta ^{}(Y;\xi ,\eta )=0`$. Then we can rewrite (66) as

$`{\displaystyle \underset{X}{}}\theta (X;\xi ,\eta )q^{[|X|+n(X)]/2}F_X`$
$`=k^{}{\displaystyle \underset{Y}{}}\theta ^{}(Y;\xi ,\eta )`$ $`q^{[|Y|+n(Y)]/2}F_Y^{},`$ (68)
and this is precisely the form described in (55).

Finally we use the Möbius relation (20) to re-express (3.6) in terms of the partial partition functions. The result has the form (56), namely

$$\underset{W}{}a(W)Z_W=k^{}\underset{U}{}b(U)Z_U^{},$$
(69)
where the coefficients are

$$a(W)=a(W;\xi ,\eta )=\underset{XW}{}\mu (W,X)\theta (X;\xi ,\eta )q^{[|X|+n(X)]/2}$$
(70)
and

$$b(U)=b(U;\xi ,\eta )=\underset{YU}{}\mu (U,Y)\theta ^{}(Y;\eta ,\xi )q^{[|Y|+n(Y)]/2}$$
(71)
To summarize, we have obtained a duality relation for every pair of compatible local partitions $`(\xi ,\eta )`$ in the form of (3.6). In some cases the relation is empty; otherwise it is given by (3.6). However, it must be emphasized that while we have obtained new duality relations for the partial partition functions $`Z_X`$, our prescription does not necessarily generate all duality relations (in fact we know that it is incomplete in the case $`N=2`$ as described below). The generation of the complete set of dualities remains an open question.

## 4 Examples

### 4.1 Two external faces, two roots each

This case was examined in detail in  (see in particular Figs. 6 and 7), and we compare our results here with those in . Label the roots $`1,2`$ on one face, and $`3,4`$ on the other. There are four local partitions on $`G`$, namely $`(1)(2)(3)(4)`$, $`(12)(3)(4)`$, $`(1)(2)(34)`$ and $`(12)(34)`$ (we follow standard notation for partitions, so, for example, $`(1)(2)(34)`$ means that there are three blocks, containing roots $`(1)`$, $`(2)`$ and $`(3,4)`$ respectively). Since there are also four local partitions on $`G^{}`$, there are 16 possible pairs $`(\xi ,\eta )`$. However, only 10 of these satisfy (37), so there are 10 possible duality relations (69). Closer examination shows that only 5 of these pairs can occur as local partitions of edge sets. We write $`1^{},2^{},3^{},4^{}`$ for the dual roots. The 5 possible pairs are

$`\xi =(1)(2)(3)(4),\eta `$ $`=`$ $`(1^{}2^{})(3^{}4^{})`$ (72)
$`\xi =(1)(2)(34),\eta `$ $`=`$ $`(1^{}2^{})(3^{})(4^{})`$
$`\xi =(12)(3)(4),\eta `$ $`=`$ $`(1^{})(2^{})(3^{}4^{})`$
$`\xi =(12)(34),\eta `$ $`=`$ $`(1^{})(2^{})(3^{})(4^{})`$
$`\xi =(1)(2)(3)(4),\eta `$ $`=`$ $`(1^{})(2^{})(3^{})(4^{})`$
There are 15 partitions of the 4 roots on $`G`$, and they are all planar, namely, each one can occur as $`\pi (S)`$ for some edge set $`S`$. By our basic result (61) each partition is associated with one of the five pairs in (72). Similarly the 15 partitions of the dual roots are each associated with one pair. In Table 1 we show these associations.

The 5 duality relations are now obtained by substituting into (66). For example, the third pair in (72) gives the identity

$`q^{5/2}`$ $`F_{(12)(3)(4)}+q^{3/2}F_{(123)(4)}+q^{3/2}F_{(124)(3)}`$
$`=`$ $`k^{}[q^{5/2}F_{(1^{})(2^{})(3^{}4^{})}+q^{3/2}F_{(1^{}3^{}4^{})(2^{})}+q^{3/2}F_{(2^{}3^{}4^{})(1^{})}]`$
This reproduces a result in .333It can be verified that Eq. (4.1) is the same as Eq. (29) of , after the substitutions $`T_X=v^{|E|/2}q^{|F|/2+|X|}F_X`$, $`T_X^{}=v^{|E|/2}q^{|F|/2+|Y|}F_Y^{}`$, and renaming roots $`11,23,34,42`$, and dual roots $`51^{},63^{},74^{},82^{}`$.  Similarly the first, second and fourth pairs in (72) reproduce three other identities in . However the last pair in (72) produces an identity which is the sum of two independent identities in .

### 4.2 The general case $`N=2`$

The general case $`N=2`$ is similar to the example above. The present method reproduces many but not all of the duality identities discovered in . In  a non-local block in a partition $`X`$ was called a bridge, and $`X`$ was called a $`k`$-bridge partition if it contained $`k`$ bridges ($`k=0,1,\mathrm{}`$). It was found that there exists one independent duality relation corresponding to every $`0`$-bridge partition, and in addition one independent relation for every $`k`$-bridge partition with $`k2`$.

Our method here reproduces all the duality relations for the $`0`$-bridge partitions, but not for the $`k`$-bridge partitions for $`k2`$. In the latter case it combines $`k`$ independent relations into a single relation. This was the case in the previous example, where the two relations for $`k=2`$ were combined into a single relation.

### 4.3 An example with $`N=3`$

As an illustration of our method we consider the case where there are three external faces, each containing two roots. Roots $`1,2`$ are on face 1, roots $`3,4`$ are on face 2 and roots $`5,6`$ are on face 3. The dual roots are $`1^{},2^{}`$ etc.

The number of local partitions of $`R`$ is 8, since there are two choices on each face. Hence the total number of pairs $`(\xi ,\eta )`$ of local partitions on $`R`$ and $`R^{}`$ is 64. The number of compatible local pairs is 27. After closer analysis, it turns out that only 15 of these pairs can be realised via edge sets. Hence our method produces 15 independent duality relations corresponding to these 15 allowed pairs.

For brevity we present just one of these relations, corresponding to the following local pair:

$$\xi =(12)(3)(4)(5)(6),\eta =(1^{})(2^{})(3^{}4^{})(5^{}6^{})$$
(74)
Instead of writing out the identity, we list in Table 2 the 13 partitions which occur on the left-hand side of (66), and the 9 dual partitions that occur on the right-hand side, along with the exponents of $`q`$.

### 4.4 Bounds on the number of relations

In general it seems to be a hard problem to determine exactly the number of compatible local partitions which can produce duality relations. However our method does provide a lower bound on this number. If $`\xi `$ is a local partition of $`R`$, then certainly $`𝒮(\xi ,\xi ^{})`$ is not empty, and hence there is always a duality relation for this compatible pair. The number of such relations is the product of the number of planar local partitions on the faces. For example, if $`N=2|R|`$ and the roots are paired on the faces, then the number of relations is at least $`2^N`$.

## 5 The Tutte polynomial

In 1955 Tutte  introduced in graph theory the notion of dichromatic, or Tutte, polynomials, which turns out to be precisely the Potts partition function. For our purposes and to conform with notations of , we define the Tutte polynomial associated with a graph $`G`$ as the two-variable polynomials

$`Q(G;t,v)`$ $`=`$ $`v^{|V|}Z(G;vt,v)`$ (75)
$`=`$ $`v^{|V|}{\displaystyle \underset{SE}{}}t^{p(S)}v^{|S|+p(S)}`$
where we have $`t=v^{}=q/v`$. Similarly, one defines as in  the rooted Tutte polynomials

$`Q_X(G;t,v)`$ $`=`$ $`v^{|V|}Z_X(G;vt,v)`$
$`Q_{X,Y}(G;t,v)`$ $`=`$ $`v^{|V|}Z_{X,Y}(G;vt,v)`$ (76)
and the associated summations

$`H_X(G;t,v)`$ $`=`$ $`v^{|V|}F_X(G;vt,v)`$
$`H_Y(G^{};v,t)`$ $`=`$ $`t^{|V^{}|}F_Y(G;vt,t)`$
$`H_{X,Y}(G;t,v)`$ $`=`$ $`v^{|V|}F_{X,Y}(G;vt,v)`$
$`H_{Y,X}(G^{};v,t)`$ $`=`$ $`t^{|V^{}|}F_{Y,X}(G;vt,t)`$ (77)
Then, the duality relation (11) for the partition function becomes

$$vQ(G;t,v)=tQ(G^{(\mathrm{D})};v,t)$$
(78)
and Lemma 1 assumes the form

$$v^{1|V|+|X|n(Y)}\left[v^{|S|}(vt)^{p(S)|X|}\right]=t^{1|V^{}|+|Y|n(X)}\left[t^{|S^{}|}(vt)^{p(S^{})|Y|}\right]$$
(79)
Likewise, the $`N=1`$ correlation duality (51) becomes

$$v^{|X|}H_X(G;t,v)=t^{|Y|}H_Y(G^{};v,t)N=1$$
(80)
and the duality (53) becomes

$$v^{|X|}t^{n(X)1}H_{X,Y}(G;t,v)=t^{|Y|}v^{n(Y)1}H_{Y,X}(G^{};v,t)$$
(81)
Furthermore, the relation (66) for general $`N`$ can be written as

$`v^{1|V|+(|R|N)/2}{\displaystyle \underset{j=1}{\overset{m}{}}}(vt)^{[|Y_i|+n(Y_i)]/2}H_{X_j}`$
$`=t^{1|V^{}|+(|R|N)/2}{\displaystyle \underset{i=1}{\overset{l}{}}}(vt)^{[|Y_i|+n(Y_i)]/2}H_{Y_i}^{}`$ (82)
These expressions reflect the symmetric roles played by the variables $`v`$ and $`t=v^{}=q/v`$.

Acknowledgments

This research was supported in part by a Research Scholarship Development Fund (RSDF) grant from Northeastern University. In addition CK was partially supported by NSF Grant DMS-9705779, and FYW was partially supported by NSF Grant DMR-9980440. We thank W. T. Lu for help in preparing the figures.

Figure Captions

Fig. 1. An $`N=3`$ example of a partition $`X`$ and the associated local partition $`X_{\mathrm{loc}}`$. (a) $`X=(135)(247)(68)`$, (b) $`X_{\mathrm{loc}}=(1)(35)(2)(4)(7)(68)`$. Cross-hatched areas denote external faces; lines denote connected edge sets.

Fig. 2. An $`N=2`$ example showing a given partition $`X=(1)(23)(4)`$ (solid line) can be the dual of two different partitions $`Y_1`$ and $`Y_2`$ (broken lines).
(a) $`Y_1=(1^{}2^{}3^{}4^{})`$, (b) $`Y_2=(1^{}2^{})(3^{}4^{})`$.


format: text

text:
# Floydian trajectories for stationary systems: a modification for bound states

## I Introduction

In a series of papers \[1, and references therein\], Floyd has developed a trajectory representation of quantum mechanics which, while related to the Bohm interpretation, is distinct in both in its trajectories and its interpretation. Whilst the Floydian method has been strongly advocated and well exemplified by its proponent, in comparison with the Bohm interpretation it seems not to have received as wide critical attention in the literature as perhaps it deserves. Here we attempt to partly redress the balance by briefly reviewing the Floydian method and then critically examining the basis of the kinematic equation of its trajectories. We refer to the Bohm interpretation to stimulate the critical analysis. Floyd has made a detailed comparison of the two approaches in .

Two novel aspects of the Floydian method are: 1) the identification of microstates of the Schrödinger equation and 2) a kinematic equation for trajectories which is derived from the Hamilton-Jacobi method of classical mechanics. The method is based upon the direct solution of the Quantum Stationary Hamilton-Jacobi Equation (QSHJE) rather than upon the solution of the corresponding Schrödinger equation and has, thus far, only been developed for 1-dimensional and pseudo-2-dimensional stationary systems. In particular, it has been applied to the finite potential step, the semi-infinite step (in 2-dimensions) , the infinite  and the finite square potential well . In early work, solutions of the QSHJE were obtained numerically  but following the discovery of the form of a general solution in 1-dimension , attention has been directed towards studying analytical solutions. The results can be extended to higher dimensions if separation of variables is permitted. However, closed form solutions are not known in the general case and whilst it has been briefly discussed in principle , the general case of higher dimensions remains untried in practice.

More recently, the work of Floyd has been set in a wider context by the work of Faraggi and Matone  who derive and generalise (to higher dimensions) the QSHJE on the basis of a postulated equivalence principle: that all physical systems are equivalent under coordinate transformations. This connection is of interest, not only because it provides a wider context for the work of Floyd but also because it motivates a connection between the quantum potential and coordinate transformations or deformations of the geometry of space.

Whilst we briefly review the Floydian method and the basis of the appearance of microstates of the Schrödinger equation, our primary focus is upon its trajectory representation. We critically assess the application of the Floydian trajectory method to bound systems. By working in a 3-dimensional coordinate representation we also highlight necessary further assumptions for its application in higher dimensions. This is in contrast to the principally 1-dimensional approach of Floyd.

A point of particular attention, is the appearance of the energy derivative of the kinetic energy in the kinematic equation of Floydian trajectories, especially in quantum systems with discrete energy spectra. Such a derivative is well defined (in fact it is unity) for classical conservative systems, since the transition between one ’stationary’ energy state of a given system to another is continuous. However, in bound quantum systems, not only is the energy spectrum discrete (complicating the definition of an energy derivative), but also any attempt to make a transition between the states in order to define a continuous derivative with respect to energy whilst at the same time continuing to satisfy the Schrödinger equation and its boundary conditions faces two problems: 1) states between the eigenstates are time dependent so resulting in a time and space dependent energy derivative 2) there is a multiplicity of such derivatives from which to choose. We explore these problems, (thus far investigated by Floyd only numerically with respect to the harmonic oscillator ), for stationary systems in general and examine their relevance to the energy variation of the QSHJE as used by Floyd. These issues do not arise in systems with continuous energy spectra where the Floydian method has been used to investigate barrier reflection and tunnelling problems.

The following sections proceed as follows. By way of establishing a comparative base, we first briefly review the equations of the Bohm trajectory method for the general case and subsequently for stationary systems. After briefly reviewing the work of Faraggi and Matone, we then review the Floydian method and the appearance of microstates contrasting the approach with that of Bohm. Deferring until later the kinematic equation of Floydian trajectories, we then derive, for its use, the variational derivative of the kinetic energy with respect to the total energy for stationary systems revealing its space-time dependence and its ambiguity. Taking a more general approach than Floyd, we then derive the kinematic equation of a Floydian trajectory from the variational derivative of the action with respect to the total energy, described by Floyd as the epoch $`\tau `$. We then use our results to demonstrate that Bohmian and Floydian trajectories are related by different assumptions about $`\tau `$ and we are led to the idea of time deformation by the quantum potential. After briefly comparing applications to different types of potential, we finally draw together our results in the conclusions.

## II The Bohm trajectory method

In the Bohm approach  the quantum Hamilton-Jacobi equation is derived, along with the equation for conservation of probability, from the Schrödinger equation and its complex conjugate by using the polar form of the wave function $`\psi =R(q,t)\mathrm{exp}(i\frac{S(q,t)}{\mathrm{}})`$. For stationary systems, $`R=R(q)`$ and $`S(q,t)=W(q)Et`$, where $`W(q)`$ is Hamilton’s characteristic function. In the Bohmian method, one obtains $`S(q,t)`$ directly from the wave function solution of the Schrödinger equation. The trajectory kinematics are derived from a classical form of the probability current density identified in the probability conservation equation. This method is therefore rooted in the probabilistic interpretation of the wave function.

In the coordinate representation, the substitution $`\psi =R(q,t)\mathrm{exp}(i\frac{S(q,t)}{\mathrm{}})`$ in the Schrödinger equation yields:

$$\frac{1}{2m}\left(S\right)^2+V(q)+Q(q)=\frac{S}{t}$$
(1)
and

$$\frac{R^2}{t}+\left(\frac{R^2}{m}S\right)=0$$
(2)
where

$$Q=\frac{\mathrm{}^2}{2mR}^2R$$

is the quantum potential. From the probability conservation equation (2), one identifies the probability current density $`j=\frac{R^2}{m}S`$. Using the classical form $`j=\rho u`$ (where $`\rho `$ is the density and $`u`$ the velocity of a flow) one then identifies $`\dot{q}=\frac{1}{m}S`$ as the streamline velocity of the probability density current which upon integration yields a trajectory. If the density of such trajectories is at any time $`R^2`$, then equation (2) ensures that it is always so. We also observe from equation (1), that with the momentum defined by $`p=S`$, this equation is the Hamilton-Jacobi equation but for the quantum potential $`Q(q)`$. One may therefore consistently assign $`p=m\dot{q}`$. The latter assignment is strictly true only for conservative systems in a cartesian coordinate representation.

In the case of stationary systems, $`S(q,t)=W(q)Et`$ and so the above equations reduce to

$$\frac{1}{2m}\left(W\right)^2+V(q)+Q(q)=E$$
(3)
and

$$\left(\frac{R^2}{m}W\right)=0$$
(4)
with the trajectory equation

$$\dot{q}=\frac{1}{m}W.$$
(5)
From the above, we see that for a stationary system of a given energy E, a Bohm trajectory is defined by a square integrable solution of the Schrödinger equation (to secure the concept of probability) and an initial condition $`x_0`$. Moreover, owing to the single-valuedness of the wave function, Bohm trajectories for a given eigenstate do not cross. For stationary systems with real wave functions, Bohm trajectories are stationary; a point of concern to Einstein who observed that, in the classical $`\mathrm{}0`$ limit, the classical motion is not obtained , p243. Floyd  has observed that the use of a trigonometric (rather than complex polar) ansatz resolves this artifact. We do not pursue this further here.

## III The quantum stationary Hamilton-Jacobi equation

### A The quantum stationary Hamilton-Jacobi equation via an equivalence principle

More recently than Floyd’s early work, Faraggi and Matone  have independently derived (and generalised to higher dimensions) the QSHJE from a postulated equivalence principle: that all physical systems are equivalent under coordinate transformations. This work extends the implications of Floyd’s work on the QSHJE to a connection between quantum mechanics and gravity . In the following, we summarise those key points of the work of Matone et al. that have a bearing on our work here:

* A key corollary of the equivalence postulate is that there must exist a trivialising coordinate transformation $`qq_0`$, defined by the relation between the respective Hamilton’s characteristic functions (reduced actions) $`W(q)=W_0(q_0)`$, which reduces any system to that of the free particle :

$$𝒲(q)𝒲_0(q_0)$$

where $`𝒲=V(q)E`$ and $`𝒲_00`$.
* The quantum potential term in the QSHJE is a direct consequence of the postulated universal existence of the trivialising coordinate transformation. In particular, the absence of the quantum potential (as in classical mechanics) prohibits the existence of such a transformation to a frame in which a system is ’at rest’. The existence of the quantum potential therefore enables the removal of this privileged frame.
* The quantum potential is expressed in terms of the Schwarzian derivative of the characteristic function or reduced action $`W_\mu `$ (see section III C) which must therefore never be a constant in order for the trivialising map to exist. To guarantee this, a general solution of the Schrödinger equation should be expressed in the form $`\psi =R_\mu (A\mathrm{exp}(iW_\mu /\mathrm{})+B\mathrm{exp}(iW_\mu /\mathrm{}))`$ rather than in the Bohmian form $`\psi =R\mathrm{exp}(iW/\mathrm{})`$ for real $`R`$, $`W`$, $`R_\mu `$, and $`W_\mu `$ with $`A`$ and $`B`$ as complex constants. In the case of real wave functions, this implies that $`A`$ and $`B`$ must be of equal magnitude and, modulo an integer factor $`2\pi `$, of opposite signed phase so that a real wave function is the sum of a Bohmian polar form and its complex conjugate. (Here and hereafter, we distinguish the characteristic function, $`W_\mu `$, as a direct solution of the QSHJE from the the reduced action, W, appearing in the phase of the wave function, as used in the Bohm formulation.)
* A general principle is that a general solution of the Schrödinger equation is always has two independent solutions. In general, they are not square integrable. However, under the special conditions for eigenfunctions, one and only one of them may be a square-integrable function on the real line. The Copenhagen interpretation selects the latter solution in order that the wave function be interpreted as a probability amplitude.
* In the case of the QSHJE, a condition for the existence of the Schwarzian derivative is that the total energy $`E`$ must take a value such that the QSHJE corresponds to a Schrödinger equation having a wave function that is a square-integrable function on the real line ie $`E`$ is an eigenvalue of the corresponding Schrödinger equation. Thus, in solving the QSHJE, quantisation arises independently of the Copenhagen interpretation demand for square-integrable solutions.
* Time does not appear explicitly in the QSHJE and so the time parameterisation is introduced, after the solutions of the QSHJE are obtained, in the form of trajectories defined using the Floydian ansatz $`t\tau =W_\mu /E`$. In classical mechanics this is equivalent to the kinematic relation $`p=m\dot{q}`$ in which time parameterisation is implicit.

The work of Faraggi et al. draws upon the requirement for satisfactory classical limits (expressed as $`\mathrm{}0`$) and through its insistence upon a non-constant reduced action, claims to have side-stepped Einstein’s objection that the Bohm interpretation does not have the correct classical limit (in the aforementioned sense) for bound systems such as the quantum harmonic oscillator. We also see that the quantum potential derived from the equivalence principle and used in Floyd’s trajectory method does not necessarily equal that of the Bohm interpretation, for the latter is based upon the polar decomposition of only square integrable wave-functions whereas the former are based on the direct solution of the QSHJE and correspond to more general solutions of the Schrödinger equation.

In the next section we focus upon the specific treatment of the QSHJE by Floyd.

### B Overview of the Floyd method

The Floydian approach  originated from an alternative ansatz to the Bohm polar form for solutions of the Schrödinger equation, which is nevertheless consistent with the second Bohm (probability conservation) equation (4). Substituted in the Schrödinger equations, the Floyd ansatz results in a single third order quantum stationary Hamilton-Jacobi (QSHJE) equation for the quantum form of Hamilton’s characteristic function $`W_\mu (q)`$. For a given energy eigenvalue of the Schrödinger equation, $`W_\mu (q)`$ is not a unique solution of the QSHJE and this leads to the idea that the square integrable wave function for bound states is not exhaustive but has associated microstates. As seen in the previous section, this conclusion has been more generally substantiated by the work of Farragi and Matone. Whilst the Floydian ansatz  is consistent with the equation for the conservation of probability, its trajectories are not distributed in accordance with the density of the wave function as in the case of the Bohm interpretation. Indeed, unlike Bohmian trajectories, the density of Floydian trajectories would subsequently not be that of the wave function. Rather, Floyd claims that the necessity of the probabilistic interpretation of the wave function in orthodox quantum mechanics arises from the incomplete specification of a particular microstate. Moreover, each trajectory corresponds to an individual microstate of the Schrödinger equation but is nevertheless sufficient to specify the corresponding wave function . We defer to a later section the derivation of the kinematic equation for Floydian trajectories and here more explicitly outline the work on the quantum stationary Hamilton-Jacobi equation and its solutions in one dimension.

### C The quantum stationary Hamilton-Jacobi equation and its solutions

The quantum stationary Hamilton-Jacobi equation (QSHJE) is given in one dimension $`q`$ by

$$\frac{(W_\mu ^{})^2}{2m}+VE=\frac{\mathrm{}^2}{4m}W_\mu ;q$$
(6)
where $`W_\mu ^{}=\frac{W_\mu }{q}`$ is the momentum conjugate to $`q`$ and $`W_\mu ;q`$ is the Schwarzian derivative of $`W_\mu `$ with respect to $`q`$. The Schwarzian derivative,

$$W_\mu ;q=\frac{W_\mu ^{\prime \prime \prime }}{W_\mu ^{}}\frac{3}{2}\left(\frac{W_\mu ^{\prime \prime }}{W_\mu ^{}}\right)^2,$$

which is unaffected by the sign of $`W_\mu `$, means that the QSHJE is a third-order nonlinear differential equation. The left hand side of equation (6) replicates the classical Hamilton-Jacobi equation, whereas the right hand side embodies the quantum effects in the Schwarzian derivative. Comparing equation (6) with equation (3) one may identify the quantum potential term as

$$Q=\frac{\mathrm{}^2}{4m}W_\mu ;q.$$
(7)
However, we recognise that the Bohm quantum potential is not always equal to the Schwarzian derivative form. The latter is only guaranteed when $`W_\mu =W`$. Equation (6) can be obtained from the Bohm formulation for stationary systems by recognising that, in one-dimension, equation (4) has a solution $`R=(Cm)/\sqrt{|W^{}|}`$, where C is a constant. (Note that this is not generally the case in higher dimensions). Substituting for $`R`$ in equation (3) one obtains the form of equation (6). This approach corresponds to substituting $`\psi =\mathrm{exp}(\pm iW_\mu /\mathrm{})/\sqrt{|W_\mu ^{}|}`$ (or a linear combination of the latter) into the Schrödinger equation. Floyd discovered that the general solution, $`W_\mu ^{}`$, of the QSHJE is given by

$$W_\mu ^{}=(2m)^{1/2}(a\varphi ^2+b\theta ^2+c\varphi \theta )^1$$
(8)
where $`[a,b,c]`$ is a set of real coefficients such that $`a,b>0`$, and $`(\varphi ,\theta )`$ must be a pair of (normalised) independent solutions of the stationary Schrödinger equation, $`\mathrm{}^2\psi ^{\prime \prime }/(2m)+(VE)\psi =0`$ for given energy $`E`$. (Neither of the the solutions $`(\varphi ,\theta )`$ is necessarily square integrable though, as mentioned in section III A, in the case of eigenstates of bound systems one and only one is.) For given $`E`$, $`W_\mu `$ is therefore prescribed by the three constants $`[a,b,c]`$. Floyd has observed  that these constants of the motion may alternatively be represented by the set $`[x_0,\dot{x_0},\ddot{x_0}]`$ of trajectory initial conditions. This is to be contrasted with the set $`[x_0]`$ required in classical mechanics to specify $`W_\mu `$ for a given $`E`$ and is a reflection of the third order nature of the QSHJE. Given the time parameterisation provided by the kinematic equation of the Floydian trajectories, each trajectory therefore specifies a unique microstate $`W_\mu `$ and vice-versa. Floydian trajectories for a given eigenstate are therefore correspondingly more various than Bohm trajectories. Indeed, unlike Bohm trajectories, Floydian trajectories (for different microstates) of an eigenstates may cross.

How do the microstates arise? Floyd has shown  that for energy eigenstates of stationary bound systems in one dimension (for which the wave function is real), equation (6) has singular values at precisely the locations where the boundary values are applied. This is because $`W_\mu ^{}0`$ as $`x\pm \mathrm{}`$ since, whilst one of the independent solutions of Schrödinger equation is square-integrable, the other must be unbounded at $`x\pm \mathrm{}`$ regardless of $`(a,b,c)`$ given their constraints. Thus, boundary conditions $`W_\mu ^{}(x\mathrm{})=0`$ admit non-unique solutions for $`W_\mu ^{}`$ and so yield an infinity of trajectories. This conclusion also applies to semi-bound systems (eg potential barriers where the total energy is less than the barrier height). We observe that the balance of energy between the kinetic energy and quantum potential terms will differ between the microstates of an eigenstate of such a bound system. Floyd  also provides evidence to suggest that an uncountable number of microstates is also possible for bound states in dimensions greater than one, though does not demonstrate that this is the case.

For unbound systems, in which the (complex) solution of the Schrödinger equation is an initial value problem, the initial conditions determine a unique solution to equation (6) for $`W_\mu ^{}`$. In this case, the same solution is obtained in the Bohm interpretation ie $`W_\mu =W`$. Thus, there are no microstates in this case and the trajectories, for given initial conditions, are unique. However, because their kinematic equations differ, the Floydian and Bohmian trajectories follow different paths through space-time .

### D General form of the wave function used for Floydian trajectories

Given a solution $`W_\mu `$ of the QSHJE for a stationary state of energy $`E`$, the corresponding wave function $`\mathrm{\Psi }`$ may written in the general bi-polar form described in section III A:

$$\mathrm{\Psi }=\psi _\mu ^++\psi _\mu ^{}$$
(9)
where

$$\psi _\mu ^\pm =A_\mu ^\pm R_\mu \mathrm{exp}(\frac{\pm iW_\mu }{\mathrm{}})\mathrm{exp}(i\frac{Et}{\mathrm{}}).$$
(10)
$`A_\mu ^\pm `$ are complex constants and $`R_\mu ^21/|W_\mu |`$ as discussed in section III C. This form is motivated by the fact that the QSHJE is invariant under a sign change of $`W_\mu `$, so that both $`\psi _\mu ^\pm `$ and any linear combination of them lead to the same QSHJE and continuity equation when substituted into the Schrödinger equation. (For a given stationary state, $`\mathrm{\Psi }`$ may, of course, also be expressed in the form of a synthesized wave  using the Bohmian ansatz $`\mathrm{\Psi }=\stackrel{~}{R_\mu }\mathrm{exp}(\frac{i\stackrel{~}{W_\mu }}{\mathrm{}})\mathrm{exp}(i\frac{Et}{\mathrm{}})`$. Selection of the proper form depends on the physical situation.) The case $`A_\mu ^{}=0`$ can represent the complex eigenfunctions of unbound stationary systems whereas $`A_\mu ^{}=A_\mu ^+`$ can represent the real eigenfunctions of bound systems. Moreover, in both instances the equivalence principle requirement that $`W_\mu constant`$ is fulfilled. Thus, the general wave function $`\mathrm{\Psi }`$ comprises two waves, $`\psi _\mu ^\pm `$, of identical form and energy, travelling in opposite directions but scaled by different different complex amplitudes $`A_\mu ^\pm `$ respectively. We refer to these two waves as the positive and negative running waves in correspondence with the sign prefix of $`W_\mu `$ in their phase factors. An important aspect for the next section is that the momentum and energy in the QSHJE can be derived from either of such a pair of running waves using $`p_\mu =S_\mu `$ and $`E=S_\mu /t`$ where $`S_\mu =W_\mu Et`$.

## IV The total energy variational derivative of kinetic energy for stationary systems

We now investigate the variation of kinetic energy with respect to the total energy for stationary systems. This variational dependence is used in the definition of the kinematic equation for Floydian trajectories and in their connection with Bohmian trajectories. Though Floyd actually worked with the variation of the quantum potential  with respect to the total energy, it is simply related to the variation of kinetic energy with respect total energy. Here we find it more convenient to work with the latter.

The use of an energy derivative implies a continuous spectrum of admissible energy eigenvalues. Indeed, Floyd provides examples of Floydian trajectories (using the energy derivative of the quantum potential) for systems with continuous energy spectra for which such a derivative is easily defined . In his work on bound states, in which microstates of the Schrödinger equation are identified, Floyd  gives examples of microstate trajectories in mechanical phase space for the quantum harmonic oscillator, based on a power series solution (in the space variable) for the modified potential (classical plus quantum potential) which has energy depend coefficients. For such systems, with discrete spectra, the concept of the energy derivative is unnatural, for the energy is not obviously a continuous variable. Indeed, in section III A we saw that a condition for the existence of the Schwarzian derivative was that the energy correspond to an eigenstate of the Schrödinger equation. It is therefore necessary to modify the methods in  to take account of the time dependence which an energy variation between eigentstates might induce. Thus, instead of Floyd’s $`Q/E`$, we adopt the notation $`\delta Q/\delta E`$ described as a variational derivative with respect to energy.

In the following sub-sections, we investigate the energy variational derivative of the kinetic energy $`\delta T/\delta E`$ from the standpoint of stationary quantum systems with distinct energy eigenstates. The derivative (justifiably described as partial) for continuous energy spectra is then obtained as a limiting case of the discrete spectrum result.

### A Stationary state variation

Let $`\mathrm{\Psi }_i`$ be a stationary solution of the Schrödinger equation

$$i\mathrm{}\frac{\mathrm{\Psi }_i}{t}=H\mathrm{\Psi }_i$$
(11)
for a given potential $`V(q)`$ and set $``$ of boundary (or initial) conditions. So that we remain within the principles of quantum mechanics and stay with a description of the same physical system, we impose upon the varied wave function, $`\mathrm{\Psi }=\mathrm{\Psi }_i+\delta \mathrm{\Psi }`$, the condition that it also satisfies the Schrödinger equation

$$i\mathrm{}\frac{\mathrm{\Psi }}{t}=H\mathrm{\Psi },$$
(12)
for the same potential $`V(q)`$ and conditions $``$. We also impose the further condition that the varied wave function preserve the space normalisation of the original wave function but accept that the varied wave function will not necessarily be a stationary solution of the Schrödinger equation.

Expressing $`\mathrm{\Psi }_i`$ and $`\mathrm{\Psi }`$ in the bi-polar form of section III D, we obtain their respective Hamilton-Jacobi equations whose difference may be written as

$$\delta \left(\frac{1}{2m}\left(S_\mu \right)^2\right)+\delta Q_\mu =\delta \left(\frac{S_\mu }{t}\right).$$
(13)
where $`S_\mu =\pm W_\mu Et`$. This equation may be expressed in the simple form

$$\delta T_\mu +\delta Q_\mu =\delta E$$
(14)
where, by identification of terms in equation (13), $`T_\mu `$ is the kinetic energy and $`E`$ the total energy. This is a general equation for variations satisfying the constraints of the Schrödinger equation regardless of boundary conditions and normalisation. Thus, under such general conditions, we may have $`E=E(q,t)`$. In the following, the variational derivative $`\delta T_\mu /\delta E`$ is obtained from the ratio of the variations $`\delta T_\mu `$ and $`\delta E`$ with respect to a particular form of infinitesimal variation of the wave function $`\delta \mathrm{\Psi }`$ which satisfies (1) the Schrödinger equation, (2) the conditions $``$ and (3) preserves the normalisation of the wave function.

All three conditions on the varied wave function $`\mathrm{\Psi }`$ are satisfied if it is expressed in the form

$$\mathrm{\Psi }=\mathrm{cos}(\theta )\mathrm{\Psi }_i+\mathrm{sin}(\theta )\mathrm{\Psi }_j$$
(15)
where $`\theta `$ is the parameter of the variation and $`\mathrm{\Psi }_j`$ is another normalised stationary solution (eigenfunction) of the Schrödinger equation for the same potential $`V(q)`$ and conditions $``$. The orthonormality of $`\mathrm{\Psi }_i`$ and $`\mathrm{\Psi }_j`$ guarantees the $`\theta `$-independence of the normalisation of $`\mathrm{\Psi }`$. That $`\mathrm{\Psi }`$ is a linear combination of solutions of the Schrödinger equation for the conditions $``$ means that it also is a solution of the same. In principle, one can construct any solution as a linear combination of any number of eigenfunctions. However, the objective in this case is to study the relative variation of the kinetic and total energies under an infinitesimal change of the wave function from one eigenstate towards another whether or not they have a continuous energy spectrum; hence the form of equation (15). Whilst this analysis could be carried out for any pair of eigenfunctions $`\mathrm{\Psi }_i`$ and $`\mathrm{\Psi }_j`$, here, we define $`\mathrm{\Psi }_i`$ and $`\mathrm{\Psi }_j`$ to be non-degenerate eigenstates corresponding respectively to adjacent different eigenvalues $`E_i`$ and $`E_j`$ in the ordered energy spectrum, so that $`E_j>E_i`$. In this way, we remove ambiguity from the definition of the variational derivative with respect to total energy and bring it into correspondence with the partial derivative $`T/E`$ in classical conservative systems.

As discussed in section III D, $`\mathrm{\Psi }_i`$, $`\mathrm{\Psi }_j`$ and thus $`\mathrm{\Psi }`$ may be decomposed into positive and negative running waves so that $`\mathrm{\Psi }=\psi _\mu ^++\psi _\mu ^{}`$ with

$$\psi _\mu ^\pm =\psi _{\mu i}^\pm \mathrm{cos}\theta +\psi _{\mu j}^\pm \mathrm{sin}\theta $$
(16)
where

$$\psi _{\mu k}^\pm =A_k^\pm R_{\mu k}\mathrm{exp}(\frac{\pm iW_{\mu k}}{\mathrm{}})\mathrm{exp}(i\frac{E_kt}{\mathrm{}})$$
(17)
are themselves the positive and negative running waves for the microstate solution $`W_{\mu k}`$ of the QSHJE for the $`k`$th eigenstate with eigenvalue $`E_k`$.

### B Analysis of kinetic energy variational derivative

In this section, for the sake of notational clarity, we initially drop the microstate suffix <sub>μ</sub> on the variables $`p`$, $`S`$, $`T`$, $`\psi `$ and $`Q`$ and re-introduce them in the final result.

Consider a variation of the wave function parameterised by $`\delta \theta `$. From equations (13) and (14), the change in the kinetic energy is

$$\delta T=\frac{p}{m}\frac{p}{\theta }\delta \theta $$
(18)
where $`p=S`$ is the conjugate momentum. Similarly, the change in the total energy is

$$\delta E=\frac{E}{\theta }\delta \theta $$
(19)
where $`E=S/t`$ is the total energy. Thus, in this work, the variational derivative of the kinetic energy with respect to the total energy is defined by

$$\frac{\delta T}{\delta E}=\frac{1}{m}\frac{pp/\theta }{E/\theta }.$$
(20)
In both the latter equations the partial derivatives with respect to $`\theta `$ are evaluated at $`\theta =0`$ ie in the unvaried state. In passing, we note from equation (14) that $`\delta Q/\delta E=1\delta T/\delta E`$.

The subsequent analysis proceeds by respectively expressing the conjugate momentum and total energy in terms of one of the running waves (we choose the positive wave) of the varied wave function (see equation (16)), so that

$`p=S=\mathrm{}\left({\displaystyle \frac{\mathrm{}}{\psi }}\psi \right)`$ (21)
$`E={\displaystyle \frac{S}{t}}=\mathrm{}\left({\displaystyle \frac{\mathrm{}}{\psi }}{\displaystyle \frac{\psi }{t}}\right)`$ (22)
in which, and hereafter, the $`\pm `$ superscript has been dropped. Therefore,

$`{\displaystyle \frac{p}{\theta }}=\mathrm{}\mathrm{}\left[{\displaystyle \frac{1}{\psi }}{\displaystyle \frac{\psi }{\theta }}{\displaystyle \frac{1}{\psi }}\psi +{\displaystyle \frac{1}{\psi }}{\displaystyle \frac{}{\theta }}\psi \right]`$ (23)
$`{\displaystyle \frac{E}{\theta }}=\mathrm{}\mathrm{}\left[{\displaystyle \frac{1}{\psi }}{\displaystyle \frac{\psi }{\theta }}{\displaystyle \frac{1}{\psi }}{\displaystyle \frac{\psi }{t}}+{\displaystyle \frac{1}{\psi }}{\displaystyle \frac{}{\theta }}{\displaystyle \frac{\psi }{t}}\right]`$ (24)
where the partial derivative terms with respect to $`\theta `$ at $`\theta =0`$ are computed from equation (16) (without the superscripts) as

$`\left({\displaystyle \frac{1}{\psi }}{\displaystyle \frac{\psi }{\theta }}\right)_{\theta =0}={\displaystyle \frac{\psi _j}{\psi _i}}`$ (25)
$`\left({\displaystyle \frac{1}{\psi }}{\displaystyle \frac{}{\theta }}\psi \right)_{\theta =0}={\displaystyle \frac{1}{\psi _i}}\psi _j`$ (26)
$`\left({\displaystyle \frac{1}{\psi }}{\displaystyle \frac{}{\theta }}{\displaystyle \frac{\psi }{t}}\right)_{\theta =0}={\displaystyle \frac{1}{\psi _i}}{\displaystyle \frac{\psi _j}{t}}.`$ (27)
With these results equations (24) become

$`\left({\displaystyle \frac{p}{\theta }}\right)_{ij}=\mathrm{}\mathrm{}\left[{\displaystyle \frac{\psi _j}{\psi _i}}\left({\displaystyle \frac{\psi _j}{\psi _j}}{\displaystyle \frac{\psi _i}{\psi _i}}\right)\right]`$ (28)
$`\left({\displaystyle \frac{E}{\theta }}\right)_{ij}=\mathrm{}\mathrm{}\left[{\displaystyle \frac{\psi _j}{\psi _i}}\left({\displaystyle \frac{1}{\psi _j}}{\displaystyle \frac{\psi _j}{t}}{\displaystyle \frac{1}{\psi _i}}{\displaystyle \frac{\psi _i}{t}}\right)\right],`$ (29)
in which the partial derivatives with respect to $`\theta `$ at $`\theta =0`$ have been labelled by the suffices $`ij`$ to show that they are evaluated at the eigenstate $`i`$ ‘in the direction of’ the eigenstate $`j`$. Expressing the running wave for the $`k`$th eigenstate in the form of equation (17) and re-introducing the microstate suffix $`\mu `$, equations (27) may be used in equation (20) to finally give

$`\left({\displaystyle \frac{\delta T}{\delta E}}\right)_{\mu ij}={\displaystyle \frac{S_{\mu i}}{m\mathrm{\Delta }E_{ji}}}`$ (30)
$`\left[(\mathrm{\Delta }S_{\mu ji})+\mathrm{}\mathrm{tan}\left({\displaystyle \frac{\mathrm{\Delta }S_{\mu ji}}{\mathrm{}}}+\mathrm{\Delta }\alpha _{\mu ji}\right)\mathrm{ln}\left({\displaystyle \frac{R_{\mu j}}{R_{\mu i}}}\right)\right],`$ (31)
where $`\mathrm{\Delta }E_{ji}=E_jE_i`$, $`\mathrm{\Delta }S_{\mu ji}=S_{\mu j}S_{\mu i}`$ and $`S_{\mu k}(q,t)=W_{\mu k}(q)E_kt`$. $`\mathrm{\Delta }\alpha _{\mu ji}=\alpha _{\mu j}\alpha _{\mu i}`$ is the difference between the phases of the running wave scale factors of the $`j`$th and $`i`$th eigentsates for microstate <sub>μ</sub>, where $`A_{\mu k}=|A_{\mu k}|\mathrm{exp}(i\alpha _{\mu k})`$.

We recall that this analysis was done for the positive running wave but observe that the result (31) is invariant to a change of sign prefix to $`W_\mu `$ (for the negative running wave) as expected if the scale factor phases also undergo a corresponding sign change. The latter condition does indeed occur when the eigenstates are real as in the case of bound states. (See section III D.) Equation (31) clearly shows that $`\left(\delta T/\delta E\right)_{\mu ij}`$ may be explicitly time dependent on account of the ’beating’ of the different frequencies of the states $`i`$ and $`j`$.

### C The kinetic energy variational derivative for unbound systems

In the limit of unbound systems (including barriers), which have continuous energy spectra, $`\mathrm{\Delta }S_{\mu ji}`$ and $`\mathrm{\Delta }E_{ji}`$ become infinitesimals and $`R_{\mu j}R_{\mu i}`$ so equation (31) becomes

$$\left(\frac{\delta T}{\delta E}\right)_{\mu ij}\left(\frac{\delta T}{\delta E}\right)_\mu =\frac{S_\mu }{m}\left(\frac{S_\mu }{E}\right),$$
(32)
where the right hand side is identical to the classical form.

In the case of a free system in a constant potential, the spectrum is continuous and $`R`$ is a constant so that the quantum potential is zero. We may then replace $`\delta T/\delta E`$ by $`T/E`$ in equation (32) and obtain precisely the classical result. Clearly, for unbound systems, the energy derivative of the kinetic energy has no explicit time dependence.

### D The kinetic energy variational derivative for bound systems

In the case of stationary bound systems with discrete spectra and real eigenfunctions, equation (31) may be applied to any one of the corresponding Floydian microstates of the Schrödinger equation. Referring to equation (8) and the text immediately thereafter, a microstate may be identified by the set $`[q_0,\dot{q_0},\ddot{q_0},E_i]`$ of trajectory initial conditions and the energy eigenvalue. Thus, to determine from equation (31) the kinetic energy variational derivative for the microstate of the $`i`$th eigenstate defined by the initial conditions $`[q_0,\dot{q_0},\ddot{q_0}]`$, we may define $`\mathrm{\Delta }W_{\mu ji}(q)=W_\mu (q,[q_0,\dot{q_0},\ddot{q_0},E_j])W_\mu (q,[q_0,\dot{q_0},\ddot{q_0},E_i])`$ holding the trajectory initial conditions constant. For such microstates, the conjugate momentum $`S_{\mu i}`$ is non-zero and the second term within the square bracket of equation (31) generates the explicit time dependence of $`\left(\delta T/\delta E\right)_{\mu ij}`$ with a frequency equal to the difference between the frequencies of the states $`i`$ and $`j`$. We therefore observe that for a given set of initial conditions there is a family of explicitly time dependent trajectories corresponding to the family of discrete eigenstates. The explicit time dependence of such bound state microstate trajectories does not appear in the work of Floyd.

## V A trajectory equation from the energy variational derivative of the action

### A The nature of the energy variational derivative

Here we study the basis of the Floydian trajectory method as defined for stationary systems. The approach, in the augmented domain of the motion $`(q,t,E)`$, is based upon the variation of the action with respect to changes in the total energy $`E`$ which are out of the constant energy plane of the motion. As such, it is restricted to stationary systems. We emphasise that the kinematic equations so obtained are for Floydian trajectories in the $`(q,t)`$ domain.

Here we make a point of clarification. Carroll  connects the evaluation of the total energy derivative of the action with the Legendre transformation from the $`(q,t)`$ domain to the $`(q,E)`$ domain. In the case of stationary systems, such a transformation corresponds to the substitution of one independent variable, $`t`$, by another, $`E`$, which is a single valued constant of the motion. Thus, the trajectory in $`(q,E)`$ space would trivially lie on a fixed energy plane labelled by $`E`$. (In the case of non-stationary systems the Legendre transformation is well defined and the motion is not restricted to such energy planes.) We emphasise that the Legendre transformation of a given motion between the $`(q,t)`$ and $`(q,E)`$ domains is to be clearly distinguished from the following study of variations of the motion in the $`(q,t,E)`$ domain. We now proceed with the latter. We once again drop the microstate suffix $`\mu `$ and the eigenstate suffices $`ij`$, using them only when clarity demands.

### B The energy variational derivative of the action and a trajectory kinematic equation

Consider the energy variational derivative of the (quantum or classical) action $`S(q,t,E)`$ which, following Floyd, we define as the epoch

$$\tau =\frac{\delta S}{\delta E}.$$
(33)
In using the variational derivative $`\delta S/\delta E`$ we are consistent with section IV but differ from Floyd, who uses the partial derivative. The use of the latter is justified only in the case of stationary classical and quantum systems which have a continuum of states of spatially uniform energy. Our use of the variational derivative acknowledges that in quantum mechanics the variation $`\delta E`$ may necessarily be a function of $`q`$ and $`t`$ in order for the Schrödinger equation and its boundary or initial conditions to remain satisfied.

We suspend Floyd’s assumption that $`\tau `$ is a constant and examine its time derivative along a trajectory in $`(q,t)`$ space in the energy plane labelled by $`E`$. Thus, from equation (33),

$$\frac{d\tau }{dt}=\frac{d}{dt}\left(\frac{\delta S}{\delta E}\right)=\left[\frac{}{t}\left(\frac{\delta S}{\delta E}\right)+\left(\frac{\delta S}{\delta E}\right)\dot{q}\right].$$
(34)
Changing the order of $`t`$ and $`\delta E`$ and using the Hamilton-Jacobi equation with $`E=S/t`$ (see equation (14)), equation (34) becomes

$$\left(\frac{\delta S}{\delta E}\right)\frac{\dot{q}}{(1d\tau /dt)}=1$$
(35)
where $`\dot{q}`$ is the total derivative of the trajectory position with respect to time.

A further equation involving $`\delta S/\delta E`$ is obtained from the energy variational derivative of the kinetic energy $`T=(S)^2/2m`$. Thus, after changing the order of $`\delta `$ and $``$ as applied to $`S`$,

$$\left(\frac{\delta S}{\delta E}\right)\frac{S}{m(\delta T/\delta E)}=1,$$
(36)
where the $`ij`$ subscripts, dropped from $`(\delta T/\delta E)_{ij}`$, are understood. Comparison of equations (35) and (36) shows that equation (35) is satisfied by

$$\dot{q}=\frac{S}{m}\frac{(1d\tau /dt)}{(\delta T/\delta E)}.$$
(37)
Whilst this solution is consistent with the one-dimensional case, for which it is the only solution, in higher dimensions it is not the only solution of equation (35). For an $`n`$-dimensional system, equation (35) defines $`\dot{q}/(1d\tau /dt)`$ as a member of an infinite set, $`𝒬_{n1}`$, of vectors passing from the origin and ending on an $`(n1)`$-dimensional plane which has a normal vector $`(\delta S/\delta E)`$ and which is located at a distance $`1/(|\left(\delta S/\delta E\right)|)`$ from the origin. For such a system, equation (36) shows that $`S/(m(\delta T/\delta E))𝒬_{n1}`$ also. By choosing $`\dot{q}/(1d\tau /dt)`$ to be this particular member of $`𝒬_{n1}`$, we obtain equation (37) and achieve consistency with the 1-dimensional case. Thus, in all except the 1-dimensional case, the form of the trajectory equation (37) is suggested, rather than prescribed, by the above variational analysis.

We note that in this derivation, no use has been made of the probability conservation equation (4). Thus, unlike Bohmian trajectories, the trajectories defined by equation (37) do not have a natural probabilistic interpretation. We discuss this matter further in the next section.

## VI Connections between classical, Floydian and Bohmian trajectories

The Floydian and Bohmian trajectories are connected by different assumptions about the time dependence of the parameter $`\tau `$. In classical mechanics, this parameter is demonstrably a constant with respect to the time $`t`$. This is not the case in quantum mechanics.

### A Classical systems

For conservative classical systems, $`\delta T/\delta E=1`$ and in cartesian coordinates $`m\dot{q}=S`$ so equation (37) implies that $`d\tau /dt=0`$ ie $`\tau `$ is indeed a constant without assumption. On the other hand, the converse is not true: Assuming that $`d\tau /dt=0`$ is a necessary though not sufficient condition for $`m\dot{q}=S`$ except in the case of one space dimension, as discussed in section V B.

### B The kinematic equation of Floydian trajectories

The kinematic equation of Floydian trajectories arises from the assumption that $`d\tau /dt=0`$. Thus, with this assumption equation (37) becomes

$$\dot{q}=\frac{S}{m}\frac{1}{(\delta T/\delta E)},$$
(38)
in which we may use the expression for $`(\delta T/\delta E)_{ij}`$ given in equation (31).

If we explicitly write the action in the stationary form $`S=WEt`$ so that $`S=W`$ and recall from equation (14) that $`\delta T/\delta E=1\delta Q/\delta E`$, then equation (38) becomes

$$\dot{q}=\frac{W}{m}\frac{1}{(1\delta Q/\delta E)},$$
(39)
which is precisely the form used by Floyd . For (Floydian) microstate trajectories defined by this kinematic equation, $`\dot{q}\frac{p}{m}`$ where $`p=W`$ is the conjugate momentum. This is different from the kinematic relation for both classical and Bohm trajectories, the latter applying directly to the reduced action $`W`$ of the wave function and not microstate solutions $`W_\mu `$.

In passing, we note that equation (39) represents a deformation of the classical kinematic relation between conjugate momentum and velocity. In making a parallel with relativity, Faraggi and Matone  capture this deformation in the form of a ’quantum mass field’

$$m_Q=m(1\delta Q/\delta E)=m\delta T/\delta E,$$

so restoring the classical form $`\dot{q}=p/m_Q`$. Only for microstate solutions of bound states, for which $`\delta T/\delta E`$ is not uniformly zero, is the ’quantum field mass’ a meaningful concept. The multiplicity of such microstates implies an equal multiplicity of ’quantum mass fields’.

### C The kinematic equation of Bohmian trajectories and quantum time

The kinematic equation of Bohmian trajectories is derived from the conservation of probability equation and the classical form for the current density. This unambiguously implies that $`\dot{q}=\frac{1}{m}S`$ in the coordinate representation so coinciding with the classical kinematic relation. (A different kinematic equation occurs in other representations of the Bohm formulation .) The classical and Bohm kinematic relations therefore coincide. Thus, rather than being derived from an assumption regarding $`\tau `$, the Bohm kinematic equation may be used to define the derivative $`d\tau /dt`$. In the the coordinate representation, equation (37) implies that

$$\frac{d\tau }{dt}=\frac{\delta Q}{\delta E}$$
(40)
showing that, in general, the ’epoch’ $`\tau `$ is not a constant for a Bohmian trajectory. Defining the ’quantum time’ as $`t_Q=t\tau `$, we find that

$$\frac{dt_Q}{dt}=(1\delta Q/\delta E)=\delta T/\delta E,$$
(41)
along the trajectory. This equation thus suggests that the quantum potential acts to deform time. The kinematic equation for Bohm trajectories with respect to $`t_Q`$ is then

$$\frac{dq}{dt_Q}=\frac{S}{m(\delta T/\delta E)}$$
(42)
which, for stationary systems, has the same form as equation (38) for Floydian trajectories with respect to the classical $`t`$.

We observe from equation (41) that for conservative classical systems, $`t_Q`$ and $`t`$ advance at the same rate and, at most, differ by an arbitrary constant shift $`\tau `$. This is consistent with section VI A and the idea that a classical particle has neither internal energy nor extension in space-time.

In the case of unbound quantum systems , there is a single microstate which can be compared with a Bohmian trajectory. In this instance, the above equation (42) suggests that, in a region of interference where the quantum potential is non-zero, the differing forms of the Floydian and Bohmian trajectories can be related by a deformation of time which is determined by the quantum potential.

For stationary quantum systems the square integrable real wave functions (such as in bound states), have $`W=S=0`$ and $`\delta Q/\delta E=1`$, so neither Bohm trajectories nor quantum time $`t_Q`$ evolve with respect to the classical time $`t`$; they stand still. In fact for such systems, equation (42) is not well defined, this being a symptom the fact that the basis of Bohm trajectories is inconsistent with that of the microstates of bound state wave functions.

### D Applications to different potentials

To illuminate the general results of the previous sections, we briefly discuss their application to different forms of potential.

In the case of unbound potentials, results are obtained for both Floydian and Bohmian trajectories and the Floydian $`W_\mu `$ and the Bohmian $`W`$ are the same. For constant potentials classical results are obtained. For the step-potential which is less than the total energy $`E`$ , there are points in the interference domain where the energy derivative of the kinetic energy passes through zero and changes sign so yielding Floydian trajectories that have infinite positive and negative velocities. This is in distinction to the positive Bohmian velocities in the same domain.

In the case of bound potentials  Floydian and Bohmian trajectories differ not only because of differences in the kinematic equations of their trajectories but also because $`W_\mu `$ and $`W`$, from which the are respectively derived, differ both in form and multiplicity. For bound potentials, Floydian trajectories are found, in this work, to have explicit time dependence which is at the beat frequency between adjacent eigenstates whereas Bohmian trajectories have zero velocity. Though not explicitly addressed here, we would also expect this principle of explicit time dependence to apply to finite bound potentials such as the finite square well  where $`E<V_0`$ (the well depth).

## VII Conclusions

The emphasis of this paper has been technical rather than philosophical and interpretational. Having briefly reviewed the Floydian trajectory method in quantum mechanics, we focused attention upon basis of the kinematic equation of Floydian trajectories.

Firstly, motivated by concerns over the definition of the total energy derivative of the kinetic energy and thus the quantum potential, especially in the case of the quantum systems with discrete energy spectra, we derived an explicitly time dependent form for the energy derivative of the kinetic energy for use in the kinematic equation of Floydian trajectories. This explicit periodic time dependence, common to all microstates of a given eigenstate and a modification of the kinematic equation used by Floyd, is at the beat frequency between two adjacent eigenstates of a bound system and is therefore itself quantised. In the limiting case of systems with continuous spectra, this modified kinematic equation becomes the original Floydian form. Our approach is based on the use a total energy variational derivative of the kinetic energy derived from a variation of the wave function constrained to be in the space spanned by the eigenfunctions of the Schrödinger equation. This is in contrast to the approach of Floyd who, treating the total energy as a parameter of each microstate, directly computes the alternative total energy partial derivative of the quantum potential without consideration of the temporal interference with other eigenstates induced by a variation in energy.

Secondly, we have examined the derivation of the kinematic equation of Floydian trajectories from the energy derivative of the action (the epoch) by first setting aside Floyd’s assumption that the epoch is a constant and working without the restriction to one dimension. We find that in general, whilst the magnitude of the velocity of the Floydian trajectory is well defined, its direction is not; the exception being in one-dimension or in systems with uncoupled degrees of freedom. This arbitrariness of direction is overcome by assuming the multi-dimensional form of the kinematic equation to be consistent with that in one dimension. We show that in classical mechanics the epoch is a constant without assumption, whereas assuming it is so in quantum mechanics leads the kinematic equation of Floydian trajectories. On the other hand, assuming the Bohm kinematic equation we find the epoch not to be constant and are led to the idea of the deformation of time by the quantum potential. In the case of unbound systems, for which only is this a consistent interpretation, we find that Floydian and Bohmian trajectories can be related by a time deformation.

Further work is required to investigate the application of the explicitly time dependent form of the total energy derivative of the kinetic energy to the various bounded potentials already analysed by Floyd.

## Acknowledgement

The author gratefully acknowledges the helpful comments provided by M. Matone and by E. R. Floyd in response to an earlier version of this paper. Some points raised by E. R. Floyd are the subject of further investigation.