Presented is an improved model of doubly-clamped microelectromechanical systems (MEMS) resonators implemented in very high-speed integrated circuit hardware description language for analogue and mixed-signal (VHDL-AMS). The model includes the effect of residual stress, which may severely shift the resonant frequency of the MEMS resonator from the analytical pre-designed value if the magnitude of intrinsic residual stress is imprecisely predicted. As the stress is not only determined by the fabrication process but also related to the structural dimensions of the resonators, in this work micro-Raman spectroscopy was utilised to measure and characterise the residual stress of a series of fabricated doubly-clamped polysilicon suspended gate field effect transistor type resonators with varying sizes. Combined with the experimentally determined residual stress, the proposed VHDL-AMS analytical model provides an error of < 3.5% resonant frequency shift with respect to the experimental result.1. Introduction: Microelectromechanical systems (MEMS) resonators are becoming increasingly attractive in timing and frequency control applications [1] because of their significant advantages of an intrinsic very high quality factor and high resonant frequency [2]. MEMS resonators also show exceptional possibilities for creating miniature-scale precision oscillators and filters because of their tiny geometric size and CMOS compatibility [3]. In addition, as the resonant frequencies are dependent on the operating conditions such as temperature and pressure, MEMS resonators show great potential as sensors such as the high sensitivity mass detector [4].Owing to the small size of the MEMS devices, microscale characterisation of residual stress is essential when considering the functionality and reliability of MEMS resonators. The residual stress can cause bending or buckling of the suspended structures of the MEMS resonators, which not only severely shifts the resonant frequency of the MEMS resonator from the analytical pre-designed value towards a desired application but also becomes the potential source leading to the failure of the MEMS resonators. If the residual stress can be accurately characterised, and eventually controlled, the design of the MEMS resonators will not be limited by residual stress. The designer may even take advantage of the intrinsic stress to increase the resonant frequency and quality factor of MEMS resonators.Residual stress is usually estimated by using the finite-element analysis (FEA) approach [5]. However, since the residual stress of a MEMS resonator is not only determined by the fabrication process but also related to the structural dimensions, there is a large margin of error in predicting stress using simulation, unless the experimental measurements are validated [6].There are several techniques for measuring and characterising residual stress in microscale structures, such as wafer curvature measurement [7], the X-ray microdiffraction technique [8] and the micro-Raman spectroscopy method [9][10][11][...