  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.
