Multiphysics Optical-Thermal and Mechanical Modeling of a CMOS-SOI-MEMS Infrared Sensor with Metasurface Absorber

1. Introduction

Thermal infrared (IR) sensors are indispensable for a wide range of applications, including night vision, thermal imaging, environmental monitoring, gas sensing, and industrial process control [1,2,3,4]. The demand for compact, high-performance, and cost-effective uncooled IR detectors has spurred significant advancements in Micro-Electro-Mechanical Systems (MEMS) fabricated on CMOS-compatible Silicon-On-Insulator (SOI) platforms. These CMOS-SOI-MEMS platforms offer exceptional thermal isolation, reduced parasitic effects, and seamless integration with on-chip electronics, enabling scalable and low-power sensing solutions [5,6,7]. By leveraging standard CMOS fabrication processes, these platforms pave the way for affordable, high-volume production of IR sensors suitable for both consumer and industrial applications.

The Thermal CMOS-SOI-MEMS (TMOS) sensor represents a breakthrough in uncooled IR detection, utilizing a suspended MOSFET as the sensing element to transduce IR-induced temperature changes into highly sensitive electrical signals [8,9,10,11]. The TMOS architecture achieves superior responsivity by maximizing the temperature rise (ΔT) on a thermally isolated micro-platform, supported by thin holding arms that minimize conductive heat losses. Key to its performance is efficient absorption of incident IR radiation across the 2–14 µm wavelength range, which is critical for applications such as CO₂ gas sensing (4.26 µm) and thermal imaging (8–14 µm). To address this, advanced nanophotonic components, such as Metasurface absorbers based on metal–insulator–metal (MIM) structures, have emerged as a powerful solution [12,13,14,15]. These Metasurfaces achieve near-unity absorption at targeted wavelengths while remaining fully compatible with CMOS fabrication, offering spectral selectivity for application-specific sensing.

In our previous work [16], we investigated how incorporating MIM absorbers affects the optical and thermal performance of CMOS-SOI-MEMS thermal IR sensors, and we detailed the simulation methodology as well as the influence of MIM geometry on the absorption spectrum.

However, integrating Metasurface absorbers into CMOS-SOI-MEMS sensors introduces complex trade-offs across optical, thermal, and mechanical domains. Enhanced absorption increases thermal capacitance, potentially slowing the sensor’s temporal response, while added mass may alter the mechanical stability of the suspended microstructure. Existing studies often focus on isolated aspects, such as electromagnetic design of absorbers or simplified thermal modeling, failing to capture the coupled interactions that govern overall sensor performance [17,18]. This lack of a holistic multiphysics framework limits the ability to optimize responsivity, speed, and structural integrity simultaneously, hindering the translation of advanced optical concepts into manufacturable devices.

This work addresses this gap by developing a comprehensive multiphysics modeling framework for a TMOS sensor enhanced with Metasurface absorbers. The optical response was modeled using Lumerical FDTD Solutions [19], where finite-difference time-domain (FDTD) simulations quantified the absorption efficiency and incident radiation distribution, serving as input heat sources for the subsequent thermal finite element analysis (FEA). The thermal and mechanical simulations were conducted using ANSYS Mechanical [20], with the thermal response validated against an analytical lumped-element RC circuit model to evaluate the sensor’s temperature rise and dynamic behavior, revealing the influence of Metasurface mass on the thermal time constant. Furthermore, mechanical modal and harmonic response analyses assessed the impact of Metasurface integration on natural frequencies and vibration modes, ensuring structural robustness under environmental excitations. As illustrated in Figure 1, this integrated optical–thermal–mechanical modeling approach provides a predictive design framework that enables systematic exploration of trade-offs toward high-responsivity, fast, and mechanically stable CMOS-compatible IR sensors for advanced applications.

2. Sensor Design and Architecture

2.1. TMOS Sensor Architecture

The CMOS-SOI-MEMS infrared sensor employs a TMOS pixel architecture, featuring a nano-machined, suspended transistor fabricated in a standard 0.13-µm CMOS-SOI process for monolithic integration with readout electronics [10]. The core sensing element is a thermally isolated metal-oxide-semiconductor field-effect transistor (MOSFET) operating in the subthreshold regime, which converts absorbed IR radiation into temperature-dependent electrical signals with high sensitivity and low power consumption.

The pixel comprises a suspended thermal platform (stage) anchored by thin holding arms, which ensure mechanical stability while minimizing thermal conductance to the substrate. Fabrication involves dry etching techniques, including reactive ion etching (RIE) for front-side dielectrics and deep reactive ion etching (DRIE) for back-side silicon removal, with the buried oxide layer serving as an etch stop. To enhance thermal isolation, the arms incorporate active silicon and polysilicon layers rather than metal interconnects, reducing conductive heat flow. Additional isolation is achieved through wafer-level vacuum packaging (WLP) via glass-frit bonding, incorporating a getter layer to sustain a high vacuum (~1 Pa) and suppress gas conduction and radiative losses [21]. Figure 2 illustrates the TMOS architecture.

2.2. Metamaterial Integration

To enhance the optical performance of the TMOS sensor, Metasurface absorbers based on metal–insulator–metal (MIM) structures are integrated atop the suspended stage. Metasurfaces, engineered nanostructures that control electromagnetic wave propagation at subwavelength scales, enable tailored absorption properties unattainable with conventional materials [17,18,22,23]. While this study evaluates two narrowband MIM configurations as exemplars, the proposed multiphysics modeling framework is versatile, applicable to any absorber architecture - whether narrowband, broadband, multilayer, or with varied structural geometries - ensuring broad relevance to diverse IR sensing applications.

The two representative MIM configurations demonstrate specific functionalities:

• Configuration 1: Optimized for peak absorption at 4.26 µm, targeting the CO₂ absorption band for non-dispersive infrared (NDIR) gas sensing.
• Configuration 2: Designed for narrowband absorption near 10 µm, suitable for biomedical sensing (e.g., protein vibrational fingerprints) and thermal imaging.

Each MIM unit cell comprises periodic layered stacks, compatible with CMOS post-processing, arranged to form a Metasurface covering the pixel’s active area. The configurations serve to illustrate the framework’s ability to quantify the optical efficiency, thermal dynamics, and mechanical stability of any absorber design. This approach evaluates absorption enhancement, thermal capacitance effects, and mass-induced frequency shifts, as detailed in subsequent sections. Figure 3 depicts the MIM unit cell configurations, including materials, thicknesses, and dimensions (in µm).

3. Optical Analysis

The optical performance of the CMOS-SOI-MEMS TMOS sensor hinges on its ability to efficiently absorb incident infrared (IR) radiation, which directly drives the thermal response. This section presents a comprehensive optical modeling framework using (FDTD) simulations to quantify absorption efficiency across the 2–14 µm band. The framework calculates incident irradiance, electromagnetic power dissipation, and absorption characteristics for baseline CMOS, titanium nitride (TiN)-enhanced, and Metasurface-enhanced configurations. While two metal–insulator–metal (MIM) Metasurface designs are evaluated as exemplars, the methodology is versatile, applicable to any absorber architecture—narrowband, broadband, multilayer, or varied geometries. Results provide heat sources for thermal analysis (Section 4).

3.1. Incident Irradiance and Radiometric Quantities

To link the electromagnetic simulations with realistic sensor operation, it is necessary to relate blackbody radiation to the optical power absorbed in the detector. The fundamental radiometric quantity is the spectral radiance $L_{\lambda }\left(T\right)$, which specifies the emitted power per unit area, per unit solid angle, and per unit wavelength of a blackbody at temperature $T$ [24,25]:

$$L_{\lambda }\left(T\right)={\displaystyle \frac{2h{c}^2}{{\lambda }^5\exp \left(\frac{hc}{\lambda k_{B}T}\right)-1}} \left[W\cdot m^{-2}\cdot s{r}^{-1}\cdot \mu m^{-1}\right]$$

(1)

where $h$ is Planck’s constant, $c$ is the speed of light, and $k$ is Boltzmann’s constant.

The spectral irradiance, $E_{\lambda }$, at the sensor plane is obtained by integrating radiance over the collection cone of the optical system:

$$E_{\lambda }\left(T\right)={\int }_{\mathsf{\Omega }}{L}_{\lambda }\left(T\right) cos\theta d\mathsf{\Omega }\left[W\cdot m^{-2}\cdot \mu m^{-1}\right]$$

(2)

where $\mathsf{\Omega }$ is the solid angle defined by the numerical aperture of the optics and $T$ is the object temperature. The factor $cos\theta$ in the irradiance equation accounts for the angle of incidence relative to the sensor’s normal axis, where $\theta$, is the polar angle, shown in Figure 4.

The total absorbed spectral power in a detector pixel of effective area $A_{\mathrm{eff}}$ is then:

$$P_{\mathrm{abs}}\left(\lambda \right)=A_{\mathrm{eff}}{\eta }_{\mathrm{abs}}\left(\lambda \right)E_{\lambda }\left(T\right) \left[W\cdot \mu m^{-1}\right]$$

(3)

where ${\eta }_{\mathrm{abs}}\left(\lambda \right)$ is the absorption efficiency obtained from the FDTD simulations.

In practice, the FDTD results were normalized to unit incident irradiance (1 W·m⁻²). Using the above relations, the normalized absorption spectra (optical efficiencies) were scaled to the blackbody-derived irradiance corresponding to a given source temperature and optical collection geometry. The resulting absorbed power distribution was used as the volumetric heat generation input for the thermal FEA (Section 4).

3.2. Power Absorbed in a Dispersive and Absorptive Material

To couple the optical and thermal domains, the volumetric heat generation in the Metasurface absorber is computed from the local electromagnetic dissipation. Starting from Poynting’s theorem, the power dissipated per unit volume in a dispersive and absorptive medium is expressed as the time-averaged product of the electric and magnetic fields with the imaginary part of the material’s permittivity and permeability [26,27,28]:

$$⟨q⟩={\displaystyle \frac{1}{2}}\omega {\epsilon }_0\mathrm{Im}\left\{\epsilon \left(\omega \right)\right\}\mid \mathbf{E}{\mid }^2+{\displaystyle \frac{1}{2}}\omega {\mu }_0\mathrm{Im}\left\{\mu \left(\omega \right)\right\}\mid \mathbf{H}{\mid }^2 \left[W\cdot m^{-3}\right]$$

(4)

where $\omega$ is the angular frequency of the incident light, $\mathbf{E}$ and $\mathbf{H}$ are the local electric and magnetic fields, and $\epsilon \left(\omega \right)$and $\mu \left(\omega \right)$are the complex, frequency-dependent permittivity and permeability of the material. For non-magnetic materials ($\mu \left(\omega \right)\approx 1$), the magnetic contribution is negligible, yielding $⟨q⟩={\displaystyle \frac{1}{2}}\omega {\epsilon }_0\mathrm{Im}\left\{\epsilon \right(\omega \left)\right\}\mid \mathbf{E}{\mid }^2$.

In this work, the three-dimensional spatial distribution of $⟨q\left(\mathbf{r}\right)⟩$ is computed using full-wave FDTD simulations of the Metasurface absorber. The resulting volumetric heat generation profile is exported to the thermal finite element model and applied as a distributed source term for subsequent thermal analysis (Section 4). This approach ensures a physically consistent mapping from electromagnetic absorption to the thermal response of the CMOS-SOI-MEMS pixel.

3.3. Optical Efficiency Analysis

The optical performance of the MEMS thermal IR sensor is dictated by the efficiency with which the suspended CMOS membrane absorbs incident infrared radiation. Since the absorbed power is directly converted into heat and subsequently sensed as a temperature rise, maximizing absorption is a critical design goal. In this study, two absorption-enhancement strategies are investigated. First, we examine the intrinsic absorption of the CMOS device layer, both in its unmodified form and with an embedded titanium nitride (TiN) layer that improves optical-to-thermal conversion. Second, we evaluate advanced metamaterial absorber designs based on metal–insulator–metal (MIM) stacks, which further tailor the spectral response and boost efficiency. The following subsections present the modeling framework and simulation results for these approaches.

3.3.1. CMOS layers Absorption and TiN-Embedded Layer

To accurately model the intrinsic absorption of the CMOS device layer, it is necessary to account for the layered structure of the membrane and the optical properties of its constituent materials. Figure 5 presents the input data based on [19,29,30] that used in the FDTD simulations and the device material layers.

The cross-sectional schematic of the CMOS–SOI stack illustrates a multilayer configuration comprising alternating dielectric and conductive films. From top to bottom, the structure includes silicon dioxide (SiO₂), an embedded titanium nitride (TiN) layer, an additional SiO₂ layer, a thin silicon nitride (Si₃N₄) film, and another SiO₂ layer. These layers are followed by the polysilicon gate and the active silicon device layer, supported on a thick buried SiO₂ layer.

The complex permittivity spectra highlight the distinct optical characteristics of each constituent material. Silicon remains nearly lossless across the 2–14 μm spectral range, ensuring transparency in the mid-infrared region. SiO₂ and Si₃N₄ exhibit pronounced phonon-polariton resonances beyond 8–9 μm, resulting in wavelength-selective absorption peaks. In contrast, TiN demonstrates metallic behavior, with a negative real permittivity and an increasing imaginary component, enabling efficient broadband light absorption.

Beyond its intrinsic optical loss, the TiN layer can also contribute to enhanced absorption by acting as an optical cavity within the dielectric stack. Its strategic placement and thickness strongly influence interference and field confinement effects, thereby modifying the local absorption distribution. The location and thickness of the TiN layer were systematically studied to optimize its contribution to the overall optical performance z [31].

These combined material and structural properties govern how the CMOS–SOI stack interacts with incident infrared radiation. The embedded TiN layer enhances absorption efficiency within the mid-IR band by introducing tunable loss and resonant field confinement, while remaining fully compatible with standard CMOS fabrication processes.

FDTD simulations were conducted using Lumerical FDTD [19] to quantify the absorption efficiency of the CMOS–SOI sensor stack over the 2–14 µm range. Two configurations were analyzed: the baseline CMOS stack without TiN and the same structure incorporating a 100 nm embedded TiN layer. The absorption efficiency, defined as the ratio of absorbed to incident optical power, was used to directly compare both designs, as shown in Figure 6.

3.3.2. Metamaterial Absorber Layers

-- While the embedded TiN layer enhances broadband absorption, further improvement in both efficiency and spectral selectivity can be achieved through the implementation of engineered metamaterial absorbers. These structures, based on metal–insulator–metal (MIM) stacks, are designed to support wavelength-selective resonant cavity modes that confine and trap infrared radiation within the sensing platform. By tailoring the geometric and material parameters, MIM absorbers can be engineered to operate as narrowband band-pass filters, selectively enhancing absorption at specific wavelengths of interest.

In this study, two representative MIM configurations - denoted MIM 1 and MIM 2—were designed and analyzed, each differing in layer thickness and overall stack geometry. The bottom metallic layer is optically opaque, ensuring zero transmission and allowing the simulation of a single unit cell. To emulate an extended Metasurface array integrated on top of the CMOS–SOI–MEMS sensor pixel, periodic boundary conditions were applied along the lateral directions, while perfectly matched layer (PML) boundaries were used along the top and bottom to eliminate spurious reflections.

Simulation results demonstrate that both MIM designs substantially enhance optical absorption within the targeted mid- and long-wavelength infrared regions, exhibiting near-unity efficiency at their respective resonance wavelengths. These resonant peaks, corresponding to the designed cavity modes, confirm the tunable narrowband filtering behavior of the MIM structures. The resulting absorption spectra for the baseline CMOS sensor, the CMOS+TiN stack, and the two MIM-enhanced configurations are presented in Figure 7.

From Figure 7, it is evident that incorporating the MIM1 structure significantly enhances the optical efficiency at ~10 µm, achieving nearly perfect absorption with an efficiency of 0.99. In contrast, without the MIM1 layers, the CMOS stack yields only 0.7564 when including the TiN layer and as low as 0.032 in the absence of TiN. Likewise, the addition of MIM2 at ~10 µm improves the efficiency to 0.9804, compared to 0.8353 for the CMOS–TiN configuration without MIM2 and 0.5740 when TiN is not present.

Furthermore, electric field and absorbed power density distributions were computed to illustrate how light is confined and absorbed within the MIM layers (see Figure 8), providing insight into the mechanisms of enhanced absorption. These results demonstrate the potential of MIM configuration to substantially improve the photon-to-heat conversion efficiency of CMOS-SOI-MEMS infrared sensors.

Using the simulated optical absorption efficiency as a function of wavelength, the spectral and total irradiance incident on the sensor from an object at temperature $T$ can be calculated. To illustrate the effect of object temperature on the irradiance received by the sensor, Figure 9 presents both the spectral and integrated irradiance for two representative temperature ranges: 300-310 K and 300-500 K. These plots provide insight into the amount of optical power delivered to the sensor under different thermal conditions, which is critical for evaluating the sensor’s response.

Using the total irradiance, the total absorbed power (Equation 3) can be computed and subsequently applied as a heat source in the thermal simulations (Section 4). This approach enables the evaluation of the thermal impact of the metamaterial layers on the performance of the CMOS-SOI-MEMS sensor pixel.

4. Thermal Analysis

The thermal performance of the TMOS sensor determines its responsivity to incident IR radiation. This section presents a thermal modeling framework that quantifies temperature rise (ΔT) and dynamic response, integrating optical heat sources from FDTD simulations (Section 3). The framework employs FEA and an analytical lumped-element RC circuit model to evaluate the impact of Metasurface absorbers on thermal behavior. Applicable to any absorber architecture—narrowband, broadband, or multilayer—the framework ensures versatility for diverse IR applications.

4.1. Governing Principles of Heat Transfer

The thermal behavior of MEMS IR sensors is dominated by conduction through the suspended support beams and solid layers, with radiation as a secondary mechanism. Convection is negligible in vacuum. The transient heat diffusion equation is [32,33]:

$$\rho c_p{\displaystyle \frac{\partial T}{\partial t}}=k{\nabla }^{2}T+Q_{gen}$$

(5.1)

$${\displaystyle \frac{\partial T}{\partial t}}=\alpha {\nabla }^{2}T+{\displaystyle \frac{Q_{gen}}{\rho c_p}}$$

(5.2)

where $\rho$ is density, $c_p \left[J\cdot k{g}^{-1}\cdot K^{-1}\right]$ is specific heat capacity, $k \left[W\cdot m^{-2}\cdot K^{-1}\right]$ is thermal conductivity, $T$ is temperature, $\alpha$ is the thermal diffusivity and defined as $\alpha \equiv {\displaystyle \frac{k}{\rho c_p}} \left[m^2\cdot s^{-1}\right]$ and $Q_{gen} \left[W\cdot m^{-3}\right]$ is the volumetric heat source.

4.2. Lumped-Circuit Equivalent Model

For microscale devices, when the thermal diffusion length $L_d=\sqrt{\alpha {\tau }_{\mathrm{t}\mathrm{h}}}$ exceeds the platform size, the stage temperature can be assumed uniform. The sensor can then be modeled as a lumped system:

$$L_d=\sqrt{\alpha {\tau }_{th}} \left[m\right]$$

(6)

If $L_d$ exceeds the dimensions of the micro-platform, temperature gradients within the stage are negligible, allowing the entire platform to be modeled as a lumped thermal capacitance $C_{\mathrm{th}}$ connected to the substrate through an overall thermal conductance $G_{\mathrm{th}}$ [16,34]:

$$C_{\mathrm{th}}{\displaystyle \frac{d\mathsf{\Delta }T\left(t\right)}{dt}}+G_{\mathrm{th}}\mathsf{\Delta }T\left(t\right)=P_{\mathrm{opt}}\left(t\right)$$

(7)

where $\mathsf{\Delta }T$ is the temperature rise of the sensing stage and $P_{\mathrm{opt}} \left[W\right]$ is the absorbed optical power. The thermal conductance of a single support beam is obtained from Fourier’s law:

$$G_{\mathrm{th}}={\displaystyle \frac{kA}{L}}$$

(8)

where $A$ is the beam cross-sectional area, and $L$ its length. For multilayer or heterogeneous beams structures, equivalent thermal conductivities can be computed for series or parallel arrangements, similar to electric conductivity [16,35,36]:

$$k_{\mathrm{eq},\mathrm{serial}}=L{\left(\sum _{i=1}^n{\displaystyle \frac{l_i}{k_i}}\right)}^{-1},k_{\mathrm{eq},\mathrm{parallel}}={\displaystyle \frac{1}{A}}\sum _{i=1}^n{k}_i{A}_i$$

(9)

with $l_i,k_i,A_i$the length, conductivity, and cross-section of each material.

The total thermal capacitance is obtained by summing contributions from all constituent materials:

$$C_{\mathrm{th}}=\sum _{i=1}^n{\rho }_i{c}_i{V}_i$$

(10)

where $V_i$is the volume of the $i$-th layer. Equivalently the thermal time constant is defined by ${\tau }_{th}=C_{th}/G_{th}$.

The thermal and physical properties of the materials relevant to the CMOS–SOI stack are summarized in Table 1.

The sensor pixel is modeled under the lumped thermal assumption, where the suspended stage and MOSFET layers (with or without MIM additions) constitute the dominant thermal capacitance, and the supporting arms act as the primary thermal resistors. Heat flows independently through each arm, forming parallel thermal paths that reduce the overall thermal resistance. Figure 10 depicts both the 3D pixel structure and its equivalent thermal circuit.

4.3. Simulation Setup and Results

The 3D FEA model was constructed using the material properties listed in Table 1. The frame-connected edges of the supporting arms were treated as isothermal boundaries, representing the bulk silicon heat sink. Heat input was applied to the sensor stage, corresponding to the absorbed optical power calculated from the optical efficiency and irradiance analysis (Section 3), defined as the difference between incident irradiance from the target and the sensor’s radiative emission at its initial temperature:

$$\begin{gathered} P_{\mathrm{opt}}={\int }_{{\lambda }_{\min }}^{{\lambda }_{\max }}{P}_{\mathrm{abs}} \left(\lambda \right)d\lambda -{\int }_{{\lambda }_{\min }}^{{\lambda }_{\max }}{P}_{emitted}\left(\lambda \right)d\lambda \\ ={\int }_{{\lambda }_{\min }}^{{\lambda }_{\max }}{A}_{\mathrm{eff}} {\eta }_{\mathrm{abs}}\left(\lambda \right)E_{\lambda }\left(T_{target}\right)d\lambda -{\int }_{{\lambda }_{\min }}^{{\lambda }_{\max }}{A}_{\mathrm{eff}} {\eta }_{\mathrm{abs}}\left(\lambda \right)E_{\lambda }\left(T_0\right)d\lambda \end{gathered}$$

(11)

where $T_0$is the initial sensor temperature and $T_{\mathrm{target}}$ is the temperature of the observed object that emits the incident radiation.

Given the lumped thermal assumption and the large thermal diffusion length relative to the stage dimensions, the detailed spatial distribution of absorbed power is not critical. Although FDTD simulations provide fine-scale absorption profiles, their mesh cannot be directly matched to the coarser FEA thermal mesh. Therefore, the total absorbed optical power $P_{\mathrm{o}\mathrm{p}\mathrm{t}}$is applied uniformly across the sensor stage, yielding comparable thermal results with significantly reduced computational complexity. Figure 11 illustrates the steady-state temperature distribution for a target object temperature of 301 K, assuming an initial sensor temperature of 300 K, for the CMOS-SOI-MEMS sensor with the TiN layer and without MIM structures, validating the lumped model assumption.

As shown in Figure 11, the stage and active MOSFET layer exhibit a nearly uniform temperature, while significant gradients occur along the support arms, validating the lumped thermal model. To assess the dynamic thermal behavior, the transient stage temperature $\mathsf{\Delta }T\left(t\right)$was simulated for four configurations: without TiN, with TiN, with TiN and MIM configuration 1, and with TiN and MIM configuration 2. Figure 12 presents the transient results of these simulations, alongside predictions from the analytical lumped-element model, enabling direct comparison with the full 3D FEA results. Thermal capacitance and conductance were computed from Equations (7)–(11) using the material parameters in Table 1, and the FEA simulations were fitted to both the analytical solution and the corresponding $G_{\mathrm{t}\mathrm{h}}$and $C_{\mathrm{t}\mathrm{h}}$values.

Figure 12 demonstrates excellent agreement between the full 3D FEA results and the analytical thermal circuit model. For the sensor without TiN, the analytical time constant is 78.65 ms, closely matching the FEA value of 81 ms. With TiN, the values are 78.93 ms (analytical) and 80.80 ms (FEA). Incorporation of MIM layers increases thermal capacitance and the corresponding time constant: MIM configuration 1 yields 94.83 ms (analytical) versus 96.8 ms (FEA), and MIM configuration 2 yields 127.50 ms (analytical) versus 129.4 ms (FEA), confirming the model’s predictive accuracy.

5. Mechanical Analysis

The dynamic mechanical integrity of MEMS is a fundamental determinant of sensor performance, stability, and reliability [37]. Specifically, the natural frequencies and mode shapes dictate the susceptibility response to external vibrations, which can cause excessive displacement, structural fatigue, or noise-inducing signal distortion. This analysis verifies that the operational bandwidth is safely decoupled from the device’s mechanical resonances.

5.1. Modal and Harmonic Analysis

The dynamic behavior of the suspended micro-platform is often conceptually modeled as a mass-spring system to provide intuitive insight, with the equivalent system illustrated in Figure 13. For the undamped, single-degree-of-freedom approximation, the governing equation is [36,38]:

$$m\ddot{x}\left(t\right)+k_{m}x\left(t\right)=0$$

(12)

With a harmonic solution $x\left(t\right)=X{e}^{j2\pi f_{n}t}$ yielding the natural frequency:

$$f_n={\displaystyle \frac{1}{2\pi }}\sqrt{{\displaystyle \frac{k_m}{m}}}$$

(13)

For the complex, multi-material 3D MEMS microstructure a FEM eigenvalue problem is necessary to capture all vibration modes:

$$\left[K\right]\left\{\varphi \right\}={\omega }_n^2\left[M\right]\left\{\varphi \right\}$$

(14)

where $\left[K\right]$ and $\left[M\right]$ are the global stiffness and mass matrices, ${\omega }_n$ are the natural angular frequencies, and $\left\{\varphi \right\}$ are the corresponding mode shapes. Each mode represents a distinct three-dimensional deformation of the suspended platform.

To assemble the stiffness $\left[K\right]$ and mass $\left[M\right]$ matrices, the FEM solver requires the material properties of each structural component: density $\rho$, Young’s modulus $E$, and Poisson’s ratio $\nu$. These parameters define the inertial and elastic response of the microstructure. The mechanical properties are depicted in Table 2.

The FEA model was constructed with a dense mesh for discretization, and fixed support boundary conditions were applied at the anchor points (Figure 14), consistent with the methodology used for similar high-performance MEMS devices. This simulates the physical constraint where the suspending arms connect to the rigid substrate frame.

The harmonic response analysis quantifies the steady-state displacement when the MEMS device is subjected to a sinusoidal external force $F\left(t\right)=F_0{e}^{j{\omega }_{m}t}$, the general equation of motion for a single-degree-of-freedom system is:

$$m\ddot{x}\left(t\right)+c\dot{x}\left(t\right)+k_{m}x\left(t\right)=F_0{e}^{j{\omega }_{m}t}$$

(15)

With the steady-state displacement amplitude:

$$\mid X\left(\omega \right)\mid ={\displaystyle \frac{F_0}{\sqrt{(k_m-m{\omega }_m^2{)}^2+(c_m\omega {)}^2}}}$$

(16)

where $c_m$ accounts for damping due to air, anchors, and internal dissipation. This frequency response function highlights resonance peaks near the natural frequencies. In practice, the MEMS sensor has many degrees of freedom, and its full dynamics are captured by the finite element formulation:

$$\left[M\right]\left\{\ddot{u}\right(t\left)\right\}+\left[C\right]\left\{\dot{u}\right(t\left)\right\}+\left[K\right]\left\{u\right(t\left)\right\}=\left\{F\right(t\left)\right\}$$

(17)

where $\left\{u\right(t\left)\right\}$ contains all nodal displacements. This multi-DOF analysis provides the accurate three-dimensional vibration characteristics, while the lumped single-DOF model offers intuitive insight into resonance behavior.

5.2. Mechanical Simulation Results and Physical Interpretation

The vibrational characteristics were investigated across three critical configurations: the baseline structure (Without MIM), and two configurations integrating the MIM absorber layer (With MIM 1 and With MIM 2). The addition of the absorber introduces mass loading that directly modifies the effective stiffness-mass balance, $\sqrt{k_m/m}$.

Table 3 summarizes the first six natural frequencies. A clear downward shift in the lower-order frequencies is observed with the inclusion and thickening of the MIM absorber layer.

Figure 15 illustrates the corresponding mode shapes, revealing a critical distinction in the dynamic response:

• Low-Order Modes ($f_1$ to $f_4$): Stage-Dominated Motion. These modes exhibit out-of-plane translation, torsional rotation, and combined deformation of the suspended platform itself (Figures 15a-d). The substantial frequency reduction (e.g., dropping from 31.03 kHz to 24.88 kHz with MIM 2) is a direct consequence of the MIM’s mass being concentrated on the suspended stage.

• High-Order Modes ($f_5$ and $f_6$): Arm-Dominated Motion. These modes, occurring at significantly higher frequencies (), primarily involve the bending and twisting of the support arms (Figure 15e,f). Crucially, their frequencies are negligibly affected by the MIM layers. This confirms that the MIM integration selectively affects the stage’s inertial properties while preserving the high stiffness and mechanical integrity of the support structure. The high arm frequencies (the fundamental frequency) are essential for ensuring thermal isolation compliance without introducing low-frequency structural vulnerabilities.

The harmonic response analysis was then performed to evaluate the displacement amplitude of the suspended structure under sinusoidal excitation. Figure 16 shows the frequency response for the three configurations. The resonance peaks correspond closely to the modal frequencies listed in Table 2, confirming the consistency of the two approaches.

The simulation results confirm that the sensor’s operational bandwidth lies well below its fundamental resonance frequency, ensuring that normal operation remains unaffected by mechanical resonance effects. This verifies that the suspended stage maintains dynamic stability under typical excitation conditions.

Critically, the lowest natural frequency of the sensor with the MIM absorber is approximately 28 kHz, which is substantially higher than the frequency range of common environmental disturbances encountered in both indoor and outdoor settings (typically below 1 kHz). This indicates strong mechanical robustness, with minimal susceptibility to destructive vibrations or excessive noise coupling during operation.

Furthermore, while the metamaterial MIM layers slightly reduce the resonance frequencies due to their added mass, they do not significantly compromise the overall mechanical stability of the device. The modest frequency shift demonstrates that the MIM layers can be integrated without adversely affecting the structural integrity of the suspended platform. Thus, the metamaterial absorber enhances optical absorption performance while preserving the sensor’s mechanical characteristics.

By maintaining sufficiently high natural frequencies, the design ensures that environmental vibrations and mechanical impacts do not interfere with thermal signal detection. This provides confidence that the MEMS IR sensor can operate reliably in real-world deployment scenarios without mechanical resonance degrading sensitivity or stability.

6. Discussion

This study develops a comprehensive multiphysics modeling framework that integrates optical, thermal, and mechanical analyses to optimize CMOS-SOI-MEMS TMOS sensors with Metasurface absorbers. The following subsections synthesize the results, elucidate trade-offs, and highlight the framework’s utility for rapid design optimization, particularly through the efficient thermal analytical model.

FDTD simulations demonstrate that Metasurface absorbers achieve near-unity absorption at targeted wavelengths (e.g., 0.99 at 4.26 µm for MIM 1, compared to 0.7564 with TiN layer and 0.032 without) (Section 3.3.2). The TiN layer enhances broadband absorption across 2–14 µm, while MIM configurations enable spectral selectivity for applications like CO₂ sensing (4.26 µm) and thermal imaging (10 µm). The optical modeling framework’s versatility, applicable to any absorber architecture-narrowband, broadband, multilayer, or varied geometries-ensures broad utility for tailored IR sensor designs.

Thermal FEA confirms the TMOS sensor’s high thermal isolation, maximizing temperature rise (ΔT) for a given IR flux (Section 4). Metasurface integration, however, increases thermal capacitance, extending the thermal time constant from 80 ms (TiN) to 95 ms (MIM 1) and 129 ms (MIM 2). While suitable for optical gas sensing, this may limit temporal bandwidth in fast-detection scenarios, necessitating careful absorber optimization. A key advancement is the analytical lumped-element RC circuit model, which exhibits excellent agreement with 3D FEA (e.g., 78.93 ms vs. 80.80 ms for TiN) (Section 4.3). This model enables rapid evaluation of diverse layer and Metasurface absorber structures—narrowband, broadband, or multilayer—without the computational overhead of full 3D FEA simulations, which are significantly more resource-intensive [34]. Mechanically, modal and harmonic analyses verify structural resilience, with natural frequencies remaining above 20 kHz despite Metasurface mass loading (e.g., 31.03 kHz to 24.88 kHz for MIM 2) (Section 5.2), ensuring stability against environmental vibrations (<1 kHz).

This integrated framework advances beyond isolated domain analyses, providing a predictive tool to navigate trade-offs in absorption efficiency, thermal response, and mechanical stability. The RC circuit model’s computational efficiency facilitates rapid design iterations, enabling optimization of any absorber architecture for CMOS-compatible IR sensors. This versatility supports applications from gas sensing to biomedical diagnostics. Future work will focus on experimental validation with fabricated prototypes to refine model accuracy and explore control strategies to mitigate ambient temperature drifts, further enhancing manufacturability and performance.

7. Conclusions

This work has presented a comprehensive multiphysics modeling framework for CMOS-SOI-MEMS thermal infrared sensors integrated with Metasurface absorbers. By coupling optical, thermal, and mechanical analyses, the framework provides a predictive and quantitative approach for optimizing sensor performance across interdependent domains.

Optical simulations using the FDTD method demonstrated that Metasurface absorbers can achieve near-unity absorption at specific wavelengths, such as 4.26 µm for CO₂ detection and 10 µm for broadband thermal imaging, while maintaining CMOS process compatibility. Thermal FEA confirmed that the TMOS structure offers excellent thermal isolation and responsivity, though Metasurface integration increases thermal capacitance and extends the time constant. The validated analytical RC model effectively reproduces full 3D FEA results, enabling rapid evaluation of absorber geometries and materials with minimal computational cost. Mechanical modal and harmonic analyses verified that the device maintains sufficient stiffness and stability, with all natural frequencies well above environmental vibration ranges.

Overall, the developed framework unifies optical absorption design, thermal performance assessment, and mechanical reliability evaluation within a single simulation environment. This approach not only accelerates design iterations but also clarifies trade-offs between responsivity, speed, and structural integrity. The presented methodology can be extended to a wide range of uncooled IR detector architectures and Metasurface configurations. Future research will focus on experimental validation through fabricated prototypes and the integration of temperature compensation and feedback mechanisms to further enhance accuracy and robustness in real-world sensing applications.