gusucode.com > PPML_v1.2 > examples/UniversalLineshapes_SciRep2016_fig3_c2.m
close all; clear all % Calculates and plots the E field % in Fig. 3c of "Universal lineshapes at the crossover between weak % and strong critical coupling in Fano-resonant % coupled oscillators" Sci. Rep. 2016; 6: 24592 % % Simone Zanotto, Firenze, feb. 2016 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % This is free software distributed under the BSD licence (see the % containing folder). % However, shall the results obtained through this code be included % in an academic publication, we kindly ask you to cite the source % website and, if applicable, the following paper: % % Simone Zanotto, Alessandro Tredicucci, "Universal lineshapes % at the crossover between weak and strong critical coupling % in Fano-resonant coupled oscillators", Sci. Rep. 2016; 6: 24592 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% addpath('PPML_root') % to be replaced by the proper path numom = 1; % number of points in energy om = linspace(150,150,numom); % photon energy (hbar omega) in meV k0 = om*(2*pi/1240); % vacuum wavenumber in inverse micron theta = 0.1; % inc. angle. NEVER SET ZERO. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Some structure parameters a = 3.2; % lattice spacing halfnpw = 10; % plane wave truncation. % WARNING! this is a computation-dependent parameter. % inaccurate simulation may result from wrong settings % material properties epsAu = -4000 + 300*1i; % Gold at 10 um %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% MQW parameters %%%% see Appl. Phys. Lett. 105, 081105 (2014) %%%% http://dx.doi.org/10.1063/1.4893730 epsbg = 10.05; % GaAs at 10 um n2deg = 0e11; % doping in cm-2 omp = sqrt(4*pi*4.8e-10^2*n2deg/(epsbg*0.067*9.1e-28*21.8e-7))*6.58e-13; % % eff. MQW bulk plasma freq, meV. :: 21.8 = Lb + Lw (nm) om12 = 150; % meV gamma12 = 3; % meV HWHM epszn = epsbg./(1-omp^2./(om12^2-om.^2-1i*2*gamma12*om)); % see i = 1; epssup = 1; epssub = 1; % super- and sub-strate permeabilities (air) L = 3; % number of internal layers % GaAs | metal pett. | GaAs f = [ 0 0.73 0 ]; % internal layer fill fraction (fraction of material B in A) d = [2 0.5 0.05 1.3 2]; % layer thicknesses epsxA = [ epsbg epsbg epsbg ]; epszA = [ epsbg epsbg epszn(i) ]; epsxB = [ epsbg epsAu epsbg ]; epszB = [ epsbg epsAu epszn(i) ]; nx = 100; nz = [50 50 10 50 50]; % calling PPML function "field_1d_tm" [x,z,Ex,Ez,Sz] = field_1d_tm(a,L,... epssup,epssub,epsxA,epszA,epsxB,epszB,f,d,... halfnpw,k0,k0*sin(theta*pi/180),nx,nz); figure(1) h = pcolor(x,z,(abs(Ez').^2)+(abs(Ex').^2)); hold on; set(h,'edgecolor','none')