In dense environments, like neutron star mergers, dynamical backcoupling between the strong and weak interactions becomes important. We therefore use non-perturbative functional methods to investigate the β -decay as the dominant dynamical process. Especially, we use both Dyson-Schwinger/Bethe-Salpeter equations and Functional-Renormalization-Group equations to control systematic errors. We present the necessary setup and outline relevant challenges.