Open access peer-reviewed chapter

Macroporous Silicon for Gas Detection

By Didac Vega and Ángel Rodríguez

Submitted: November 8th 2017Reviewed: March 12th 2018Published: November 5th 2018

DOI: 10.5772/intechopen.76439

Downloaded: 748


Macroporous silicon (MPS) has been shown to be a promising material in many areas of technical interest. In particular, MPS has been applied for electronic devices and microfluidic applications. One of the most promising features of MPS is that it enables the development of optical applications using simple and cost-effective technology, compatible with MEMS fabrication processes and suitable for mass production. This chapter describes the application of MPS structures fabricated using electrochemical etching (EE) for the detection of gases of environmental concern in the wavelength range comprising 4 μ m to 15 μ m , such as C O 2 . Vertical-modulated MPS structures are reported, whose photonic bandgaps can be placed at different wavelengths depending on the application needs. These structures have been applied to the quantification of C O 2 , and these results are summarised here. Detection is performed by the direct measure of absorption, obtaining promising results with short optical paths.


  • macroporous silicon
  • photonic crystal
  • electrochemical etching
  • non-dispersive infrared
  • gas sensing

1. Introduction

Macroporous silicon (MPS) has been shown to be a versatile material with a broad spectrum of promising applications [1]. MPS was first described by Lehmann in the early 1990s [2, 3, 4] and has since attracted great interest among researchers. Of the initial works in MPS development, it is also worth mentioning those from Zhang [5], Propst [6], and Parkhutik [7]. Of particular interest is the orderedform of MPS. In such form, pores are etched in a silicon substrate arranged in a uniform pattern. This results in an engineered material whose dielectric constant εxyzhas now a characteristic spatial distribution. Given a regular, periodic εdistribution, the propagation of an electromagnetic wave through the media is affected. In particular, with an adequately designed pattern, it is possible to create an interference pattern giving rise to photonic bandgaps(PBGs) and other optical characteristics for the specific structure. Such structures are thus termed photonic crystals1 (PCs) [8]. A PBG defines a specific frequency range and, possibly, direction in the crystal, in which an electromagnetic wave propagation is forbidden through the material [9].

Thanks to the existence of PBGs, photonic crystals have been used for advanced applications in optical communications [10], photovoltaics [11], photonics [12], light emission [13], anti-reflection/blackbody emitters [14], and gas sensing [15, 16]. Particularly, in this field, the peculiar functionalities of PCs make these structures very attractive for their use in chemical sensing: gas detection [15, 17, 18] and bio-sensing [19]. The main advantage of using PCs for these areas is their potential to design and fabricate very compact and cheap photonic devices [20, 21]. Furthermore, the possibility of integration into large-scale circuits or microelectronic fabrication processes [21, 22] opens up vast opportunities for novel devices.

In particular, in this chapter, we focus on the use of MPS PCs for the optical detection and quantification of gases. More specifically, the characteristic absorption lines, or bands, of a gas in the infrared—that is, the fingerprintof the gas—are used to quantify the presence of the gas of interest. MPS demonstrates that is a suitable material and technology to operate in the mid-infrared (MIR), and is useful for the detection of environmental concern gases, such as carbon dioxide (CO2), methane (CH4) or nitrogen oxides (NOx).


2. Sensing application: gas detection

Sensor demand for everyday applications is rapidly growing. The areas of use are many and multidisciplinary [23, 24]—to name a few: environmental [25, 26], safety [27], security [25], health, transportation, and wearables. Market research shows that gas sensor segment has strong growth [24] (forecast to achieve $765 M in 2020 [23]). Furthermore, climate change concerns is making governments and other agencies push the research and deployment of smart sensor gridsfor large area monitoring [28], resulting in new opportunities for emerging technologies such as ultra-compact PC sensors.

2.1. Detection strategies: optical

Optical gas detection provides very desirable advantages over other methods. In first place, pure optical methods like spectroscopy have exceptionally fast response times to changes in mixture concentration. Furthermore, the spectroscopic optical systems are highly selective: they permit identifying a target gas by its spectral signatureor fingerprint2. Optical systems can be made extremely sensitive, though this generally requires very long optical paths.

Traditional optical-based measuring is based on the direct measurement of optical power. One of its main drawbacks is that the equipment is large and expensive. Indeed, one has to trade-off space for detection limit. Some applications require several centimetres or even metres of optical path to reach the required sensitivity. Furthermore, spectroscopic systems are energy limited, thus they require stable and high power light sources such as lasers or thermal radiators. Spectroscopic systems also require additional complex mechanical and electronic equipment for the signal conditioning and processing. On the other hand, non-spectroscopic optical systems impose less strict requirements on some of these aspects. However, reducing the optical path length has not found a good solution until recently with the advent of PCs.

Other optical systems exist that use alternative detection strategies. Nevertheless, these alternative systems lose some of the more desirable traits of optical detection, like response time and selectivity. Special mention has to be done with respect to terahertz gas detection. This technology is closely related to IR optical detection: typical working frequencies range from 0.1to 10THzwhich correspond to wavelengths from 3mmto 30μm[29]. At these frequencies, gas molecules have mainly vibrational resonances which are then used in the same manner as for IR optical detection: by absorption measurement. Terahertz technology has been successfully applied for the detection of air pollutants and health applications [30] and long range detection. However, terahertz technology has some serious drawbacks: the principal being the generation of terahertz waves [29].For commercial and industrial applications, non-dispersive infrared detectionseems to be the most promising optical solution. This is discussed latter in more detail.

2.2. Gases of environmental concern

The growing global consciousness in environmental preservation and climate change has driven the research and development of sensing devices. In particular, monitoring the environment for pollution control [25, 31] is one of the most important applications. Gas sensors are also significant for health [32] and indoor air quality assessment [33].

Air pollutants and greenhouse gases are primarily related to the exhaust gases of combustion processes. Major air pollutants are carbon monoxide (CO), ground-level ozone (O3), nitrogen dioxide (NO2), and sulphur dioxide (SO2). Greenhouse gases are carbon dioxide (CO2), methane (CH4), and nitrous oxide (N2O) [34]. There are many other toxic air pollutants such as nitric oxide (NO), ammonia (NH3), and hydrocarbons.

These gases have simple molecular composition with strong light absorption in the medium infrared wavelength range. The absorption spectrum of a gas at these light frequencies is caused by the different vibrational and rotational modes of the atomic bonds in the molecules [35]. The absorption coefficient spectra for several of the environmental concern gases are plotted in Figure 1. Absorption data of several gases can be found on the freely available HiTRAN database [35]. Each vibrational-rotational mode of a gas corresponds to an absorption line. These modes are narrow and close together3. However, at atmospheric pressure line mixingoccurs, and the instrument resolution further broadens the spectra, transforming the absorption-lines into an absorption regionor band, as seen in Figure 1.

Figure 1.

Absorption coefficients of some of the most relevant health and environmental concern gases. The coefficients are calculated from HiTRAN data using anapparatus functionofsinc2xwith an aperture resolutionΔf=4cm1.

From the data shown in Figure 1, it is clear that CO2is potentially the largest contributor to greenhouse effect. In addition, its absorption band is clearly separated from other strong absorbers, making this gas easy to detect and quantify in an unknown mixture.

2.3. Non-dispersive infrared detection of gases

Non-dispersive infrared (NDIR) detection relies on the fact that certain absorption lines of gases are “isolated” and their wavelengths have little overlap with other gases, as seen in Figure 1. This is profited to simplify the design of a gas detector reducing cost, complexity, space requirements, and power [36]. Identification and quantification of a certain gas can be done looking at a narrowband region of the spectrum. This can be done using optical filters or selective light sources, or a combination. At the other end, the detector (a photodiode: PD) will give a measure of the optical power received, directly corresponding to the concentration of the gas. When light passes through a gas mixture, some wavelengths will be absorbed following the Beer–Lambert (B-L) law


where I0is the optical intensity of the source (or the reference value), Iis the received optical power, lis the optical path, cis the gas concentration, and ais the absorption coefficient of the gas; source intensity and absorption coefficient depend on wavelength λ. An optical filter thus selects the absorption lines to measure. This allows for a multigas detection scheme using different filters for each gas [37]. Alternatively, the absorption can be measured as follows: A=1T=1I/I0(where Ais the absorption, and Tthe transmission). This is valid under the assumption that no light reflection occurs at the gas interfaces4.

The simplest NDIR system is depicted in Figure 2a. A sample of the unknown gas mixture is placed in a gas cell. Light is then passed through the cell from a source to a detector. From Eq. (1), it is clear that optical path length lis critical for the performance of the system. The necessary optical length lto achieve a desired detection limit can be several meters. Such long paths are highly impractical. Typical commercial products employ special cells (e.g., a Harriot cell) or optical systems using several mirrors to extend the effective optical path keeping a small volume.

Figure 2.

Conventional optical gas sensor using mirrors to attain a long optical path. A PC-based gas sensor can be made much more compact exploiting the special properties of PCs such as slow light. Reprinted from [18], with permission of SPIE.

To remedy the need of complex, bulky, and fragile optical systems to achieve the long optical paths required, a PC can be used as proposed in Figure 2b. PCs-based sensors are projected to require very small footprints, a few centimetres at most. The idea is to take advantage of the special features that PCs exhibit. One of such is the existence of propagation modes with extremely slow group velocity vg(also called slow lightmodes,) such that interaction time of the propagating electromagnetic wave is increased. Another possibility is the inclusion of defects in the PC to create resonant states. These are localised [38, 39] and can have very high-quality factors (Q-factor) prolonging the interaction time of light with the gas. Other alternative sensing methods can be used where PCs can enhance the system response, e.g., photothermal [40, 41] or refractive index change. On the other hand, PCs can be simply used as highly selective filters.


3. Macroporous silicon applied to gas detection

Macroporous silicon is a structuredmaterial which has some very desirable properties. In particular, MPS has been extensively applied for the fabrication of PCs. In the first place, silicon is widely used in the industry, and many techniques and tooling exist. Additionally, MPS can benefit of established fabrication processes such as those used for micro electro-mechanical systems (MEMS). For dedicated applications, there are available fabrication technologies suitable for mass production, like electrochemical etching (EE). Given these advantageous traits, MPS has been proposed as a revolutionary material in many application areas. One of such areas is gas sensing. There has been some very intense work in order to develop the technology of MPS as gas sensing devices or as enhancers.

3.1. Seminal works

One of the first uses of MPS for gas sensing is in the work by Geppert [18]. From the analysis of the photonic band structure5 of the PCs, they observe that some bands are very flat at certain wave vectors, thus the group velocity vgreduces to values about 0.02cto 0.05c(c, the speed of light). The slow propagating light has more time to interact with the media making possible to reduce the physical path length while maintaining an acceptable interaction time. It is also observed that some bands have the maximum field intensity located in the high dielectric region (dielectric bands), while other bands have the maximum in the low εregion (air bands). Thus, air bands must be used for the maximum interaction.

Some of the described devices are three-dimensional (3-d) PCs made of MPS. For these devices, light is shone normal to the surface, propagating parallel to the pores’ axis. The devices were used to detect ammonia (NH3) and sulphur hexafluoride (SF6). Transmitted light was measured through the PC and a small cavity of 300μmdepth filled with gas. They show promising results, but they argue that the small figures obtained are due to the low coupling efficiency of the PC. Furthermore, no enhancement was observed by the presence of the PC. This is also attributed to the large change in the effective refractive index6 (neff) as for high contrast materials the transmittance T1/neff, thus reducing the coupling efficiency. In [18], it is also discussed a simple way to compensate for deviations in the bandgap due to fabrication tolerances: tilting the sample 15°, a shift as large as 35cm1of the PBG (about 4%) can be achieved.

They also present an alternative structure, refined in later works [17, 42], which uses two-dimensional (2-d) MPS structures. These are easier to fabricate and tolerances are better. As the structure is two-dimensional, the photonic band structure exists only for in-planepropagation. Thus, light must be coupled to the sideof the PC. An example of such arrangement is shown in Figure 3a. Additionally, small vgimplies large neff, diminishing coupling efficiency. To avoid this problem, a special taperedsection is placed on the interfaces of the PC. This tapered section gradually adapts the air nairto the neffof the mode. This section can be done by deforming the lattice (elongating one axis) [18] or by leaving additional bulk material at the interfaces (promoting surface waves) [17]. The measurement results of these samples, see Figure 3b, show that in this case the PC does enhance moderately the detection of CO2(assumed 0to 100%change). The devices have a 15dB/mminsertion loss, and variability in the enhancement factor is high. These are attributed to uncertainties during fabrication: a fluctuation of 1%of the pore diameter is reported.

Figure 3.

Optical setup for the measurement of gas using 2-d a MPS PC (a). (b) Shows results of measurement ofCO2. The top panel shows the absolute values, while the lower panel shows the relative variation of the measured value. Cycles 2 and 4 correspond to the presence ofCO2, while the other cycles are with reference atmosphere. Reprinted from [17], with the permission of AIP Publishing.

3.2. Structure design

From the possible strategies when using MPS structures as gas sensors, 2-d crystals impose the need to inject the light from the sides. In this way is easy to obtain long optical paths along the samples, but light coupling is complicated. However, from a practical standpoint, a system based on the normal incidence on the sample, using 3-d PCs, is easier to assemble, align, and calibrate7. In this chapter we expose results based in this option with PCs produced with macroporous silicon.

The designed MPS silicon structures are simulated using a simplified model by the finite-difference time-domain (FDTD) numerical method. Some MPS structures will be used as reflectors, while others will be used as filters. Reflection design requires a PBG encompassing just the absorption region of the desired gas. On the other hand, transmission design needs a PGB as wide as possible, with some crystal defects to block all light but the corresponding to the gas absorption line. The largest PBGs are obtained with an opal like PC [43]. As the samples will be illuminated from the top8, the horizontal pore arrangement is not as important. Our samples pore disposition is a square array. The EE of silicon does not allow getting perfect spherical shapes (opals), but the actual pore profile has a sinusoidal-like shape, or skewedsinusoidal shape as seen in Figure 4. For the simulations, this profile is approximated by a cylinder and an ellipsoid, or by cylindrical slices, see Figure 4.

Figure 4.

Sample model. The unitary cell can be approximated by a “ball and stick” model consisting of a cylinder and an ellipsoid, or by a slice method, using short cylinders.

To place the PBG at the CO2absorption line, the vertical period is calculated to be around 700nm. As a rule of thumb, in a silicon PC, optical bandgaps can be found at about four times the pitch. Therefore, to have a PBG at λ0=4.25μm, l01.1μm. Considering the MPS porosity9 p¯is about 40%, the effective refractive index neff=nSip¯=1.48[44], where nSi=3.47[45]. This gives that the vertical period has to be lSi=l0/neff=718.0nm.

The FDTD simulations show that in the vertical direction there is a narrow PBG (7.1%). A chirpedmodulation smoothly varying the vertical period can double the PBG width [46]. For the samples with included defects, initially the photonic response for a single defect was studied. As a further step, the use of several defects has also been analysed. In particular, the coupling between resonant cavities creates additional resonance modes [47]. Using two cavities, the fundamental resonant mode splits into two, making it possible to filter two separate frequencies using a simple PC design, to perform a differential measurement. One of the resonances selects the gas absorption line, while the other is placed in a “safe” wavelength were it will not be influenced by the media. Such configuration would allow ultra-compact filters or selective emitters.

3.3. Fabrication

For MPS fabrication, the EE method is one of the most versatile, and one of the preferred for three-dimensional structure definition. A general description of the method is the dissolution of silicon by an electrically promoted reduction–oxidation chemical process (redox reaction). The particular implementation of EE used for the fabrication of orderedMPS uses a pre-treatment of the silicon substrate to define the sites (nucleation points) where localised etching of the material occurs. This method was first described by Lehmann in the early 1990s [3, 4]. Succinctly, the silicon substrate to be processed is placed in a bath of hydrofluoric acid (HF). Silicon is unreactive in normal conditions with HF; however, passing a current activates the redox reaction [48]. As detailed in previous works, the reaction exchanges holesat the semiconductor-HFinterface. Either n-type or p-type doped silicon substrates can be used; however, for the controlled growth of macro-pores10, n-type substrates are advantageous, as holes are minority carriers. To etch the substrate holes must then be generated somehow. The most common way is by using photogeneration(light assisted EE). This provides great control on the etching process.

The etching setup places the patterned face to be etched in contact with the electrolyte, opposing to the cathode. The back-face (i.e., the other side) is contacted electrically to close the electrochemical circuit. The illumination can be provided from either the front or back surfaces. For practicality reasons, it is simplest to irradiate from the backside. Our etching setup follows these principles, plus the IR light is provided by IR LEDS with λLED=950nm, and HFis continuously pumped through the etching cell at a controlled temperature.

The pore growth and morphology depend on several factors of the etching process. The pore front (i.e., the tipof the pores) advances at a speed (etching velocity) relatively unaffected by the electrical conditions; however, the etch rate (i.e., the amount of dissolved silicon) is proportional to the electrical current flowing through the electrochemical cell. Thanks to this behaviour, it is possible to use the EE method to fabricate 3-d MPS structures by the modulation of the current and potential during the etching process [49]. More detailed information on the etching process can be found in the works of Lehmann [2], Zhang [5], and Kolasinski [48].

The MPS photonic crystals developed in this work have all been fabricated using the light assisted EE of silicon using n-type substrates. Having pores arranged in an orderedfashion generates internal interference patterns for an electromagnetic wave propagating through the porosified media. These interference patterns will depend on the direction of propagation and frequency; hence, the materialisation of a PC. To arrange the pores in a desired pattern, the method by Lehmann [4] is used. This method prepares the wafers before etching using lithography to define the pore sites. Then nucleation points are created by a short silicon anisotropic etching at the sites. This creates small pyramidal pits in the surface, which promote the electrochemical etching at these spots. EE also allows the formation of very large aspect-ratio (AR, ratio of pore diameter to pore length) pores in relatively short time. For example, AR=1:465has been reported [12].

After the porosification, some samples were post-processed to create a membrane. The membrane is done by anisotropic etching of the back face. For gas measurements, an open membrane is the best option: gas can flow freely through the MPS sample, solving any issues with gas trapping or residence time for the sensors.

3.3.1. Fabrication quality

In general, the EE of silicon produces good quality porous silicon. What is “good quality” with respect to MPS? In a minimal sense, macroporous silicon is of “good quality” if pores are of the same shape and have grown uniformly. Of course this does not provide much information whether the “quality” is acceptable for the intended application. In particular, for photonic applications, the requirements can be very strict. Fabrication imperfections arise due to numerous reasons: wafer crystalline defects, crystal alignment, and local dopant distribution; but also the etching process itself can account for some variability, and lithography errors will also cause flaws in the grown pores. The common fabrication defects one can encounter in a MPS structure fabricated by EE are summarised in Table 1.

Table 1.

Defects found in macroporous silicon.

From Table 1, the defects that are of greater concern are the ones classified as “large optical influence.” Here surface roughnessis not considering as a defect. It is a consequence of minute variations inherent in any fabrication process. The amount of perturbation does not alter the pore morphology or cause long-range effects and will not give unexpected results.

Pore deathis the premature ending of a pore resulting in a shortened pore with respect to its neighbours. This can be single, isolated pores or can affect several pores. Current starvationis the most common cause. During etching, there is a minimal current under which there are not enough carriers to sustain stable growth. In this situation, current will concentrate on some pores resulting in wider, distorted, and shorter pores. The solution is increasing the etching current. Also, shorter structures are less susceptible to suffer from pore death.

Pore branchingis generally a consequence of pore death but can also happen if the etching parameters are not adequate: if applied potential is too small, the space charge region (SCR) will be small and unintended holes could find their way to spots in the pore wall far from the pore tip and secondary pores grow from that points.

Pore axis alignmentdepends on the crystalline direction of the Si substrate. Pores preferentially develop following the principal crystalline directions of the silicon substrate; however, the pores also have a tendency to grow following the current paths. If the crystal is severely misaligned, the mismatch between the crystal direction and the current direction will cause that the pores grow in several directions simultaneously [50]. The pore cross-sectionaffects the symmetry of the PC. In the stable growth phase, cross-section approximates a circle; however, under certain condition of applied potential and electrolyte additives, the cross section can change its shape.

Spikingis the appearance of random streaks of nanometric branches, spreading out in random directions and with irregular shape. Excessive voltage during the etching promotes spiking. Moderate to low spiking has little influence on the optical performance. However, it changes the porosity, which may shift the wavelengths of the optical features. Unstable growthis basically caused by substrate inhomogeneity, though the dynamic behaviour of the etching system may also influence [51]. In general, unstable growth is mild and can be observed with pores not evolving in a straight line but slightly wandering. Another form of unstable growth is the long-rangevariation of pore diameter.

Profile fidelityrefers to the closeness of the actual pore profile to the intended profile. It is considered a defect because it is an unintended result which affects severelythe optical response of the etched sample. Several test runs with progressive refinements are needed before the near-optimal conditions and waveforms are found. After EE, the pore’s walls may be covered with a layer of microporous silicon. This layer will depend on the electrolyte concentration and present additives. It has no noticeable optical influence in the MIR, and it can be left; however, it is very easy to remove by a quick dip in HF.

The abovementioned discussion was made considering the local effect of perturbations in the pore shape and outcome of the etching process. Nonetheless, the uncertainties of etching, the inhomogeneity of the substrate, and the flaws in lithography can be considered on a larger scale or long-range. Pores at distant locations11 are also affected by the etching system. The most common long-range effect is the variation of pore diameter of up to a few percent [42]. Also, the etching process generates H2bubbles which must be removed otherwise they adhere to the surface and create bigspots where the etching stops. More in-depth analysis can be found in the works by Föll and co-workers [52, 53].

3.4. Absorption and losses of macroporous silicon

One important concern is determining how absorption and fabrication tolerances will affect the performance of the detection system. Intrinsic silicon absorbs light at wavelengths shorter than 950nm, approximately, however when silicon is doped with impurities, strong absorption may be observed in the MIR range. Several papers report on the absorption dependence with doping level in silicon [54]. Interestingly, for the doping concentrations used in the fabrication of optical gas sensors working in the middle infrared (n<5·1017cm3), it is found that material absorption is negligible [54, 55].

On the other hand, fabrication tolerances definitively do have an impact on optical performance. For instance, in [17], it is reported that a 1%variation of the pore radius attenuates 15dB/cmthe transmitted light. This is further studied in [42], global pore diameter fluctuation of 3.5%was caused by the spatial variation of dopant. They further claim that, for a better than 90%transmission in a 1mmthick PC, the pore position should change less than 0.3%and diameter has to be within 0.5%. A systematic study of the effect of in-plane disturbances in electrochemically etched MPS is presented in [56]. Some interesting conclusions extracted from this work are that ellipsoid perturbation up to ±10%has little impact on the PBG if the porosity is maintained near its optimum, and that a small variation of the air fraction of just ±4%is enough to shift a narrow PBG (Δf=5%) out of the designed central gap frequency f0. Similar conclusion was found in the work by our group [57], where the effect of vertical period variation in 3-d MPS was analysed. Wide PBGs were designed for the samples used in this work, but noticeable narrowing and shift could be observed with periodicity variations above ±5%. Furthermore, it was found that the PBG became more transparent allowing light transmission higher than a 10%in the forbidden band.


4. Measurements and results

Samples of MPs photonic crystals were fabricated with a lattice pitch of 700 nm and modulated pore profile. The fabrication conditions were the ones described earlier in Section 3.3—using the EE method. The modulation profiles for the 3-d PCs were programmed to generate a “strong”12 profile having a 700nmvertical period, the same as the lattice pitch. The samples were initially characterised to obtain their optical spectrum in the MIR range using a Vertex 70 Fourier-Transform Infrared Spectrometer (FT-IR) from Bruker Optics. Also some of the structures were cleaved and latter inspected by SEM to determine the pore morphology and actual etched profile. After the samples were analysed, the valid samples were used in a dedicated gas measuring system also using FT-IR spectrometers to measure the performance of the PCs in sensing carbon dioxide.

4.1. Fabrication, sample characterisation, and morphological study

The fabricated samples were etched using a silicon substrate with 001crystal orientation. For the 700 nm period MPS structures, substrates of the appropriate resistivity were used. The pores are arranged in a square pattern. This pattern is transferred to the substrate surface using nano-imprint lithography (NIL). The etching temperature was set to 10°C. The resulting PCs are to be used with light coming from the top, aligned to the pores’ axis, with normal or quasi-normal angle of incidence (θifrom 0°to 13°). The actual modulation waveforms depend whether the samples will be used for reflectance or transmittance. For the samples used in transmission gas measurements, membranes were made. However, the structure resilience is critical given the thin porous layer, so membranes are notopen. Indeed, the porosified depth of the MPS crystals ranges from 20to 50μm. Such thin layers are very fragile. For this reason, membranes did not reach the pore tips.

4.1.1. Morphological analysis

The morphological analysis of some samples shows, as seen in Figure 5, that the etched profiles of the MPS structures have some imperfections. The observed flaws are mainly small variations in pore radius between adjacent pores, pore “wiggling”, skew, and pore death. On the surface of the pores, a rough finish due to microporous silicon can be appreciated (see Figure 5). Pore radius variation is small—less than 10%—for the samples used here. This has a small impact on the optical response of the PC as has demonstrated by our works [57] and others [58, 59]. Skew is generally not an issue. Some extreme cases have shown up to 5°of crystal misalignment from the 001direction (see Figure 5b), but otherwise pores grow fine. Some pore “wiggling” can be observed in the fabricated samples (see Figure 5a). As a consequence, there is some decrease of the optical performance of the PC. In general, the degree of wiggling in our samples is small.

Figure 5.

Cross section detail of a MPS photonic crystal (a) and (b), and the current waveform used to generate (a).

The fabrication process for the MPS samples has been optimised to obtain structures virtually free from dead pores. In general, optical performance is preserved if dead pores are few, or if pores die late during etching. To avoid pore death in the initial samples, the lower current limit was increased, so dynamic range was reduced and modulation could not be as strong as desired (compare Figure 5a and b).

The most significant defect found in the samples here used was vertical period variation13. As seen from Figure 5, for neighbouring pores, the beginning and ending of each vertical period is slightly different. As shown in our previous work [57], this variation has a noticeable influence on the PC optical response. In general, the 700 nm samples fabricated with our equipment show about 5%variation in periodicity which results in about 30%PBG narrowing from the ideal and nearly 10%minimum transmittance [57]14.

Two types of PC structures were fabricated for this work. Samples used for reflectance measurements are regular 3-d structures with a continuous sinusoidal-like profile for the pores. Such samples are shown in Figure 5 along with an example of current waveform used to generate them. The cross section images reveal that the modulation is slightly skewed to the beginning of the modulation. The vertical period is measured to be 700nmwhich is the desired length. Total etched depth is 25μmand the mean pore diameter is d¯=500μm. The programmed “modulation index” 15 is about min=0.7, but results in metch0.1. Despite this, this structure has good enough optical performance.

Comparing the etched pore profile from the input waveform (Figure 5c), it can be seen that during the high plateau portion of the profile, once the pore has reached its maximum diameter, it slowly starts narrowing as the pore front advances. Presently, this can only be “corrected” by trial and error and judicious changes to the current and potential waveforms. For example, a second profile was designed with smoother transitions and “pre-skewing” resulting however in the PC of Figure 5b. Better modulation index was obtained (metch0.2) but the skewing did not improve. In addition, the vertical period increased to 850nm, which corresponds to a PBG centred at λ0=5.03μm.

Samples used for transmittance measurements include one or more defects, as shown in Figure 6. These structures have been fabricated using a refined profile. As seen on the SEM cross-section micrograph, the modulation skew is still present16, but the vertical period length is better adjusted. Otherwise, the quality of the fabricated structures is maintained. The modulation index has been increased, with pores having a smaller diameter at the necks. This greatly improves the PBG (both width and transmittance blocking) of the PC compared to the first structures depicted in Figure 5. The vertical period of the PC used for these samples is p¯v1μm. These samples have at least one planar defect in the regular PC profile to define a resonant cavity. These cavities have a diameter ddef=230nm, and length was varied around ldef=2.1μm[60], to place the resonance at λ=4.25μm.

Figure 6.

Macroporous silicon structures used in transmission measurements. (a) Is a single defect PC, while (b) shows a two defect MPS structure. Panel (b) © 2017 IEEE. Reprinted, with permission, from [60].

4.1.2. Optical response

Samples were optically characterised after fabrication to ensure the adequacy to sensing CO2. This was done in a FT-IR spectrometer in reflection mode with light incident at a quasi-normal angle of θi=13°. Also some samples were analysed in transmission mode with light illuminating in normal angle of incidence. The wavelength range of the analysis extends the MIR and some of the NIR regions from λ=1μmto 20μm.

Samples used for reflectance gas measurements, such as the one in Figure 5a, have a single PBG in the optical response. The morphological analysis shows that, as the modulation index is small, the PBG will be relatively narrow. This is confirmed in the measurement shown in Figure 7. The results confirm that as expected by the criteria given above, the PBG of a 3-d PC with a vertical period of 700nmis centred at λ0=4.25μm. PC quality and modulation index limit the reflectance to about 30%. In spite of these shortcomings, the obtained response is good enough to be used for the sensing of carbon dioxide.

Figure 7.

Measured spectra of the PC shown inFigure 5a. The absorption lines ofCO2are marked with arrows, the length giving a qualitative impression of their strength. It can be clearly seen how the PBG of the MPS structure corresponds to theλ=4.25μmabsorption line.

The samples for transmission measurements have different responses according to the number of PC defects placed. As seen in Figure 8a, b, placing two defects gives rise to two resonant states, and three states if three defects are used. Coupling between resonant cavities induce the appearance of several resonances in the optical spectra even if all cavities are of the same dimensions. The placement of the defect along the PC depth also influences the coupling efficiency and quality factor of the resonance. The creation of a membrane also improves the transmitted signal, compare Figure 8c and d where the transmitted peak is almost 10%to the 4%of the non-membrane sample.

Figure 8.

Characterisation of transmission PC samples. (a) Is a two defect sample which shows two resonance states inside the PBG. (b) Is for a three defect PC (thick line) compared to a single defect crystal. (c) Is the characterisation of the sample inFigure 6b. (d) Shows the response of a two-defect membrane PC with a peak centred in theCO2absorption line at4.25μm, and a reference peak at3.8μm. Panels (a) and (b) © 2017 IEEE. Reprinted, with permission, from [60]. Panels (c) and (d) from [61] are licensed under CC BY-NC-ND 3.0.

4.2. Gas measurement system

After the devices were characterised, they were measured under different atmospheres consisting of CO2diluted in pure nitrogen. The amount of carbon dioxide was controlled using a two gas mixer in which one line is pure nitrogen as the carrier gas, and the process gas (CO2is this case) is connected to the second line. Each line is controllable separately: shutting, pressure, and flow. Both gas lines are then mixed in a manifold and the mixture is then output to the gas cell. The gas cell is sealed with one input port and one exhaust port. The flow control is performed using mass flow controllers (MFCs) of different range: the carrier gas controller (MFC1) is 2l/minfull scale (FS), and the process gas controller (MFC2) is 200ml/minFS. Both MFCs are calibrated for N2, so for the process gas a correction factor is used. For the measurements, a constant flow rate of V̇=400ml/minfor the mixture is used. In these conditions, the minimum gas concentration goes from 7500ppmto 50%. Given the performance data of the MFCs, the concentration uncertainty is given by Bc2c23.15·103+0.25·103c112, where c=V̇2/V̇1+V̇2is the gas concentration, and V̇1and V̇2are the flow rates of carrier and gas, respectively. Note that V̇=V̇1+V̇2.

The gas mixture is fed to a gas cell purposely built where the MPS photonic crystal is placed. The gas flows continuously through the cell. For the measurements, a broadband infrared light source is used. This light is directed to the cell and then collected by a photodetector (PD). In the work here presented, the MPS structures were characterised by FT-IR spectrometry. The expected response of an autonomous NDIR measuring device can be then extrapolated from the spectroscopic gathered data in the characterisation of the PCs.

4.2.1. Principle of detection and method

The basic idea here proposed for an NDIR system is to use a MPS photonic crystal with a PBG wide enough to comprise one of the gas absorption regions. In the case of CO2, this region extends from λ4.20to 4.35μm(see Section 2 – above). Ideally, all light outside the PBG should be discarded and not arrive at the detector. Once the light has been filtered, the photodetector gives a reading proportional to the incoming optical power, and further filtered by the detector’s photo-response. In the following discussion, it is assumed that the PD response is flat in the absorption region of the gas.

From the characterisation measurements, it is straightforward to approximate the expected response of the NDIR system. The spectra can be any of the reflectance, transmittance or absorbance, as these are proportional to the power of the electromagnetic wave: Erλ2=RλEiλ2, and so on. Now it is possible to define the quantity SR=Rλdλas the normalised optical power received at the PD from the reflectance data. SRis thus proportional to the gas concentration and estimates the actual output of a PD. Equivalent power quantities STand SAcan be defined for the transmittance and absorbance spectra respectively. A second measuring criteria is simply evaluating the reflectance spectrum at a given wavelength Rλgas=Rgas(equivalent quantities are defined for transmittance and absorbance). Such wavelength is chosen to be where absorption is the greatest. The relation of this quantity with an actual output of a PD is less clear but serves to evaluate the performance of the gas cell and PC.

4.3. Reflection CO2measurement

Reflection measurement of CO2uses a PC as a selective reflector. The reflection spectra are obtained from which the absorbance spectra are calculated as Aλ=1Tλ=1RλR0λ; where Rλis the measured spectrum at some gas concentration, and R0λis the reference spectrum (only N2). The “optical power” is then S=Aλdλ. Using a reference spectrum removes any effect of the optical system, such as reflectance at optical interfaces.

4.3.1. Setup

The measurement setup consists of the gas mixer connected to a specifically built gas cell. A schematic view of the cell design is shown in Figure 9 alongside the actual built device. The cell consists of two plates separated with an O-ring such that the cell is made airtight. The top plate has one port for gas input and another for output, and a central opening where a potassium bromide (KBr) window is placed and held in place by a sealant. When assembled, the two plates are held by several screw ensuring mechanical stability and airtightness, with the PC underneath the IR window. The gap from the sample to the window is about 0.5mm, resulting in a total17 light path of 1mm.

Figure 9.

Gas cell for reflection measures. A schematic section view is shown and the actual fabricated cell is pictured next to it. The sample is held in place by locking it with a spring screwed in a suitable position of the cell base grid. The cell is made airtight by an O-ring, and sealant around the IR window.

The cell is then placed in the spectrometer to make the gas concentration readings. In particular, a Renishaw Raman microscope equipped with Smiths IlluminatIR FT-IR module is used to take reflectance measurements at normal incidence, θi=0°. The aperture sizes used in the optical setup was 50μm×50μm. The wavelength range extends from k=1000cm1(10μm) to 3500cm1(2.86μm). Resolution is 4cm1.Each spectrum at a given concentration was averaged over 256scans.

4.3.2. Results

Reflectance spectra were measured for CO2at concentrations from 3to 43%, the results summarised in Figure 10. Lower concentrations were tested but for concentrations below 3%changes are very small and difficult to discriminate. This effect is then observed in the “optical power” signal Sat the output (Figure 10b) calculated as mentioned above. Figure 10b also shows that the Beer–Lambert law is closely followed by the measured data for concentrations above 10%. From the fitted B-L model, it is possible to extract the exponential term l·a3.65(where lis the optical path length, and athe absorption coefficient). From the HiTran data for carbon dioxide, at λ=4.25μmand 4cm1instrument resolution, aCO250cm1. This gives that the effective path length is leff0.75mm. Measurement uncertainty is dominated by the mixer uncertainty, as seen in Figure 10b, as power measures are averaged over a large number of scans. This results in a fairly constant uncertainty across the full measurement range about ±2%.

Figure 10.

Measurement results for reflectance setup. The spectral response (a) fromk=1000cm1(10μm) to3500cm1(2.86μm) has been normalised to a pureN2atmosphere and shows the effect of absorption due toCO2in theλ=4.25μmband. The output response (b) against concentration shows good sensitivity for concentrations above5%. The dashed line is the ideal Beer–Lambert response.


5. Conclusion

Macroporous silicon is a versatile material that has shown to be a good candidate for the obtainment of PCs for gas sensing applications. Using fabrication methods such as EE opens up the possibility of obtaining high-quality photonic filters, in large quantities, and economically competitive. Furthermore, this fabrication technique is very flexible allowing creating customised designs with little effort. Using silicon as the base material has other benefits, such as the reutilisation of the existing manufacturing tooling and the reuse of process flows. Moreover, EE is compatible with microfabrication technology and might be incorporated in VLSI designs to build complete sensing devices. This will result in more compact and integrated system design lowering the bill of materials, costs, and improving manufacturability.

Here we demonstrate one possible way to use MPS PCs for gas sensing: as selective filters. Carbon dioxide has been detected and quantified using NDIR reflectance measurements. It has been found the MPS crystal has also an effect in the measured absorption. This is due to the very nature of the PC, slowing the group velocity of the incident light and enhancing the interaction time—increasing the effective optical path length—with the gas mixture. The inclusion of resonant cavities further enhances light absorption by inducing resonant states and spatially confining the IR radiation.

The PCs here shown, prove that a compact CO2sensor using MPS technology is feasible, achieving a detection range near to that of commercial optical devices based on IR PD/LED. The devices shown here have room for improvement in particular regarding the fabrication, and progress is steadily being made in this area.



This work has been funded by the Spanish Ministerio de Economía y Competitividadwith research grants TEC2010-18222 and TEC2013-48147-C6-2-R. The authors would also like to acknowledge the contribution and help from Daniel Segura García and David Cardador Maza in the preparation of this chapter.


Conflict of interest

The authors declare that there is no conflict of interest.


  • Nonetheless, MPS is not the only material nor technology in which PCs can be devised and fabricated.
  • The chemical composition of a gas compound has specific vibrational (atomic bonds; from MIR to VIS) and electronic (electron excitation; from VIS to X-Ray) resonances that result in particular absorption frequencies unique for such gas.
  • An ideal vibration-rotation mode has a single frequency. However, the actual line profile of a mode depends on external factors such as the gas pressure, velocity, temperature, etc. Typical line profiles are Gaussian (Doppler broadening), Lorentzian (pressure broadening), and Voigt (mixture of the previous). In standard conditions, ideal gas line profile separation is around Δ k = 2 cm − 1 , and full-width half-maximum approximately FWHM = 0.2 cm − 1 .
  • The construction of a gas cell as well as other optical elements existing in the light path introduces several reflections. However, all of them are accounted when measuring the reference value I 0 .
  • The band structure is the reciprocal of the dispersion diagram, thus the slope of the band represents the group velocity v g = ∂ ω / ∂ k , where ω is the wave frequency and k the wave vector in the reciprocal lattice of the PC, for a plane wave propagating through the media [62].
  • There exist several methods to find the effective refractive index of a PC. The simplest one for low frequencies is calculating the average of the different materials that compose the PC. To obtain the n eff at higher frequencies, a theoretical study of the photonic band structure is needed, from it, n eff can be derived from the group velocity.
  • Indeed, using 3-d PCs, the light can be launched into free space and coupled into the crystal from the top or bottom surface, which are much larger than the sides of the MPS structure.
  • Concretely, the “top” surface is the surface from where the pores are etched. This surface has the initial pattern of the pore sites. For prime quality wafers, this face is polished and the incident light will have little scattering.
  • Porosity is defined as the ratio of air to silicon volume in the unit-cell.
  • A porous media is termed microporous for pore diameters less than 5 nm , mesoporous for 50 nm > d pore > 5 nm , and macroporous if pores are larger than 50 nm wide.
  • For example, a scratch on the back-surface or shadow of the illumination will cause the shape of such scratch or shadow to be transferred to the pores grown at the front surface: the affected area will have smaller diameter pores or, in extreme cases, dead pores and branching.
  • This profile tries to obtain a spherical shape as close as possible.
  • That is for the current way the PCs are being used: with light coming from the “top” along the pores’ axis.
  • These performance figures were achieved with the latter fabricated samples, with an optimised EE process.
  • The input waveform is a square signal, so the concept of modulated index in AM is extended here as m = max r t ∕ min r t , where r t is the radius of the poresuperfluous.
  • However for this particular instance, the asymmetric modulation was actually designed.
  • This figure may change due to uncertainty. For example MPS sample thickness or tightening of the screws can change the gap as much as 0.2 mm . It must be remarked that this is a proof-of concept cell.

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Didac Vega and Ángel Rodríguez (November 5th 2018). Macroporous Silicon for Gas Detection, Applications of Silicon Photonics in Sensors and Waveguides, Lakshmi Narayana Deepak Kallepalli, IntechOpen, DOI: 10.5772/intechopen.76439. Available from:

chapter statistics

748total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

New Approach to Mach-Zehnder Interferometer (MZI) Cell Based on Silicon Waveguides for Nanophotonic Circuits

By Trung-Thanh Le

Related Book

First chapter

Semiconductor Optical Amplifiers

By M. Haridim, B.I. Lembrikov, Y. Ben-Ezra

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us