Atomic nuclei are one of the most multifaceted systems in the universe. The approach of studying these complex systems through the use of global nuclear models has been a cornerstone of LANL and proven insightful to a wide variety of applications. We report on recent progress made by combining nuclear structure calculations, e.g. from Quasi-particle Random Phase Approximation (QRPA), using the Finite-Range Droplet Model (FRDM) with statistical Hauser-Feshbach (HF). The QRPA+HF approach allows for the description of gamma-ray competition at each stage of particle emission. We apply this model to neutron-rich nuclei which yields branching ratios for beta-delayed neutron emission and beta-delayed fission. We further explore a low energy M1 scissors mode correlated to nuclear deformation and show how this improves the match between theoretical calculations and evaluated data. We extend the results of these models across the chart of nuclides and comment on the potential impact to astrophysical applications.