\documentclass{standalone}
\input{palette}
\usepackage{pgfplots}
\pgfplotsset{compat=default}
% With this compatability setting, pgfplots runs in backwards compatability mode using pgfplots version pre-1.3. It's safe to use a more current version, but axis and tick labels may move and need to be repositioned.

%\usepackage{txfonts}
%\usepackage[mathscr]{euscript}
% The above font packages were used to generate figures for the article, but are here commented out to ensure better compability.

\usetikzlibrary{pgfplots.groupplots}

\newcommand{\mytick}[2]{\draw[thick,#2] (#1, 0.01) -- (#1, -0.15);}

\begin{document}
\begin{tikzpicture}
\def\xclas{1.58}
\def\xclasT{1.58*2.1/1.604}
\begin{axis}[name=inset1,
	footnotesize,
	xshift=.17\textwidth,
	yshift=0.028\textwidth,
	width=0.3\textwidth,
	height=0.225\textwidth,
	xmin=1.58,
	xmax=1.61,
	ymin=0.4,
	ymax=1.2,
    ytick = {0.3, 0.5, ..., 0.9, 1.1},
    yticklabel style={xshift=0.05cm},
    xticklabel style={yshift=0.05cm},
	]
\addplot [draw = black, color = black, mark=diamond*, densely dashed, mark options={solid}] table [x=<h>, y=<SvN>, y expr=\thisrowno{2}/ln(2), col sep=space] {EntanglementData/EntanglementEntropy_Delta025_L050_h.dat};
\addplot [draw = virmiddle, color = virmiddle, mark=*, solid, mark options={solid}] table [x=<h>, y=<SvN>, y expr=\thisrowno{2}/ln(2), col sep=space] {EntanglementData/EntanglementEntropy_Delta025_L100_h.dat};
\addplot [draw = viryellowdark, color = viryellowdark, mark=x, dashed, mark options={solid}] table [x=<h>, y=<SvN>, y expr=\thisrowno{2}/ln(2), col sep=space] {EntanglementData/EntanglementEntropy_Delta025_L150_h.dat};
\addplot [draw = virdark, color = virdark, mark=o, densely dotted, mark options={solid}] table [x=<h>, y=<SvN>, y expr=\thisrowno{2}/ln(2), col sep=space] {EntanglementData/EntanglementEntropy_Delta025_L200_h.dat};
\end{axis}

\begin{axis}[name=plot1,
    width=0.9\textwidth,
    height=0.33 \textwidth,
    xmin=0,
    xmax=2.6,
    ymin=0,
    ylabel=$S_\mathrm{vN}/\ln 2$,
    ylabel style = {yshift=-0.2cm},
    xlabel=$h_x/J$,
    xlabel style = {at = {(axis description cs: 0.95,1.4025)}},
    yticklabel style = {/pgf/number format/fixed, /pgf/number format/zerofill, /pgf/number format/precision=1},
    ytick pos=left,
    xtick pos=upper,
    ytick align=outside,
    ymax=2.0,
    ytick = {0.0, 0.5, 1.0, ..., 2.0},
    xticklabel pos=upper,
    xtick = {0.0,0.2, 0.4, ..., 2.3},
    xticklabel style = {/pgf/number format/fixed, /pgf/number format/zerofill, /pgf/number format/precision=1,
    yshift=-.05cm},
    legend style = {anchor = south east,at={(axis description cs: 1, 0.09)},
    	inner sep=0.5pt},
    legend cell align={left},
    title = {(a)}, title style = {at = {(axis description cs:0.97,0.715)}},
]
\addplot [draw = black, color = black, mark=diamond*, densely dashed, mark options={solid}] table [x=<h>, y=<SvN>, y expr=\thisrowno{2}/ln(2), col sep=space] {EntanglementData/EntanglementEntropy_Delta025_L050_h.dat};
\addplot [draw = virmiddle, color = virmiddle, mark=*, solid, mark options={solid}] table [x=<h>, y=<SvN>, y expr=\thisrowno{2}/ln(2), col sep=space] {EntanglementData/EntanglementEntropy_Delta025_L100_h.dat};
\addplot [draw = viryellowdark, color = viryellowdark, mark=x, dashed, mark options={solid}] table [x=<h>, y=<SvN>, y expr=\thisrowno{2}/ln(2), col sep=space] {EntanglementData/EntanglementEntropy_Delta025_L150_h.dat};
\addplot [draw = virdark, color = virdark, mark=o, densely dotted, mark options={solid}] table [x=<h>, y=<SvN>, y expr=\thisrowno{2}/ln(2), col sep=space] {EntanglementData/EntanglementEntropy_Delta025_L200_h.dat};
\addplot [domain=\xclas:3.5, virdark]{ ln(2)/ln(2) };
\addplot +[mark=none,solid,virdark] coordinates {(\xclas, 0) (\xclas, 2.5)};
\legend{\textcolor{black}{DMRG $L=50$}, \textcolor{virmiddle}{DMRG $L=100$}, \textcolor{viryellowdark}{DMRG $L=150$}, \textcolor{virdark}{DMRG $L=200$}};
\end{axis}
%
%
\begin{axis}[name=plot2,at={(plot1.south west)},anchor=north west,
width=0.9\textwidth,
height=0.31\textwidth,
xmin=0,
xmax=2.6,
ymin=0,
xticklabels = {,,},
ylabel=$\tau_n$,
ylabel style = {yshift=-0.2cm},
yticklabel style = {/pgf/number format/fixed, /pgf/number format/zerofill, /pgf/number format/precision=1, },
tickpos=left,
ytick align=outside,
ymin=0.0,
ymax=1.0,
ytick = {0.0, 0.2, 0.4, ..., 0.8},
xtick = {0.0,0.2, 0.4, ..., 2.3},
xtick pos=upper,
legend style = {anchor = east,at={(axis description cs: 1, 0.49)}, 
	inner sep=0.5pt},
legend cell align={left},
title = {(b)}, title style = {at = {(axis description cs:0.97,0.67)}},
]
% DMRG tau1
\addplot[restrict x to domain=0:2.2][fill= white, draw = virdark, mark=*, dashed, mark options={solid}] table [x=<h>, y=<MySt>, col sep=space, y expr = 1.0-4.0*((\thisrowno{2}*\thisrowno{2})+(\thisrowno{6}*\thisrowno{6}))] {DMRG_Data/StaggeredMagnetization_h_mixedm.dat};
% Experimental tau1
\addplot[restrict x to domain=0:3.5][draw = virdark, mark=diamond*, dashdotted, mark options={solid}] table [x=<B>, y=<tau1>, col sep=space, x expr=1.604/2.1*\thisrowno{0}] {EntanglementData/Experiment_OneTangle.dat};
% DMRG tau2
\addplot[restrict x to domain=0:2.2][draw = virmiddle, color = virmiddle, mark=*, solid, mark options={solid}] table [x=<h>, y=<tau2>, col sep=space, y expr = \thisrowno{2}]{EntanglementData/TwoTangle_h.dat};

\legend{\textcolor{virdark}{$\tau_1$ DMRG}, \textcolor{virdark}{$\tau_1$ Cs$_2$CoCl$_4$}, \textcolor{virmiddle}{$\tau_2$ DMRG}};
\addplot +[mark=none,solid,virdark] coordinates {(\xclas, 0) (\xclas, 1.5)};
\end{axis}
%
\begin{axis}[at={(plot2.south west)},
anchor=north west,
width=0.9\textwidth,
height=0.37\textwidth,
xmin=0,
xmax=2.6*2.1/1.604,
ymin=0,
ymax=5,
xtick = {0.0, 0.5, 1.0, ..., 2},
xtick pos=top,
xticklabels={},
xlabel=,
ylabel=,
yticklabels={},
ytick style={draw=none}]
\end{axis}
%
\begin{axis}[at={(plot2.south west)},
anchor=north west,
width=0.9\textwidth,
height=0.37\textwidth,
xlabel=$B (T)$,
xlabel style={yshift=.1cm,xshift=.1cm},
xmin=0,
xmax=2.6*2.1/1.604,
ymin=0,
ymax=5,
xticklabel style = {/pgf/number format/fixed, /pgf/number format/zerofill, /pgf/number format/precision=1,
yshift=.05cm},
xtick = {0.0, 0.5, 1.0, ..., 3},
xtick pos=bottom,
ylabel=$f_\mathcal{Q}$,
ylabel style = {yshift=-0.2cm},
yticklabel style = {/pgf/number format/fixed, /pgf/number format/zerofill, /pgf/number format/precision=0},
ytick pos=left,
ytick = {0,1,2,3,4},
ytick align=outside,
y tick label style = {black},
legend style = {anchor = north east,at={(axis description cs: .923, 1)}, fill=none,
inner sep=0.5pt
},
legend cell align={left},
title = {(c)}, title style = {at = {(axis description cs:0.97,0.75)}}, xlabel style = {at = {(axis description cs: 0.95,0.16)}}]

\addplot[black, mark=diamond, dotted, mark options={solid}] table [x=<Field(T)>, y=<QFI>, col sep=space, y expr = \thisrowno{2}] {EntanglementData/DMRG_QFI_vs_T_NoPolarizationFactor.dat};
\addplot[virdark, mark=diamond*, solid, mark options={solid}] table [x=<Field(T)>, y=<QFI>, col sep=space, y expr = \thisrowno{2}] {EntanglementData/DMRG_QFI_vs_T.dat};
\addplot[virmiddle, mark=*, dashed, mark options={solid}] table [x=<Field(T)>, y=<QFI>, col sep=space, x expr = \thisrowno{0}, y expr = \thisrowno{1}] {EntanglementData/Experiment_QFI_vs_T.dat};
\addplot[viryellowdark, mark=square*, densely dashed, mark options={solid}] table [x=<Field(T)>, y=<QFI>, col sep=space, x expr = \thisrowno{0}, y expr = \thisrowno{1}] {EntanglementData/Experiment_Corrected_QFI_vs_T.dat};


\addplot [domain=0.0:\xclasT, black, dashed]{ 3 };
\addplot +[mark=none,solid,virdark] coordinates {(\xclasT, 0) (\xclasT, 6.05)};
\legend{\textcolor{black}{DMRG w/o PF},\textcolor{virdark}{DMRG w. PF},\textcolor{virmiddle}{INS Cs$_2$CoCl$_4$},\textcolor{viryellowdark}{Corrected INS}};
\end{axis}
%
\end{tikzpicture}

\end{document}