Modeling the ferroelectric phase transition in barium titanate with DFT accuracy and converged sampling