{"id":235766,"date":"2022-09-20T08:37:20","date_gmt":"2022-09-20T08:37:20","guid":{"rendered":"\/forum\/forums\/topic\/how-to-set-the-oam-parameters-of-incident-light-in-lumerical-fdtd\/"},"modified":"2022-09-20T08:37:20","modified_gmt":"2022-09-20T08:37:20","slug":"how-to-set-the-oam-parameters-of-incident-light-in-lumerical-fdtd","status":"closed","type":"topic","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/how-to-set-the-oam-parameters-of-incident-light-in-lumerical-fdtd\/","title":{"rendered":"How to set the OAM parameters of incident light in Lumerical FDTD?"},"content":{"rendered":"

Hi Everyone,<\/p>\n

Right now, I want to set up the orbital angular momentum (OAM) modes of the incident light in Lumerical FDTD. Does anyone know how to set the OAM mode of incident light?<\/p>\n

It looks like that we can’t set up this in the software itself, but we might be able to set it up in the lumerical script section. Below is an example of setting some parameters of the incident light, which I got online. It might a hint\/reference to you.<\/p>\n

Really hope someone can tackle this issue. Thanks.<\/p>\n

# scriptfile:usr_custom_source.lsf<\/div>\n
 <\/div>\n
# define position vectors<\/div>\n
 <\/div>\n
x=linspace(-10e-6,10e-6,401);<\/div>\n
y=linspace(-10e-6,10e-6,401);<\/div>\n
z=0; #this is z z-normal source<\/div>\n
 <\/div>\n
X=meshgridx(x,y);<\/div>\n
Y=meshgridx(x,y);<\/div>\n
 <\/div>\n
lambda0=0.5e-6;<\/div>\n
f=c\/lambda0;<\/div>\n
w=2*pi*f;<\/div>\n
k=2*pi\/lambda0;<\/div>\n
 <\/div>\n
#calculate k-space field profiles<\/div>\n
NA=0.2;<\/div>\n
 <\/div>\n
kx=linspace(-k,k,200);<\/div>\n
ky=linspace(-k,k,200);<\/div>\n
Kx=meshgridx(kx,ky);<\/div>\n
Ky=meshgridy(kx,ky);<\/div>\n
phi=atan2(Ky,Kx);<\/div>\n
theta=real(acos(sqrt(1-Kx^2\/k^2-Ky^2\/k^2)));<\/div>\n
envelope=exp(-0.5*(Kx^2+Ky^2)\/(NA*k)^2);<\/div>\n
 <\/div>\n
Exk=cos(phi)*cos(theta)*envelope;<\/div>\n
 <\/div>\n
Eyk=sin(phi)*cos(theta)*envelope;<\/div>\n
 <\/div>\n
Ezk=sin(theta)*envelope;<\/div>\n
 <\/div>\n
Ex=-1i*czt(Exk,kx,ky,x,y);<\/div>\n
 <\/div>\n
Ey=-1i*czt(Eyk,kx,ky,x,y);<\/div>\n
 <\/div>\n
Ez=-1i*czt(Ezk,kx,ky,x,y);<\/div>\n
 <\/div>\n
Hxk=-sqrt(eps0\/mu0)*sin(phi)*envelope;<\/div>\n
Hyk=sqrt(eps0\/mu0)*cos(phi)*envelope;<\/div>\n
Hx=-1i*czt(Hxk,kx,ky,x,y);<\/div>\n
Hy=-1i*czt(Hyk,kx,ky,x,y);<\/div>\n
Hz=0*Ex;<\/div>\n
 <\/div>\n
# sacle fild so E^2=1<\/div>\n
 <\/div>\n
E2=abs(Ex)^2+abs(Ey)^2+abs(Ez)^2;<\/div>\n
scaleFactor=sqrt(max(E2));<\/div>\n
Ex=Ex\/scaleFactor;<\/div>\n
Ey=Ey\/scaleFactor;<\/div>\n
Ez=Ez\/scaleFactor;<\/div>\n
Hx=Hx\/scaleFactor;<\/div>\n
Hy=Hy\/scaleFactor;<\/div>\n
Hz=Hz\/scaleFactor;<\/div>\n
 <\/div>\n
# package field data into the EM fields dataset<\/div>\n
 <\/div>\n
EM=rectilineardataset(“EM fields”,x,y,z);<\/div>\n
EM.addparameter(“lambda”,c\/f,”f”,f);<\/div>\n
EM.addattribute(“E”,Ex,Ey,Ez);<\/div>\n
EM.addattribute(“H”,Hx,Hy,Hz);<\/div>\n
 <\/div>\n
# save dataset to .mat file<\/div>\n
matlabsave(“import_data.mat”,EM);<\/div>\n","protected":false},"template":"","class_list":["post-235766","topic","type-topic","status-closed","hentry","topic-tag-define-source","topic-tag-LumericalFDTD-1","topic-tag-photonics","topic-tag-source-setting-1"],"aioseo_notices":[],"aioseo_head":"\n\t\t\n\t\n\t\n\t\n\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t