Diffusion GYACOMO (OBSOLETE: LATEST VERSION ON gitlab.epfl.ch) (Miller) Browse

gyacomo/wk/g2k3_transport_scaling.m
fa01fb8b1798Miller

g2k3_transport_scaling.m
View Options

			%% G ~ g^2/k^3 rule

			K_N = [1.60 1.70 1.80 1.90 2.00 2.10 2.20 2.30 2.40 2.50];
			%DG
			g_DG = [0.08 0.20 0.31 0.42 0.55 0.61 0.70 0.79 0.87 0.96];
			k_DG = [0.21 0.31 0.37 0.42 0.47 0.47 0.47 0.52 0.52 0.47];
			G_DG = g_DG.^2./k_DG.^3./K_N;
			%SG
			g_SG = [0.08 0.19 0.30 0.40 0.49 0.59 0.68 0.76 0.85 0.93];
			k_SG = [0.26 0.37 0.42 0.42 0.47 0.47 0.52 0.52 0.52 0.47];
			G_SG = g_SG.^2./k_SG.^3./K_N;
			%LR
			g_LR = [0.12 0.18 0.25 0.32 0.39 0.47 0.55 0.64 0.72 0.81];
			k_LR = [0.31 0.31 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37];
			G_LR = g_LR.^2./k_LR.^3./K_N;
			%LD
			g_LD = [0.05 0.11 0.19 0.27 0.35 0.44 0.54 0.63 0.73 0.83];
			k_LD = [0.26 0.26 0.31 0.31 0.31 0.31 0.31 0.31 0.31 0.31];
			G_LD = g_LD.^2./k_LD.^3./K_N;
			% nu=0
			g_CL = [0.22 0.27 0.32 0.37 0.44 0.51 0.59 0.67 0.76 0.85];
			k_CL = [0.73 0.73 0.73 0.73 0.63 0.52 0.52 0.52 0.52 0.42];
			G_CL = g_CL.^2./k_CL.^3./K_N;
			%% Plot

			figure;
			% plot(K_N,G_SG); hold on;
			% plot(K_N,G_DG); hold on;
			% plot(K_N,G_LR); hold on;
			% plot(K_N,G_LD); hold on;
			% plot(K_N,G_CL,'--k'); hold on;
			k = 0.35;
			plot(K_N,g_SG.^2/k.^3); hold on;
			plot(K_N,g_DG.^2/k.^3); hold on;
			plot(K_N,g_LR.^2/k.^3); hold on;
			plot(K_N,g_LD.^2/k.^3); hold on;
			plot(K_N,g_CL.^2/k.^3,'--k'); hold on;
			% plot(K_N,(G_DG+G_SG+G_LR+G_LD)/4);