  Large-scale models of nuclear structure are currently the only way to provide
consistent datasets for the many properties of thousands of exotic nuclei that
are required by nucleosynthesis simulations. In [W.Ryssens et al., Eur. Phys.
J. A 58, 246 (2022)], we recently presented the new BSkG2 model based on an
energy density functional of the Skyrme type. Relying on a flexible
three-dimensional coordinate representation of the nucleus, the model takes
into account both triaxial deformation and time-reversal symmetry breaking.
BSkG2 achieves a state-of-the-art global description of nuclear ground state
(g.s.) properties and reproduces in particular the known masses with a
root-mean-square (rms) deviation of 678 keV. Moving beyond g.s. properties, the
model also reproduces all empirical values for the primary and secondary
barriers as well as isomer excitation energies of actinide nuclei with rms
deviations below 500 keV, i.e. with unprecedented accuracy. Here we discuss in
detail the extension of our framework to the calculation of the fission
barriers of 45 actinide nuclei, including odd-mass and odd-odd systems. We
focus in particular on the impact of symmetry breaking which is key to the
accuracy of the model: we allow systematically for axial, reflection and
time-reversal symmetry breaking. The effect of the latter on the fission
properties of odd-mass and odd-odd nuclei is small, but we find that allowing
for shapes with triaxial or octupole deformation, as well as shapes with both,
is crucial to achieving this accuracy. The numerical accuracy of our coordinate
space approach, the variety of nuclear configurations explored and the
simultaneous successful description of fission properties and known masses
makes BSkG2 the tool of choice for the large-scale study of nuclear structure.
