I wanted to plot E2(T) fot a constant composition, or constant chemical potential, T-scans, but I keep getting E2=0.
How do I turn it on?
Sorry I just first want to check for a possible simple mistake:
E2 should be in column 6 in both canonical and grandcanonical modes.
It is the susceptibility (column 7) that turns off in canonical mode (because there are no composition fluctuations).
If that is not the problem, try increasing the -sigdig parameter to output more digits (in case the fluctuations are too small).