The serpent Monte Carlo code and the serpent-ants two-step calculation chain are used to model the hot zero power physics tests described in the BEAVRS benchmark. The predicted critical boron concentrations, control rod group worths, and isothermal temperature coefficients are compared between serpent and serpent-ants as well as against the experimental measurements. Furthermore, radial power distributions in the unrodded and rodded core configurations are compared between serpent and serpent-ants. In addition to providing results using a best practices calculation chain, the effects of several simplifications or omissions in the group constant generation process on the results are estimated. Both the direct and two-step neutronics solutions provide results close to the measured values. Comparison between the measured data and the direct serpent Monte Carlo solution yields root-mean-square (RMS) differences of 12.1 mg/kg, and for boron, control rod worths and temperature coefficients, respectively. The two-step serpent-ants solution reaches a similar level of accuracy with RMS differences of 17.4 mg/kg, and . The match in the radial power distribution between serpent and serpent-ants was very good with the RMS and MAXimum (MAX) for pin power errors being 1.31% and 4.99%, respectively, in the unrodded core and 1.67% (RMS) and 8.39% (MAX) in the rodded core.