1. Introduction
The existence of a formal analogy among the equations that govern wave propagation, allows to extend, in some cases, the results to phenomena whose study originates in areas of different physical disciplines. This is the case, for example, of the phenomenon of the so-called superluminal propagation of evanescent waves, whose name directly recalls the field of electromagnetism but which also has close analogies in quantum physics and, as demonstrated in this paragraph, even in acoustics. The unifying concept is that of the propagation of an impulse or a wave packet in a medium with a strong or anomalous dispersion: a packet of acoustic, electromagnetic or quantum wave pobability density waves can therefore be considered equivalent.
Let us therefore imagine not a monochromatic wave but a wavepacket limited in time and space with a central carrier frequency and an envelope with modulated amplitude which, therefore, introduces components at different frequencies through Fourier analysis: the so-called phase velocity
can be identified as the velocity of the crests of the carrier frequency, but already Lord Rayleigh [
1] identified the fact that the packet envelope moves with the group velocity
(at first order) which has, in the case of anomalous dispersion, peculiar characteristics e.g. can become negative when the frequency decreases rather than increases with increasing
k.
Anomalous dispersion was first studied for mechanical oscillators [
2] and later, by Sommerfeld and Brillouin, in materials that absorb light in which the group velocity can be greater than c (velocity of light in vacuum) or even negative within the absorption [
3]. The phase velocity can be greater than c in many cases, for example inside waveguides but this does not create problems because, representing the velocity of a continuous series of crests, this does not carry the signal information.
On the other end, if the group velocity is identified as the velocity at which a signal (information) travels, the problem of the limit of the principle of causality set by the constancy of the velocity of light in the vacuum of special relativity immediately arises, but this identification is possible only in cases of "normal" dispersion where the deformation of the packet is slight and the group velocity remains below the velocity c.
Sommerfeld and Brillouin then continued to define signal velocity as that of the point of rise of the signal intensity equal to half its steady-state value, which can be demonstrated to proceed at the group velocity, but they introduced the concept of "velocity of the wave front" ideally defined as that of a discontinuity in the signal of the step function type with infinite derivative and it is this last velocity, certainly carrying information, which is subject to the limitation of the velocity
c just as the velocity, to be defined, with which the energy is transported [
3,
4].
All these velocities, in a non-dispersive medium coincide and are equivalent, for elastic waves, to the bulk transverse waves velocity of that medium while, for electromagnetic waves, is the velocity of light in that medium or c if it is a vacuum. However, it must be noted that, the definition of a wave front with infinite derivative (vertical slope) is only ideal, requiring an infinite bandwidth in the frequency components which is unachievable in any real signal generator. In reality, physical signal generators only produce finite spectra due to their natural inertia in reaching the steady-state amplitude (impossibility of a vertical rise of the wave front) and the necessarily finite content of energy in the signal because an infinite frequency bandwidth would require, due to the Plank relation , an infinite energy.
The limitation of the frequency bandwidth for a physical detector highlights the limits of reasoning based on classical physics. In the latter case, a "classical" detector can detect a theoretically small quantity of energy as desired, whereas a physical detector needs at least a quantum of energy
. Regardless of the interpretation, when we talk about propagation at superluminal velocity we are referring to a clear and measurable effect that involves the group velocity and therefore requires the use of a wave packet. This effect has been predicted theoretically and measured experimentally in different conditions: for example in anomalous dispersion zones near an absorption line [
5,
6,
7], in nonlinear [
8] and linear gain lines [
9,
10,
11], in a active plasma [
12], optically [
13] and finally in a tunnel effect barrier [
14,
15,
16].
This last physical situation in which superluminal wavepackets are found, is the one we want to take into consideration here because there is a close analogy between the quantum tunneling effect of a particle through a potential barrier and the crossing of a zone, limited in space, in which propagation is prohibited for both electromagnetic and acoustic guided wave packets.
The so-called Hartman effect [
17] has been theoretically analyzed [
18] and experimentally measured [
16,
19] for which the crossing time of a quantum potential barrier [
20] initially increases with the dimensions of the barrier and then rapidly reaches a constant value, which is of the order of magnitude of the inverse of the carrier frequency of the wave packet, independent of the dimensions of the barrier. It is immediate that, by simply defining the velocity of crossing the barrier as the ratio between its dimensions and the time taken to cross it, the velocity grows with the dimensions of the barrier when the time reaches a constant value, until it reaches superluminal values.
It is necessary to point out, however, in this phenomenon that the time considered is that of detecting the half height of the rise of the wave front, therefore we are talking about observation of the packet envelope which should correspond to the group velocity, on the other hand the superluminal velocity defined in this way is not a real velocity that can be defined at every point, since, inside the barrier where the waves are evanescent, there is no defined propagation and corresponding velocity. Other possible definitions of the tunneling crossing time and the corresponding velocity are possible although the interpretations are not univocal [
21,
22].
2. Analogy between the Quantum Tunneling Effect and Evanescent, Electromagnetic or Acoustic Guided Waves
Recall that the Schrödinger equation leads to a negative kinetic energy in the case of tunneling since the potential energy
V is greater than the total energy of the particle
E: if we write the equation for a one-dimensional barrier
Its monochromatic solution leads to evanescent modes within the barrier (where V>E) of the type
where
q is such that, having defined the free de Broglie wavenumber as
, that becomes
if a potential
is present, is defined as
inside the barrier wher
.
So, defining for the barrier a characteristic
, then it is possible to write down the Schrödinger equation as
The Schrödinger equation is then completely analogous to the Helmholtz equation for guided waves. This equation is applicable to both guided electromagnetic waves and guided acoustic waves, where represents the electric field inside the waveguide in the former case and the polarized displacement in the horizontal transverse direction, relative to the typical waveguide axis, of the so-called SH (Shear Horizontal) in the latter case of acoustic guided waves.
For guided waves, k is then the wavenumber of the bulk waves of the fields in the material medium considered while is the cut-off wavenumber of the guide such that where is the order of the possible propagation modes inside a guide of thickness b.
In a wave guide it is therefore possible to simulate a potential barrier by narrowing the guide along a section d in such a way that, there, and therefore the propagation wave number in the direction of the guide axis , becomes imaginary giving rise to evanescent waves inside the narrowed part of the guide. For the electromagnetic waves where v is the bulk velocity of the wave in the material while for the ultrasonic SH waves with the velocity of the transverse bulk waves in the medium considered.
The quantum tunneling effect happens when the particle, represented by a wave packet of the probability density wave function with central frequency corresponding to the De Broglie wavelength, passes through an area in which a forbidden potential barrier is present, therefore it has a precise analogy with the propagation of a packet of electromagnetic or acoustic waves within a tapered guide in such a way that the waves are evanescents in the narrowed, limited area representing the barrier.
The equations describing the quantum tunneling effect can be transformed in the ones describing the propagation, in a waveguide, trough a forbidden barrier with the substitution
where
v is the bulk electromagnetic wave velocity or the bulk ultrasonic wave transverse velocity in the material that constitutes the waveguide.
3. A Possible Definition of Traversal Time: The Phase Time
Let’s consider the solution of a one-dimensional Schrödinger equation for a given value of energy
E, in presence of a potential barrier
V. The solution can be expressed in three parts,
,
and
respectively the wave function, on the left, inside and on the rigth of the barrier:
Considering a rectangular shaped potential barrier placed along the
x axes between positions
and
, the boundary conditions of continuity of the wave function and of its derivative leads to a system for the wave amplitudes to solve:
Eliminating the amplitudes
B,
C e
D, it is possible to find the complex transmission coefficient
such that
and an analogous expression for the reflection coefficient
.
It is given a wavepacket strickly picked around a value
k, inpinging on the barrier; its Fourier components will be of the form
exp
. The transmitted wavepacket will be thus described by the expression
where, for a particle,
.
The phase time is the time associated with a recognizable element of the packet such as, for example, the peak to which we assign a position
. To follow the peak, we can use the stationary phase method: we rewrite the exponential in the integral as
where
represents the total changing of phase of the individual components, with a phase shift due to propagation with the path inside the barrier, included. The method, therefore, consists in neglecting the contribution to the integral in the regions in which
varies rapidly with
k, since the rapid oscillations of the function tend to give a null contribution to the integral. The main contribution to the integral comes instead from the regions around the extremal points of
, those where the first derivative with respect to
k vanishes [
14]. This leads to the equation
thus the process of tunneling through a barrier of length
d leads to a spatial delay
and the traversal phase time
is defined as the ratio of the spatial delay due to the overall phase shift to the group velocity
[
23].
From equation (
7)
so the phase displacement
The phase time for the barrier crossing is then
The time defined in this way is certainly an important reference concept, but it is necessary to mention why this definition does not conclude the theoretical discussion, indeed other definitions of traversal time have been hypothesized and the problem is still open.
3.1. Problems with Phase Time
Let’s consider the reflection phase time
(defined in the same way as
but now for the reflected particle) and the dwell time
(Dwell time) [
24] defined as the average time spent by the particle in a certain spatial interval
, which may possibly include the expression barrier having the expression
that is the ratio between the number of the particles in the interval that includes the barrier and the flux of the incident probability density
.
Then, necessarily, the condition
where
R and
T are respectively the reflexion and transmission coefficients, should be satisfied. In our case, instead, there is an additional sinusoidal term [
14,
21,
25]:
where
is the phase displacement of the reflected wave. This term represents the interference between the reflected wave function and the one incident on the barrier. Since, to define the velocity of crossing the barrier we need the actual traversal time taken by the particle (or by the wave packet) to cross it, the phase transmission time
as defined above, cannot represent the actual traversal time because it contains an additional interference term in which the contributions of transmission and reflection are not separable.
4. Modes in Ultrasonic Waveguides
It’s given an ultrasonic rectangular waveguide, characterized by an ideally infinite length (direction z) along which the waves propagate, an ideally infinite width (direction y) and a finite thickness b (direction x) comparable with the wavelength. The normal stress free condition at the surface of the waveguide rules the propagation of the possible modes. With respect to the waveguide it is possible to distinguish three polarization of the ultrasonic waves. The longitudinal one L and, of the two shear, the one polarized in the vertical direction () along the thickness and the other polarized in the horizontal direction () along the width of the waveguide.
When a shear horizontal wave is reflected at the interface constitute by the waveguide surfaces, the boundary conditions are fully satisfied by a reflection in another wave, while the shear vertical is partially reflected and partially transformed in a longitudinal L and so the longitudinal generates at reflection both the and L kind. So two kind of modes can propagate along the waveguide, the pure shear modes, horizontally polarized, and modes that are a combination of longitudinal ad shear vertical polarization called Lamb modes.
Let’s focus on the
modes that have a dispersion relation that is the analogous of the electromagnetic and quantum cases and has an analytic expression for the function
while the dispersion relation for Lamb waves must be solved numerically to obtain the function
. From a mathematical point of view it is possible to define a displacement vector
, a scalar potential
and a vector potential
such that
Through Christoffel equation for isotropic medium [
26] it is possible to demonstrate that the shear waves depend only on
with a wave equation
where
is the shear bulk velocity
and
is the shear modulus [
27].
In the rectangular waveguide the
waves reflect on itself back and forth between the sufaces. Their shear wavenumber vector of modulus
, has transversal component
along
y direction and component
along the direction of propagation
z.
waves have the only displacement component along the horizontal direction
x so the vector potential
has the only component along the vertical
y direction such that
and
The boundary conditions impose the component of the stress, normal to the surfaces
to be null at
. This leads to
and to the ondition of
transverse resonance given by
So a number
n of
modes exist that satisfy the boundary conditions and that have a wavenumber of propagation
such that
This dispersion relation is the analogous of that of the electromagnetic waveguide. Excluding the mode, that is simply a shear bulk wave propagating parallel to the z direction, all the other modes have a cut-off frequency below which, the mode is evanescent.
5. Phase Time for Ultrasonic Waves
The transmission phase time
, as defined in section (
Section 3), is an important reference concept that can be defined also in the acoustic field. Let’s consider its expression for a potential barrier simulated by a zone in which the waves become evanescent (e.g. an area with less thickness than a waveguide such that the frequency is under the cut-off limit) for
waves. A waveguide of thickness
b is considered in which the propagation of a particular
mode is possible with dispersion dispersion equation (
21)
where the wave vector
, of the propagating wave along the waveguide, is real and the propagating
mode has expression
In the narrowed section of the waveguide of length
d, if the thickness
is such that
, the correspondent
becomes imaginary and it is substituted by
and the wave becomes evanescent with expression
Defining, then,
the traversal phase time for
waves has an expression analogous to (
12)
where the group velocity is
.
6. Results
6.1. Theorical Results for SH Waves
In a concrete example, we can consider a waveguide with thickness
mm tapered in an area with a smaller thickness
mm, of length
d; the velocity of the transverse waves is set to
m/s typical of aluminum. In this case it is found that the cut-off frequency of the evanescent wave zone is
MHz. In
Figure 1 the phase time is graphed as a function of frequency for the first mode
and a barrier
mm long. To be noted the barrier resonances above the cutoff frequency and the monotonic decrease in traversal time below the cutoff frequency.
If we then consider, for example, a frequency
MHz (wavelength
mm) just below the cut-off frequency
MHz in
Table 1 are indicated the values of the transversal phase time
and the corresponding traversal velocity
defined as the ratio between the length of the barrier (evanescent zone) and the traversal phase time, corresponding to different lengths
d of the barrier. It is possible to note the so called Hartman effect for which the traversal time tends to a constant limit value, independent by the barrier length, resulting in rapid increasing in traversal velocity; the effect happens when, increasing the barrier length
d, it becomes
opaque i.e.
. Remembering that the velocity of bulk transverse waves
m/s represents the analogue of the velocity of light in the medium that constitutes the electromagnetic wave guide, it is possible to note that, with these parameters, the acoustic analogous of the apparent superluminal behavior, is already reached with a barrier length corresponding to few wavelenghts of the signal.
At lower frequencies then, i.e. further inside the potential well, the Hartman effect is reached even for shorter barrier lengths and the value of the traversal time limit decreases slightly and it is within the order of magnitude of a fraction of a microsecond.
This is graphed in
Figure 2 where the phase time is in function of the barrier length for several frequencies under the cut-off. The plateau of constant phase time, showing the Hartman effect, is reached, for lower frequencies at shorter barrier lengths.
Defining the traversal velocity as the ratio between barrier length and phase time increment during the crossing, in
Figure 3 it is possible to see at which barrier length, the value of this velocity goes above the limit velocity of shear bulk waves in the material, for different frequencies. The lower the frequency, the shorter barrier length is necessary to cross above this limit.
6.2. Some Experimental Results in Literature
Traversal time measurements in the tunneling effect in quantum mechanics pose several problems for issues related to non-invasive measurements of quantum objects [
20]. In electromagnetism and acoustics, measurements of the delay of the packet envelope were instead carried out for different frequencies and in different situations in which evanescent waves were involved.
Table 2 lists some experimental results in photonics (tunneling between two FTIR total internal reflection prisms, photonic chain of dielectric layers, tapered waveguide), in guided microwaves and in guided Lamb waves.
Table 2 shows experimental results that support a hypothesis made by Nimtz and Stahlhofen [
35] for which the traversal time is of the order of magnitude of the reciprocal of the frequency used. This is found for both electromagnetic waves and acoustic waves. This fact, combined with the Hartman effect of independence of the traversal time from the length of the barrier for "opaque" barriers (beyond a minimum length) leads to the possibility of an envelope of the wave packet traveling at very large velocities, even superluminal for electromagnetic waves and beyond the maximum elastic velocity of the considered medium for elastic waves.
We remember that in any case in particular conditions, the velocity of the envelope does not correspond to the velocity of information transfer, therefore the causality principle of special relativity remains intact.
7. Conclusions
Following the formal analogy among quantum, electromagnetic and acoustic waves, the problem of determination of the traversal time of a potential barrier by evanescent acoustic waves has been addressed. The analogy has been demonstrated for the (shear horizontal) modes in a rectangular waveguide that have the same analytical form of dispersion relation as quantum and electromagnetic waves; then with an appropriate tapering of the waveguide a potential barrier has been simulated. The article focused in particular on the determination of the so called phase time i.e. the temporal delay due to a phase displacement for a wavepacket crossing a region of forbidden propagation, putting in evidence the so called Hartman effect also for acoustic modes.
The Hartman effect happens for opaque barriers when the length d is such that . The typical traversal time, in that case, is found for waves and compared with other results in literature obtained in different fields, confirming the conjecture that the traversal phase time is proportional to the inverse of the typical frequency of the phenomenon independently from the length of the opaque barrier.
This leads to an increasing of traversal velocity defined as barrier length divided by the phase time, well above the velocity limit relative to the phenomenon in object, which, for acoustic shear guided waves, is the shear velocity of bulk waves in the material of which the waveguide is made and that is the analogous of the ligth c velocity in electromagnetic and quantum fields.
Conflicts of Interest
The author declare no conflict of interest.
References
- Lord Rayleigh, Nature (London) XXIV, 1881, 382.
- Lord Rayleigh, Philos. Mag. XLVIII, 1889, 151.
- Brillouin, L. Wave propagation and group velocity, Academic Press, New York, 1960.
- Winful, H.G. Group delay, stored energy, and the tunnelling of evanescent electromagnetic waves. Phys. Rev. E 2003, 68. [Google Scholar]
- Garrett, C.G.B.; McCumber, D.E. Phys. Rev. A 1970, 1, 305. [CrossRef]
- Chu, S.; Wong, S. Linear pulse propagation in an absorbing medium. Pys. Rev. Lett. 1982, 48, 738. [Google Scholar] [CrossRef]
- Akulshin, A.M.; Barreiro, S.; Lezama, A. Pys Rev. Lett. 1999, 83, 4277. [CrossRef]
- Basov, N. G et al. Sov. Phys. JETP 1966, 23, 16. [Google Scholar]
- Chiao, R.Y. Phys Rev. A 1993, 48, R34. [CrossRef]
- Steinberg, A.M.; Chiao, R.Y. Phys. Rev. A 1994, 49, 2071. [CrossRef]
- Li-Gang, W.; Nian-Hua, L.; Qiang, L.; Shi-Yao, Z. Effect of coherence on the superluminal propagation of light pulses through anomalously dispersive media with gain. Europhys. Lett. 2002, 60, 834. [Google Scholar]
- Fisher, D.L.; Tajima, T.; Downer, M.C.; Siders, C.W. Phys. Rev. E 1995, 51, 4860.
- Steinberg, A. M.; Kwiat, P.G.; Chiao, R.Y. Measurement of the single-photon tunnelling time. Phys. Rev. Lett. 1993, 71, 708–711. [Google Scholar] [CrossRef]
- Hauge, E.H.; Støvneng, J.A. Tunneling times: a critical review. Rev. Mod. Phys. 1989, 61, 917–936. [Google Scholar] [CrossRef]
- Ranfagni, A.; Mugnai, D.; Fabeni, P.; Pazzi, G.P. Delay-time measurements in narrowed waveguides as a test of tunnelling. Appl. Phys. Lett. 1991, 58, 774. [Google Scholar] [CrossRef]
- Enders, A.; Nimtz, G. Evanescent-mode propagation and quantum tunneling. Phys. Rev. E 1993, 48, 632–634. [Google Scholar] [CrossRef] [PubMed]
- Hartman, T.H. Tunneling of a wave packet. J. Appl. Phys. 1962, 33, 3427. [Google Scholar] [CrossRef]
- Ghatak, A.; Banerjee, S. Appl. Opt. 1989, 28, 1960.
- Spielmann, Ch.; Szipöcs, R.; Stingl, A.; Krausz, F. Tunneling of optical pulses through photonic band gaps. Phys. Rev. Lett. 1994, 73, 2308. [Google Scholar] [CrossRef]
- Muga, J.G.; Leavens, C.R. Arrival time in quantum mechanics. Phys. Reports 2000, 338, 353–438. [Google Scholar] [CrossRef]
- Büttiker, M. Larmor precession and the traversal time for tunnelling. Phys. Rev. B 1983, 27, 6178–6188. [Google Scholar] [CrossRef]
- Winful, H.G. Tunneling time, the Hartman effect and superluminality: a proposed resolution of an old paradox. Physics Reports 2006, 436, 1–69. [Google Scholar] [CrossRef]
- Landauer, R.; Martin, Th. Barrier interaction time in tunnelling. Rev. of Mod. Phys. 1994, 66, 217–227. [Google Scholar] [CrossRef]
- Leavens, C.R.; Aers, G.C. Dwell time and phase times for transmission and reflection. Phys. Rev. B 1989, 39, 1202–1206. [Google Scholar] [CrossRef] [PubMed]
- Landauer, R.; Martin, Th. Time delay in wave packet tunnelling. Solid State Commun. 1992, 84, 115–117. [Google Scholar] [CrossRef]
- Slawinski, M, A. Waves and rays in elastic continua, World Scientific, 2010.
- Salencon J. Handbook of continuum mechanics: general concepts, thermoelasticity, springer science & business media, 2001.
- Balcou, Ph.; Dutriaux, L. Phys. Appl. Rev. Lett. 1997, 78, 851. [CrossRef]
- Mugnai, D.; Ranfagni, A.; Ronchi, L. The queston of tunneling time duration: a new experimental test at microwave scale. Phys, Lett. A 1998, 247, 281. [Google Scholar] [CrossRef]
- Carey, J.J.; Zawadzka, J.; Jaroszynski, D.; Wynne, K. Phys. Rev. Lett. 2000, 84, 1431. [CrossRef] [PubMed]
- Haibel, A.; Nimtz, G. Universal tunnelling time in photonic barrier. Ann. Phys. 2001, 10 Ed. 8, 707–712. [Google Scholar] [CrossRef]
- Nimtz, G.; Enders, A.; Spieker, H. J. Phys. I. 1994, 4, 565.
- Enders, A.; Nimtz, G. J. Phys. I. 1992, 2, 1693.
- Alippi, A.; Germano, M.; Bettucci, A.; Farrelly, F.A.; Muzio, G. Traversal time of acoustic plate waves through a tunneling section. Phys. Rev. E 1998, 57, 4907–4910. [Google Scholar] [CrossRef]
- Nimtz, G.; Stahlhofen, A. A. Universal tunneling time for all fields. Ann. Phys. 2008, no.6, 17, 374–379. [Google Scholar] [CrossRef]
|
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).