Lattice Field Theory can be used to study finite temperature first-order phase transitions in new, strongly-coupled gauge theories of phenomenological interest. Metastable dynamics arising in proximity of the phase transition can lead to large, uncontrolled numerical errors when analysed with standard methods. In this contribution, we discuss a prototype lattice calculation in which the first-order deconfinement transition in the strong Yang-Mills sector of the standard model is analysed using a novel lattice method, the logarithmic linear relaxation algorithm. This method provides a determination of the density of states of the system with exponential error suppression.