gusucode.com > PPML_v1.2 > examples/UniversalLineshapes_SciRep2016_fig3_c3.m

    close all;
clear all
          
% Calculates and plots the double-peaked spectrum (solid points) 
% 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 = 201;                  % number of points in energy
om = linspace(100,200,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 = 3e11;     % 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 


% main cycle over photon energy
for i = 1:size(om,2)
    i

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     = [0      0.5      0.05          1.3      0]; % layer thicknesses
epsxA = [       epsbg    epsbg         epsbg     ]; 
epszA = [       epsbg    epsbg         epszn(i)  ];   
epsxB = [       epsbg    epsAu         epsbg     ];     
epszB = [       epsbg    epsAu         epszn(i)  ];   

% calling PPML function "RTA_1d_tm"
[RR,TT,AA] = RTA_1d_tm(a,L,...
   epssup,epssub,epsxA,epszA,epsxB,epszB,f,d,...
   halfnpw,k0(i),k0(i)*sin(theta*pi/180));
T(i) = TT; % transmitted intensity
end

figure(1)
plot(om,T,'k.')