\documentclass[twocolumn]{article}
\setlength{\textwidth}{160mm}
\setlength{\oddsidemargin}{-0mm}
\setlength{\textheight}{220mm}
\setlength{\topmargin}{-8mm}
\usepackage{graphicx,color,hyperref}
\graphicspath{{./fig/}{./png/}}
\thispagestyle{empty}
\begin{document}

\begin{figure}[h!]\begin{center}
    \includegraphics[width=\columnwidth]{fig/EEKM-EEGW_CME+QCD.eps}
\end{center}\caption[]{EEGW vs EEKM with all CME and QCD runs.}\label{fig:EEGW-EEKM}
\end{figure}

In \autoref{fig:EEGW-EEKM}, Grey and red: \cite{Brandenburg:2021aln}; Blue: QCD runs \cite{Bran+Clar+He+Kahn21}; Black: new runs.

\begin{figure}[h!]\begin{center}
    \includegraphics[width=\columnwidth]{fig/RunA.eps}
\end{center}\caption[]{Magnetic and GW energy vs. viscosity.}\label{fig:energy-visc}
\end{figure}

\autoref{fig:energy-visc}: Run A from CME paper. Red is magnetic energy and blue is GW.

\begin{figure}[h!]\begin{center}
    \includegraphics[width=\columnwidth]{fig/EEKM-EEGW_CME+QCD_r2.eps}
\end{center}\caption[]{EEGW vs EEKM in regime II.}\label{fig:regimeII}
\end{figure}

In \autoref{fig:regimeII} I plotted the runs from Roper Pol et al. (2020) in grey, CME runs in regime II in red, and new runs in black with the extrapolated run in orange, which are all in regime II. The grey line is the quadratic relation.

\begin{figure}[h!]\begin{center}
    \includegraphics[width=\columnwidth]{fig/EEKM-EEGW_CME+extrap.eps}
\end{center}\caption[]{}\label{fig:no_qcd}
\end{figure}

I removed the QCD runs in \autoref{fig:no_qcd} (\autoref{fig:EEGW-EEKM} blue points are wrong; need to normalize the same way as EWPT runs).

\bibliographystyle{h-physrev5}
\bibliography{Chirality-BBN}
\end{document}