Gaussian beam incident a Negative refractive index plate

中文/English

Reference


1. 'Left-Handed' Materials Propel Spectrum Physics


2. Lumerical FDTD Solution


3.Left-Hand material , Metalmaterial 。

J. B. Pendry, “Negative Refraction Makes a Perfect Lens,” Phys. Rev. Lett., vol. 85, no. 18, pp. 3966–3969, Oct. 2000


4.Negative index metamaterials


To achieve a negative refractive index, we can excite a Gaussian beam incident into a double negative refractive (DNR) materials (permeability = -1 & permittivity = -1) . Here we can use the "Drude model"


Simulation parameters:

DNR plate thickness = 600 nm

Incident wavelength = 200 nm

Gaussian beam size = 200 nm




Define Parameters

1. Selected the simulation dimension.

2. Define total x, y, z length that equal to CAD, and then set a appropriate resolution.

3. Press【建立網格 Cad to Grids (Create)】button to create simulation grids.





A. Define the boundary conditions

B. Press the【創建 (Create)】button => Create the total grid size (Include the boundary condition & add space)


The oblique incident wave is linear polarization from bottom side(theta =45).





1. Checked the simulation material that will be used

2. Input the builtin function

3. Press the run button to execute the code

4. Checked out the geometries

5. If correct! then Press 【輸出 output】button to output the *.csv files.


the builtin function geometries variable are

gdx is (grid) dx value, gdy = dy. gdz = dz;

ib= X grids, jb= Y grids, kb= Z grids;

icenter= ib/2; jcenter= jb/2; kcenter=kb/2;



%=================矩形 Brick===========================
xstart=1*gdx;
xend=ib*gdx;
ystart=1*gdy;
yend=1*gdy;
zstart=kcenter*gdz-300e-9;
zend=kcenter*gdz+300e-9;
nindex=1^2;
sigma=0;
choice='EM_Model1'; %E_Iso,PEC,M_Iso,PMC,E_Model1,M_Model1,EM_Model1
gridtype=-1; %
Iso_Brick(choice,gridtype,nindex,sigma,xstart,xend,ystart,yend,zstart,zend)
%=================矩形 Brick===========================


double negative refractive (DNR) materials (permeability = -1 & permittivity = -1)

To use the Drude model, if we define epsilon_infinity=1, , gd=0, then epsilonr=1-wd^2/w^2; and then obtain wd=sqrt(2)*w

if let the incident wavelength = 200 nm, can get wd=sqrt(2)*pi/3e8/200e-9 = 1.332864881447510e+016。

take this value into both E-Models & M_Models, Using the Drude model.





1. Set the Spectrum Analysis & Wavelength range




1. Set the simulation time setp

2. Save parameters

3. Calculation


★(Result Analysis):