Discovering Bohr’s Yin-Yang Diagram in Quantum Tunneling Dynamics

: On 17 October 1947, Niels Bohr was made a knight of the Order of the Elephant by the King of Denmark in view of his outstanding achievements and contributions to science. Bohr designed his own coat of arms that featured a pattern of Yin and Yang (Tai Chi symbol) to symbolize the wave–particle complementarity. However, Bohr’s Yin-Yang diagram (YYD) was neither drawn based on the principles of quantum mechanics, nor did it originate from the traditional Taoist YYD. Scientists still have doubts about the legitimacy of using YYD as the icon of the wave–particle complementarity, because the YYD belonging to quantum mechanics itself is unknown so far. This paper reports the YYDs existing in quantum mechanics and justifies the role of YYD in the wave–particle duality by showing that any system, whether classical or quantum, has an ideal YYD as long as it satisfies Bohr’s principle of complementarity (BPC). The deviation of a deformed YYD from the ideal YYD indicates the extent to which a real system satisfies BPC. This paper constructs the quantum YYD by the complex quantum trajectory of a particle tunneling via a step barrier, which displays the continuous transition between the wave behavior and the particle behavior. It appears that the YYD designed by Bohr in his coat of arms resembles the YYD generated by tunneling motion, not only in appearance but also in the governing equation.


Introduction
Before Niels Bohr, the principle of complementarity had appeared in many intellectual disciplines, dating as far back as Yin-Yang thought in ancient China more than 3000 years ago.Yin and Yang elements appear in pairs, called the complementary pairs, such as the moon and the sun, female and male, dark and bright, cold and hot, passive and active, and so on.Complementary elements are opposite in nature yet rely on each other to exist.Bohr was the first person to introduce the principle of complementarity into the world of physics and regarded it as a universal principle.Bohr's principle of complementarity (BPC) holds that objects have certain pairs of complementary properties which cannot all be observed or measured simultaneously.Using one particular piece of apparatus, only one of the features could be made manifest at the expense of the other.
Taoist Yin-Yang diagram (YYD), which is also known as Tai Chi symbol in China [1,2], is a graphical demonstration of the complementary Yin-Yang pair by means of areas in black and red to display their interlaced changes within a period (see Figure 1a).Why Bohr used this ancient symbol to represent his principle of complementarity [3] has remained a mystery for nearly eighty years.Except for philosophical or even religious connections, physicists have so far found no scientific clues to substantiate Bohr's original ideas [4].
The key reason why this mystery has not been solved for many years is because physicists and even Bohr himself have no idea how to create the YYD from the principle of complementarity.Taoist YYD is a graphic representation of the ancient Chinese version of the principle of complementarity [5].However, Taoist YYD has no unified form, and principle of complementarity.Taoist YYD is a graphic representation of the ancient Chinese version of the principle of complementarity [5].However, Taoist YYD has no unified form, and Chinese scholars of various dynasties have proposed different methods for the drawing of Taoist YYD.Different complementary pairs will derive different forms of YYD, so there is indeed no YYD that can uniquely represent BPC.Although the forms of YYD are quite diverse, these different forms of YYD have common characteristics.This paper proposes the common rules in the construction of various YYDs.The ancient Taoist YYD and the quantum YYDs of wave-particle duality and tunneling effect can all be constructed by the proposed common rules.The resulting YYDs are all concrete realizations of the idea of complementarity, and any one of them can be used as a symbol of BPC.As shown in Figure 1, the areas occupied by red and black in Taoist YYD serve as the quantitative indices  and  for the state A (Yang) and state B (Yin), respectively.Using the two indices  and  , the two properties of BPC, i.e., joint completeness and mutual exclusion, can be expressed quantitatively.
YYD is a graphical indicator that reflects the degree to which a system satisfies the principle of complementarity, and it is present in almost every system in nature, whether classical or quantum.What Bohr could not have imagined was that he did not need to borrow Taoist YYD to convey his principle of complementarity, because the principle of complementarity is the theoretical foundation of Taoist YYD.This study reverses Bohr's idea by applying the definition of YYD to reproduce Taoist YYD, which is an ideal YYD satisfying BPC, and using the same definition to construct the YYDs belonging to quantum mechanics itself.The results obtained here reveal that the YYD designed by Bohr in his coat of arms is different from Taoist YYD but is close to the quantum YYD based on tunneling effect.As shown in Figure 1, the areas occupied by red and black in Taoist YYD serve as the quantitative indices p A and p B for the state A (Yang) and state B (Yin), respectively.Using the two indices p A and p B , the two properties of BPC, i.e., joint completeness and mutual exclusion, can be expressed quantitatively.
YYD is a graphical indicator that reflects the degree to which a system satisfies the principle of complementarity, and it is present in almost every system in nature, whether classical or quantum.What Bohr could not have imagined was that he did not need to borrow Taoist YYD to convey his principle of complementarity, because the principle of complementarity is the theoretical foundation of Taoist YYD.This study reverses Bohr's idea by applying the definition of YYD to reproduce Taoist YYD, which is an ideal YYD satisfying BPC, and using the same definition to construct the YYDs belonging to quantum mechanics itself.The results obtained here reveal that the YYD designed by Bohr in his coat of arms is different from Taoist YYD but is close to the quantum YYD based on tunneling effect.
Taoist YYD contains two states: Yin and Yang, whose magnitudes wax and wane over time, but their sum remains constant, hence forming a pattern like two fish rotating around each other within a circle.The Yin-Yang states in Taoist YYD are a complementary pair possessing three properties.The first property is that they are jointly complete, that is, the two states constitute the entire system.The second property is that they are mutually exclusive, that is, there is no intersection or interference between the two states.The combination of the above two properties indicates that the Yin-Yang states satisfy BPC [7], which can be expressed as p A + p B = 1 in terms of the quantitative measures p A and p B (see Figure 1b).The third property of the Yin-Yang pair is that the pair remains complementary over a complete period, i.e., p A (t) + p B (t) = 1, 0 ≤ t ≤ T. The YYD composed of the Yin-Yang states having the above three properties is said to be an ideal YYD.The traditional Taoist YYD is a typical example of ideal YYDs, and this paper reveals two ideal YYDs emerging from quantum mechanics, one from two-path interference experiments and the other from quantum tunneling dynamics.
The complementary pair (A, B) that Bohr considered includes: (i) position and mo- mentum, (ii) spin on different axes, (iii) wave nature and particle nature of a quanton, (iv) value of a field and its change, (v) entanglement and coherence, and (vi) photon polarization.Among these complementary pairs, the most discussed is the particle nature and the wave nature of a quanton, whose relationship with BPC has been tested by many experiments in recent decades.In these experiments, the quantitative measures for the complementary pair are given by (p A , p B ) = D 2 , V 2 , where D is the distinguishability of the which-way information, and V is visibility of wave fringe measured by the Mach-Zehnder interferometer (MZI) [8,9].The manifested wave and/or particle property of the tested quantons depend on the detecting devices used in the experiments.In the all-or-nothing cases [10][11][12], quantons behave as a particle (p A = 1 and p B = 0) or as a wave (p A = 0 and p B = 1), while in the intermediate cases [13][14][15][16], quantons behave simultaneously as a particle and as a wave, whose quantitative measures satisfy the inequality p A + p B ≤ 1.It was further verified [17] that the equality p A + p B = 1 holds if the tested quantons have no internal degree of freedom, for which the conditions of joint completeness and mutual exclusion are both met.
In recent years, the advancement of quantum beam splitters [18] has allowed the MZI to be set in a state of quantum superposition of "on" and "off" such that the wave and particle behaviors of a single photon can coexist simultaneously, with a continuous morphing between them.Unlike the classical MZI, where the wave nature and the particle nature are treated independently, the wave and particle natures of a single quanton in the setup of quantum MZI become quantum states [19,20] so that there can be interference between them.The interference between the wave state and the particle state in the quantum MZI implies that they are not mutually exclusive and cannot form a qualified complementary pair (A, B), for which the violation of the complementarity, 1 < p A + p B ≤ 2, has been experimentally observed [21,22].
Although BPC has been tested by numerous experiments, these experiments still cannot tell us why Bohr used the YYD to symbolize his principle.As a first point, this paper constructs the ideal YYD and the deformed YYD from the experimental data (p A , p B ) = D 2 , V 2 .According to the difference between the ideal YYD and the deformed YYD, one can judge the extent to which the experimental results satisfy the principle of complementarity.Interestingly, the ideal YYD constructed by the two-path interference experiment is closer to the traditional Taoist YYD but different from the YYD designed by Bohr in his coat of arms.This finding prompts us to explore another complementary pair in order to more closer represent Bohr's YYD.
The finding of a more suitable complementary pair is inspired by the quantum tunneling motion of a quanton via a step barrier, wherein one finds that the changes in the quanton's real velocity .
x R and imaginary velocity .
x I synchronously reflect the alternating changes in its particle nature and wave nature.It was pointed out that the wave-particle duality is a direct result of quantum motion occurring in complex space [23][24][25].The validity of the complex-trajectory interpretation of wave-particle duality is further supported by the recent findings of weak measurements, which reveal that every physical observable has a real part as well as an imaginary part [26,27].Similar to the entanglement, the imaginarity has been recognized as an indispensable resource responsible for quantum advantages [28][29][30], which has ruled out the real-valued standard formalism of quantum theory [31].When a quanton exhibits particle behavior, the real component x R alone is sufficient to describe its motion, as treated in classical mechanics.But when a quanton's motion exhibits wave behavior, the imaginary component x I is also needed to describe its motion.The results obtained reveal that the quantum YYD constructed by the comple- x 2 I is closer to Bohr's YYD than the YYD constructed by the complementary pair (p In what follows, first the general YYD and the ideal YYD are defined in Section 2 and it is shown that the traditional Taoist YYD is the oldest example of the ideal YYD satisfying BPC.With (p A , p B ) = D 2 , V 2 , Section 3 constructs the quantum YYD by the data of the two-path interference experiments, from which the similarity between Taoist YYD and the interference-based YYD is shown.Section 4 introduces the mapping between the complex tunneling trajectory and the quantum YYD for a particle's tunneling motion via a step barrier.The particle's tunneling trajectories are solved by quantum Hamilton mechanics [32] in Section 5 to yield the particle's complex position x R + ix I and complex velocity .
x R + i .
x I by which the tunneling-based YYD is then generated.The obtained tunneling-based YYD is a graphical recorder of the particle's tunneling dynamics, which moves together with the particle, recording the velocity of the particle, and presenting it graphically, as demonstrated in Section 6.The tunneling-based YYD is not always an ideal YYD, and how its shape evolves with the intensity of the tunneling effect is discussed in Section 7. Section 8 compares the four different YYDs considered in this paper and clarifies their relationship with BPC, which brings to the following conclusions: complementarity is a universal principle, as Bohr said, and the YYD constructed in this paper is a universal image found all over in the nature.

Definitions for General YYD and Ideal YYD
Traditional Taoist YYD is an ideal YYD satisfying BPC, which is a special case of the general YYD.In this Section, the general YYD is defined first and then it is explained under what conditions the ideal YYD can be induced.A general YYD is defined for any periodic system, which is composed of state A (Yang) and state B (Yin) with quantitative measures p A and p B normalized between 0 and 1.The growth and decline of the Yin-Yang states are conveyed through the interlacing changes of p A (in red) and p B (in black).As shown in Figure 2, the shape of YYD is uniquely determined by its inner and outer curves.
The radial coordinates of the inner and outer curves are called inner and outer radii, respectively.In terms of the quantitative measures p A (t) and p B (t) recorded at each moment t, the inner radius of a general YYD is defined by and its outer radius is defined by r = p A (θ) + p B (θ), where the azimuth angle θ = 2πt/T is measured clockwise from the vertical axis and T is the period of time, as shown in Figure 2. It is noted that the periodic property of the system implies p A (0) = p A (2π) and p B (0) = p B (2π).
The ideal YYD is a special case of the general YYD, which further meets the continuity condition and the complementarity condition.The continuity condition requires that the two sections of the inner curve defined by Equation (1) must be continuous at the origin, i.e., p A (π) = p B (2π) = 0, or equivalently, the state A must occupy the whole system at θ = 2π (θ = 0), and the state B must occupy the whole system at θ = π.The complementarity condition requires that the states A and B must satisfy BPC: p A (θ) + p B (θ) = 1 over the entire interval, or equivalently, the outer curve has to be a unit circle, r = p A (θ) The radial coordinates of the inner and outer curves are called inner and outer radii, respectively.In terms of the quantitative measures  () and  () recorded at each moment , the inner radius of a general YYD is defined by and its outer radius is defined by  =  () +  (), where the azimuth angle  = 2/ is measured clockwise from the vertical axis and  is the period of time, as shown in Figure 2. It is noted that the periodic property of the system implies  (0) =  (2) and  (0) =  (2).
The ideal YYD is a special case of the general YYD, which further meets the continuity condition and the complementarity condition.The continuity condition requires that the two sections of the inner curve defined by Equation (1) must be continuous at the origin, i.e.,  () =  (2) = 0, or equivalently, the state A must occupy the whole system at  = 2 ( = 0), and the state B must occupy the whole system at  = .The complementarity condition requires that the states A and B must satisfy BPC:  () +  () = 1 over the entire interval, or equivalently, the outer curve has to be a unit circle,  =  () +  () = 1, 0 ≤  ≤ 2.
In ancient China, the most representative complementary pair (, ) was day and night.In this case,  represents the proportion of daytime in a day, and  represents (A, B) a qualified complementary pair of a periodic system satisfying the condition p A (t k ) + p B (t k ) = 1 at every moment t k within the period 0 ≤ t k ≤ T. The inner curve of Taoist YYD is defined by the polar coordinates (r k , θ k ), where r k = p A (t k ) for 0 ≤ t k ≤ T/2, and r k = p B (t k ), for T/2 ≤ t k ≤ T, and the outer curve is defined by r k = p A (t k ) + p B (t k ).The quantitative measures p A (t k ) and p B (t k ) used in the traditional Taoist YYD are the proportions of the states A and B in the entire system at time t k , just like the proportions of day and night on a certain day t k of the year.
In ancient China, the most representative complementary pair (A, B) was day and night.In this case, p A represents the proportion of daytime in a day, and p B represents the proportion of night.These two proportions vary on different days of the year, but they always meet the complementarity condition p A + p B = 1.By recording the proportions of day and night for each day of the year in the polar coordinates, with white area denoting the daytime and black area denoting the nighttime, the original YYD in ancient China was then formed [5].However, the original YYD determined by sunshine time is not an ideal YYD, because it does not necessarily satisfy the continuity condition.The continuity condition requires that in the year there be exactly one day with 24 h daytime (at the summer solstice, θ = 0) and exactly one another day with 24 h nighttime (at the winter solstice, θ = π).Except for those who live at a certain latitude on the earth, people cannot actually draw an ideal YYD based on local sunshine hours.
The first ideal YYD, called Taoist YYD, was constructed by Chen Tuan (871AC-989AC), a philosopher in the early Song Dynasty of China.He drew the YYD according to the sequence of the sixty-four hexagrams appeared in the Book of Changes [10] (also called Yijing or I Ching), which is the oldest Chinese classics formed around the ninth century BC.Here, the Taoist YYD is reproduced by the BPC. Figure 2 shows a schematic diagram of constructing Taoist YYD from the measured p A (t) and p B (t) under the situation that the system is divided into four parts, and one period is divided into eight intervals.In Taoist YYD, the quantitative measures p A (t k ) and p B (t k ), k = 0, 1, • • • , 8, represent the Physics 2024, 6 969 proportions of the states A and B in the entire system at time t k .For example, at time t 1 , state A occupies three-quarters of the system, and state B occupies one-quarter of the system, so one gets p A = 3/4 and p B = 1/4.The situation considered in Figure 2 can be extended to the general case, where the system has N/2 partitions and the period has N partitions.Now, the measure p A (p B ) represents the ratio of the number of partitions containing A (B) state to the total N/2 partitions of the system.Therefore, In order to satisfy the complementarity condition, the measure of B is set to p B (t k ) = 1 − p A (t k ).Substituting the p A (t k ) and p B (t k ) defined by Equation (2) into Equation (1) yields the polar coordinates (r k , θ k ) describing the inner curve of Taoist YYD: where the index k and the azimuth angle θ k have the linear relation k/(N/2) = θ k /π; by noting that when the integer k increases sequentially from 0 to N/2, the azimuth angle θ k increases linearly from 0 to π.It is noted that the continuity condition p A (π) = p B (2π) = 0 is also satisfied by p A and p B defined by Equation (2).When N is set to 8, the eight points (r k , θ k ) defined by Equation (3) are just those shown in Figure 2. When N is set to 64, the curve formed by connecting the sixty four coordinate points (r k , θ k ) defined in Equation (3)  recovers the traditional Taoist YYD, which was first constructed according to the order of the sixty four hexagrams in the Book of Changes.The ancient Taoist YYD is reproduced here by the principle of complementarity, i.e., by requiring that the complementary pair (A, B) must be jointly complete and mutually exclusive.As the number of partitions N approaches infinity, the discrete points defined by Equation (3) become a continuous polar curve This is the inner curve of the ideal Taoist YYD formed by two sections of Archimedean spiral, which has a general expression r = a + bθ.If Bohr's original intention was to use Taoist YYD to convey his principle of complementarity, his correct choice should be the pattern in Figure 2, which is Taoist YYD drawn exactly according to his principle of complementarity.However, the YYD designed by Bohr in his coat of arms, as shown in Figure 3, is based on a popular drawing method, which is composed of a unit circle and two inscribed circles with polar equation given by Apparently, the inner curve of Bohr's YYD only passes through the first and third quadrants and is significantly different from Taoist YYD defined by Equation (4).In what follows, the definition (1) of the general YYD to derive the YYDs belonging to quantum mechanics itself.
Figure 3, is based on a popular drawing method, which is composed of a unit circle and two inscribed circles with polar equation given by Apparently, the inner curve of Bohr's YYD only passes through the first and third quadrants and is significantly different from Taoist YYD defined by Equation (4).In what follows, the definition (1) of the general YYD to derive the YYDs belonging to quantum mechanics itself.

Constructing Quantum YYD by Two-Path Interference Experiments
Among various natural phenomena, wave-particle duality best embodies the complementary relationship.In the past several decades, many experiments have been proposed to test whether the wave-particle duality satisfies BPC.However, all experiments on wave-particle duality so far still cannot explain the relationship between Bohr's YYD and wave-particle duality.In this Section, it is shown that regardless of the test results of BPC, a YYD can always be constructed according to Equation (1) by the experimental data ( ,  ) = ( ,  ) from the MZI, where  =  is a measure of the particle nature regarding the distinguishability of the which-way information, and  =  is a measure of the wave nature regarding the visibility of wave fringe.
An experimental setup to construct the YYD for wave-particle duality is shown in Figure 4a [16,33], wherein a single photon pulse is split by a variable beam splitter (VBS) with adjustable reflectivity in the range 0 ≤  ≤ 1, and a piezoelectric transducer (PZT) is employed to adjust the relative phase  between the two paths.The two paths of photon are recombined at the beam merger (BM), which is a fixed equal-probability beam splitter, and finally the photon is detected by the detectors  and  .The measured values of  and  depend on the reflectivity  of the VBS.If  = 0 or  = 1, the which-path information is known deterministically, and the photon manifests the property of a pure particle with  = 1 and  = 0.If 0 <  < 1, the interference fringe appears behind the beam merger and yields a nonzero wave visibility , which achieves its maximum value  = 1 at  = 1/2.

Constructing Quantum YYD by Two-Path Interference Experiments
Among various natural phenomena, wave-particle duality best embodies the complementary relationship.In the past several decades, many experiments have been proposed to test whether the wave-particle duality satisfies BPC.However, all experiments on wave-particle duality so far still cannot explain the relationship between Bohr's YYD and wave-particle duality.In this Section, it is shown that regardless of the test results of BPC, a YYD can always be constructed according to Equation ( 1) by the experimental data (p A , p B ) = D 2 , V 2 from the MZI, where p A = D 2 is a measure of the particle nature regarding the distinguishability of the which-way information, and p B = V 2 is a measure of the wave nature regarding the visibility of wave fringe.
An experimental setup to construct the YYD for wave-particle duality is shown in Figure 4a [16,33], wherein a single photon pulse is split by a variable beam splitter (VBS) with adjustable reflectivity in the range 0 ≤ R 1 ≤ 1, and a piezoelectric transducer (PZT) is employed to adjust the relative phase ϕ between the two paths.The two paths of photon are recombined at the beam merger (BM), which is a fixed equal-probability beam splitter, and finally the photon is detected by the detectors D 1 and D 2 .The measured values of D and V depend on the reflectivity R 1 of the VBS.If R 1 = 0 or R 1 = 1, the which-path information is known deterministically, and the photon manifests the property of a pure particle with D = 1 and V = 0.If 0 < R 1 < 1, the interference fringe appears behind the beam merger and yields a nonzero wave visibility V, which achieves its maximum value Following Berthold-Georg Englert [9], the distinguishability D has two different notions: a-priori distinguishability and a-posteriori distinguishability, where the former is also called predictability P, which refers to a which-way information obtained using an unbalanced interferometer with different particle flux along the two paths.The a-priori distinguishability and visibility are especially useful in the construction of the ideal YYD, because it requires analytical expressions for the distinguishability and the visibility.In the case of a MZI without path loss, the a-priori distinguishability, P, and the a-priori visibility, V, can be expressed analytically as a function of R 1 [33]: Physics 2024, 6 971 Substituting p A = P 2 and p B = V 2 into Equation ( 1), the polar equation of the YYD corresponding to the two-path interference experiment reads where the azimuth angle, θ, is related to the reflectivity, R 1 , as θ = 2πR 1 .The continuity condition p A (π) = p B (2π) = 0 and the complementarity condition p A (θ) + p B (θ) = P 2 + V 2 = 1 are satisfied so that the resulting YYD is an ideal YYD as shown in Figure 4b.
Physics 2024, 6, FOR PEER REVIEW 8 . Following Berthold-Georg Englert [9], the distinguishability  has two different notions: a-priori distinguishability and a-posteriori distinguishability, where the former is also called predictability , which refers to a which-way information obtained using an unbalanced interferometer with different particle flux along the two paths.The a-priori distinguishability and visibility are especially useful in the construction of the ideal YYD, because it requires analytical expressions for the distinguishability and the visibility.In the case of a MZI without path loss, the a-priori distinguishability, , and the a-priori visibility, , can be expressed analytically as a function of  [33]: Substituting  =  and  =  into Equation ( 1), the polar equation of the YYD corresponding to the two-path interference experiment reads There are visible differences in the mathematical form between Taoist YYD (4) and the interference-based YYD (7) in that the former is a linear function of θ/π, while the latter is a quadratic function of θ/π.The quadratic form of the interference-based YYD comes from the feature that quantum probability is the squared magnitude of the wave function.Despite of their difference, both Taoist YYD and the interference-based YYD satisfy the condition of BPC: p A (θ) + p B (θ) = 1, 0 ≤ θ ≤ 2π, so both can be regarded as representative symbols of BPC.
The ideal YYD shown in Figure 4b is constructed from a perfect MZI without path loss; however, path loss in a MZI is inevitable, which causes the actual YYD to deviate from the ideal YYD.In what follows in this Section, the two kinds of nonideal YYDs are considered, one from the violation of the continuity condition and the other from the violation of the complementarity condition.When path loss occurs, the measured values of V and P are different from their ideal values (6) and read [33] respectively, where L 1 and L 2 with values between 0 and 1 are the losses of the two paths and are simulated in the experiment by the two beam splitters BS1 and BS2 with reflectivity equal to L 1 and L 2 , as shown in Figure 4a.Again by substituting p A = P 2 and p B = V 2 from Equation (8) into Equation ( 1), one obtains the deformed YYD caused by the path loss as shown in Figure 4c,d.It can be seen that this type of path loss does not produce the deformation of the outer curve of the YYD, that is, the complementarity condition p A + p B = P 2 + V 2 = 1 is still satisfied, but the inner curve of the YYD is not continuous.The continuity condition requires p A (π) = p B (2π) = 0, but for this case, The discrepancy between p A (π) and p B (2π) denoted by δp A increases with the path loss L 2 as shown in Figure 4c,d.
Figure 4e,f demonstrates the YYD deformation caused by the path loss under another configuration of the experiment, where the setups in the input and output ports are interchanged such that the input port is a fixed equal-probability beam splitter, while the output port is a variable beam splitter with adjustable reflectivity in the range 0 ≤ R 2 ≤ 1.The predictability P and the visibility V for this configuration can be expressed by [33] Opposite to the case (8), as can be seen from Figure 4e,f, the type (9) of path loss causes the outer curve of the YYD to deviate from the unit circle, i.e., the complementarity condition p A + p B = 1 is not satisfied over the entire interval, but the continuity conditions p A (π) = P 2 (1/2) = 0 and p B (2π) = V 2 (1) = 0 still hold.As the path loss increases, the deformation of YYD intensifies and deviates further from the ideal YYD.In short, both the violation of the continuity condition (Figure 4c,d) and the violation of complementarity condition (Figure 4e,f) can deform the ideal YYD (Figure 4b), and the deformation tells us graphically the extent to which BPC is violated.

The Mapping between Quantum YYD and Complex Tunneling Velocity
The YYD generated by the two-path interference experiments shows a static picture of YYD, as found in Section 3.This Section considers a complementary pair (A, B), whose quantitative measures (p A (t), p B (t)) are time dependent so that the pattern of YYD changes with time.This dynamic picture of the YYD can be visualized by the quantum tunneling process, where the particle's real and imaginary velocity components .
x I (t) constitute a time-dependent complementary pair.The concepts of trajectory and velocity given to a tunneling particle are derived from the complex trajectory interpretation of quantum mechanics [32,[34][35][36][37][38][39][40], which is different from the probability interpretation, as compared in Figure 5, but they have a common source: the wavefunction ψ(x).Probability interpretation treats the particle's position x ∈ R as a real variable and employs |ψ(x)| 2 as the probability of finding the particle at the position x.On the other hand, the complex trajectory interpretation uses the same wavefunction, ψ(x) but with x ∈ C to generate the complex quantum trajectory x(t) = x R (t) + ix I (t).employs |()| as the probability of finding the particle at the position .On the other hand, the complex trajectory interpretation uses the same wavefunction, () but with  ∈ ℂ to generate the complex quantum trajectory () =  () + i ().Here, the wavefunction is considered describing a particle's tunneling motion via a step barrier with height,  greater than the particle's energy, .Such a wavefunction () can be found analytically from the Schrödinger equation and then used in the quantum Hamilton mechanics [32] to produce the particle's complex position, () =  () + i (), and complex velocity,  () =  () + i (), during the tunneling process.As soon as () and  () are determined, the quantitative measures ( ,  ) = ( ,  ) and ( ,  ) = ( ,  ) can be used in Equation ( 1) to construct the tunneling-based YYDs.
Figure 6 shows how the complex tunneling velocities are mapped into the quantum YYD.The Figure 6, left, shows the positions  of the particle on the tunneling trajectory at nine moments  ,  = 0, 1, … , 8. A trajectory point  with velocity  ( ) + i ( ) is mapped into the point  on the YYD with radial coordinate  =  ( ) +  ( ) and azimuth coordinate  = 2 /, where  is the total tunneling time.As time  increases from 0 to , the area swept by the radial line segment  forms the Yang region (in red) of the YYD, while the area swept by the radial line segment  forms the Yin region (in black) of the YYD.(b) Trajectory interpretation uses the same wavefunction ψ(x) with x ∈ C to generate the complex trajectory x(t) = x R (t) + ix I (t), which is a continuous trajectory traversing the barrier.E denotes the particle energy.
Here, the wavefunction is considered describing a particle's tunneling motion via a step barrier with height, V 0 greater than the particle's energy, E. Such a wavefunction ψ(x) can be found analytically from the Schrödinger equation and then used in the quantum Hamilton mechanics [32] to produce the particle's complex position, x(t) = x R (t) + ix I (t), and complex velocity, .
x R (t) + i .
x I (t), during the tunneling process.As soon as x(t) and .
x(t) are determined, the quantitative measures (p x 2 I can be used in Equation ( 1) to construct the tunneling-based YYDs. Figure 6 shows how the complex tunneling velocities are mapped into the quantum YYD.The Figure 6, left, shows the positions P i of the particle on the tunneling trajectory at nine moments t i , i = 0, 1, . . ., 8. A trajectory point P i with velocity .
x R (t i ) + i .
x I (t i ) is mapped into the point P ′ i on the YYD with radial coordinate r i =  Among the nine trajectory points, P 0 is the entry point whose velocity .
x 0 is pure real, P 4 is the turning point whose velocity .
x 4 is pure imaginary, and P 8 is the exit point whose velocity is opposite to the entrance velocity .
x 8 = − .x 0 .Given the trajectory point P i with the accompanying velocity .
x R (t i ) + i .
x I (t i ), the inner and outer curves of the YYD can be determined as shown in the Figure 6, right.As time t evolves continuously, the tunneling particle follows a continuous complex trajectory x(t) = x R (t) + ix I (t) and, according to Equation (1), the inner curve of the velocity-based quantum YYD can be expressed in polar coordinates: and the outer curve is given by r = p A + p B = . x The azimuth coordinate θ is related to time t as θ = 2πt/T, 0 ≤ t ≤ T. In the time range of 0 ≤ t < T/2, the inner radius of the YYD is given by the squared real velocity r = p A = .
x 2 R , while in the time range of T/2 ≤ t < T, the inter radius is given by the squared imaginary velocity r = p B = .Among the nine trajectory points,  is the entry point whose velocity  is pure real,  is the turning point whose velocity  is pure imaginary, and  is the exit point whose velocity is opposite to the entrance velocity  = − .Given the trajectory point  with the accompanying velocity  ( ) =  ( ) + i ( ), the inner and outer curves of the YYD can be determined as shown in the Figure 6, right.As time  evolves continuously, the tunneling particle follows a continuous complex trajectory () =  () + i () and, according to Equation (1), the inner curve of the velocity-based quantum YYD can be expressed in polar coordinates: and the outer curve is given by  =  +  =  () +  ().The azimuth coordinate  is related to time  as  = 2/, 0 ≤  ≤ .In the time range of 0 ≤  < /2, the inner radius of the YYD is given by the squared real velocity  =  =  , while in the time range of /2 ≤  < , the inter radius is given by the squared imaginary velocity  =  =  .

Computing Complex Tunneling Trajectory
To construct the YYD by Equation ( 10), one needs the particle's tunneling velocity  () =  () + i (), which can be derived from the complex-extended Schrödinger where ℏ is the reduced Planck constant,  is the particle mass, and the coordinate  is defined on the complex plane  =  + i ∈ ℂ.Using the transformation (, ) =  ( , ) ℏ ⁄ , one transforms the Schrödinger equation ( 11) into the quantum Hamilton-Jacobi equation (QHJE): x R , and imaginary velocity, .
x I , (left) is mapped into a point P ′ i on the YYD (right) in such a way that the Yang (in red) region is formed by the time evolution of the squared real velocity,

Computing Complex Tunneling Trajectory
To construct the YYD by Equation (10), one needs the particle's tunneling velocity .
x R (t) + i .
x I (t), which can be derived from the complex-extended Schrödinger equation, iℏ where ℏ is the reduced Planck constant, m is the particle mass, and the coordinate x is defined on the complex plane x = x R + ix I ∈ C. Using the transformation Ψ(t, x) = e iS(t,x)/ℏ , one transforms the Schrödinger Equation ( 11) into the quantum Hamilton-Jacobi equation (QHJE): where p = ∂S/∂x is the canonical momentum conjugate to x, S(t, x) is the quantum action function, and H(t, x, p) is the quantum Hamiltonian defined as The quantum Hamiltonian H(t, x, p) (13) contains three items: (i) the particle's kinetic energy p 2 /2m, (ii) the applied potential V(t, x), and (iii) the internal quantum potential Q(t, x), defined as: Physics 2024, 6

975
It can be seen that the quantum potential is state-dependent and is uniquely determined by the wavefunction Ψ(t, x).If the applied potential V(t, x) is independent of time, that is, V(t, x) = V(x), then the wavefunction Ψ(t, x) can be further separated as The particle's momentum, p, according to the momentum definition in QHJE, is then given by Substituting the separated form of S(t, x) into QHJE (12) yields Equation ( 17) expresses the law of energy conservation, indicating that the total energy, E, of the particle represented by the Hamiltonian H(t, x, p) is a constant.The energy conservation law (17) is actually an alternative expression of the time-independent Schrödinger equation, as can be shown by substituting the momentum p (16) and the quantum potential Q (14) into Equation (17).
According to the Hamiltonian H(t, x, p) (17), the quantum Hamilton equations of motion read . .
where the sum of V and Q is defined as the total potential, The wave behavior of the particle is attributed to the quantum potential Q.If Q is constant, Equations ( 19) and ( 20) become classical equations of motion and the particle exhibits particle behavior.The quantum potential Q has a close relationship with the quantum probability.From Equation ( 21), one can see that except for a constant bias E, the magnitude of the total potential, |V T | = |V + Q|, is inversely proportional to the probability density function, |ψ(x)| 2 .Therefore, the smaller the probability density function |ψ(x)| 2 the higher the corresponding total potential is.Where ψ(x) is equal to zero, the total potential is infinitely high, so a particle with limited energy cannot reach it, and the probability of finding the particle there is zero.
For a given quantum state ψ(x), the velocity and acceleration of a particle moving in this state are computed by Equations ( 19) and (20), respectively: ..
The particle's trajectory x(t) in the state ψ(x) is obtained by integrating Equation ( 22) and the quantum force exerted on the particle at different positions is given by Equation (23).
The applied potential V(x) in Equation ( 18) to be considered here is the step barrier: where the height, V 0 , of the step barrier is greater than the energy, E, of the incident particle in order to explore quantum tunneling motion.The solution of the Schrödinger Equation ( 18) with V(x) (24) is obtained readily as where The range x R < 0 is the region of free motion, while x R ≥ 0 is the region where quantum tunneling motion takes place to generate the associated YYD.With the wavefunction ψ(x) given by Equation ( 25), the particle's velocity in the region x R ≥ 0 is obtained by Equation (22) as where the four constants A, B, C, and D are to be determined from the given initial conditions.To avoid the influence of physical units on the shape of the quantum YYD, the dimensionless variables √ 2mE/ℏ x → x and (2E/ℏ)t → t are adopted here to rewrite Equations ( 26) and (27) as where the parameter n = √ V 0 /E − 1 is a measure of the tunneling intensity.When n = 0, i.e., V 0 = E, the tunneling effect disappears completely.The four constants in Equations ( 28) and ( 29) are determined by the given initial conditions x(0 + ) = x(0 − ) = 0 and .
x(0 − ) = 1, from which In the region Re(x) < 0, the particle is free from the action of the applied potential and has a constant velocity, .
x = 1 given by Equation ( 29) with B = 0.By substituting ψ(x) = Ae ik 1 x for Re(x) < 0 into Equation ( 21), one obtains V T = V + Q = 0, which together with Equation (30) shows that the particle exhibits particle behavior with constant velocity, as shown in the region outside the step potential in Figure 5.
In the region Re(x) ≥ 0, the particle undergoes tunneling process with wave behavior caused by the quantum potential, Q.Once a particle enters the tunneling process, its motion extends from the real axis to the complex plane, and its complex trajectory can be solved analytically from Equation (28) as Physics 2024, 6 977 where the two constants α and β are determined from the given initial conditions as Due to |α| = 1, the constant α can be expressed as α = e iθ α with θ α = 2tan −1 (1/n).In terms of θ α and n, the particle's complex trajectory x(t) can be solved from Equation ( 31) as an explicit function of time Regarding the particle's velocity, Equation ( 31) is substituted into Equation ( 28) to yield an analytical expression for . x(t): β 2 e −2in 2 t + 4α = .
x R (t) + i .
x I (t). ( where the plus sign is taken for the time interval 0 ≤ t < T/2, and the minus sign is for T/2 ≤ t < T, where T is the total tunneling time to be determined.The quantum YYD shown in Figure 6 is drawn using . x R (t) and .
x I (t) computed from Equation (33).The remaining task is to determine the tunneling time T = t out , the turning point x tp , and the exit point x out of the tunneling trajectory, as shown in Figure 7.
gether with Equation (30) shows that the particle exhibits particle behavior with constant velocity, as shown in the region outside the step potential in Figure 5.
In the region Re() ≥ 0, the particle undergoes tunneling process with wave behavior caused by the quantum potential, .Once a particle enters the tunneling process, its motion extends from the real axis to the complex plane, and its complex trajectory can be solved analytically from Equation (28) as where the two constants  and  are determined from the given initial conditions as Due to || = 1, the constant  can be expressed as  =  with  = 2tan (1/ ).In terms of  and , the particle's complex trajectory () can be solved from Equation (31) as an explicit function of time Regarding the particle's velocity, Equation ( 31) is substituted into Equation ( 28) to yield an analytical expression for  (): where the plus sign is taken for the time interval 0 ≤  < /2, and the minus sign is for /2 ≤  < , where  is the total tunneling time to be determined.The quantum YYD shown in Figure 6 is drawn using  () and  () computed from Equation (33).The remaining task is to determine the tunneling time  =  , the turning point  , and the exit point  of the tunneling trajectory, as shown in Figure 7.The turning point  is the position where the particle reaches the maximum depth of penetration.Taking the absolute values of both sides of Equation (31), the relationship between  and  is obtained as follows: The turning point x tp is the position where the particle reaches the maximum depth of penetration.Taking the absolute values of both sides of Equation ( 31), the relationship between x R and x I is obtained as follows: The maximum value of x R occurs at cos(θ α + 2nx I ) = 1, and the related coordinate x I solution is Substituting the above x I (36) into Equation ( 35) yields the maximum penetration depth: When the particle reaches the maximum penetration depth, its position is the turning point x tp , so the complex coordinates of x tp can be expressed as Next, let us compute the time t tp and the velocity .
x tp when the particle reaches the turning point t tp .Squaring both sides of Equation ( 31) and using the turning condition, θ α + 2nx I = 0, one obtains t tp as The velocity .
x tp at the turning point is then determined by substituting the position x tp (39) into Equation ( 28) to give .
This result verifies what Figures 5 and 7 show, that the particle has only imaginary velocity component at the turning point.
Finally, let us calculate the position x out and time t out at the exit point.From Figure 7, one observes that the real coordinates x R of the entry point and exit point are both zero.By using this condition in Equation (35), one obtains from which x I = 0 or x I = −θ α /n, where the former corresponds to the position of the entrance, and the latter corresponds to the position of the exit.Therefore, the position coordinate of the particle at the exit is Comparing Equation (42) with Equation ( 38), one has Im(x out ) = 2Im x tp = −θ α /n, that is, the imaginary coordinate of the exit point is twice that of the turning point.Furthermore, the time for the particle to reach the exit point can be obtained by substituting x out (42) into Equation (31) as Comparing T with t tp (35), one finds that the time for the particle to reach the turning point is exactly half of the overall tunneling T. This result explains why in Figure 6, when the YYD is drawn halfway at P ′ 4 , the particle just reaches the turning point P 4 .The particle's velocity at the exit point is found by substituting x out = iθ α /n into Equation ( 28) as .
x out from the step barrier is the same as the x in , but in the opposite direction.The characteristics of the entire tunneling trajectory are summarized in Figure 7.

Time Evolution of Quantum YYD
The obtained tunneling-based YYD in Figure 6 is actually a graphical recorder of the particle's tunneling dynamics, which stores the particle's velocity at every moment so far.Figure 8 shows the time evolution of the YYD at nine moments, t 0 , t 1 , • • • , t 8 , corresponding to nine trajectory points in Figure 6.When time increases from 0 to T/2, the real velocity of the particle (in red) decreases from the maximum value at P ′ 0 to zero at P ′ 4 , while at the same time, the imaginary velocity of the particle (in black) increases from zero to the maximum value.In the time range of T/2 ≤ t < T, one sees the opposite trend that the real velocity of the particle increases from zero at P ′ 4 to the maximum value at P ′ 8 , and at the same time, the imaginary velocity decreases from the maximum value to zero.Only when the quantum tunneling process is completed one can see the complete quantum YYD, wherein the intertwined red and black pattern expresses the alternating changes of the real velocity and the imaginary velocity during the tunneling process.
Physics 2024, 6, FOR PEER REVIEW 16 as the incident velocity  , but in the opposite direction.The characteristics of the entire tunneling trajectory are summarized in Figure 7.

Time Evolution of Quantum YYD
The obtained tunneling-based YYD in Figure 6 is actually a graphical recorder of the particle's tunneling dynamics, which stores the particle's velocity at every moment so far.Figure 8 shows the time evolution of the YYD at nine moments,  ,  , ⋯ ,  , corresponding to nine trajectory points in Figure 6.When time increases from 0 to /2, the real velocity of the particle (in red) decreases from the maximum value at  to zero at  , while at the same time, the imaginary velocity of the particle (in black) increases from zero to the maximum value.In the time range of /2 ≤  < , one sees the opposite trend that the real velocity of the particle increases from zero at  to the maximum value at  , and at the same time, the imaginary velocity decreases from the maximum value to zero.Only when the quantum tunneling process is completed one can see the complete quantum YYD, wherein the intertwined red and black pattern expresses the alternating changes of the real velocity and the imaginary velocity during the tunneling process.Physics 2024, 6 980 Figure 9 shows the velocity distributions at nine trajectory points along the tunneling trajectory.It appears that the quantum YYD moves together with the particle, recording the velocity of the particle, and then presenting it graphically.The quantum YYD allows us to examine how the particle nature and the wave nature interchanges during the tunneling process.The imaginary velocity of a particle is a measure of its wave nature such that the greater the imaginary velocity of a particle is, the higher the uncertainty of its motion in the real space is, and the greater the wave nature it has [23].As shown in Figure 9, when the particle continues to tunnel towards the barrier, the particle's real velocity (in red) gradually decreases while its imaginary velocity (in black) increases.This means that as the particle tunnels deeper, its wave nature becomes larger.When the particle reaches the turning point P 4 , its imaginary velocity reaches the maximum, while its real velocity is zero, which means that the particle behaves completely like a wave at the turning point.In the process of returning from point P 4 to the exit point P 8 , the real velocity increases while the imaginary velocity decreases, indicating that the particle nature gradually recovers towards the exit, at which the imaginary velocity of the particle is zero and the wave nature of the particle disappears completely.After leaving the step potential, the particle returns to the free space, and its motion is again described by Equation ( 29) with A = 0, which yields .
x = −1 and indicates a classical motion along the negative real axis.
Physics 2024, 6, FOR PEER REVIEW 17 Figure 9 shows the velocity distributions at nine trajectory points along the tunneling trajectory.It appears that the quantum YYD moves together with the particle, recording the velocity of the particle, and then presenting it graphically.The quantum YYD allows us to examine how the particle nature and the wave nature interchanges during the tunneling process.The imaginary velocity of a particle is a measure of its wave nature such that the greater the imaginary velocity of a particle is, the higher the uncertainty of its motion in the real space is, and the greater the wave nature it has [23].As shown in Figure 9, when the particle continues to tunnel towards the barrier, the particle's real velocity (in red) gradually decreases while its imaginary velocity (in black) increases.This means that as the particle tunnels deeper, its wave nature becomes larger.When the particle reaches the turning point  , its imaginary velocity reaches the maximum, while its real velocity is zero, which means that the particle behaves completely like a wave at the turning point.In the process of returning from point  to the exit point  , the real velocity increases while the imaginary velocity decreases, indicating that the particle nature gradually recovers towards the exit, at which the imaginary velocity of the particle is zero and the wave nature of the particle disappears completely.After leaving the step potential, the particle returns to the free space, and its motion is again described by Equation ( 29) with  = 0, which yields  = −1 and indicates a classical motion along the negative real axis.The tunneling-based quantum YYD moves with the particle and records the particle's velocity.From the entry point to the turning point, the real velocity component (in red) decreases while the imaginary velocity component (in black) increases, indicating that the particle nature decreases and the wave nature increases.When reaching the turning point, the particle has only an imaginary velocity component (pure wave nature).From the turning point to the exit point, the red area increases while the black area decreases, indicating that the particle nature increases and the wave nature decreases.When the particle reaches the exit, the particle is left with only real velocity component and exhibits pure particle nature.The tunneling-based quantum YYD moves with the particle and records the particle's velocity.From the entry point to the turning point, the real velocity component (in red) decreases while the imaginary velocity component (in black) increases, indicating that the particle nature decreases and the wave nature increases.When reaching the turning point, the particle has only an imaginary velocity component (pure wave nature).From the turning point to the exit point, the red area increases while the black area decreases, indicating that the particle nature increases and the wave nature decreases.When the particle reaches the exit, the particle is left with only real velocity component and exhibits pure particle nature.
Physics 2024, 6 981 In quantum mechanics, the wave-like behavior of a tunneling particle is described by the probability density function |ψ(x)| 2 as shown in Figure 5a.As soon as the particle imaginary velocity is uniquely determined by |ψ(x)| 2 , as derived below in this Section, the imaginary velocity serves as an alternative indicator of the particle's wave behavior.The relation between the particle's velocity .
x and the wavefunction ψ(x) is given by Equation ( 16) as To separate the complex velocity .
x into its real part and imaginary part, let us express wavefunction ψ(x) as ψ(x) = |ψ|e iθ ψ , where θ ψ ∈ R is the phase of ψ(x).Using this expression of ψ(x) in Equation ( 44 from which the imaginary velocity is .
One can see that the imaginary velocity .
x I is uniquely determined by the probability density function |ψ(x)| 2 and can be interpreted as the change rate of |ψ(x)| per unit magnitude of |ψ(x)|.Therefore, when the particle's imaginary velocity .
x I reaches the maximum, it is the place where the wave magnitude |ψ(x)| changes most drastically, that is, where the wave behavior of the particle is most significant.Corresponding to the quantum tunneling process, this is exactly where the particle reaches the turning point, as shown in Figure 9.
The driving force behind the tunneling dynamics comes from the total potential V T = V 0 + Q, as shown by Equation (23).Since the applied potential V 0 is a constant, the quantum potential Q is the dominant potential governing the particle's tunneling motion.It is precisely because of the effect of quantum potential Q that the particle is forced to turn back after penetrating the step barrier for a certain distance.The deeper the particle penetrates the step barrier, the greater resistance it encounters.The maximum depth to which a particle penetrates is where it tolerates the greatest resistance.
One determine the total potential V T by substituting the wavefunction ψ(x) from Equation (25) into Equation ( 21) to obtain the following result: where V T is expressed in a dimensionless form.Once the total potential is obtained, the quantum force exerted on the particle during the tunneling process can be found as If expressed in terms of the particle's velocity .
x, F Q can be rewritten in a simple form as Physics 2024, 6 982 With the velocity .
x tp (40), one can determine the quantum force exerted on the particle at the turning point as As expected, the quantum force experienced by the particle at the turning point is in the negative x R direction, preventing the particle from continuing to penetrate in the positive x R direction.Equation (50) shows that as n = √ V 0 /E − 1 increases, the resistance encountered by the particle becomes greater, resulting in a smaller penetration depth (37).
From the perspective of probability, the higher the step barrier V 0 is, the lower the probability of tunneling is.The transition from the probability interpretation to the trajectory interpretation of the tunneling phenomenon requires the assistance of the quantum potential Q(x). Figure 10 shows the complex trajectories of the particle moving over the total potential surface V T = Q(x) + V 0 in the two cases of n = 1 and n = 5, indicating that when n is larger, the total potential surface becomes steeper, and the depth that the particle can penetrate is shallower.The turning point P 4 is the place with the highest total potential in the overall tunneling trajectory, while the entry point P 0 and the exit point P 8 are the places with the lowest total potential.
Physics 2024, 6, FOR PEER REVIEW 19 As expected, the quantum force experienced by the particle at the turning point is in the negative  direction, preventing the particle from continuing to penetrate in the positive  direction.Equation (50) shows that as  =  / − 1 increases, the resistance encountered by the particle becomes greater, resulting in a smaller penetration depth (37).
From the perspective of probability, the higher the step barrier  is, the lower the probability of tunneling is.The transition from the probability interpretation to the trajectory interpretation of the tunneling phenomenon requires the assistance of the quantum potential ().Figure 10 shows the complex trajectories of the particle moving over the total potential surface  = () +  in the two cases of  = 1 and  = 5, indicating that when  is larger, the total potential surface becomes steeper, and the depth that the particle can penetrate is shallower.The turning point  is the place with the highest total potential in the overall tunneling trajectory, while the entry point  and the exit point  are the places with the lowest total potential.Figure 10.The complex trajectories of a tunneling particle on the total potential surface.When a particle undergoes quantum tunneling towards the real axis,  , the height of the quantum potential increases with the tunneling depth.When the particle reaches the turning point ( ), the quantum potential reaches its maximum.The higher the step barrier height,  , the steeper the tunnel trajectory and the shallower the tunneling depth.When  approaches infinity, the tunneling depth of the particle is zero, that is, the particle is immediately reflected when it encounters the surface of the step barrier.The red arrows along the trajectory denote the moving directions of the particle.

Growth of Quantum YYD with Tunneling Intensity
The quantum YYD based on tunneling velocity ( ,  ) = ( ,  ) is not always an ideal YYD, and its shape depends on the intensity of the quantum tunneling effect.It is found here that the stronger the quantum tunneling effect the clearer and more complete the quantum YYD is.Once the intensity of the tunneling effect approaches zero, the quantum YYD disappears.The intensity of the tunneling effect is related to the relative magnitude between the incident particle's energy  and the step barrier height,  .This relative magnitude is measured by the dimensionless parameter  =  / − 1.The larger the value of  is, the smaller the particle's energy is compared to the barrier height, and a larger tunneling effect is required for the particle to penetrate the barrier.
What makes us curious is how the intensity of the tunneling effect affects the quantum YYD. Figure 11 shows the quantum YYDs corresponding to different values of .One can see that the quantum YYD is born from  = 0, initially like a melon seed, and gradually grows into a perfect appearance as n increases.As  ≫ 1, the quantum YYD no longer changes, and the corresponding polar Equation (10) converges to Figure 10.The complex trajectories of a tunneling particle on the total potential surface.When a particle undergoes quantum tunneling towards the real axis, x R , the height of the quantum potential increases with the tunneling depth.When the particle reaches the turning point (P 4 ), the quantum potential reaches its maximum.The higher the step barrier height, V 0 , the steeper the tunnel trajectory and the shallower the tunneling depth.When V 0 approaches infinity, the tunneling depth of the particle is zero, that is, the particle is immediately reflected when it encounters the surface of the step barrier.The red arrows along the trajectory denote the moving directions of the particle.

Growth of Quantum YYD with Tunneling Intensity
The quantum YYD based on tunneling velocity (p A , p B ) = .
x 2 I is not always an ideal YYD, and its shape depends on the intensity of the quantum tunneling effect.It is found here that the stronger the quantum tunneling effect the clearer and more complete the quantum YYD is.Once the intensity of the tunneling effect approaches zero, the quantum YYD disappears.The intensity of the tunneling effect is related to the relative magnitude between the incident particle's energy E and the step barrier height, V 0 .This relative magnitude is measured by the dimensionless parameter n = √ V 0 /E − 1.The larger the value of n is, the smaller the particle's energy is compared to the barrier height, and a larger tunneling effect is required for the particle to penetrate the barrier.
Physics 2024, 6 983 What makes us curious is how the intensity of the tunneling effect affects the quantum YYD. Figure 11 shows the quantum YYDs corresponding to different values of n.One can see that the quantum YYD is born from n = 0, initially like a melon seed, and gradually grows into a perfect appearance as n increases.As n ≫ 1, the quantum YYD no longer changes, and the corresponding polar Equation (10) where T = π/n 2 is the dimensionless tunneling time.The satisfaction of the complementarity condition p A + p B = cos 2 n 2 t + sin 2 n 2 t = 1 and the continuity condition p A (π) = p B (2π) = 0 show that the quantum YYD described by Equation ( 51) is an ideal YYD, whose inner curve is composed of two sections of cardioid.For small value of n, the corresponding YYDs as shown in Figure 11   The numbers refer to the nine trajectory points in Figure 6.
Comparing Equation (51) with Equation ( 5), it was quite surprising to find that the ideal YYD generated by the tunneling motion resembles Bohr's YYD designed in his coat of arms.If Bohr drew YYD based on Equation (51) instead of Equation ( 5), then the coat of arms he designed would completely reflect his principle of complementarity not only in image but also in theory.
Figure 11 shows that when the parameter  =  / − 1 ≫ 1, the tunneling-based YYD converges to a steady-state pattern, known as the ideal YYD.Since the ideal YYD is drawn based on the particle's tunneling trajectory, it is expected that the tunneling trajectory also converges to a steady-state pattern as  ≫ 1 .Let us verify that this steady-state trajectory is a semicircle on the complex plane.
When n is small, the real velocity component .
x R (particlelike behavior) is dominant in the YYD, leading to a deformed YYD with a long strip shape concentrating near θ = 0.When n increases, the imaginary velocity component .
x I (the wavelike behavior) becomes comparable to .
x R , and the resulting YYD is fuller.As n ≫ 1, the velocity-based YYD approaches the ideal YYD whose inner curve is composed of two sections of cardioid and whose outer curve is a unit circle.The numbers refer to the nine trajectory points in Figure 6.
Physics 2024, 6 984 Comparing Equation (51) with Equation ( 5), it was quite surprising to find that the ideal YYD generated by the tunneling motion resembles Bohr's YYD designed in his coat of arms.If Bohr drew YYD based on Equation (51) instead of Equation ( 5), then the coat of arms he designed would completely reflect his principle of complementarity not only in image but also in theory.
Figure 11 shows that when the parameter n = √ V 0 /E − 1 ≫ 1, the tunneling-based YYD converges to a steady-state pattern, known as the ideal YYD.Since the ideal YYD is drawn based on the particle's tunneling trajectory, it is expected that the tunneling trajectory also converges to a steady-state pattern as n ≫ 1.Let us verify that this steadystate trajectory is a semicircle on the complex plane.
Figure 12 shows the tunneling trajectories of the particle on the complex plane solved from Equation (28) under several different values of n.The main feature of this series of trajectories indicates that as the value of n increases, the covering range of the particle's position x becomes narrower, as expressed by the relation |nx| ≈ 1/n.Therefore, when n is large, one obtains the result |nx| ≪ 1.With this observation, the exponential function e nx can be estimated by its first-order expansion 1 + nx, which is then used in Equation (28) to yield the particle's velocity for n ≫ 1 as where note that |nx| ≪ 1 and 1 + ix ≈ 1.With the initial condition x(0) = 0, the solution of Equation ( 48) can be found readily as Physics 2024, 6, FOR PEER REVIEW 21 Figure 12 shows the tunneling trajectories of the particle on the complex plane solved from Equation (28) under several different values of .The main feature of this series of trajectories indicates that as the value of  increases, the covering range of the particle's position  becomes narrower, as expressed by the relation || ≈ 1/.Therefore, when  is large, one obtains the result || ≪ 1.With this observation, the exponential function  can be estimated by its first-order expansion 1 + , which is then used in Equation (28) to yield the particle's velocity for  ≫ 1 as where note that || ≪ 1 and 1 + i ≈ 1.With the initial condition  0 = 0, the solution of Equation ( 48) can be found readily as From Equation (53), one can obtain the position, velocity, and time when the particle reaches the turning point and the exit point.These results are identical to those summarized in Figure 7, when  is large.For example, the overall tunneling time  is the time that the real coordinate  () returns to zero, and from Equation ( 53), one has  = / , which is identical to Equation (43).In addition, the penetration depth  is the maximum value of the real part of (), which is equal to 1/ from Equation (53), and this result is identical to Equation (37) as  is large.The tunneling trajectory expressed by Equation ( 53) is a semicircle on the complex plane with radius 1/ and center (0, −1/ ): It is based on this ideal tunneling trajectory that one can construct the ideal quantum YYD.Differentiating Equation (53) with respect to time, one obtains the real and imaginary velocities of the particle as Using Equation (55) in Equation (10), one obtains the ideal YYD as given by Equation (51).
The YYD derived previously is based on the velocity distribution of a tunneling particle, so a natural question is whether the quantum YYD can also be constructed by the position distribution of a tunneling particle.The answer is yes, and surprisingly, From Equation (53), one can obtain the position, velocity, and time when the particle reaches the turning point and the exit point.These results are identical to those summarized in Figure 7, when n is large.For example, the overall tunneling time T is the time that the real coordinate x R (t) returns to zero, and from Equation (53), one has T = π/n 2 , which is identical to Equation (43).In addition, the penetration depth δ p is the maximum value of the real part of x(t), which is equal to 1/n 2 from Equation (53), and this result is identical to Equation (37) as n is large.

985
The tunneling trajectory expressed by Equation ( 53) is a semicircle on the complex plane with radius 1/n 2 and center 0, −1/n 2 : It is based on this ideal tunneling trajectory that one can construct the ideal quantum YYD.Differentiating Equation (53) with respect to time, one obtains the real and imaginary velocities of the particle as .
The YYD derived previously is based on the velocity distribution of a tunneling particle, so a natural question is whether the quantum YYD can also be constructed by the position distribution of a tunneling particle.The answer is yes, and surprisingly, when the value of n is large, the YYDs obtained by these two different methods are the same.To search for the complementary pair (A, B) for the position-based YYD, let us consider Equation (54), which under the coordinate transformation x R = n 2 x R and x I = n 2 x I + 1 can be expressed as x 2 R + x 2 I = 1.Hence, the new pair (A, B) = (x R , x I ) with the quantitative measures becomes complementary by noting p A + p B = x 2 I + x 2 R = 1.Therefore, the outer curve of the position-based YYD is r = p A + p = 1, which is a unit circle, and the corresponding inner curve defined by Equation (1) becomes where x R (t) and x I (t) are given by Equation (53).Equation (51) and Equation (57) are equivalent, which implies that the velocity-based YYD and the position-based YYD become identical as n ≫ 1.
For a general value of n not satisfying n ≫ 1, the symmetric center of the tunneling trajectory is not at (x R , x I ) = 0, −1/n 2 , but at 0, −(1/n)tan −1 (1/n) , as can be seen from Figure 7. Hence, for a general position-based YYD, the quantitative measures p A and p B in Equation (56) are replaced by where the time history of the complex position x(t) = x R (t) + ix I (t) is solved from Equation ( 28) for a general value of n.
Figure 13 shows how the position-based YYD evolves with the tunneling intensity n.The YYD emerges from the origin with n = 0, initially resembling a walnut in appearance, and grows into an ideal YYD as n increases.Regardless of whether the velocity or position data are used to construct the YYD, one finds that when the value of n is large, the two methods result in the same ideal YYDs, whose outer curve approaches a perfect circle, and whose inner curve approaches a cardioid.
Figure 13 shows how the position-based YYD evolves with the tunneling intensity .The YYD emerges from the origin with  = 0, initially resembling a walnut in appearance, and grows into an ideal YYD as  increases.Regardless of whether the velocity or position data are used to construct the YYD, one finds that when the value of n is large, the two methods result in the same ideal YYDs, whose outer curve approaches a perfect circle, and whose inner curve approaches a cardioid.

Conclusions
In this paper, it is found that Bohr's YYD in quantum tunneling dynamics and discussed four special forms of YYD relating to BPC.Bohr's YYD in his coat of arms is a popular YYD consisting of a unit circle and two inscribed circles, as shown in Figure 14a.Except for being like the tunneling-based YYD in appearance, this popular YYD has no link to physical systems.The second type of YYD is the traditional Taoist YYD, which is an ancient realization of BPC as shown in Figure 14b.The third type of YYD in Figure 14c visualizes a particle's tunneling motion by displaying the time evolution of the particle's position and velocity.The fourth type of YYD in Figure 14d is drawn based on the data of the two-path interference experiments for testing the wave-particle duality.Taoist YYD has been reproduced in this article by BPC, which, in turn, justifies Bohr's use of Taoist YYD as the icon of his principle.The only pity is that the YYD designed by Bohr in his cost of arms is not Taoist YYD drawn based on his principle.However, Bohr did not need to borrow Taoist YYD to symbolize the wave-particle duality, because the wave-particle duality has its own representative YYDs, i.e., the tunneling-based YYD and the interference-based YYD, as revealed in this article.The investigation in this paper shows that Bohr's YYD and its quantum version, the tunneling-based YYD, belong to the same category, while Taoist YYD and its quantum version, the interference-based YYD, belong to another category.The types of YYD are not limited to those shown in Figure 14.Actually, every system in nature has its accompanying YYD that can be constructed according to Equation (1).Any complementary pair A and B in a system, whether classical or quantum, can be employed to construct the system's YYD, if their quantitative measures  and  are available within a period.In ancient China, the complementary pair could be the daytime and nighttime in a day.In classical mechanics, the complementary pair could be the kinetic energy and the gravitational potential energy of a particle during vertical projectile motion or could be the kinetic energy and the elastic potential energy of a spring-mass system.In quantum mechanics, the complementary pair could be the particle nature and the wave nature of a photon in a two-path interference experiment or could be the real and imaginary velocities of a tunneling particle.If, as Bohr said, complementarity is a universal principle, the YYD constructed in this paper is a universal image found all over in nature.
When the complementary pair satisfies the conditions of continuity and complementarity, the obtained YYD turns out to be an ideal one, whose outer curve is a perfect The types of YYD are not limited to those shown in Figure 14.Actually, every system in nature has its accompanying YYD that can be constructed according to Equation (1).Any complementary pair A and B in a system, whether classical or quantum, can be employed to construct the system's YYD, if their quantitative measures p A and p B are available within a period.In ancient China, the complementary pair could be the daytime and nighttime in a day.In classical mechanics, the complementary pair could be the kinetic energy and the gravitational potential energy of a particle during vertical projectile motion or could be the kinetic energy and the elastic potential energy of a spring-mass system.In quantum mechanics, the complementary pair could be the particle nature and the wave nature of a photon in a two-path interference experiment or could be the real and imaginary velocities of a tunneling particle.If, as Bohr said, complementarity is a universal principle, the YYD constructed in this paper is a universal image found all over in nature.
When the complementary pair satisfies the conditions of continuity and complementarity, the obtained YYD turns out to be an ideal one, whose outer curve is a perfect circle with continuous inner curve.One can then judge the degree to which the system satisfies BPC based on the deviation of the system's YYD from the ideal YYD.

Figure 1 .
Figure 1.(a) Bohr designed his own coat of arms that featured a symbol of Yin and Yang, and a motto in Latin "contraria sunt complementa" (opposites are complementary) [6].(b) Two states A and B of a system are said to be complementary, if they are jointly complete and mutually exclusive, or equivalently,  +  = 1, where  and  are the quantitative measures for A and B, respectively.

Figure 1 .
Figure 1.(a) Bohr designed his own coat of arms that featured a symbol of Yin and Yang, and a motto in Latin "contraria sunt complementa" (opposites are complementary) [6].(b) Two states A and B of a system are said to be complementary, if they are jointly complete and mutually exclusive, or equivalently, p A + p B = 1, where p A and p B are the quantitative measures for A and B, respectively.

Figure 2 .
Figure 2. A schematic diagram of constructing Taoist YYD according to the principle of complementarity.(, ) a qualified complementary pair of a periodic system satisfying the condition  ( ) +  ( ) = 1 at every moment  within the period 0 ≤  ≤ .The inner curve of Taoist YYD is defined by the polar coordinates ( ,  ), where  =  ( ) for 0 ≤  ≤ /2, and  =  ( ), for /2 ≤  ≤ , and the outer curve is defined by  =  ( ) +  ( ).The quantitative measures  ( ) and  ( ) used in the traditional Taoist YYD are the proportions of the states A and B in the entire system at time  , just like the proportions of day and night on a certain day  of the year.

Figure 2 .
Figure 2. A schematic diagram of constructing Taoist YYD according to the principle of complementarity.(A,B) a qualified complementary pair of a periodic system satisfying the condition p A (t k ) + p B (t k ) = 1 at every moment t k within the period 0 ≤ t k ≤ T. The inner curve of Taoist YYD is defined by the polar coordinates (r k , θ k ), where r k = p A (t k ) for 0 ≤ t k ≤ T/2, and r k = p B (t k ), for T/2 ≤ t k ≤ T, and the outer curve is defined by r k = p A (t k ) + p B (t k ).The quantitative measures p A (t k ) and p B (t k ) used in the traditional Taoist YYD are the proportions of the states A and B in the entire system at time t k , just like the proportions of day and night on a certain day t k of the year.

Figure 3 .
Figure 3. Bohr's YYD designed in his coat of arms is composed of a unit circle and two inscribed circles.The polar equation of the inner curve is given by  = ± cos .

Figure 3 .
Figure 3. Bohr's YYD designed in his coat of arms is composed of a unit circle and two inscribed circles.The polar equation of the inner curve is given by r = ±cos θ.

Figure 4 .
Figure 4.The quantum YYDs obtained from the two-path interference experiments.(a) The schematic diagram of the Mach-Zehnder interferometer used to measure the predictability  and the visibility .See text for details.(b) The ideal quantum YYD constructed by using the quantitative measures ( ,  ) = ( ,  ) without path loss in the interferometer.L1 and L2 denote the losses of the two paths, and R1 and R2 are the adjustable reflectivities of the beam splitters in the input and output ports, respectively.(c,d) The deformed YYDs caused by the path loss  ., whose outer curves remain the same as the ideal YYD, but whose inner curves have a discontinuity  at  = .(e,f) The deformed YYDs with ( ,  ) measured from the other configuration of the experiment by interchanging the input port and the output port in (a).The resulting YYD has a deformation in the outer curve without causing the discontinuity in the inner curve.

Figure 4 .
Figure 4.The quantum YYDs obtained from the two-path interference experiments.(a) The schematic diagram of the Mach-Zehnder interferometer used to measure the predictability P and the visibility V. See text for details.(b) The ideal quantum YYD constructed by using the quantitative measures (p A , p B ) = P 2 , V 2 without path loss in the interferometer.L 1 and L 2 denote the losses of the two paths, and R 1 and R 2 are the adjustable reflectivities of the beam splitters in the input and output ports, respectively.(c,d) The deformed YYDs caused by the path loss L 2 ., whose outer curves remain the same as the ideal YYD, but whose inner curves have a discontinuity δp A at θ = π.(e,f) The deformed YYDs with (p A , p B ) measured from the other configuration of the experiment by interchanging the input port and the output port in (a).The resulting YYD has a deformation in the outer curve without causing the discontinuity in the inner curve.

Figure 5 .
Figure 5. Comparing two quantum interpretations of tunneling phenomenon.(a) Probability interpretation uses the wavefunction, () with  ∈ ℝ to describe the particle's wave behavior within the step barrier and treats |()| as the probability of finding the particle at position .(b) Trajectory interpretation uses the same wavefunction () with  ∈ ℂ to generate the complex trajectory () =  () + i (), which is a continuous trajectory traversing the barrier.E denotes the particle energy.

Figure 5 .
Figure 5. Comparing two quantum interpretations of tunneling phenomenon.(a) Probability interpretation uses the wavefunction, ψ(x) with x ∈ R to describe the particle's wave behavior within the step barrier and treats |ψ(x)| 2 as the probability of finding the particle at position x.(b) Trajectory interpretation uses the same wavefunction ψ(x) with x ∈ C to generate the complex trajectory x(t) = x R (t) + ix I (t), which is a continuous trajectory traversing the barrier.E denotes the particle energy.

x 2 R 2 I 2 R 2 I
(t i ) + .x(t i ) and azimuth coordinate θ i = 2πt i /T, where T is the total tunneling time.As time t i increases from 0 to T, the area swept by the radial line segment .xforms the Yang region (in red) of the YYD, while the area swept by the radial line segment .xforms the Yin region (in black) of the YYD.

11 Figure 6 .
Figure 6.The mapping of the complex tunneling velocity into the quantum YYD.A trajectory point  with real velocity,  , and imaginary velocity,  , (left) is mapped into a point  on the YYD (right) in such a way that the Yang (in red) region is formed by the time evolution of the squared real velocity,  , while the Yin (in black) region is formed by the squared imaginary velocity,  .See text for details.

Figure 6 .
Figure 6.The mapping of the complex tunneling velocity into the quantum YYD.A trajectory point P i with real velocity, .

x 2 R 2 I
, while the Yin (in black) region is formed by the squared imaginary velocity, .x .See text for details.

Figure 7 .
Figure 7.The analytical expression of the tunneling trajectory.The time, position, and velocity of the three trajectory points (entry point  , turning point  , and exit point  ) are expressed analytically in terms of the parameter  =  / − 1.The red arrows along the trajectory denote the moving directions of the particle.See text for details.

Figure 7 .
Figure 7.The analytical expression of the tunneling trajectory.The time, position, and velocity of the three trajectory points (entry point x in , turning point x tp , and exit point x out ) are expressed analytically in terms of the parameter n = √ V 0 /E − 1.The red arrows along the trajectory denote the moving directions of the particle.See text for details.

Figure 8 .
Figure 8.The time evolution of the YYD in the quantum tunneling process.The areas swept by the radius  =  =  and  =  =  are plotted at nine instances  = /8,  = 0, 1, ⋯ , 8. The Yang (in red) region is the area swept by  =  (), showing the distribution of particle's real velocity, while the Yin (in black) region is the area swept  =  (), showing the distribution of particle's imaginary velocity.

Figure 8 . 2 R 2 I
Figure 8.The time evolution of the YYD in the quantum tunneling process.The areas swept by the radius r = p A = .x 2 R and r = p B = .x 2 I are plotted at nine instances t = kT/8, k = 0, 1, • • • , 8. The Yang (in red) region is the area swept by r = p A (θ), showing the distribution of particle's real velocity, while the Yin (in black) region is the area swept r = p B (θ), showing the distribution of particle's imaginary velocity.

Figure 9 .
Figure 9.The tunneling-based quantum YYD moves with the particle and records the particle's velocity.From the entry point to the turning point, the real velocity component (in red) decreases while the imaginary velocity component (in black) increases, indicating that the particle nature decreases and the wave nature increases.When reaching the turning point, the particle has only an imaginary velocity component (pure wave nature).From the turning point to the exit point, the red area increases while the black area decreases, indicating that the particle nature increases and the wave nature decreases.When the particle reaches the exit, the particle is left with only real velocity component and exhibits pure particle nature.

Figure 9 .
Figure 9.The tunneling-based quantum YYD moves with the particle and records the particle's velocity.From the entry point to the turning point, the real velocity component (in red) decreases while the imaginary velocity component (in black) increases, indicating that the particle nature decreases and the wave nature increases.When reaching the turning point, the particle has only an imaginary velocity component (pure wave nature).From the turning point to the exit point, the red area increases while the black area decreases, indicating that the particle nature increases and the wave nature decreases.When the particle reaches the exit, the particle is left with only real velocity component and exhibits pure particle nature.

Figure 11 .
Figure 11.The evolution of the velocity-based YYD with the intensity of the tunneling effect,  =  / − 1.When  is small, the real velocity component  (particlelike behavior) is dominant in the YYD, leading to a deformed YYD with a long strip shape concentrating near  = 0.When  increases, the imaginary velocity component  (the wavelike behavior) becomes comparable to  , and the resulting YYD is fuller.As  ≫ 1, the velocity-based YYD approaches the ideal YYD whose inner curve is composed of two sections of cardioid and whose outer curve is a unit circle.The numbers refer to the nine trajectory points in Figure 6.

Figure 11 .
Figure 11.The evolution of the velocity-based YYD with the intensity of the tunneling effect, n= √ V 0 /E − 1.When n is small, the real velocity component

Figure 12 .
Figure 12.The change of the particle's tunneling trajectory on the complex plane with the parameter .For large , the particle's position  within the barrier is more restricted, as expressed by the relation || ≪  .As  increases, the tunneling trajectory of the particle eventually approaches a perfect semicircle.

Figure 12 .
Figure12.The change of the particle's tunneling trajectory on the complex plane with the parameter n.For large n, the particle's position x within the barrier is more restricted, as expressed by the relation |x| ≪ n −1 .As n increases, the tunneling trajectory of the particle eventually approaches a perfect semicircle.

Figure 13 .
Figure 13.The growth of the position-based YYD with the tunneling intensity n.On plotting the position-based YYD, the quantitative measures are given by p A = n 2 x I + ntan −1 (1/n) 2 andp B = n 2 x R 2 .As n increases from 0, the YYD grows from the origin, initially like the shape of a walnut.As n ≫ 1, the position-based YYD becomes identical to the velocity-based YYD.The numbers refer to the nine trajectory points in Figure6.

Figure 14 .
Figure 14.The interconnection of four YYDs and their relationship with BPC.(a) Bohr's YYD is based on a popular construction to symbolize the principle of complementarity.(b) Taoist YYD is proved to be an ancient realization of BPC, which justifies the role of YYD in Bohr's principle.(c) The tunneling-based YYD visualizes a particle's tunneling motion via a step barrier by displaying the time evolution of the particle's velocity.(d) The interference-based YYD is drawn by the data from two-path interferometers to test BPC.According to the inner curves of the four YYDs, Bohr's YYD and its quantum version, the tunneling-based YYD, belong to the same category, while Taoist YYD and its quantum version, the interference-based YYD, belong to another category.

Figure 14 .
Figure 14.The interconnection of four YYDs and their relationship with BPC.(a) Bohr's YYD is based on a popular construction to symbolize the principle of complementarity.(b) Taoist YYD is proved to be an ancient realization of BPC, which justifies the role of YYD in Bohr's principle.(c) The tunneling-based YYD visualizes a particle's tunneling motion via a step barrier by displaying the time evolution of the particle's velocity.(d) The interference-based YYD is drawn by the data from two-path interferometers to test BPC.According to the inner curves of the four YYDs, Bohr's YYD and its quantum version, the tunneling-based YYD, belong to the same category, while Taoist YYD and its quantum version, the interference-based YYD, belong to another category.