OAM轨道角动量的产生 matlab

动量是物理学中一个基本物理量，结构光除了携带线性动量，也可以携带角动量。角动 量又分为自旋角动量（spin angular momentum, SAM）和轨道角动量（OAM）。利用 OAM 光 束的空间维度作为信息载体，可以独立的传输信息。1992 年 Allen 等人通过实验得出 OAM 由具有螺旋相位的波前光束携带，相位因子为 exp(1i*l*phi ) ，l 代表螺旋波阵面的螺旋模式指 数即拓扑荷(topological charge)，phi代表方位角，涡旋光的相位呈螺旋分布，中心处的电场强 度为零，因此涡旋光束的分布在横截面上呈“甜甜圈”状。由于无限螺旋模式数和不同阶数 之间的正交性使得 OAM 成为一种备受关注的光学自由度，用于提高光学和量子的信息 容量，多个 OAM 模式可以实现复用/解复用，将 OAM 全息与高分辨率超表面结合实现数字 全息，已经成为目前研究的热点，为超高容量全息设备和系统开辟新的路径。

利用OAM的第一步是要学会生成涡旋光，而在实现OAM全息的过程中，不需要考虑光的强度，只需要利用相位，但采样的流程中需要考虑到OAM的强度对采样间隔的影响。

一般利用LG光束

实现的matlab代码如下：

N = 200;
lambda = 632e-9;    %波长为632nm
k = 2*pi/lambda;    %波数
w0 = 3e-3;          %束腰尺寸
x = linspace(-3*w0,3*w0,N);     y = x;
[X,Y] = meshgrid(x,y);
[theta,r] = cart2pol(X,Y);
Z_R = pi*w0^2/lambda;      %瑞利长度
z = 0;
w_z = w0*sqrt(1+(z/Z_R)^2);%光束在z位置的半径
figure;
p = 2; 
m=1；        %携带拓扑荷为1
 E3 = sqrt(2*factorial(p)/pi/(p+factorial(abs(m))))*(1/w_z)*(sqrt(2)*r/w_z).^abs(m).*exp(-r.^2/w_z^2).*laguerre(p,abs(m),2*r.^2/w_z^2).*exp(-1i*m*theta).*exp(-1i*k*z).*exp(-1i*k*r.^2*z/2/(z^2+Z_R^2))*exp(-1i*(2*p+abs(m)+1)*atan(z/Z_R));
    I3 = E3.*conj(E3);  I3 = I3/max(max(I3));

效果如下 所示：可根据需要获得所需拓扑荷的强度分布E3