Toward Reliable $0νββ$-decay Nuclear Matrix Elements: Exploring the potential of measuring $γγ$-transitions

In this work, a physics process known since quite long ago, double-gamma decay ($γγ$), has been revisited from a new perspective: providing valuable insights into neutrinoless double-beta decay ($0νββ$) nuclear matrix elements. At the same time that an eager experimental search of $0νββ$ decay is underway, the nuclear and particle physics communities have made huge progress during the last years. The main goal of this thesis has been to investigate one of those approaches which is the computation of nuclear observables related to $0νββ$ decay as a way to help in determining and reducing the theoretical uncertainties in $0νββ$-decay NMEs. This way, we have proposed that the measurement of the double magnetic dipole $γγ$ decay from the double isobaric analog state in the $ββ$-decay final nucleus could establish the value and reduce the uncertainty in $0νββ$-decay NMEs, because the NMEs of the two processes are very well correlated. We have explored the validity of this approach to predict and quantify NMEs uncertainties in the case where data is available for the nuclear observable related to $0νββ$ decay, that is for the Standard Model allowed $2νββ$ decay. A further objective has been to start with the theoretical characterization of the first steps toward the measurement of the proposed $γγ$ double magnetic dipole decay from the double isobaric analogue state. In particular, we have studied the main decaying modes that can compete with this process: single gamma decay and proton emission, vi and we have calculated the corresponding branching ratios. In addition, we have also given the first steps in the study of the relation between $γγ$ and $0νββ$-decay NMEs in the ab initio valence space in-medium similarity renormalization group.