Introduction
The value of the Young's modulus of silicon is often required for engineering designs using microelectromechanical 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 peerreviewed 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. 229238, 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>").


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.
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}=c_{ijkl}ε_{kl} , or ε_{ij}=s_{ijkl}σ_{kl} .
The general case of a fourthrank 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: 10^{9} Pa S: 10^{12} Pa 
C_{11}  C_{12}  C_{44}  S_{11}  S_{12}  S_{44} 
Silicon <100>  165.6  63.9  79.5  7.68  2.14  12.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 "offaxis", 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 Case  Symbol  Appropriate E value 
Axial stress/strain, incl. beam bending <110> directions ("X or Y axis") <100> directions ("45° offaxis")  E  169 GPa 130 GPa 
Plate bending, ν = 0.064 for <110> in (100)  E'_{110}  170 GPa 
Thin film biaxial stress/strain, (100) plane  B_{100}  180 GPa 
Shear load, twisting a <110> ("X" or "Y") beam  G  50.9 GPa 
Stress Concentrations  E_{V}  166 GPa 
Hydrostatic loads (Bulk modulus)  B  97.8 GPa 
Temperature Coefficient of E  TCE  ~ 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 
E_{x} = E_{y} = 169 GPa, E_{z} = 130 GPa ν_{xz} = ν_{yz} = 0.36, ν_{zx} = ν_{zy} = 0.28, ν_{xy} = ν_{yx} = 0.064 G_{yz} = G_{zx} = 79.6 GPa, G_{xy} = 50.9 GPa 
Cubic stiffness matrix values (principal axes) 
C_{11} = C_{22} = C_{33} = 165.6 GPa C_{12 / 21} = C_{13 / 31} = C_{32 / 23} = 63.9 GPa C_{11} = C_{22} = C_{33} = 79.6 GPa 
Cubic stiffness tensor values (principal axes) 
c_{1111} = c_{2222} = c_{3333} = 165.6 GPa c_{1122} = c_{1133} = 63.9 GPa c_{2211} = c_{2233} = 63.9 GPa c_{3311} = c_{3322} = 63.9 GPa c_{1313} = c_{3113} = c_{1331} = c_{3131} = 79.6 GPa c_{2323} = c_{3223} = c_{2332} = c_{3232} = 79.6 GPa c_{1212} = c_{2112} = c_{1221} = c_{2121} = 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 singlecrystallinesilicon cantilever resonators...", Appl. Phys. Lett. 83, 3081 (2003).Elasticity values for polycrystalline silicon ("polysilicon") can vary significantly with deposition methods and conditions. For estimation purposes, polysilicon can be treated isotropically, using E_{poly} = 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. 229238, 2010 (DOI: 10.1109/JMEMS.2009.2039697).
A complete list of references is available in that paper. Some selected references are:
 "Foundations of MEMS",
2006
C. Liu
Book published by Pearson Education, Inc., Upper Saddle River, NJ
ISBN: 0131472860
 "Microsystem Design",
2001
S. D. Senturia
Book published by Kluwer Academic Publishers, Boston
ISBN: 0792372468
 "Physical Properties of Crystals: Their Representation by Tensors and Matrices",
1985
J. F. Nye
Book published by Oxford University Press, Oxford
ISBN: 0198511655
 "Electronic effects in the constants of ntype silicon", 1967
J. J. Hall
Physical Review, 161, (3), pp. 756
 "Youngs' Modulus, Shear Modulus and Poisson's Ratio in Silicon and Germanium", 1965
J. J. Wortman and R. A. Evans
Journal of Applied Physics, 36, (1), pp. 153156