Chapter
Properties of Silicon
1.1.4 Electrical Properties
Silicon is a group IV element in the periodic table and is a semiconductor with a bandgap of 1.12 eV, which means that pure silicon at room temperature is almost an insulator. By doping with group III or group V elements, the resistivity of silicon can be varied over a wide range.
1.1.4.1 Introduction—Dopants and Impurities in Silicon
Semiconductors are solid materials that have electrical conductivities in between those of conductors and those of insulators. The physical reason causing a material to behave as a conductor, semiconductor, or insulator lies in the availability, or lack thereof, of free current carriers in the material. Semiconductors are characterized by the narrow bandgap between the valence bands, occupied by electrons, and the conduction band, in which electrons move freely according to applied electrical fields. Intrinsic (i.e., pure) semiconductors act as insulators at room temperatures, but their behavior changes dramatically with temperature, and, more to the point, with small impurities present in the crystal. Very small amounts of electrically active impurities can totally alter the electrical properties of semiconductors such as silicon. This is because the electrically active impurities either easily donate valence electrons (donors) or accept them, creating holes (acceptors). These electrons or holes are free (i.e., not bound to individual atoms). Their movement due to applied electrical fields carries electrical currents, giving rise to the term charge carriers used to denote them.
The electrical properties of semiconductor materials such as single-crystal silicon are thus defined by the impurity concentrations present in the silicon lattice. Impurities are introduced into the starting materials during crystal growth and modified during device processing by additional doping of the silicon material with electrically active impurities. In intentional doping of silicon, impurity atoms from group III (acceptors) and group V (donors) are used. The dopants used in crystal pulling, and the conductivities reached with specific impurity concentrations, are described in more detail in Sections 3.1 and 3.3Section 3.1Section 3.3, Dopants and impurities in silicon crystals.
Manipulation of the electrical properties in the structures created during MEMS device manufacturing follows practices established in semiconductor device manufacturing. Techniques used include both very traditional methods, such as deep diffusions of dopants, which have been abandoned in mainstream semiconductor processes, and current standard techniques such as ion implantation and epitaxial deposition. While semiconductor-grade starting materials are largely free of other electrically active impurities, the incorporation of unintentional contamination into silicon during processing can have significant effects on the electrical properties of the manufactured devices.
Unintentional doping of silicon includes the introduction of unwanted donors or acceptors to the crystal lattice from the processing environment. These impurities can be either group III/V-type atoms misplaced, or other contaminants, such as some transition metals.
The generation of electrically active donors also takes place within the single-crystal CZ silicon itself, without the introduction of additional impurities. CZ silicon always includes a few parts per million atoms of interstitial oxygen atoms, originating from the quartz crucible used to hold the melt during crystal pulling (see also Chapter 3). At certain temperatures, ranging from 400°C to 550°C, these interstitial atoms create conglomerates of several oxygen atoms within the lattice. Such silicon–oxygen microclusters are known as thermal donors (TDs), as they donate free electrons to the conduction band, influencing the electrical properties accordingly [22,23]. The concentration is, however, usually below 1015×cm−3, and thus has only very marginal effects on other than high-resistivity silicon. These donors are not stable at temperatures above 600°C, and even for applications requiring high-resistivity silicon, their effect can be suppressed by quick cooling down over the generation temperature range. This method is called thermal donor anneal (TDA), and while it is effective, it does not prevent the generation of new TD if temperatures in the critical range are used later in the device processing. For further details, see Section 3.6.
The use of semiconductors is based on the fact that the charge carrier concentrations are also influenced by any electric fields that are present. This can take place intentionally, as, for example, in transistors, but also due to electric fields generated by surface effects such as charging. These effects are more pronounced in high-resistivity silicon, but they bear consideration in general.
In the vast majority of cases, the electronics required for the realization of a MEMS-based sensor consist of the basic building blocks of semiconductor electronics used since the birth of the silicon-based semiconductor industry and described in the basic handbooks of the industry, such as Sze [24]. In MEMS applications requiring very high resistivity, in fields such as RF and optical, special considerations apply. The very high resistivity materials used are obviously strongly affected by even the smallest concentrations of unintended charge carriers. These unintended charge carriers can be introduced into the material by methods such as TD generation (described earlier), oxidation, and contamination of the surface if the donors/acceptors are not subsequently evaporated. In very high resistivity silicon, these effects can be severe, in some cases even leading to type reversal.
1.1.4.2 Piezoresistive Effect in Silicon
1.1.4.2.1 General Piezoresistive Effect
The change in the resistance of metal devices due to an applied mechanical load was first discovered by Lord Kelvin in 1857. With the large-scale use of single-crystal silicon in the making of semiconductor circuits, a much stronger piezoresistive effect was discovered in silicon [25]. This discovery forms the basis for practically all piezoresistive MEMS applications.
1.1.4.2.2 Strain
When a sample is subjected to physical force, the force yields a change in length, dL, that follows the well known Hooke’s law F=kΔL, where k is a material constant. The stress is defined as the applied force per unit area. Stress (σ) is thus given by:
where εL denotes the differential deformation, ΔL/L, known as strain. Most solids exhibit elastic behavior for small stress loads. In elastic deformation the strain is proportional to the applied stress and the material follows Hooke’s law, which, in the simple case of uniaxial stress, can be written as:
where the scalar Y stands for the modulus of elasticity, also known as Young’s modulus. This material parameter for crystalline materials such as silicon is not omnidirectional, but varies according to crystal planes.
For strains above certain thresholds, Hooke’s law is no longer valid. This is caused by non-elastic, or irreversible (i.e., plastic) deformation.
The piezoresistive effect is quantified using the gauge factor G, describing the relationship between the applied strain and the change in resistivity. The piezoresistive effect, and thus also the gauge factor of semiconductor materials, is several magnitudes larger than the geometrical effect observed in metals. The effect is observed in several semiconductor materials, such as germanium, polycrystalline silicon, and single-crystal silicon.
The large piezoresistive effect in silicon is due to the fact that, instead of a stress-dependent change in geometry, the effect is first and foremost due to the stress-dependent resistivity of the material. The strain applied to the lattice causes the anisotropic band structure to deform, giving rise to anisotropic change in the mobility of the current carriers.
1.1.4.2.3 Stress in Anisotropic Materials
The elastic behavior of isotropic materials can be defined with just two elastic constants, commonly used parameters being the Young’s modulus Y and Poisson’s ratio ν, but modified constants known as bulk modulus K and shear modulus G are also sometimes used. In the general anisotropic case, the stress state of the material, is defined by a matrix known as the stress tensor, a 3×3 matrix in three dimensions.
Crystallographic symmetry allows the reduction of stress tensors in silicon to six-vector notation. In static equilibrium, the situation is further simplified to just six independent variables, and the stress tensor can be written as [26]:
Typically, the geometries used in MEMS applications reduce the situation, as in the case of a piezoresistive sensing layer created on a relatively thick silicon substrate. For such a layer, the vertical (z) stresses can be equated to zero, and the stress is reduced to:
This expression can be further simplified if we limit our discussion to stresses parallel, and stresses perpendicular, to the direction of current. For practical sensor systems, this limitation can easily be applied.
Thus, the forces can be denoted as σl (parallel stress) and σt (perpendicular or transverse stress).
1.1.4.2.4 Strain Effect on Resistivity
In metals the conduction/valence bands are partially filled with charge carriers, and small changes in band shape do not affect most of the carriers. In semiconductors the shape and dimension of the bandgap has a much larger effect, and shifts in energy bands due to applied stress affect the mobility of charge carriers. The result is resistivity change due to applied stress and gauge factors up to two magnitudes larger than those observed in metals. Because the relative change in resistivity is proportional to strain, it can be written as:
For strained anisotropic material, the resistivity is thus no longer scalar and is described through tensor mathematics, similar to the mathematics described earlier for stresses in anisotropic material. Thus, the dependence of the current density J on the electric field E is denoted as:
Details of the tensor mathematics explaining the anisotropic resistivity are described in, for example, Thomsen [26].
In typical applications the MEMS devices are constructed in such a way that the current in the piezoresistors is either parallel or perpendicular to the direction of the stress. Also, the vast majority of applications are made on (100) silicon wafer substrates, with the substrate oriented in such a way that the resistors are aligned on the [110] direction on p-type and the [100] direction in n-type, as depicted in Figure 1.13. These directions yield the maximum positive and negative piezo coefficients (i.e., the maximum relative resistivity change for a given applied force).

Figure 1.13. Typical piezoresistor alignments used on MEMS-standard (100) silicon wafers, aligned to utilize maximum piezo coefficients for rectangular membranes of the type used in pressure sensors. (a) P-type silicon, Dicing along <110>, (b) P-type silicon, Dicing along <100>, (c) N-type silicon, Dicing along <100>.
For a system in which the coordinates are chosen to coincide with the crystallographic axes, the relationship between the stress tensor and the anisotropic change in resistivity simplifies due to crystallographic symmetries of silicon, and can be defined as:
where Π11, Π12 and Π44 represent three independent piezo coefficients, listed in Table 1.2.
Table 1.2. Independent Piezo Coefficients for Crystalline Bulk Silicon [25,27]
10−11 1/Pa | n-Si ND=4×1014 [25] | p-Si NA 2×1015 [25] | p-Si NA 1×1018 [27] | p-Si NA 1×1019 [27] |
---|---|---|---|---|
Π11 | −102 | +6,6 | ||
Π12 | +53 | −1,1 | ||
Π44 | −14 | +138 | +103 | +81 |
In most applications the change in resistivity is limited to stresses parallel and perpendicular to the resistive path. In such cases, the situation is greatly simplified [28]. The tensor element becomes a scalar equation, and can be written as:
where the common notions of longitudinal piezo coefficient πl and transversal piezo coefficient πt are used. As an example of a common situation, let us consider the case of a p-type piezoresistor constructed to utilize the maximum piezo sensitivity as described in Figure 1.13(a). For such a resistor, which is uniaxially stressed and lying in the (110) direction, the shear strain is zero and the longitudinal piezo coefficient is simplified to:
The corresponding transversal piezo coefficient then becomes
For piezoresistors in the <100> plane of an n-silicon substrate, the maximum piezo coefficient is found along the (100) crystal axis, and the piezo coefficients are reduced to:
and
It is immediately observed that for both n- and p-type piezoresistors, uniaxially stressed at respective optimum locations for sensitivity, the longitudinal and transverse piezo coefficients have opposite signs. The advantages offered by this fact in a Wheatstone bridge have contributed to the dominance of bridge circuits in piezoresistive sensing based on silicon.
The alignment of piezoresistors to be constructed is governed by the piezoresistive coefficients listed in Table 1.2. Almost universally, the optimal solution is to align the sensing resistors according to maximum piezoresistive effect, as described by the examples depicted in Figure 1.13. The choices of resistor direction, dicing direction, and anisotropic etching are determined by the crystal orientations. In addition, the specification of the crystallographic location of the substrate wafer flat is based on the chosen crystal alignment. The primary flat has historically been located at the <110> direction, and dicing along and perpendicular to the primary flat is depicted in Figure 1.13(a) and (c).
1.1.4.2.5 Linearity
In comparison to many other physical phenomena, the piezoresistive effect is a fairly linear effect. To be exact, the linear equations for piezo coefficients should, however, be replaced with higher order polynomials. For resistors aligned on the <110> crystal axis, for maximum piezoresistive effect in p-type silicon, a quadratic correction to the transversal piezoresistive coefficient is sufficient for a wide range of doping concentrations and applied stresses. The quadratic correction for p-type silicon along the <110> crystal axis is of the form [28]:
yielding a quadratic correction of less than 1% for stresses up to 10 Mpa.
1.1.4.2.6 Effect of Temperature and Doping
In addition to the piezoresistive effect, the silicon resistors exhibit strong temperature dependence. For piezoresistors with low to moderate doping, the resistivity changes can reach one magnitude for every 100°C. Change in offset is caused by leakage currents, which are strongly dependent on temperature according to Laermer [29]. While in some cases the temperature dependence can be utilized or ignored, in the vast majority of cases this must be compensated for. Combined with the effect of doping on the piezoresistive effect, the temperature effect can be corrected with the piezoresistance correction factor P(N,T). Generally, a higher dopant concentration in silicon reduces the temperature effect, whereas the higher the temperature, the lower is the effect of doping level. Thus, designing of the resistor requires a trade-off between sensitivity and temperature stability (Figure 1.14).

Figure 1.14. Piezo correction factor P(N,T) in n-Si as a function of doping level, for various temperatures.
Redrawn and modified from Kanda [30].The modern IC technology can provide a mathematical solution to the compensation of the piezoresistors, by introducing signal conditioning of the linearity and temperature effect to the sensing setup. This approach provides additional accuracy and flexibility at the cost of some cost and complexity. While the approach is commonly used to improve nonidealities in surface micromachined sensors, the performance of sensing elements manufactured from bulk, single-crystal silicon is usually close to ideal when temperature effects are accounted for.
Basic solutions for correcting for temperature variations include the use of a reference (i.e., a non-strained resistor which is at the same temperature). The difference of signal between the strained and non-strained sensors provides the desired signal due to strain only. The other basic technique is to utilize a temperature independent circuit—in practice, a bridge circuit.
In a Wheatstone bridge circuit application, the four resistors are doped to the silicon and arranged in such a way on the sensing membrane that they react to strain in pairs. In addition to having the advantages in sensitivity, discussed earlier, the bridge circuit can be used to compensate for temperature effects. The standard setup of a Wheatstone bridge in (100) silicon substrate consists of one pair of resistors aligned to give maximum positive change in resistivity (longitudinally stressed if p-type, transversely stressed if n-type), while the other pair yields the maximum negative change (a perpendicular alignment with respect to the first pair of resistors). This bridge circuit provides an output signal which is independent of temperature, assuming that the whole membrane is subjected to equal temperature changes. The circuit output is thus, to a first-degree approximation, directly proportional to the change in resistivity due to the applied stress. The main disadvantage of piezoresistive gauging using bridge circuits lies in the current requirements in the milliampere range, which cause limitations in the use of these circuits in power-sensitive applications.
1.1.4.2.7 Example of a Piezoresistive Sensor Design
Pressure sensing has utilized piezoresistive sensing based on MEMS technology for years and is the first widely adopted application of the technology. This is due to the fact that thin micromachined silicon membranes with a piezoresistor bridge implemented on top of them provide a rather ideal solution to the common problem of sensing pressure over a known, rather limited range of pressures, over long periods of time and countless cycles of pressure changes, for low cost and no maintenance.
A very basic sensing setup is as follows:
Starting point:
Silicon substrate, N-type, (100) orientation, device dicing along <110> direction.
Membrane defined by etching, using an alkali etchant from the backside of the wafer. The long etching is carried out in concentrated KOH or TMAH, requiring a thick protective layer of oxide, or nitride. The traditional bulk micromachining has utilized backside patterning for these cavity-etching processes.
Cavity sealed through wafer bonding. The sealing options range from glues and glass frits to fusion bonding of silicon wafers.
Piezoresistors are boron implanted on top of the wafer, at the edges of the membrane, where the stress maxims are located when the cavity pressure differs from the outside pressure. Boron implantation of the ohmic contacts can usually also be performed at this stage.
The basic device is finished with one or two layers of metals.
The dopant profile of a simple implanted and diffused piezoresistor of the type described is illustrated in Figure 1.15. The denuded zone responsible for the electrical isolation is shown without bias; with reverse bias, it can be extended to meet the requirements.

Figure 1.15. Basic diffused piezoresistor dopant depth profile. The wafer surface is to the left, and the piezoresistor is hatched. The transition zone under the piezoresistor acts as the electrical isolation. The hatched line (– · · –) on the left denotes free hole concentration, and the hatched line (— · —) on the right side of the transition denotes free electron density, approaching the starting material ND density in the depth of the wafer. The graph is calculated using IC-simulation software ICECREM 4.3 for WINDOWS™ by Fraunhofer Gesellschaft for Integrated Systems and Device Technology (IISB), Erlangen, Germany [31], and the carrier concentrations refer to zero bias.
The basic device outlined earlier follows the geometry outlined in Figure 1.16. This device has the piezoresistors located at the stress maximums of the rectangular membrane and along the crystal axis yielding greatest change in resistivity for a given stress. Also, the shear stresses of the membrane are minimized, theoretically to zero, at the locations of the resistors.

Figure 1.16. Basic MEMS-based sensor application using piezoresistors in <100> silicon. (a) Side view (cross-section along the <110> plane). (b) Top view, with the piezoresistors located in a Wheatstone bridge, along the (110) direction on the top of the silicon membrane.
1.1.4.2.8 Surface Effects
In many cases, the carrier concentration depth profile of the piezoresistor can be very important. Especially in accelerometers, the piezoresistive MEMS-sensor is often loaded so as to cause compressive strain on one side of the membrane and tensile strain on the other. In such a case the effective piezo coefficient, and thus the sensitivity, is very strongly influenced by the ratio of carrier concentrations between the membrane surfaces. Also, those piezoresistors whose stress is purely compressive or tensile are affected by the vertical carrier profiles, as the greatest strains are located at the very surfaces of the membrane. On the other hand, the carrier concentrations tend to change more rapidly the closer they are to the surface. Because the carrier concentrations near the device surfaces can also be affected by charging of the dielectric layers used to isolate and seal the micromachined structures, a careful consideration of the design principles used in analog transistor design is necessary to ensure the reliable operation of the piezoresistive element.