We present a novel random field model based on maximum entropy distributions and fourth order moments to capture the steady state distribution of a passive scalar in isotropic homogeneous turbulence. We show that the model's parameters are able to be efficiently learned through a pseudolikelihood approach. Empirical results show that the model performs better than empirical and Gaussian baselines in measuring important statistics of the passive scalar.