1. Introduction
In classical physics position and momentum are two conjugate variables that determine motion. It can be said by analogy that amplitude and phase are used in a similar role in quantum mechanics. However, the analogy is not perfect. Pauli [
1] asked whether the wave function can be constructed from the knowledge of a set of amplitudes only. Lamb [
2] argued that from a set of values of the amplitude of the wave function and their rates of change, the wave function including its phase will be found uniquely. Counter examples were later given [
3,
4] and now it turns out that the knowledge of amplitude and certain information about the analytical values of the wave function are required together to construct the quantum states. In the study of Guimaraes and Baseia [
5] situations with defined phases for stationary or moving type of fields were produced.
Rayleigh [
6] showed in the field of classical waves that through the interference between the research wave and the researcher wave, the magnitude and phase of acoustic waves can be determined individually, that is, by finding the minimum and zero values. Mandel and Wolf [
7] noted with reference to the ‘complex analytic signal’ (an electromagnetic field with positive frequency components) that the position of the zeros (from which the phase can be determined) and the intensity represent two groups of information intertwined by the analytical value of the wave.
Interference of optical waves is clearly a phase phenomenon. In classical systems this results from the superposition of positive and negative real wave amplitudes. Phase interference can be used in optical systems to obtain a type of quantum measurement, known as non-destructible measurements [
8] (chapter 19). This is done to make a measurement that does not change any value of the system at the expense of other values that may change. In optics, the phase is the one that can act as a test to determine the intensity (or number of photons). The phase may change during measurement, while the number of photons does not change [
7,
9]. To conclude, we will note the theoretical demonstration, presented in [
10], which shows that any operator (discrete, finite dimensional) can be constructed using only optical means.
The question of determining the phase of a field (classical or quantum, such as a wave function) from a modulus (absolute value, amplitude) of the field along an actual parameter (for which a single experimental determination is possible), is known as the ‘phase problem’ [
7]. The relations derived in the next section represent a formal description for determining the phase given the amplitude, and vice versa.
In 1984 Berry described his discovery of time-independent phase changes in multi-component states known as a geometric/topological phase or Berry phase [
11] based, among others, on Aharonov and Bohm who discovered the topologically acquired phase [
12] named after them. In their discovery, they showed that when an electron moves along a closed path along which the magnetic field is zero, it acquires an observed phase change that is proportional to the “vector potential”. The topological aspect, i.e., that the path is within a multi-linked part of space (in physical terms, the closed path cannot be contracted without encountering a magnetic field), has also been shown to be of great importance [
13,
14], mainly through extensions and applications of the phase-change concept of Aharonov–Bohm [
15] which led to a number of developments in many fields of physics [
16]. The term “open path phase” denotes a fully non-cyclical development [
17,
18]. This term, unlike the value of the Berry phase, is not a measurable constant, but is accessible, in part, by experiments.
As already noted, the Berry phase and the open path phase indicate changes in the phases of the components of the state, rather than the overall phase change of the wave function, which belongs to the “dynamic phase” [
11,
19]. In some cases the presence of more than one component in the state function is a topological effect. This determination is based on Longuet–Higgins’ theorem “Topological Test for a Intersections” [
20], which states that if a given wave function of an electronic state changes sign when it moves around a loop in a nuclear configuration space, then the state must be degenerate with another state at some point in a loop.
To summarize, regarding the effects of the phase in complex states, we will look at the two ways in which it is possible to arrive at a complex description of a phenomenon that takes place in the real world:
First, the time-dependent wave function is necessarily complex due to the form of the Schrödinger equation that is time-dependent, which includes the square root of minus one i. Second, there are also defined functions that do not include the i (such as the Schrödinger equation that is independent of time). Here, too, the wave equation can be complex by having some of the variables obtain complex values. This allows for the removal of possible ambiguities that arise in the solution at a singular point that can be infinite. In addition to this it can often be useful to refer to a number of physical parameters that appear in theory as a complex quantity and that the wave function will include analytical values in relation to them. This formal procedure may even include basic constants such as and so on.
Yahalom and Englman [
21] presented an analytical formulation of a one dimensional scattering process of a microscopic wave packet, which is returned from an infinite potential barrier. It has been found that under conditions suitable for the electron, there is a reciprocal relationship (Kramers–Kronig) between the phase and the amplitude of the wave function of the propagating particle. These interrelationships show the analytical part (phase) uniquely from within the modulus.
The physical basis for this relationship was clarified in [
22,
23], as stemming from the lower delineation of energies. When the analytical conditions are not fully met, for example, due to being zero points in the wave function, the phase can still be calculated from the amplitude, which is given as a function of the composite time with the introduction of the conditions of Toll [
24]. When the underlying analytical conditions are fully met, then they are sufficient to calculate the amplitude as a function of time, i.e., along the actual t-axis.
Following Yahalom and Englman’s research, the present study will seek to expand their analysis and present a reciprocal relationship between the modulus-log and the phase of superposition of the two energy levels of the particle in an infinite potential well at the center of which there is a finite barrier. We remark that the boundary conditions in a well are different, however, than the boundary conditions in the case of an infinite barrier described in [
21].
3. Temporal Kramers–Kronig Relations Theory
The wave function is a complex function and therefore the following relationship is satisfied:
where
is the amplitude and
is the phase. If proper analytical conditions are met the phase and amplitude are related by the temporal Kramers–Kronig relations, which are reminiscent in form with the relations obtained by Kramers and Kronig for the real and imaginary parts of the dielectric function [
25]. However, there are also differences, in particular the classical Kramers–Kronig relations are derived in the complex frequency plane, while our results are derived in the complex time plane. To avoid confusion, we use the term
temporal Kramers–Kronig relations whenever it is appropriate. Another major difference is that we deal with the quantum wave function and not with the electromagnetic dielectric function as was done originally by Kramers and Kronig, hence the physical object under consideration is entirely different. An analysis of
in the lower half of the
t-plane is required for using the Cauchy theorem to link the real and the imaginary parts of
on the real axis. For all complex times
z within a closed contour
C in the lower half of the
t-plane, Cauchy theorem gives:
Provided that is analytical inside the contour.
We assume that the contour
C can be chosen to include the real t-axis and a large infinite semicircle in the lower half of the plane. In case the logarithm of the wave function vanishes on the circle half at infinity, the Cauchy integral can be written as follows:
Taking the limit as a complex time approaching the real axis from below, we write
in (5):
Taking the limit
we may contain the point inside the trajectory by drawing a tiny semicircle over point
t (see
Figure 1). The denominator can be formally written as:
where
P is the principal value. The delta function is used to describe the contribution from the small semicircle that goes in against the clock direction over the pole at
. Using (5) and simply rearranging turns (6) into:
The real and imaginary parts of this equation are:
The temporal Kramers–Kronig relations expresses the relation between the real part and the imaginary part of the logarithm of the wave function. As stated, the relation expressed in Equation (
9) is true only in the case where the function is log-analytic in the lower half of the plane. This means that the function will have neither singular points nor zeros in the lower half of the plane.
4. Checking Zeros and Singular Points in the Lower Half of the Complex Time
Our wave function (the general solution for the Schrödinger equation that is time-independent) is given in formula (1) where
and also holds:
We want to test analyticity in the lower half of the complex time
so we will write:
We will write the function (1) as follows:
When
, it is clear that the function has no singular points in the lower half of the plane and therefore we will concentrate on identifying the zeros if any. We will define:
We will calculate a limit of
when the imaginary part of time approaches to minus infinity:
We will now define an asymptotic wave function:
We need a log analytic function, i.e., a function that does not accept values of zero or infinity in the lower half of the complex time plane. Moreover, it is required that the function approaches unity on the infinity circle so that the logarithm vanishes on this circle. Thus we will use the asymptotic function. We have not yet proved that this function is log analytical and this can only be proved with respect to individual cases where the temporal Kramers–Kronig relations exists. We will first deal with the superposition of two eigenfunction
When
. In order for the wave function to be log analytical so that
is not diverging at any point, we must find the condition under which there are no zeros in the complex plane of time. This is a prerequisite for the temporal Kramers–Kronig relations to be satisfied. That is, since
then we require that
. Assuming that the placement coefficients are equal
the condition for the absence of zeros is:
Therefore if condition (17) is met at each and every point in the lower part of the complex plane.
6. Combining a Double Potential Well Model
We consider a potential well as depicted in
Figure 2. A particle is moving to the right from the left hand side of a finite-sized barrier located in the center of the well at
. The wave function of the particle is represented as the superposition of the two energy levels of the particle in a double well. Quantum mechanics predicts that the particle will be returned after encountering the barrier with high probability. However there is also a probability of tunneling through the barrier and towards the right side of the potential well as shown in
Figure 2:
The eigenfunctions of the double well are given in Equation (
31):
When
and
. The expression for energy as a function of the wave-number
k is given in the formula
and will be related to
q according to equation:
Following Equation (
30) and using the chain rule we obtain the following inequality as a condition for the Kramers Kronig relations to hold:
Since the derivative
is positive it will not affect the derivative of the wave function and therefore we obtain the condition:
In order to establish the domain without zero points we will use Equation (
35). In the domain left to the barrier we have according to Equation (
31):
Since
x is positive then we are left with the cosine function which is negative provided that:
, we will note that it is not possible to reach the upper limit and therefore it is irrelevant. The momentum is given by the relation:
Therefore the result is:
or
If those condition are met, the function is a log analytical and satisfies the temporal Kramers–Kronig relations, meaning that the amplitude can be calculated using the phase and the phase can be calculated by using the amplitude. When the condition is not met it is not possible to calculate the phase from the amplitude and the amplitude from the phase, so we have uncertainty about the second magnitude given the first magnitude. For this reason the above inequality, is denoted the principle of log analytical certainty. Thus, it can also be said that there is uncertainty in the relationship between the phase and amplitude in the field in which the following condition is met:
or
According to inequalities (39a) and (39b) the larger the momentum (or energy) of the particle, the smaller the x area where the uncertainty between the phase and the amplitude takes place, and the greater the area where certainty exists. Things can also be presented differently if we look at the distance
d between the barrier and the measuring point (see
Figure 3):
This condition means that as
p increases (i.e., as energies increase) so, the area of certainty adjacent to the barrier increases and then the area of uncertainty decreases. It is also seen that as the width of the well increases,
d (the range of certainty) increases. We will multiply the two sides of inequality (41) by
and use Equation (
37) to obtain:
The
product is constant under single potential well conditions (without a barrier), whereas in the case of a double well it increases moderately with the increase in
a (for example, with doubling
a the product increases by 4.55%, and by taking four time of
a times the product increases by 6.95%). The above inequality is reminiscent in its form of Heisenberg’s principle of uncertainty
. It can be seen, however, that Heisenberg’s inequality does not define a range in which the certainty or uncertainty exists, but rather defines a minimum magnitude to which the product of standard deviations of momentum and position of the particle wave function is equal. So as the standard deviation of the momentum increases so we will talk about a smaller initial standard deviation of the position, whereas in the present case inequality describes the condition for being in an area where there is an uncertainty in the relationship between phase and amplitude. We see that near the barrier we have certainty in the sense that the amplitude and phase can be reconstructed from each other while at a distance exceeding the value derived the phase can not be reconstructed from the amplitude and vice versa. In the domain within the barrier (the second domain) it can be seen according to
Figure 4 that there is no uncertainty since the derivative will always be negative. That is, at each point there is a phase-amplitude relationship.
7. Areas of Certainty and Uncertainty beyond the Barrier
The expression for the energy eigenfunction in the third domain (the domain to the right of the barrier) is:
We will define
and simplify the function as follows:
where:
Taking a derivative of the function in Equation (
44) by Equation (
35):
In Equations (48) and (49) we will neglect the terms
and
. This can be justified on the basis that the domain of certainty is close the barrier. It is obvious that for a very small
(very close the barrier) the negligence of
and
is justified. For slightly larger
we notice that the first term in Equation (
48) which is
if we assume that
a and
b are the same order of magnitude. Now, if
and
are the same order of magnitude it is clear that if
is smaller than
a (or
b) even if not extremely small the terms
and
can be neglected in comparison to
and
, respectively. We examined a number of concrete cases and saw that the neglection is indeed justified. As the log analytic condition requires that:
Dividing by a positive cosine and a positive
the following inequality is obtained as a condition for log analyticity:
Note that the cosine is negative when
. It can be seen that the log analytical principle does not hold and this is because inequality is reversed. Thus:
Therefore, the principle of certainty in the third domain of the well is:
The condition for uncertainty involves reversing the sign of inequality:
It can be said that these inequalities are similar to Heisenberg’s uncertainty principle. From this it can be deduced that in the third domain that provided that the higher the eigenenergies the larger the value on the inequality right side. It will also follow that under the same conditions the area of certainty increases and therefore the area of uncertainty decreases.
A concrete example is given in
Figure 4. Where we illustrate that the areas of certainty and uncertainty as derived from our numerical calculations. The following parameters of the double well are assumed: Depth
eV, width of each of the two wells
Å, and barrier width
Å.
9. Testing Phase-Amplitude Relations
We will now examine whether the temporal Kramers–Kronig relations hold. We will look at the first domain of the well and consider the uncertainty relation described in the Equation (39). Let us take a point closest to the barrier (because this is the domain where the phase-amplitude relations should be satisfied), for example,
Å. We will use Equation (
56) and write the principal value explicitly:
is selected. We compare a graph for the amplitude logarithm derived by a standard straight forward calculation to the temporal Kramers–Kronig integral of the phase we obtain that the graphs are the same as can easily be verified in
Figure 6. The graph of the amplitude for a longer period of time is shown in
Figure 7.
Now we will do an Inverse transform to find the phase according to Equation (
57) and use the principal value as we did in Equation (
58), then we will draw a phase graph as a function of time for straight forward phase calculation and for the temporal Kramers–Kronig calculated phase and obtain the graphs for both cases shown in
Figure 8.
Now we study the barrier, we know that the barrier domain is certain, that is, there are no zeros at all in the complex time plane (see
Figure 9). We will use the transform (56) and draw graphs in
Figure 10.
Now we turn to Equation (
57) which is an inverse transform and draw a graph for the phase shown in
Figure 11.
Now we will turn to the third domain of the well in which the uncertainty in described by Equation (
55). In this domain we choose a point close to the barrier that satisfies the conditions of certainty, for example,
Å. We will refer to Equation (
56) and draw graphs of the amplitude and its logarithm in
Figure 12 and
Figure 13.
Now we turn to Equation (
57) and draw a graph depicted in
Figure 14.
From the graphs we have obtained it can be seen that indeed the relations between phase and amplitude exists according to Equations (56) and (57) but only at points where there are no zeros in the complex time plane. For completeness, we will make a temporal Kramers–Kronig evaluation at points that do not meet the conditions of certainty (points where there is no phase-amplitude relations). For example, we take
Å and draw a graph of the amplitude logarithm as a function of time depicted in
Figure 15.
A graph of the phase as a function of time in
Figure 16.
We plot the temporal dependence of the amplitude of the wave function for
Å and observe that the amplitude has zeros on the time axis and hence there is no phase-amplitude relations according to
Figure 17.
10. General Properties of a Shannon’s Entropy in the Case of Tunneling
We shall now discuss briefly the concept of the phase of a wave function and its relationship to Shannon’s entropy associated with the wave function. Let us investigate the time evolution of the entropy [
26]. We shall find it beneficial to describe the evolution of a wave function
using the amplitude-phase description:
The amplitude and phase can be used to define Madelung’s density and velocity field as follows:
Those quantities satisfy the continuity equation:
and the Euler equation:
in the above
Q is the quantum potential. Those real equations are mathematically equivalent to Schrödinger’s complex equation for a single quantum particle. In terms of the probability density the entropy takes the simple form:
Taking a temporal derivative of the above quantity we obtain:
Using the continuity Equation (
61) results in:
The above expression can be integrated by parts and using Gauss theorem we obtain:
If the surface
is chosen such that there is no probability flux on this surface we are left with:
According to the theory of random variables the right hand side is an expectation value [
27], hence:
A correlation of two random variable
A and
B is defined as [
27]:
We may generalize the above to vector quantities
and
such that:
Thus the entropy production rate is essentially the correlation between the gradient of phase and log amplitude of the wave function. We also notice that for independent random variables [
27]:
Hence,
if one may treat
as independent of
we may write:
Now if in addition the phase is chaotic such that
the entropy production will be rather low. However, if on the other hand the temporal Kramers–Kronig relations hold those quantities are strongly correlated. It follows that one may write the entropy production rate equation using Equation (
9) in the form:
Now, according to the theory of random processes, an autocorrelation of a random process
is given by [
27] a correlation of two time samples of the process:
Thus the entropy production rate is given in terms of the phase autocorrelation as:
We recall that autocorrelation even of white noise is not null and gives a delta function expression in time.
Now the connection between entropy increase and tunneling discussed in [
26] can be also understood in the framework of the temporal Kramers–Kronig relations which are prominent in and in the vicinity of the barrier and imply correlation between phase and amplitude (as described in the above heuristic discussion). We also suggested in [
26] (see also [
28]) that entropy increase may be connected to the emergence of quantum chaos.