diff --git a/Carlen_Holzmann_Jauslin_Lieb_2020.tex b/Carlen_Holzmann_Jauslin_Lieb_2020.tex
index 3e2a478..67398df 100644
--- a/Carlen_Holzmann_Jauslin_Lieb_2020.tex
+++ b/Carlen_Holzmann_Jauslin_Lieb_2020.tex
@@ -35,13 +35,13 @@
 \email{lieb@princeton.edu}
 
 \begin{abstract}
-In 1963, a {\it Simple Approach} was developed to study the ground state energy of an interacting Bose gas.
+In 1963, a {\it Simplified Approach} was developed to study the ground state energy of an interacting Bose gas with a purely repulsive potential.
 It consists in the derivation of an Equation, which is not based on perturbation theory, and which gives the exact expansion of the energy at low densities.
 This Equation is expressed directly in the thermodynamic limit, and only involves functions of $3$ variables, rather than $3N$.
-Here, we revisit this approach, and show that the Equation yields accurate predictions for various observables for {\it all} densities.
-Specifically, in addition to the ground state energy, we have shown that the Simple Approach gives predictions for the condensate fraction, two-point correlation function, and momentum distribution.
-We have carried out a variety of tests by comparing the predictions of the Equation with Quantum Monte Carlo calculations, and have found remarkable agreement.
-We thus show that the Simple Approach provides a new theoretical tool to understand the behavior of the many-body Bose gas, not only in the small and large density ranges, which have been studied before, but also in the range of intermediate density, for which little is known.
+Here, we revisit this approach, introduce two more equations and show that these yields accurate predictions for various observables for {\it all} densities for repulsive potentials with positive Fourier transform.
+Specifically, in addition to the ground state energy, we have shown that the Simplified Approach gives predictions for the condensate fraction, two-point correlation function, and momentum distribution.
+We have carried out a variety of tests by comparing the predictions of the Equations with Quantum Monte Carlo calculations for exponential interaction potentials as well as a different, finite range potential of positive type, and have found remarkable agreement.
+We thus show that the Simplified Approach provides a new theoretical tool to understand the behavior of the many-body Bose gas, not only in the small and large density ranges, which have been studied before, but also in the range of intermediate density, for which much less is known.
 \end{abstract}
 
 \maketitle
@@ -51,7 +51,7 @@ We thus show that the Simple Approach provides a new theoretical tool to underst
 Bose gases are one of the foundational objects in the statistical mechanics of quantum systems, and have been the focus of much scrutiny, dating back to the early days of quantum mechanics\-~\cite{Le29}.
 Nevertheless, there are still several important problems to be solved, in the case of interacting Bose gases, in which the correlations between particles make the analysis very difficult.
 In this case, observables may be computed by either performing numerical computations using finite-size approximations and extrapolations, or by devising effective theories which capture some of the correlations between particles, while remaining integrable.
-In this paper, we present an effective theory which goes back to 1963\-~\cite{Li63}, and which we have found gives astonishingly accurate predictions in the thermodynamic limit at {\it all} densities that have been verified numerically by Quantum Monte Carlo (QMC) computations.
+In this paper, we present an effective theory which goes back to 1963\-~\cite{Li63}, and which we have found gives accurate predictions in the thermodynamic limit at {\it all} densities that have been verified numerically by Quantum Monte Carlo (QMC) computations.
 This remarkable agreement leads us to suggest that this may be a new way of understanding and analyzing the quantum many-body problem.
 \bigskip
 
@@ -97,9 +97,9 @@ Note that, whereas Hartree theory is accurate at asymptotically large densities,
 
 \indent
 Therefore, the Bose gas is described by Bogolubov theory at low density, and Hartree theory or the MSA at high density.
-In this paper, we will discuss another effective theory for the ground state of the repulsive Bose gas with a positive type potential, which is highly accurate at all densities.
+In this paper, we will discuss another effective theory for the ground state of the repulsive Bose gas with a positive type potential, which is highly accurate at all densities, which is {\it exact} at low and high densities, and highly accurate at all intermediate densities.
 In other words, it is a physically descriptive interpolation between Bogolubov and Hartree theory.
-To justify our claim that it is in good {\it quantitative} agreement  the physics all all densities, we rely on with QMC simulations of the Bose gas for intermediate densities.
+To justify our claim that it is in good {\it quantitative} agreement with the physics at all densities, we rely on with QMC simulations of the Bose gas for intermediate densities.
 This equation was originally introduced in 1963\-~\cite{Li63}, and studied for the high density Jellium\-~\cite{LS64}, and in one dimension\-~\cite{LL64}.
 There has been no research progress since then.
 The merit of this equation is twofold.
@@ -209,7 +209,7 @@ In the present paper, we discuss some more quantitative results, with more of a
 We will consider potentials that are of positive type, with a special focus on exponential potentials of the form $\alpha e^{-|\mathbf x|}$.
 We have found that the prediction for the energy is very accurate for {\it all} densities, see Figure\-~\ref{fig:energy}.
 In the case $\alpha=1$, the relative error compared to the QMC simulation is as small as $0.1\%$, and is comparable to the error made by a Bijl-Dingle-Jastrow function Ansatz \cite{Bi40,Di49,Ja55}, see Figure\-~\ref{fig:cmp_energy}, even though the solution of the Big Equation is much easier to compute numerically than the Bijl-Dingle-Jastrow optimizer.
-The prediction for the condensate fraction is less accurate, though still remarkably good for small values of $\alpha$, see Figure\-~\ref{fig:condensate0.5}.
+The prediction for the condensate fraction is less accurate in the intermediate density regime, though still remarkably good for small values of $\alpha$, see Figure\-~\ref{fig:condensate0.5}.
 For larger $\alpha$, the Big Equation is off the mark, see Figure\-~\ref{fig:condensate16}, although the qualitative features of the condensate fraction are still well reproduced.
 We have also carried out similar computations for the hard core potential, for which we also find good agreement, see Figure\-~\ref{fig:hardcore}.
 
@@ -238,10 +238,16 @@ For the Big and Simple Equations discussed in this paper, we have found that thi
   \sqrt{\rho a_0}\ll|\mathbf k|\ll1
 \end{equation}
 which is another confirmation of the accuracy of the effective equation at small densities.
-However, if $\sqrt\rho\gtrsim1$, then the Tan regime does not exist, and the picture in terms of strongly coupled few-particle configurations inherent to the analysis of unitary Bose gases\-~\cite{CW11,SBe14} breaks down, as the Bose gas transitions to a strongly correlated liquid.
+However, if $\sqrt\rho\gtrsim1$, then the universal Tan regime does not exist, and the picture in terms of strongly coupled few-particle configurations inherent to the analysis of unitary Bose gases\-~\cite{CW11,SBe14} breaks down, as the Bose gas transitions to a strongly correlated liquid.
 This is confirmed for the prediction of the Big Equation, see Figure\-~\ref{fig:tan}.
 \bigskip
 
+\indent
+As further evidence of the breakdown of universality in the intermediate density regime, we have also compared the ground state energy for two very different potentials, which have the same scattering length and the same integral.
+We have found that the energy for these two potentials is significantly different in the intermediate density regime, see Figure\-~\ref{fig:compare_pots}.
+For these two potentials, we have also found that the Quantum Monte Carlo data fits very well with the prediction of the Big Equation.
+\bigskip
+
 \indent
 The rest of the paper is structured as follows.
 In section\-~\ref{sec:approx}, we detail the approximation needed to get from the many-body Bose gas to the Full Equation, and then discuss the approximations leading to the Big, Medium and Simple Equations.
@@ -261,6 +267,9 @@ We start from the many-body Hamiltonian: denoting the number of particles by $N$
 We confine the $N$ particles in a cubic box $\Lambda$  of volume $V$, and impose periodic boundary conditions.
 Later on, we will take the thermodynamic limit $N,V\to\infty$, $\frac NV=\rho$ fixed.
 
+\indent
+In the derivation presented here, we will rely on the translation invariance of the Hamiltonian, which does not allow us to study a system with a trapping potential at this time.
+
 \indent Let  $E_N$ denote the ground state energy and let $\psi_N(\mathbf x_1,\cdots,\mathbf x_N)$  denote the ground state wave function so that
 \begin{equation}
   H\psi_0(\mathbf x_1,\cdots,\mathbf x_N)=E_N\psi_N(\mathbf x_1,\cdots,\mathbf x_N)
@@ -366,7 +375,7 @@ At high densities, since the system approaches a mean-field regime, one might al
 The Full Equation we have derived is quite difficult to study, even numerically.
 As was discussed in Section\-~\ref{sec:intro}, we will introduce further approximations to simplify the equation.
 The first approximation is to neglect the $\frac12u(\mathbf z)u(\mathbf y-\mathbf x)$ term in\-~(\ref{L}), which is the most difficult term, from a computational point of view. 
-We expect that, at low densities, this term is expected to be of order $\rho^{3/2}$ uniformly $\mathbf x$, whereas the leading order term in $L$ should be of order $\rho$.
+We expect that, at low densities, this term is expected to be of order $\rho^{3/2}$ uniformly in $\mathbf x$, whereas the leading order term in $L$ should be of order $\rho$.
 This leads us to the Big Equation defined in\-~(\ref{bigeq}).
 This equation is easier to solve numerically than the Full Equation, because in Fourier space, it involves only two convolution operators, whereas the Full Equation contains three, which makes it computationally heavier.
 Nevertheless, this equation is still difficult to study analytically, so we make further approximations
@@ -385,7 +394,7 @@ On account of (\ref{EN}), the function $S(\mathbf x)$ defined in (\ref{K}) satis
   \int d\mathbf x\ S(\mathbf x) = \frac{2\tilde e}{\rho}
 \end{equation}
 which is just another way of stating (\ref{erel}).
-There are two different length scales in the problem: the first is the scattering length of the potential $a_0$ and the interparticle distance $\rho^{-1/3}$.
+There are two different length scales in the problem: the first is the scattering length of the potential $a_0$ and the second is the interparticle distance $\rho^{-1/3}$.
 At sufficiently low densities we will have 
 \begin{equation}
   a_0 \ll  \rho^{-1/3}
@@ -429,7 +438,7 @@ mixed-estimator bias occurring for observables different from the ground state e
 In principle, the mixed-estimator bias can be controlled either by systematic improvement of
 the trial wave function \cite{RMH18} or by different projection Monte Carlo methods, e.g. Reptation 
 Monte Carlo \cite{BM99}. For the system under consideration, the mixed estimator bias
-of the pair-product wave function was found to be sufficiently accurate, 
+of the pair-product wave function was found to be sufficiently small, 
 the overall precision being limited rather by
 the finite system size of the QMC calculations.
 \bigskip
@@ -473,6 +482,7 @@ For $\alpha=16$ this is even clearer, and one sees that the Medium Equation is m
   \caption{
     The energy as a function of density for the potential $e^{-|\mathbf x|}$ (top) and $16e^{-|\mathbf x|}$ (bottom).
     We compare the predictions of the Big, Medium and Simple Equations to a QMC simulation.
+    For comparison, we also plot the Lee-Huang-Yang (LHY) energy\-~(\ref{lhy}).
   }
   \label{fig:energy}
 \end{figure}
@@ -506,10 +516,11 @@ The approximations leading to the Big, Simple and Medium Equations reduce the nu
 In doing so, we lose some information, and, in particular, we do not obtain a prediction for the many-body wavefunction $\psi_0$.
 Therefore, computing observables other than the ground state energy is not entirely straightforward.
 To compute the condensate fraction, we first express it in terms of the energy of an auxiliary system, from which we derive an approximation following the prescriptions in section\-~\ref{sec:approx}.
-Specifically, the condensate fraction of the many-body ground state $\psi_0$ is in terms of the projector $P_i\psi_0:=\int\frac{d\mathbf x_i}V\psi_0$ onto the condensate wavefunction (which is the constant function):
+Specifically, the {\it non}-condensed fraction of the many-body ground state $\psi_0$
 \begin{equation}
   \eta_0:=1-\frac1N\sum_{i=1}^N\left<\psi_0\right|P_i\left|\psi_0\right>
 \end{equation}
+is expressed in terms of the projector $P_i\psi_0:=\int\frac{d\mathbf x_i}V\psi_0$ onto the condensate wavefunction (which is the constant function):
 which we re-express in terms of the modified Hamiltonian
 \begin{equation}
   H_\mu=-\frac12\sum_{i=1}^N\Delta_i+\sum_{1\leqslant i<j\leqslant N}v(\mathbf x_i-\mathbf x_j)-\mu\frac1N\sum_{i=1}^NP_i
@@ -571,6 +582,7 @@ The location of the maximum of the non-condensed fraction (or the minimum of the
   \caption{
     The non-condensed fraction as a function of the density for the potential $\frac12e^{-|\mathbf x|}$.
     We compare the predictions of the Big, Medium and Simple Equations to a QMC simulation.
+    The prediction of Bogolubov theory\-~(\ref{eta0}) is also plotted for comparison (Bog).
   }
   \label{fig:condensate0.5}
 \end{figure}
@@ -590,12 +602,12 @@ We first note that this can be rewritten in a way that makes the translation inv
 \end{equation}
 which we can rewrite as a functional derivative of the ground state energy per-particle $e_0$:
 \begin{equation}
-  C_2(\mathbf x)=2\rho\frac{\delta e_0}{\delta v(\mathbf x)}
+  C_2(\mathbf x)=2\rho^2\frac{\delta e_0}{\delta v(\mathbf x)}
   .
 \end{equation}
 The prediction $\tilde C_2$ of the Big and Medium Equations for the two-point correlation function are therefore defined by differentiating $\tilde e$ in\-~(\ref{erel}) with respect to $v$:
 \begin{equation}
-  \tilde C_2(\mathbf x):=2\rho\frac{\delta\tilde e}{\delta v(\mathbf x)}
+  \tilde C_2(\mathbf x):=2\rho^2\frac{\delta\tilde e}{\delta v(\mathbf x)}
   .
   \label{C2}
 \end{equation}
@@ -632,6 +644,7 @@ The prediction of the Big Equation remains quite accurate, when compared to the
 The presence of this local maximum in $g_2$ shows that, in the probability distribution $\psi_0$, there is a larger probability of finding pairs of particles that are separated by a certain fixed distance.
 This indicates the appearance of a new physical length scale at intermediate densities, and indicates that the system exhibits a non-trivial physical behavior in this regime.
 Note that this behavior was observed for the stronger potential $16e^{-|\mathbf x|}$, but seems to be absent for $e^{-|\mathbf x|}$.
+Note, also, that, as will be discussed next, this maximum is also present in the two-point correlation $C_2$, and is, therefore, the manifestation of a physical phenomenon.
 \bigskip
 
 \begin{figure}
@@ -650,6 +663,7 @@ At low densities, the prediction of the Big Equation agrees rather well with the
 The Simple and Medium Equations are not as accurate.
 At larger densities, the Simple and Medium Equations are quite far from the QMC computation, and the Big Equation is not as accurate as in the case of $\tilde g_2$, but it does reproduce some of the qualitative behavior of the QMC computation.
 In particular, there is a local maximum in the two-point correlation function, which occurs at a length scale that is close to that observed for $\tilde g_2$.
+This suggests the emergence of a non-trivial phase, which resembles a liquid.
 At small $\mathbf x$, $\tilde C_2$ is negative, which is clearly not physical, and those values should be discarded.
 
 \begin{figure}
@@ -739,6 +753,49 @@ This is particularly true in the region in which $|\mathbf k|^{-4}$ transitions
   \label{fig:tan}
 \end{figure}
 
+
+\subsection{Non-universal behavior at intermediate densities}
+\indent
+The low density asymptotics of the energy, given by the Lee-Huang-Yang formula\-~(\ref{lhy}), only depend on the potential through the scattering length.
+At high density\-~(\ref{ehigh}), they only depend on the potential through $\int d\mathbf x\ v(\mathbf x)$.
+In this sense, the low and high density behavior of the Bose gas is {\it universal}.
+In this section, we show that, at intermediate densities, the energy does not only depend on the scattering length and the integral of the potential, thus suggesting that the behavior of the Bose gas at intermediate densities is {\it not universal}.
+
+\indent
+To that end, we have compared the predictions of the Big Equation for the energy for two potentials that have the same scattering length, and the same integral.
+The first potential, $v_{32}^{(0)}$, is defined in the next section, see\-~(\ref{vn}), and the second is an exponential potential
+\begin{equation}
+  \Phi_{\alpha,\beta}(\mathbf x):=\alpha e^{-\beta|\mathbf x|}
+\end{equation}
+where $\alpha$ and $\beta$ are chosen in such a way that the scattering length and integral of $\Phi$ are equal to those of $v_{32}^{(0)}$.
+The scattering length of $v_{32}^{(0)}$ was computed numerically and found to be $\approx0.5878$, and its integral is $\frac{64\pi^2}9$.
+The scattering length of $\Phi_{\alpha,\beta}$ is
+\begin{equation}
+  \frac1\beta\left(\log\frac\alpha{\beta^2}+2\gamma+2\frac{K_0(2\sqrt{\frac\alpha{\beta^2}})}{I_0(2\sqrt{\frac\alpha{\beta^2}})}\right)
+\end{equation}
+where $\gamma$ is the Euler constant and $K_0$ and $I_0$ are modified Bessel functions.
+The integral of $\Phi_{\alpha,\beta}$ is $\frac{8\pi\alpha}{\beta^3}$.
+We thus find that, in order to make the scattering length and integral of $v_{32}^{(0)}$ and $\Phi_{\alpha,\beta}$ coincide, we must choose $\alpha\approx907.2$ and $\beta\approx6.874$.
+
+\indent
+The prediction of the energy for these two potentials is plotted in Figure\-~\ref{fig:compare_pots}.
+We find that, as expected, the energies coincide at low and high density, but they differ significantly in the intermediate density regime.
+We have confirmed this fact by QMC computations, and found good agreement of the QMC data with our prediction for both potentials.
+
+\begin{figure}
+  \hfil\includegraphics[width=8cm]{compare_pots.pdf}
+  \caption{
+    The prediction of the energy by the Big Equation for the potentials $v_{32}^{(0)}$ and $\Phi(\mathbf x)\equiv\alpha e^{-\beta|\mathbf x|}$ with $\alpha\approx 907.2$ and 
+$\beta \approx 6.873$.
+    The potentials are chosen to have the same scattering length, $a_0\approx 0.5878$, as well
+as the same value for their integrals, so they coincide at low and at high densities.
+    They differ signinficantly at intermediate densities.
+    We compare each curve to a few QMC points, which fit well.
+    We also plot the Lee-Huang-Yang (LHY) energy\-~(\ref{lhy}).
+  }
+  \label{fig:compare_pots}
+\end{figure}
+
 \section{Hard-core potential}\label{sec:hardcore}
 
 \indent
@@ -747,18 +804,21 @@ We have investigated two directions to get around this restriction.
 
 \indent
 The first, and most straightforward, is to consider the hard-core potential as a limit of soft core potentials.
+Obviously, this approach will not be accurate at densities approaching close-packing, but as we will see, is rather accurate at smaller densities.
 As was mentioned in section\-~\ref{sec:intro}, it is preferable to only use potentials of positive type (that is, non-negative potentials with a non-negative Fourier transform).
 With this in mind, we consider the sequence of potentials
-\begin{equation}
-  v^{(0)}_n(|\mathbf x|):=\alpha_n\mathds 1_{|\mathbf x|<\frac12}\ast\mathds 1_{|\mathbf x|<\frac12}
-\end{equation}
-that is
 \begin{equation}
   v^{(0)}_n(|\mathbf x|)
-  =\mathds 1_{|\mathbf x|<1}
+  :=\Theta(1-|\mathbf x|)
   \alpha_n\frac{2\pi}{3}(|\mathbf x|-1)^2(|\mathbf x|+2)
 \end{equation}
-where $\mathds 1_{|\mathbf x|<\frac12}$ is the indicator function that $|\mathbf x|<\frac12$, and $\alpha_n\to\infty$ ($v^{(0)}_n$ has positive type because it is the convolution of a function with itself).
+where $\Theta(x)$ is the Heaviside function, which is equal to $1$ for $x>0$ and $0$ otherwise, and $\alpha_n\to\infty$.
+This potential can also be written as
+\begin{equation}
+  v^{(0)}_n(|\mathbf x|)=
+  \alpha_n\int d\mathbf y\ \Theta({\textstyle\frac12-|\mathbf y|})\Theta({\textstyle\frac12-|\mathbf x-\mathbf y|})
+\end{equation}
+which shows that it is of positive type because it is the convolution of the function $\Theta(\frac12-|\mathbf x|)$ with itself.
 In addition, we fix the scattering length of the potential to 1, by rescaling space: denoting the scattering length of $v^{(0)}_n$ by $a_n$, we take the potential to be
 \begin{equation}
   v_n(\mathbf x):=v^{(0)}_n\left({\textstyle\frac{|\mathbf x|}{a_n}}\right)
@@ -786,12 +846,15 @@ For smaller densities, for the Simple Equation, we see that the predictions made
 
 \begin{figure}
   \hfil\includegraphics[width=8cm]{hardcore_energy.pdf}
+  \hfil\includegraphics[width=8cm]{hardcore_compare.pdf}
   \hfil\includegraphics[width=8cm]{hardcore_condensate.pdf}
   \caption{
-    The energy (top) and non-condensed fraction (bottom) as a function of the density for the hard core potential.
+    The energy (top), relative error in the energy $\frac{\tilde e-e_{\mathrm{QMC}}}{e_{\mathrm{QMC}}}$ (middle), and non-condensed fraction (bottom) as a function of the density for the hard core potential.
     The circles were computed by solving the hard core Simple Equation for $|\mathbf x|>1$ (simple hc).
     The lines were computed by approximating the hard core potential by the potential $v_{512}(\mathbf x)$, see\-~(\ref{vn}).
     We compare the predictions of the Big, Medium and Simple Equations to QMC results reported\-~\cite{GBC99}.
+    The prediction of Bogolubov theory\-~(\ref{eta0}) is also plotted for comparison (Bog).
+    The right edge of the plots correspond to the close-packing density $\rho_{\mathrm{cp}}=\sqrt2$\-~\cite{Ha05}.
   }
   \label{fig:hardcore}
 \end{figure}
@@ -811,7 +874,7 @@ It is quite easy to find a counter-example if $v$ is not of positive type.
 For instance, if $v(\mathbf x)=0$ for all $|\mathbf x|<1$, then, consider a wavefunction $\psi$ that is smooth and supported on $|\mathbf x_1|,\cdots,|\mathbf x_N|<\frac 12$.
 Since all particles are at a distance that is $<1$, the potential energy of such a wavefunction is 0, and its kinetic energy is $O(N)$.
 Thus, the energy per particle is of order 1, which, for large $\rho$, is $\ll\frac\rho2\int d\mathbf x\ v(\mathbf x)$.
-(Note that a non-trivial, non-negative potential with $v(\mathbf x)$ cannot be of positive type if $v(0)=0$, since the maximumof a positive type function is attained at $0$.)
+(Note that a non-trivial, non-negative potential with $v(\mathbf x)$ cannot be of positive type if $v(0)=0$, since the maximum of a positive type function is attained at $0$.)
 \bigskip
 
 \indent
@@ -825,6 +888,7 @@ This is further confirmed by the computations for the hard core potential, in wh
   \caption{
     The non-condensed fraction as a function of the density for the potential $16e^{-|\mathbf x|}$.
     We compare the predictions of the Big, Medium and Simple Equations to a QMC simulation.
+    The prediction of Bogolubov theory\-~(\ref{eta0}) is also plotted for comparison (Bog).
   }
   \label{fig:condensate16}
 \end{figure}
@@ -832,10 +896,19 @@ This is further confirmed by the computations for the hard core potential, in wh
 \section{Conclusions}\label{sec:conclusions}
 
 \indent
-In this paper we show the astonishing agreement in the predictions of the ground state energy, condensate fraction and correlation function of the interacting Bose gas given by the {\it simplified approach} developed in 1963\-~\cite{Li63} with the values obtained by Quantum Monte-Carlo calculations.
+In this paper we show the good agreement in the predictions of the ground state energy, condensate fraction and correlation function of the repulsive Bose gas given by the {\it simplified approach} developed in 1963\-~\cite{Li63} with the values obtained by Quantum Monte-Carlo calculations, for the potentials $e^{-|\mathbf x|}$ and $16e^{-|\mathbf x|}$.
 The simplified approach was thought to be accurate  only at low densities, in complete agreement with other analyses of the time.
 Here, we show that it is accurate at {\it all} densities.
-This establishes a new paradigm for many body bosonic physics.
+This establishes a new approach to many body bosonic physics.
+Combining this analysis with the exact results in\-~\cite{CJL20,CJL20b} leads us to conjecture that the simplified approach is accurate for any repulsive potential of positive type with a scattering length and an integral that is not too large.
+\bigskip
+
+\indent
+We have discussed three different approximations, the Big, Medium and Simple Equations.
+The Big Equation is the most accurate, but also the most difficult to solve.
+The Medium Equation is obtained by neglecting terms of higher order in $u$, which makes it much more easy to compute with, while remaining rather close to the Big Equation.
+The Simple Equation is then obtained by approximating $g_2(x)v(x)$ by a Dirac-delta function.
+This drastically simplifies the equation, but is also less accurate at intermediate densities (while the low and high densities are still asymptotically exact).
 \bigskip
 
 \indent
@@ -844,6 +917,8 @@ The method provides a promising avenue to approach singular potentials, such as
 In addition, this allows us to approach various physical questions, such as Bose-Einstein condensation, even in the intermediate density regime, away from the dilute and dense limits.
 
 \begin{acknowledgements}
+  We thank two anonymous referees for many helpful comments.
+  E.H.L. thanks the Institute for Advanced study for its  hospitality.
   U.S.~National Science Foundation grants DMS-1764254 (E.A.C.),  DMS-1802170 (I.J.) are gratefully acknowledged.
 \end{acknowledgements}
 
diff --git a/Changelog b/Changelog
index d467512..ab62ec5 100644
--- a/Changelog
+++ b/Changelog
@@ -1,3 +1,18 @@
+v0.2:
+
+  * Added: Comparison of the energy for two different potentials that share
+  the same scattering length and integral (new figure).
+
+  * Added: Lee-Huang-Yang prediction in figures.
+
+  * Fixed: Wrong exponent in definition of correlation function.
+
+  * Added: More precise comparison of the big equation prediction for the hard
+  core potential (new figure).
+
+  * Changed: Miscellaneous minor fixes and clarifications.
+
+
 v0.1:
 
   * Changed: Improved definition of the two-point correlation function for the
diff --git a/bibliography.bib b/bibliography.bib
index bfd0835..3d7f517 100644
--- a/bibliography.bib
+++ b/bibliography.bib
@@ -34,6 +34,8 @@
 2Fphysreva.60.5129},  year = 1999,  month = {dec},  publisher = {American Physical Society ({APS})},  volume = {60},  number = {6},  pages = {5129--5132},  author = {S. Giorgini and J. Boronat and J. Casulleras},  title = {Ground state of a homogeneous Bose gas: A diffusion Monte Carlo calculation},  journal = {Physical Review A} , pages={}}
 @article{GS09, doi = {10.1007/s10955-009-9718-0},  url = {https://doi.org/10.1007
 2Fs10955-009-9718-0},  year = 2009,  month = {apr},  publisher = {Springer Science and Business Media {LLC}},  volume = {135},  number = {5-6},  pages = {915--934},  author = {Alessandro Giuliani and Robert Seiringer},  title = {The Ground State Energy of the Weakly Interacting Bose Gas at High Density},  journal = {Journal of Statistical Physics} }
+@article{Ha05, doi = {10.4007/annals.2005.162.1065}, url = {https://doi.org/10.4007
+2Fannals.2005.162.1065}, year = 2005, month = {nov}, publisher = {Annals of Mathematics, Princeton U}, volume = {162}, number = {3}, pages = {1065--1185}, author = {Thomas Hales}, title = {A proof of the Kepler conjecture}, journal = {Annals of Mathematics} }
 @article{HCe16, doi = {10.1103/physrevb.94.035126},  url = {https://doi.org/10.1103
 2Fphysrevb.94.035126},  year = 2016,  month = {jul},  publisher = {American Physical Society ({APS})},  volume = {94},  number = {3},  author = {Markus Holzmann and Raymond C. Clay and Miguel A. Morales and Norm M. Tubman and David M. Ceperley and Carlo Pierleoni},  title = {Theory of finite size effects for electronic quantum Monte Carlo calculations of liquids and solids},  journal = {Physical Review B} }
 @article{Ja55, doi = {10.1103/physrev.98.1479},  url = {https://doi.org/10.1103
@@ -44,14 +46,14 @@
 2Fphysrevlett.119.143401},  year = 2017,  month = {oct},  publisher = {American Physical Society ({APS})},  volume = {119},  number = {14},  author = {Catherine E. Klauss and Xin Xie and Carlos Lopez-Abadia and Jos{\'{e}} P. D'Incao and Zoran Hadzibabic and Deborah S. Jin and Eric A. Cornell},  title = {Observation of Efimov Molecules Created from a Resonantly Interacting Bose Gas},  journal = {Physical Review Letters} }
 @article{KMe06, doi = {10.1038/nature04626},  url = {https://doi.org/10.1038
 2Fnature04626},  year = 2006,  month = {mar},  publisher = {Springer Science and Business Media {LLC}},  volume = {440},  number = {7082},  pages = {315--318},  author = {T. Kraemer and M. Mark and P. Waldburger and J. G. Danzl and C. Chin and B. Engeser and A. D. Lange and K. Pilch and A. Jaakkola and H.-C. N\"agerl and R. Grimm},  title = {Evidence for Efimov quantum states in an ultracold gas of caesium atoms},  journal = {Nature} }
-@incollection{Kr02, doi = {10.1142/9789812777072_0007},  year = 2002,  month = {aug},  publisher = {World Scientific},  pages = {265--328},  author = {E. Krotscheck},  title = {Theory of Correlated Basis Functions},  booktitle = {Series on Advances in Quantum Many-Body Theory}, volume = {4}, editor = {A. Fabrocini, S. Fantoni, E. Krotscheck}}
+@book{Kr02, doi = {10.1142/9789812777072_0007},  year = 2002,  month = {aug},  publisher = {World Scientific},  pages = {265--328},  author = {E. Krotscheck},  title = {Theory of Correlated Basis Functions},  series = {Series on Advances in Quantum Many-Body Theory}, volume = {4}, editor = {A. Fabrocini, S. Fantoni, E. Krotscheck}}
 @article{LHY57, doi = {10.1103/physrev.106.1135},  url = {https://doi.org/10.1103
 2Fphysrev.106.1135},  year = 1957,  month = {jun},  publisher = {American Physical Society ({APS})},  volume = {106},  number = {6},  pages = {1135--1145},  author = {T. D. Lee and Kerson Huang and C. N. Yang},  title = {Eigenvalues and Eigenfunctions of a Bose System of Hard Spheres and Its Low-Temperature Properties},  journal = {Physical Review} }
 @article{Le29, doi = {10.1007/bf01340138},  url = {https://doi.org/10.1007
 2Fbf01340138},  year = 1929,  month = {nov},  publisher = {Springer Science and Business Media {LLC}},  volume = {56},  number = {11-12},  pages = {778--789},  author = {W. Lenz},  title = {Die Wellenfunktion und Geschwindigkeitsverteilung des entarteten Gases},  journal = {Zeitschrift f\"ur Physik} }
 @article{Li63, doi = {10.1103/physrev.130.2518},  url = {https://doi.org/10.1103
 2Fphysrev.130.2518},  year = 1963,  month = {jun},  publisher = {American Physical Society ({APS})},  volume = {130},  number = {6},  pages = {2518--2528},  author = {Elliott H. Lieb},  title = {Simplified Approach to the Ground-State Energy of an Imperfect Bose Gas},  journal = {Physical Review} }
-@book{Li65, author = {Elliott H. Lieb}, title = {The Bose fluid}, booktitle = {in Lectures in Theoretical Physics}, year = {1965}, pages = {175--224}, publisher = {of Colorado Press} }
+@book{Li65, author = {Elliott H. Lieb}, title = {The Bose fluid},  series = {Lectures in Theoretical Physics}, volume={VIIC}, year = {1965}, pages = {175--224}, publisher = {University of Colorado Press} }
 @article{LL64, doi = {10.1103/physrev.134.a312},  url = {https://doi.org/10.1103
 2Fphysrev.134.a312},  year = 1964,  month = {apr},  publisher = {American Physical Society ({APS})},  volume = {134},  number = {2A},  pages = {A312--A315},  author = {Elliott H. Lieb and Werner Liniger},  title = {Simplified Approach to the Ground-State Energy of an Imperfect Bose Gas. {III}. Application to the One-Dimensional Model},  journal = {Physical Review} }
 @article{LS64, doi = {10.1103/physrev.133.a899},  url = {https://doi.org/10.1103
diff --git a/figs/compare_pots.fig/Makefile b/figs/compare_pots.fig/Makefile
new file mode 100644
index 0000000..9938c4e
--- /dev/null
+++ b/figs/compare_pots.fig/Makefile
@@ -0,0 +1,36 @@
+PROJECTNAME=compare_pots
+
+SIMPLEQ=simpleq
+
+DATS=bigeq_tent.dat bigeq_exp.dat
+PDFS=$(addsuffix .pdf, $(PROJECTNAME))
+TEXS=$(addsuffix .tikz.tex, $(PROJECTNAME))
+
+all: $(PDFS)
+
+$(PDFS): $(DATS)
+	gnuplot $(patsubst %.pdf, %.gnuplot, $@) > $(patsubst %.pdf, %.tikz.tex, $@)
+	pdflatex -jobname $(basename $@) -file-line-error $(patsubst %.pdf, %.tikz.tex, $@)
+
+bigeq_tent.dat:
+	julia -p 8 $(SIMPLEQ)/main.jl -p "N=12;P=8;J=10;maxiter=100;v_a=32;minlrho=-4;nlrho=100;maxlrho=4;eq=bigeq" -U tent -M medeq energy_rho > $@
+
+bigeq_exp.dat:
+	# rescale rho by beta^3~324.860
+	julia -p 8 $(SIMPLEQ)/main.jl -p "N=12;P=8;J=10;maxiter=100;v_a=19.1968;minlrho=-6.512;nlrho=100;maxlrho=1.488;eq=bigeq" -M medeq energy_rho > $@
+
+install: $(PDFS)
+	cp $^ $(INSTALLDIR)/
+
+clean-aux:
+	rm -f $(addsuffix .tikz.tex, $(PROJECTNAME))
+	rm -f $(addsuffix .aux, $(PROJECTNAME))
+	rm -f $(addsuffix .log, $(PROJECTNAME))
+
+clean-dat:
+	rm -f $(DATS)
+
+clean-tex:
+	rm -f $(PDFS)
+
+clean: clean-aux clean-tex
diff --git a/figs/compare_pots.fig/compare_pots.gnuplot b/figs/compare_pots.fig/compare_pots.gnuplot
new file mode 100644
index 0000000..0f8ac58
--- /dev/null
+++ b/figs/compare_pots.fig/compare_pots.gnuplot
@@ -0,0 +1,39 @@
+set ylabel "$\\frac{e}{\\rho}$" norotate
+set xlabel "$\\rho$"
+
+set xtics 1e-4, 100, 10000
+set xtics add ("$10^{-4}$" 0.0001, "$10^{-2}$" 0.01, "$1$" 1.0, "$10^2$" 100, "$10^4$" 10000)
+unset mxtics
+
+set xrange [1e-5:1e4]
+
+#set ytics 0.8, 0.05, 1.1
+set yrange [0:40]
+
+# default output canvas size: 12.5cm x 8.75cm
+set term lua tikz size 8,6 standalone
+
+set key top left box linetype rgbcolor"#999999" width 0.3 height 0.3 spacing 1.3
+
+
+# set linestyle
+set style line 1 linetype rgbcolor "#FF3399" linewidth 2
+set style line 2 linetype rgbcolor "#DC143C" linewidth 2
+set style line 3 linetype rgbcolor "#DC143C" linewidth 2
+set style line 4 linetype rgbcolor "#4B0082" linewidth 2 dashtype "."
+set style line 5 linetype rgbcolor "#DAA520" linewidth 2
+set style line 6 linetype rgbcolor "#999999" linewidth 1.5 dashtype (1,5)
+
+set pointsize 1
+
+set logscale x
+
+a=0.597754
+beta=6.87436
+
+plot \
+ "bigeq_tent.dat" using 1:($2/$1) with lines ls 4 title "big $v_{32}^{(0)}$",\
+ "bigeq_exp.dat" using ($1*beta**3):($2/$1/beta) with lines ls 5 title "big $\\Phi$",\
+ "2021-02-16+tent32.dat" u 1:($2/$1) ls 1 title "QMC $v_{32}^{(0)}$",\
+ "2021-02-16+19.dat" u ($1*beta**3):($2/$1/beta) ls 2 title "QMC $\\Phi$",\
+ 2*pi*a*(1+128./15/sqrt(pi)*sqrt(x*a**3)) ls 6 title "LHY"
diff --git a/figs/condensate.fig/condensate05.gnuplot b/figs/condensate.fig/condensate05.gnuplot
index 66a642e..f03eb92 100644
--- a/figs/condensate.fig/condensate05.gnuplot
+++ b/figs/condensate.fig/condensate05.gnuplot
@@ -8,7 +8,7 @@ unset mxtics
 #set ytics 0.6, 0.1
 set mytics 2
 
-set xrange [4e-7:100]
+set xrange [4e-8:100]
 
 set yrange [:0.06]
 
@@ -24,14 +24,18 @@ set style line 2 linetype rgbcolor "#DC143C" linewidth 2
 set style line 3 linetype rgbcolor "#32CD32" linewidth 2 dashtype "-"
 set style line 4 linetype rgbcolor "#4B0082" linewidth 2
 set style line 5 linetype rgbcolor "#DAA520" linewidth 2
+set style line 6 linetype rgbcolor "#999999" linewidth 1.5 dashtype (1,5)
 
 set pointsize 1
 
 set logscale x
 
+a=0.7666858699084325
+
 plot \
  "simpleq0.5.dat" using 1:($3<1e-5 ? $2 : 1/0) with lines ls 1 title "simple",\
  "mueq0.5.dat" using 1:($3<1e-5 ? $2 : 1/0) with lines ls 3 title "med",\
  "bigeq0.5.dat" using 1:($3<1e-5 ? $2 : 1/0) with lines ls 5 title "big",\
- "2020-08-27+0.5.dat" using 1:(1-$3) with points ls 2 title "QMC"
+ "2020-08-27+0.5.dat" using 1:(1-$3) with points ls 2 title "QMC" ,\
+ 8./(3*sqrt(pi))*sqrt(x*a**3) ls 6 title "Bog"
 
diff --git a/figs/condensate.fig/condensate16.gnuplot b/figs/condensate.fig/condensate16.gnuplot
index b581781..488d036 100644
--- a/figs/condensate.fig/condensate16.gnuplot
+++ b/figs/condensate.fig/condensate16.gnuplot
@@ -10,6 +10,7 @@ set mytics 2
 
 set xrange [1e-7:100]
 
+set yrange [:1.2]
 
 # default output canvas size: 12.5cm x 8.75cm
 set term lua tikz size 8,6 standalone
@@ -23,14 +24,18 @@ set style line 2 linetype rgbcolor "#DC143C" linewidth 2
 set style line 3 linetype rgbcolor "#32CD32" linewidth 2 dashtype "-"
 set style line 4 linetype rgbcolor "#4B0082" linewidth 2
 set style line 5 linetype rgbcolor "#DAA520" linewidth 2
+set style line 6 linetype rgbcolor "#999999" linewidth 1.5 dashtype (1,5)
 
 set pointsize 1
 
 set logscale x
 
+a=3.9270207371831565
+
 plot \
  "simpleq16.dat" using 1:($3<1e-5 ? $2 : 1/0) with lines ls 1 title "simple",\
  "mueq16.dat" using 1:($3<1e-5 ? $2 : 1/0) with lines ls 3 title "med",\
  "bigeq16.dat" using 1:($3<1e-5 ? $2 : 1/0) with lines ls 5 title "big",\
- "2020-08-27+16.dat" using 1:(1-$3) with points ls 2 title "QMC"
+ "2020-08-27+16.dat" using 1:(1-$3) with points ls 2 title "QMC" ,\
+ 8./(3*sqrt(pi))*sqrt(x*a**3) ls 6 title "Bog"
 
diff --git a/figs/energy.fig/energy1.gnuplot b/figs/energy.fig/energy1.gnuplot
index 9667dba..52fc376 100644
--- a/figs/energy.fig/energy1.gnuplot
+++ b/figs/energy.fig/energy1.gnuplot
@@ -24,14 +24,18 @@ set style line 2 linetype rgbcolor "#DC143C" linewidth 2
 set style line 3 linetype rgbcolor "#32CD32" linewidth 2 dashtype "-"
 set style line 4 linetype rgbcolor "#4B0082" linewidth 2
 set style line 5 linetype rgbcolor "#DAA520" linewidth 2
+set style line 6 linetype rgbcolor "#999999" linewidth 1.5 dashtype (1,5)
 
 set pointsize 1
 
 set logscale x
 
+a=1.2543564105910647
+
 plot \
  "simpleq1.dat" using 1:($3<1e-5 ? ($2/$1) : 1/0) with lines ls 1 title "simple" ,\
  "mueq1.dat" using 1:($3<1e-5 ? ($2/$1) : 1/0) with lines ls 3 title "med",\
  "bigeq1.dat" using 1:($3<1e-5 ? ($2/$1) : 1/0) with lines ls 5 title "big",\
- "2020-10-15+energy.dat" using 1:($2/$1) with points ls 2 title "QMC"
+ "2020-10-15+energy.dat" using 1:($2/$1) with points ls 2 title "QMC" ,\
+ 2*pi*a*(1+128./(15*sqrt(pi))*sqrt(x*a**3)) ls 6 title "LHY"
 
diff --git a/figs/energy.fig/energy16.gnuplot b/figs/energy.fig/energy16.gnuplot
index 4fd8de6..ee6448e 100644
--- a/figs/energy.fig/energy16.gnuplot
+++ b/figs/energy.fig/energy16.gnuplot
@@ -24,14 +24,18 @@ set style line 2 linetype rgbcolor "#DC143C" linewidth 2
 set style line 3 linetype rgbcolor "#32CD32" linewidth 2 dashtype "-"
 set style line 4 linetype rgbcolor "#4B0082" linewidth 2
 set style line 5 linetype rgbcolor "#DAA520" linewidth 2
+set style line 6 linetype rgbcolor "#999999" linewidth 1.5 dashtype (1,5)
 
 set pointsize 1
 
 set logscale x
 
+a=3.9270207371831565
+
 plot \
  "simpleq16.dat" using 1:($3<1e-5 ? ($2/$1) : 1/0) with lines ls 1 title "simple",\
  "mueq16.dat" using 1:($3<1e-5 ? ($2/$1) : 1/0) with lines ls 3 title "med",\
  "bigeq16.dat" using 1:($3<1e-5 ? ($2/$1) : 1/0) with lines ls 5 title "big",\
- "2020-08-27+16.dat" using 1:($2/$1) with points ls 2 title "QMC"
+ "2020-08-27+16.dat" using 1:($2/$1) with points ls 2 title "QMC" ,\
+ 2*pi*a*(1+128./(15*sqrt(pi))*sqrt(x*a**3)) ls 6 title "LHY"
 
diff --git a/figs/energy_relative_error.fig/Makefile b/figs/energy_relative_error.fig/Makefile
index 8c75f0f..d4dfb0e 100644
--- a/figs/energy_relative_error.fig/Makefile
+++ b/figs/energy_relative_error.fig/Makefile
@@ -1,6 +1,6 @@
 PROJECTNAME=cmp_energy16 cmp_energy1
 
-SIMPLEQ=simpleq
+SIMPLEQ=~/Work/Research/2018+bose_gas/cmp/simpleq
 
 DATS= simpleq16.dat mueq16.dat bigeq16.dat simpleq1.dat mueq1.dat bigeq1.dat
 PDFS=$(addsuffix .pdf, $(PROJECTNAME))
@@ -13,42 +13,28 @@ $(PDFS): $(DATS)
 	pdflatex -jobname $(basename $@) -file-line-error $(patsubst %.pdf, %.tikz.tex, $@)
 
 bigeq16.dat:
-	for rho in 1e-3 5e-3 1e-2 2e-2 1e-1 1e+0; do \
-	  echo -n $$rho "  " >> $@-tmp ;\
-	  julia $(SIMPLEQ)/main.jl -p "eq=bigeq;N=12;P=8;J=10;maxiter=100;v_a=16;rho=$$rho;minlrho=-6;nlrho=50" -M medeq energy >> $@-tmp ;\
-	done
+	julia $(SIMPLEQ)/main.jl -p "eq=bigeq;N=12;P=8;J=10;maxiter=100;v_a=16;minlrho=-6;nlrho=50;rhos=1e-3,5e-3,1e-2,2e-2,1e-1,1e+0" -M medeq energy > $@-tmp
 	paste $@-tmp 2020-08-27+16.dat > $@
 	rm $@-tmp
 simpleq16.dat:
-	for lrho in -3 -2.3010299956639813 -2 -1.6989700043360187 -1 0; do \
-	  julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=16;minlrho=$$lrho;nlrho=1" -M simpleq energy_rho >> $@-tmp ;\
-	done
+	julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=16;rhos=1e-3,5e-3,1e-2,2e-2,1e-1,1e+0" -M simpleq energy_rho > $@-tmp
 	paste $@-tmp 2020-08-27+16.dat > $@
 	rm $@-tmp
 mueq16.dat:
-	for lrho in -3 -2.3010299956639813 -2 -1.6989700043360187 -1 0; do \
-	  julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=16;minlrho=$$lrho;nlrho=1" -M mueq energy_rho >> $@-tmp ;\
-	done
+	julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=16;rhos=1e-3,5e-3,1e-2,2e-2,1e-1,1e+0" -M mueq energy_rho > $@-tmp
 	paste $@-tmp 2020-08-27+16.dat > $@
 	rm $@-tmp
 
 bigeq1.dat:
-	for rho in 1e-6 1e-4 1e-3 5e-3 1e-2 1.5e-2 2e-2 5e-2 1e-1 1.0 1e1 5e1; do \
-	  echo -n $$rho "  " >> $@-tmp ;\
-	  julia $(SIMPLEQ)/main.jl -p "N=12;P=8;J=10;maxiter=100;v_a=1;rho=$$rho;minlrho=-6;nlrho=50;eq=bigeq" -M medeq energy >> $@-tmp ;\
-	done
+	julia $(SIMPLEQ)/main.jl -p "N=12;P=8;J=10;maxiter=100;v_a=1;minlrho=-6;nlrho=50;eq=bigeq;rhos=1e-6,1e-4,1e-3,5e-3,1e-2,1.5e-2,2e-2,5e-2,1e-1,1.0,1e1,5e1" -M medeq energy > $@-tmp
 	paste $@-tmp 2020-10-15+energy.dat > $@
 	rm $@-tmp
 simpleq1.dat:
-	for lrho in -6 -4 -3 -2.3010299956639813 -2 -1.8239087409443189 -1.6989700043360187 -1.3010299956639813 -1 0 1 1.6989700043360187; do\
-	  julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=1;minlrho=$$lrho;nlrho=1" -M simpleq energy_rho >> $@-tmp ;\
-	done
+	julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=1;rhos=1e-6,1e-4,1e-3,5e-3,1e-2,1.5e-2,2e-2,5e-2,1e-1,1.0,1e1,5e1" -M simpleq energy_rho > $@-tmp
 	paste $@-tmp 2020-10-15+energy.dat > $@
 	rm $@-tmp
 mueq1.dat:
-	for lrho in -6 -4 -3 -2.3010299956639813 -2 -1.8239087409443189 -1.6989700043360187 -1.3010299956639813 -1 0 1 1.6989700043360187; do\
-	  julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=1;minlrho=$$lrho;nlrho=1" -M mueq energy_rho >> $@-tmp ;\
-	done
+	julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=1;rhos=1e-6,1e-4,1e-3,5e-3,1e-2,1.5e-2,2e-2,5e-2,1e-1,1.0,1e1,5e1" -M mueq energy_rho >> $@-tmp ;\
 	paste $@-tmp 2020-10-15+energy.dat > $@
 	rm $@-tmp
 
diff --git a/figs/hardcore.fig/Makefile b/figs/hardcore.fig/Makefile
index 1cf8737..d8e1dfb 100644
--- a/figs/hardcore.fig/Makefile
+++ b/figs/hardcore.fig/Makefile
@@ -1,7 +1,7 @@
-PROJECTNAME=hardcore_energy hardcore_condensate
+PROJECTNAME=hardcore_energy hardcore_compare hardcore_condensate
 SIMPLEQ=simpleq
 
-DATS=simpleq-energy.dat simpleq_soft-energy.dat mueq_soft-energy.dat bigeq_soft-energy.dat simpleq-condensate.dat simpleq_soft-condensate.dat mueq_soft-condensate.dat bigeq_soft-condensate.dat
+DATS=simpleq-energy.dat simpleq_soft-energy.dat mueq_soft-energy.dat bigeq_soft-energy.dat simpleq-compare.dat simpleq_soft-compare.dat mueq_soft-compare.dat bigeq_soft-compare.dat simpleq-condensate.dat simpleq_soft-condensate.dat mueq_soft-condensate.dat bigeq_soft-condensate.dat
 
 PDFS=$(addsuffix .pdf, $(PROJECTNAME))
 TEXS=$(addsuffix .tikz.tex, $(PROJECTNAME))
@@ -21,6 +21,23 @@ mueq_soft-energy.dat:
 bigeq_soft-energy.dat:
 	julia -p 8 $(SIMPLEQ)/main.jl -p "N=12;P=8;eq=bigeq;maxiter=100;maxlrho=0.15;nlrho=64;v_a=512;v_b=1.2494804800401416" -U tent -M medeq energy_rho_init_prevrho > $@
 
+simpleq-compare.dat:
+	julia -p 8 $(SIMPLEQ)/main.jl -p "N=12;P=8;J=5;rhos=1e-6,5e-6,1e-5,5e-5,1e-4,5e-4,1e-3,5e-3,1e-2,5e-2,1e-1,0.166,0.244" -M simpleq-hardcore energy_rho > $@-tmp
+	paste $@-tmp GBC99-energy.dat > $@
+	rm $@-tmp
+simpleq_soft-compare.dat:
+	julia $(SIMPLEQ)/main.jl -p "maxiter=100;v_a=328.77174503097945;v_b=1.2479228505627133;rhos=1e-6,5e-6,1e-5,5e-5,1e-4,5e-4,1e-3,5e-3,1e-2,5e-2,1e-1,0.166,0.244" -U tent -M simpleq energy_rho > $@-tmp
+	paste $@-tmp GBC99-energy.dat > $@
+	rm $@-tmp
+mueq_soft-compare.dat:
+	julia $(SIMPLEQ)/main.jl -p "maxiter=100;v_a=328.77174503097945;v_b=1.2479228505627133;rhos=1e-6,5e-6,1e-5,5e-5,1e-4,5e-4,1e-3,5e-3,1e-2,5e-2,1e-1,0.166,0.244" -U tent -M mueq energy_rho > $@-tmp
+	paste $@-tmp GBC99-energy.dat > $@
+	rm $@-tmp
+bigeq_soft-compare.dat:
+	julia -p 8 $(SIMPLEQ)/main.jl -p "N=12;P=8;eq=bigeq;maxiter=100;v_a=328.77174503097945;v_b=1.2479228505627133;rhos=1e-6,5e-6,1e-5,5e-5,1e-4,5e-4,1e-3,5e-3,1e-2,5e-2,1e-1,0.166,0.244" -U tent -M medeq energy_rho_init_prevrho > $@-tmp
+	paste $@-tmp GBC99-energy.dat > $@
+	rm $@-tmp
+
 simpleq-condensate.dat:
 	julia -p 8 $(SIMPLEQ)/main.jl -p "N=12;P=8;J=5;maxlrho=0.15;nlrho=24" -M simpleq-hardcore condensate_fraction_rho > $@
 simpleq_soft-condensate.dat:
diff --git a/figs/hardcore.fig/hardcore_compare.gnuplot b/figs/hardcore.fig/hardcore_compare.gnuplot
new file mode 100644
index 0000000..ef67c1e
--- /dev/null
+++ b/figs/hardcore.fig/hardcore_compare.gnuplot
@@ -0,0 +1,36 @@
+set ylabel "energy relative error"
+set xlabel "$\\rho$"
+
+set xtics 1e-6, 100, 1
+set xtics add ("$10^{-6}$" 0.000001, "$10^{-4}$" 0.0001, "$10^{-2}$" 0.01, "$1$" 1.0)
+unset mxtics
+
+set xrange [1e-9:sqrt(2)]
+
+set ytics 1e-5, 10, 10
+set ytics add ("$10^{-5}$" 1e-5, "$10^{-4}$" 1e-4, "$10^{-3}$" 1e-3, "$10^{-2}$" 1e-2, "$10^{-1}$" 1e-1, "$10^{0}$" 1, "$10^1$" 10)
+#unset mytics
+
+# default output canvas size: 12.5cm x 8.75cm
+set term lua tikz size 8,6 standalone
+
+set key top left box linetype rgbcolor"#999999" width 0.3 height 0.3 spacing 1.3
+
+# set linestyle
+set style line 1 linetype rgbcolor "#4169E1" linewidth 2 dashtype "."
+set style line 2 linetype rgbcolor "#DC143C" linewidth 2
+set style line 3 linetype rgbcolor "#32CD32" linewidth 2 dashtype "-"
+set style line 4 linetype rgbcolor "#4B0082" linewidth 2
+set style line 5 linetype rgbcolor "#DAA520" linewidth 2
+set style line 6 linetype rgbcolor "#999999" linewidth 2
+
+set pointsize 1
+
+set logscale xy
+
+plot \
+ "simpleq-compare.dat" using 1:($3<1e-2 ? 2*abs($5/2-$2)/$5 : 1/0) with lp ls 4 title "simple hc",\
+ "simpleq_soft-compare.dat" using 1:($3<1e-2 ? 2*abs($5/2-$2)/$5 : 1/0) with lp ls 1 title "simple",\
+ "mueq_soft-compare.dat" using 1:($3<1e-2 ? 2*abs($5/2-$2)/$5 : 1/0) with lp ls 3 title "med",\
+ "bigeq_soft-compare.dat" using 1:($3<1e-2 ? 2*abs($5/2-$2)/$5 : 1/0) with lp ls 5 title "big"
+
diff --git a/figs/hardcore.fig/hardcore_condensate.gnuplot b/figs/hardcore.fig/hardcore_condensate.gnuplot
index 974c9fd..7cd2fd2 100644
--- a/figs/hardcore.fig/hardcore_condensate.gnuplot
+++ b/figs/hardcore.fig/hardcore_condensate.gnuplot
@@ -24,16 +24,19 @@ set style line 2 linetype rgbcolor "#DC143C" linewidth 2
 set style line 3 linetype rgbcolor "#32CD32" linewidth 2 dashtype "-"
 set style line 4 linetype rgbcolor "#4B0082" linewidth 2
 set style line 5 linetype rgbcolor "#DAA520" linewidth 2
-set style line 6 linetype rgbcolor "#999999" linewidth 2
+set style line 6 linetype rgbcolor "#999999" linewidth 1.5 dashtype (1,5)
 
 set pointsize 1
 
 set logscale x
 
+a=1.
+
 plot \
- "simpleq-condensate.dat" using 1:($3<1e-2 ? $2 : 1/0) ls 6 title "simple hc",\
+ "simpleq-condensate.dat" using 1:($3<1e-2 ? $2 : 1/0) ls 4 title "simple hc",\
  "simpleq_soft-condensate.dat" using 1:($3<1e-2 ? $2 : 1/0) with lines ls 1 title "simple" ,\
  "mueq_soft-condensate.dat" using 1:($3<1e-2 ? $2 : 1/0) with lines ls 3 title "med",\
  "bigeq_soft-condensate.dat" using 1:($3<1e-2 ? $2 : 1/0) with lines ls 5 title "big",\
- "GBC99-condensate.dat" using 1:(1-$2) ls 2 title "QMC [15]"
+ "GBC99-condensate.dat" using 1:(1-$2) ls 2 title "QMC [15]" ,\
+ 8./(3*sqrt(pi))*sqrt(x*a**3) ls 6 title "Bog"
 
diff --git a/figs/hardcore.fig/hardcore_energy.gnuplot b/figs/hardcore.fig/hardcore_energy.gnuplot
index 9967232..d478518 100644
--- a/figs/hardcore.fig/hardcore_energy.gnuplot
+++ b/figs/hardcore.fig/hardcore_energy.gnuplot
@@ -23,16 +23,19 @@ set style line 2 linetype rgbcolor "#DC143C" linewidth 2
 set style line 3 linetype rgbcolor "#32CD32" linewidth 2 dashtype "-"
 set style line 4 linetype rgbcolor "#4B0082" linewidth 2
 set style line 5 linetype rgbcolor "#DAA520" linewidth 2
-set style line 6 linetype rgbcolor "#999999" linewidth 2
+set style line 6 linetype rgbcolor "#999999" linewidth 1.5 dashtype (1,5)
 
 set pointsize 1
 
 set logscale xy
 
+a=1.
+
 plot \
- "simpleq-energy.dat" using 1:($3<1e-2 ? $2/$1 : 1/0) ls 6 title "simple hc",\
+ "simpleq-energy.dat" using 1:($3<1e-2 ? $2/$1 : 1/0) ls 4 title "simple hc",\
  "simpleq_soft-energy.dat" using 1:($3<1e-2 ? $2/$1 : 1/0) with lines ls 1 title "simple",\
  "mueq_soft-energy.dat" using 1:($3<1e-2 ? $2/$1 : 1/0) every ::::95 with lines ls 3 title "med",\
  "bigeq_soft-energy.dat" using 1:($3<1e-2 ? $2/$1 : 1/0) with lines ls 5 title "big",\
- "GBC99-energy.dat" using 1:($2/(2*$1)) ls 2 title "QMC [15]"
+ "GBC99-energy.dat" using 1:($2/(2*$1)) ls 2 title "QMC [15]" ,\
+ 2*pi*a*(1+128./(15*sqrt(pi))*sqrt(x*a**3)) ls 6 title "LHY"