论文复现|tightly focused circularly polarized ring Airy beam

发布时间:2024年01月16日

请尊重原创的劳动成果
如需要转载,请后台联系

前言

采用MATLAB复现一篇论文里面的插图,涡旋光束的聚焦的仿真方式有很多种,这里采用MATLAB进行仿真,当然也有其他的很多方式,不同的方式各有千秋。

论文摘要

本文证明了在紧聚焦透镜作用下聚焦圆偏振环艾里光束(RAB)的双聚焦特性。计算了紧聚焦环空光束两个焦点处的辐射力,数值结果表明粒子可能在两个焦点处被纵向和横向捕获。通过改变相应的参数,我们可以控制两个陷波的特性。捕获力随 NA 和缩放参数 w 的增大而增大;适当的初始半径 r0 是增强任一捕获的必要条件。随着 w 的增大或 r0 的减小,两个陷阱可以靠得更近。为了实现双光学陷阱,我们应该选择较小的衰减参数 a 和较大的 NA,否则双光学陷阱会退化为单光陷阱。此外,由于布朗运动和散射力的影响,粒子的尺寸应在一个特殊的范围内。

复现内容

在这里插入图片描述
在这里插入图片描述
上面的第一张图是论文里面的插图,下面两张图是采用MATLAB和论文中给出的公式进行编写复现的,可以看到,复现的结果非常好。
此外还修改了部分参数进行进一步的分析:
在这里插入图片描述



Ep = k*f*exp(1i*phi0)/(2*sqrt(2)).*sum(E1,3);
Er = -1i*k*f*exp(1i*phi0)/(2*sqrt(2)).*sum(E2,3);
Ez = k*f*exp(1i*phi0)/(sqrt(2)).*sum(E3,3);

I = Ep.*conj(Ep)+Er.*conj(Er)+Ez.*conj(Ez);
normalizeI = I/max(max(I));
% writematrix(normalizeI,'normalizeI.xls');
I_Ex = Ep.*conj(Ep)/max(max(I));
I_Ey = Er.*conj(Er)/max(max(I));
I_Ez = Ez.*conj(Ez)/max(max(I));

%% =============================================
% 绘图操作
figure('Name','vortex beam','Position',[471 326 1058 652])
pcolor(z(:,:,1),r(:,:,1),I);
set(gca,'fontname','times new roman','fontsize',18, ...
    'XGrid','off','YGrid','off','GridLineStyle','-.','XMinorTick','off','YMinorTick','off','XMinorGrid','off','YMinorGrid','off', ...
    'TickDir','none','Box','off');
set(gca,'XTick',1.98e-3:0.01e-3:2.02e-3,'XTicklabel',[1.98:0.01:2.02]);  %修改坐标轴的范围
set(gca,'YTick',-10e-6:5e-6:10e-6,'YTicklabel',[-10:5:10]);   %修改坐标轴的范围
xlabel('z(mm)');
ylabel('r(\mum)');
% caxis([0,1]);
shading interp;
colormap jet;
ck = colorbar("location",'northoutside','Ticks',0:0.5e14:2e14,'TickLabels',[0:0.5:2]);
set(get(ck,'title'),'string','Intensity(10^{14}V^{2}/m^{2})','fontsize',15);

上面只是展示了复现的部分代码,如果需要全部完整的代码,请直接在微信公众号搜索:“XD悟理”,找到第六期推文就可以拿到完整的代码。不要任何其他操作

文章来源:https://blog.csdn.net/weixin_44449744/article/details/135614125
本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。