SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Nano Express

Tunneling times of acoustic phonon packets through a distributed Bragg reflector

Zorayda Lazcano1, Pedro Luis Valdés Negrín2, Diosdado Villegas12, Jesus Arriaga1* and Rolando Pérez-Álvarez3

Author Affiliations

1 Instituto de Física, Benemérita Universidad Autónoma de Puebla, Apartado Postal J-48, CP 7250 Puebla, México

2 Departamento de Física, Universidad Central ‘Marta Abreu’ de Las Villas, CP 54830 Santa Clara, Cuba

3 Universidad Autónoma del Estado de Morelos, Ave. Universidad 1001, CP 62209 Cuernavaca, México

For all author emails, please log on.

Nanoscale Research Letters 2014, 9:449  doi:10.1186/1556-276X-9-449

The electronic version of this article is the complete one and can be found online at: http://www.nanoscalereslett.com/content/9/1/449

Received:10 May 2014
Accepted:12 August 2014
Published:29 August 2014

© 2014 Lazcano et al.; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.


The longwave phenomenological model is used to make simple and precise calculations of various physical quantities such as the vibrational energy density, the vibrational energy, the relative mechanical displacement, and the one-dimensional stress tensor of a porous silicon distributed Bragg reflector. From general principles such as invariance under time reversal, invariance under space reflection, and conservation of energy density flux, the equivalence of the tunneling times for both transmission and reflection is demonstrated. Here, we study the tunneling times of acoustic phonon packets through a distributed Bragg reflector in porous silicon multilayer structures, and we report the possibility that a phenomenon called Hartman effect appears in these structures.

Distributed Bragg reflector; Tunneling time; Hartman effect


Phonons, the quanta of lattice vibrations, manifest themselves practically in all electrical, thermal, and optical phenomena in semiconductors and other material systems. The reduction of the size of electronic devices below the acoustic phonon mean free path creates a new situation for phonon propagation and interaction, opening up an exciting opportunity for engineering phonon spectrum in nanostructured materials [1]. Since the early work of Narayanamurti et al. [2], important progress has lately emerged in the development of nanowave phononic devices including, e.g., mirrors, cavities, and monochromatic sources.

How long does it take for a particle to tunnel through a potential barrier? This is a question that has occupied physicists for decades and one for which there is still no definitive answer [3]. The Hartman effect (HE) states that the tunneling time becomes independent of the barrier length [4]. The independence of tunneling time on barrier length would imply arbitrarily large and indeed superluminal velocities for tunneling wave packets, if tunneling was in fact a propagation phenomenon.

Phonon tunneling studies have also revealed phenomena related to the HE. Recent experiments on tunneling acoustic waves have reported the breaking sound barrier [5,6]. Yang et al. found that inside a phononic band, the group velocity increases linearly with the sample thickness, a rather remarkable effect that is a signature of tunneling in quantum mechanics [6,7]. Villegas et al. have discussed the physical conditions under which the tunneling time for long-wavelength phonons through semiconductor heterostructures is independent on the system’s size, i.e., the effect equivalent to the HE for electrons [8]. Experimental studies of HE in the context of the nanophononics have been carried out [9,10]. In these works, very short transit times in the stop bands have been measured, one acoustic equivalent of HE of electron tunneling through potential barriers.

During the last decade, interest in achieving all-silicon-based opto- and microelectronics was highly stimulated by the discovery of the unique optical properties of porous silicon [11]. Porous silicon is known as a versatile material with applications in light emission, sensing, and photonic crystal devices. It is well known that the introduction of artificial spatial periodicity in the elastic properties of a system results in Brillouin zone folding. Such folding is often accompanied by the appearance of bandgaps in the phonon frequency spectrum [12,13]. In the last few years, this interest has been translated to porous silicon-based phononic systems [14-17]. Here, we study the tunneling times of acoustic phonon packets through a distributed Bragg reflector in porous silicon multilayer structures.

The paper is organized as follows: The ‘Methods’ section provides some fundamentals in both the long-wavelength model and the transfer matrix method. The main theoretical findings are presented in the ‘Results and discussion’ section. In particular in this section, from general principles such as invariance under time reversal, invariance under space reflection, and conservation of energy density flux, the equivalence of the tunneling times for both transmission and reflection is demonstrated. At the end of the paper, the main conclusions are given.


Long-wavelength model

The one-dimensional energy density [18] is defined as

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M1">View MathML</a>


where the first term in (1) represents the kinetic energy density, the second one represents the interaction energy density of the phonon field with itself, and the third one represents the strain energy density that accounts for the dispersive character of the oscillations. These terms depend on the atomic relative displacements u, the linear mass density ρ, the phonon frequency at the center of the Brillouin zone ωΓ, the one-dimensional strain tensor u/z, and the stress tensor σ, which is equal to

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M2">View MathML</a>


being β a parameter that accounts on the behavior of the bulk phonon dispersion relation. From Equations 1 and 2, we can obtain the one-dimensional equations of motion,

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M3">View MathML</a>



<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M4">View MathML</a>


The continuity equation for the energy density is defined as

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M5">View MathML</a>


where the energy density flux, j, is given by

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M6">View MathML</a>


By substituting Equation 2 into Equation 6, we obtain the following expression for the energy density flux:

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M7">View MathML</a>


A very convenient method for calculating the propagation of harmonic waves through a system consisting of a finite number of uniform layers is the transfer matrix method [19,20]. The transfer matrix, which propagates the wave amplitudes across a j-th homogeneous layer is given by

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M8">View MathML</a>


where Zj=ρjvj is the acoustic impedance, ρj the mass density, vj the velocity of sound, and dj the layer width.

Results and discussion

Equivalence of the transmission and reflection times in the tunneling of long-wavelength phonons

In this paper, for the reason of complement, we give a demonstration of the equivalence of the transmission and reflection times (τt=τr) in the tunneling of long-wavelength phonons [21]. Consider the general phonon scattering process as shown schematically in Figure 1. Let us consider that the phonon propagation is normal to the layer interfaces and adopt the continuum model valid for long-wavelength oscillations. In this section, we analyze the properties of the transfer matrix in the tunneling process of long-wavelength phonons.

thumbnailFigure 1. General stationary scattering configuration. The general stationary scattering configuration in one dimension. An arbitrary barrier is confined to the interval (−a,a).

Invariance under time reversal

Let us observe that if we take the complex conjugate of (3) and if t is replaced by −t, we get

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M9">View MathML</a>


provided only that β(z) and ωΓ(z) are real functions. Observe that this equation has the same form as (3). Therefore, if u(z,t) is a solution of (3), then u(z,−t) is also a solution. u(z,−t) is often referred as the time-reversed solution. The behavior of the wave equation exhibited by (9) is called invariance under time reversal. For the stationary state, invariance under time reversal implies that if u(z) is a stationary-state wave function, then u(z) is also one.

The general solution of the wave equation for the system depicted in Figure 1 is given by

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M10">View MathML</a>


where <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M11">View MathML</a> is the wave number in the left (right) part of the system. Using the transfer matrix, we can relate the wave amplitudes in the left and right parts according to

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M12">View MathML</a>



<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M13">View MathML</a>


On the other hand, the time-reversed solution of the wave equation is

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M14">View MathML</a>


Hence, comparing Equations 10, 11, and 13, we get

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M15">View MathML</a>



<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M16">View MathML</a>


The properties of the transfer matrix in the tunneling of long-wavelength phonons can be obtained from Equations 11 and 14

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M17">View MathML</a>


Conservation of energy density flux

To demonstrate the conservation of the energy density flux, we need to calculate this quantity in the left and right parts of the system. Using Equations 7 and 10, we obtain

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M18">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M19">View MathML</a>


where the superscript L (R) refers to the left (right) of the region (−a,a), Im denote imaginary part and <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M20">View MathML</a> and <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M21">View MathML</a>. The conservation of the energy density flux is expressed as jL=jR, or equivalently

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M22">View MathML</a>


where Mt and M denote the transpose and conjugate matrix, respectively, and

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M23">View MathML</a>


Using Equation 19 and the condition (16), one further condition on matrix M can be added

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M24">View MathML</a>


Invariance under space reflection

Now, let us consider the invariance under space reflection of Equation 3. Because β2(z), ω2(z), and <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M25">View MathML</a> are even functions of z, another solution of the wave equation is obtained by replacing z by −z. We can immediately write the general solution of the wave equation as

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M26">View MathML</a>


and obtain the following relation

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M27">View MathML</a>


Substituting (23) into (14) gives us the condition

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M28">View MathML</a>


where I is the identity matrix. From this last equation, we obtain

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M29">View MathML</a>


The transmission (reflection) amplitude <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M30">View MathML</a> for a phonon incident perpendicular to the barrier from right to left (the inverse process is described by D=1,C=ArR,A=0,B=AtR) and the transmission (reflection) amplitude <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M31','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M31">View MathML</a> for a phonon incident from left to right (the direct process is described by <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M32','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M32">View MathML</a>, <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M33','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M33">View MathML</a> are given by the following equations:

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M34','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M34">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M35','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M35">View MathML</a>


in terms of the modulus and phase. By applying conditions (16) and (21) to the previous relations, we obtain

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M36">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M37','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M37">View MathML</a>


Finally, the transmission and reflection coefficients are given by

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M38','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M38">View MathML</a>


where t, r, and i denote transmitted, reflected, and incident, respectively. From Equations 16, 27, 28, and 30, it is very easy to prove that <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M39','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M39">View MathML</a> and <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M40','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M40">View MathML</a>. From Equation 28, we obtain αtR(k3,k1)=αtL(k1,k3), and from Equations 16, 26, 27, 28, and 29, we can deduce the relation

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M41">View MathML</a>


If we consider that both materials in the left and right parts of Figure 1 are identical, i.e., (k1=k3k), then αtR(k3,k1)=αtL(k1,k3)≡αt(k) and βrR(k3,k1)=βrL(k1,k3)≡βr(k). By using this relations, the phase is

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M42">View MathML</a>


Finally, using Equation 32 is straightforward to prove that the transmission and reflection phase times given by

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M43','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M43">View MathML</a>



<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M44">View MathML</a>


are equal,

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M45','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M45">View MathML</a>


Falck et al. have obtained a similar result for symmetric scattering potentials [22].

Tunneling of acoustic phonons through a distributed Bragg reflector

In this section, we apply the concepts discussed previously to study a distributed Bragg reflector (DBR) formed by a finite number of periods, N, based on porous silicon (PSi). The DBR is obtained by stacking periodically N times two constituent layers A and B with different porosity. The thicknesses of A and B layers are denoted by dA and dB, respectively. In PSi, the mass density is related with the porosity P via ρ=ρ0(1−P), being ρ0 the mass density of bulk silicon. The propagating velocity of the longitudinal waves through PSi is related with the porosity as v=v0(1−P)k, being v0 the longitudinal velocity of the wave in bulk silicon, and k is a parameter. The acoustic impedance of layer A(B) is given by ZA=ρAvA(ZB=ρBvB). The dependence on porosity of these two parameters, velocity and mass density, requires a very accurate control of the etching process. We grow the layers forming the DBR according to the methodology reported in [23]. The thicknesses of the PSi layers are dA=2 μm and dB=1.65 μm, and their respective porosities equal to PA=0.47 and PB=0.67. The parameters k=0.56, ρ0=2,330 kg/m 3, and v0=8,440 m/s were determined by fitting the experimental results as is explained in [23]. In Figure 2, the theoretical (black solid line) and experimental (blue solid line) transmission coefficients (in dB), =jt/ji, are plotted as a function of the acoustic phonon frequencies (in GHz). The theoretical transmission spectrum has been modeled using the transfer matrix method and includes the effect of the sample, transducers, and liquid coupling the transducers to the sample. All these components are included writing their corresponding transfer matrix like appearing in (8) but using their respective parameters. The experimental transmission has been measured using a vector network analyzer (VNA) using two piezoelectric transducers according to the experimental setup described in [23]. We can observe that in the range of frequencies reported here, the DBR shows two acoustic gaps, corresponding to the first and second acoustic gap, and centered around ≅0.71 and ≅1.42 GHz, respectively. In the inset, we show the calculated phonon transmittance for frequencies around the zone center of the first minigap. In particular, note that at frequencies fL≅0.54 GHz and fR≅0.88 GHz, corresponding to the edges of the gap, =0 dB. We can observe a noticeable difference in the fundamental bandgap between the experimental and the theoretical spectrum. The experimental bandgap has a depth of approximately 50 dB which is less than the modeled value of approximately 200 dB. However, this is due to the experimental limitation on the setup used, and it is only attributed to the noise of our VNA.

thumbnailFigure 2. Transmission coefficients. Theoretical and experimental transmission coefficients (in dB) plotted as a function of the acoustic phonon frequencie (in GHz) for the DBR, with N=18. The inset shows some details of the theoretical transmission coefficient around the first acoustic stop band.

The behavior of the vibrational energy density (1) as a function of the normalized distance, along the axis of the DBR, is shown in Figure 3a. We observe that the vibrational energy density is a one piecewise constant function. Let us now investigate the behavior of at the resonances. At the frequencies fL≅0.54 GHz (blue solid line) and fR≅0.88 GHz (red solid line), the vibrational energy density exhibits the maximum value in the center region of the DBR. However, these modes are characterized by complete transmission of the vibrational energy through the system. The one-dimensional stress tensor σ and the relative displacements u are shown in Figure 3b,c, respectively. In these figures, we have included a schematic representation of the structure, representing with dark (light) regions the layers with low (high) porosity. We observed from Figure 3c that the relative displacement for the mode appearing at fL≅0.54 GHz (fR≅0.88 GHz) is localized in the low (high) porosity regions, which is consistent with the agreement used to refer at the valence (dielectric) and conduction (air) band used in semiconductors (photonic crystals). From these figures, it is clearly observed that the phonon modes propagate through the structure.

thumbnailFigure 3. The vibrational energy density, the one-dimensional stress tensor and the relative displacements.(a) The vibrational energy density (in a.u.), (b) the one-dimensional stress tensor σ (in a.u.), (c) and the relative displacements u (in a.u.) in the structure, plotted as a function of the normalized distance along the axis of the DBR. The distance is normalized by the factor lc=dA+dB.

The phonon modes with energies within the first minigap appear marked by the black solid line in Figure 3b,c. We observe that at the frequency fC≅0.70 GHz, the amplitude of the wave shows an abrupt decay along the axis of the DBR. This qualitative behavior is very similar for the other physical quantities.

Tunneling times

The dwell time, tD, for acoustic phonons inside the DBR is calculated by the following equation [21]:

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M49','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M49">View MathML</a>


where H is the total energy in the interval (0,L), being L the total length of the DBR. This energy is obtained by integrating Equation 1, i.e.,

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M50','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M50">View MathML</a>


The free time, tf, which is the time associated to the transit of the sound pulse along a distance L with velocity, vf, is simply defined by

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/449/mathml/M51','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/449/mathml/M51">View MathML</a>


being vf the phase velocity.

In Figure 4, the tunneling times as a function of the phonon frequency (in GHz) are depicted for the acoustic phonons around the second acoustic gap. In Figure 4a, we can observe that the dwell time (36) increases near the gap. The behavior at such frequencies is the same as if the phonon was trapped for a long time in the spatial region occupied for the structure before being transmitted. In Figure 4b, we observe that, in the stop bands, the transmission time (33) is shorter than free time (38), though a large noise is observed due to the coexistence of the liquid modes. The very short transit times in the stop bands are the acoustic equivalent of the Hartman effect [4] of electrons tunneling through potential barriers [8,9].

thumbnailFigure 4. Tunneling times.(a) The dwell time tD (in μs) and (b) the transmission time τt and free time tf (in μs) plotted as a function of the acoustic phonon frequencies (in GHz). (c) Theoretical and experimental transmission coefficients (in dB) plotted as a function of the acoustic phonon frequencies (in GHz).


In this paper, we studied tunneling times of acoustic phonon packets through a distributed Bragg reflector made of porous silicon layers. Under the assumption that the long-wavelength approximation is valid, and from general principles of symmetry and conservation, we report an explicit demonstration of the equivalence of the transmission and reflection times in the tunneling of long-wavelength phonons. Calculations of the vibrational energy density and the vibrational energy stored within the structure allows a better visualization of the physical phenomena occurring in this system. The description of the stress and strain fields complements the energetic description. We report the possibility that a phenomenon called Hartman effect appears in porous silicon multilayer structures, an acoustic equivalent of Hartman effect of electrons tunneling through potential barriers. The results of this study could be useful for the design of acoustic devices.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

All authors contributed to the analytical calculations. ZL performed the experimental measurements and fabrication of the porous silicon DBR. All authors contributed to the writing of the manuscript. All authors read and approved the final manuscript.


DV and RPA acknowledge the hospitality of Benemérita Universidad Autónoma de Puebla, México. DV acknowledges the support of S.R.E. (México). JA acknowledges Conacyt, México for the partial support under grant no. 167939.


  1. Balandin AA: Nanophononics: phonon engineering in nanostructures and nanodevices.

    J Nanosci Nanotechnol 2005, 5:1015-1022. PubMed Abstract | Publisher Full Text OpenURL

  2. Narayanamurti V, Stormer HL, Chin MA, Gossard AC, Wiegmann W: Selective transmission of high-frequency phonons by a superlattice: the “dielectric” phonon filter.

    Phys Rev Lett 1979, 43:2012-2016. Publisher Full Text OpenURL

  3. Winful HG: Delay time and the Hartman effect in quantum tunneling.

    Phys Rev Lett 2003, 91:260401. PubMed Abstract | Publisher Full Text OpenURL

  4. Hartman TE: Tunneling of a wave packet.

    J of Appl Phys 1962, 33:3427-3433. Publisher Full Text OpenURL

  5. Robertson W, Ash J, McGaugh J: Breaking the sound barrier: tunneling of acoustic waves through the forbidden transmission region of a one-dimensional acoustic band gap array.

    Am J of Phys 2002, 70:689. Publisher Full Text OpenURL

  6. Page J, Suxia Y, Cowan M, Zhengyou L, Chan C, Ping S: 3D Phononic cristals. In Wave Scattering in Complex Media: from Theory to Applications. Edited by van Tiggelen B, Skipetrov S. Amsterdam: Kluwer Academic; 2003:283-307. OpenURL

  7. Yang S, Page J, Liu Z, Cowan M, Chan CT, Sheng P: Ultrasound tunneling through 3D phononic crystals.

    Phys Rev Lett 2002, 88:104301. PubMed Abstract | Publisher Full Text OpenURL

  8. Villegas D, de León-Pérez F, Pérez-Álvarez R: Tunneling time for phonons: dependence on the system’s size.

    Phys Status Solidi (b) 2005, 242:1767-1770. Publisher Full Text OpenURL

  9. Huynh A, Lanzillotti-Kimura N, Jusserand B, Perrin B, Fainstein A, Pascual-Winter M, Peronne E, Lemaître A: Subterahertz phonon dynamics in acoustic nanocavities.

    Phys Rev Lett 2006, 97:115502. PubMed Abstract | Publisher Full Text OpenURL

  10. Lanzillotti-Kimura ND: Dinámica de fonones acústicos en multicapas nanométricas.

    PhD thesis Instituto Balseiro, SC Bariloche, Argentina 2009. OpenURL

  11. Ghulinyan M, Otton CJ, Gaburro Z, Bettotti P, Pavesi L: Porous silicon free-standing coupled microcavities.

    Appl Phs Lett 2003, 82:1550-1552. Publisher Full Text OpenURL

  12. Parsons LC, Andrews GT: Observation of hypersonic phononic crystal effects in porous silicon superlattices.

    Appl Phys Lett 2009, 95:241909. Publisher Full Text OpenURL

  13. Aliev GN, Goller B, Kovalev D, Snow PA: Hypersonic acoustic mirrors and microcavities in porous silicon.

    Appl Phys Lett 2010, 96:124101. Publisher Full Text OpenURL

  14. Da Fonseca RJM, Saurel JM, Foucaran A, Massone E, Talierco E, Camassel J: Acoustic microscopy investigation of porous silicon.

    Thin Solid Films 1995, 225:155. OpenURL

  15. Andrews GT, Zuk J, Kiefte H, Clouter MJ, Nossarzewska-Orlowska E: Elastic characterization of a supported porous silicon layer by Brillouin scattering.

    Appl Phys Lett 1996, 69:1217. Publisher Full Text OpenURL

  16. Thomas L, Aliev GN, Snow PA: Hypersonic rugate filters based on porous silicon.

    Appl Phys Lett 2010, 97:173503. Publisher Full Text OpenURL

  17. Aliev GN, Goller B, Snow PA: Porous silicon bulk acoustic wave resonator with integrated transducer.

    Nanoscale Res Lett 2012, 7:378. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  18. de León-Pérez F, Pérez-Alvarez R: Phonon propagation in nonpolar semiconductor heterostructures.

    Phys Rev B 2001, 63:245304. OpenURL

  19. Kavokin A, Malpuech G, Di Carlo A, Lugli P: Photonic Bloch oscillations in laterally confined Bragg mirrors.

    Phys Rev B 2000, 61:4413. Publisher Full Text OpenURL

  20. Karabutov AA, Kosevichc Y, Sapozhnikov OA: Bloch oscillations of an acoustic field in a layered structure.

    Acoust Phys 2013, 59:158-169. OpenURL

  21. Villegas D, de León-Pérez F, Pérez-Alvarez R: Tunneling time of long-wavelength phonons through semiconductor heterostructures.

    Phys Rev B 2005, 71:035322. OpenURL

  22. Hauge EH, Falck JP, Fjeldly TA: Transmission and refletion times for scattering of wave packets off tunneling barriers.

    Phys Rev B 1987, 36:4203. OpenURL

  23. Lazcano Z, Arriaga J, Aliev GN: Experimental and theoretical demonstration of acoustic Bloch oscillations in porous silicon structures.

    J Appl Phys 2014, 115:154505. Publisher Full Text OpenURL