added line graphs

2D-Conductivity/conductivity_calc_surface.m

@@ -32,46 +32,46 @@ f_vals = f_vals .* (2*pi); % rads-1
 
 % Carrier Density
 %%%%%%%
-% carrier_vals = logspace(0, MAX_Y, Y_TOTAL); % m-2
-% 
-% % below turns turns carrier densities into Fermi energies
-% fermi_vals = zeros(1, length(carrier_vals));
-% for carr=1:length(carrier_vals)
-%     fermi_vals(carr) = fermi_from_carrier_density(carrier_vals(carr), ev_to_j(t));
-% end
-% 
-% % CALCULATE SHEET CONDUCTIVITY
-% cond = zeros(length(f_vals),... % frequency
-%              length(fermi_vals),... % fermi
-%              2); % intra/inter
-% for freq=1:length(f_vals)
-%     for y=1:length(fermi_vals)
-%         
-%         cond(freq, y, :) = sheet_conductivity(f_vals(freq),... % omega (rads-1)
-%                                               fermi_vals(y),... % fermi_level (J)
-%                                               300,... % temp (K)
-%                                               5e-12); % scatter_lifetime (s)
-%     end
-% end
+carrier_vals = logspace(0, MAX_Y, Y_TOTAL); % m-2
 
-% Temperature
-%%%%%%%
-temp_vals = linspace(0, 2230, Y_TOTAL); % K
+% below turns turns carrier densities into Fermi energies
+fermi_vals = zeros(1, length(carrier_vals));
+for carr=1:length(carrier_vals)
+    fermi_vals(carr) = fermi_from_carrier_density(carrier_vals(carr), ev_to_j(t));
+end
 
 % CALCULATE SHEET CONDUCTIVITY
 cond = zeros(length(f_vals),... % frequency
-             length(temp_vals),... % fermi
+             length(fermi_vals),... % fermi
              2); % intra/inter
 for freq=1:length(f_vals)
-    for y=1:length(temp_vals)
+    for y=1:length(fermi_vals)
         
         cond(freq, y, :) = sheet_conductivity(f_vals(freq),... % omega (rads-1)
-                                              fermi_from_carrier_density(1.3e13*10000, ev_to_j(t)),... % fermi_level (J)
-                                              temp_vals(y),... % temp (K)
+                                              fermi_vals(y),... % fermi_level (J)
+                                              300,... % temp (K)
                                               5e-12); % scatter_lifetime (s)
     end
 end
 
+% Temperature
+%%%%%%%
+% temp_vals = linspace(0, 2230, Y_TOTAL); % K
+% 
+% CALCULATE SHEET CONDUCTIVITY
+% cond = zeros(length(f_vals),... % frequency
+%              length(temp_vals),... % fermi
+%              2); % intra/inter
+% for freq=1:length(f_vals)
+%     for y=1:length(temp_vals)
+%         
+%         cond(freq, y, :) = sheet_conductivity(f_vals(freq),... % omega (rads-1)
+%                                               fermi_from_carrier_density(1.3e13*10000, ev_to_j(t)),... % fermi_level (J)
+%                                               temp_vals(y),... % temp (K)
+%                                               5e-12); % scatter_lifetime (s)
+%     end
+% end
+
 % Scatter Lifetime
 %%%%%%%
 % scatt_vals = logspace(-11, -14, Y_TOTAL); % s-1
@@ -98,7 +98,7 @@ if DISPLAY_HZ % divide radians back to hertz
     f_vals = f_vals ./ (2*pi);
 end
 
-y_vals = temp_vals;
+y_vals = carrier_vals;
 
 % cond = sign(cond).*log10(abs(cond));
 
@@ -119,13 +119,13 @@ grid;
 colorbar;
 axis tight;
 set(gca, 'xscale', 'log')
-% set(gca, 'yscale', 'log')
+set(gca, 'yscale', 'log')
 % set(gca, 'zscale', 'log')
 
 set(gca, 'ColorScale', 'log')
 
-% ylabel('Net Carrier Density (m^{-2})');
-ylabel('Temperature (K)');
+ylabel('Net Carrier Density (m^{-2})');
+% ylabel('Temperature (K)');
 % ylabel('Scatter Lifetime (s)');
 
 zlabel('Conductivity (S)');
@@ -152,13 +152,13 @@ grid;
 colorbar;
 axis tight;
 set(gca, 'xscale', 'log')
-% set(gca, 'yscale', 'log')
+set(gca, 'yscale', 'log')
 % set(gca, 'zscale', 'log')
 
-% set(gca, 'ColorScale', 'log')
+set(gca, 'ColorScale', 'log')
 
-% ylabel('Net Carrier Density (m^{-2})');
-ylabel('Temperature (K)');
+ylabel('Net Carrier Density (m^{-2})');
+% ylabel('Temperature (K)');
 % ylabel('Scatter Lifetime (s)');
 
 zlabel('Conductivity (S)');

2D-Conductivity/conductivity_calculations.m

@@ -18,7 +18,7 @@ F_TOTAL = 1e2; % number of points to generate
 % EXCITATION_TYPE = 'inter';
 EXCITATION_TYPE = 'all';
 
-TWO_SERIES = true; % for comparing two dopants
+MULTIPLE_SERIES = true; % for comparing two dopants
 
 %%%%%%%%%%%%%%%%%%%%%%%%%%%
 %%     CALCULATE
@@ -31,19 +31,27 @@ x_vals = x_vals .* (2*pi); % rads-1
 cond = zeros(length(x_vals), 2);
 for x=1:length(x_vals)
     cond(x, :) = sheet_conductivity(x_vals(x),... % omega (rads-1)
-                                    fermi_from_carrier_density(1.3e13*10000, ev_to_j(3)),... % fermi_level (J)
+                                    fermi_from_carrier_density(1.3e13*100*100, ev_to_j(3)),... % fermi_level (J)
                                     300,... % temp (K)
                                     1e-12); % scatter_lifetime (s)
 end
 
-if TWO_SERIES
+if MULTIPLE_SERIES
     cond2 = zeros(length(x_vals), 2);
     for x=1:length(x_vals)
         cond2(x, :) = sheet_conductivity(x_vals(x),... % omega (rads-1)
-                                        fermi_from_carrier_density(2.2e13*10000, ev_to_j(3)),... % fermi_level (J)
+                                        fermi_from_carrier_density(2.2e13*100*100, ev_to_j(3)),... % fermi_level (J)
                                         300,... % temp (K)
                                         1e-12); % scatter_lifetime (s)
     end
+    
+%     cond3 = zeros(length(x_vals), 2);
+%     for x=1:length(x_vals)
+%         cond3(x, :) = sheet_conductivity(x_vals(x),... % omega (rads-1)
+%                                         fermi_from_carrier_density(1.3e13*10000, ev_to_j(3)),... % fermi_level (J)
+%                                         300,... % temp (K)
+%                                         1e-12); % scatter_lifetime (s)
+%     end
 end
 
 if DISPLAY_HZ % divide radians back to hertz
@@ -54,11 +62,16 @@ end
 %%       RENDER
 %%%%%%%%%%%%%%%%%%%%%%%%%%%
 
-RE_COLOUR = 'r';
+RE_COLOUR = 'r-';
 IM_COLOUR = 'r--';
-RE_COLOUR2 = 'b';
+MAG_COLOUR = 'r:';
+RE_COLOUR2 = 'b-';
 IM_COLOUR2 = 'b--';
-LW = 1.5;
+MAG_COLOUR2 = 'b:';
+RE_COLOUR3 = 'b';
+IM_COLOUR3 = 'b--';
+MAG_COLOUR3 = 'b:';
+LW = 2;
 
 figure(1);
 hold on;
@@ -66,10 +79,16 @@ hold on;
 if strcmp(EXCITATION_TYPE, 'intra')
     plot(x_vals, real(cond(:, 1)), RE_COLOUR, 'LineWidth', LW);
     plot(x_vals, imag(cond(:, 1)), IM_COLOUR, 'LineWidth', LW);
+    plot(x_vals, abs(cond(:, 1)), MAG_COLOUR, 'LineWidth', LW);
     
-    if TWO_SERIES
+    if MULTIPLE_SERIES
         plot(x_vals, real(cond2(:, 1)), RE_COLOUR2, 'LineWidth', LW);
         plot(x_vals, imag(cond2(:, 1)), IM_COLOUR2, 'LineWidth', LW);
+        plot(x_vals, abs(cond2(:, 1)), MAG_COLOUR2, 'LineWidth', LW);
+        
+%         plot(x_vals, real(cond3(:, 1)), RE_COLOUR3, 'LineWidth', LW);
+%         plot(x_vals, imag(cond3(:, 1)), IM_COLOUR3, 'LineWidth', LW);
+%         plot(x_vals, abs(cond3(:, 1)), MAG_COLOUR3, 'LineWidth', LW);
     end
     title('2D Intraband Sheet Conductivity');
 
@@ -77,10 +96,16 @@ if strcmp(EXCITATION_TYPE, 'intra')
 elseif strcmp(EXCITATION_TYPE, 'inter')
     plot(x_vals, real(cond(:, 2)), RE_COLOUR, 'LineWidth', LW);
     plot(x_vals, imag(cond(:, 2)), IM_COLOUR, 'LineWidth', LW);
+    plot(x_vals, abs(cond(:, 2)), MAG_COLOUR, 'LineWidth', LW);
     
-    if TWO_SERIES
+    if MULTIPLE_SERIES
         plot(x_vals, real(cond2(:, 2)), RE_COLOUR2, 'LineWidth', LW);
         plot(x_vals, imag(cond2(:, 2)), IM_COLOUR2, 'LineWidth', LW);
+        plot(x_vals, abs(cond2(:, 2)), MAG_COLOUR2, 'LineWidth', LW);
+        
+%         plot(x_vals, real(cond3(:, 2)), RE_COLOUR3, 'LineWidth', LW);
+%         plot(x_vals, imag(cond3(:, 2)), IM_COLOUR3, 'LineWidth', LW);
+%         plot(x_vals, abs(cond3(:, 2)), MAG_COLOUR3, 'LineWidth', LW);
     end
     title('2D Interband Sheet Conductivity');
     
@@ -88,10 +113,16 @@ elseif strcmp(EXCITATION_TYPE, 'inter')
 else
     plot(x_vals, real(sum(cond, 2)), RE_COLOUR, 'LineWidth', LW);
     plot(x_vals, imag(sum(cond, 2)), IM_COLOUR, 'LineWidth', LW);
+    plot(x_vals, abs(sum(cond, 2)), MAG_COLOUR, 'LineWidth', LW);
     
-    if TWO_SERIES
+    if MULTIPLE_SERIES
         plot(x_vals, real(sum(cond2, 2)), RE_COLOUR2, 'LineWidth', LW);
         plot(x_vals, imag(sum(cond2, 2)), IM_COLOUR2, 'LineWidth', LW);
+        plot(x_vals, abs(sum(cond2, 2)), MAG_COLOUR2, 'LineWidth', LW);
+        
+%         plot(x_vals, real(sum(cond3, 2)), RE_COLOUR3, 'LineWidth', LW);
+%         plot(x_vals, imag(sum(cond3, 2)), IM_COLOUR3, 'LineWidth', LW);
+%         plot(x_vals, abs(sum(cond3, 2)), MAG_COLOUR3, 'LineWidth', LW);
     end
     title('2D Sheet Conductivity');
 end
@@ -100,8 +131,9 @@ set(gca,'Xscale','log')
 % set(gca,'Yscale','log')
 axis tight
 
-if TWO_SERIES
-    legend('Real TTF', 'Imaginary TTF', 'Real CoCp2', 'Imaginary CoCp2');
+if MULTIPLE_SERIES
+    legend('Re(TTF)', 'Im(TTF)', '|TTF|', 'Re(CoCp_2)', 'Im(CoCp_2)', '|CoCp_2|');
+%     legend('Re(1x10^{8}m^{-2})', 'Im(1x10^{8}m^{-2})', '|1x10^{8}m^{-2}|', 'Re(1x10^{15}m^{-2})', 'Im(1x10^{15}m^{-2})', '|1x10^{15}m^{-2}|', 'Re(1.3x10^{17}m^{-2})', 'Im(1.3x10^{17}m^{-2})', '|1.3x10^{17}m^{-2}|');
 else
     legend('Real', 'Imaginary');
 end