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 planes of 450-nm-diameter SD fiber with oblique endface. 980-nm fundamental mode (i.e., HE11 mode) with quasi-y polarization is excited at z = 0 μm and propagates along SD fiber (whose length occupies the range from z = 0 μm to z = 4 μm). 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 xz plane (y = 0 μm). ⟨Fx⟩, ⟨Fy⟩, and ⟨Fz⟩, respectively in xy plane with z = 1.3 μm (second row, i.e., (d)–(f)), z = 1.92 μm (third row, i.e., (g)–(i)), z = 2.74 μm (fourth second row, i.e., (j)–(l)), z = 3.255 μm (bottom row, i.e., (m)–(o)). The dashed line represents the profile of the SD fiber.