Unified Field-integral Thermodynamics of Bose Mixtures: Stability and Critical Behavior

We establish a unified thermodynamic framework for Bose mixtures at finite temperatures based on the functional field integral, within which the decision on whether to discard the anomalous densities, when determining the density configuration and stability matrix, yields distinct theories. Beyond the existing Hartree-Fock approximation and Ota-Giorgini-Stringari theory, retaining the anomalous densities throughout will provide a completely self-consistent thermodynamic description, requiring the combination of the Hartree-Fock-Bogoliubov approximation and the representative statistical ensemble. Comparing three approaches for predicting magnetic susceptibility, we highlight the role of anomalous densities in stabilizing superfluid mixtures. We further unveil that Feshbach coupling can either expand the regime of atomic and molecular superfluids, or induce a phase transition to a pure molecular superfluid, depending on their density ratio. Importantly, we show that thermal fluctuations will trigger a phase transition from stable to unstable mixtures, where anomalous densities can serve as distinct signatures for experimental observation.