SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

This article is part of the series Porous Semiconductors - Science and Technology 2012.

Open Access Nano Express

Demonstration of photon Bloch oscillations and Wannier-Stark ladders in dual-periodical multilayer structures based on porous silicon

J Octavio Estevez1, Jesús Arriaga1, Antonio Mendez-Blas1, Edgar Reyes-Ayona1, José Escorcia2 and Vivechana Agarwal2*

Author affiliations

1 Instituto de Física, Universidad Autónoma de Puebla, Puebla, A.P. J-48, México

2 Centro de Investigación en Ingeniería y Ciencias Aplicadas, UAEM, Av. Universidad 1001, CuernavacaMéxico, Morelos 62210, Col. Chamilpa

For all author emails, please log on.

Citation and License

Nanoscale Research Letters 2012, 7:413  doi:10.1186/1556-276X-7-413


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


Received:30 April 2012
Accepted:26 June 2012
Published:23 July 2012

© 2012 Estevez 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/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Theoretical demonstration and experimental evidence of photon Bloch oscillations and Wannier-Stark ladders (WSLs) in dual-periodical (DP) multilayers, based on porous silicon, are presented. An introduction of the linear gradient in refractive indices in DP structure, which is composed by stacking two different periodic substructures N times, resulted in the appearance of WSLs. Theoretical time-resolved reflection spectrum shows the photon Bloch oscillations with a period of 130 fs. Depending on the values of the structural parameters, one can observe the WSLs in the near infrared or visible region which may allow the generation of terahertz radiation with a potential applications in several fields like imaging.

Keywords:
Photon Bloch oscillations; Wannier-Stark ladders; Dual-periodical multilayers; Porous silicon; Photonic crystals

Background

The analogs between electron transport and propagation of the optical waves in dielectric structures opened the possibility of the implementation of Bloch oscillations for electromagnetic waves in photonic crystals (PC) [1,2]. The photonic analog of the above-mentioned effect appears when a PC is subjected to a slowly varying refractive index or a geometric parameter modulation, resulting in a linear tilting of the band structure. Such ‘chirped’ PCs give rise to a set of equidistant frequency levels [3], i.e., the optical counterpart of the Wannier-Stark ladders (WSLs) in semiconductor superlattices. Recently, different methods have been adopted to tilt the photonic band for the observation of WSLs and photon Bloch oscillations (PBOs) [4-7]. In confined Bragg mirrors, the band structure modifications are due to the gradual change of the lateral confinement [4]. On the other hand, in geometrically chirped PCs, the band structure modulation arises from a gradual increase in the thicknesses of the layers [5,6]. Furthermore, in graded-index optical superlattices, the index gradient comes from a linear modification in the refractive indices of the layers [7,8], and so on. The use of above-mentioned PCs with gradient in optical thickness has evoked special interests in many applications due to their novel properties [9]. Such multilayer structures originate a new type of Fabry-Perot cavity where the reflectors are replaced by nonpropagating regions associated with the local periodicity of the structure. If the linear gradient is considered, the distance between the band edges where the PBOs occur can be maintained constant. Hence, the period of PBOs remains constant with the change in the frequency of the incident wave. In case of nonlinear gradient, the distance between the band edges and, therefore, the period of the resulting PBOs can be tuned by changing the frequency of the incident light. For example in 2005 Lousse and Fan [8] reported the tunable terahertz Bloch oscillations in the chirped photonic crystals, with the potential applications in several fields, like biomedical sensing. Such useful photonic structures can be fabricated with different materials. Recently, one-dimensional photonic superlattices made of porous silicon (PSi) have allowed the demonstration of optical analogs [10] of electronic phenomena [11,12], such as PBOs, Zener tunneling, and Anderson localization [6,7,13,14].

Porous silicon provides good flexibility in the design of optical devices due to its easy fabrication technique [15-17] and tunable optical properties. PSi can be obtained by electrochemical etching of doped silicon wafers, which allows the fabrication of several types of one-dimensional (1D) porous silicon photonic bandgap structures, such as distributed Bragg reflectors [18], omnidirectional mirrors [19-22], Fabry-Perot optical microcavities [23,24], waveguides [25], rugate filters [26], and optical biosensors [27-31].

In the present work, we demonstrate the theoretical and experimental evidence of WSLs, using dual-periodical multilayer structures (shown schematically in Figure 1), with a linear gradient in refractive indices, based on porous silicon. Theoretical evidence of the presence of PBOs in such structures is also presented.

thumbnailFigure 1. Dual-periodical multilayer structure. Schematic view of the 1D dual-periodic structure showing the layer parameters, where n and m represent the period numbers of a and b in A and B substructures, respectively; n1 and n2 are high and low refractive indexes of alternating dielectric layers in a and b; the layer thickness is d1 and d2 for a, and d3 and d4 for b, respectively.

Dual-periodic structures

The optical properties of dual-periodical (DP) structures have been theoretically and experimentally reported by several groups [32,33]. Recently, Pérez et al. [34] reported DP structures from PSi multilayers. Dual-periodic structure (Figure 1) is composed of two substructures, A and B, repeating alternatively in the sequence <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/413/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/413/mathml/M1">View MathML</a>. The An and Bmare in turn composed of two different periodic units, a and b, respectively, where subscripts, n and m, are the number of periods for a and b in the A and B substructures, respectively. Both a and b consist of a pair of layers with high and low refractive indices. The thickness of the double layer a is da = d1 + d2d1 and d2 being the thicknesses for the layers with the high (n1) and low (n2) refractive indices, respectively. Similarly, the double layer b has thickness db = d3 + d4; d3 and d4 being the thicknesses for layers with the high and low refractive index as well. In particular, the following sequence was used: A2B4A2B4A2B4A2B4A2B4A2B4= <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/413/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/413/mathml/M2">View MathML</a> for the infrared region. If the substructure B is considered as a defective layer, the frequency intervals where the resonances of the transmission peaks appear can be reduced by increasing the number of periods a in the substructure A. On the other hand, if substructure A is a defective layer, the frequency intervals of the resonances can be increased by reducing the number of periods b in the substructure B. When identical AnBmstructures are coupled, a degenerate mode repulsion arises. Each degenerate optical resonance splits up and a miniphotonic band forms [32-34]. Due to the periodicity of the structure, the miniphotonic bands are separated by photonic band gaps in which propagation is prohibited. Moreover, when N>1, there are N−1 defect layers; therefore, N−1 resonance modes and N−1 transmission peaks will appear in the spectra. By adjusting the structural parameters, it is possible to tune the number, frequency, and full width at half-maximum (FWHM) of the resonance modes, opening the possibility to fabricate optical filters based on porous silicon multilayers. Such DP photonic structures are very promising in the field of optoelectronics, optical communications, and optical biosensors [34].

Furthermore, 1D translational symmetry of the system should be broken by introducing a small gradient in the refractive indices along the depth of the DP structure to obtain PBOs in periodic 1D photonic crystals. The gradient in the refractive indices results in a spatial tilting of the miniband and photonic band gaps in which the resonances, due to defects in DP structure, change slightly while preserving the mode coupling. In this way, the extended photonic states are turned into a discrete sequence of energy levels with level spacing △E, which is an optical equivalent of a WSL in frequency domain. The refractive index gradient in layers is given by △n = (<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/413/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/413/mathml/M3">View MathML</a> - <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/413/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/413/mathml/M4">View MathML</a>)/<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/413/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/413/mathml/M5">View MathML</a>, where the subscripts z1and zm are the first and the m-th layer along the depth within the sample. This gradient is the optical counterpart of the external electric field used in electronic superlattices.

Methods

The structures were fabricated by electrochemical etching of boron-doped silicon wafers with resistivity from 0.007 to 0.013 Ω· cm and (100) orientation. The substrates were etched in an electrolyte consisting of HF (40%) and ethanol (99.98%) in the volumetric ratio of 1:1. To estimate the refractive index corresponding to a given current density, the effective medium Bruggeman’s model is used. To measure the porosity of the layers, gravimetrical method was employed. The reflectivity spectra were measured using p-polarized transverse-magnetic (TM) modes light at an incidence angle of 20° using VARIAN-CARY 5000 spectrometer (Varian Inc., NC, USA ). Two different sets of structural parameters were used to fabricate the PSi structures to observe WSLs and PBOs in the near infrared region. The experimental reflectivity results were compared with the theoretical simulations. To obtain <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/413/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/413/mathml/M6">View MathML</a> structure with resonance transmission modes in the near infrared region d1 = 105, d2 = 315, d3 = 216, and d4 = 103 nm were considered. Their refractive indices were taken as n1= 2.2 and n2 = 1.4, which correspond to the porosities of 48% and 76%, and are obtained with current densities (J) of 29 and 134 mA/cm2, respectively. For obtaining a linear refractive index gradient of 16% in <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/413/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/413/mathml/M7">View MathML</a> structure (named as G(A2B4)6), n1 changes from 2.2 to 2.55 (37% porosity, J = 5.2 mA/cm2) and n2 from 1.4 to 1.57 (69% porosity, J = 89.0 mA/cm2).

Theoretical model

Light propagation in dual-periodical multilayers

Photon propagation in one-dimensional structures has been modeled by simple transfer matrix method considering p-polarized light [35]. If we consider an electromagnetic (EM) wave propagating in the structure with propagation constant k=k + kz, there are two independent EM modes: TM and transverse electric (TE). The electric and magnetic field for the TE and TM mode are perpendicular to the plane defined by the wave vector and the direction of periodicity. Using the transfer matrix method, we can relate the amplitudes of the fields <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/413/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/413/mathml/M8">View MathML</a> and <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/413/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/413/mathml/M9">View MathML</a> in the j-th cell of the system to the amplitudes of the field in the (j + 1)-th cell according to [36],

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

(1)

where <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/413/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/413/mathml/M11">View MathML</a> and <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/413/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/413/mathml/M12">View MathML</a> is the amplitude of the wave in layer j with polarization μ= sptraveling to the right and to the left, respectively. For the case considered in this work, the total transfer matrix of the system can be written as a product of matrices of the type [36]:

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

(2)

where q = <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/413/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/413/mathml/M14">View MathML</a> for p-polarization and q = kjz for s-polarization; ϕj=kjzdjkj=(ω/c)njkjz being the component of the wave vector along the growth direction of the system in the j-th layer given by <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/413/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/413/mathml/M15">View MathML</a>; and <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/413/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/413/mathml/M16">View MathML</a>, its complex refractive index. A real approach over the transfer matrix was implemented by considering the absorption and dispersion dependence. Reflectivity calculations of the system are given in terms of the matrix elements of the total transfer matrix <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/413/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/413/mathml/M17">View MathML</a> according to <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/413/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/413/mathml/M18">View MathML</a>. Using transfer matrix calculations, it is possible to compute the electric field distribution inside the structure. The electric field in the structure along the z direction for a certain wavelength can be expressed as follows:

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

(3)

where <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/413/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/413/mathml/M20">View MathML</a> are the elements of the transmission matrix from the first to the k-th interface. As can be seen from this expression, the reflectivity coefficient of the whole structure must be calculated previously before evaluating the field amplitude inside the structure. Hence, the intensity is simply the square of the electric field, |E|2 which is normalized to the amplitude of the incident signal. The amplitude of the time-resolved reflection can be calculated using the following expression:

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

(4)

where r(ω) is the reflectivity coefficient, and g(ω) is the incident pulse with a Gaussian spectral function in frequency domain:

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

(5)

where E0 is the central photon energy of the pulse and δ is the pulse width. The parameter δcontrols the pulse duration, which should be sized to a value close to the Bloch oscillations periods. The expression in Equation 4 is valid only for times longer than the duration of the incident pulse [6].

Porous silicon as an effective medium

As PSi layer consists of two distinct components (air and silicon); its complex dielectric function (εPSi) has an intermediate value between silicon (εSi) and air (εair), i.e., a volume fraction fp of silicon and a volume fraction (1 - fp) of pores (where refractive index of air = 1). Bruggeman effective medium approximation was used to estimate the effective optical parameters of PSi [37,38]:

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

(6)

From this equation we obtain εPSiin terms of εSi and the porosity fP:

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

(7)

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

(8)

The refractive index n(λ) and the extinction coefficient k(λ) can be obtained by the Cauchy model, which is useful for dielectric materials (with exponential absorption), far from the absorption bands [39].

Results and discussion

Figure 2a shows the square of the electric field inside the DP <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/413/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/413/mathml/M26">View MathML</a> structure, without any refractive index gradient (△n = 0), calculated by transfer matrix methods using Equation 3 for infrared region. The scattering states map provides the electric field intensity at each position inside the structure when plane electromagnetic waves of unit amplitude impinge onto the structure to 20° of incidence. A change from dark to bright corresponds to the increase in the light intensity. Although, the input intensity has been normalized to unity, in some of the regions, the intensity inside the structures can be more than one due to internal resonances. No field localization is observed among the different zones in the DP structure. Composed by some delocalized states in space, each bright band is an optical analog of electronic energy band caused by the periodic potential in electronic superlattices. In Figure 2b the optical reflectivity measurements performed on DP structure, fabricated for infrared region of the electromagnetic spectrum, are presented. The dotted lines in Figure 2b correspond to the results obtained by numerical simulations. The measured reflectance spectrum shows the presence of five resonance peaks at photons energy 0.837, 0.853, 0.870, 0.885, and 0.896 eV, with FWHM of 6, 8.2, 8, 7.1, and 4 meV, respectively. The FWHMs of the resonance peaks can be tuned by increasing or decreasing the number of layers in either of the substructures.

thumbnailFigure 2. Scattering states map and reflectance spectra of (A2B4)6structure. (a) Calculated electric field intensity distribution of TM modes in (A2B4)6 structure with no refractive index gradient i.e., △n= 0 (flat miniband situation). A change of color from dark to bright corresponds to the increase in the light intensity. (b) Measured and simulated reflectance spectra of the structure. The schematic illustration of the corresponding DP structure is shown just above the scattering states map.

The decreasing or the increasing number of layers of the substructure, acting as a defect in the system (substructure B for the samples considered here), results in the widening or narrowing of the peaks and consequently, the value of FWHM [34].

Figure 3a shows the square of the electric field inside the G<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/413/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/413/mathml/M27">View MathML</a> structure with △n = 16% in the same range of frequencies as in Figure 2. One can clearly observe the light confinement due to the inclination of the minibands and photonic bandgaps (PBGs) (Figure 3a). The occurrence of WSL as a series of narrow resonant peaks is clearly observed in the reflectivity spectra of the structure (Figure 3b). When the incident pulse enters in G<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/413/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/413/mathml/M28">View MathML</a> structure in the frequency region between two local gaps, most of the light is reflected back on the lower band gap, and only a small percentage which crosses the gap elastically (Zener tunneling) feeds the horizontal resonances associated with the WSLs. Hence, the confined light starts to oscillate within the inclined allowed miniband with a regular period of τB = h/△E. These are PBOs which represent the time domain counterpart of the photonic WSLs. Time domain oscillations can be detected by measuring the time-resolved reflection of the sample. The center of the Wannier-Stark resonance defines, in space, the center of the PBO. The distance between successive frequency levels can be obtained by the following Fabry-Perot formula: △E = hc/<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/413/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/413/mathml/M29">View MathML</a>, where d is the local distance between the gap edges at the given input frequency and <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/413/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/413/mathml/M30">View MathML</a> is the average refractive index. The measured values of △E, from Figure 3b, are found to be in agreement with the theoretical simulations presented in Figure 3a (as a function of film depth). The structural parameters are modulated to obtain the energy spacing △E ≈ 27 meV, and therefore, τB = 148 fs.

thumbnailFigure 3. Scattering states map and reflectance spectra of G(A2B4)6structure. (a) Calculated electric field intensity distribution of TM modes in G(A2B4)6 structure with a linear gradient in refractive index (△n = 16%). Inclined white lines are shown as indicators of the miniphotonic band inclination. One can observe the formation of the photonic WSLs confined between two minigaps. (b) Measured and simulated reflectance spectra of the structure and the input Gaussian pulse (red dotted line) in frequency domain.

Figure 4 shows the theoretical time-resolved reflection of the G<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/413/mathml/M31','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/413/mathml/M31">View MathML</a> structure calculated using Equation 4 for a Gaussian pulse with E0= 809 meV (δ = 10.5 meV), as shown in Figure 3b, and E0= 830 meV (δ= 11 meV). One can observe the PBOs with an oscillation period of 111 and 130 fs(Figure 4), compared to 148 fs (△E = 27 meV) measured from the reflectance spectrum of the corresponding structure (Figure 3b). These oscillations are observable due to the tunneling of photons through the lower inclined gap. The form and the energy of the pulse are found to slightly influence the period of the PBOs, e.g., an increase in the pulse energy by 21 meV results in a increase of PBO period by 19 fs.

thumbnailFigure 4. Time-resolved reflection spectra in G(A2B4)6structure. Calculated time-resolved reflection spectra of the G(A2B4)6 structure for incident pulse energy of (top image) E0= 809 meV (bottom image) E0= 830 meV.

Conclusions

Successful demonstration of WSLs and consequently the PBOs in one-dimensional dual-periodical PSi structures for near infrared range is presented. The possibility of observing the WSLs in reflectance spectrum is evaluated by selecting the adequate parameters with a linear gradient in refractive indices. Such photonic structures can be very promising in the observation of Bloch oscillations in the different regions of the electromagnetic spectrum.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

JOE prepared the PSi multilayers. JA, ER-A, and JE performed the treatment of experimental data and simulations. AMB and VA proposed the experimental conditions and measured the structures. JOE, JA,and VA prepared the manuscript initially and participated in its design and coordination. All authors read and approved the final manuscript.

Acknowledgements

This work has been partially supported by CONACyT under project 128953 (for VA) and 101486 (for AMB).

References

  1. Yablonovitch E: Inhibited spontaneous emission in solid-stated physics and electronics.

    Phys Rev Lett 1987, 58:2059-2062. PubMed Abstract | Publisher Full Text OpenURL

  2. Joannopoulos JD, Meade RD, Winn JN: Photonic Crystals: Molding the Flow of Light. New Jersey: Princeton University Press; 1995. OpenURL

  3. Monsivais G, del Castillo-Mussot M, Claro F: Stark-Ladder resonances in the propagation of electromagnetic waves.

    Phys Rev Lett 1990, 64:1433-1436. PubMed Abstract | Publisher Full Text OpenURL

  4. Kavokin A, Malpuech G: Photonic Bloch oscillations in laterally confined Bragg mirrors.

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

  5. Malpuech G, Kavokin A, Panzarini G, DiCarlo A: Theory of photon Bloch oscillations in photonic crystals.

    Phys Rev B 2001, 63:035108. OpenURL

  6. Agarwal V, del Río JA, Malpuech G, Zamfirescu M, Kavokin A, Coquillat D, Scalbert D, Vladimirova M, Gil B: Photon Bloch oscillations in porous silicon optical superlattices.

    Phys Rev Lett 2004, 92:097401. PubMed Abstract | Publisher Full Text OpenURL

  7. Sapienza R, Costantino P, Wiersma D, Ghulinyan M, Oton CJ, Pavesi L: Optical analogue of electronic Bloch oscillations.

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

  8. Lousse V, Fan S: Tunable terahertz Bloch oscillations in chirped photonic crystals.

    Phys Rev B 2005, 72:075119. OpenURL

  9. Barthelemy P, Ghulinyan M, Gaburro Z, Toninelli C, Pavesi L, Wiersma DS: Optical switching by capillary condensation.

    Nature Photonics 2007, 1:172-175. Publisher Full Text OpenURL

  10. Ghulinyan M, Gelloz B, Ohta T, Pavesi L, Lockwood DJ, Koshida N: Stabilized porous silicon optical superlattices with controlled surface passivation.

    Appl Phys Lett 2008, 93:061113. Publisher Full Text OpenURL

  11. Bloch F: Über die Quantenmechanik der Elektronen in Kristallgittern.

    Z Phys 1928, 52:555. OpenURL

  12. Zener C: A theory of electrical breakdown of solid.

    Dielectrics Proc R Soc 1934, 145(London A):523. OpenURL

  13. Ghulinyan M, Oton CJ, Gaburro Z, Pavesi L: Zener tunneling of light waves in an optical superlattice.

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

  14. Bertolotti J, Gottardo S, Wiersma DS: Optical necklace states in Anderson localized 1D systems.

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

  15. DeLouise LA, Miller BL: Quantitative assessment of enzyme immobilization capacity in porous silicon.

    Anal Chem 2004, 76:6915-6920. PubMed Abstract | Publisher Full Text OpenURL

  16. Pavesi L, Turan R: Silicon Nanocrystals. Weinheim: WILEY-VCH; 2010. OpenURL

  17. Muñoz-Sandoval E, Agarwal V, Escorcia-García J, Ramírez-González D, Martínez-Mondragón MM, Cruz-Silva E, Meneses-Rodríguez D, Rodríguez-Manzo JA, Terrones H, Terrones M: Architectures from aligned nanotubes using controlled micropatterning of silicon substrates and electrochemical methods.

    Small 2007, 3:1157. PubMed Abstract | Publisher Full Text OpenURL

  18. Mazzoleni C, Pavesi L: Application to optical components of dielectric porous silicon multilayers.

    Appl Phys Lett 1995, 67:2983-2985. Publisher Full Text OpenURL

  19. Bruyant A, Lerondel G, Reece PJ, Gal M: All-silicon omnidirectional mirrors based on one-dimensional photonic crystals.

    Appl Phys Lett 2003, 82:3227-3229. Publisher Full Text OpenURL

  20. Estevez JO, Arriaga J, Méndez Blas A, Agarwal V: Omnidirectional photonic bandgaps in porous silicon based mirrors with a Gaussian profile refractive index.

    Appl Phys Lett 2008, 93:191915. Publisher Full Text OpenURL

  21. Estevez JO, Arriaga J, Méndez Blas A, Agarwal V: Enlargement of omnidirectional photonic bandgap in porous silicon dielectric mirrors with a Gaussian profile refractive index.

    Appl Phys Lett 2009, 94:061914. Publisher Full Text OpenURL

  22. Ariza-Flores D, Gaggero-Sager LM, Agarwal V: Omnidirectional photonic bandgap in dielectric mirrors: a comparative study.

    J Phys D Appl Phys 2012, 45:015102. Publisher Full Text OpenURL

  23. Pavesi L, Mazzoleni C: Controlled photon emission in porous silicon microcavities.

    Appl Phys Lett 1995, 67:3280-3282. Publisher Full Text OpenURL

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

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

  25. Apiratikul P, Rossi AM, Murphy TE: Nonlinearities in porous silicon optical waveguides at 1550 nm.

    Opt Express 2009, 17:3396-3406. PubMed Abstract | Publisher Full Text OpenURL

  26. Salem MS, Sailor MJ, Sakka T, Ogata YH: Electrochemical preparation of a rugate filter in silicon and its deviation from the ideal structure.

    J Appl Phys 2007, 101:063503. Publisher Full Text OpenURL

  27. Lin VS-Y, Motesharei K, Dancil K-PS, Sailor MJ, Ghadiri MR: A porous silicon-based optical interferometric biosensor.

    Science 1997, 278:840-843. PubMed Abstract | Publisher Full Text OpenURL

  28. Agarwal V, del Río JA: Tailoring the photonic band gap of a porous silicon dielectric mirror.

    Appl Phys Lett 2003, 82:1512-1514. Publisher Full Text OpenURL

  29. Palestino G, Agarwal V, Aulombard R, Pérez E, Gergely C: Biosensing and protein fluorescence enhancement by functionalized porous silicon devices.

    Langmuir 2008, 24:13765-13771. PubMed Abstract | Publisher Full Text OpenURL

  30. Escorcia-Garcia J, Agarwal V, Parmananda P: Noise mediated regularity of porous silicon nanostructures.

    Appl Phys Lett 2009, 94:133103. Publisher Full Text OpenURL

  31. Estevez JO, Arriaga J, Méndez-Blas A, Robles-Chaírez MG, Contreras-Solorio DA: Experimental realization of the porous silicon optical multilayers based on the 1-s sequence.

    J Appl Phys 2012, 111:013103. Publisher Full Text OpenURL

  32. Qin Q, Lu H, Zhu SN, Yuan CS, Zhu YY, Ming NB: Resonance transmission modes in dual-periodical dielectric multilayer films.

    Appl Phys Lett 2003, 82:4654-4656. Publisher Full Text OpenURL

  33. Yamilov AG, Herrera MR, Bertino MF: Slow-light effect in dual-periodic photonic lattice.

    J Opt Soc Am B 2008, 25:599-608. Publisher Full Text OpenURL

  34. Pérez K, Estevez JO, Méndez-Blas A, Arriaga J: Localized defect modes in dual-periodical multilayer structures based on porous silicon.

    J Opt Soc Am: B 2012, 29:538-542. Publisher Full Text OpenURL

  35. Yeh P: Optical Waves in Layered Media. New Jersey: Wiley VCH; 1998. OpenURL

  36. Arriaga J, Saldaña XI: Band structure and reflectivity of omnidirectional Si-based mirrors with a Gaussian profile refractive index.

    J Appl Phys 2006, 100:044911. Publisher Full Text OpenURL

  37. Pap AE, Kordás K, Vähäkangas J, Uusimäki A, Leppävuori S, Pilon L, ans Szatmári S: Optical properties of porous silicon. Part III: comparison of experimental and theoretical results.

    Opt Mater 2006, 28:506-513. Publisher Full Text OpenURL

  38. Khardani M, Bouaïcha M: Bruggeman effective medium approach for modelling optical properties of porous silicon: comparison with experiment.

    Phys Stat Sol (C) 2007, 4:1986-1990. Publisher Full Text OpenURL

  39. Bosch S, Ferré-Borrull J, Sancho-Parramon J: A general-purpose software for optical characterization of thin films: specific features for microelectronic applications.

    Solid-State Electronics 2001, 45:703-709. Publisher Full Text OpenURL