Bose Einstein Distribution (3D)

Surface plot of the Bose-Einstein distribution over the complex plane. This plot shows an important feature of the Matsubara formalism developed for QFT at non-zero temperature. It’s a method to calculate expectation values of operators in a canonical ensemble evolved by an imaginary time. In momentum space, this leads to the replacement of continuous frequencies by discrete imaginary (Matsubara) frequencies. Away from the imaginary axis, the distribution becomes approximately flat, particularly at sufficiently low temperatures. Pulled from arxiv:1712.09863.

Edit and compile if you like:

% Surface plot of the Bose-Einstein distribution over the complex plane.
% This plot shows an important feature of the Matsubara formalism developed for QFT at non-zero temperature. It's a method to calculate expectation values of operators in a canonical ensemble evolved by an imaginary time. In momentum space, this leads to the replacement of continuous frequencies by discrete imaginary (Matsubara) frequencies.
% Away from the imaginary axis, the distribution becomes approximately flat, particularly at sufficiently low temperatures.
% Pulled from https://arxiv.org/abs/1712.09863.

\documentclass{standalone}

\usepackage{pgfplots,mathtools}
\pgfplotsset{compat=newest}

\let\Re\relax
\DeclareMathOperator{\Re}{Re}
\let\Im\relax
\DeclareMathOperator{\Im}{Im}

\begin{document}
\begin{tikzpicture}
  \begin{axis}[
      xlabel=$\Re(p_0)$,
      ylabel=$\Im(p_0)$,
      zlabel=$n_\text{B}(p_0)$,
      x label style={at={(0.35, 0)}},
      y label style={at={(0.95, 0.15)}},
      shader=flat,
      tickwidth=0pt
    ]

    \def\nB{(e^(2*x) - 2*e^x*cos(deg(y)) + 1)^(-1/2)}

    \addplot3[surf, opacity=0.5, domain=1:10, y domain=-10:10]{\nB};

    \addplot3[surf, opacity=0.5, domain=-10:-2, y domain=-10:10]{\nB};

    \addplot3[surf, opacity=0.5, domain=-2:1, y domain=-10:10, restrict z to domain*=0:2]{\nB};

  \end{axis}
\end{tikzpicture}
\end{document}