Suppose the coupling between two spins is $C_{i,j}<0$, then the classical partition function is given by $$Z=\sum_{\{s_i\}}e^{\sum_{i,j}s_iK_{ij}s_j+h\sum_{i}s_i}$$ where $K_{ij}=-\beta C_{ij}$ and $h=\beta H$ After Hubbard-Stratonovich transformation and some simple manipulations we get $$Z=\frac{1}{\sqrt{\mathrm{det}(4\pi K)}}\int D\phi e^{-\sum_{i,j}\phi_iK_{ij}\phi_j+h\sum_{i}\phi_i+\sum_{i}\mathrm{ln}(\mathrm{cosh}(2\sum_j K_{ij}\phi_j))}$$ where the spin configurations have been summed over.
When trying to apply low $T$ perturbation, $K_{ij}$ is large and positive. Then normally people just assume that strong fluctuations will be suppressed, or more precisely it is assumed that $|\phi_i|<<1$ and that the spatial profile of the field is smooth. Then based on these assumptions, the standard Landau mean-field expansion with a ${\phi}^4$ term can be obtained.
My question is: After the transformation, the variables $\phi_i$ are supposed to be randomly fluctuating numbers that can take arbitrary values. However,the matrix $K$ is obviously not positive definite and as a result some large $\phi$ can contribute significantly to the integration instead of being suppressed. Furthermore, based on this same argument, there seems to be no reason for the fluctuating field to be smooth. Where does my argument go wrong?
