贝塞尔光镊颗粒物观测系统的设计研究

赵威伦 蔡宸 赵罡 赵春生†

北京大学物理学院大气与海洋科学系, 北京100871; †通信作者, E-mail: zcs@pku.edu.cn

摘要 搭建贝塞尔光镊颗粒物观测系统。通过圆锥透镜产生贝塞尔光, 利用贝塞尔光施加的光压力和反向气流施加的阻力, 稳定地悬浮气溶胶颗粒物。结合弹性光散射信号, 实现对气溶胶单颗粒粒径和折射率的测量。该系统还可用于研究不同环境条件下气溶胶的吸湿性、挥发性以及折射率等物理化学性质。关键词 贝塞尔光; 光镊; 气溶胶单颗粒

第56卷 第6期 2020年11月

上述设计不能在气相环境中捕获颗粒物。并且, 相比于对被捕获颗粒物进行观测, 更侧重对颗粒物进行微操作, 因此在气溶胶研究中的作用十分有限。贝塞尔光施加在颗粒物上的光压力无法让颗粒物受力平衡, 要在气相中稳定地悬浮颗粒物, 需要引入其他外力。2006年, Summers 等[25]改进 Arlt等[23]的设计, 在贝塞尔光传播方向引入反向的气体稳定层流。这种巧妙的设计既让颗粒物因同时受到方向相反的光压力和气流阻力而保持平衡, 又为颗粒物提供气相环境, 使贝塞尔光镊可以在气相中捕

[26]获颗粒物。此外, Carruthers 等 2010年发明利用两束反向的贝塞尔光来平衡颗粒物的方法, 也在气相中成功地悬浮颗粒物。

实现在气相中稳定悬浮颗粒物后, 就要对颗粒物进行测量。2009年, Meresman 等[27]利用 CCD相机(Charge Coupled Device Camera), 观测到颗粒物的散射光强分布, 借助米散射理论[28], 从散射光强分布反演颗粒物的粒径, 首次实现对贝塞尔光镊悬浮颗粒物的测量。2013年, Walkers 等[11]将衰荡腔与贝塞尔光镊结合, 实现对单颗粒消光截面的测量,并根据米散射理论[28], 从散射光强反演得到粒径。反演粒径时, 需要预先确定颗粒物的折射率。2015年之前, 对折射率的处理要么是用给定值, 要么是将折射率参数化成相对湿度的函数, 这些处理方法都需要已知颗粒物的化学组分。2015 年, Cotterell等[29]基于最大皮尔森相关的思想编写算法, 直接从散射光强同时得到粒径和折射率。实际环境中气溶胶颗粒物成分复杂, 这种不需要已知颗粒物的化学组分或相对湿度, 就能得到粒径和折射率的方法在气溶胶研究中具有现实意义。

目前贝塞尔光镊的应用场景比高斯光镊少, 主要应用在气溶胶热动力学(如吸湿性和挥发性[30)、

[29]气溶胶光学(如气溶胶消光效率 和折射率参数化 [31])方面, 应进一步拓展贝塞尔光镊的应用范围。贝塞尔光镊可以悬浮的颗粒物粒径接近积聚模态, 积聚模态气溶胶在实际大气中对环境及辐射都有重要影响, 因此利用贝塞尔光镊表征颗粒物具有实际意义。鉴于我国目前还没有贝塞尔光镊相关技术, 我们自主搭建贝塞尔光镊系统, 实现对气溶胶单颗粒的捕捉和测量。

1

北京大学学报(自然科学版)

贝塞尔光镊测量的原理

如图1(a)所示, 激光(Gaussian laser)沿 z轴正方1032

向传播, 垂直于z轴的平面为横截面。激光在横截面的光强分布为二维高斯分布, 中心强, 周围弱。激光通过圆锥透镜Axicon 后, 波矢方向改变的角度为 θ, 形成贝塞尔光[22]。圆锥透镜是一种特殊的透镜, 形态为圆锥体, 圆锥体的轴与z轴一致, 与z轴正交的平面在图1(a)中记为平面和平面′, 在图1(b)中记为平面12。贝塞尔光在zmax区域内可被观测到, 它的横截面由中心的亮核与一系列围绕中心亮核的同心圆环构成(图1(b)), 每个圆环上的光功率与中心亮核相同[32]。随着光束的传播, 贝塞尔光横截面的形态会发生变化, 在靠近圆锥透镜的某个位置为一个中心亮核和一个圆环(图2(a))。随着远离圆锥透镜, 圆环数目增加。当与圆锥透镜的距离足够远时, 形如图1(b)的结构消失(图2(t))。在一定的距离内, 中心亮核的形态不发生变化, 即中心亮核在一定的距离内保持“准无衍射”, 这段距离称为传播距离[33], 在图1(a)中标记为 zmax 和 z′max。当中心亮核内有颗粒物时, 中心亮核的光场将对颗粒物施加光压力[34], 在图 1(c)中用 Foptical 表示。当中心亮核的尺寸与颗粒物尺寸相当时, 光压力大到不可忽略。气溶胶颗粒尺寸为μm级, 可将凸透镜组中心亮核的尺寸缩小到μm级(如图1(a)所示)。凸透镜L1和L2相距 f1+f2, f1 和 f2分别是凸透镜L1和L2的焦距, f1>f2, 这样的透镜组配置称为望远镜。在z′max区域内将产生新的贝塞尔光, 中心亮核尺寸缩小 f1/f2 倍。

悬浮颗粒物, 要求其受力平衡。仅有贝塞尔光中心光核(Central core of Bessel beam)施加的光压力不能使颗粒物平衡, 需要引入其他外力。本文参考

[25] [34] Summers 等 和 Preston 等 的方案, 如图 1(c)所示, 颗粒物处于z′max 区域时, 受到贝塞尔光施加的正 z方向的光压力。引入沿负z方向流动的气体层流(Laminar gas flow), 对颗粒物施加负z方向的斯托克斯阻力, 在图1(c)中用 Fdrag表示。适当地调整光强和气流流速, 当 Foptical=fdrag 时, 颗粒物受力平衡, 被稳定地悬浮。

2贝塞尔光镊测量系统的搭建

本研究搭建的贝塞尔光镊测量系统如图3(a)所示, 图3(b)展示实景。采用 532 nm激光作为光源,激光依次经平面镜M1和 M2反射后, 照射在凸透镜L1上。平面镜M1和M2的作用为准直经过凸透镜L1和L2的激光, 通过调整平面镜M1和M2的位

赵威伦等 贝塞尔光镊颗粒物观测系统的设计研究

置或仰角, 使激光依次经过凸透镜L1 和 L2的中心。凸透镜L1和L2组成一望远镜, 焦距分别为f1= 75 mm, f2=150 mm, 相距 f1+f2=225 mm, 其作用是扩大激光横截面宽度w (扩大倍数为 f2/f1=2), 进而增加贝塞尔光的传播距离 zmax[35]:

z  w /2tan  ,  (  /2  /2)( n nm  1),其中, na和 nm分别是圆锥透镜和介质的折射率, 本文中介质为空气; α是圆锥透镜Axicon 的顶角, 本文中α = 178°。

激光经过扩束后, 依次经过平面镜M3和 M4反射, 照射在圆锥透镜 Axicon 上。平面镜M3 和M4的作用是准直经过Axicon 的激光, 使激光沿着圆锥透镜的轴传播。本文使用的圆锥透镜Axicon的顶角为178°。

北京大学学报(自然科学版)第56卷 第6期 2020年11月

激光经过圆锥透镜Axicon 后, 在 Axicon 与凸透镜L3之间形成贝塞尔光。贝塞尔光的尺寸可用

[35]贝塞尔光中心亮核的半径rb=2.4048/k sinθ 衡量,其中k为波矢。当贝塞尔光用来捕获颗粒物时, 要求 rb在微米级, 本文的圆锥透镜Axicon产生的贝塞尔光尺寸过大, 不能满足要求, 于是利用两个凸透镜L3和L4组成的望远镜来缩小rb。贝塞尔光经过平面镜M5和M6反射后, 进入凸透镜L3和L4组成的望远镜中。平面镜M5和M6的作用是使光场中心可以分别经过凸透镜L3和L4的中心。凸透镜L3和L4的焦距分别为f3=300 mm, f4=50.8 mm, 相距f3+f4=350.8 mm, rb被缩小为原来的f4/f3=1/5.9。新的贝塞尔光在凸透镜L4上方的样品池(Sample Cell)中形成, 用于捕获微粒。

平面镜M7的作用是改变光路, 使新生成的贝塞尔光竖直向上, 方便观测颗粒物。凸透镜L4与样品池之间有一透明载玻片, 用来防止气溶胶样品污染凸透镜L4。样品池是颗粒物被悬浮的地方, 是两端连通的中空结构, 用于引导并保持层流, 是捕获颗粒物的关键部位。为了保持层流, 国外的样品池常为长方体结构, 如1 cm × 1 cm × 10 cm [26]。本文的样品池为圆柱体结构, 长10 cm, 横截面直径为1 cm。与长方体相比, 光滑的圆柱体结构没有直角,更有利于层流的保持。

在样品池的另一端通入气流, 本文实验中为氮气(N2)。具有一定流量的氮气被分成两条支路, 分别经过流量控制器1 和 2(图 3(a)中 Flowmeter1和Flowmeter2)。经过 Flowmeter1的气流不经过水浸润, 始终保持干燥, 称为“干气”; 经过 Flowmeter2的气流经过水浸润, 含有一定量的水汽, 称为“湿气”。可以通过两个流量控制器分别控制“湿气”与“干气”的比例, 再将两路气流混合, 从而调节进入样品池中气流的相对湿度。

样品池的另一个上游是超声雾化器(图 3(a)中Nebulizer), 用来将样品溶液雾化, 产生一定量的多分散气溶胶颗粒。雾化产生的颗粒物顺着气路被氮气引入样品池, 在其中受到方向相反的两个力(贝塞尔光施加的竖直向上的光压力和气流施加的竖直向下的斯托克斯阻力)的作用, 颗粒物的重力可忽略。当颗粒物到达光压力与气流阻力平衡的位置时, 就被稳定地悬浮。对样品池环境的监测通过温湿探头实现。国外使用单一探头, 将其置于样品池内, 位于颗粒物的

1034下游[30]。使用单一探头可能导致测量偏差, 将其置于样品池内也可能对颗粒物所处的环境造成扰动。为了避免这些问题, 本研究在样品池外靠近样品池的上游和下游分别放置温‒湿探头(图 3(a)中 T/RH probe1和T/RH probe2), 颗粒物附近的温‒湿度可由两者插值得到, 探测结果比单一探头更可靠。

3气溶胶单颗粒的捕获及其散射函数的初步测量

当气溶胶颗粒物被稳定地悬浮后, 可以通过CCD相机对其进行观察。如图3(a)所示, CCD相机前应当加上足够的消光片, 防止因颗粒物散射光过强而损伤相机, 也可加上适当的透镜组来调整相机接收的视场。适当地调整CCD前的透镜组后, 可以通过CCD相机看到如图4(a)所示明暗相间的颗粒物散射光强度分布, 其中心左侧为气溶胶前向散射的光, 右侧为后向散射的光。贝塞尔光具有自重构的性质[24], 当其遇到一个颗粒物时, 在颗粒后部一段距离处又会形成新的贝塞尔光。这段距离约为Dpk/2kz, 其中Dp是颗粒物的直径, kz是波矢k在光轴上的投影。贝塞尔光的自重构性使得同时捕获多个颗粒物成为可能[21], 图4(b)显示利用本文搭建系统同时捕获的5个颗粒物组成的微粒链。这种微粒链不稳定, 一般最后只会剩下一个微粒。本文利用凸透镜组收集散射光。由于凸透镜的直径有限, 无法收集0°~180°的所有散射光, 收集到的散射角范围与凸透镜的直径以及凸透镜到被悬浮颗粒物之间的距离有关。如图3(a)所示, 被悬浮颗粒物在凸透镜上的投影应位于凸透镜的中心, 凸透镜的直径为d, 颗粒物与凸透镜的距离为x, 则收集到的散射光角度范围为 d /2  d /2 1 1   /2  tan  ,  x   x 每个角度对应的散射光强Iobs(ϕ)的相对大小可由图4(a)沿 y方向平均得到, 本文收集到的散射角度范围为约70°到约110°。经过y方向平均后, 图 4(a)中颗粒物的散射光强度 Iobs(ϕ)如图 5中五角星所示,散射由强变弱的部分为70°~74°, 79°~83°, 88°~93°和 99°~106°, 分别对应图5中区域 I, III, V和 VII;由弱变强的部分为 75°~78°, 84°~87°, 94°~98°和107°~110°, 分别对应图5中区域II, IV, VI和VIII。

根据前面得到的 Iobs(ϕ), 可以反演得到气溶胶 /2  tan

的颗粒物为硫酸铵与水的混合物, 硫酸铵的复折射率为 1.52, 水的复折射率为1.33, 因此被悬浮颗粒的折射率在1.33~1.52之间。基于 Carruthers 等[36]的工作, 可以假定被悬浮颗粒的粒径在1~4 μm之间。假定被悬浮颗粒物是均匀内混合的球形, 被颗粒物散射的光无偏振, 则根据米散射理论[28], 可由此计算不同粒径(Dp)、不同复折射率(m)的颗粒物在不同散射方向 ϕ的散射强度 Imie(dp, m, ϕ), 1 μm

在查算表 Imie(dp, m, ϕ)中找到与观测散射光强Iobs(ϕ)相近的最优米散射光强Imie(dp*, m*, ϕ)的思路如下所述。考虑到凸透镜的边缘对散射光强接收的影响, 在图5的8个区域中, 优先找到90°附近区域III, IV, V和VI中与Iobs(ϕ)趋势一致的 Imie(dp, m, ϕ),即 Imie(dp, m, ϕ)和 Iobs(ϕ)在区域 III, IV, V和VI都是先减小, 再增加, 再减小, 再增加。如果满足要求的(Dp, m)只有1组, 则这组(Dp, m)就定为被悬浮颗粒的粒径和折射率(D *, m*); 如果这时满足要求的

p (Dp, m)大于1 组, 则再增加区域限定条件; 如果没有满足要求的(Dp, m), 则逐个减小区域限定条件。本文案例中满足要求的(Dp, m)大于1组, 故增加区域限定条件, 直到满足区域限定条件的(Dp, m)只有1组为止。最终, 当区域限定条件增加到区域I, II, III, VI, V和 VI, 共6个区域时, 满足要求的(Dp, m)

1035

[16]

Journal of Quantitative Spectroscopy & Radiative Transfer, 2009, 110(14): 1293‒1306 Davis E J. A history of single aerosol particle levitation. Aerosol Science and Technology, 1997, 26(3): 212‒254 Al Zaitone B A, Tropea C. Evaporation of pure liquid droplets: comparison of droplet evaporation in an acoustic field versus glass-filament. Chemical Engineering Science, 2011, 66(17): 3914‒3921 Power R M, Reid J P. Probing the micro-rheological properties of aerosol particles using optical tweezers. Reports on Progress in Physics Physical Society, 2014, 77(7): 074601 Power R M, Simpson S H, Reid J P, et al. The transition from liquid to solid-like behaviour in ultrahigh viscosity aerosol particles. Chemical Science, 2013, 4(6): 2597‒2604 Mason B J, Walker J S, Reid J P, et al. Deviations from plane-wave mie scattering and precise retrieval of refractive index for a single spherical particle in an optical cavity. Journal of Physical Chemistry A, 2014, 118(11): 2083‒2088 Walker J S, Carruthers A E, Orr-ewing A J, et al. Measurements of light extinction by single aerosol particles. Journal of Physical Chemistry Letters, 2013, 4(10): 1748‒1752 Davies J F, Miles R E H, Haddrell A E, et al. Temperature dependence of the vapor pressure and evaporation coefficient of supercooled water. Journal of Geophysical Research Atmospheres, 2015, 119(18): 10931‒10940 Davies J F, Miles R E H, Haddrell A E, et al. Influence of organic films on the evaporation and condensation of water in aerosol. Proceedings of the National Academy of Sciences of the United States of America, 2013, 110(22): 8807‒8812 Cai C, Stewart D J, Preston T C, et al. A new approach to determine vapour pressures and hygroscopicities of aqueous aerosols containing semi-volatile organic compounds. Physical Chemistry Chemical Physics, 2014, 16(7): 3162‒3172 Lienhard D M, Bones D L, Zuend A, et al. Measurements of thermodynamic and optical properties of selected aqueous organic and organic–inorganic mixtures of atmospheric relevance. Journal of Physical Chemistry A, 2012, 116(40): 9954‒9968 Dennis-smither B J, Marshall F H, Miles R E H, et al.

[25]

Volatility and oxidative aging of aqueous maleic acid aerosol droplets and the dependence on relative humidity. Journal of Physical Chemistry A, 2014, 118 (30): 5680‒5691 Dennis-smither B J, Miles R E H, Reid J P. Oxidative aging of mixed oleic acid/sodium chloride aerosol particles. Journal of Geophysical Research Atmospheres, 2012, 117(D20): 1‒13 Lu J W, Rickards A M J, Walker J S, et al. Timescales of water transport in viscous aerosol: measurements on sub-micron particles and dependence on conditioning history. Physical Chemistry Chemical Physics, 2014, 16(21): 9819‒9830 Bones D L, Reid J P, Lienhard D M, et al. Comparing the mechanism of water condensation and evaporation in glassy aerosol. Proceedings of the National Academy of Sciences of the United States of America, 2012, 109(29): 11613‒11618 Haddrell A E, Davies J F, Miles R E H, et al. Dynamics of aerosol size during inhalation: hygroscopic growth of commercial nebulizer formulations. International Journal of Pharmaceutics, 2014, 463(1): 50‒61 Carruthers A E, Reid J P, Orr-ewing A J. Longitudinal optical trapping and sizing of aerosol droplets. Optics Express, 2010, 18(13): 14238‒14244 Mcgloin D, Dholakia K. Bessel beams: diffraction in a new light. Contemporary Physics, 2005, 46(1): 15‒ 28 Arlt J, Garces-chavez V, Sibbett W, et al. Optical micromanipulation using a Bessel light beam. Opt Commun, 2001, 197(4/5/6): 239‒245 Garces-chavez V, Mcgloin D, Melville H, et al. Simultaneous micromanipulation in multiple planes using a self-reconstructing light beam. Nature, 2002, 419: 145‒147 Summers M D, Reid J P, Mcgloin D. Optical guiding of aerosol droplets. Optics Express, 2006, 14(14): 6373‒6380 Carruthers A E, Orrewing A J, Reid J P. Optical trapping and sizing of aerosol droplets using counterpropagating Bessel beams. 2010, 7762(1): 77610V‒ 77620V