基于光纤光镊的长距离微粒可控操纵

基于光纤光镊的长距离微粒可控操纵 Controllable manipulation of long-distance microparticles based on fiber optic tweezers first-author 姜春雷（1977-），男，黑龙江青冈人，教授，博士生导师，2006年于大庆石油学院获得硕士学位，2017年于哈尔滨工业大学获得博士学位，主要从事微纳光纤传感器、激光超精密测量与超显微成像等方面的理论及应用研究。 E-mail： jiangchunlei_nepu@163.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965645&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965647&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965646&type= 22.01332855 32.00399780 姜春雷 （1977-），男，黑龙江青冈人，教授，博士生导师，2006年于大庆石油学院获得硕士学位，2017年于哈尔滨工业大学获得博士学位，主要从事微纳光纤传感器、激光超精密测量与超显微成像等方面的理论及应用研究。 E-mail： jiangchunlei_nepu@163.com jiangchunlei_nepu@163.com corresp E-mail：eq0687@126.com E-mail：eq0687@126.com 陈　朋（1981-），男，安徽临泉人，博士，讲师，2010年，2015于哈尔滨工业大学分别获得硕士和博士学位，主要从事精密测量、外差激光干涉测量方面的理论及应用研究。E-mail：eq0687@126.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965648&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965650&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965649&type= 22.01332855 32.00399780 陈　朋 （1981-），男，安徽临泉人，博士，讲师，2010年，2015于哈尔滨工业大学分别获得硕士和博士学位，主要从事精密测量、外差激光干涉测量方面的理论及应用研究。E-mail： eq0687@126.com eq0687@126.com 东北石油大学 电气信息工程学学院，黑龙江 大庆 163318 东北石油大学 电气信息工程学学院，黑龙江 大庆 163318 College of Electrical Information Engineering， Northeast Petroleum University， Daqing 163318， China College of Electrical Information Engineering， Northeast Petroleum University， Daqing 163318， China 在液体中灵活操纵微粒或细胞，特别是将细胞或微粒运输到指定位置，已经被证明在细胞分析、疾病诊断、药物递送等方面有着至关重要的作用。针对细胞或微粒非接触光学捕获的灵活性受限于光纤光镊操纵距离短的问题，提出了一种结构简单且具有长距离非接触可控操纵微粒的新型光纤光镊。利用加热和拉伸技术制作了类锥形平口光纤探针，980 nm的激光经过光纤探针后会对微粒产生大的散射力，将微粒逐渐推离光纤端口，同时借助反向流体阻力，在不移动光纤探针的情况下通过调节激光器光源的输出功率，对轴向位置上直径为6 µm的聚苯乙烯微粒可以进行长达102.2 µm的可控往返操纵，应用有限元法仿真了光镊的光场强度分布，并采用麦克斯韦应力张量法分析了光镊对微粒的作用力。实验和仿真结果表明，所提出的类锥形平口光镊是可行的。 The flexible manipulation of particles or cells in fluids， particularly the transport of cells or particles to a specified location， has proven to be of vital importance in cell analysis， disease diagnosis， drug delivery， etc. However， the flexibility of contact-free optical capture of cells or particles is clearly insufficient owing to the short operating distance of optical fiber tweezers. Therefore， a novel optical fiber tweezer with a simple structure and controllable long-distance contact-free control of particles is proposed. In this study， a special cone-like flat port fiber probe is fabricated via heating and stretching techniques. A 980 nm laser passing through the fiber probe exerts a large scattering force on the particles， which gradually pushes the particles away from the fiber port. Simultaneously， with the aid of reverse fluid resistance， the output power of the laser light source is adjusted without moving the fiber probe. A polystyrene particle with a diameter of 6 µm in an axial position can be controlled round-trip with a distance of 102.2 µm. The intensity distribution of the optical tweezers is simulated using the finite element method， and the force exerted by the optical tweezers on the particles is analyzed using the Maxwell stress tensor method. The feasibility of the proposed cone-like flat port optical tweezers is verified via both experiments and simulations. 光纤光镊 光学捕获 有限元法 麦克斯韦应力张量法 optical fiber tweezers optical trapping finite element method Maxwell stress tensor method 1 引 言 1 微纳米颗粒或生物细胞的光学操纵一直是研究热点 ［ 1 1 - 3 3 1-3 ］ 。自从美国贝尔实验室的Ashkin利用单束激光引入高数值孔径物镜形成三维光学势阱，实现对微粒的非接触三维空间俘获后，光镊技术开创了光学微观操控研究的新纪元。然而，传统的光镊利用散装光学元件通过相位调制产生多个梯度光学势阱，系统结构复杂，且大数值孔径物镜的工作距离较低，严重限制了对微粒或细胞的光学操纵。为了克服这些限制，研究人员采用光纤光镊对细胞或微粒进行操纵 ［ 4 4 - 5 5 4-5 ］ 。与双光束光纤镊子 ［ 6 6 - 8 8 6-8 ］ 和微纳米光纤的光学操作 ［ 9 9 - 10 10 9-10 ］ 相比，基于单光纤尖端的锥形光纤光镊紧凑灵活，具有更大的应用潜力。但锥形光纤光镊输出光强烈聚焦在光纤尖端附近，导致粒子或细胞的操纵距离受到限制 ［ 11 11 - 14 14 11-14 ］ ，如同光学显微镜中工作距离短的物镜，使得光纤光镊的使用变得不灵活，并且由于光纤光镊的操作距离短，细胞靠近光纤端口，可能会对生物样品造成损伤。因此，具有长距离操纵的光纤光镊对于实现真正的非接触光学操纵具有重要意义。 在过去的十几年，研究人员对如何产生具有长距离操纵的光纤光镊展开了大量研究。2012年，李等利用亚微米锥形光纤探针实现了高度灵活的粒子捕获、驱动和精确排列 ［ 15 15 ］ 。2013年，Liu等利用模分复用技术实现了单光纤捕获系统中被捕获微粒的轴向捕获位置动态调整 ［ 16 16 ］ 。通过改变基模光束（LP 01 ）和低阶模光束（LP 11 ）两种模式光束的功率比，控制捕获的直径为6 μm酵母细胞沿光轴方向移动3 μm。2015年，张等提出一种梯度折射率光纤镊子 ［ 17 17 - 19 19 17-19 ］ ，在单模光纤和梯度折射率光纤之间引入空气腔，通过调节激光的输出功率、空气腔的长度以及流速来改变操作距离。空气腔是由两种不同的光纤对齐拼接而成，在调整空气腔的长度时有些许复杂。2019年，Liu等利用倏逝波稳定地将纳米颗粒和细胞捕获在由天然生物细胞组装而成的光学传送带表面 ［ 20 20 ］ 。通过调节注入激光的相对功率，粒子或细胞可以沿着生物传送带实现双向传输。引入两个反向传播光束的双光束光阱可能对粒子或细胞的传输距离存在限制。2021年，Mumati等提出一种用波长2 µm的掺铥光纤激光器进行光热泳来俘获微粒的方法 ［ 21 21 ］ 。由于激光束被水溶液强烈吸收，在焦点周围产生局部温度梯度，最终颗粒沿温度梯度迁移并聚集在距焦点125 µm的地方。如果应用在生物细胞和组织运输领域，这种热效应 ［ 22 22 ］ 可能会对生物细胞和组织造成损伤。 为了实现结构简单的长距离可控操纵微粒或细胞的光纤光镊，本文利用改造的类锥形平口光纤探针在980 nm的激光输出功率较低的情况下，提供一个大的散射力将粒子推离光纤端口，同时利用流体的阻力，从而可以对微粒进行长距离可控往返操纵。采用有限元分析法建立模型仿真类锥形平口光纤的光场分布，并结合麦克斯韦应力张量法推导出光阱力公式，分析捕获微粒的受力情况，通过仿真验证了该方法的可行性。 2 原 理 1 在光纤端口，微球主要受到两个力的作用：梯度力和散射力。梯度力是微球中电偶极矩在不均匀电磁场中受到的力，它正比于光强梯度的平方，指向光场强度的最大处，使微球被稳定俘获到光纤表面。光在散射过程中与光子交换动量获得散射力，其方向沿着光的传播方向，使微球沿着光束的传播方向运动。梯度力和散射力的合力被称为捕获力。捕获力与光的波长、粒子属性及粒子尺寸等因素密切相关。本文使用有限元分析法，通过解麦克斯韦方程对矢量光场的分布和变化进行精确描述。电磁场作用在体积为 V 粒子的光阱力 <math id="M1"><mi>F</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965652&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965651&type= 2.53999996 2.79399991 为： <math display="block" id="M2"><mi>F</mi><mo>=</mo><msub><mo>∫</mo><mrow><mi>V</mi></mrow></msub><mrow><mo>∇</mo><mo>⋅</mo><mi mathvariant="bold-italic">T</mi><mi mathvariant="normal">d</mi><mi>V</mi><mo>=</mo><msub><mo>∫</mo><mrow><mi>S</mi></mrow></msub><mrow><mo stretchy="false">(</mo><mi mathvariant="bold-italic">T</mi><mo>⋅</mo><mi mathvariant="bold-italic">n</mi><mo stretchy="false">)</mo><mi mathvariant="normal">d</mi><mi>S</mi></mrow></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965655&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965653&type= 46.65133286 6.60400009 . 通过高斯定理（或散度定理）可知，作用在体积 <math id="M3"><mi>V</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965659&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965657&type= 3.04800010 2.87866688 粒子的作用力可换算为作用其闭合曲 <math id="M4"><mi>S</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965717&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965716&type= 2.20133328 2.79399991 上的力。本文采用麦克斯韦应力张量法 ［ 23 23 ］ ，真空中麦克斯韦应力张量的表达式为： <math display="block" id="M5"><mi mathvariant="bold-italic">T</mi><mo>=</mo><msub><mrow><mi>ε</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mi mathvariant="bold-italic">E</mi><mi mathvariant="bold-italic">E</mi><mo>+</mo><msub><mrow><mi>μ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mi mathvariant="bold-italic">H</mi><mi mathvariant="bold-italic">H</mi><mo>-</mo><mfenced separators="|"><mrow><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><msub><mrow><mi>ε</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msup><mrow><mi mathvariant="bold-italic">E</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><msub><mrow><mi>μ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msup><mrow><mi mathvariant="bold-italic">H</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfenced><mi mathvariant="bold-italic">I</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965665&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965664&type= 67.39466858 9.05933285 ， 其中： <math id="M6"><msub><mrow><mi>ε</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965667&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965666&type= 2.62466669 3.72533321 是真空中介电常数， <math id="M7"><msub><mrow><mi>μ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965669&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965668&type= 3.13266683 4.57200003 是真空中磁感应系数。麦克斯韦应力张量 <math id="M8"><mi mathvariant="bold-italic">T</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965685&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965684&type= 2.62466669 2.79399991 有9个分量： <math display="block" id="M9"><mi mathvariant="bold-italic">T</mi><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><msub><mrow><mi mathvariant="bold-italic">a</mi></mrow><mrow><mn mathvariant="normal">11</mn></mrow></msub></mtd><mtd><msub><mrow><mi mathvariant="bold-italic">a</mi></mrow><mrow><mn mathvariant="normal">12</mn></mrow></msub></mtd><mtd><msub><mrow><mi mathvariant="bold-italic">a</mi></mrow><mrow><mn mathvariant="normal">13</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi mathvariant="bold-italic">a</mi></mrow><mrow><mn mathvariant="normal">21</mn></mrow></msub></mtd><mtd><msub><mrow><mi mathvariant="bold-italic">a</mi></mrow><mrow><mn mathvariant="normal">22</mn></mrow></msub></mtd><mtd><msub><mrow><mi mathvariant="bold-italic">a</mi></mrow><mrow><mn mathvariant="normal">23</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi mathvariant="bold-italic">a</mi></mrow><mrow><mn mathvariant="normal">31</mn></mrow></msub></mtd><mtd><msub><mrow><mi mathvariant="bold-italic">a</mi></mrow><mrow><mn mathvariant="normal">32</mn></mrow></msub></mtd><mtd><msub><mrow><mi mathvariant="bold-italic">a</mi></mrow><mrow><mn mathvariant="normal">33</mn></mrow></msub></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965673&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965672&type= 30.14133072 13.54666710 . 这9个分量具体为： <math display="block" id="M10"><mfenced open="{" close="" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><msub><mrow><mi mathvariant="bold-italic">a</mi></mrow><mrow><mn mathvariant="normal">11</mn></mrow></msub><mo>=</mo><msub><mrow><mi>ε</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msubsup><mrow><mi mathvariant="bold-italic">E</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>+</mo><msub><mrow><mi>μ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msubsup><mrow><mi mathvariant="bold-italic">H</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>-</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mfenced separators="|"><mrow><msub><mrow><mi>ε</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msup><mrow><mi mathvariant="bold-italic">E</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msub><mrow><mi>μ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msup><mrow><mi mathvariant="bold-italic">H</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfenced></mtd></mtr><mtr><mtd><msub><mrow><mi mathvariant="bold-italic">a</mi></mrow><mrow><mn mathvariant="normal">12</mn></mrow></msub><mo>=</mo><msub><mrow><mi>ε</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">E</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">E</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>+</mo><msub><mrow><mi>μ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">H</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">H</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi mathvariant="bold-italic">a</mi></mrow><mrow><mn mathvariant="normal">13</mn></mrow></msub><mo>=</mo><msub><mrow><mi>ε</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">E</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">E</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><mo>+</mo><msub><mrow><mi>μ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">H</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">H</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi mathvariant="bold-italic">a</mi></mrow><mrow><mn mathvariant="normal">21</mn></mrow></msub><mo>=</mo><msub><mrow><mi>ε</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">E</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">E</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>+</mo><msub><mrow><mi>μ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">H</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">H</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi mathvariant="bold-italic">a</mi></mrow><mrow><mn mathvariant="normal">22</mn></mrow></msub><mo>=</mo><msub><mrow><mi>ε</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msubsup><mrow><mi mathvariant="bold-italic">E</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>+</mo><msub><mrow><mi>μ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msubsup><mrow><mi mathvariant="bold-italic">H</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>-</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mfenced separators="|"><mrow><msub><mrow><mi>ε</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msup><mrow><mi mathvariant="bold-italic">E</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msub><mrow><mi>μ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msup><mrow><mi mathvariant="bold-italic">H</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfenced></mtd></mtr><mtr><mtd><msub><mrow><mi mathvariant="bold-italic">a</mi></mrow><mrow><mn mathvariant="normal">23</mn></mrow></msub><mo>=</mo><msub><mrow><mi>ε</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">E</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">E</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>+</mo><msub><mrow><mi>μ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">H</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">H</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi mathvariant="bold-italic">a</mi></mrow><mrow><mn mathvariant="normal">31</mn></mrow></msub><mo>=</mo><msub><mrow><mi>ε</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">E</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">E</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>+</mo><msub><mrow><mi>μ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">H</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">H</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi mathvariant="bold-italic">a</mi></mrow><mrow><mn mathvariant="normal">32</mn></mrow></msub><mo>=</mo><msub><mrow><mi>ε</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">E</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">E</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>+</mo><msub><mrow><mi>μ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">H</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">H</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi mathvariant="bold-italic">a</mi></mrow><mrow><mn mathvariant="normal">33</mn></mrow></msub><mo>=</mo><msub><mrow><mi>ε</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msubsup><mrow><mi mathvariant="bold-italic">E</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>+</mo><msub><mrow><mi>μ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msubsup><mrow><mi mathvariant="bold-italic">H</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>-</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mfenced separators="|"><mrow><msub><mrow><mi>ε</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msup><mrow><mi mathvariant="bold-italic">E</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msub><mrow><mi>μ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msup><mrow><mi mathvariant="bold-italic">H</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfenced></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965675&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965674&type= 62.39933777 63.83866882 ， 其中： <math id="M11"><msub><mrow><mi>ε</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965677&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965676&type= 2.62466669 3.72533321 是介质中介电常数， <math id="M12"><msub><mrow><mi>μ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965679&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965678&type= 3.13266683 4.57200003 是介质中磁感应系数， <math id="M13"><mi mathvariant="bold-italic">E</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965741&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965738&type= 2.62466669 2.79399991 是电场矢量， <math id="M14"><mi mathvariant="bold-italic">H</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965753&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965750&type= 3.04800010 2.79399991 是磁场矢量。数字1，2，3分别对应笛卡尔坐标系中的 x ， y ， z 三个分量。 <math id="M15"><mi mathvariant="bold-italic">T</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965685&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965684&type= 2.62466669 2.79399991 的每个分量可表示为： <math display="block" id="M16"><mtable><mtr><mtd><mfenced open="<" close=">" separators="|"><mrow><msub><mrow><mi mathvariant="bold-italic">T</mi></mrow><mrow><mi>i</mi><mi>j</mi></mrow></msub></mrow></mfenced><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mi mathvariant="normal">R</mi><mi mathvariant="normal">e</mi><mfenced open="[" close="" separators="|"><mrow><msub><mrow><mi>ε</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msub><mrow><mi>ε</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">E</mi></mrow><mrow><mi mathvariant="bold-italic">i</mi></mrow></msub><msubsup><mrow><mi mathvariant="bold-italic">E</mi></mrow><mrow><mi mathvariant="bold-italic">j</mi></mrow><mrow><mi mathvariant="bold-italic">*</mi></mrow></msubsup><mo>-</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><msub><mrow><mi>δ</mi></mrow><mrow><mi>i</mi><mi>j</mi></mrow></msub><mfenced separators="|"><mrow><msub><mrow><mi>ε</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msub><mrow><mi>ε</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msup><mrow><mi mathvariant="bold-italic">E</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfenced><mo>+</mo></mrow></mfenced></mtd></mtr><mtr><mtd><mfenced open="" close="]" separators="|"><mrow><msub><mrow><mi>μ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msub><mrow><mi>μ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">H</mi></mrow><mrow><mi mathvariant="bold-italic">i</mi></mrow></msub><msubsup><mrow><mi mathvariant="bold-italic">H</mi></mrow><mrow><mi mathvariant="bold-italic">j</mi></mrow><mrow><mi mathvariant="bold-italic">*</mi></mrow></msubsup><mo>-</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><msub><mrow><mi>δ</mi></mrow><mrow><mi>i</mi><mi>j</mi></mrow></msub><mfenced separators="|"><mrow><msub><mrow><mi>μ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msub><mrow><mi>μ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msup><mrow><mi mathvariant="bold-italic">H</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfenced></mrow></mfenced><mo>,</mo></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965687&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965686&type= 64.34667206 19.98133469 其中： <math id="M17"><msub><mrow><mi mathvariant="bold-italic">T</mi></mrow><mrow><mi>i</mi><mi>j</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965689&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965688&type= 3.97933316 3.72533321 表示作用在垂直于 <math id="M18"><mi>j</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965699&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965698&type= 1.10066664 3.64066648 轴单位面积上的力在 <math id="M19"><mi>i</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965694&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965693&type= 1.10066664 2.79399991 轴上的分量。 <math id="M20"><msub><mrow><mi mathvariant="bold-italic">n</mi></mrow><mrow><mi>j</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965697&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965695&type= 2.87866688 3.72533321 是垂直于 <math id="M21"><mi>j</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965699&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965698&type= 1.10066664 3.64066648 轴的微粒外表面 <math id="M22"><mi>S</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965717&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965716&type= 2.20133328 2.79399991 的向外法向量。其中： <math display="block" id="M23"><msub><mrow><mi>δ</mi></mrow><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>=</mo><mfenced open="{" close="" separators="|"><mrow><mtable><mtr><mtd><mn mathvariant="normal">0</mn><mo>,</mo><mi>i</mi><mo>≠</mo><mi>j</mi></mtd></mtr><mtr><mtd><mn mathvariant="normal">1</mn><mo>,</mo><mi>i</mi><mo>=</mo><mi>j</mi></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965705&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965703&type= 20.74333191 10.32933331 . 因此，粒子在轴向上受到的光学力 <math id="M24"><msub><mrow><mi>F</mi></mrow><mrow><mi mathvariant="normal">a</mi><mi mathvariant="normal">o</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965709&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965707&type= 4.57200003 3.72533321 可表示为 ［ 24 24 ］ ： <math display="block" id="M25"><msub><mrow><mi>F</mi></mrow><mrow><mi mathvariant="normal">a</mi><mi mathvariant="normal">o</mi></mrow></msub><mo>=</mo><msub><mo>∫</mo><mrow><mi>S</mi></mrow></msub><mrow><mfenced separators="|"><mrow><mfenced open="<" close=">" separators="|"><mrow><msub><mrow><mi mathvariant="bold-italic">T</mi></mrow><mrow><mi mathvariant="bold-italic">M</mi></mrow></msub></mrow></mfenced><mo>⋅</mo><mi mathvariant="bold-italic">n</mi></mrow></mfenced><mi mathvariant="normal">d</mi><mi>S</mi></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965713&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965712&type= 32.85066605 6.60400009 ， 其中：对粒子周围的封闭曲面 <math id="M26"><mi>S</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965717&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965716&type= 2.20133328 2.79399991 进行积分， <math id="M27"><mi mathvariant="bold-italic">n</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965723&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965720&type= 2.11666679 2.79399991 为曲面法向量， <math id="M28"><mfenced open="<" close=">" separators="|"><mrow><msub><mrow><mi mathvariant="bold-italic">T</mi></mrow><mrow><mi mathvariant="bold-italic">M</mi></mrow></msub></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965729&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965726&type= 7.78933382 4.06400013 为与时间无关的麦克斯韦应力张量，其计算方法如下： <math display="block" id="M29"><mfenced open="<" close=">" separators="|"><mrow><msub><mrow><mi mathvariant="bold-italic">T</mi></mrow><mrow><mi mathvariant="bold-italic">M</mi></mrow></msub></mrow></mfenced><mo>=</mo><mi mathvariant="bold-italic">D</mi><msup><mrow><mi mathvariant="bold-italic">E</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup><mo>+</mo><mi mathvariant="bold-italic">H</mi><msup><mrow><mi mathvariant="bold-italic">B</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup><mo>-</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mfenced separators="|"><mrow><mi mathvariant="bold-italic">D</mi><msup><mrow><mi mathvariant="bold-italic">E</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup><mo>+</mo><mi mathvariant="bold-italic">H</mi><msup><mrow><mi mathvariant="bold-italic">B</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup></mrow></mfenced><mi mathvariant="bold-italic">I</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965735&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965732&type= 62.65333176 7.87400007 ， 式中： <math id="M30"><mi mathvariant="bold-italic">E</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965741&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965738&type= 2.62466669 2.79399991 是电场矢量， <math id="M31"><mi mathvariant="bold-italic">D</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965747&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965744&type= 2.87866688 2.79399991 是电位移矢量， <math id="M32"><mi mathvariant="bold-italic">H</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965753&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965750&type= 3.04800010 2.79399991 是磁场矢量， <math id="M33"><mi mathvariant="bold-italic">B</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965758&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965755&type= 2.53999996 2.79399991 是磁通量密度矢量，上标“ <math id="M34"><mi mathvariant="normal">*</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965764&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965762&type= 1.77800000 2.79399991 ”表示共轭复数， <math id="M35"><mi mathvariant="bold-italic">I</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965771&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965769&type= 1.43933344 2.79399991 是各向同性张量。 普通抛物线形或锥形光纤在捕获微粒时，由于端口类似凸透镜，从而形成强汇聚场，微粒通常被捕获在光纤端口或者附近。如果对粒子或细胞进行长距离可控操纵，微粒可以双向运输，需要一个特定的光纤探针沿光轴提供一个大的光散射力，而不是将细胞捕获在光纤端口的光梯度力。本文提出了一个类锥形平口状的光纤探针，该光纤端口两侧结构类似于抛物线形，避免了光束在通过端口后向四周发散，降低了能量损失。其次端口呈现平口状，无法形成有效的强汇聚场，从而在轴向上对微粒始终表现为远离光纤探针尖端的光学力。为了验证抛物线形光纤与类锥形平口状光纤对微粒捕获的影响，本文通过有限元分析建立了一个仿真模型来分别对这两种不同结构的光纤进行仿真分析。两种不同光纤探针的仿真结果如 图1 图1 所示。 仿真中，水、光纤探针、聚苯乙烯微粒的折射率分别设置为1.33，1.46，1.59。聚苯乙烯微粒的半径设置为3 μm，激光的波长设置为980 nm，光纤端口的输入功率 P in =5 W/m。其中， 图1 图1 （a）为抛物线形光纤的光场分布仿真和微粒在轴向上的受力分析； 图1 图1 （b）为类锥形平口光纤的光场分布仿真和微粒在轴向上的受力分析。从 图1 图1 可以看出，光经过抛物线形光纤探针作用后的光场分布在探针端口附近，聚苯乙烯微粒在靠近光纤端口处受到一个负的光学力 F ao ，其作用是将粒子吸引并束缚在光纤端口，导致粒子无法进行长距离可控操纵，与本文实验要求不符。而光经过类锥形平口光纤探后汇聚的光场分布在离光纤端口约10 μm处位置，微粒在轴向位置上受到一个正的光学力 F ao ，其作用是将粒子推离光纤端口，使得微粒存在双向运输的可能性。结合以上仿真结果，类锥形平口状光纤更符合本文的实验要求。 本文采用拉制的类锥形平口状光纤探针进行实验。 图2 图2 为实验中使用的光纤探头及其结构参数。它是由一种普通单模光纤G.652D（康宁SMF-28e）经过加热和拉伸制作而成。首先，使用光纤剥离器将光纤最外层保护套和涂覆层剥离，再把光纤放置在酒精灯火焰下加热至光纤的熔点，加热区域沿光纤方向约为3~4 mm，加热约5 s后，以大约为0.03 mm/s的速度水平拉伸纤维，拉伸3 s后，拉伸速度突变为0.05 mm/s，使光纤断裂。由于熔融石英材料表面的张力作用，光纤形成类锥形平口状端口。该方法具有很高的重复性，有助于大规模的生产及应用。 实验采用中心波长 <math id="M36"><mi>λ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965794&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965792&type= 1.77800000 2.79399991 =980 nm的激光器作为光源，其光源输出功率 P 可在0~20 mW内调节。在这个波长范围内，水和微粒的吸收都可以忽略。首先，将波长为980 nm激光器的光直接耦合到由普通单模光纤G.652D（康宁SMF-28e）拉制而成的光纤探头中，然后将光纤探头固定在精密三维平台上，从而实现光纤探针的精确定位。用不锈钢毛细管（内径约为0.9 mm，壁厚约为0.1 mm，长度约为120 mm）包裹光纤探头，防止光纤探头断裂和弯曲。随后，将光纤探头浸入直径为6 µm的聚苯乙烯（Polystyrene）微粒的悬浮液中，而装有聚苯乙烯微粒的悬浮液则放在精密三维台上。实验采用集成电荷耦合器件（Charge Coupled Device，CCD）摄像机、光学变焦镜头（Optem zoom 70XL）、物镜（20X，NA=0.40）以及增距镜（2×）的组合来观察细胞的光学捕获，并通过计算机进行实时监控，包括图像采集和视频记录。 图3 图3 为聚苯乙烯（PS）微粒长距离可控操纵实验装置。 实验利用溶液的蒸发力以及分子间的引力产生一个与光学力平衡的流体阻力 F v ，实验模型如 图4 图4 （a）所示。在载玻片的4个角落分别粘有少量的蓝凝胶（用于固定载玻片和盖玻片），然后蓝凝胶上放置一层盖玻片，在盖玻片与载玻片之间注入配置好的聚苯乙烯（PS）溶液。需要注意的是，放置光纤探头处的端口高度要低于相对一侧端口的高度，其余两侧保持平行。其目的是减小放置光纤探针端口处溶液与空气的接触面积，从而降低该端口的蒸发力。放置光纤探头的一侧高度约为1.5 mm，相对一侧的高度约为3.0 mm。流体流向的控制原理如 图4 图4 （b）所示，图中 F e 代表左右两侧的蒸发力。因为放置光纤探针的左侧溶液与空气的接触面积要小于相对一侧溶液与空气的接触面积，所以放置光纤探头的左侧蒸发力要小于右侧的蒸发力，同时由于溶液与气体接触的表面层中的分子比液体内部稀疏，其分子引力小于溶液内部的分子引力，从而表面层中的分子受到指向溶液内部力的作用 ［ 25 25 - 27 27 25-27 ］ 。在张力和蒸发力的作用下，水会带动粒子向左侧运动，从而提供一个与光学力 F ao 方向相反的流体阻力 F v ，通过两者之间的平衡，在光轴上实现微粒的捕获。对于球形细胞或微粒来说，流体阻力 F v 可根据斯托克斯定律 ［ 28 28 ］ 表示为： <math display="block" id="M37"><msub><mrow><mi>F</mi></mrow><mrow><mi>v</mi></mrow></msub><mo>=</mo><mn mathvariant="normal">6</mn><mi mathvariant="normal">π</mi><mo>⋅</mo><mi>η</mi><mo>⋅</mo><mi>r</mi><mo>⋅</mo><mi>v</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965810&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965807&type= 24.04533195 4.57200003 ， 其中： <math id="M38"><mi>η</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965816&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965814&type= 1.86266661 3.64066648 是溶液的黏度系数， <math id="M39"><mi>r</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965822&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965819&type= 1.52400005 2.79399991 是微粒半径， <math id="M40"><mi>v</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965828&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965825&type= 1.86266661 2.79399991 是流体的流速。 基于光纤光镊的聚苯乙烯微粒可控操纵工作原理如 图4 图4 （c）所示。当细胞或微粒处于光纤探头端口处，在横向上，横向梯度力 F tg 负责将微粒或细胞约束在光轴方向上，使微粒不发生偏离，从而实现三维空间捕获。在轴向上会受到轴向光学力 F ao 的作用，而轴向光学力 F ao 是由散射力 F as 和轴向梯度力 F ag 共同作用的合力，它主要表现为将微粒或细胞推离光纤端口。其中，横向梯度力 F tg 和轴向光学力 F ao 主要通过激光器光源的输出功率 P 来控制，并随着 P 的增大而增大，沿光传播方向逐渐减小。而流体阻力 F v 的流向主要通过调整放置光纤探头一侧的高度和相对一侧的高度，它与轴向光学力 F ao 之间的平衡决定了微粒的轴向位置。 3 实验结果及分析 1 3.1　光纤探针A的实验 2 本文使用光纤探头A在不同光源输出功率下对直径为6 μm的聚苯乙烯微粒进行可调节操纵，其光学捕获过程的实验显微图像如 图5 图5 所示。其中，通过观察视频中光纤附近参考微粒通过单位距离内的用时来计算当前的流速 v ，共计测量了5组记录，最后取其平均值计算出当前溶液中粒子的平均运动速度 v =14.7 μm/s。本文将微粒中心与光纤端口的距离定义为微粒的操纵距离 L 。当 t =0 s时，如 图5 图5 （a）所示，此时光源输出功率 P 为2.3 mW，聚苯乙烯微粒在 L =11.1 μm的位置被稳定捕获，并作为微粒的初始捕获位置，逐渐增大光源输出功率并观察微粒移动情况。当光源输出功率 P 为3.5，4.5，5.6 mW时，如图 5 5 （b）~ 5 5 （d）所示，粒子的操纵距离 L 分别为13.3，28.9和36.7 μm；当光源输出功率 P 为8.0，10.2，12.5 mW时，如图 6 6 （e）~ 6 6 （f）所示，操纵距离 L 分别为48.3，57.2和62.8 μm。根据实验数据可得，当微粒处于光纤端口附近时，光学力变化梯度较为明显，微粒的移速较快；随着微粒逐渐远离光纤端口，微粒的移动趋势减缓，最终当 t =13 s时，如 图6 图6 （h）所示，测得最大操纵距离 L =65.6 μm，此时激光器输出功率 P =15.0 mW。随后降低激光源的输出功率，发现操纵距离 L 随光源输出功率 P 的降低而逐渐减小，直至恢复初始距离。由此可知，通过调节激光器的输出功率可以操控聚苯乙烯微粒进行往返运动。实验的具体调节过程可参考附件视频，该视频是在相对稳定的状态下单独录制的。 3.2　光纤探针B的实验 2 虽然在第一次实验中最大操纵距离 L 仅被调节到65.56 μm，但从理论上猜想，通过优化光纤探头的几何形状同样可以提高操纵距离 L 。因此，采用光纤探针进行对比实验。对比实验演示了光纤探头B长距离操控聚苯乙烯微粒进行往返运动，结果如 图6 图6 所示。在 t =0 s时，如 图6 图6 （a）所示，初始操纵距离 L =23.3 μm，此时激光器的输出功率 P =2.3 mW。由此发现，在光源输出功率相同的情况下，经过优化后的光纤探头B的初始操纵距离约为光纤探头A的一倍。随后，调节激光光源的输出功率 P 为3.5，4.5，5.6，6.8 mW，如图 6 6 （b）~ 6 6 （e）所示，光纤探头B对微粒的操纵距离 L 分别为32.2，38.9，43.3，66.7 μm。然后继续调节激光光源的输出功率 P 为8.0，9.1，10.2，13.7 mW，如图 6 6 （f）~ 6 6 （i）所示，操纵距离 L 分别为80.0，86.7，94.4，97.8 μm。最终当激光器光源的输出功率 P =15.9 mW时，如 图6 图6 （j）所示，测得光纤探头B的最大操纵距离 L =102.2 μm。根据实验数据可得，光纤探头B的最大操纵距离远大于光纤探头A的最大操纵距离，其约为光纤探头A的1.56倍，进而可以通过调整光纤探头端口的尺寸来调整粒子的操纵距离。 3.3　分析讨论 2 两次实验结果得到激光器输出功率 P 与聚苯乙烯微粒操控距离 L 之间的关系，如 图7 图7 所示。从 图7 图7 可以看出，聚苯乙烯微粒的可控距离 L 随激光器输出功率 P 的提高而增大，并且光纤探头B操纵微粒的距离远大于光纤探头A。由 图2 图2 可知，光纤探头B的光输出端口宽度大约是光纤探头A的2.3倍，光场能够沿光轴覆盖更长的距离，从而使同样尺寸的微粒在相同流速下获得更长的操作距离。为了验证猜想的合理性，采用二维有限元方法分别对光纤探头A 和光纤探头B进行了数值仿真，仿真结果如 图8 图8 所示。由 图8 图8 可知，光通过光纤探头A端口后，在轴向上所形成的光场区域在距离光纤端口约20 μm的位置，光轴的覆盖长度约为30 μm，微粒在靠近光纤探头A端口处受到一个将粒子推离端口的正向光学力 F ao ， F ao 沿光的传播方向逐渐减弱。而光源通过光纤探头B端口后，在轴向上所形成的光场在距离光纤端口约45 μm处，沿光轴的覆盖长度约为50 μm，光纤探针B的光场覆盖长度以及与光纤端口的距离远远大于光纤探针A，所以光纤探针B对同样大小微粒的可控操纵距离大于光纤探针A的操纵距离。结合实验结果和以上仿真分析，本文设计的类锥形平口光纤可以对微粒进行长距离可控操控，通过调整光纤探头端口的尺寸可以调整操纵距离。 4 结 论 1 为了解决光纤光镊捕获细胞或粒子距离短的问题，提出了一种结构简单且具有长距离非接触可控操纵微粒的方法。该方法利用特殊的类锥形平口光纤探头对微粒产生大的散射力，同时与设计的简易流体阻力装置结合在一起，制作了两种不同尺寸端口的类锥形平口光纤并进行对比实验。其中，光纤探头A对微粒的最大操纵距离约为65.6 μm，而在相同光源功率下，光纤探头B对同样大小聚苯乙烯微粒的最大操纵距离可以达到102.2 μm。为了验证实验结果的合理性，本文分别采用有限元法和麦克斯韦应力张量法仿真分析了类锥形平口光纤光镊的光场强度分布以及对微粒的受力情况。仿真结果和实验两方面都验证了所提出类锥形平口光纤光镊实现长距离操纵微粒的可行性。在生物和医学领域，该方法对于生物样品的灵活操控、靶向给药和生物化学物质的协同观察等具有重要意义。由于本实验设计的流体阻力装置目前只能控制粒子的流向，未来将该方法与微流控芯片结合，那么在降低流速的情况下，微粒的可控操纵距离会进一步提高。 参考文献 GRIER D G . A revolution in optical manipulation ［J］. Nature ， 2003 ， 424 （ 6950 ）： 810 - 816 . doi: 10.1038/nature01935 http://dx.doi.org/10.1038/nature01935 A revolution in optical manipulation CHIOU P Y ， OHTA A T ， WU M C . Massively parallel manipulation of single cells and microparticles using optical images ［J］. Nature ， 2005 ， 436 （ 7049 ）： 370 - 372 . doi: 10.1038/nature03831 http://dx.doi.org/10.1038/nature03831 Massively parallel manipulation of single cells and microparticles using optical images YANG A H J ， MOORE S D ， SCHMIDT B S ， et al . Optical manipulation of nanoparticles and biomolecules in sub-wavelength slot waveguides ［J］. Nature ， 2009 ， 457 （ 7225 ）： 71 - 75 . doi: 10.1038/nature07593 http://dx.doi.org/10.1038/nature07593 Optical manipulation of nanoparticles and biomolecules in sub-wavelength slot waveguides 赵明雪 ， 朱刘昊 ， 邵彤彤 ， 等 . 非对称动态涡旋光束对酵母菌细胞的光操纵特性 ［J］. 液晶与显示 ， 2021 ， 36 （ 6 ）： 841 - 847 . doi: 10.37188/CJLCD.2021-0023 http://dx.doi.org/10.37188/CJLCD.2021-0023 非对称动态涡旋光束对酵母菌细胞的光操纵特性 ZHAO M X ， ZHU L H ， SHAO T T ， et al . Optical manipulation characteristics of yeast cells by asymmetric dynamic vortex beams ［J］. Chinese Journal of Liquid Crystals and Displays ， 2021 ， 36 （ 6 ）： 841 - 847 . （in Chinese） . doi: 10.37188/CJLCD.2021-0023 http://dx.doi.org/10.37188/CJLCD.2021-0023 Optical manipulation characteristics of yeast cells by asymmetric dynamic vortex beams 林兴磊 ， 付文升 ， 邹永刚 ， 等 . 多局域空心光阱及其光场调控 ［J］. 光学 精密工程 ， 2021 ， 29 （ 2 ）： 251 - 258 . doi: 10.37188/OPE.20212902.0251 http://dx.doi.org/10.37188/OPE.20212902.0251 多局域空心光阱及其光场调控 LIN X L ， FU W SH ， ZOU Y G ， et al . Multi-bottle beam optical trap and its light field regulation ［J］. Opt. Precision Eng. ， 2021 ， 29 （ 2 ）： 251 - 258 . （in Chinese） . doi: 10.37188/OPE.20212902.0251 http://dx.doi.org/10.37188/OPE.20212902.0251 Multi-bottle beam optical trap and its light field regulation GUCK J ， ANANTHAKRISHNAN R ， MAHMOOD H ， et al . The optical stretcher： a novel laser tool to micromanipulate cells ［J］. Biophysical Journal ， 2001 ， 81 （ 2 ）： 767 - 784 . doi: 10.1016/s0006-3495(01)75740-2 http://dx.doi.org/10.1016/s0006-3495(01)75740-2 The optical stretcher： a novel laser tool to micromanipulate cells JESS P R T ， GARCÉS-CHÁVEZ V ， SMITH D ， et al . Dual beam fibre trap for Raman micro-spectroscopy of single cells ［J］. Optics Express ， 2006 ， 14 （ 12 ）： 5779 - 5791 . doi: 10.1364/oe.14.005779 http://dx.doi.org/10.1364/oe.14.005779 Dual beam fibre trap for Raman micro-spectroscopy of single cells LIU Y X ， YU M . Investigation of inclined dual-fiber optical tweezers for 3D manipulation and force sensing ［J］. Optics Express ， 2009 ， 17 （ 16 ）： 13624 - 13638 . doi: 10.1364/oe.17.013624 http://dx.doi.org/10.1364/oe.17.013624 Investigation of inclined dual-fiber optical tweezers for 3D manipulation and force sensing BRAMBILLA G ， MURUGAN G S ， WILKINSON J S ， et al . Optical manipulation of microspheres along a subwavelength optical wire ［J］. Optics Letters ， 2007 ， 32 （ 20 ）： 3041 - 3043 . doi: 10.1364/ol.32.003041 http://dx.doi.org/10.1364/ol.32.003041 Optical manipulation of microspheres along a subwavelength optical wire ZHANG Y ， LI B J . Particle sorting using a subwavelength optical fiber ［J］. Laser & Photonics Reviews ， 2013 ， 7 （ 2 ）： 289 - 296 . doi: 10.1002/lpor.201200087 http://dx.doi.org/10.1002/lpor.201200087 Particle sorting using a subwavelength optical fiber LIU Z H ， GUO C K ， YANG J ， et al . Tapered fiber optical tweezers for microscopic particle trapping： fabrication and application ［J］. Optics Express ， 2006 ， 14 （ 25 ）： 12510 - 12516 . doi: 10.1364/oe.14.012510 http://dx.doi.org/10.1364/oe.14.012510 Tapered fiber optical tweezers for microscopic particle trapping： fabrication and application MOHANTY S K ， MOHANTY K S ， BERNS M W . Manipulation of mammalian cells using a single-fiber optical microbeam ［J］. Journal of Biomedical Optics ， 2008 ， 13 （ 5 ）： 054049 . doi: 10.1117/1.2983663 http://dx.doi.org/10.1117/1.2983663 Manipulation of mammalian cells using a single-fiber optical microbeam YUAN L B ， LIU Z H ， YANG J ， et al . Twin-core fiber optical tweezers ［J］. Optics Express ， 2008 ， 16 （ 7 ）： 4559 - 4566 . doi: 10.1364/oe.16.004559 http://dx.doi.org/10.1364/oe.16.004559 Twin-core fiber optical tweezers GONG Y ， YE A Y ， WU Y ， et al . Graded-index fiber tip optical tweezers： numerical simulation and trapping experiment ［J］. Optics Express ， 2013 ， 21 （ 13 ）： 16181 - 16190 . doi: 10.1364/oe.21.016181 http://dx.doi.org/10.1364/oe.21.016181 Graded-index fiber tip optical tweezers： numerical simulation and trapping experiment XIN H B ， XU R ， LI B J . Optical trapping， driving and arrangement of particles using a tapered fibre probe ［J］. Scientific Reports ， 2012 ， 2 ： 818 . doi: 10.1038/srep00818 http://dx.doi.org/10.1038/srep00818 Optical trapping， driving and arrangement of particles using a tapered fibre probe LIU Z H ， WANG L ， LIANG P B ， et al . Mode division multiplexing technology for single-fiber optical trapping axial-position adjustment ［J］. Optics Letters ， 2013 ， 38 （ 14 ）： 2617 - 2620 . doi: 10.1364/ol.38.002617 http://dx.doi.org/10.1364/ol.38.002617 Mode division multiplexing technology for single-fiber optical trapping axial-position adjustment GONG Y ， HUANG W ， LIU Q F ， et al . Graded-index optical fiber tweezers with long manipulation length ［J］. Optics Express ， 2014 ， 22 （ 21 ）： 25267 - 25276 . doi: 10.1364/oe.22.025267 http://dx.doi.org/10.1364/oe.22.025267 Graded-index optical fiber tweezers with long manipulation length GONG Y ， ZHANG C L ， LIU Q F ， et al . Optofluidic tunable manipulation of microparticles by integrating graded-index fiber taper with a microcavity ［J］. Optics Express ， 2015 ， 23 （ 3 ）： 3762 - 3769 . doi: 10.1364/oe.23.003762 http://dx.doi.org/10.1364/oe.23.003762 Optofluidic tunable manipulation of microparticles by integrating graded-index fiber taper with a microcavity ZHANG C L ， GONG Y ， LIU Q F ， et al . Graded-index fiber enabled strain-controllable optofluidic manipulation ［J］. IEEE Photonics Technology Letters ， 2016 ， 28 （ 3 ）： 256 - 259 . doi: 10.1109/lpt.2015.2494583 http://dx.doi.org/10.1109/lpt.2015.2494583 Graded-index fiber enabled strain-controllable optofluidic manipulation LIU X S ， WU Y ， XU X H ， et al . Bidirectional transport of nanoparticles and cells with a bio-conveyor belt ［J］. Small ， 2019 ， 15 （ 50 ）： 1905209 . doi: 10.1002/smll.201905209 http://dx.doi.org/10.1002/smll.201905209 Bidirectional transport of nanoparticles and cells with a bio-conveyor belt MAMUTI R ， FUJI T K ， KUDO T . Opto-thermophoretic trapping of micro and nanoparticles with a 2 µm Tm-doped fiber laser ［J］. Optics Express ， 2021 ， 29 （ 23 ）： 38314 - 38323 . doi: 10.1364/oe.440866 http://dx.doi.org/10.1364/oe.440866 Opto-thermophoretic trapping of micro and nanoparticles with a 2 µm Tm-doped fiber laser PENG X L ， CHEN Z H ， KOLLIPARA P S ， et al . Opto-thermoelectric microswimmers ［J］. Light： Science & Applications ， 2020 ， 9 ： 141 . doi: 10.1038/s41377-020-00378-5 http://dx.doi.org/10.1038/s41377-020-00378-5 Opto-thermoelectric microswimmers 刘炳辉 ， 杨立军 ， 王扬 . 镀膜光纤探针近场捕获的模拟与实验 ［J］. 光学 精密工程 ， 2011 ， 19 （ 10 ）： 2355 - 2365 . doi: 10.3788/ope.20111910.2355 http://dx.doi.org/10.3788/ope.20111910.2355 镀膜光纤探针近场捕获的模拟与实验 LIU B H ， YANG L J ， WANG Y . Simulation and experiments of near-field trapping using metal-coated optical fiber probe ［J］. Opt. Precision Eng. ， 2011 ， 19 （ 10 ）： 2355 - 2365 . （in Chinese） . doi: 10.3788/ope.20111910.2355 http://dx.doi.org/10.3788/ope.20111910.2355 Simulation and experiments of near-field trapping using metal-coated optical fiber probe SCHMIDT B S ， YANG A H ， ERICKSON D ， et al . Optofluidic trapping and transport on solid core waveguides within a microfluidic device ［J］. Optics Express ， 2007 ， 15 （ 22 ）： 14322 - 14334 . doi: 10.1364/oe.15.014322 http://dx.doi.org/10.1364/oe.15.014322 Optofluidic trapping and transport on solid core waveguides within a microfluidic device 罗昊 . 基于微纳米结构表面的浸润性及液滴蒸发行为的研究 ［D］. 西安 ： 西北大学 ， 2016 . 基于微纳米结构表面的浸润性及液滴蒸发行为的研究 LUO H . Wetting and Evaporation Dynamics on Micro/nano Structured Surfaces ［D］. Xi'an ： Northwest University ， 2016 . （in Chinese） Wetting and Evaporation Dynamics on Micro/nano Structured Surfaces 李喜乐 ， 许文然 ， 张攀 . 表面张力法判断亲水性的研究 ［J］. 信息记录材料 ， 2018 ， 19 （ 5 ）： 17 - 20 . 表面张力法判断亲水性的研究 LI X L ， XU W R ， ZHANG P . Research on the determination of hydrophilicity by surface tension method ［J］. Information Recording Materials ， 2018 ， 19 （ 5 ）： 17 - 20 . （in Chinese） Research on the determination of hydrophilicity by surface tension method FAROOQ U . 基于微流控技术的声表面波片上实验室的研究 ［D］. 杭州 ： 浙江大学 ， 2019 . 基于微流控技术的声表面波片上实验室的研究 FAROOQ U . Study of Surface Acoustic wave Lab on Chip based Microfluidics ［D］. Hangzhou ： Zhejiang University ， 2019 . （in Chinese） Study of Surface Acoustic wave Lab on Chip based Microfluidics LIU X S ， HUANG J B ， LI Y C ， et al . Optofluidic organization and transport of cell chain ［J］. Journal of Biophotonics ， 2017 ， 10 （ 12 ）： 1627 - 1635 . doi: 10.1002/jbio.201600306 http://dx.doi.org/10.1002/jbio.201600306 Optofluidic organization and transport of cell chain 10.37188/OPE.20223013.1555 Q631 A 1004-924X（2022）13-1555-09 姜春雷,水华胜,陈朋等.基于光纤光镊的长距离微粒可控操纵[J].光学精密工程,2022,30(13):1555-1563. 姜春雷,水华胜,陈朋等.基于光纤光镊的长距离微粒可控操纵[J].光学精密工程,2022,30(13):1555-1563. JIANG Chunlei,SHUI Huasheng,CHEN Peng,et al.Controllable manipulation of long-distance microparticles based on fiber optic tweezers[J].Optics and Precision Engineering,2022,30(13):1555-1563. JIANG Chunlei,SHUI Huasheng,CHEN Peng,et al.Controllable manipulation of long-distance microparticles based on fiber optic tweezers[J].Optics and Precision Engineering,2022,30(13):1555-1563. 黑龙江省自然科学基金资助项目（No.LH2021F008） 2022-03-15 2022-04-12 2022-07-10 2022-08-07 光学精密工程 13 30