This paper presents a new approach to the simulation of the thermal reflow of e-beam-exposed polymethyl methacrylate (PMMA) taking into account its nonuniform viscosity profile. This approach is based on numerical "soapfilm" modeling of the surface evolution, processed by the free software "Surface Evolver" in area normalization mode. The PMMA viscosity profile is calculated via the simulation of the exposed PMMA number average molecular weight distribution using the Monte-Carlo method and empirical formulas. The relation between the PMMA viscosity and the mobility of PMMA surface vertices was determined via the thermal reflow simulation for uniform PMMA gratings using analytical and numerical approaches in a wide viscosity range. The agreement between reflowed profiles simulated with these two approaches emphasizes the applicability of "soapfilm" modeling in the simulation of polymer thermal reflow. The inverse mobility of PMMA surface vertices appeared to be proportional to the PMMA viscosity with a high precision. The developed approach enables thermal reflow simulations for complex nonuniform structures, which allows the use of predictable reflow as a stage of 3D microfabrication.
Keywords: PMMA; grayscale e-beam lithography; nonuniform viscosity profile; thermal reflow.