SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Nano Express

Localization of acoustic modes in periodic porous silicon structures

Zorayda Lazcano*, Octavio Meza and Jesús Arriaga

Author Affiliations

Instituto de Física, Benemérita Universidad Autónoma de Puebla, 18 Sur y San Claudio, Edif. 110-A, Puebla 72570, México

For all author emails, please log on.

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


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


Received:12 May 2014
Accepted:4 July 2014
Published:21 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.

Abstract

The propagation of longitudinal acoustic waves in multilayer structures based on porous silicon and the experimental measurement of acoustic transmission for the structures in the gigahertz range are reported and studied theoretically. The considered structures exhibit band gaps in the transmission spectrum and these are localized modes inside the band gap, coming from defect layers introduced in periodic systems. The frequency at which the acoustic resonances appear can be tuned by changing the porosity and/or thickness of the defect layer.

Keywords:
Porous silicon; Localized defects; Phononic crystals; Acoustic transmission

Background

The study of acoustic and elastic wave propagation in phononic crystals (PCs) [1-3] have been studied theoretically [4] and experimentally [5] in recent years. In analogy with the photonic band gap materials, emphasis in phononic crystals has been on achieving large acoustic band gaps within which propagation of sound is forbidden.

The interest about phononic band gap materials is in both, pure and applied physics. Fundamental physics has a special interest concerned with the localization phenomena of sound and vibrations in PCs. Researchers have prospected numerous applications based on cavity structures built around PCs, such as wave filters, waveguides, and splitters [6-9]. Furthermore, it is possible to design cavities for coherent (single-wavelength) phonon generation and control, to attain phonon amplification and ‘lasing’ in the called ‘saser’, one of the most important potential applications [10-12].

Periodic solid-state structures exhibit transmission stop bands for waves at certain frequencies. By placing one or more defects into a perfect phononic crystal, acoustic cavities are created inside the system. The presence of these defects, produces localization of elastic or acoustic modes inside the phononic band gap. These localized modes are the acoustic analog of donor or acceptor states produced inside the band gap of semiconductors. In analogy with electronic systems, one can consider these acoustic states to control the sound propagation through the structure.

If a defect is introduced into a periodic structure, the translational symmetry is broken and highly localized defect modes within the band gaps are created [6,8,13,14]. Point, linear, and planar defect states have been theoretically investigated in one-dimensional (1D), two-dimensional (2D), and three-dimensional (3D) phononic crystals [3,15,16]. In 1D structures, a microcavity can be a spacer layer of thickness λ/2 enclosed by two Bragg reflectors [17]. In 2002, Trigo et al. proposed phonon cavities in structures consisting of two semiconductor superlattices enclosing a spacer layer, showing that acoustical phonons can be confined in such layered structures if the spacer thickness is an integer multiple of the acoustic half-wavelength at the center of one of the superlattice-folded minigaps. These acoustic cavities are semiconductor multilayers in the nanometer scale and are fabricated by molecular beam epitaxy (MBE), which is a sophisticated and expensive technique that requires ultra-high vacuum system and a very tight control on the growth parameters, and modulate the thicknesses is easier than to modulate the elastic properties of the layers. Contrasting, porous silicon (PS) multilayer fabrication is relatively easy and considerably less costly, besides that this material allows to modulate both the thicknesses and the elastic properties of each layer.

PS is known as a versatile material with applications in light emission, sensing, and photonic devices [18]. The possibility of producing acoustic band gaps in PS was proposed in 2005 [19], and detailed calculations of predicted bandwidths were subsequently published [20]. Recently, experimental results of Brillouin light scattering suggested the existence of zone-folded phonons and phononic band gaps in PS multilayers [21]. G. N. Aliev et al. in 2010 reported experimental measurements of transmission for longitudinal acoustic waves and observation of hypersonic band gaps in a periodic PS structure with an acoustic microcavity consisting of a spacer layer of thickness λ/2 [22].

By taking advantage of the possibility to modulate the elastic properties of PS layers, and considering that it is possible to create localized modes by introducing a defect layer with different acoustic properties into a periodic structure, in this paper, we investigate the propagation of longitudinal acoustic waves in multilayer structures based on PS, that exhibit resonant cavity modes in frequencies of gigahertz (GHz), consisting of defect layers intentionally introduced in periodic structures. The design and material parameters that allow to create these localized acoustic modes is discussed, and experimental results of the measured acoustic transmission in PS samples fabricated by electrochemical etching are presented.

Methods

Theoretical models

The multilayer PS structures studied here have thicknesses in micrometer range and the procedure used to fabricate them creates mesoporous silicon with an average pore diameter of 20 to 50 nm. On the other hand, in our experiments, the typical longitudinal wavelengths excited throughout the samples are 3 to 7 μm depending on porosity. Accordingly, each of the individual layers in the structures is assumed to be homogeneous.

The longitudinal acoustic wave equation in the continuum limit for a solid inhomogeneous along the z direction (but homogeneous along the x and y directions) is given by [23],

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

(1)

where ρj is the mass density, <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/419/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/419/mathml/M2">View MathML</a> and u(z,t) is the atomic displacement. Here, j is an index identifying each layer. The limits of the elastic continuum description of wave propagation in ordered media depends on the dimensions of the system compared with the wavelength. When the dimensions approach nanometer-length scales, atomistic treatments using first principles or semi-empirical methods may become necessary [24]. However, in our case, the thicknesses of the layers are in the micrometer range and each layer can be considered as a homogeneous layer; thus, the model described before is assumed valid.

In a solid, the acoustic waves can be longitudinal or transversal. In this letter, only longitudinal waves propagating through PS are considered because in our experiments, the waves are coupled to the samples through a liquid at normal incidence. The mass density ρ is a function of the porosity and is described by ρ=ρ0(1−P) where ρ0=2.330 g/cm 3 is the density of bulk silicon and P the porosity. The acoustic velocity dependence on porosity is given empirically by vL=vL0(1−P)k, being vL0 the longitudinal velocity of sound in bulk silicon along the (100) crystallographic direction and k≥0.5 is a constant [25-28]. In general, the parameter k depends on PS morphology which in turn depends on the doping level of the Si substrate [25,26].

The general solution of Equation 1 takes the form,

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

(2)

corresponding to waves propagating to the right (A+) and left (A). At each interface, this solution must satisfy the boundary conditions related to the continuity of the atomic displacement and stress,

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

(3)

and

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

(4)

respectively. Here, dj denotes the position of the j-th interface between j and j+1 layers. The frequency ω is related to its wave vector via ω=kjvj, with vj the sound speed in the j-th layer and ω=2πf, being f is the frequency in s−1. Using the transfer matrix method (TMM) [26], we can relate the amplitudes of the fields <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/419/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/419/mathml/M6">View MathML</a> and <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/419/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/419/mathml/M7">View MathML</a> in the layer j of the system with the amplitudes of the wave in the j+1 layer according to

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

(5)

The transfer matrix Tj appearing in the previous equation propagates the amplitudes through a layer with thickness dj, mass density ρj, and sound longitudinal velocity vLj, and is given explicitly by,

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

(6)

If we consider a structure formed by N layers, the total transfer matrix representing the structure is obtained by multiplying, in the appropriate order, a series of N transfer matrices, each one given by a matrix of the type appearing in Equation 6. The obtained matrix relates the displacement vector at the beginning of the structure with that at the end, and represents a 2 × 2 set of equations that can be fully solved.

With the above formalism, one can derive the acoustic eigenenergies and eigenvectors. The reflectivity and transmission can also be calculated as the square modulus of <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/419/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/419/mathml/M10">View MathML</a> and <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/419/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/419/mathml/M11">View MathML</a>, respectively, by imposing the boundary conditions <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/419/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/419/mathml/M12">View MathML</a> and <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/419/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/419/mathml/M13">View MathML</a> for a wave traveling from right to left. Here 0 and N label the first and last layer of the structure, respectively.

Attenuation can be included by taking the wave vector kj complex, such that Kj=kjαi, where αi is attenuation coefficient. The form of the attenuation coefficient depends on the physical process causing loss and we assume that the Akhiezer model is dominant in a semiconducting material. This gives α=ηω2/2ρv3, where η is the viscosity. However, it is known that introducing acoustic attenuation into the model leads to important effects as the shrinking of gaps, only for frequencies higher than 180 GHz [29]; therefore, no absorptive behavior is considered in our model since no important effects are obtained if they are included.

Furthermore, the position and width of the band gap are critical parameters for devices that reflect or localize the acoustic waves [30]. Band structures of many kinds of periodic phononic crystals have been reported [31-33]. The most commonly studied acoustic band gaps in 1D PCs are the Bragg type, appearing at an angular frequency ω of the order of vL(T)/d (vL(T) refers to the longitudinal (transverse) elastic wave velocity and d is the lattice constant).

An acoustic Bragg mirror can be made by repeating n times a basic block of two materials with different acoustic properties. The mismatch in the acoustic impedance between these layers results in waves that are reflected and interfere, giving an acoustic band gap around a central frequency fB. For normal incidence, this frequency is given by

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

(7)

being m the order of the stop band, d1 and d2 are the layer thicknesses, and Z1 and Z2 are the acoustic impedances of layers 1 and 2, respectively. The acoustic impedance Z is given by ρv, with v as the sound velocity and ρ as the mass density. The condition ρ1d1/Z1=ρ2d2/3Z2 optimize the stop-band width and reflectivity, corresponding in an infinite stack, to the first minigap at the Brillouin zone center. The reflectivity at the center of the stop-band depends on the acoustic impedance mismatch between the two materials Z2/Z1, and for n pairs of layers is given by [17,22],

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

(8)

In [34], the authors considered periodic semiconductor structures of GaAs/AlAs to introduce microcavities as spacer layers of thickness λ/2. However, for a 10-period GaAs/AlAs mirror, RB∼0.880, while RB∼0.996 if n=20. For a PS structure, a porosity variation of 15 % between the constituent layers of 52 % and 67 % of porosity, leads to RB∼0.997 for n=6. Thus, by modulating the porosity of the PS structures, very high reflectivity values can be achieved. This is an essential condition to obtain narrow transmission bands into the stop bands corresponding to the cavity modes.

To demonstrate the localization in time domain, we consider the propagation of a Gaussian pulse through the structure. The Gaussian pulse is described by g(f)= exp(−4π[(ff0)/σ]2), were f0 is the central frequency and σ the pulse width. In response to the incident pulse, the time and spatial variations of the displacement field u(z,t) inside the sample can be calculated according to the scattering state method as [35],

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

(9)

where u(z,f) is the displacement field distribution at each frequency, which is obtained by the transfer matrix method.

Experimental details

Samples were electrochemically etched from boron-doped (100)-oriented Si substrates with a resistivity of 0.007 to 0.013 Ωcm. Room-temperature anodization was performed using a 1:1 solution of HF (40 %) and ethanol (99.98 %). The acoustic transmission measurements reported here were done using a Vector Network Analyzer (VNA). Each sample was placed between two ZnO-based piezoelectric transducers with a central frequency of 1.1 GHz and an operation bandwidth of 500 MHz. The transducers consist of a piezoelectric layer driving waves into a silicon pillar with a thickness of 520 μm. To couple the transducers to the specimen, In-Ga eutectic was used. The transducer front surface was aligned parallel to the sample surface using two orthogonal microscopes so that the acoustic waves impinge normally into the PS layers. The transducers were connected to the VNA ports and transmission parameters were measured as function of frequency, more details of the experimental set-up can be found in [36]. Dependence of sound velocity with porosity for Si wafers and experimental conditions used here were determined previously and reported in [36], the values of vL0 and the parameter k used here are 8.44 km/s and 0.56, respectively.

For the fabrication of PS multilayers, we consider the inclusion of ‘etch stops’ or ‘etch breaks’ where the current is interrupted to stop the etching of the Si wafer in order to prevent any depletion of HF [37]. The introduction of these etching breaks is necessary to obtain layers with constant porosity with depth [38]. Because our samples include very thick layers, with large mismatch porosities between them, the number and length of the etch breaks are important to obtain homogeneous structures. We found that etch breaks of 4 s with a ratio (etch break time)/(etching time) from 3.3 for low porosities (52 %) to 7.3 for high porosities (67 %) are enough to minimize any chirp in the layers.

Results and discussion

Thicknesses of layers were measured by optical microscopy, and the layer porosities were determined from optical reflectance spectra by fitting our experimental measurements and comparing them with our theoretical simulations for each sample. The acoustic transmission and field intensity distribution have been modeled using the transfer matrix method described before and taking into account the effect of the sample (PS-Si substrate), transducers (Si pillars), and In-Ga eutectic liquid used to couple the transducers to the sample.

Three PS multilayer samples are considered here to show the effect of localization inside the structures. All of them consist of layers a and b repeating alternatively, and a defect layer, c, in the middle of the structure. The sequence used for structures was ababababababcbabababababa=(ab)6c(ba)6. In the first sample (1) porosities and thicknesses of layers a and b are Pa=53%, da=1.15 μm, Pb=67%, and db=1.10 μm, respectively. Here, layer c has the same thickness and porosity of layer a, and therefore, this sample is completely periodic. The porosities and thicknesses of the layers were chosen to obtain the fundamental stop band within the bandwidth of the acoustic transducers, and satisfying Equation 7. A scheme of structure 1 is displayed in the top panel of Figure 1. The central panel of Figure 1 (solid line) shows the measured acoustic transmission spectrum of the PS periodic structure with a total thickness of approximately 27 μm. The band gap in the transmission spectrum observed around 1.15 GHz and ranged from 0.94 to 1.38 GHz is the first-order acoustic stop band of the mirror, corresponding to m=1 in Equation 7. This fundamental stop band shows an attenuation of approximately 50 dB with a fractional bandwidth of 37 %. The dashed line is the result of calculations using TMM. Good agreement between modeled and measured spectra is seen. The fine features of the spectrum are not noise but the longitudinal modes of the Si pillars of the transducers and the Si substrate of the sample. The fundamental band gap has a depth of approximately 50 dB which is less than the modeled value of approximately 100 dB. However, this is an experimental limitation on the setup used and its noise floor. Thus, we can conclude that the stop band has a depth of at least 50 dB. The bottom panel of Figure 1 shows the squared displacement field corresponding to the central frequency of the gap, 1.15 GHz. The dashed line represents the material acoustic impedance and is useful to identify the position in the sample. As can be seen, the displacement field is not localized, as is expected.

thumbnailFigure 1. Acoustic transmission and distribution of the displacement field for the periodic case, sample 1. (Top) Scheme of the periodic structure consisting of 12.5 periods of layers a and b. (Middle) Acoustic transmission spectra, measured in solid line and calculated in dashed line. The measured transmission, recorded on a logarithmic scale, is normalized to its maximum and corrected by an envelope function of the transducer response. (Bottom) In solid line, squared phonon displacement corresponding to the central frequency of the gap. The dashed line represents the material acoustic impedance and serve to identify the position in the sample.

Now, based on the concepts mentioned before about cavities, we will show how the intentional introduction of a defect layer between a pair of mirrors can lead to formation of an acoustic cavity mode within the stop band. For this purpose, we consider two structures: sample 2 and sample 3.

In sample 2, porosities and thicknesses of layers a, b, and c are: da=1.15 μm, Pa=52%, db=1.00 μm, Pb=65%, dc=1.15 μm and Pc=74%, respectively. The defect (layer c) corresponds to a layer with the same thickness, as the periodic case, but higher porosity (lower impedance), as is shown schematically at the top of Figure 2. In the middle of Figure 2 are shown the acoustic transmission spectra, measured experimentally (solid line) and calculated theoretically (dashed line). The introduction of the defect layer results in well-localized transmission modes at 1.01 and 1.27 GHz, within the fundamental stop band ranged from 1.02 to 1.47 GHz and with a fractional bandwidth of 35 %, as it can be seen in the transmission spectrum. At the bottom of the Figure 2 is shown (in solid line) the displacement field distribution as a function of the position in the sample, corresponding to the cavity modes, the first (thick line) and second (thin line) modes at 1.01 and 1.27 GHz, respectively. It can be seen that the amplitude of the acoustic displacement is maximum around the defect layer. The dashed line is the material acoustic impedance.

thumbnailFigure 2. Acoustic transmission and distribution of the displacement field for sample 2. (Top) Scheme of a structure consisting of two mirrors with six periods of layers a and b enclosing a defect layer of higher porosity between them. (Middle) Measured acoustic wave transmission spectrum through the sample (solid line). The dashed curve is the calculated spectrum (see text for details). (Bottom) In solid line, squared phonon displacement corresponding to the first (thick line) and second (thin line) modes at 1.01 and 1.27 GHz, respectively. The dashed line represent the layer acoustic impedance.

Sample 3, represented schematically at the top of Figure 3, contains a defect consisting of a layer with lower porosity (higher impedance) at the center of the structure. Here, thickness and porosities are: da=0.89 μm, Pa=65.5%, db=1.12 μm, Pb=53%, dc=0.89 μm, Pc=42%, for layers a, b, and c, respectively. The defect layer (c) keeps the periodicity in thickness but the porosity changes. As it can be clearly seen in measured transmission spectrum shown in Figure 3, this results in an acoustic cavity mode at 1.15 GHz within the fundamental stop band ranging from 1.02 to 1.44 GHz (34 % fractional bandwidth). The corresponding displacement field distribution for this cavity mode is shown at the bottom of the same figure (thick line) and demonstrates that the displacement field is maximum around this cavity in the same way as the second mode in sample 2. For demonstration purposes, we have calculated the displacement field for 1.46 GHz and the results are shown in Figure 3 using a thin line. Localization effects cannot be observed.In Figures 1, 2, and 3, good agreement between modeled and measured spectra is observed, and the slight differences between theoretical and experimental acoustic transmissions are due to features of porous silicon layers which are not considered here, as the roughness at the interfaces, as well as intrinsic error coming from the measured procedure, and not to absorption properties, as was explained before.

thumbnailFigure 3. Acoustic transmission and distribution of the displacement field for sample 3. (Top) Scheme of a structure of two mirrors with six periods of layers a and b enclosing a defect layer of lower porosity. (Middle) Measured (solid line) and calculated acoustic transmission spectra (see text for details). (Bottom) In solid line, squared phonon displacement corresponding to the cavity mode frequency (thick line) at 1.15 GHz, and for a frequency of 1.46 GHz (thin line). The dashed line represent the layer acoustic impedance.

In Figure 4, we show the time-resolved displacement field u(z,t), corresponding to the time evolution of a Gaussian pulse in the samples calculated using Equation 9. Figure 4a,b corresponds to the time and spatial variations of the displacement field inside sample 2, using f0=1.01 GHz in Figure 4a and 1.27 GHz in Figure 4b. These values correspond to the frequencies where the first and the second cavity modes appear, respectively. Figure 4c shows the displacement field inside of sample 3 for f0=1.15 GHz, the frequency of the corresponding cavity mode. Figure 4d corresponds to sample 3 using f0=1.46 GHz. We use a pulse with σ=200 MHz for all cases. In Figure 4a, it can be seen that the displacement field is in the center of the PS structure, corresponding to the defect layer. Figure 4b corresponds to the same structure but the incident Gaussian pulse has a frequency equal to 1.27 GHz; in this case, localization is around the defect layer and not inside it. In Figure 4c, as Figure 4b, the localization is around the defect layer and not inside it, because this corresponds to the second mode of sample 3. Calculations for sample 3 for the peak at 1.46 GHz appears in Figure 4d, as can be seen, localization of the displacement field is not observed.

thumbnailFigure 4. Time evolution of acoustic Gaussian pulses. Calculations of ln(u(z,t)) for acoustic Gaussian pulses centered at the frequency f0 indicated in each figure and σ=200 MHz for all cases. (a,b) Sample 2. (c,d) Sample 3. Zero in x-axis is placed at the surface of the PS sample.

In order to estimate the displacement field intensity u(z,t)2, within the defect layer as a function to the time, we integrate the displacement field on the defect layer using,

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

(10)

where Φ(t) is the displacement field intensity contained in the defect as a function of the time. Figure 5a shows Φ(t) for sample 2 for two Gaussian pulses centered in the frequencies f0 indicated there, as expected the first mode has higher displacement field intensity in the defect layer because the acoustic wave is localized in the center of the PS structure, on the contrary to the second mode, see Figure 5a. In the case of sample 3, the localization is less than sample 2 for the two Gaussian pulses considered, that is, the Φ(t) amplitude for sample 2, see Figure 5b, in the first mode is around 30 times more the Φ(t) amplitude in sample 3 for one incident Gaussian pulse with frequency equal to 1.15 GHz. Finally, localization is not observed in sample 3 for a Gaussian pulse with a frequency of 1.46 GHz, as is expected, see Figure 5b.

thumbnailFigure 5. Displacement field intensity as a function of time. Theoretical calculations for the displacement field intensity u(z,t)2 in the defect as a function of time for (a) sample 2 and (b) sample 3, for frequencies indicated in each figure.

The modeled transmittance of the periodic case (sample 1) and for the two cavity structures (samples 2 and 3), obtained by the TMM, shows a good match with the experimental results. The localized acoustic resonances can be tuned at different frequencies (within the acoustic band gap) by changing the porosity of the defect layer.

Moreover, for commercial acoustic mirrors which are components of solidly mounted resonators and filters [39], a low-acoustic-impedance material such as SiO 2 is layered with high-impedance materials such as tungsten or molybdenum. Following Equation 8, for the layer pair of molybdenum and silica, where acoustic impedances are 66.2 MRayl and 13.1 MRayl, respectively, the fixed impedance ratio is 5.1, and the same impedance ratio can be obtained using PS layers of 30 % and 75 %, so, by modulating the porosity, very high reflectivity values can be achieved. Besides, by increasing the number of pairs of layers enclosing the defect layer leads to narrowing of the transmission band together with deepening of the stop band.

Conclusions

We have demonstrated theoretically by using the TMM and experimentally by acoustic transmission measured directly, the formation of acoustic cavity modes in GHz frequencies by introduction of defects into periodic structures based on PS. Acoustic resonances can be tuned at different frequencies by changing the porosity of the defect. And we proved that these resonant modes appear due to the localization of the field into the defect. The acoustic mirrors and cavity structures based on PS have a performance which is at least comparable with that devices based on semiconductor superlattices. This study could be useful for the design of acoustic devices, such as highly selective frequency filters with applications in GHz range.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

JA had the original idea of the study. ZL performed the experiments and measurements. OM and ZL did the numerical calculations. All authors contributed in the writing of the manuscript. All authors read and approved the final manuscript.

Acknowledgements

The authors acknowledge CONACyT for support under project No. 167939.

References

  1. Kushwaha MS, Halevi P, Dobrzynski L, Djafari-Rouhani B: Acoustic band structure of periodic elastic composites.

    Phys Rev Lett 1993, 71:2022. PubMed Abstract | Publisher Full Text OpenURL

  2. Kushwaha MS, Halevi P, Martínez G: Theory of acoustic band structure of periodic elastic composites.

    Phys Rev B 1994, 49:2313. Publisher Full Text OpenURL

  3. Sigalas MM, Economou EN: Attenuation of multiple-scattered sound.

    Europhys Lett 1996, 36:241. Publisher Full Text OpenURL

  4. Kushwaha MS, Halevi P: Giant acoustic stop bands in two-dimensional periodic arrays of liquid cylinders.

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

  5. Sánchez-Pérez JV, Caballero D, Martínez-Sala R, Rubio C, Sánchez-Dehesa J, Meseguer F, Llinares J, Gálvez F: Sound attenuation by a two-dimensional array of rigid cylinders.

    Phys Rev Lett 1998, 80:5325. Publisher Full Text OpenURL

  6. Sigalas MM: Elastic wave band gaps and defect states in two-dimensional composites.

    J Acoust Soc Am 1997, 101:1256. Publisher Full Text OpenURL

  7. Sigalas MM: Defect states of acoustic waves in a two-dimensional lattice of solid cylinders.

    J Appl Phys 1998, 84:3026. Publisher Full Text OpenURL

  8. Kafesaki M, Sigalas MM, Garca N: Frequency modulation in the transmittivity of wave guides in elastic-wave band-gap materials.

    Phys Rev Lett 2000, 85:4044. PubMed Abstract | Publisher Full Text OpenURL

  9. Khelif A, Djafari-Rouhani B, Vasseur JO, Deymier PA: Transmission and dispersion relations of perfect and defect-containing waveguide structures in phononic band gap materials.

    Phys Rev B 2003, 68:024302. OpenURL

  10. Lacharmoise P, Fainstein A, Jusserand B, Thierry-Mieg V: Optical cavity enhancement of light–sound interaction in acoustic phonon cavities.

    Appl Phys Lett 2004, 84:3274. Publisher Full Text OpenURL

  11. Fokker PA, Dijkhuis JI, de Wijn HW: Stimulated emission of phonons in an acoustical cavity.

    Phys Rev B 1997, 55:2925. OpenURL

  12. Camps I, Makler SS, Pastawski HM, Foa Torres LEF: GaAs-AlxGa1−xAs double-barrier heterostructure phonon laser: a full quantum treatment.

    Phys Rev B 2011, 64:125311. OpenURL

  13. Wu FG, Liu JY, Liu ZY, Liu YY: Splitting and tuning characteristics of the point defect modes in two-dimensional phononic crystals.

    Phys Rev E 2004, 69:066609. OpenURL

  14. Khelif A, Choujaa A, Benchabane S, Djafari-Rouhani B, Laude V: Guiding and bending of acoustic waves in highly confined phononic crystal waveguides.

    Appl Phys Lett 2004, 84:4400. Publisher Full Text OpenURL

  15. Psarobas IE, Stefanou N, Modinos A: Phononic crystals with planar defects.

    Phys Rev B 2000, 62:5536. Publisher Full Text OpenURL

  16. Wang ZG, Lee SH, Kim CK, Park CM, Nahm K, Nikitov SA: Acoustic wave propagation in one-dimensional phononic crystals containing Helmholtz resonators.

    J Appl Phys 2008, 103:064907. Publisher Full Text OpenURL

  17. Trigo M, Bruchhausen A, Fainstein A, Jusserand B, Thierry-Mieg V: Confinement of acoustical vibrations in a semiconductor planar phonon cavity.

    Phys Rev Lett 2002, 89:227402. PubMed Abstract | Publisher Full Text OpenURL

  18. Bisi O, Ossicini S, Pavesi L: Porous silicon: a quantum sponge structure for silicon based optoelectronics.

    Surf Sci Rep 2000, 38(1–3):5. OpenURL

  19. Kiuchi A, Gelloz B, Kojima A, Koshida N: Possible operation of periodically layered nanocrystalline porous silicon as an acoustic band crystal device.

    In Group-IV Semiconductor Nanostructures: 29 Nov - 2 Dec, Boston Edited by Tsybeskov L, Lockwood DJ, Delerue C, Ichikawa M.

    2004 Materials Research Society Symposia Proceedings Series vol. 832 2005:F371–F376

    OpenURL

  20. Reinhardt A, Snow PA: Theoretical study of acoustic band-gap structures made of porous silicon.

    Phys Status Solidi A 2007, 204:15281535. OpenURL

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

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

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

  23. Landau LD, Lifshitz M: Theory of Elasticity. Bristol, Pergamon; 1975. OpenURL

  24. Ramprasad R, Shi N: Scalability of phononic crystal heterostructures.

    Appl Phys Lett 2005, 87:111101. Publisher Full Text OpenURL

  25. Phani KK, Niyogi SK, Maitra AK, Roychaudhury M: Strength and elastic modulus of a porous brittle solid: an acousto-ultrasonic study.

    J Mater Sci 1986, 21:4335. Publisher Full Text OpenURL

  26. Maitra AK, Phani KK: Ultrasonic evaluation of elastic parameters of sintered powder compacts.

    J Mater Sci 1994, 29:4415. Publisher Full Text OpenURL

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

    Thin Solid Films 1995, 225:155-158. OpenURL

  28. Fonseca RJM, Saurel JM, Foucaran A, Camassel J, Massone E, Taliercio T: Acoustic investigation of porous silicon layers.

    J Mat Sci 1995, 30:3539. Publisher Full Text OpenURL

  29. He J, Sapriel J, Azoulay R: Acoustic attenuation and optical-absorption effects on light scattering by acoustic phonons in superlattices.

    Phys Rev B 1989, 40:1121. Publisher Full Text OpenURL

  30. Weaver RL: Anomalous diffusivity and localization of classical waves in disordered media: the effect of dissipation.

    Phys Rev B 1993, 47:1077. Publisher Full Text OpenURL

  31. Liu Z, Zhang X, Mao Y, Zhu YY, Yang Z, Chan CT, Sheng P: Locally resonant sonic materials.

    Science 2000, 289:1734-1736. PubMed Abstract | Publisher Full Text OpenURL

  32. Hirsekorn M: Small-size sonic crystals with strong attenuation bands in the audible frequency range.

    Appl Phys Lett 2004, 84:3364. Publisher Full Text OpenURL

  33. Sainidou R, Stefanou N, Modinos A: Widening of phononic transmission gaps via Anderson localization.

    Phys Rev Lett 2005, 94:205503. PubMed Abstract | Publisher Full Text OpenURL

  34. Lanzillotti Kimura ND, Fainstein A, Balseiro CA, Jusserand B: Phonon engineering with acoustic nanocavities: Theoretical considerations on phonon molecules, band structures, and acoustic Bloch oscillations.

    Phys Rev B 2007, 75:024301. OpenURL

  35. Malpuech G, Kavokin A, Panzarini G, Di Carlo A: Theory of photon Bloch oscillations in photonic crystals.

    Phys Rev B 2001, 63:035108. OpenURL

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

  37. Thönissen M, Berger MG, Billat S, Arens-Fischer R, Krüger M, Lüth H, Theiss W, Hillbrich S, Grosse P, Lerondel G, Frotscher U: Analysis of the depth homogeneity of p-PS by reflectance measurements.

    Thin Solid Films 1997, 297:92. Publisher Full Text OpenURL

  38. Lazcano Z, Arriaga J: High quality porous silicon multilayer structures for infra-red applications.

    Prog Electromagn Res 2013, 1:1404. OpenURL

  39. Satoh Y, Nishihara T, Yokoyama T, Ueda M, Miyashita T: Development of piezoelectric thin film resonator and its impact on future wireless communication systems.

    Jpn J Appl Phys 2005, 44:2883.

    Part 1

    Publisher Full Text OpenURL