On the generalized Hartman effect presumption in semiconductors and photonic structures

We analyze different examples to show that the so-called 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 [2][3][4][5][6] 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 non-resonant 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 *Correspondence: pereyrapedro@gmail.com 2 Departamento de Ciencias Básicas, UAM-Azcapotzalco, C.P. 02200, Mexico, Federal District, Mexico 3 Institute für Experimental and Applied Physics, University of Regensburg, Regensburg, Bavaria, Germany Full list of author information is available at the end of the article transmitted through a double-barrier system of Bragg gratings. In this reference, non-conclusive 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 single-barrier 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. http://www.nanoscalereslett.com/content/8/1/145 For the inference of a generalized Hartman effect to be meaningful for multi-barriers, 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 non-resonant 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 so-called 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 multi-barrier 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 single-cell 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 n-cell 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

Electron tunneling
If we consider electrons through superlattices with unit cell length c = a + b, we will have 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 well-known Hartman effect. The energy E is kept fixed and λ = 2π / √ 2m A E is the de Broglie wavelength.
With α R and α I growing exponentially with the barrier width b, one can easily show from Equation 2 that for large b, the non-resonant tunneling time approaches that for a single barrier, i.e., τ n (E) ≈ τ 1 (E) as This is the well-known 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 http://www.nanoscalereslett.com/content/8/1/145 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 non-resonant 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.
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 so-called 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 non-resonant 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.
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 well-known 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 E n ν is found for different number of cells as long as the ratio ν/n is constant. This means that E n 1 = E 2n 2 = E 3n 3 = . . . . From Equation 5, it is also evident the linear dependence of τ n on n.

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 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.

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 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 u = 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].
(13) http://www.nanoscalereslett.com/content/8/1/145 Using parameters of Longhi et al. [10] for n 0 , n 1 , k B , and L o , the non-resonant 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 incident-field 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.