HighPerformance Computing Strategies for Complex Engineering Optimization Problems
View this Special IssueResearch Article  Open Access
Modeling the BoronDoping Silicon Beam by a Multilayer Model
Abstract
The borondoping silicon beam commonly used in microdevices exhibits a nonuniform material property along its thickness or width because of the gradient of boron concentration induced by diffusion process. The constant of rigidity, one of the most important parameters of microbeam, needs to be accurately calculated and designed in the development of high precise sensors and actuators. Current design methods, mainly depending on the analytical solutions derived under the assumption of a uniform material property or some commercial software for a varied property, are not adequate and time consuming to calculate the constant of rigidity of borondoping silicon beam. A multilayer model is proposed in this paper to replace the continuous solid model by dividing the beam into separated layers glued together. The finite element lamination method is utilized to acquire the equivalent Young modulus and moment of inertia of cross section of multilayer model. The equivalent values are calculated from doublelayer structures to multilayer ones based on the small deformation theory and the material mechanics theory. The proposed method provides an effective method to design the stiffness or frequency of microdevice and its results are validated by COMSOL simulation.
1. Introduction
Boron diffusion is the process commonly used in silicon based MEMS for acquiring a certain electrical property in microelectronics [1] or a process requirement in bulk silicon micromachining [2, 3]. It is one of the most important steps in bulk silicon dissolved wafer process (BSDWP) which is extensively used in fabrication of microaccelerometers [4, 5], microgyroscopes [6, 7], microswitches [8, 9], and microgears [10]. The microstructures fabricated by BSDWP, however, exhibit a varied property from top surface to bottom surface. When the movable or deformable parts are involved in devices, the borondoping silicon is always chosen as the structural material for the flexible components, like beams [11], flexure hinges [12, 13], and so forth. These components play the role of springs in microsensors to connect the movable parts to the substrate. Their stiffness or rigidity, therefore, has significant importance to the design of movable structures in MEMS, because it directly determines the sensitivity, stability, natural frequency, and measurement range of microdevices. Microbeam design always focuses on the dependence of the stiffness on its structure dimensions and material property. The analytical formula of stiffness of folded beam which is used for flexible supporting of mass of microaccelerometers has been given in details by Senturia [14]. The formula was derived by replacing the folded beam by a fixedfixed beam with the load on its middle position. Wong et al. have given out the theoretical analysis of stiffness constant of roundfolded beams in MEMS accelerometers [15]. The mentioned formulae are based on the assumption the beams have a uniform material property or Young’s modulus at least. However, the microstructures always exhibit a varied property after the silicon undergoes a series of doping and deposition processes, specifically in the selfstop etching process where a heavy borondoping step exists [3].
The doping process including boron ion implantation and diffusion will lead to a changed property of doping silicon from the surface to the interior. This paper proposes a finite element method to study the stiffness of borondoped microfolded beams of inplane capacitive microaccelerometers. This method utilizes an approximated multilayer model to replace the continuous solid beam. Each layer of the model has a uniform property and Young’s modulus of one layer is different from another, and then mechanics equations and boundary conditions are established analytically to analyze the stiffness of the approximated model under the transverse force vertical to stratification plane. In Section 2, the folded beam form in microaccelerometer is introduced, and its material distribution of borondoped silicon is recognized as an extended Gaussian function along the thickness. Young’s modulus corresponding to the concentration of boron is identified by the characteristics of selfstop etching of silicon. Section 3 gives out the equivalent model by dividing the cross section of beam into a series of rectangles, and its corresponding effective modulus, effective moment of inertia, and bending stiffness are derived by material mechanics equations. The followed simulation and experiments indicate a good agreement between proposed model and simulation. The application condition of this model and conclusion can been found in Section 4.
2. Material Property of Beam
The microfabrication, mainly consisting of surface silicon and bulk silicon technologies, always involves diffusion, deposition, ion implantation and etching, and so forth. Those processing steps generally concern the rearrangement or redistribution of materials, especially in the borondoping silicon of BSDWP [3]. The studied folded beams in this paper are fabricated by BSDWP and applied in microaccelerometers for supporting the movable mass block which can transfer the input acceleration into the shift of differential capacitance. The whole sensor of accelerometer is given in Figure 1(a), and one of its supporting beams is shown in Figure 1(b). The sensitive axis of sensor is parallel to axis, and the deformation plane of the folded beams is parallel to plane under the symmetrical layout of mass block and supporting beams. Therefore, the load along with is called the inplane load, while the load perpendicular to plane is called the outofplane load. The folded beam can be replaced by a fixedfixed one with the load applied on the center according to [14]. Figure 1(c) gives its cross section whose material property is represented by gradually changing colors. The bottom part of the beam has the highest concentration of 4 boron and the top part has the lightest according to [3, 16] in which the authors described the whole process flow of BSDWP in detail.
The boron concentration in silicon is assumed to be an extended half Gaussian function of the depth according to the solid diffusion mechanism [17]. The maximum concentration occurs at surface and then gradually decreases. Therefore, only two points of concentration and diffusion depth are needed to predict the boron profile in silicon solid. This condition can be achieved by two separate selfstopped etching processes using different etchants [18]. One etchant is (HF : HNO_{3} : CH_{3}COOH 1 : 3 : 8); the stop point is at a boron atom concentration of cm^{−3}. The other etchant is 30%–40% KOH; the stop point is at the concentration of cm^{−3}. Therefore, the Gaussian functions () of the boron profile with the boundary conditions are given as cm^{−3} and cm^{−3}, where is the coordinate along depth and and are the distances of stop points from the surface, respectively.
The relative concentration, defined as the ratio of boron atoms to silicon atoms per unit volume, can be expressed as with the boundary conditions and .
Compared with undoped silicon, heavily borondoped silicon has a 0.8% higher thermal expansion coefficient [19] and a 20%–30% higher Young’s modulus [20, 21]. Since it is related to the concentration of boron, Young’s modulus can be assumed as where (=) is the maximum relative concentration determined by the solid solubility of boron in silicon, (=200–220 GPa) is Young’s modulus to , and (=169 GPa) is Young’s modulus of silicon. The values of , therefore, can be calculated under the condition that Young’s modulus is proportional to the boron concentration.
3. Equivalent Model and Simulation
The constant of rigidity under the outofplane load is considered in this paper in order to evaluate the cross sensitivity of microaccelerometers. The cross section of beam is depicted in Figure 2(a); the density of red represents the concentration of boron in silicon, which decreases with the depth from the top surface according to the Fick’s diffusion laws. Figure 2(b) is a simplified model of borondoping silicon with four separated layers divided by three contour lines. When the beam is divided into a series of layers by the method applied to cross section, a multilayer beam model is gained to replace the continuous structure by such a finite element lamination (FEL) procedure.
According to the material mechanics theory [22], the doublelayer beam in Figure 3(a) with different Young’s modulus and same width can by replaced by the one in Figure 3(b) with the same Young’s module but a different width under the transferring condition that the total strains in both beams are equal to each other under the same load. It is concluded that “to replace the original material structure by another equivalent material structure the area must be multiplied by the modular ratio ”; therefore, the equivalent width in Figure 3(b) with the same thickness is
The equivalent Young modulus of doublelayer beam can be acquired by the same condition of material replacement. The relationship can be expressed as where is the equivalent Young modulus which equals
When a series of layers are considered, the doublelayer FEL procedure should be extended to a multilayers model of which the equivalent Young modulus is
Another important parameter involved in the evaluation of outofplane stiffness is the moment of inertia of cross section which can be calculated by the following expression: where is moment of inertia of cross section, is the moment of inertia of the th section relative to , is the moment of inertia of the th section relative to its own neutral axis, is the area of the th section, is the coordinate of neutral axis of the th section, and is the coordinate of neutral axis in the coordinate system defined in Figure 2; its position has been determined in [3]: where and are the width and thickness of the th layer, respectively.
Replacing Young’s modulus of the 2nd to th layer by that of the first layer, the corresponding width of the 2nd to th layer, , becomes . Thus the parameters mentioned above can be expressed as
Substituting (9)–(11) into (7), the moment of inertia of cross section of borondoping silicon beam can be expressed as
According to [14], the outofplane rigidity of folded beam can be written as
The proposed FEL method avoids the timeconsumed modeling and complicated meshing process and provides an analytical solution to the outofplane rigidity of borondoping silicon beams based on the material characteristics of diffusion mechanism. This solution procedure can be programmed by any computer language and the model parameters can be adjusted easily through parametric design for many kinds of beams. Figure 4 gives the calculation procedure of the FEL method. Because each layer after laminating is assumed to be with a uniform material property, the thickness of layers must be precisely controlled to acquire a more accurate model. This requirement is realized by defining the difference of Young’s modulus between two adjacent layers less than a set error, .
The results of FEL method are compared with those gained in commercial software COMSOL Multiphysics to verify its availability and effectiveness. The boron profile is determined by two diffusion depth points, and , which are gained by the selfstop etching process and their values are measured as 42 μm and 83 μm, respectively [3]. The parameters of simulated beams are shown in Table 1. Young’s modulus of monocrystalline silicon is 169 Gpa, while that of the heavy borondoping silicon is 200 Gpa. The dimensions of beam are equal to those of the designed accelerometers. The rigidity constants of uniform beams with 169 Gpa and 200 GPa are 1.372 kN/m and 1.624 kN/m, respectively. The constant becomes 1.481 kN/m when the modulus variation induced by boron implantation is considered. These constants acquired by software are shown in Table 2, in which the relative errors are less than 7%. Figure 5 shows the curve of relationship of rigidity constant and layer number from 100 to 1000, which indicates that a larger number, , results in a more precise constant comparing with that of COMSOL.


The FEL method shows an effective solution to the constants of rigidity of microbeams with a nonuniform material distribution. Table 3 lists the results from COMSOL and FEL method when the thickness of beams varies with process condition. The relative error decreasing with increasing thickness results from the ignored effects of hinges which, is not considered in simplification of microbeams. When the external force is applied on the folded beam, the single beam connecting to the anchors undergoes both bending and torsional deflection. The later will generate an additional displacement of the loaded point, of which the value is inversely proportional to the torsional stiffness. As the thickness of beam is increasing, the torsional stiffness increases linearly with thickness, while the bending stiffness increases cubically with thickness. Consequently, a thicker beam leads a higher torsional stiffness but a much higher bending stiffness, which means that the contribution of torsional displacement is declined and a more precise rigidity constant of folded beam is acquired with a smaller relative error.

4. Conclusion
To overcome the disadvantages of traditional methods, a finite element lamination model is developed to calculate the constant of rigidity of borondoping silicon beam which shows a nonuniform material property along its thickness. The method divides the continuous beam into a series of layers which are glued rigidly together, and each layer has the same properties which are determined by the characteristics of borondoping silicon material. The equivalent beams created by proposed method can be recognized as a multilayer structure. Its equivalent Young’s modulus and moment of inertia can be easily acquired by the theory of material mechanics. The constant of rigidity, therefore, can be derived by the small deformation theory. The results show a good agreement with those of commercial software COMSOL Multiphysics. The availability and effectiveness of this method can find its applications in MEMS fields, where the material nonuniformity and dimensional uncertainty are dominating, for a shortened modeling time and a higher calculating precise.
The presented method now is just suitable and powerful enough for the beams with a uniform structural dimension and a variable material property in any manner. Besides, the proposed method is not superior to FEM simulation in the time consumed. But it can be programmed in any computer language after reasonable simplification of borondoped beam which frequently appears in MEMS design. This convenience can save the engineers much time in learning complex software and sharply reduce the requirements of professional knowledge.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This research was partially supported by National Natural Science Foundation of China (Grant no. 51175437) and Fundamental Research Funds for the Central Universities (Grant no. ZYGX2011J085). The authors would like to acknowledge Shi Zhigui and Du Lianming for their assistance in fabrication and experiments.
References
 S. Solmi, F. Baruffaldi, and R. Canteri, “Diffusion of boron in silicon during postimplantation annealing,” Journal of Applied Physics, vol. 69, no. 4, pp. 2135–2142, 1991. View at: Publisher Site  Google Scholar
 Y. B. Gianchandani and K. Najafi, “A bulk silicon dissolved wafer process for microelectromechanical devices,” Journal of Microelectromechanical Systems, vol. 1, no. 2, pp. 77–85, 1992. View at: Publisher Site  Google Scholar
 W. Zhou, H. Yu, B. Peng, H. Shen, X. He, and W. Su, “Modeling the microstructure curvature of Borondoped silicon in bulk micromachined accelerometer,” Materials, vol. 6, no. 1, pp. 244–254, 2013. View at: Publisher Site  Google Scholar
 B. P. Van Drieenhuizen, N. I. Maluf, I. E. Opris, and G. T. A. Kovacs, “Forcebalanced accelerometer with mG resolution, fabricated using silicon fusion bonding and deep reactive ion etching,” in Proceedings of the International Conference on SolidState Sensors and Actuators (TRANSDUCERS '97), vol. 2, pp. 1229–1232, Chicago, Ill, USA, June 1997. View at: Google Scholar
 W. Zhou, B. L. Li, B. Peng, W. Su, and X. P. He, “Prediction of gap asymmetry in differential micro accelerometers,” Sensors, vol. 12, no. 6, pp. 6857–6868, 2012. View at: Publisher Site  Google Scholar
 H. Xie and G. K. Fedder, “A DRIE CMOSMEMS gyroscope,” in Proceedings of the 1st IEEE International Conference on SensorsIEEE Sensors 2002, vol. 2, pp. 1413–1418, June 2002. View at: Google Scholar
 H. Xie and G. K. Fedder, “Fabrication, characterization, and analysis of a DRIE CMOSMEMS gyroscope,” IEEE Sensors Journal, vol. 3, no. 5, pp. 622–631, 2003. View at: Publisher Site  Google Scholar
 M. Tang, A. Q. Liu, A. Agarwal, Q. X. Zhang, and P. Win, “A new approach of lateral RF MEMS switch,” Analog Integrated Circuits and Signal Processing, vol. 40, no. 2, pp. 165–173, 2004. View at: Publisher Site  Google Scholar
 J. Li, Q. X. Zhang, and A. Q. Liu, “Advanced fiber optical switches using deep RIE (DRIE) fabrication,” Sensors and Actuators A, vol. 102, no. 3, pp. 286–295, 2003. View at: Publisher Site  Google Scholar
 K. Fu, A. J. Knobloch, F. C. Martinez et al., “Design and fabrication of a siliconbased MEMS rotary engine,” in Proceedings of ASME International Mechanical Engineering Congress and Exposition (IMECE2001/MEMS '01), pp. 3447–3452, New York, NY, USA, November 2001. View at: Google Scholar
 Y. Sun, D. Fang, and A. K. Soh, “Thermoelastic damping in microbeam resonators,” International Journal of Solids and Structures, vol. 43, no. 10, pp. 3213–3229, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Y. Bellouard, A. A. Said, and P. Bado, “Integrating optics and micromechanics in a single substrate: a step toward monolithic integration in fused silica,” Optics Express, vol. 13, no. 17, pp. 6635–6644, 2005. View at: Publisher Site  Google Scholar
 G. Xu and L. Qu, “Some analytical problems of high performance flexure hinge and micromotion stage design,” in Proceedings of the IEEE International Conference on Industrial Technology (ICIT '96), pp. 771–775, Shanghai, China, December 1994. View at: Google Scholar
 S. D. Senturia, Microsystem Design, Kluwer Academic, New York, NY, USA, 2001.
 W. C. Wong, I. A. Azid, and B. Y. Majlis, “Theoretical analysis of stiffness constant and effective mass for a roundfolded beam in MEMS accelerometer,” Journal of Mechanical Engineering, vol. 57, no. 6, pp. 517–525, 2011. View at: Publisher Site  Google Scholar
 X. P. He, W. Su, B. Peng, and W. Zhou, “Electrostatic compensation method in frequency robustness design of micro accelerometer,” Journal of Micro/Nanolithography, Microfabrication, and Microsystems, vol. 11, no. 4, Article ID 043001, 2012. View at: Publisher Site  Google Scholar
 P. V. Zant, Microchip Fabrication: A Practical Guide to Semiconductor Processing (5th Ed.), McGrawHill, New York, NY, USA, 2004.
 H. L. Tang and H. Z. Ling, “Research on deep boron diffusion etchstop layer silicon wafer process,” Master ElectronicsElectronic Component, vol. 26, pp. 59–61, 2007. View at: Google Scholar
 B. S. Berry and W. C. Pritchet, “Internal stress and internal friction in thinlayer microelectronic materials,” Journal of Applied Physics, vol. 67, no. 8, pp. 3661–3668, 1990. View at: Publisher Site  Google Scholar
 K. Najafi and K. Suzuki, “Measurement of fracture stress, young's modulus, and intrinsic stress of heavily borondoped silicon microstructures,” Thin Solid Films, vol. 181, no. 12, pp. 251–258, 1989. View at: Publisher Site  Google Scholar
 J. T. Borenstein, P. Greiff, J. B. Sohn, and M. S. Weinberg, “Characterization of membrane curvature in micromachined silicon accelerometers and gyroscopes using optical interferometry,” in Proceedings of the SPIE Micromachining and Microfabrication Process Technology II, vol. 2879, pp. 116–125, October 1996. View at: Google Scholar
 E. J. Hearn, Mechanics of Materials I (Third Edition), ButterworthHeinemann, Oxford, UK, 2000.
Copyright
Copyright © 2014 Huijun Yu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.