\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=plot1,
    width=0.65\textwidth,
    height=0.6\textwidth,
    xmin=0,
    xmax=3,
    ymin=0,
    ylabel=$f_\mathcal{Q}$,
    ylabel style = {yshift=-0.2cm},
    xlabel=$B$ (T),
    yticklabel style = {/pgf/number format/fixed, /pgf/number format/zerofill, /pgf/number format/precision=1},
    ytick pos=left,
    xtick pos=both,
    ytick align=outside,
    ymax=5,
    ytick = {0, 1, 2, 3, 4, 5},
    xticklabel pos=lower,
    xtick = {0.0,0.5, 1.0, ..., 2.5},
    xticklabel style = {/pgf/number format/fixed, /pgf/number format/zerofill, /pgf/number format/precision=1},
    legend style = {anchor = north east,at={(axis description cs: 0.985, 1)}},
    legend cell align={left},
    title = {(a)}, title style = {at = {(axis description cs:0.05,0.87)}},
]
\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[black, 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};

%Components DMRG QFI
\addplot[viryellowdark, mark=pentagon*, solid, mark options={solid}] table [x=<Field(T)>, y=<QFI>, col sep=space, y expr = \thisrowno{2}] {EntanglementData/DMRG_QFI_xx_vs_T.dat};
\addplot[virmiddle, mark=triangle*, solid, mark options={solid}] table [x=<Field(T)>, y=<QFI>, col sep=space, y expr = \thisrowno{2}] {EntanglementData/DMRG_QFI_yy_vs_T.dat};
\addplot[virdark, mark=square, solid, mark options={solid}] table [x=<Field(T)>, y=<QFI>, col sep=space, y expr = \thisrowno{2}] {EntanglementData/DMRG_QFI_zz_vs_T.dat};

\addplot [domain=0.0:3.5, 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{viryellowdark}{DMRG $S^{xx}$ w/o PF},\textcolor{virmiddle}{DMRG $S^{yy}$ w/o PF},\textcolor{virdark}{DMRG $S^{zz}$ w/o PF}};
\end{axis}
%
%
\begin{axis}[name=plot2,at={(plot1.north east)},anchor=north west,
width=0.65\textwidth,
height=0.6\textwidth,
xmin=0,
xmax=3.0,
ymin=0,
ylabel=,
yticklabel style = {/pgf/number format/fixed, /pgf/number format/zerofill, /pgf/number format/precision=1, },%xshift=0.8cm},
tickpos=left,
ytick align=inside,
ymin=0.0,
ymax=5.0,
ytick = {0.0, 1.0, 2.0, ..., 5.0},
xtick = {0.0,0.5, 1.0, ..., 2.5},
xticklabel style = {/pgf/number format/fixed, /pgf/number format/zerofill, /pgf/number format/precision=1, },%xshift=0.8cm},
xlabel=$B$ (T),
xtick pos=both,
yticklabels = {,,},
legend style = {anchor = north east,at={(axis description cs: 0.985, 1)}},
legend cell align={left},
title = {(b)}, title style = {at = {(axis description cs:0.05,0.87)}},
]
\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};

%Components experimental QFI
\addplot[viryellowdark, mark=pentagon*, solid, mark options={solid}] table [x=<Field(T)>, y=<QFI>, col sep=space, y expr = \thisrowno{1}] {EntanglementData/Experiment_Corrected_QFI_xx_vs_T.dat};
\addplot[virmiddle, mark=triangle*, solid, mark options={solid}] table [x=<Field(T)>, y=<QFI>, col sep=space, y expr = \thisrowno{1}] {EntanglementData/Experiment_Corrected_QFI_yy_vs_T.dat};
\addplot[virdark, mark=square, solid, mark options={solid}] table [x=<Field(T)>, y=<QFI>, col sep=space, y expr = \thisrowno{1}] {EntanglementData/Experiment_Corrected_QFI_zz_vs_T.dat};

\addplot [domain=0.0:3.5, black, dashed]{ 3 };
\addplot +[mark=none,solid,virdark] coordinates {(\xclasT, 0) (\xclasT, 6.05)};


\legend{\textcolor{virmiddle}{INS Cs$_2$CoCl$_4$},\textcolor{viryellowdark}{Corrected INS},
	\textcolor{viryellowdark}{Exp. $S^{xx}$ w/o PF},\textcolor{virmiddle}{Exp. $S^{yy}$ w/o PF},\textcolor{virdark}{Exp. $S^{zz}$ w/o PF}};
\end{axis}
%
\begin{axis}[name=plot3, at={(plot2.north east)},
anchor=north west,
width=0.65\textwidth,
height=0.2871\textwidth,
xmin=0,
xmax=3,
ymin=0,
ymax=2,
xticklabel style = {/pgf/number format/fixed, /pgf/number format/zerofill, /pgf/number format/precision=1},
xtick = {0.0, 0.5, 1.0, ..., 3},
xtick pos=both,
xticklabels={},
ylabel=,
ylabel style = {yshift=-0.2cm},
yticklabel style = {/pgf/number format/fixed, /pgf/number format/zerofill, /pgf/number format/precision=1},%, xshift=0.7cm},
ytick pos=right,
ytick = {0,1,2,3,4},
ytick align=outside,
y tick label style = {black},
legend style = {anchor = north east,at={(axis description cs: 0.985, 1)}},
legend cell align={left},
title = {(c) $S^{xx}$}, title style = {at = {(axis description cs:0.1,0.)}}, ]


\addplot[viryellowdark, mark=pentagon, solid, mark options={solid}] table [x=<Field(T)>, y=<QFI>, col sep=space, y expr = \thisrowno{2}] {EntanglementData/DMRG_QFI_xx_vs_T.dat};
\addplot[viryellowdark, mark=pentagon*, solid, mark options={solid}] table [x=<Field(T)>, y=<QFI>, col sep=space, y expr = \thisrowno{1}] {EntanglementData/Experiment_Corrected_QFI_xx_vs_T.dat};

\legend{\textcolor{viryellowdark}{DMRG},\textcolor{viryellowdark}{Exp.}};
\end{axis}


\begin{axis}[name=plot4, at={(plot3.south west)},
anchor=north west,
width=0.65\textwidth,
height=0.2871\textwidth,
xmin=0,
xmax=3,
ymin=0,
ymax=2.0,
xticklabel style = {/pgf/number format/fixed, /pgf/number format/zerofill, /pgf/number format/precision=1},
xtick = {0.0, 0.5, 1.0, ..., 3},
xtick pos=both,
xticklabels={},
ylabel=,
ylabel style = {yshift=-0.2cm},
yticklabel style = {/pgf/number format/fixed, /pgf/number format/zerofill, /pgf/number format/precision=1},
ytick pos=right,
ytick = {0,1},
ytick align=outside,
y tick label style = {black},
legend style = {anchor = north east,at={(axis description cs: 0.72, 1)}},
legend cell align={left},
title = {(d) $S^{yy}$}, title style = {at = {(axis description cs:0.1,0.)}}, ]

\addplot[virmiddle, mark=triangle, solid, mark options={solid}] table [x=<Field(T)>, y=<QFI>, col sep=space, y expr = \thisrowno{2}] {EntanglementData/DMRG_QFI_yy_vs_T.dat};
\addplot[virmiddle, mark=triangle*, solid, mark options={solid}] table [x=<Field(T)>, y=<QFI>, col sep=space, y expr = \thisrowno{1}] {EntanglementData/Experiment_Corrected_QFI_yy_vs_T.dat};



\legend{\textcolor{virmiddle}{DMRG},\textcolor{virmiddle}{Exp.}};
\end{axis}


\begin{axis}[name=plot5, at={(plot4.south west)},
anchor=north west,
width=0.65\textwidth,
height=0.2871\textwidth,
xlabel=$B$ (T),
xmin=0,
xmax=3,
ymin=0,
ymax=2.0,
xticklabel style = {/pgf/number format/fixed, /pgf/number format/zerofill, /pgf/number format/precision=1},
xtick = {0.0, 0.5, 1.0, ..., 3},
xtick pos=both,
ylabel=,
ylabel style = {yshift=-0.2cm},
yticklabel style = {/pgf/number format/fixed, /pgf/number format/zerofill, /pgf/number format/precision=1},
ytick pos=right,
ytick = {0,1},
ytick align=outside,
y tick label style = {black},
legend style = {anchor = north east,at={(axis description cs: 0.985, 1)}},
legend cell align={left},
title = {(e) $S^{zz}$}, title style = {at = {(axis description cs:0.1,0.0)}}, ]

\addplot[virdark, mark=square, solid, mark options={solid}] table [x=<Field(T)>, y=<QFI>, col sep=space, y expr = \thisrowno{2}] {EntanglementData/DMRG_QFI_zz_vs_T.dat};
\addplot[virdark, mark=square*, solid, mark options={solid}] table [x=<Field(T)>, y=<QFI>, col sep=space, y expr = \thisrowno{1}] {EntanglementData/Experiment_Corrected_QFI_zz_vs_T.dat};


\legend{\textcolor{virdark}{DMRG},\textcolor{virdark}{Exp.}};
\end{axis}
%
\end{tikzpicture}

\end{document}