h = figure(1);
set(h,'Position',[100 678 1200 360]);
index = 0;
prefix = 'Data/TwoPlatesLoose_Charged_L25_index';
prefix_PMF = 'Data/TwoPlatesTight_Charged_PMF_L25_index';
fname = sprintf([prefix, '%1.1d.mat'], index);
fname_PMF = sprintf([prefix_PMF, '%1.1d.mat'], index);
load(fname);
geom_d = energy(1,:);
VdW_d = energy(2,:);
elec_d = energy(3,:);
load(fname_PMF);
LJ_pair = energy(1,:);
elec_pair = energy(2,:);
elec_d_outside = energy(3,:);
geom_infty = energy(4,1);
VdW_infty = energy(4,2);
elec_infty = energy(4,3);
geom = geom_d - geom_infty;
VdW = VdW_d - VdW_infty + LJ_pair;
elec = elec_d + elec_d_outside - elec_infty + elec_pair;
total = geom + VdW + elec;
subplot(1,3,1); plot(4:1:24,total,'b-+');hold on;
prefix = 'Data/TwoPlatesTight_Charged_L25_index';
prefix_PMF = 'Data/TwoPlatesTight_Charged_PMF_L25_index';
fname = sprintf([prefix, '%1.1d.mat'], index);
fname_PMF = sprintf([prefix_PMF, '%1.1d.mat'], index);
load(fname);
geom_d = energy(1,:);
VdW_d = energy(2,:);
elec_d = energy(3,:);
load(fname_PMF);
LJ_pair = energy(1,:);
elec_pair = energy(2,:);
elec_d_outside = energy(3,:);
geom_infty = energy(4,1);
VdW_infty = energy(4,2);
elec_infty = energy(4,3);
geom = geom_d - geom_infty;
VdW = VdW_d - VdW_infty + LJ_pair;
elec = elec_d + elec_d_outside - elec_infty + elec_pair;
total = geom + VdW + elec;
subplot(1,3,1); plot(4:1:24,total,'r-o');hold on;
xlim([4,24]); xlabel('d (\AA)','fontsize',20,'interpreter', 'latex');
ylim([-300,100]); ylabel('G_{tot}^{PMF}( k_BT )','fontsize',20,'rot',90);
h_legend = legend('0.0e/0.0e loose','0.0e/0.0e tight');
set(h_legend,'Location','SouthEast','FontSize',18);
set(gca,'xtick',[4,6,8,10,12,14,16,18,20,22,24],'fontsize',15);
index = 3;
prefix = 'Data/TwoPlatesLoose_Charged_L25_index';
prefix_PMF = 'Data/TwoPlatesTight_Charged_PMF_L25_index';
fname = sprintf([prefix, '%1.1d.mat'], index);
fname_PMF = sprintf([prefix_PMF, '%1.1d.mat'], index);
load(fname);
geom_d = energy(1,:);
VdW_d = energy(2,:);
elec_d = energy(3,:);
load(fname_PMF);
LJ_pair = energy(1,:);
elec_pair = energy(2,:);
elec_d_outside = energy(3,:);
geom_infty = energy(4,1);
VdW_infty = energy(4,2);
elec_infty = energy(4,3);
geom = geom_d - geom_infty;
VdW = VdW_d - VdW_infty + LJ_pair;
elec = elec_d + elec_d_outside - elec_infty + elec_pair;
total = geom + VdW + elec;
subplot(1,3,2); plot(4:1:24,total,'b-+');hold on;
prefix = 'Data/TwoPlatesTight_Charged_L25_index';
prefix_PMF = 'Data/TwoPlatesTight_Charged_PMF_L25_index';
fname = sprintf([prefix, '%1.1d.mat'], index);
fname_PMF = sprintf([prefix_PMF, '%1.1d.mat'], index);
load(fname);
geom_d = energy(1,:);
VdW_d = energy(2,:);
elec_d = energy(3,:);
load(fname_PMF);
LJ_pair = energy(1,:);
elec_pair = energy(2,:);
elec_d_outside = energy(3,:);
geom_infty = energy(4,1);
VdW_infty = energy(4,2);
elec_infty = energy(4,3);
geom = geom_d - geom_infty;
VdW = VdW_d - VdW_infty + LJ_pair;
elec = elec_d + elec_d_outside - elec_infty + elec_pair;
total = geom + VdW + elec;
subplot(1,3,2); plot(4:1:24,total,'r-o');hold on;
xlim([4,24]); xlabel('d (\AA)','fontsize',20,'interpreter', 'latex');
ylim([0,350]); ylabel('G_{tot}^{PMF}( k_BT )','fontsize',20,'rot',90);
h_legend = legend('0.2e/0.2e loose','0.2e/0.2e tight');
set(h_legend,'Location','SouthWest','FontSize',18);
set(gca,'xtick',[4,6,8,10,12,14,16,18,20,22,24],'fontsize',15);
index = 4;
prefix = 'Data/TwoPlatesLoose_Charged_L25_index';
prefix_PMF = 'Data/TwoPlatesTight_Charged_PMF_L25_index';
fname = sprintf([prefix, '%1.1d.mat'], index);
fname_PMF = sprintf([prefix_PMF, '%1.1d.mat'], index);
load(fname);
geom_d = energy(1,:);
VdW_d = energy(2,:);
elec_d = energy(3,:);
load(fname_PMF);
LJ_pair = energy(1,:);
elec_pair = energy(2,:);
elec_d_outside = energy(3,:);
geom_infty = energy(4,1);
VdW_infty = energy(4,2);
elec_infty = energy(4,3);
geom = geom_d - geom_infty;
VdW = VdW_d - VdW_infty + LJ_pair;
elec = elec_d + elec_d_outside - elec_infty + elec_pair;
total = geom + VdW + elec;
subplot(1,3,3); plot(4:1:24,total,'b-+');hold on;
prefix = 'Data/TwoPlatesTight_Charged_L25_index';
prefix_PMF = 'Data/TwoPlatesTight_Charged_PMF_L25_index';
fname = sprintf([prefix, '%1.1d.mat'], index);
fname_PMF = sprintf([prefix_PMF, '%1.1d.mat'], index);
load(fname);
geom_d = energy(1,:);
VdW_d = energy(2,:);
elec_d = energy(3,:);
load(fname_PMF);
LJ_pair = energy(1,:);
elec_pair = energy(2,:);
elec_d_outside = energy(3,:);
geom_infty = energy(4,1);
VdW_infty = energy(4,2);
elec_infty = energy(4,3);
geom = geom_d - geom_infty;
VdW = VdW_d - VdW_infty + LJ_pair;
elec = elec_d + elec_d_outside - elec_infty + elec_pair;
total = geom + VdW + elec;
subplot(1,3,3); plot(4:1:24,total,'r-o');hold on;
xlim([4,24]); xlabel('d (\AA)','fontsize',20,'interpreter', 'latex');
ylim([-300,200]); ylabel('G_{tot}^{PMF}( k_BT )','fontsize',20,'rot',90);
h_legend = legend('- 0.2e/0.2e loose','- 0.2e/0.2e tight');
set(h_legend,'Location','SouthEast','FontSize',18);
set(gca,'xtick',[4,6,8,10,12,14,16,18,20,22,24],'fontsize',15);