Theoretical analysis of optical force density distribution inside subwavelength-diameter optical fibers*

Project supported by the National Natural Science Foundation of China (Grant Nos. 11604230 and 11434017), the Guangdong Provincial Innovative and Entrepreneurial Research Team Program, China (Grant No. 2016ZT06C594), and the National Key Research and Development Program of China (Grant No. 2018YFA 0306200).
Zhang Yun-Yuan1, Yu Hua-Kang1, †, Wang Xiang-Ke1, Wu Wan-Ling1, Gu Fu-Xing2, Li Zhi-Yuan1, ‡
       

(color online) Cross-sectional plots of time-averaged optical force-density components ⟨Fx⟩, ⟨Fy⟩, and ⟨Fz⟩ through central plane of 450-nm-diameter SD fiber with flat endface. A 980-nm fundamental mode (i.e., HE11 mode) with quasi-y polarization is excited at z = 0 μm and propagates in positive z-direction along the SD fiber. The plots of calculated optical force density correspond to incident optical power Pinc = 3.06 × 10−16 W. Color scale bar shows force density in units of μN/m3. (a) ⟨Fx⟩, (b) ⟨Fy⟩, (c) ⟨Fz⟩, respectively in central xy plane with z = 2 μm. (d) ⟨Fx⟩, (e) ⟨Fy⟩, (f) ⟨Fz⟩, respectively in central xz plane with y = 0 μm. (g) ⟨Fx⟩, (h) ⟨Fy⟩, (i) ⟨Fz⟩, respectively in central yz plane with x = 0 μm. Dashed lines represent profile of SD fiber.