网站地图
下载:
-
在实际激光产生中，由于设备的原因会出现光束像散和扭曲，且一般是不可控的。为了描述这类光束，1993年，SIMON和MUKUNDA构造了一种部分相干轴对称高斯-谢尔模光束[1]，在其结构中引入了扭曲相位，被称为扭曲高斯-谢尔模(twisted Gaussian-Schell model, TGSM)光束。扭曲相位具有固定的手性或旋向性[2]，且携带有轨道角动量(orbital angular momentum, OAM)[3]，使得光斑在传输过程中绕轴旋转，这一新颖的特性引起了科研人员对TGSM光束的研究兴趣。1994年, FRIBERG在利用两组柱透镜组成的光学系统实验产生了TGSM光束[4]，并证明了这种典型的扭曲光束可分解为互不相关的椭圆高斯光束的叠加[5]。此后，科研人员长期致力于研究扭曲光束在鬼成像[6]、粒子捕获[7]、自由空间光通信[8]、照明光源[9]等方面的应用价值。
相干和相位是光波的基本特性，它们之间是相互影响的。过去很长一段时间对扭曲相位的研究局限在了TGSM光束这一种模型，这是由于扭曲相位是用一个与光场横截面两点位置相关的相位因子描述，新型扭曲光束的产生依赖于对光场横截面两点位置相干结构的精确调控。得益于GORI在2007年提出的构造特殊关联部分相干光束的方法[10]，空间关联结构光束的调控及产生的研究因此迎来爆发式增长，许多具有特殊关联结构的光束模型被提出，例如多高斯-谢尔模(muti-Gaussian-Schell model, MGSM)光束[11]，余弦高斯关联谢(cosine-Gaussian-Schell model, CGCSM)光束[12]，拉盖尔-高斯关联谢尔模光束(Laguerre-Gaussian correlation Schell model, LGCSM)光束等[13]，其中，MGSM光束以其在远场的平顶特性引起了广泛的关注，研究者对其在大气湍流[14]、随机介质[15]、海洋湍流[16]中的传输特性进行了研究，其对抑制光束在大气湍流中的闪烁[17-18]和粒子捕获[19]等方面也有重要作用。近年来，得益于BORGHI等人提出的轴对称谢尔模光束是否可携带扭曲相位的判断条件[20]，以及GORI等人提出的一种构建无对称约束扭曲光束的方法[21]，新型扭曲光束不断地被提出并进行了研究[22-30]，如径向偏振扭曲光束[22-24]、扭曲阵列光束[25-27]、扭曲涡旋高斯-谢尔模光束[28]、扭曲拉盖尔-高斯-谢尔模光束[29]，以及扭曲椭圆多高斯-谢尔模光束[30]。
随着移动互联网、物联网和云计算等技术的高速发展，人们对于通信容量和速率的要求越来越高。扭曲光束携带有轨道角动量，为高速光纤通信提供了新的信息调制维度，将极大地提升系统的信道容量和频谱利用率。基于模式分解理论，部分相干光可看作多模式相干光的非相干叠加，多模光纤可用于传输多种模式的光束，常见的多模光纤为阶跃折射率和梯度折射率(gradient index, GRIN)光纤。梯度折射率光纤具有自聚焦特性，脉冲色散低于阶跃折射率光纤，在光学制造[31]、自成像[32]、光束整形[33]等领域有广阔的应用前景。近些年来，研究者分别对电磁高斯-谢尔模光束[34-35]、电磁多高斯-谢尔模光束[36]、部分相干涡旋光束[37]、拉盖尔-高斯光束[38]在梯度折射率光纤中的传输特性进行了研究。然而，作为一种新型光束，扭曲多高斯光束在梯度折射率光纤中的传输特性鲜被研究。本文中基于Mercer展开理论[39]和交叉谱密度函数的非负性原则，证明了多高斯-谢尔模关联结构可携带扭曲相位，得到了一种椭圆扭曲多高斯-谢尔模光束，推导了其通过ABCD光学系统后交叉谱密度函数的解析式，并详细研究了其在梯度折射率光纤中传输时光强和相干度的演化特性。
-
部分相干光束的2阶统计特性可用交叉谱密度函数来表示，多高斯-谢尔模光束在源平面上的交叉谱密度函数(cross-spectral density, CSD)可以表示为[11]:
$ W\left(\boldsymbol{r}_1, \boldsymbol{r}_2\right)=\sqrt{S\left(\boldsymbol{r}_1\right) S\left(\boldsymbol{r}_2\right)} \eta\left(\boldsymbol{r}_1, \boldsymbol{r}_2\right) $
(1) 式中, S(r)表征强度轮廓分布函数，r1≡(x1, y1)，r2≡(x2, y2)是源平面上的任意两点的空间位置坐标，η(r1, r2)表征光束相干结构，多高斯-谢尔模关联结构可表示为:
$ \eta\left(\boldsymbol{r}_1, \boldsymbol{r}_2\right)=\frac{1}{C_0} \sum\limits_{m=1}^M\left(\begin{array}{l} M \\ m \end{array}\right) \frac{(-1)^{m-1}}{m} \exp \left[-\frac{\left(\boldsymbol{r}_2-\boldsymbol{r}_1\right)^2}{2 m \delta^2}\right] $
(2) 式中, M表征多高斯关联函数的模数，为归一化参数，表示二项式系数，σ为光束束腰宽度，δ为光束的相干长度。当M=1时，(1)式转化为高斯-谢尔模光束的交叉谱密度函数表达式。
BORGHI等人证明了光束是否能携带扭曲相位只与相干结构有关，而与光束的振幅无关[20]。因此，为了便于分析，可将(1)式携带扭曲相位，并忽略掉强度可表示为：
$ \bar{W}\left(\boldsymbol{r}_1, \boldsymbol{r}_2\right)=\eta\left(\left|\boldsymbol{r}_1-\boldsymbol{r}_2\right|\right) \exp \left(-\mathrm{i} \mu \boldsymbol{r}_1 \times \boldsymbol{r}_2\right) $
(3) 式中, μ表示扭曲因子。根据部分相干光束模式叠加理论，将(3)式用Mercer展开表示为[39]:
$ \bar{W}\left(\boldsymbol{r}_1, \boldsymbol{r}_2\right)=\sum\limits_{j=0, 1 / 2, 1, \cdots} \sum\limits_{l=-j}^j \varLambda_{j+l} \varPhi_{j, l}{ }^*\left(\boldsymbol{r}_1\right) \varPhi_{j, l}\left(\boldsymbol{r}_2\right) $
(4) 其特征函数Φj, l为正交，Φj, l*表示特征函数的复共轭，特征值Λj+l为非负，当部分相干光束携带扭曲相位后，其特征函数形式可为拉盖尔-高斯型[39]:
$ \begin{gathered} \varPhi_{j, l}=\sqrt{\frac{\mu}{\pi}}\left[\frac{(j-|l|) !}{(j+|l|) !}\right]^{1 / 2}(r \sqrt{\mu})^{2|l|} \times \\ \quad \exp (\mathrm{i} 2 l \theta) \mathrm{L}_{j-|l|}^{2|l|}\left(\mu r^2\right) \exp \left(-\frac{\mu r^2}{2}\right) \end{gathered} $
(5) 式中, j=0, 1/2, 1, …; l=－j, －j+1, …, j; Lj－|l|2|l|为拉盖尔多项式，l代表拓扑荷数，j－|l|为径向指数，2|l|为角向指数，两指数均为非负整数，为方便接下来的证明，(4)式中的空间坐标改写为极坐标形式r≡(r, θ); exp(i2lθ)代表每个模式均带有涡旋相位。从Mercer展开形式可看出，每个模式均携带有涡旋相位，揭示了携带轨道角动量的扭曲光束与涡旋光束之间的内在联系。接下来只需验证W(r1, r2)模式展开后所有特征值是否都是非负实数，且特征值满足Λs=Λj+l，特征值序列{Λs}s=0∞可由相干结构函数积分得到[20]:
$ \varLambda_s=\int \eta(|\boldsymbol{r}|) \mathrm{L}_s\left(\mu r^2\right) \exp \left(-\frac{\mu r^2}{2}\right) \mathrm{d}^2 \boldsymbol{r} $
(6) 式中, r=r2－r1，Ls表示角向指数为0、径向指数为s的拉盖尔多项式，将(2)式代入(6)式计算得到:
$ \begin{gathered} \varLambda_s=\frac{1}{C_0} \sum\limits_{m=1}^M\left(\begin{array}{l} M \\ m \end{array}\right) \frac{(-1)^{m-1}}{m} \times \\ \int_0^{2 \pi} \mathrm{d} \theta \int_0^{\infty} \exp \left[-\frac{r^2}{2 m \delta^2}-\frac{\mu r^2}{2}\right] \mathrm{L}_s\left(\mu r^2\right) r \mathrm{~d} r \end{gathered} $
(7) 最终积分得到特征值序列表达式为:
$ \varLambda_s=\frac{1}{C_0} \sum\limits_{m=1}^M\left(\begin{array}{l} M \\ m \end{array}\right)(-1)^{m-1} 2 \pi \delta^2 \frac{\left(1-\mu m \delta^2\right)^s}{\left(1+\mu m \delta^2\right)^{s+1}} $
(8) 对(8)式分析后可知，当μ≤1/mδ2时，该特征值为非负实数，此时椭圆多高斯-谢尔模光束可以携带扭曲相位。为便于之后的分析，将强度轮廓函数写为椭圆高斯函数，就得到了一类新型椭圆扭曲多高斯-谢尔模(twisted multi-Gaussian-Schell model, TMGSM)光束，其交叉谱密度函数可以表示为:
$ \begin{gathered} W\left(\boldsymbol{r}_1, \boldsymbol{r}_2\right)=\frac{1}{C_0} \sum\limits_{m=1}^M\left(\begin{array}{l} M \\ m \end{array}\right) \frac{(-1)^{m-1}}{m} \times \\ \exp \left(-\frac{r_{1 x}^2+r_{2 x}^2}{4 \sigma_x^2}-\frac{r_{1 y}{ }^2+r_{2 y}{ }^2}{4 \sigma_y^2}\right) \times \\ \exp \left[-\frac{\left(\boldsymbol{r}_1-\boldsymbol{r}_2\right)^2}{2 m \delta^2}\right] \exp \left(-\mathrm{i} \mu \boldsymbol{r}_1 \times \boldsymbol{r}_2\right) \end{gathered} $
(9) 式中, r1x，r1y，r2x，r2y分别为r1，r2在x，y方向上的分量; σx，σy表示x，y方向上的束腰宽度。在傍轴近似下利用Collins公式，椭圆TMGSM光束在像散ABCD光学系统中的传输公式可表示为:
$ \begin{gathered} W\left(\boldsymbol{\rho}_1, \boldsymbol{\rho}_2, z\right)=\frac{1}{\lambda^2 B^2} \iint W\left(\boldsymbol{r}_1, \boldsymbol{r}_2\right) \times \\ \exp \left[-\frac{\mathrm{i} k}{2 B}\left(A \boldsymbol{r}_1^2-2 \boldsymbol{r}_1 \boldsymbol{\rho}_1+D \boldsymbol{\rho}_1^2\right)\right] \times \\ \exp \left[\frac{\mathrm{i} k}{2 B}\left(A \boldsymbol{r}_2^2-2 \boldsymbol{r}_2 \boldsymbol{\rho}_2+D \boldsymbol{\rho}_2^2\right)\right] \mathrm{d}^2 \boldsymbol{r}_1 \mathrm{~d}^2 \boldsymbol{r}_2 \end{gathered} $
(10) 式中, ρ1=(ρ1x, ρ1y)和ρ2=(ρ2x, ρ2y)是输出平面上任意两个位置向量; A, B, C, D为光学系统传递矩阵的元素; λ为光束波长，k=2π/λ是波数。将(9)式代入(10)式，经过积分，就得到了输出平面的交叉谱密度函数解析式:
$ \begin{gathered} W\left(\boldsymbol{\rho}_1, \boldsymbol{\rho}_2, z\right)=\sum\limits_{m=1}^M\left(\begin{array}{l} M \\ m \end{array}\right) \frac{(-1)^{m-1}}{m} \times \\ \frac{k^2}{4 B^2 C_0 \sqrt{N_1 N_3\left(N_4-\varOmega_2\right)\left(N_5-\varOmega_1\right)}} \exp \left[\frac{\mathrm{i} D k}{2 B}\left(\boldsymbol{\rho}_2^2-\boldsymbol{\rho}_1^2\right)\right] \times \\ \exp \left[-\frac{u_1^2}{4 N_1}-\frac{v_1^2}{4 N_3}+\frac{\left(\mathrm{i} v_2-\varPi_2\right)^2}{4\left(N_4-\varOmega_2\right)}+\frac{\left(\mathrm{i} u_2-\varPi_3\right)^2}{4\left(N_5-\varOmega_1\right)}\right] \end{gathered} $
(11) 其中，使用了一系列参数替换：
$ \begin{gathered} N_1=\frac{1}{4 \sigma_x^2}+\frac{1}{2 m \delta^2}+\frac{\mathrm{i} A k}{2 B} ; N_2=\frac{1}{4 \sigma_x^2}+\frac{1}{2 m \delta^2}-\frac{\mathrm{i} A k}{2 B} ; \\ N_3=\frac{1}{4 \sigma_y^2}+\frac{1}{2 m \delta^2}+\frac{\mathrm{i} A k}{2 B}; \end{gathered} $
$ \begin{gathered} N_4=\frac{1}{4 \sigma_y^2}+\frac{1}{2 m \delta^2}-\frac{\mathrm{i} A k}{2 B} ; N_5=N_2-\frac{\eta_1^2}{4\left(N_4-\varOmega_2\right)} ; \\ u_1=\frac{k \rho_{1 x}}{B} ; u_2=\frac{k \rho_{2 x}}{B} ; v_1=\frac{k \rho_{1 y}}{B} ; v_2=\frac{k \rho_{2 y}}{B} ; \end{gathered} $
$ \begin{gathered} \varOmega_1=\frac{1}{4 N_1 m^2 \delta^4}-\frac{\mu^2}{4 N_3} ; \varOmega_2=\frac{1}{4 N_3 m^2 \delta^4}-\frac{\mu^2}{4 N_1} ; \\ \varPi_1=\frac{i u_1}{2 N_1 m \delta^2}-\frac{v_1 \mu}{2 N_3} ; \end{gathered} $
$ \begin{gathered} \varPi_2=\frac{\mathrm{i} v_1}{2 N_3 m \delta^2}+\frac{u_1 \mu}{2 N_1} ; \varPi_3=\varPi_1+\frac{\left(\varPi_2-\mathrm{i} v_2\right) \eta_1}{2\left(N_4-\varOmega_2\right)} \\ \eta_1=\frac{\mathrm{i} \mu}{2 N_3 m \delta^2}-\frac{\mathrm{i} \mu}{2 N_1 m \delta^2} \end{gathered} $
(12) 接下来，本文中将研究椭圆TMGSM光束在梯度折射率光纤中的传输特性。梯度折射率光纤的折射率分布可由下式表示:
$ n^2=\left\{\begin{array}{l} n_0{ }^2\left(1-\beta^2 \rho^2\right), \left(\rho^2 \leqslant R_0{ }^2\right) \\ n_0{ }^2\left(1-\beta^2 R_0{ }^2\right), \left(\rho^2>R_0{ }^2\right) \end{array}\right. $
(13) 其中,
$ \beta=\frac{1}{R_0}\left[1-\frac{n_1^2}{n_0^2}\right]^{\frac{1}{2}} $
(14) 式中，R0为纤芯半径，ρ为空间任意点距纤芯中心的距离，β为梯度折射率系数，n0为光纤中心折射率，n1为光纤包层折射率，梯度折射率光纤的ABCD传输矩阵可表示为:
$ \left[\begin{array}{ll} A & B \\ C & D \end{array}\right]=\left[\begin{array}{cc} \cos (\beta z) & \frac{\sin (\beta z)}{n_0 \beta} \\ -n_0 \beta \sin (\beta z) & \cos (\beta z) \end{array}\right] $
(15) 在这里用一种特定的梯度折射率光纤，纤芯为掺杂锗的二氧化硅(掺杂质量分数为0.079的氧化锗)，包层为二氧化硅制成，折射率由塞米尔方程计算:
$ n^2(\omega)=1+\sum\limits_{p=1}^3 \frac{B_p \omega_p{ }^2}{\omega_p^2-\omega^2} $
(16) 式中，ωp=2πc/λp，Bp，λp均为特定材料下的塞米尔系数，c为光速，ω为λ=632.8nm时对应的角频率。对纯二氧化硅，塞米尔系数B1=0.6961663, B2=0.4079426, B3=0.8974794，λ1=0.0684043μm, λ2=0.1162414μm, λ3=9.896161μm。对掺杂质量分数为0.079锗的二氧化硅，其塞米尔系数B1=0.7136824, B2=0.4254807, B3=0.8964226，λ1=0.0617167μm, λ2=0.1270814μm，λ3=9.896161μm。取光纤纤芯半径为25μm，通过(13)式、(14)式和(16)式可以得到n0=1.46977, n1=1.45702, β=5.25726mm－1，从而确定ABCD传输矩阵的具体参数。部分相干光束的光强和相干度可由交叉光谱密度函数表示为：
$ S(\boldsymbol{\rho})=W(\boldsymbol{\rho}, \boldsymbol{\rho}) $
(17) $ \eta\left(\boldsymbol{\rho}_1, \boldsymbol{\rho}_2\right)=\frac{W\left(\boldsymbol{\rho}_1, \boldsymbol{\rho}_2\right)}{\sqrt{S\left(\boldsymbol{\rho}_1\right)} \sqrt{S\left(\boldsymbol{\rho}_2\right)}} $
(18) -
本节中，将对椭圆TMGSM光束在梯度折射率光纤中传输的光强和相干度(degree of coherence, DOC)进行模拟和分析，除特殊说明外，初始参数设定为：λ=632.8nm，k=2π/λ，n0=1.46977，n1=1.45702，σx=100/k，σy=25/k，δ=50/k, β=5.25726mm－1，μ=4×10－5k2，M=5，L=π/β=597.5722μm。图 1a和图 1b是椭圆TMGSM光束在梯度折射率光纤中传输时的x-z归一化光强剖面图和y-z归一化光强剖面图。结果表明，初始的椭圆光斑在传输过程中发生周期为L的变化，在z=0处σx>σy，而在z=0.5tL(t=1, 3, 5…)的焦平面处σx < σy，这种现象是由梯度折射率光纤的自聚焦特性决定的。在图 1b中可以清楚地看到反常规的聚焦现象，即y方向的光斑宽度在聚焦过程中逐渐增大，这是由于光斑旋转造成的。此外，当光束传输到焦平面时，光强分布呈平顶分布，这种现象是由于多高斯关联结构造成的[11]。
Figure 1. Normalized intensity distribution of elliptically twisted multi-Gaussian-Schell model beams propagating in a gradient index fiber
图 2是不同多高斯模数的椭圆TMGSM光束在梯度折射率光纤中传输不同距离的归一化光强密度图。除M外，其它取值与初始参数一致。当M=1时，椭圆TMGSM光束转化为TGSM光束。从图中可以发现，光束光斑随传输距离的增加发生顺时针旋转，传输到焦平面z=0.5L处刚好转过π/2。同时，光斑的椭圆度也在传输过程中不断降低。观察图 2f、图 2l和图 2r，发现随M的增加不会影响光斑的旋转速度，但光斑在焦平面处的平顶分布区域增大，同时光束发散也随之增大，在焦平面处得到了较大的聚焦平顶光斑。
Figure 2. Normalized light intensity diagrams of elliptically twisted multi-Gaussian-Schell model beams propagating at different distances in gradient index fibers with different M
图 3是不同扭曲因子影响下，椭圆TMGSM光束在梯度折射率光纤中传输不同距离时的归一化光强密度图。除扭曲因子μ以外，其它参数均取初始参数，且扭曲因子取值满足μ≤1/(mδ2)的条件限制。结果表明: 当μ=0时，光束退化为MGSM光束，光斑在传输过程中只存在聚焦过程，不会发生旋转; 当μ>0时，光斑发生旋转，但旋转速度并不均匀，在刚开始传输时，随着扭曲因子增大，光斑的旋转速度加快，而传输到接近焦平面时，扭曲因子小的反而旋转速度更快，在焦点z=0.5L处，无论扭曲因子取何值，光斑总是刚好旋转了π/2。此外，观察图 3e、图 3k和图 3q发现，在接近焦平面处，光斑的椭圆度也随着扭曲因子的增大而增大。
Figure 3. Normalized light intensity diagrams of elliptically twisted multi-Gaussian-Schell model beams propagating at different distances in gradient index fibers under different twist factors μ
接下来研究了椭圆TMGSM光束在梯度折射率光纤中传输时相干度的演化特性。图 4a和图 4b是椭圆TMGSM光束在梯度折射率光纤中传输时的x-z相干度剖面图和y-z相干度剖面图, 其它参数与图 1一致。图 4显示，在梯度折射率光纤中传输时，相干度也会发生周期为L的变化。初始的相干度分布相同的情况下，由于光束强度的各向异性对相干性的影响，在传输过程中将导致相干度逐步呈现各向异性，由圆形高斯分布演化为椭圆高斯分布。
Figure 4. Degree of coherence distribution of elliptically elliptically twisted multi-Gaussian-Schell model beams propagating in a gradient index fiber
图 5是不同模数M下的椭圆TMGSM光束在梯度折射率光纤中传输不同距离的相干度密度图, 其它参数与图 2一致。从图中可发现，随着传输距离的增加，光束相干度分布由圆形高斯分布演化为椭圆高斯分布，椭圆度逐渐增大，且相干度分布逐渐逆时针旋转，传输到焦平面z=0.5L处刚好转过π/2。观察图 5a~图 5c得出，随着M的增加，相干度分布轮廓变小，周围开始出现暗环。与光强变化规律类似，增大M不会影响相干度分布的旋转速度。图 6是椭圆TMGSM光束在梯度折射率光纤中传输不同距离时在不同扭曲因子μ影响下的相干度密度图，参数同图 3。观察图 6a、图 6g和图 6m可得出扭曲因子的大小对初始相干度分布没有影响。图 6表明，当μ=0时，随着传输距离的增加，相干度分布不发生旋转; 当μ>0时，相干度分布发生旋转，但旋转速度并不均匀，与光强相反，在刚开始传输时，随着扭曲因子增大，相干度分布的旋转速度变慢，而传输到接近焦平面时，旋转速度加快，在焦平面z=0.5L处，无论扭曲因子取何值，相干度分布也是刚好旋转了π/2。
Figure 5. Degree of coherence density diagrams of elliptically twisted multi-Gaussian-Schell model beams propagating at different distances in gradient index fiber under different M
Figure 6. Degree of coherence density diagrams of elliptically twisted multi-Gaussian-Schell model beams propagating at different distances in gradient index fiber under different twist factors μ
-
本文中基于Mercer模式展开理论和交叉谱密度函数的非负性原则，严格证明了当扭曲因子满足μ≤1/(mδ2)时，多高斯-谢尔模关联结构可携带扭曲相位，得到了一种椭圆TMGSM光束，推导了其通过ABCD光学系统时交叉谱密度函数的解析式，并研究了其在梯度折射率光纤中传输时模数、扭曲因子对光强和相干度变化特性的影响。结果表明: 椭圆TMGSM光束在梯度折射率光纤中传输时光强和相干度均呈现周期性变化；光强分布呈现为由椭圆高斯分布到平顶分布再恢复为椭圆高斯分布的周期性变化，模数M越大，光强分布的平顶区域越大，相干度分布轮廓越小；由于扭曲相位的存在，光强和相干度分布均在传输过程中发生旋转，同时，旋转速度是非线性的并与扭曲因子的大小有关，且旋转角度在焦平面处限定为π/2。本文中的研究结果在光纤通信、聚焦成像、光学捕获等方面具有潜在的应用前景。
扭曲多高斯光束在梯度折射率光纤中的传输特性
Propagation characteristics of twisted multi-Gaussian beams in gradient index fibers
-
摘要: 为了得到一种新型椭圆扭曲多高斯-谢尔模(TMGSM)光束, 采用Mercer模式展开的方法进行了理论分析和验证, 证明了多高斯-谢尔模关联结构可携带扭曲相位, 详细研究了其在梯度折射率光纤中传输时的光强和相干度演化。结果表明, 椭圆TMGSM光束在梯度折射率光纤中传输时, 光强和相干度分布随着传输距离的增加发生周期性旋转, 并在0.5L(L为周期)的整数倍处偏转π/2, 其旋转角速度呈非线性变化且与扭曲因子的大小有关; 增大多高斯模数, 焦平面处光强分布的平顶区域增大, 相干度分布轮廓变小。此研究结果在光纤通信、聚焦成像、光学捕获等方面具有潜在的应用前景。
-
关键词:
- 激光光学 /
- 扭曲相位 /
- Mercer模式展开理论 /
- 多高斯关联 /
- 梯度折射率光纤
Abstract: In order to obtain a new elliptically twisted multi-Gaussian-Schell model beam, a Mercer model expansion method was adopted, the theoretical analysis and verification were carried out. It is proved that the twisted phase can be carried in the multi-Gaussian-Shell model correlation structure, and the intensity and coherence evolution of the beam propagating in gradient index fibers was studied. The results show that when elliptically twisted multi-Gaussian-Schell model beams are propagating in the graded index fibers, the distribution of intensity and coherence rotates periodically with the increase of transmission distance, and the deflection is π/2 degrees at the integer multiples of 0.5L (L is a cycle). The rotation angular velocity is nonlinear and related to the size of the twist factor. With the increase of multi-Gaussian modulus, the flat top region of the intensity distribution at the focal plane increases, and the coherence distribution contour becomes smaller. The research results have potential applications in optical fiber communication, focused imaging, optical capture and so on. -
Figure 1. Normalized intensity distribution of elliptically twisted multi-Gaussian-Schell model beams propagating in a gradient index fiber
a—z-x profile b—z-y profile
Figure 2. Normalized light intensity diagrams of elliptically twisted multi-Gaussian-Schell model beams propagating at different distances in gradient index fibers with different M
a~f—M=1 g~l—M=5 m~r—M=10
Figure 3. Normalized light intensity diagrams of elliptically twisted multi-Gaussian-Schell model beams propagating at different distances in gradient index fibers under different twist factors μ
a~f—μ=0 g~l—μ=4×10-5k2 m~r—μ=8×10-5k2
Figure 4. Degree of coherence distribution of elliptically elliptically twisted multi-Gaussian-Schell model beams propagating in a gradient index fiber
a—z-x profile b—z-y profile
Figure 5. Degree of coherence density diagrams of elliptically twisted multi-Gaussian-Schell model beams propagating at different distances in gradient index fiber under different M
a—M=1 b—M=5 c—M=10
-
[1] SIMON R, MUKUNDA N. Twisted Gaussian Schell-model beams[J]. Journal of the Optical Society of America, 1993, A10(1): 95-109. [2] SIMON R, MUKUNDA N. Twist phase in Gaussian-beam optics[J]. Journal of the Optical Society of America, 1998, A15(9): 2373-2382. [3] SERNA J, MOVILLA J M. Orbital angular momentum of partially coherent beams[J]. Optics Letters, 2001, 26(7): 405-407. doi: 10.1364/OL.26.000405 [4] FRIBERG A T, TERVONEN E, TURUNEN J. Focusing of twisted Gaussian Schell-model beams[J]. Optics Communications, 1994, 106(4/6): 127-132. [5] FRIBERG A T, TERVONEN E, TURUNEN J. Interpretation and experimental demonstration of twisted Gaussian Schell-model beams[J]. Journal of the Optical Society of America, 1994, A11(6): 1818-1826. [6] DING G L, LV B D. Propagation of twisted Gaussian Schell-model beams through a misaligned first-order optical system[J]. Journal of Modern Optics, 2009, 48(10): 1617-1621. [7] ZHAO Ch L, CAI Y J, KOROTKOVA O. Radiation force of scalar and electromagnetic twisted Gaussian Schell-model beams[J]. Optics Express, 2009, 17(24): 21472-21487. doi: 10.1364/OE.17.021472 [8] WANG F, CAI Y J, EYYUBOGLU H T, et al. Twist phase-induced reduction in scintillation of a partially coherent beam in turbulent atmosphere[J]. Optics Letters, 2012, 37(2): 184-186. doi: 10.1364/OL.37.000184 [9] TONG Zh S, KOROTKOVA O. Beyond the classical Rayleigh limit with twisted light[J]. Optics Letters, 2012, 37(13): 2595-2597. doi: 10.1364/OL.37.002595 [10] GORI F, SANTARSIERO M. Devising genuine spatial correlation functions[J]. Optics Letters, 2007, 32(24): 3531-3533. doi: 10.1364/OL.32.003531 [11] SAHIN S, KOROTKOVA O. Light sources generating far fields with tunable flat profiles[J]. Optics Letters, 2012, 37(14): 2970-2972. doi: 10.1364/OL.37.002970 [12] MEI Zh R, KOROTKOVA O. Cosine-Gaussian Schell-model sources[J]. Optics Letters, 2013, 38(14): 2578-2580. doi: 10.1364/OL.38.002578 [13] MEI Zh R, KOROTKOVA O. Random sources generating ring-shaped beams[J]. Optics Letters, 2013, 38(2): 91-93. doi: 10.1364/OL.38.000091 [14] DU Sh C, YUAN Y S, LIANG C H, et al. Second-order moments of a multi-Gaussian Schell-model beam in a turbulent atmosphere[J]. Optics & Laser Technology, 2013, 50: 14-19. [15] ZHANG Y Y, ZHAO D M. Scattering of multi-Gaussian Schell-model beams on a random medium[J]. Optics Express, 2013, 21(21): 24781-24792. doi: 10.1364/OE.21.024781 [16] ZHU W T, TANG M M, ZHAO D M. Propagation of multi-Gaussian Schell-model beams in oceanic turbulence[J]. Optik—International Journal for Light and Electron Optics, 2016, 127(8): 3775-3778. doi: 10.1016/j.ijleo.2016.01.047 [17] YUAN Y Sh, LIU X L, WANG F, et al. Scintillation index of a multi-Gaussian Schell-model beam in turbulent atmosphere[J]. Optics Communications, 2013, 305: 57-65. doi: 10.1016/j.optcom.2013.04.076 [18] AVRAMOV-ZAMUROVIC S, NELSON C, GUTH S, et al. Flatness parameter influence on scintillation reduction for multi-Gaussian Schell-model beams propagating in turbulent air[J]. Applied Optics, 2016, 55(13): 3442-3446. doi: 10.1364/AO.55.003442 [19] LIU X Y, ZHAO D M. Optical trapping Rayleigh particles by using focused multi-Gaussian Schell-model beams[J]. Applied Optics, 2014, 53(18): 3976-3981. doi: 10.1364/AO.53.003976 [20] BORGHI R, GORI F, GUATTARI G, et al. Twisted Schell-model beams with axial symmetry[J]. Optics Letters, 2015, 40(19): 4504-4507. doi: 10.1364/OL.40.004504 [21] GORI F, SANTARSIERO M. Devising genuine twisted cross-spectral densities[J]. Optics Letters, 2018, 43(3): 595-598. doi: 10.1364/OL.43.000595 [22] WU G F. Propagation properties of a radially polarized partially coherent twisted beam in free space[J]. Journal of the Optical Society of America, 2016, A33(3): 345-350. [23] PENG X F, LIU L, YU J Y, et al. Propagation of a radially polarized twisted Gaussian Schell-model beam in turbulent atmosphere[J]. Journal of Optics, 2016, 18(12): 125601. doi: 10.1088/2040-8978/18/12/125601 [24] FU W Y, CAO P F. Second-order statistics of a radially polarized partially coherent twisted beam in a uniaxial crystal[J]. Journal of the Optical Society of America, 2017, A34(9): 1703-1710. [25] WAN L P, ZHAO D M. Twisted Gaussian Schell-model array beams[J]. Optics Letters, 2018, 43(15): 3554-3557. doi: 10.1364/OL.43.003554 [26] ZHANG J P, WANG J, HUANG H K, et al. Propagation characteristics of a twisted Cosine-Gaussian correlated radially polarized beam[J]. Applied Sciences, 2018, 8(9): 1485-1498. doi: 10.3390/app8091485 [27] WANG J, HUANG H K, CHEN Y K, et al. Twisted partially coherent array sources and their transmission in anisotropic turbulence[J]. Optics Express, 2018, 26(20): 25974-25988. doi: 10.1364/OE.26.025974 [28] STAHL C S D, GBUR G. Twisted vortex Gaussian Schell-model beams[J]. Journal of the Optical Society of America, 2018, A35(11): 1899-1906. [29] PENG X F, LIU L, WANG F, et al. Twisted Laguerre-Gaussian Schell-model beam and its orbital angular moment[J]. Optics Express, 2018, 26(26): 33956-33969. doi: 10.1364/OE.26.033956 [30] WANG H Y, PENG X F, LIU L, et al. Twisted elliptical multi-Gaussian Schell-model beams and their propagation properties[J]. Journal of the Optical Society of America, 2020, A37(1): 89-97. [31] SONG S H, PARK S, OH C H, et al. Gradient-index planar optics for optical interconnections[J]. Optics Letters, 1998, 23(13): 1025-1027. doi: 10.1364/OL.23.001025 [32] PONOMARENKO S A. Self-imaging of partially coherent light in graded-index media[J]. Optics Letters, 2015, 40(4): 566-568. doi: 10.1364/OL.40.000566 [33] KUNKEL W M, LEGER J R. Gradient-index design for mode conversion of diffracting beams[J]. Optics Express, 2016, 24(12): 13480-13488. doi: 10.1364/OE.24.013480 [34] ROYCHOWDHURY H, AGRAWAL G P, WOLF E. Changes in the spectrum, in the spectral degree of polarization, and in the spectral degree of coherence of a partially coherent beam propagating through a gradient-index fiber[J]. Journal of the Optical Society of America, 2006, A23(4): 940-948. [35] ZHU Sh J, LIU L, CHEN Y H, et al. State of polarization and propagation factor of a stochastic electromagnetic beam in a gradient-index fiber[J]. Journal of the Optical Society of America, 2013, A30(11): 2306-2313. [36] SAHIN S, ZHANG M H, CHEN Y H, et al. Transmission of a polychromatic electromagnetic multi-Gaussian Schell-model beam in an inhomogeneous gradient-index fiber[J]. Journal of the Optical Society of America, 2018, A35(9): 1604-1611. [37] WANG J, YANG Sh, GUO M J, et al. Change in phase singularities of a partially coherent Gaussian vortex beam propagating in a GRIN fiber[J]. Optics Express, 2020, 28(4): 4661-4673. doi: 10.1364/OE.386167 [38] GAO Y R, GAO X Y, LIN R, et al. Periodic properties of Laguerre-Gaussian correlated Schell-model beams in a gradient-index fiber[J]. Optik—International Journal for Light and Electron Optics, 2021, 228: 165755. doi: 10.1016/j.ijleo.2020.165755 [39] WOLF E. New theory of partial coherence in the space-frequency domain Part Ⅰ: Spectra and cross spectra of steady-state sources[J]. Journal of the Optical Society of America, 1982, 72(3): 343-351. doi: 10.1364/JOSA.72.000343 -
点击查看大图
图(6)
计量
- 文章访问数: 3563
- HTML全文浏览量: 2556
- PDF下载量: 24
- 被引次数: 0
