The Young's Modulus of Silicon

A Simple Explanation

By M.A. Hopcroft, Ph.D.

Introduction

The value of the Young's modulus of silicon is often required for engineering designs using micro-electro-mechanical systems (MEMS) technology. However, silicon is an anisotropic crystalline material and so the answer to the question, "What is the Young's modulus of silicon?" is, "it depends, and it can vary by up to 45%". Fortunately, it is straightforward to determine the correct answer for any situation. A condensed answer to the question and explanation is given here.
The information here has been published with significant additional material in a peer-reviewed journal. Please cite this reference when using this information:

M. A. Hopcroft, W. D. Nix, and T. W. Kenny, "What is the Young's Modulus of Silicon?", IEEE Journal of Microelectromechanical Systems, vol. 19, Issue 2, pp. 229-238, 2010 (DOI: 10.1109/JMEMS.2009.2039697).

Skip the explanation and go straight to the answer...

Crystal Orientation

Silicon is an anisotropic crystalline material with cubic symmetry (Fig. 1). The orientation of crystal structures is described by Miller indices, which are integer groups such as “hkl”. The integer groups can describe a direction vector, indicated by square brackets ("[hkl]"), or a plane perpendicular to the direction, indicated by parentheses ("(hkl)"). The cubic crystal directions are shown in Figure 2. In a cubic crystal, directions separated by 90° are equivalent, and these "families" of equivalent directions are indicated by corner brackets ("<hkl>").
Silicon crystal structure
Directions in a cubic crystal
Figure 1 - Silicon crystal structure Figure 2 - Cubic crystal directions

Silicon wafers are described by the crystal orientation of their surface and by the location of the wafer flats. The most common wafer type is called a "(100) wafer", meaning that the surface of the wafer is parallel to a (100) plane of the silicon in the wafer (Fig. 3). The major flat of these wafers is parallel to the <110> direction of the silicon crystal structure. Therefore, for a (100) silicon wafer, the "X" and "Y" directions most typically used in the layout of a MEMS design correspond the the <110> directions of the silicon crystal, as shown in Figure 4.
Crystal orientation of <100> silicon wafer X-Y axes on typical wafer
Figure 3 - Crystal orientation of "(100)" silicon wafer Figure 4 - Typical wafer orientation for MEMS design

Anisotropic Elasticity

Hooke's law describes the linear, isotropic relationship between stress (σ) and strain (ε) using compliance (S) and stiffness (C):
 σ=Cε  , or  ε=Sσ  .
The mechanical stiffness of materials under uniaxial loading is called the Young's modulus, and is typically represented by the symbol E in engineering texts, so Hooke's law is often written as σ=Eε  .
Hooke's law can be extended to anisotropic materials with a tensor formulation that includes orientation:
 σij=cijklεkl  , or   εij=sijklσkl  .
The general case of a fourth-rank tensor requires 81 terms. Fortunately, due to the properties of cubic symmetry, the anisotropic properties of a cubic crystal can be described with a 6x6 matrix using simplfied tensor notation and only 3 independent values are required in this matrix. The values in this matrix are determined with respect to a specific set of XYZ axes. The mathematical convention for these matrix transforms is to use the "principal crystallographic axes", which for a cubic crystal are the <100> directions. Therefore, unless otherwise specified, values for the elasticity constants are for the <100> axes, not the <110> axes typically used in MEMS design.
The best values for the elastic constants for silicon from the scientific literature are given in Table 1.

Table 1: Elastic Constants of Silicon
C: 109 Pa
S: 10-12 Pa
C11C12C44 S11S12S44
Silicon <100> 165.663.979.5 7.68-2.1412.6

Elasticity for Design Calculations

When doing design calculations, isotropic approximations are typically used, in order to make the calculation tractable. For an anisotropic material such as silicon, the correct value for the relevant crystal orientation must be used. For example, a silicon cantilever which is drawn parallel to the X axis in a MEMS layout tool for fabrication on a typical (100) silicon wafer would have its axis of bending in a <110> direction. Therefore, material properties for <110> should be used for design calculations, such as estimating the deflection of the cantilever. In this case, the appropriate value of E is 169 GPa. However, if the cantilever was fabricated "off-axis", aligned with a <100> direction, then the appropriate value of E for calculation is 130 GPa. The most common values of E used in design calculations are given in Table 2.

Table 2: Common Values of E for Silicon
Load CaseSymbolAppropriate E value
Axial stress/strain, incl. beam bending
  <110> directions ("X or Y axis")
  <100> directions ("45° off-axis")
E

 
 
169 GPa
130 GPa
Plate bending, ν = 0.064 for <110> in (100)E'110170 GPa
Thin film biaxial stress/strain, (100) planeB100180 GPa
Shear load, twisting a <110> ("X" or "Y") beamG50.9 GPa
Stress ConcentrationsEV166 GPa
Hydrostatic loads (Bulk modulus)B97.8 GPa
Temperature Coefficient of ETCE~ -60 ppm/°C

Elasticity for Finite Element Modeling

All modern finite element modeling (FEM) software has the capability to perform calculations using the anisotropic formulation of elasticity. The trick is to enter the appropriate values into the software. In some cases, the software will require the elastic constants in simplified tensor notation. In other cases, the user may be required to enter the elastic constants which have been transformed from the <100> axes to the to the axes of the silicon wafer crystal orientation. Refer to the documentation of your finite element software for details. The values for some common cases are given in Table 3.

Table 3: Common Elasticity Values for Silicon for FEM
Orthotropic stiffness values
for (100) silicon wafers
Ex = Ey = 169 GPa, Ez = 130 GPa
νxz = νyz = 0.36, νzx = νzy = 0.28, νxy = νyx = 0.064
Gyz = Gzx = 79.6 GPa, Gxy = 50.9 GPa
Cubic stiffness matrix values
(principal axes)
C11 = C22 = C33 = 165.6 GPa
C12 / 21 = C13 / 31 = C32 / 23 = 63.9 GPa
C11 = C22 = C33 = 79.6 GPa
Cubic stiffness tensor values
(principal axes)
c1111 = c2222 = c3333 = 165.6 GPa
c1122 = c1133 = 63.9 GPa
c2211 = c2233 = 63.9 GPa
c3311 = c3322 = 63.9 GPa
c1313 = c3113 = c1331 = c3131 = 79.6 GPa
c2323 = c3223 = c2332 = c3232 = 79.6 GPa
c1212 = c2112 = c1221 = c2121 = 79.6 GPa

Additional Notes

The elasticity values given here are correct for bulk silicon. Silicon structures with dimensions less than 100 nm will experience a significant reduction in E due to surface effects. For more information, see: Li, et al., "Ultrathin single-crystalline-silicon cantilever resonators...", Appl. Phys. Lett. 83, 3081 (2003).

Elasticity values for poly-crystalline silicon ("polysilicon") can vary significantly with deposition methods and conditions. For estimation purposes, polysilicon can be treated isotropically, using Epoly = 160 GPa, νpoly = 0.22 .

References

The material presented here is condensed from:

M. A. Hopcroft, W. D. Nix, and T. W. Kenny, "What is the Young's Modulus of Silicon?", Journal of Microelectromechanical Systems, vol. 19, Issue 2, pp. 229-238, 2010 (DOI: 10.1109/JMEMS.2009.2039697).

A complete list of references is available in that paper. Some selected references are:


M.A. Hopcroft, copyright 2013