This technical seminar details finite element method (FEM) modeling protocols for simulating laser-induced heating and thermal expansion in optomechanical systems, specifically micro-scale resonators. The presentation delineates the physical distinction between radiation pressure and bolometric (photothermal) forces, emphasizing how finite absorption in dielectrics leads to volumetric heating, structural deformation, and subsequent shifts in resonance frequency. Through COMSOL Multiphysics, the instructor demonstrates the setup of heat transfer modules, the definition of Gaussian heat sources, and the analysis of both stationary temperature gradients and time-dependent thermal relaxation. Furthermore, the session explores the origins of intrinsic stress in silicon nitride thin films resulting from coefficient of thermal expansion (CTE) mismatches during high-temperature fabrication processes. The session concludes by establishing the necessity of frequency-domain analysis to distinguish photothermal back-action from radiation pressure in high-precision sensing applications.
Review Group Recommendation
The ideal audience for this material includes MEMS (Micro-Electromechanical Systems) Design Engineers, Computational Physicists specializing in Nanophotonics, and Graduate Researchers in Optomechanics. These professionals possess the requisite background in partial differential equations and solid-state physics to apply these simulation techniques to device fabrication and sensing characterization.
Technical Summary: Laser Heating and Thermal Expansion in Resonators
0:52 Volumetric Thermal Effects: Beyond radiation pressure, laser absorption in dielectrics (e.g., silicon nitride) triggers volumetric heating. This leads to thermal expansion, which modifies the physical dimensions and resonance frequency of the cavity.
2:44 Competing Optomechanical Forces: Photothermal (bolometric) forces and radiation pressure often compete. While radiation pressure is nearly instantaneous, photothermal effects involve a retardation (diffusion time) that can lead to either cooling or amplification of mechanical motion.
5:13 Governing Equations for Thermal Transport: The simulation relies on the diffusion equation, requiring parameters such as thermal capacitance ($C_p$), thermal conductivity ($k$), and the absorption coefficient ($\alpha$).
8:17 Motivation for Photothermal Simulation: High-power optical probing is required for ground-state cooling of mechanical oscillators, but excessive power induces significant temperature rises (e.g., from 5K to 40K), necessitating accurate thermal modeling to validate experimental hypotheses.
10:48 Modeling Beam Heating: The instructor demonstrates a double-clamped beam simulation. Key parameters include the circulating power in the cavity and the "sampling length" (the ~5–7 $\mu m$ interaction region where the beam sits in the cavity's near-field).
13:44 Mode Volume Fraction: In near-field interactions, only a small fraction (approximately $10^{-3}$) of the total circulating optical power is typically absorbed by the adjacent mechanical beam.
24:36 Defining User-Induced Heat Sources: COMSOL’s "Definitions" node is used to implement analytical Gaussian functions ($W/m^3$) to represent the laser spot's intensity profile across the beam's longitudinal axis.
32:53 Time-Dependent Thermal Dynamics: A transient study reveals the thermal relaxation time—the duration required for phonons to diffuse out of the system. In the provided example, the thermal response cut-off is approximately 20 Hz; modulation faster than this prevents the temperature from following the power fluctuations.
37:30 Origins of Thermal Stress: High intrinsic stress (~1 GPa) in LPCVD silicon nitride is partially attributed to CTE mismatch. When cooling from $1000^\circ C$ to room temperature, the silicon substrate shrinks less than the nitride film, effectively "pulling" the beam and inducing tensile stress.
45:42 Quantifying Fabricated Stress: FEM analysis indicates that thermal expansion mismatch accounts for approximately 500 MPa of stress, leaving the remainder of the 1 GPa total to other amorphous material properties.
49:06 Bolometric Force Frequency Response: In Whispering Gallery Mode (WGM) resonators, the bolometric force exhibits a flat response at low frequencies but drops off at a specific pole (cutoff frequency). Engineers must ensure the mechanical resonance frequency of interest is significantly higher than this cutoff to isolate pure radiation pressure effects.