function []=figureS14a() close all figure() hold on %% Data parameters conversion_factor300=1/1.00703256; % unit length definition for N=300 conversion_factor1000=1/1.00663684; % unit length definition for N=1000 conversion_factor8000=1/1.00648915; % unit length definition for N=8000 sample_indices=[3,5]; Ts=[0.101,0.050]; % temperatures %% Fit lines for sample_index_dummy=1:2 % color sample_index=sample_indices(sample_index_dummy); if sample_index==3 colorB=(4-1)/(11-1); end if sample_index==5 colorB=(6-1)/(11-1); end colorR=1-colorB; color=[colorR,0,colorB]; %% Loading and plotting data for N=1000 % Structural PTS data_name=sprintf('./PTS%d.dat',sample_index); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all'); NCP=round(size(numbers,1)/3); cavity_sizes=numbers(3*(1:NCP)-2); PTSs=numbers(3*(1:NCP)-1); PTS_errors=numbers(3*(1:NCP)); bx=conversion_factor1000*cavity_sizes; by=PTSs; be=PTS_errors; % Plot errorbar(bx,by,be,'Color',color,'Marker','.','MarkerSize',20,'LineStyle','-') %% Loading and plotting data for N=300 % Structural PTS data_name=sprintf('./PTS_SMALL%d.dat',sample_index); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all'); NCP=round(size(numbers,1)/3); cavity_sizes=numbers(3*(1:NCP)-2); PTSs=numbers(3*(1:NCP)-1); PTS_errors=numbers(3*(1:NCP)); bx=conversion_factor300*cavity_sizes; by=PTSs; be=PTS_errors; % Plot errorbar(bx,by,be,'Color',color,'Marker','.','MarkerSize',20,'LineStyle','--') %% Loading and plotting data for N=8000 % Structural PTS data_name=sprintf('./PTS_LARGE%d.dat',sample_index); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all'); NCP=round(size(numbers,1)/3); cavity_sizes=numbers(3*(1:NCP)-2); PTSs=numbers(3*(1:NCP)-1); PTS_errors=numbers(3*(1:NCP)); bx=conversion_factor8000*cavity_sizes; by=PTSs; be=PTS_errors; % Plot errorbar(bx,by,be,'Color',color,'Marker','.','MarkerSize',20,'LineStyle',':') end xlim([0,8]) ylim([0,1]) xlabel('$R$','Interpreter','latex'); ylabel('$Q$','Interpreter','latex') print(gcf, '-dpdf', 'figureS14a.pdf'); end