-
基于宽带光源微光学元件(如集成导光板)在衍射色散方面的设计需求, 本文构建了宽带光源微光学元件衍射理论分析模型, 探讨分析了衍射光谱的色度规律特性, 提出并定义了能准确定量衡量衍射光束色散程度的色散量C, 同时明确给出了零色散的边界判据点. 通过对研制的矩形位相光栅进行测试分析, 所得的光谱色度特性规律与理论分析结果相一致, 实验结果验证了本文提出的色散度判据参量C和零色散边界点的正确性. 本文提出的宽带光源色散度判据参量C、零色散边界判据点, 不仅为集成导光板结构参数的设计提供指导, 而且也能为其他宽带微光学元件的设计过程中探讨色散特性时提供指导.With the development of the microstructure fabrication process and the integration of micro-optical elements, diffractive micro-optical elements are widely used in broadband light sources, such as the integrated light guide plate (ILGP). And with the structural feature size of the ILGP decreasing from tens of microns to microns and even sub-microns, the diffraction dispersion phenomenon will inevitably become a prominent problem in research and design of non-dispersion elements. Nevertheless, under the broadband light source illumination, the analysis of the dispersion characteristic of diffraction spectrum of the microstructure array has not been reported in detail. Therefore a theoretical model of micro-optical element with a typical one-dimensional rectangular phase grating (RPG) and a widely used white LED source is established in this paper. The dispersion characteristic of the diffraction spectrum is studied, that is, with the increase of period of the RGP or the cone angle of incident beam, the dispersion of diffraction spectrum weakens. Dispersion parameter C and its formula are proposed, which can precisely measure the chromatic dispersion degree of the diffraction spectrum. Furthermore, the boundary criterion point of non-dispersion C = 0.3 is given explicitly. It is explored that no matter whether the cone angle of incident beam or the RPG period increases, the non-dispersion output light can be obtained only by matching the two parameters to make the dispersion parameter C less than 0.3. Then an RPG sample, of which the structural parameters are consistent with the designed ones, is fabricated by using micro-nano processing technology. By changing the cone angle of incident beam, the luminance and the chromaticity coordinates of the diffraction beam are tested. The analyses of the test results of the fabricated RPG sample show that the spectrum dispersion regularity is in accord with the theoretical analysis. The consistency verifies the correctness of dispersion parameter C, its formula and the non-dispersion boundary criterion point. The dispersion parameter C and non-dispersion boundary criterion point presented in this paper provide a guidance for analyzing the dispersion characteristics when the structural parameters of the integrated light guide plate and other broadband micro-optical element are designed.
1. 引 言
尘埃等离子体由电子、离子、中性粒子和尘埃颗粒组成, 在空间物理[1,2]、半导体工业领域[3]、聚变[4]以及实验室[5-9] 中受到广泛研究. 尘埃颗粒与等离子体中的电子和离子频繁碰撞而携带电荷, 通过颗粒间的相互作用形成强耦合系统. 在实验室研究中, 这一特性使尘埃等离子体成为模拟许多微观物理过程的理想平台, 如结晶[10-12]、相变[13-15]、孤波[16,17]和马赫锥[18,19]等. 尘埃颗粒通过粒子-粒子以及粒子-波等相互作用形式从等离子体中获取能量, 这些尘埃颗粒的动能可以高达上百个eV[20,21]. 然而, 通常情况下尘埃颗粒悬浮在极板鞘层附近进行随机运动, 其携带的能量不能被有效地收集并利用. 因此, 如何将尘埃颗粒的随机运动整流成定向运动成为解决这一问题的关键.
费曼棘轮可以将粒子的随机运动转化为定向运动, 且不违反热力学第二定律[22,23]. 费曼棘轮通常由具有一定不对称性的周期性棘轮单元组成, 由于费曼棘轮原理的普适性与重要性, 关于它的研究已经延伸到量子棘轮[24,25]、有机电子棘轮[26]以及光学棘轮[27]等学科领域. 这些系统中具有的某种对称性破缺机制是费曼棘轮运行的关键. 这些棘轮实验主要集中在微纳尺度的棘轮结构、超导下或者液相环境下, 实现条件普遍比较苛刻且实验难度较大. 近期, 研究者们利用尘埃等离子体系统易于实验和便于观察等特点, 将费曼棘轮原理引入到尘埃等离子体物理领域, 设计出尘埃等离子体圆棘轮, 实现了对尘埃颗粒的可控性整流, 即将尘埃颗粒的随机运动整流为定向运动[28]. 系统中的不对称性来源于水平极板上两个圆形树脂棘轮围成的棘轮通道, 带电尘埃颗粒在棘轮通道中受到棘轮势的作用, 形成的棘轮效应使得尘埃颗粒可以沿棘轮通道定向输运. 由于树脂棘轮位于极板上的鞘层内, 其形成的是悬浮电位, 在实验中很难精确控制其电位进而有效调控棘轮势与颗粒输运. 如果使用金属棘轮并将其置于极板上, 金属棘轮可以与极板同电位, 通过调节放电条件可准确控制金属棘轮电位, 并且可使其具有更低的电位水平, 由此可实现对棘轮势在较大范围内的精确调控, 此时, 尘埃颗粒输运可能会有新现象, 这方面的实验研究还未见报道. 本文通过设计尘埃等离子体金属直棘轮实验装置, 研究了微米级单分散尘埃颗粒在不同放电条件下的整流规律以及不同尺寸尘埃颗粒在相同放电条件下的整流现象, 进一步通过等离子体流体模拟获得等离子体参量的二维分布, 并基于Langevin方程对不同尺寸颗粒的运动进行数值模拟, 在相同条件下分析不同尺寸的尘埃颗粒在尘埃等离子体金属直棘轮中的整流机制.
2. 实验装置
图1为尘埃等离子体金属直棘轮实验装置示意图. 在尺寸为30 cm×20 cm×15 cm的真空室内水平放置两个极板, 上极板氧化铟锡(ITO)玻璃接地, 下极板不锈钢板接13.56 MHz的射频电源和匹配器. 下极板的上表面由30个宽为5 mm和高为3 mm的锯齿槽构成(如图1所示依次编号), 其沿x方向具有不对称性和周期性变化特点. 同时, 在下极板上面沿x方向平行设置两个金属条, 其下表面具有与下极板相同尺寸的锯齿槽, 并与下极板完全啮合. 两金属条之间距离为10 mm, 因此, 两个金属条和下极板构成沿x方向的棘轮通道. 利用容性耦合氩气辉光放电产生均匀的等离子体后, 棘轮通道表面形成鞘层, 其沿x方向具有不对称性与周期性分布特征. 利用颗粒注射器在棘轮通道中间“0”号锯齿槽处持续不断地缓慢投入微米大小的聚苯乙烯微球(尘埃颗粒). 这些尘埃颗粒在等离子体中自由下落时快速充电, 在受到下极板鞘层电场的作用下一般悬浮于棘轮通道上方4—8 mm处, 在一定的实验条件下能够沿棘轮通道形成定向输运, 当棘轮通道内所有锯齿槽均被尘埃颗粒填满时, 尘埃颗粒会流出棘轮通道并流入棘轮通道两端外侧的收集槽中. 通过调节放电条件或使用不同大小的尘埃颗粒, 可以控制尘埃颗粒的输运方向与速度. 为了方便观测, 在真空室外两侧相向沿棘轮通道方向施加两束5 mW的激光束(532 nm)用于照明尘埃颗粒, 尘埃颗粒的运动状态利用位于真空室上方的摄像机进行记录. 为了研究金属直棘轮中尘埃颗粒的整流机制, 分别使用直径为10, 15和 24 µm的单分散聚苯乙烯微球进行实验, 气压p变化范围是20—50 Pa, 功率P变化范围是10—50 W.
图 1 尘埃等离子体金属直棘轮实验示意图. 下极板的上表面由30个锯齿槽构成(图中只显示7个), 结合两个金属条, 形成沿x方向的棘轮通道. 单分散尘埃颗粒由棘轮通道中部的“0”号锯齿槽投入, 在特定条件沿棘轮通道会产生定向运动Fig. 1. Schematic diagram of the dusty plasma metal straight ratchet. The upper surface of the lower electrode is composed of 30 sawteeth (shown only 7 sawteeth), combined with two metal strips, a ratchet channel along the x direction is formed. Monodispersed dust particles are introduced into the “0” sawtooth at the middle of the ratchet channel, which exhibit directional flow along the ratchet channel under specific conditions.3. 实验结果
利用尘埃等离子体金属直棘轮能够实现尘埃颗粒的整流. 为了展示整流现象, 实验中从棘轮通道中间“0”号锯齿槽处持续地撒入直径为15 µm的单分散尘埃颗粒, 并观察尘埃颗粒的运动情况. 尘埃颗粒在等离子体中自由下落时快速充电, 在落入棘轮通道后受下极板及两侧金属条鞘层电场作用悬浮并被约束在棘轮通道的中间位置. 当单个锯齿槽中只有少量尘埃颗粒时, 尘埃颗粒在锯齿槽中呈现单链分布, 如果继续加入尘埃颗粒, 由于颗粒间Yukawa静电排斥作用, 尘埃颗粒在单个锯齿槽内形成中间双链两端单链的分布, 此时单个锯齿槽中的尘埃颗粒达到饱和, 约有nmax=17个, 此过程如图2(a), (c)和(e)中局部放大插图所示. 这个饱和颗粒数nmax的大小与尘埃颗粒半径和放电条件有一定的关系, 但是影响不大, 只有2—3个的差别. 此时, 再继续加入尘埃颗粒, 在颗粒间相互作用力和带电颗粒之间的集体效应的推动下, 尘埃颗粒分布由“0”号锯齿槽向其相邻的锯齿槽扩展, 如图2(a), (b)所示, 此时尘埃颗粒并没有形成持续的单方向输运(整流). 如果持续不断地撒入尘埃颗粒, 棘轮的整流效应开始显现出来, 在非平衡离子拖拽力的作用下[28], 尘埃颗粒由“0”号锯齿槽依次流入“1”, “2”, “3”, ⋯号锯齿槽进而形成持续的正向输运, 即正流, 如图2(c)—(f)所示. 特别要注意的是, “–1”号锯齿槽在初始撒颗粒时会涌入少量颗粒, 在尘埃颗粒形成持续的整流后, “–1”号锯齿槽内颗粒数逐渐饱和, 其左侧的“–2”号等锯齿槽内的颗粒数为 0, 如图2(e)所示. 图2(b), (d), (f)分别给出了随着颗粒数的增加颗粒数N在棘轮通道内的分布与演化过程. 此外, 当停止向棘轮通道中撒入尘埃颗粒时, 棘轮通道内的尘埃颗粒将在各自的平衡位置上趋于静止, 这是由于尘埃颗粒的动能小于沿棘轮通道x方向的棘轮势能, 尘埃颗粒被束缚在棘轮势阱中, 此时尘埃颗粒不能形成沿棘轮通道持续的输运, 只有持续不断地向棘轮通道中撒入尘埃颗粒时, 棘轮通道内尘埃颗粒的平衡状态才会被打破, 尘埃颗粒才会沿着棘轮通道形成持续的正流. 尘埃颗粒的定向输运速度与往“0”号锯齿槽中投入颗粒数的速率有关, 在本次实验中, 尘埃颗粒形成稳定正流后的平均速度约为0.53 mm/s. 因此, 当放电气压p=35 Pa, 放电功率P=40 W时, 尘埃等离子体金属直棘轮中颗粒形成正流.
图 2 尘埃等离子体金属直棘轮中尘埃颗粒的正流 (a), (c), (e)随着投入尘埃颗粒数的增加, 尘埃颗粒在棘轮通道中的空间分布逐渐向右扩展, 形成正向整流, 图中N代表棘轮通道中尘埃颗粒总数量; (b), (d), (f)对应的不同锯齿槽位置处颗粒数n的分布情况, 放电气压为 , 放电功率为p=35 Pa , 颗粒直径为15 µm. 图中尘埃颗粒做了亮度增强处理P=40 W Fig. 2. Positive flow of dust particles in the dusty plasma metal straight ratchet: (a), (c), (e) Number of introduced dust particles increases, the spatial distribution of dust particles in the ratchet channel gradually expands to the right, forming a positive flow, N represents the total number of dust particles in the ratchet channel; (b), (d), (f) the corresponding distributions of particle number n along the ratchet channel, gas pressure , discharge powerp=35 Pa , particle diameter is 15 µm. Brightness of dust particles in images is enhanced.P=40 W 通过改变放电功率可以改变尘埃颗粒的整流速度. 当保持放电气压p=35 Pa不变, 逐渐降低放电功率会导致尘埃颗粒整流速度的减缓, 当放电功率降至P=25 W时, 尘埃颗粒将无法形成定向输运. 在新的一次实验中(P=25 W), 从棘轮通道中间“0”号锯齿槽处缓慢地撒入尘埃颗粒. 刚开始时, 由于撒入的颗粒数量较少且受到棘轮电势的约束, 尘埃颗粒会积聚在“0”号锯齿槽内. 随着撒入颗粒数N的增加, 由于带电尘埃颗粒之间的排斥作用, “0”号锯齿槽不能再容纳下这些尘埃颗粒, 它们将在集体效应的作用下以近似相等的速度均匀地向“0”号锯齿槽两侧流动, 直到这些尘埃颗粒均匀分布并铺满整个棘轮通道, 如图3(a), (c), (e)所示. 从图3(b), (d), (f)所示的尘埃颗粒在棘轮通道中的分布曲线可以看出, 随着尘埃颗粒数量的增加, 颗粒的分布曲线逐渐向棘轮通道两端同时上扬, 并且在“0”号锯齿槽(虚线位置)附近呈现出明显的对称性, 这表明尘埃颗粒沿棘轮通道输运的净位移为零. 此时, 尽管尘埃颗粒表现出双向输运现象, 但这一现象并非源自棘轮势的影响, 而只是由颗粒之间的相互排斥作用引起的.
图 3 尘埃等离子体金属直棘轮中尘埃颗粒的无净位移运动 (a), (c), (e)随着投入尘埃颗粒数的增加, 尘埃颗粒在棘轮通道中的空间分布逐渐向两侧同时扩展, 图中N代表棘轮通道中尘埃颗粒总数量; (b), (d), (f)对应的不同锯齿槽位置处颗粒数n的分布情况. 放电气压为 , 放电功率为p=35 Pa , 颗粒直径为15 µm, 图中尘埃颗粒进行亮度增强处理P=25 W Fig. 3. None net flow of dust particles in the dusty plasma metal straight ratchet: (a), (c), (e) Number of introduced dust particles increases, the spatial distribution of dust particles in the ratchet channel gradually expands to both sides simultaneously, N represents the total number of dust particles in the ratchet channel; (b), (d), (f) the corresponding distributions of the particle number n along the ratchet channel, gas pressure , discharge powerp=35 Pa , particle diameter is 15 µm, brightness of dust particles in images is enhanced.P=25 W 通过改变放电功率还可以进一步控制尘埃颗粒的整流方向. 在新的一次实验中, 保持放电气压p = 35 Pa不变, 将放电功率降到P=15 W时, 从棘轮通道中间“0”号锯齿槽处逐渐撒入尘埃颗粒. 随着撒入颗粒数量的增加, 颗粒会从“0”号锯齿槽依次流入“–1”, “–2”, “–3”, ⋯号锯齿槽进而形成持续的逆向输运, 即逆流, 如图4(a), (c)和(e)所示. 在经过约50 s后, 尘埃颗粒在棘轮通道中形成稳定的逆流. 此时棘轮通道中尘埃颗粒的分布不会再随时间的推移而发生显著变化. 在此过程中, 尘埃颗粒的数量分布曲线向棘齿通道左侧逐渐上扬直到饱和, 如图4(b), (d), (f)所示. 此时, 如果停止投入尘埃颗粒, 这种持续的尘埃颗粒整流过程将中断, 这与图2所示的尘埃颗粒正流过程类似. 在本次实验中, 尘埃颗粒沿棘轮通道形成稳定逆流后的平均速度约为–0.12 mm/s. 通过对比图2—图4发现, 在保持放电气压p=35 Pa不变的情况下, 随着放电功率的逐渐降低, 尘埃颗粒沿棘轮通道的整流速度表现出由正向逐渐降低至零并进一步转为逆向, 发生了流的反转现象.
图 4 尘埃等离子体金属直棘轮中尘埃颗粒的逆向整流 (a), (c), (e)随着投入尘埃颗粒数的增加, 尘埃颗粒在棘轮通道中的空间分布逐渐向左扩展, 形成逆向整流. 图中N代表棘轮通道中尘埃颗粒总数量. (b), (d), (f)对应的不同锯齿槽位置处颗粒数n的分布情况. 放电气压为 , 放电功率为p=35 Pa , 颗粒直径为15 µm. 图中尘埃颗粒做了亮度增强处理P=15 W Fig. 4. Negative flow of dust particles in the dusty plasma metal straight ratchet: (a), (c), (e) Number of introduced dust particles increases, the spatial distribution of dust particles in the ratchet channel gradually expands to the left, forming a negative flow, N represents the total number of dust particles in the ratchet channel; (b), (d), (f) the corresponding distributions of the particle number n along the ratchet channel. Gas pressure , discharge powerp=35 Pa , particle diameter is 15 µm. Brightness of dust particles in images is enhanced.P=15 W 尘埃颗粒整流及其反转的现象不仅与放电功率的变化有关, 还与放电气压密切相关. 通过对放电气压和放电功率进行较大范围的参数扫描, 将其对尘埃颗粒整流行为的影响绘制成相图, 如图5所示. 在该相图中, 黑色曲线为尘埃颗粒呈现正、逆流的临界线, 即在这条黑色临界线上对应的放电功率与气压条件下尘埃颗粒沿棘轮通道的净位移为零(没有定向输运), 在临界线上方区域的实验条件下尘埃颗粒形成正流, 反之形成逆流. 这表明, 尘埃颗粒的整流可以在一个较为宽泛的放电功率和气压条件下实现, 并且, 在较高气压和较高放电功率的条件下尘埃颗粒更容易发生正流现象, 而在较低气压和较低放电功率的条件下尘埃颗粒则更容易发生逆流现象. 因此, 图5表明, 通过调节放电功率和气压能够精确调控尘埃颗粒的整流方向. 此外, 尘埃颗粒的投入速度对尘埃颗粒整流速度的大小有一定的影响, 由于实验中很难实现匀速定量的投入尘埃颗粒, 导致尘埃颗粒的整流速度不太稳定, 所以目前很难在相图中同时给出整流速度大小与放电功率和气压的关系.
图 5 不同放电气压和功率实验条件下尘埃颗粒整流的相图. 黑色曲线表示尘埃颗粒在棘轮通道中呈现正流和逆流的临界线, 颗粒直径为15 µmFig. 5. Phase diagram of the directional flow of dust particles under different gas pressure and power experimental conditions. The black curve represents the critical line separating the positive and negative flows of dust particles in the ratchet channel, particle diameter is 15 µm.以上实验结果展示了直径为15 µm的单分散尘埃颗粒在不同放电条件下的整流现象. 由于尘埃颗粒在棘轮通道中的悬浮高度与其直径密切相关, 因此, 不同直径的尘埃颗粒的悬浮高度不同, 所处的等离子体环境不同, 进而其整流行为也不相同. 通过对不同直径(10, 15和24 µm)的尘埃颗粒在相同装置中的整流行为进行系统研究, 结果表明这些尺寸的尘埃颗粒在棘轮通道中均可实现整流, 但其正流和逆流发生转变的临界放电条件存在一定差异, 这些差异并不影响尘埃颗粒定向输运的定性结论. 图6给出了在放电气压p=35 Pa, 放电功率P=25 W的实验条件下, 这3种尺寸尘埃颗粒在棘轮通道中的整流情况. 实验结果表明, 在这一放电条件下, 直径为10 µm和24 µm的尘埃颗粒在棘轮通道中分别形成正流(平均速度约0.39 mm/s)和逆流(平均速度约–0.11 mm/s), 分别如图6(a), (c)所示, 而直径为15 µm的尘埃颗粒在棘轮通道中没有形成定向输运(净位移为零, 图6(b)). 通过测量进一步发现, 随着尘埃颗粒直径的增大, 其在棘轮通道中的悬浮高度分别由z = 5.9 mm降低到4.7 mm和4.1 mm, 因此, 尘埃颗粒的悬浮高度影响了其输运结果. 此外, 在图5所示相图的放电参数范围内, 直径在 5—35 µm范围内的尘埃颗粒都可以进行整流.
图 6 同一放电条件下不同尺寸尘埃颗粒在棘轮通道内的运动情况 (a)直径为10 µm的尘埃颗粒随着投入颗粒数目的增加逐渐形成正流; (b)直径为15 µm的尘埃颗粒在棘轮通道中没有形成定向输运; (c)直径为24 µm的尘埃颗粒随着撒入颗粒数目的增加逐渐形成逆流, 放电气压 , 功率p=35 Pa . 图中尘埃颗粒做了亮度增强处理P=25 W Fig. 6. The flow of dust particles with different sizes in the ratchet channel under same discharge conditions: (a) Dust particles with a diameter of 10 µm gradually form a positive flow as the number of introduced dust particles increases; (b) dust particles with a diameter of 15 µm do not form a directional flow in the ratchet channel; (c) dust particles with a diameter of 24 µm gradually form a negative flow as the number of introduced particles increases, gas pressure , discharge powerp=35 Pa . Brightness of dust particles in images is enhanced.P=25 W 4. 数值模拟
尘埃颗粒在棘轮通道中的定向输运与其周围的鞘层等离子体参量密切相关. 为了揭示金属直棘轮中尘埃颗粒整流机制, 首先要掌握等离子体参量的分布情况. 由于棘轮通道的空间尺度非常小, 常用的朗缪尔探针不能准确测量此种情况下的电子密度等参量. 因此, 本文基于尘埃等离子体棘轮模型[29,30]通过数值模拟的方法来获得棘轮通道内的等离子体参量. 这里给出主要的模型计算流程.
由于实验中不同尺寸的尘埃颗粒悬浮在不同高度上(z方向), 并沿棘轮通道x方向进行整流, 因此选择沿棘轮通道中心(y = 0处)的x-z平面进行二维数值模拟. 首先, 利用COMSOL软件建立金属直棘轮射频辉光放电流体模型, 以获得下极板锯齿槽表面鞘层的等离子体参量, 由于金属直棘轮呈现周期性排列结构, 模拟中只选取一个锯齿单元作为模拟区域, 并在x方向上采用周期性边界条件, 放电条件为气压p=35 Pa, 功率P=25 W, 其余各项参数均与实验中保持一致. 随后再利用Matlab软件进一步分析尘埃颗粒在鞘层中的动力学过程. 尘埃颗粒在鞘层中受到的主要作用力包括: 重力 Fg、离子拖拽力Fi、中性粒子阻力Fn、颗粒之间的Yukawa作用力FY、鞘层电场力Fe和涨落力Fη, 其Langevin运动方程为
m¨rk=Fg+Fi+Fn+FY+Fe+Fη,k=1,2,⋯,n, (1) 其中, n为颗粒数, m为颗粒质量, 详细的受力公式参见文献[28, 29]. 根据上述方程, 进一步模拟实验中不同尺寸尘埃颗粒在棘轮通道中的整流过程. 由于尘埃颗粒整流实验中使用的气压较高, 带电颗粒间相互作用势能很强, 处于强耦合状态, 其热运动非常弱, 因此在动力学模拟中将涨落力忽略. 锯齿槽中尘埃颗粒的数量为19个, 沿水平方向上等间隔分布, 初始高度z0 = 7 mm, 初始速度为零. 根据模拟得到的等离子体参量, 利用OML理论[31]可以计算出尘埃颗粒的带电量. 以直径为15 µm的尘埃颗粒为例, 其带电量Q约为−105e, 在棘轮通道中的整流过程中, 其带电量基本保持不变.
下极板锯齿槽表面鞘层中垂直方向电场Ez = –∂U/∂z的等值线表现出明显的不对称性和周期性特征, 如图7(a)所示, 这是由锯齿形电极结构所呈现的不对称与周期性变化决定的. 尘埃颗粒在鞘层中的悬浮高度主要依靠重力、鞘层电场力和竖直z方向上的离子拖拽力Fiz的共同作用, 且满足mg + Fiz=QEz, 其中Fiz约为尘埃颗粒重力的3%, 可以忽略. 由于尘埃颗粒在棘轮通道内水平输运时其带电量Q基本不变, 因此, 尘埃颗粒的悬浮高度由Ez决定, 即同一尺寸的单分散尘埃颗粒将沿其对应的电场等值线悬浮在鞘层中. 对于实验中使用的3种尺寸的尘埃颗粒, 它们分别沿各自对应的Ez等值线悬浮在鞘层中, 其平均悬浮高度分别为5.4, 4.3和3.9 mm, 较小的颗粒悬浮高度较高, 较大的颗粒悬浮高度较低, 如图7(a)所示, 这与实验观测结果一致.
图 7 两个锯齿槽上表面的鞘层电场 的等值线分布(a)以及在这些等值线上对应的电势变化(b). 图(a)中鞘层电场的竖直分量Ez 的等值线呈现不对称性分布, 不同尺寸的尘埃颗粒沿不同的Ez 的等值线悬浮于鞘层中; 图(b)中的电势分布具有棘轮势特征, 因此, 不同尺寸的尘埃颗粒悬浮于鞘层不同的高度, 受到的棘轮势具有不同的特征, 进而使得其形成沿棘轮通道不一样的输运结果. 图中红、绿、蓝球分别表示直径为10, 15和24 µm的尘埃颗粒, 放电气压Ez , 功率p=35 Pa P=25 W Fig. 7. The contour distribution of the sheath electric field on the surface of two sawteeth (a) and the corresponding potential variation along these contour curves (b). (a) The contour curves of the vertical component of the sheath electric fieldEz show asymmetric distributions, dust particles of different sizes are suspended along different contour curves within the sheath; (b) the potential distribution has the characteristic of a ratchet potential. Therefore, dust particles of different sizes are suspended at different heights within the sheath, experiencing ratchet potentials with different characteristics, which leads to distinct transport outcomes along the ratchet channel. The red, green and blue balls in the curves indicate dust particles with diameter of 10, 15 and 24 µm, respectively. Gas pressureEz , discharge powerp=35 Pa .P=25 W 鞘层中这些不同Ez等值线上的电势分布具有其各自的分布特点, 如图7(b)所示. 这些电势线的分布具有一定的不对称性与周期性特征, 是典型的棘轮势, 可通过双正弦函数来反映[22]:
U(x)=ΔU[sin(2πx/L)−14sin(4πx/L)]+U0, (2) 式中, L代表一个锯齿单元的长度, ΔU代表棘轮势振幅, U0代表棘轮势偏置, ΔU和U0分别为
ΔU=(Ul−Uh)/2, (3) U0=(Ul+Uh)/2, (4) 其中, Uh为锯齿尖处的电势, Ul为锯齿槽处的电势. ΔU的符号与大小分别决定了棘轮势的不对称取向与势阱高度.
棘轮势决定了尘埃颗粒在棘轮通道中的输运过程. 这是因为棘轮势产生的电场Ex=−∂U/∂x会引起x方向的离子流, 从而对尘埃颗粒施加非平衡离子拖拽力Fix. 棘轮势的不对称性使非平衡的离子拖拽力在单个棘轮势中的积分fix=∫L0Fixdx可不为零, 其量级约为10−13 N, 足以克服同等量级的中性气体阻力[5,32], 进而推动尘埃颗粒形成定向输运. 棘轮势的不对称性取向决定了尘埃颗粒在棘轮通道中的输运方向. 当尘埃颗粒悬浮在鞘层较高位置时, Ul>Uh, 棘轮势具有正取向(ΔU>0), 尘埃颗粒呈现正流, 如图7(b); 当尘埃颗粒处于较低位置处, Ul<Uh, 棘轮势具有负取向(ΔU<0), 尘埃颗粒呈现逆流; 特别是, 在某一中等高度处, 电势呈水平直线分布(ΔU=0), 不对称性消失, 尘埃颗粒不能形成定向输运(净位移为零). 图7数值计算结果表明, 在相同的放电条件下, 不同尺寸的尘埃颗粒悬浮在鞘层不同高度上, 受到的棘轮势取向不同, 进而形成不同的定向输运, 直径10 µm的尘埃颗粒正流速度平均约0.37 mm/s, 直径24 µm的尘埃颗粒逆流速度平均约–0.13 mm/s, 这与图6实验结果基本一致. 此外, 由于尘埃颗粒在鞘层中的悬浮高度还与放电条件有关, 因此, 还可以通过调节放电气压与功率来控制单分散尘埃颗粒的悬浮高度, 进而控制其输运方向, 即图2—图5所示实验结果.
5. 结 论
本文首先通过设计尘埃等离子体金属直棘轮实验装置, 实现了在较大放电参数范围内对尘埃颗粒的可控性整流. 其次, 在相同的放电条件下对不同尺寸的尘埃颗粒进行整流实验, 观察到较小的尘埃颗粒呈现正流, 较大的尘埃颗粒呈现逆流. 最后, 通过对尘埃等离子体金属直棘轮进行二维数值模拟, 研究不同尺寸尘埃颗粒在相同实验条件下的整流机制. 研究结果表明: 不同尺寸的尘埃颗粒在棘轮通道中的悬浮高度存在差异, 小颗粒悬浮在较高位置, 而大颗粒悬浮在较低位置. 由于这种悬浮高度的差异, 不同尺寸尘埃颗粒在其各自悬浮高度上受到的棘轮势的不对称取向相反, 这种棘轮势取向差异直接决定了不同尺寸的尘埃颗粒沿棘轮通道朝相反方向流动. 本文研究结果从根本上阐明了尘埃等离子体金属直棘轮中尘埃颗粒整流的机制, 为进一步研究双分散尘埃颗粒分离奠定了理论基础, 同时促进了等离子体物理和非平衡统计物理的交叉融合, 对解决芯片制造工艺中的尘埃污染问题也具有重要意义.
[1] Xu P, Huang H X, Wang K, Ruan S C, Yang J, Wan L L, Chen X X, Liu J Y 2007 Opt. Express 15 809
Google Scholar
[2] Xu P, Hong C Q, Cheng G X, Zhou L, Sun Z L 2015 Opt. Express 23 6773
Google Scholar
[3] Huang H X, Ruan S C, Yang T, Xu P 2015 Nano-Micro Lett. 7 177
Google Scholar
[4] 黄海漩, 徐平, 阮双琛, 杨拓, 袁霞, 黄燕燕 2015 物理学报 64 154212
Google Scholar
Huang H X, Xu P, Ruan S C, Yang T, Yuan X, Huang Y Y 2015 Acta Phys. Sin. 64 154212
Google Scholar
[5] 徐平, 袁霞, 杨拓, 黄海漩, 唐少拓, 黄燕燕, 肖钰斐, 彭文达 2017 物理学报 66 124201
Google Scholar
Xu P, Yuan X, Yang T, Huang H X, Tang S T, Huang Y Y, Xiao Y F, Peng W D 2017 Acta Phys. Sin. 66 124201
Google Scholar
[6] 徐平, 唐少拓, 袁霞, 黄海漩, 杨拓, 罗统政, 喻珺 2018 物理学报 67 024202
Google Scholar
Xu P, Tang S T, Yuan X, Huang H X, Yang T, Luo T Z, Yu J 2018 Acta Phys. Sin. 67 024202
Google Scholar
[7] Thomschke M, Reineke S, Lüssem B, Leo K 2012 Nano Lett. 12 424
Google Scholar
[8] Siitonen S, Laakkonen P, Vahimaa P, Kuittinen M, Tossavainen N 2006 Appl. Opt. 45 2623
Google Scholar
[9] Zhao X N, Hu J P, Lin Y, Xu F, Zhu X J, Pu D L, Chen L S, Wang C H 2016 Sci. Rep. 6 28319
Google Scholar
[10] Xue C X, Cui Q F 2010 Opt. Lett. 35 986
Google Scholar
[11] Tsukamoto H, Nishiyama M 2006 Jpn. J. Appl. Phys. 45 6678
Google Scholar
[12] Xu P, Huang Y Y, Zhang X L, Huang J F, Li B B, Ye E, Duan S F, Su Z J 2013 Opt. Express 21 20159
Google Scholar
[13] Xu P, Huang Y Y, Su Z J, Zhang X L, Luo T Z, Peng W D 2015 Opt. Express 23 4887
Google Scholar
[14] Xu P, Yan Z L, Wan L L, Huang H X 2004 Proceedings of SPIE Holography Diffractive Optics and Applications Ⅱ Beijing, China, November 8−11, 2004 p66
[15] Park S R, Kwon O J, Shin D, Song S H, Lee H S, Choi H Y 2007 Opt. Express 15 2888
Google Scholar
[16] Yang X P, Yan Y B, Jin G F 2005 Opt. Express 13 8349
Google Scholar
[17] Caputo R, Sio L D, Jak M J J, Hornix E J, Boer D K G, Cornelissen H J 2007 Opt. Express 15 10540
Google Scholar
[18] Xu M, Urbach H P, Boer D K G 2007 Opt. Express 15 5789
Google Scholar
[19] 张以谟 2008 应用光学 (北京:电子工业出版社) 第7章
Zhang Y M 2008 Applied Optics (Beijing: Publishing House of Electronics Industry) Chapter 7 (in Chinese)
[20] 吕乃光 2006 傅里叶光学 (北京:机械工业出版社)第70−113页
Lü N G 2006 Fourier Optics (Beijing: China Machine Press) pp70−113 (in Chinese)
[21] 苏显渝, 李继陶 1999 信息光学 (北京:科学出版社) 第44页
Su X Y, Li J T 1999 Information Optics (Beijing: Science Press) p44 (in Chinese)
-
图 2 零级、± 1级三衍射级次的相对光强分布 (a)−(e)入射光束锥角θ和光栅周期d分别分别为(7.5°, 4 μm), (7.5°, 8 μm), (7.5°, 40 μm), (14.5°, 4 μm), (102.5°, 4 μm); 图中红色、绿色、蓝色线分别代表三原色的红光、绿光、蓝光; 零级、+1级、–1级分别用粗横线、竖短线、细横线表示
Fig. 2. Relative light intensity distributions of zero, positive and negative one order of diffraction beams of three primary colors, where the RPG period d and the cone angle of incident beam θ of (a)−(e) are (7.5°, 4 μm), (7.5°, 8 μm), (7.5°, 40 μm), (14.5°, 4 μm) and (102.5°, 4 μm) respectively. The red, green, blue line represents the red, green, blue light of three primary colors respectively. The zero, positive and negative one order beam is represented by the thick horizontal, vertical short and thin horizontal line respectively.
图 4 色散量C与矩形位相光栅周期d, LED入射光束锥角θ的关系 (a)矩形位相光栅周期d; (b) LED入射光束锥角θ; 其中红圆圈表示零色散的边界点
Fig. 4. Influences of grating period d and incident light cone angle θ on the dispersion parameter C. (a), (b) is the calculated relationship curve between C and d, or θ respectively, where the red circles represent the zero-dispersion boundary points.
表 1 不同入射光束锥角θ在观察平面中点所对应的衍射光谱的亮度值和色坐标
Table 1. Luminance and chromaticity coordinate of the center diffraction spectrum with different θ.
θ 102.68º 64.01º 45.24º 34.71º Luminance L (cd/m2) 6105.6 5161.2 4018.2 3262.1 Chromaticity Coordinate x = 0.2978y = 0.2828 x = 0.3024y = 0.2798 x = 0.3074y = 0.2779 x = 0.3025y = 0.2770 -
[1] Xu P, Huang H X, Wang K, Ruan S C, Yang J, Wan L L, Chen X X, Liu J Y 2007 Opt. Express 15 809
Google Scholar
[2] Xu P, Hong C Q, Cheng G X, Zhou L, Sun Z L 2015 Opt. Express 23 6773
Google Scholar
[3] Huang H X, Ruan S C, Yang T, Xu P 2015 Nano-Micro Lett. 7 177
Google Scholar
[4] 黄海漩, 徐平, 阮双琛, 杨拓, 袁霞, 黄燕燕 2015 物理学报 64 154212
Google Scholar
Huang H X, Xu P, Ruan S C, Yang T, Yuan X, Huang Y Y 2015 Acta Phys. Sin. 64 154212
Google Scholar
[5] 徐平, 袁霞, 杨拓, 黄海漩, 唐少拓, 黄燕燕, 肖钰斐, 彭文达 2017 物理学报 66 124201
Google Scholar
Xu P, Yuan X, Yang T, Huang H X, Tang S T, Huang Y Y, Xiao Y F, Peng W D 2017 Acta Phys. Sin. 66 124201
Google Scholar
[6] 徐平, 唐少拓, 袁霞, 黄海漩, 杨拓, 罗统政, 喻珺 2018 物理学报 67 024202
Google Scholar
Xu P, Tang S T, Yuan X, Huang H X, Yang T, Luo T Z, Yu J 2018 Acta Phys. Sin. 67 024202
Google Scholar
[7] Thomschke M, Reineke S, Lüssem B, Leo K 2012 Nano Lett. 12 424
Google Scholar
[8] Siitonen S, Laakkonen P, Vahimaa P, Kuittinen M, Tossavainen N 2006 Appl. Opt. 45 2623
Google Scholar
[9] Zhao X N, Hu J P, Lin Y, Xu F, Zhu X J, Pu D L, Chen L S, Wang C H 2016 Sci. Rep. 6 28319
Google Scholar
[10] Xue C X, Cui Q F 2010 Opt. Lett. 35 986
Google Scholar
[11] Tsukamoto H, Nishiyama M 2006 Jpn. J. Appl. Phys. 45 6678
Google Scholar
[12] Xu P, Huang Y Y, Zhang X L, Huang J F, Li B B, Ye E, Duan S F, Su Z J 2013 Opt. Express 21 20159
Google Scholar
[13] Xu P, Huang Y Y, Su Z J, Zhang X L, Luo T Z, Peng W D 2015 Opt. Express 23 4887
Google Scholar
[14] Xu P, Yan Z L, Wan L L, Huang H X 2004 Proceedings of SPIE Holography Diffractive Optics and Applications Ⅱ Beijing, China, November 8−11, 2004 p66
[15] Park S R, Kwon O J, Shin D, Song S H, Lee H S, Choi H Y 2007 Opt. Express 15 2888
Google Scholar
[16] Yang X P, Yan Y B, Jin G F 2005 Opt. Express 13 8349
Google Scholar
[17] Caputo R, Sio L D, Jak M J J, Hornix E J, Boer D K G, Cornelissen H J 2007 Opt. Express 15 10540
Google Scholar
[18] Xu M, Urbach H P, Boer D K G 2007 Opt. Express 15 5789
Google Scholar
[19] 张以谟 2008 应用光学 (北京:电子工业出版社) 第7章
Zhang Y M 2008 Applied Optics (Beijing: Publishing House of Electronics Industry) Chapter 7 (in Chinese)
[20] 吕乃光 2006 傅里叶光学 (北京:机械工业出版社)第70−113页
Lü N G 2006 Fourier Optics (Beijing: China Machine Press) pp70−113 (in Chinese)
[21] 苏显渝, 李继陶 1999 信息光学 (北京:科学出版社) 第44页
Su X Y, Li J T 1999 Information Optics (Beijing: Science Press) p44 (in Chinese)
计量
- 文章访问数: 8037
- PDF下载量: 51