7 参考文献
[1] V. R.Almain,C.A.Barrios, R. R. Puccinia and M. Lipson,All optical control of
light on a silicon chip,Nature, vol. 431, pp. 1081-1084, 2004. [2] Arthur Ashkin,How it all begain,nature photonics,vol5,June 2011.
[3] M. NietoVesperinas, P. C. Chaumet and A. Rahmani, Near-field photonic forces, Phil. Trans. R. Soc. Lond. Ser. A, vol. 362, pp. 719-737, 2004.
[4] K.C.Neuman and A. Nagy, Single-molecule force spectroscopy: optical tweezers,
magnetic tweezers and atomic force microscopy,Nat. Methods 5,491-505 (2008). [5] Chu, S. Laser manipulation of atoms and particles. Science 253, 861–866 (1991). [6] Pernice, W. H. P., Li, M. & Tang, H. X. Theoretical investigation of the transverse
optical force between a silicon nanowire waveguide and a substrate. Opt. Express 17, 1806–1816 (2009).
[7] Peter T. Rakich and Milo. Popovi.Engineering optical forces in waveguides and cavities based on optical response.Opt. Express,1130-1145(2009).
[8] Notomi,M.&Mitsugi,S.Wavelength conversion via dynamic refractive index tuning of a cavity.Phys. Rev. A 73, 051803(R) (2006).
[9] Li, M. et al. Harnessing optical forces in integrated photonic circuits. Nature 456, 480–484 (2008).
[10] M. Ohkawa, M. Izutsu and T. Sueta, Integrated optic pressure sensor on silicon substrate,Appl. Opt., vol. 28, pp. 5153-5157, 1989.
[11] Dorsel, A.,McCullen.Optical bistability and mirror confinement induced by
radiation pressure. Phys. Rev. Lett. 51, 1550-1553 (1983). [12] Povinelli, M. L. et al. Evanescent-wave bonding between optical waveguides. Opt.
Lett. 30, 3042–3044 (2005).
[13] Favero, I. & Karrai,K. Optomechanics of deformable optical cavities. Nature Photon.3, 201–205 (2009). [14] Dries Van Thourhout and Joris Roels,Optomechanical device actuation through the optical gradient force,nature photonics,vol4,June 2010.
[15] A. Ashkin, Optical trapping and manipulation of neutral particles using lasers,Proc. Natl.Acad.Sci.USA,vol. 94,pp.4853–4860,1997.
[16] O. Arcizet, P.F. Cohadon, Radiation pressure cooling and optomechanical instability of a micromirror,Nature, vol. 444, pp.71-74, 2006.
[17] S. Lin, E.Schonbrun,and K. Crozier, Optical Manipulation with Planar Silicon
Microring Resonators,Nano Lett., vol. 10, pp. 2408-2411,2010.
[18] S. A. Hamilton, B. S. Robinson, 100 GB/s optical time-division multiplexed
networks, J. Light-wave Techno.20, 2086-2100 (2002). [19] P. T. Rakich,General Treatment of Optical Forces and Potentials in Mechanically
Variable Photonic Systems, Opt. Express, vol. 17, pp. 18116-18135, 2009.
[20] Barclay, P. E., Srinivasan, Nonlinear response of silicon photonic crystal microcavities excited via an integrated waveguide and fiber taper.Opt. Express 13, 801-820 (2005).
[21] M., Uys, H. & Meystre, P. Optomechanical trapping and cooling of partially
reflective mirrors. Phys. Rev. A 77, 033819 (2008). [22] Eichenield, A picogram and nanometer-scale photonic-crystal optomechanical
cavity. Nature 459, 550-555 (2009). [23] Hohberger Metzger, C. & Karrai, K. Cavity cooling of a micro-lever. Nature 432,
1002-1005 (2004).
[24] Anetsberger, G. et al. Near-i eld cavity optomechanics with nanomechanical
oscillators. Nature Phys. 5, 909-914 (2009).
[25] Kippenberg, T. J. & Vahala, K. J. Cavity opto-mechanics. Opt. Express 15, 17172-17205 (2007)
[26] Rakich, P., Corralling and spectral bonding of optical resonances through
optically induced potentials. Nature Photon. 1, 658-665 (2007).
[27] Pai, A.,P. & Vinet, Radiation pressure induced instabilities in laser
interferometric detectors of gravitational waves. Eur. Phys. J. D 8, 333-346 (1999).
[28] Povinelli, M. L. et al . High-Q enhancement of attractive and repulsive optical
forces between coupled whispering-gallery-mode resonators. Opt. Express 13,
8287-8295 (2005).
[29] Notomi, M., Optomechanical wavelength and energy conversion in high-Q
double-layer cavities of photonic crystal slabs. Phys. Rev. Lett. 97,023903 (2006).
[30] Gigan, S. et al. Self-cooling of a micromirror by radiation pressure. Nature
444, 67 -70 (2006).
8 附录
wave_in=1600;wave_n=1601.72; %wave_in 代表初始调试波长 wave_n代表微环初始本证谐振波长
del=1.72; Tl代表初始差值
while(del>=1e-6) %循环总体部分
r=30*10^-6; %r代表微环的半径
q=0.75; %代表微环固定部分所占的比例
k=0.1883*10^-10; %微环的总体弹性系数
f=@(x)force(wave_in,wave_n,q,r,k,x); %表示在输入参数为输入波长、本证波长、比例、半径、弹性系数下的光力表达式 具体内容参见fuction函数
x=fzero(f,60); %求解形变量
n1=2.2079;wave_nn=1601.72; %微环初始有效折射率和微环本证波长
a=(4*pi^2*r*n1)./wave_nn; %光在微环中的行程相位
n_effc=2.2+0.108*exp(-0.0125.*(200-x)); %变形导致的有效折射率的变化值
b=(4*q*pi^2*r*2.2089)+(4*(1-q)*pi^2*r*n_effc); %求解变形后的谐振波长
wave_y=b./a;
wave_in=wave_in+0.01; %调试波长不断的增加 wave_n=wave_y; %本证波长变化
del=wave_n-wave_in; %循环控制条件
[x,wave_in,wave_n] %输出形变量和最终波长 End
function f=force(wave_in,wave_n,q,r,k,x) %求解形变量的方程 delta=wave_n-wave_in; %两者之间的差值 q1=-0.2*delta+1; %耦合到微环里的比例 p=1.5625*10^-3; %直波导中的功率
B=16.2859; %微环谐振器的放大系数 power_ring=B*p*q1; %微环中总体的能量 length=(1-q)*2*pi*r; %悬空微环的长度 c=3*10^8; %真空中的光速
n_eff=0.1080*0.0125*10^9*exp(-0.0125.*(200-x)); %光力关于距离的导数
f=((n_eff*power_ring*length)/c)-k*x; %光力作用下的微环最终表达式
基于微环谐振器光力的研究
1 绪论
1.1 光力的产生
早在1871年，James就已经预测了光可以产生力[1]。随着现代实验技术的发展，现在我们已经认识到光可以产生两种力：散射力和梯度力[2]。散射力我们可以将它形象的认为是一束光子打在物体上，光子携带的动量会转移给物体，从而在物体表面作用一个力，这个力的方向与光路的方向平行，大小与光的频率和能量密度有关，我们把这种力称之为散射力。光散射力在微腔、微环、干涉仪等方面的应用已经被广泛研究。光梯度力，与动态变化的电磁场有关系，我们可以将其认为是场与场之间的力[3]。放在不断变化电磁场中的物体会被极化，形成一个偶极子，处在电场中的偶极子会受到电场力的作用，大小与电场的强度有关，方向指向场强更大的方向。光梯度力相较于光散射力有着更大的数值，且更容易控制方向，因此在实际中有着更加广泛的应用。众所周知的“光镊”的原理就是依据光梯度力，操作毫米纳米级层次的物体[4][5]。
1.2 光波导的原理
要想控制光沿着固定的方向和路径传播，就需要用合适的媒介来传播光。用来运输光的媒介我们称之为光波导。光波导是由光透明介质(如石英玻璃)构成的传输光频电磁波的导行结构。光波导的传输原理是在不同折射率的介质分界面上，电磁波的全反射现象使光波局限在波导及其周围有限区域内传播。只有特定角度、频率和模式的光才能在光波导中传输。为适应不同的用途，光波导有着各种各样的形状。当光在波导中传播时，电磁场并非全部局限在波导结构中，在波导周围有限区域内会有一个快速消散到空气中的消散场，我们把这个消散的场域称之为消散波。消散波是一个近场耦合，是许多应用的原理基础。 最近几年，由于对性能更好价格更低廉的新型光纤通信器件迫切需求以及波导制作技术和平面工艺水平的不断提高，微环谐振器在理论和实验方面得到了快速的发展，并成为构建和实现集成光子学功能器件的重要的基础光波导单元[6].环状波导的谐振效应使其具有独特的波长选择、高品质因子等特性。目前，微环谐振器结构在研究中已被大量地应用于制作激光调频器、光波导分插复用器、生物化学传感器、调制器、光开关等。
1.3 光力驱动器件的意义
通入到光波导中的光会通过消散波与介质层耦合，使介质层极化，从而产生光梯度力。在微纳米层次下，这种力可以使微结构弯曲，从而产生各种各样有趣的现象。目前，利用光梯度力驱动的设备已被广泛的研究。比如单根波导与基底之间的装置，两根波导之间相互吸引排斥以及波导与各种微腔之间的光机械关系已被深入研究[7][8]。如图1-1所示，在很多光波导结构中可以产生光梯度力[9]。当光波导中有光传播时，消散波会和周围的结构产生耦合，就会有光梯度力作用在波导结构上，使光波导结构发生变形。这一奇特的性质被用在了光驱动器件上，在光开关、光延时、光通信等方面均有着重要的用途。
(a)两根平行波导之间产生力 (b)微环谐振器和波导之间产生力
图1-1 光波导结构上作用的光梯度力
利用光力驱动器件，是实现全光控集成光路的一个重要组成部分。全光控制具有快速直接、功耗低、集成度高等众多优点，是实现下一代光通信技术的基石。研究相关领域，是实现我国前沿研究的重要部分。
1.4 光驱动器件的研究现状
目前在该领域已有大量研究，很多光驱动器件被提出和制造出来，如图1-2所示[10]
，其不同的工作特性也被深入研究。
(a)两根平行波导 (b)垂直耦合谐振器
图1-2 光驱动器件
然而由于该领域还在进一步的研究和探讨中，很多光驱动器件工作的内部机理尚未完全弄清。许多文献只是提出阐述了器件某方面优异的性能，未给出详细的解释。本文针对微环和光波导这一具体器件，其具体结构和相应参数如图1-3所示，
图1-3 本文研究的波导结构
来初步探讨其机械和光学性质的机理。目前在一领域，有关波导与微环之间光耦合关系的研究已经比较多，但是关于光波导与微环之间由于光耦合，微环与基底产生光力导致微环弯曲的研究却非常的少，相关的研究主要集中在波导之间的作用弯曲。而且目前的很多研究是探究在光散射力作用下光机械的相互反作用[11]，有关光梯度力作用下的作用机理目前尚不明确。尤其是有关机械变形方面，相关文献均未给出具体的细节。此外，光波导与微环之间的作用关系非常复杂，牵涉到光耦合、光机械耦合以及光机械反作用等重要问题，对这一领域的深入研究有利于加深对光力和光机械作用原理的深入理解。
1.5 本论文的主要工作
本文主要针对直波导和微环谐振器，从其内部机理入手，详细探讨微环变形对波导力学和光学性能的影响，给出相关现象的理论解释和计算方法。主要内容包括以下几个部分：绪论。概述光力的发展及相关知识的简介。说明了本文选题的背景及意义，最后对本论文的工作及论文章节的安排做了简要介绍。微环的理论基础。首先从基本的原理入手，研究了光在波导中传输的基本条件，然后对微环建立模型，研究了表征微环光传输特性的指标，得到了波导和微环之间的能量耦合
关系式，为后续章节的计算提供了理论基础。光力的计算和位移变形的求解。首先得到光力的理论表达式，然后探究影响光力的因子，在得到具体光力表达式的情况下利用能量原理求解最大变形量。对微环双稳态现象的解释。利用前面讨论的基本知识和并结合光和机械之间相互耦合反馈理论对双稳态现象给出解释。关于光波导结构尺寸的优化。光波导结构尺寸的变化不仅会影响谐振腔中光学性质的改变，而且对微环的机械应力变化也有着很大的影响，综合考虑波导尺寸对这两方面的影响，利用研究的结果对双稳态现象提出有利的建议。总结和展望。对本文得出的结论进行总结，指出在该课题研究中存在的一些问题，并指出下一步研究的方向。
2 微环的理论基础
2.1 引言
微环属于集成光学领域中的一种重要的光波导器件。它独特的光传输特性使它在光通信网络、传感器和谐振器等方面得到了大量的应用。本章首先介绍光波导基本的传播原理和模式耦合理论，然后采用近似分析法对全通型微环进行分析，得到表征微环特性的参量，最后分析参量随条件变化的关系。
2.2 微环谐振器的工作原理
微环实际上属于众多光学谐振器中的一种，它是通过反馈的方式将光送入到一个闭合的回路，当光波在回路中传播一个周期后的相位差是2?的整数倍时，由于谐振效应则回路中形成稳定的谐振模式。当谐振器的腔长较大的时候，在谐振器中会存在许多的谐振模式，每一个不同的模式则对应于不同的谐振波长，因此对于一个谐振器来说，会有着许多不同的谐振波长。微环谐振器一般包括全通型微环和上下型微环两种，本文研究的微环属于全通型微环。当光波经过微环和光波导的耦合区时，由于消散波的耦合，一部分光波耦合进入微环，一部分光波直接通过光波导，耦合进入微环的光在微环中传播一个周期到达耦合区时，会再次发生耦合，一部分光耦合到直波导中，另一部分继续在环中传播。微环的传输特性取决于光波在微环中的传输、微环和直波导的耦合因素，主要的原理为干涉效应和反馈作用。通过合理设计耦合区间的间距可以控制微环的耦合系数来精确设计所需要的谱线。实际的谐振器除了环形还可以是盘形、跑道型、多边形和球形等各种各样的形式，根据实际情况可选择不同的形式。
2.3 光在微环中的基本传播原理
在对微环进行建模之前，需要了解光在微环中的基本传播原理。基本的光波导结构是由一个较高折射率内核和较低折射率覆盖层组成的，如图2-1所示，其中
n0?n1。由基本的光学原理可以知道，光在光波导中传播必须满足全反射，即
n1sin??/2????n0。由这个关系式我们得到最大入射角?max?sin?1n12?n02。
xn0??n1z 图2-1 基本的光传播示意图
但是实践发现，光线仅仅满足这个条件还不能在光波导中传播，光线的相位面传播还需要满足一定的条件，经过相关推导，我们可以得到这样的关系式：
am??tan?kn1sin??22?2????1 (2-1) ?2sin??n12?n02式中k是波矢量，a是波导宽度，m是整数值，代表相应的模式，??是22n1中心波导和覆盖层之间的折射率之差。
对应于每一个折射角，就会对应一个m值。我们将m称之为模式。当m变化的时候，代表着不同模式的传导波，m?0代表光波导基本的传播模式。不同的传导波有着不同的光传输特性，会对光力的产生及其特性产生不同的影响。
为了更加精确的描述光，我们使用麦克斯韦方程组来描述光场变化的基本规律。我们一般采用电场强度E、磁场强度H、电通量密度D、磁通量密度B等四个矢量来描述。麦克斯韦方程组的微分形式可以写作成如下的形式： ??H??D?J （2-2） ?t ??E???B （2-3） ?t
5.2.2 截面参数的变化
截面参数对微环的影响比较复杂，截面参数不仅能影响微环的光学性质，还影响着微环的力学性能。在光学性质方面表现为会影响波导内光学模式的传输，导致有效折射率的变化，从而影响光力的变化，而同时力学性能的改变意味着微环弹性的变化，为解决这个问题，利用COMSOL软件对微环的力学性能进行研究，在给定力下的微环的弯曲量与截面高度的关系如图5-3所示：
图5-3波导截面高度与最大位移的关系
从图中可以看出，随着波导截面高度的增加，微环变形量逐渐减小，这是因为截面高度的增加，导致微环抗弯截面系数增加，因而使形变量减小。
截面长度和微环最大位移的关系如图5-4所示，从图中可以看出，随着波导截面长度的增加，微环变形量也是逐渐减小的，其原因也是因为抗弯截面系数的增加，但相对来说，截面高度的变化对微环形变量的影响更大，这是因为在抗弯截面系数的表达式中高度项的指数比较高。
图 5-4 波导截面长度与最大位移的关系
波导截面参数的变化对光学的影响也非常大，由于时间和精力有限，本文未对这一部分做出详细的研究，根据相关文献[26-27]，存在一个最优波导截面尺寸，使得产生的光力最大，通过考虑到以上力学性能的关系，就可以得到一个微环形变量最大的截面尺寸，从而可以得到更大范围的双稳态现象。
5.2.3 结构的材料属性对对微环变形的影响
微环结构的材料属性包括杨氏模量、密度等，这些属性同样会影响到微环的传输特性。目前，有关新型材料的研究也是一个热点，在光子晶体、新型等离子体方面都有着很多的研究。限于本文的研究水平，本文从杨氏模量和密度这两个条件入手，探究这两个条件对微环特性的影响。
利用COMSOL软件，在上文讨论的基础上，通过改变杨氏模量和密度，得到如图5-5和图5-6所示的关系。
图5-5杨氏模量和最大变形的关系
图5-6密度和最大形变量的关系
从图中可以看出，密度和杨氏模量的增加都会使微环的变形量减小，因而在选用材料的时候，尽量选用样式模量和密度都较大的材料。
5.3 本章总结
本章在上面工作的基础上，对波导相关的尺寸参数进行优化，首先研究了微环悬空长度对力学性能和光学性能的影响，给出对微环悬空长度的建议：在满足其他条件的情况下，微环悬空长度可适度增加。然后研究了截面高度和长度对微环力学性能的影响，最后简单探究了一下杨氏模量和密度对微环性能的影响，并提出了相应的建议。
6 总结与展望
光通信以其时空带宽高、抗电磁干扰能力强、互连密度大、功耗低以及集成度高等优点在光信号传输、光交换和光计算中被广泛研究与应用。光子电路作为光通信的一个重要器件已日益成为人们研究的重点。依靠于集成电路成熟的制造工艺和方法，集成光路得到了快速的发展。目前的集成光路的控制方法多集中在光电、光热、光化学、机械等方面，这些方法在实现光通信方面起到了一定的作用，但是却有着自身无法克服的缺点。全光控制的概念被提出，全光控制具有快速直接、功耗低、集成度高等众多优点。如何实现全光控制，是目前研究的热点[28-29]。 本文所研究的光梯度力作用下的微环谐振器，对于实现全光控制有着重要的作用。微环结构本身产生的光梯度力使得微环本身的传输特性发生改变，光学性能和力学性能相互耦合，相关文献已经发现[30]，由于光梯度力的作用，微环本身会产生双稳态现象。虽然这一现象已被发现，但是关于双稳态中具体的一些参数如何影响最终结果以及是否可以对结果进行优化等方面缺乏相应的工作。因而本文从光力和微环特性入手，通过相应的理论计算和仿真探究相关参数的影响并对结果进行优化，最后对微环双稳态的使用提出了相关的建议。本文的研究内容主要集中在以下几个方面:
1.首先研究光在波导中传播的基本原理，利用模式耦合理论得到波导之间耦合的传输矩阵，然后对全通型微环建立数学模型，通过研究光在微环中传输的特性，得到微环中功率、相位的表达式，对其静态和动态特性进行求解，最后得到表征微环特性的性能指标参数。这一部分主要是为了下文的计算提供理论基础。 2.微环变形的分析。首先通过能量转化过程得到光力的表达式，为了求解光力，还需要得到有效折射率随间距的变化关系，通过理论分析求解全矢量光场，得到折射率与间距的关系，仿真的结果与理论推算基本吻合。然后对微环的力学性能进行分析，通过仿真求解出微环的弹性系数，利用能量转化求解出微环变形的位移量。求得的位移量为75nm。最后，对微环振动性能进行了简单的分析，并讨论了相关参数对各阶振型的影响。 3.双稳态现象的解释。首先利用前面的知识得到微环谐振波长偏移量随间距变化的关系和环内耦合功率与控制波长与谐振波长的关系，然后将调制过程分成两种情况，对于每一种情况，给出过程的系统框图，并编写相关的程序，计算出最终的谐振波长值。计算得到的最终谐振波长值为1603.4nm。将计算结果与实验数据做比较，二者基本吻合。 4.波导尺寸的优化。对微环的悬空长度、截面的高度和长度以及材料的属性进行了研究。结果发现，增加微环的悬空长度可以有效地增加微环的偏移量，使得最终的谐振波长值增加，同时使得滞回线之间的差值增大，双稳态现象会更加明显。截面的高度和长度参数存在一个最优值使得光力最大，也意味着会使双稳态更加明显。材料的杨氏模量和密度与微环变形量成正比，因而在选取材料时尽量选取较大杨氏模量和密度的材料。
本论文在微环谐振器的具体数值计算以及相关参数的优化方面做了一些工作，取得了阶段性的成果，由于时间有限，还有一些工作没有进行深入的研究，希望在下一步的工作中能完成相关内容，并取得较好的成果。相关建议如下： 探究微环动态响应性能与光学性能之间的关系，目前相关的研究集中在光散射力
的范围，有关光梯度力的光机械耦合性能的研究还比较少，同时，由于光力做功，还会使机械阵子的频率得以加速或者冷却，这是一个非常有趣的研究课题，那么在下一步的研究中会逐渐接触这些研究。 有关双稳态现象的应用。本文针对双稳态现象给出了数值计算和参数优化，如何利用这些现象来更好地设计应用具体的器件，是一个值得关注的课题。本文打算利用这个现象来设计一个压力传感器，当外力作用在上面使得微环发生弯曲时，通过读取直波导光学性质的改变来检测较小的力。这只是一个简单的设想，在以后的研究中还会继续对其应用性进行深入的研究。
本文主要做的理论计算和软件仿真，缺乏相应的实验过程和实验操作，因而在以后的研究中，需要尽快掌握有关纳米制造一些方面的技术，做一些自己的器件，使得自己全面发展。
图3-5 微环动态响应示意图
仿真得到的微环振动频率在1-10MHz之间，与相关文献上的数据是吻合的，从而证明了仿真的有效性。在得到微环振动频率的响应之后，通过改变微环悬空部分的长度、微环横截面的高度和宽度来探究振动频率随参数变化的关系。
(a)当悬空长度增加时，如图3-6所示，各阶的振动频率会随着相应减小，总体来说，阶次越高，振动频率衰减的就越快。
图3-6悬空长度和振动频率的关系
(b)当截面长度增加时，截面长度和振动频率的关系如图3-7所示：
图3-7截面长度和振动频率的关系
从图中可以看出，当截面长度从200nm增加到400nm的时候，三阶和四阶振动频率增加明显，近似为线性关系，一阶、二阶和五阶振动频率无明显变化；当截面长度增加为400nm以上时，一阶、二阶、三阶和四阶无明显变化，而五阶振动频率出现大幅度增加。
(c)当截面高度变化时，截面高度和频率的变化关系如图3-8所示：
图3-8截面高度和振动频率的关系
从图中可以看出，微环振动频率相对于截面高度的变化更加复杂，在100nm到400nm之间，一阶和二阶振动频率近似线性增加，三阶和四阶振动频率存在局部
最大值，在400nm之后，二阶振动频率平稳不再增加，一阶和三阶四阶振动频率线性增加。
3.4 总结
本章针对光力的计算和微环形变量的问题进行了探讨，首先得到微腔光力的表达式，然后对表达式中的有效折射率进行理论和软件仿真，取得了理想的结果。针对光力是变力，难以求解微环形变量这一问题，利用能量守恒和仿真模拟法得到微环的最终形变量，并对结果进行了取舍和分析。最后对微环的动态响应频率进行了简单探讨，并分析了悬空长度、截面的高度和宽度对振动频率的影响关系。
4 微环双稳态现象的探讨
4.1 引言
微环在光梯度力的作用下发生弯曲，弯曲导致微环的光传输性质发生变化，光性质变化会使光力的大小发生变化，力大小变化会反作用于微环弯曲量，因此，对于微环变形来讲，这是一个光和机械变形强耦合的过程。如果逐渐增加或者减小来调制控制波导中光的波长，那么在微环中就会产生双稳态现象。所谓双稳态现象，是指对于同一个输入会有两种不同的稳态输出。对于本文来说，当调制条件不同时，对于同一种波长，会有不同的输出。利用以上章节所讨论的知识，本章针对这一问题给出解释。
4.2 微环谐振频率的改变
由第二章微环传输特性的讨论可以知道，微环有效折射率的改变会使微环谐振频率发生变化。由第三章关于光力的讨论可知，微环有效折射率随间距的变化近似服从负指数规律，综合这两种因素，我们就可以得到谐振频率随间距变化的关系。 已知光在微环中传播所满足的角度关系式为：
??k?neff?2?r?2?neff?2?r??42rneff ??2m? （4-1）
?其中m是整数。从式子中可以看出，对于特定的m值，有效折射率的改变会使谐振波长发生变化，有效折射率越大，谐振波长就会越大。
对于本文研究的微环，有34的微环固定，其折射率不变，14的微环悬臂，其有效折射率随间距变化而变化，将上式改写一下，我们得到：
?2r3?2rn1eff?n2eff?2m? （4-2） ??k?neff?2?r? ??其中n1eff是悬臂部分的有效折射率，n2eff是固定部分的有效折射率。 由第三章讨论可知，悬臂部分有效折射率随间距变化的关系为：
n2eff?n0?ne （4-3） ??g )1xp(
将该式代入到中，则得到谐振频率随间距变化的关系：
?2r3?2rn?(n0?n1exp(??g))?2m? （4-4） ??k?neff?2?r? ?1eff?利用matlab可以得到如图4-1所示的微环变形量和谐振波长的关系：
图4-1变形量和谐振波长的关系
从图中可以看出，当微环变形量超过100nm之后，谐振波长的变化就非常迅速，但实际上，当微环谐振波长的变化值超过微环自由光谱范围（FSR）的时候，就不能使用上述的关系式，因为对应的本证谐振波长会发生变化，因而，在微环变形量满足100nm以下才能使用本关系式。通过第二章的计算，我们得到最大变形量为75nm，是满足该关系式的，可以使用。
4.3 耦合光功率的计算
通过第二章对微环的分析，我们知道当控制波频率与微环谐振频率一样时，环内耦合的功率最大，二者差值越大，环内耦合的功率就会越小。它们之间的关系近似满足：
?PinPring?22wc?wr???????4? （4-5） 微环谐振频率随间距变化的关系在上一节我们已经研究过了，他们之间的关系满
足：
?2r3?2rn1eff?(n0?n1exp(??g))?2m? （4-6） ??那么由此，我们就可以得到环内耦合功率随间距变化的关系。它存在一个极值，当二者波长一致时，耦合功率最大，差值越大，耦合功率越小，且以谐振波长所在的直线为对称轴。 由此可见，微环的弯曲既会影响谐振波长的改变，还会影响环内耦合功率的改变，因此在微环弯曲时要综合考虑环内能量和有效折射率的关系。
4.4 双稳态的计算
当采用不同的方法调制输入到直波导中的控制波时，微环中会有不同的现象产生，下面分情况讨论这些现象，并给出解释。
1）当控制波长从1600nm逐渐增加的情况：
其整体过程是：控制波的能量耦合到环中，产生光力，光力使环弯曲，环弯曲导致环中有效折射率变化，使环的谐振频率变化，谐振频率的变化导致光耦合能量变化，从而反作用光力。这是一个相互耦合的过程。为了解决这一个问题，我们假定在给定控制波的波长时，环中能量不发生变化，当波长增加时，下一个时刻点的环谐振频率以上一时刻变化后的为准。 在给定控制波长的条件下，一方面因为本身微环弯曲量较小，从实验和仿真我们都可以看到，其变化量在几十个纳米的范围，微环谐振频率的变化值不明显，另一方面，耦合功率的关系式也决定了小数值变化对微环能量影响较小，因此，我们在给定波长的条件下不考虑环内能量变化。
那么问题就转变成当给定控制波长时，环中能得到一个确切的功率，这个环内功率在下一个控制波长到来之前一直保持稳定，在这个功率作用下，微环弯曲，谐振波长变化，然后在继续一下个循环。因此，当波长以0.01纳米的速度逐渐增加时，我们把这个过程看成是一个个给定波长下的迭代，当波长差为零的时候，迭代停止，此时我们可以得到微环谐振器的偏移量以及迭代结束时的控制波长。 通过以上章节的讨论，本文将这个过程用matlab语言编写成一个程序，程序的系统框图如图4-2所示：
图4-2 循环过程的系统框图
Matlab的源程序放在附录处。
通过调整相应参数，得到在调制波长从1600nm逐渐增加时，最后稳定波长为1603.4nm。这与实验数据图4-3所反映的结果是基本符合的。
图4-3控制波长增加的数据图
结果分析：图中表示当控制波长从左向右逐渐增加的时候，对应的微环三个谐振波长也会相应的增加，当控制波长超过1603.7nm时，微环的谐振波长没有发生变化，表示此时微环没有发生谐振。
2）当控制波长从1603纳米逐渐减少的情况：
其原理过程与波长逐渐增加的过程相似，区别在于当波长减少的时候，微环弯曲导致其谐振频率仍然是增加的，对于波长差会很快的衰减到零，因此相对控制波长增加的情况来说，其最终谐振波长会减小，且谐振器谐振频率值会变小。理论和实验数据图4-4吻合的非常好。
图4-4控制波长减小的数据图
结果分析：图中表示当控制波长从右向左逐渐减小的时候，对应的微环三个谐振波长也会相应的增加，当控制波长低于1601.7nm时，微环的谐振波长恢复到原来的位置，实验数据显示，在1602nm处达到了极值。
3）综合以上两个调制过程，最终得到图4-5所示的调制图。
图4-5 滞回线图
图中红线表示当调制波长从1600nm增加到1605nm微环谐振偏移量的变化。蓝线表示当调制波长从1605nm减小到1600nm时微环偏移量的变化。综合前面的计算结果和理论分析，我们发现了以下规律：
(一)微环发生最大谐振值时对应的控制波长等于微环未发生形变前的谐振波长和偏移量之和。这一点从计算过程中可以明显的看出来，当控制波长与谐振波长相等时，微环发生最大形变，此后微环会迅速的恢复到初始位置。
(二)红线对应的最大控制波长值比蓝线对应的最大控制波长值大。这是由控制波长的调制过程决定的。当调制波长逐渐增加的时候，微环的谐振波长也在增加，初始值谐振波长比调制波长大，而调制波长的增加幅度比谐振波长的大，因而在经过若干次迭代之后二者最终达到一个平衡值。当调制波长逐渐减小的时候，微环的谐振波长依然在增加，因而二者会很快达到平衡值。在图像上的表现为红线调制的偏移量比蓝线调制的偏移量大。
4.5 本章小结
本章首先研究了微环的传输特性得到谐振波长偏移量随着间距变化的关系，然后得到环内耦合功率与控制波长和谐振波长的关系，最后利用以上的关系对微环双稳态进行计算，将得到的计算值与实验值进行比较，并给出在调制过程中的一些规律。计算结果和实验数据基本吻合，达到了较为理想的结果。
5 光波导结构尺寸的优化
5.1 引言
前面几章综合论述了微环在光力作用下的特性，可以看出，目前而言光力的大小还处在纳米级的层次，如何增强光学力的量级以更好的促进光学力的应用，是很多研究机构目前正在研究的课题。相关研究已经证明[25]，改变光波导的结构尺寸能够改变光学力的大小，同时结构尺寸的变化会对微腔的光学性能和机械性能产生很大的影响。因此本章尝试探究改变光波导的结构尺寸参数来优化波导，以期望对更好的利用双稳态现象提出一些建议。
5.2 光波导结构变化的影响
光波导的结构尺寸包括悬空部分的比例、微环横截面的高度和长度以及微环的材料属性。这些尺寸的变化会首先影响到微环中光传播的性质，继而导致相关光学参数的改变，光学参数的改变会影响到光力的产生和大小，同时尺寸的变化还会影响到微环力学性能的变化，因而对最后的结果影响比较复杂。本章节在上面章节所取得的成果下，对参数进行优化。
5.2.1 悬臂长度的变化
本文针对研究的微环有四分之一的部分悬空，当悬空部分增加时，微环的弹性系数会降低，通过软件模拟，我们得到了图5-1所示的悬空长度和最大位移的关系。
图 5-1悬空长度和最大位移的关系
从图中可以看出，当微环悬空长度逐渐增加时，其最大位移量也增加，近似满足
线性关系，由F?Kx可知对于给定作用力下其刚度矩阵变小，最大位移量增加。那么由此我们得到了图5-2所示的悬空长度和微环弹性系数的关系。
图 5-2 悬空长度和微环弹性系数的关系
悬空长度增加，弹性系数减小，通过计算得到其最终的谐振波长会增加，其具体数值为k每减小0.01，变形量会增加10nm左右。
同时悬臂长度的增加导致微环固定部分与悬空部分的比例减小，由
??k?neff?2?r?4??2r?n1eff?4(1??)?2r?n2eff?2m? （5-1）
其中?是固定部分所占的比例，可知，当?减小时，微环谐振波长对有效折射率
的变化更加敏感，也意味着其谐振波长的改变值更大。
当??0.75时，得到最终谐振波长为1603.4nm，当??0.7时，得到最终谐振波长为1604.3nm，因此增加悬空长度可以有效的增大变形量，同时得到一个较大的最终谐振波长。但对于实际的器件来说，悬空长度太大不利于封装和制作，并且本文没有考虑到实际的稳定情况，因此在条件许可的情况下，建议适度增大悬空长度的比例。
综合以上的讨论，我们得到如下结论： [1] 增加悬空部分的比例，会使微环谐振器的变形量更明显，微环谐振波长变化
值更大。
[2] 增加悬空部分比例，会使控制波长逐渐增加时的最终谐振波长值变大，控制
波长逐渐减小时的最终谐振波长变小，在滞后回线上意味着滞后之差变大，双稳态现象更加明显。 [3] 增加悬空部分的比例，会使微环谐振波长偏移量增加，在滞后回线上意味着
峰值增加。