Abstract
We analyze different examples to show that the socalled generalized Hartman effect is an erroneous presumption. The results obtained for electron tunneling and transmission of electromagnetic waves through superlattices and Bragg gratings show clearly the resonant character of the phase time behavior so that a generalized Hartman effect is not expected to occur. A reinterpretation of the experimental results in double Bragg gratings is proposed.
Background
The Hartman [1] effect is known as the independence of the tunneling time on the barrier width as this parameter gets large. It has been shown that the experimental evidences of this effect on the transmission times of photons and electromagnetic pulses [26] are compatible with phase time calculations [7]. The Hartman effect has been investigated in various ways by extending the system not only for a single barrier but also for double [8,9] and multiple barrier [10,11] structures. Olkhovsky, Recami, and Salesi came out with the idea that for nonresonant tunneling through two potential barriers, the tunneling time (which is a phase time) is independent not only of the barrier width but also of the barrier separation [8]. The approximations introduced in this reference to obtain the unknown coefficients, led these authors to unphysical results like the generalized Hartman effect. This has been called the generalized Hartman effect (GHE). The twobarrier problem can be solved without approximations, see for example, in the work of Pereyra [12]. An experiment to check this effect was performed by Longhi et al. [10] where optical pulses of 1,550 nm wavelength were transmitted through a doublebarrier system of Bragg gratings. In this reference, nonconclusive and inappropriately presented results for five different values of the gratings separation were reported. Most of the theoretical conclusions were based on questionable formulas and unnecessarily involved calculations. For example, Equation 2 (used in Equations 3 and 4) of [8] is not the actual transmission coefficient through a double Bragg grating. A criticism on the mathematical rigor on GHE is also given by S. Kudaka and S. Matsumoto [13,14]. It is easy to check from a straightforward calculation, or from the precise and general formulas published in [7] as quoted below, that the phase time for a double barrier (DB) with separation L has the structure
with T_{2} and T the double and singlebarrier transmission coefficients, respectively, k the wave number, ω the frequency and A_{i}, A_{r}, F, and G simple functions of the potential parameters (P. Pereyra and H. P. Simanjuntak, unpublished work). Despite this clear dependence on L, involved and contradictory arguments lead to establish that τ becomes independent of L[8,10,11]. In the following we will consistently use a for the separation between barriers.
For the inference of a generalized Hartman effect to be meaningful for multibarriers, double superlattices (SLs) or double Bragg gratings (BG), one would of course need to keep the physical parameters [like the energy (wavelength) of the particle (wave)] fixed as the length between barriers is increased. The tunneling and transmission times behavior should be taken with care when one tries to find a Hartman effect due to barrier separation in multibarrier systems [8,11] since, in general, the density of resonance energies grows rapidly as the separation increases. It is well known that the nonresonant gaps in the band structure of a SL or a BG become resonating when these systems are divided and separated; and the separation is increasingly varied. This was already recognized in [15] (for double SL) and in [10] (for double BG). On the other hand, it is well known that the tunneling time follows the resonant band structure [7,16]. Thus, it is not possible to keep increasing the separation between barriers and superlattices without crossing resonances. For this reason, visualized here with specific examples for electrons and electromagnetic waves, the existence of a generalized Hartman effect is a rather questionable issue. For these examples we perform first principle calculations using the actual transmission coefficient of the system (such as that of double BG in the experiment in [10]) so that we can justify completely that the socalled generalized Hartman effect is erroneous.
To study the Hartman effect and to criticize the presumption of a generalized Hartman effect in superlattices, Bragg gratings, and multibarrier systems, we will use the theory of finite periodic system that allows straightforward calculation of the phase time. For electron tunneling, we shall assume periodic and sectionally constant potentials with cells of length ℓ_{c}=a+b and a barrier of width b and strength V_{o} in the middle. For electromagnetic waves, each cell consisting of dielectrics 1 and 2 will contain a dielectric 2 of length b in the middle. In this case ϵ_{i}, n_{i}, and μ_{i} (with i=1,2) are the corresponding permittivities, refractive indices, and permeabilities; the regions outside the SL are assumed to be air. For Bragg gratings, the refractive indices are periodic.
Methods
If we have a Gaussian wave packet (of electrons or electromagnetic waves) through a SL of length nℓ_{c}−a, the centroid phase time (which is taken here as the tunneling or transmission time) is given by [7,17,18]
Here α=α_{R}+iα_{I} is the (1,1) element of the singlecell transfer matrix M; U_{n}(α_{R}) are the Chebyshev polynomials of the second kind evaluated at α_{R}; and α_{n} is the (1,1) element of the ncell transfer matrix M_{n}. This is given by [16]
At resonance, where U_{n−1}=U_{2n−1}=0, we have [16]
The expression for the tunneling or transmission time simplifies as
The tunneling time in Equation 2 is exact and general and valid for arbitrary number of cells, barrier width, and barrier separation. Thus, one can check the existence or not of a (generalized) Hartman effect at will. For concrete examples, we consider superlattices like (GaAs/Al_{0.3}Ga_{0.7}As)^{n}/GaAs, with electron effective mass m_{A}=0.067 m in GaAs layers, m_{B}=0.1 m in Al_{0.3}Ga_{0.7}As layers (m is the bare electron mass) and V_{o}=0.23 eV, and Bragg gratings with periodic refractive index.
Results and discussion
Electron tunneling
If we consider electrons through superlattices with unit cell length ℓ_{c}=a+b, we will have
with and . When m_{A}, m_{B} and V_{o} are taken as fixed parameters, we choose a=100 Å and b=30 Å.
For a single barrier, n=1, the tunneling time τ_{1} plotted in Figure 1 as a function of the reduced barrier width b/λ shows the wellknown Hartman effect. The energy E is kept fixed and is the de Broglie wavelength.
Figure 1. The singlebarrier transmission coefficient and the tunneling time. The singlebarrier transmission coefficient 1/α^{2} (gray lines) and the tunneling time τ_{1} (dark lines) as functions of the reduced barrier width b/λ, when the electron energies are E=0.122516 eV, E=0.15 eV and E=0.2 eV. In the tunneling time curves, the Hartman effect is evident.
With α_{R} and α_{I} growing exponentially with the barrier width b, one can easily show from Equation 2 that for large b, the nonresonant tunneling time approaches that for a single barrier, i.e., τ_{n}(E)≈τ_{1}(E) as
This is the wellknown Hartman effect. Since this quantity becomes also independent of the barrier separation [8,11]a, it has been taken as the analytical evidence of a generalized Hartman effect. However, such an approximation that leads to the independence on a and n is obtained by taking the limit of large b first that is strictly speaking infinite, which makes the first barrier the only one that matters for the incoming wave to penetrate while the rest of the SL is immaterial. This was also pointed out by Winful [9]. However, Winful [9] used an approximation: The transmission of the double square barrier potential to model the transmission through the double BG. Here, we present calculations using the actual transmission coefficient through the double BG. As mentioned before, for the generalized Hartman effect to be meaningful, it should not matter whatever limit we take first whether on a, b, or n. It turns out that a nonresonant energy region becomes resonant as the separation a increases (see the discussion on the double Bragg gratings in section ‘Hartman effect in two Bragg gratings systems’).
The situation is completely different for resonant tunneling through a SL with large but finite barrier width b where Equation 5 shows that the tunneling time becomes τ_{n}(E)∝be^{2qb} (since α_{R} and α_{I} behave as e^{qb} for large b). Thus, relatively small barrier width would be needed to study the effect of the barrier separation and the number of barriers on the tunneling time. The tunneling time for a relatively small barrier width is shown in Figure 2 for an electron (with energy E=0.15 eV) through SLs which number of cells are n=3,4, and 6.
Figure 2. The tunneling time τ_{n} as a function of the reduced barrier width. The tunneling time τ_{n} as a function of the reduced barrier width b/λ for electrons (with energy E=0.15 eV) through superlattices with n=3,4, and 6.
Looking at α_{R} and α_{I}, that are oscillating functions in a, it is clear that it is not possible to have the tunneling time to be independent of the barrier separation a, by keeping the barrier width and number of cells fixed. Therefore, the socalled generalized Hartman effect is at least dubious. The tunneling time behavior that will be found below for the double BG is easy to understand here. Starting with a certain barrier separation a, a nonresonant phase time becomes a resonant one as a is increased, while the other parameters are kept fixed. This is shown in Figure 3 where the tunneling time is plotted as a function of the reduced barrier separation, a/λ, for fixed b, n, and electron energy E. This result shows that in this kind of systems, the presumption of a generalized Hartman effect is incorrect.
Figure 3. The tunneling time τ_{6} as a function of reduced barrier separation and fixed barrier width. The tunneling time τ_{6} as a function of reduced barrier separation a/λ for fixed barrier width b, number of cells n=6 and electron energy E=0.15 eV with the corresponding de Broglie wavelength λ.
The Hartman effect as a consequence of varying the number of cells was already discussed in [7]. In Figure 4 we show three qualitatively different examples on the behavior of the tunneling time as a function of n. In Figure 4a for energies in the gap (E=0.15 eV and E=0.2 eV), the saturation of the tunneling time exhibits the wellknown Hartman effect. In Figure 4b, the energy lies at the edge of a resonant region. The phase time τ_{n} resonates for multiples of n=21. This behavior is clearly understood if we consider Equations 4 and 5. Equation 4 implies that the same resonance energy is found for different number of cells as long as the ratio ν/n is constant. This means that . From Equation 5, it is also evident the linear dependence of τ_{n} on n.
Figure 4. The tunneling time τ_{n} as the number of cellsn in a SL is varied. (a) Saturation of τ_{n} for electron energies E=0.15 eV and E=0.2 eV in the gap. (b) The energy is close to a resonant bandedge. In this case, more resonances appear as n is increased with the energy fixed. No Hartman effect can be inferred from this figure.
The Hartman effect and the electromagnetic waves
Electromagnetic waves have been used for discussions on the Hartman effect [9]. For a superlattice L(H/L)^{n} made of alternating layers with refractive indices n_{L} and n_{H}, the phase time (PT) for each frequency component of a Gaussian wave packet through a SL of length nℓ_{c}−a is also obtained from Equation 2 with k_{L,H}=ωn_{L,H}/c and with [7]
To see the effect of varying the size of the SL on the PT, one has to be sure that such variation will still keep the wavelength inside a photonic band gap. It was shown that by increasing the number of cells, for fixed thicknesses of layers and wavelength in a gap, the PT exhibits [7] the observed Hartman effect [2,3]. However, this condition will not be possible by varying arbitrarily the thicknesses of the layers. The reason is that there is only a small range of thicknesses that one can use to keep the chosen wavelength to lie in a gap before going out of it and may even reach resonances, as shown in Figure 5 where the PT oscillates (with a band structure) and grows as a function of the reduced thicknesses a/λ and b/λ. This is analogous to the electron tunneling time shown in Figure 3.
Figure 5. The phase times τ_{n} as functions of the reduced thicknesses. The phase times τ_{n} as functions of the reduced thicknesses a/λ (a) and b/λ (b) for a superlattice (L/H/L)^{n} made of fused silica (L) (n_{L}=1.41) and titanium oxide (H) (n_{H}=2.22) for λ=702 nm.The oscillating behavior of τ_{n} with clear band structure is far away of any Hartman effect.
Hartman effect in two Bragg gratings systems
We now consider the system that was taken in [10] as thought to support the idea of a generalized Hartman effect: the double Bragg gratings (DBG). Independent of the approximate method used in that paper, we find that assuming sin(k_{B}a)=0 (the only way to obtain the reduced expressions of Table 1 in [10]) and still considering a as a variable are incongruous. Moreover, the idea that the PT becomes independent of a is incompatible with the Equation (4b) in their work, where a linear dependence on a is reported. In the DBG, the gratings of length L_{o} and refractive index n(z)=n_{0}+n_{1} cos(2k_{B}z) are separated by a distance a. The values of a considered in the experiment are indicated by arrows in Figure 6. The BG wave equation
Figure 6. The phase time as a function of the Bragg gratings separation. (a) The phase time as a function of the separation a between two Bragg gratings, for incident λ=1,542 nm, k_{B}=6.1074/μm, n_{0}=1.452, n_{1}/n_{0}=1.8×10^{−4}, and L_{o}=8.5 mm. (b, c) The PT is plotted as a function of ω, for a=42 mm. The phase time in (b) is the same as that in (c) but plotted from 0 to 10 ns to compare with Figure 2 in [10]. Arrows indicate the as in [10].
when ignoring the (n_{1}/n_{0})^{2} term for n_{1}/n_{0}≪1 (as in [10]), becomes the Mathieu equation, in which solutions ψ_{1}(z)=Se(u,v;k_{B}z+Π/2) and ψ_{2}(z)=So(u,v;k_{B}z+Π/2) are Mathieu functions [19] with and . The real and imaginary parts of the (1,1) element of the transfer matrix are
with W the Wronskian and
Here θ_{1}=θ(L_{o},0), θ_{2}=θ(2L_{o}+a,L_{o}+a) analogously for χ_{1,2}, μ_{1,2}, ν_{1,2}, with
Using parameters of Longhi et al. [10] for n_{0},n_{1}, k_{B}, and L_{o}, the nonresonant gap becomes resonant as the gratings separation increases. Though details are beyond the purpose of this paper, we plot in Figure 6 the PT as a function of the separation a for incidentfield wavelength λ=1542 nm, and as a function of the frequency ω, for a=42 mm. Recall that in [10], λ≃1,550 nm was considered. While the PT appears completely in graph (c), in (b) it is plotted in a different range to compare with the experiment. The resonant behavior of the PT with a and the absence of any generalized Hartman effect are evident. Similar results are obtained when λ=2Π/k_{B}.
Conclusion
We have shown that the presumption of generalized Hartman effect for tunneling of particles and transmission of electromagnetic waves is not correct.
Competing interests
The authors declare that they have no competing interest.
Authors’ contributions
HPS initialized the work. Both authors carried out the calculations. Both authors read and approved the final manuscript.
Acknowledgements
The authors would like to thank Professor Norman H. March for comments and suggestions on the manuscript.
References

Hartman TE: Tunneling of a wave packet.
J Appl Phys 1962, 33:3427. Publisher Full Text

Steinberg AM, Kwiat PG, Chiao RY: Measurement of the singlephoton tunneling time.
Phys Rev Lett 1993, 71:708. PubMed Abstract  Publisher Full Text

Spielmann Ch, Szipöcs R, Stingl A, Krausz F: Tunneling of optical pulses through photonic band gaps.
Phys Rev Lett 1994, 73:2308. PubMed Abstract  Publisher Full Text

Ranfagni A, Mugnai D, Fabeni P, Pazzi GP: Delaytime measurements in narrowed waveguides as a test of tunneling.

Leavens CR, Aers GC: Dwell time and phase times for transmission and reflection.
Phys Rev B 1989, 39:1202. Publisher Full Text

Pereyra P: Closed formulas for tunneling time in superlattices.
Phys Rev Lett 2000, 84:1772. PubMed Abstract  Publisher Full Text

Olkhovsky VS, Recami E, Salesi G: Superluminal tunneling through two successive barriers.
Europhys Lett 2002, 57:8792002. Publisher Full Text

Winful HG: Apparent superluminality and the generalized Hartman effect in doublebarrier tunneling.

Longhi S, Laporta P, Belmonte M, Recami E: Measurement of superluminal optical tunneling times in doublebarrier photonic band gaps.

Pereyra P: Fundamentals of Quantum Physics. Heidelberg: Springer Verlag; 2012.

Kudaka S, Matsumoto S: Questions concerning the generalized Hartman effect.
Phys Lett A 2011, 375:3259. Publisher Full Text

Kudaka S, Matsumoto S: Reply to Comments on: ‘questions concerning the generalized Hartman effect’.
Phys Lett A 2012, 376:1403. Publisher Full Text

Nimtz G, Haibel A, Vetter RM: Pulse reflection by photonic barriers.

Pereyra P, Castillo E: Theory of finite periodic systems: general expressions and various simple and illustrative examples.

Simanjuntak HP, Pereyra P: Evolution and tunneling time of electron wave packets through a superlattice.

Pereyra P, Simanjuntak HP: Time evolution of electromagnetic wave packets through superlattices: evidence for superluminal velocities.

Morse PM, Feshbach V: Methods of Theoretical Physics. Part I. New York: McGraw Hill; 1953.