Search

We investigate the experimental feasibility of detecting second-order double-magnetic dipole ($\gamma\gamma$-$M1M1$) decays from double isobaric analog states (DIAS), which have recently been found to be strongly correlated with the nuclear matrix elements of neutrinoless $\beta\beta$ decay. Using the nuclear shell model, we compute theoretical branching ratios for $\gamma\gamma$-$M1M1$ decays and compare them with other competing processes, such as single-$\gamma$ decay and proton emission, which represent the dominant decay channels. We also estimate the potential competition from internal conversion and internal pair creation, which can influence the decay dynamics. Additionally, we propose an experimental strategy based on using LaBr$_3$ scintillators to identify $\gamma\gamma$-$M1M1$ transitions from the DIAS amidst the background of the competing processes. Our approach emphasizes the challenges of isolating the rare $\gamma\gamma$-$M1M1$ decay and suggests ways to enhance the experimental detection sensitivity. Our simulations suggest that it may be possible to access experimentally $\gamma\gamma$-$M1M1$ decays from DIAS, shedding light on the neutrinoless $\beta\beta$ decay nuclear matrix elements., Comment: 8 pages, 4 figures, 5 tables

Simulating physical systems with variational quantum algorithms is a well-studied approach, but it is challenging to implement in current devices due to demands in qubit number and circuit depth. We show how limited knowledge of the system, namely the entropy of its subsystems, its entanglement structure or certain symmetries, can be used to reduce the cost of these algorithms with entanglement forging. To do so, we simulate a Fermi-Hubbard one-dimensional chain with a parametrized hopping term, as well as atomic nuclei ${}^{28}$Ne and ${}^{60}$Ti with the nuclear shell model. Using an adaptive variational quantum eigensolver we find significant reductions in both the maximum number of qubits (up to one fourth) and the amount of two-qubit gates (over an order of magnitude) required in the quantum circuits. Our findings indicate that our method, entropy-driven entanglement forging, can be used to adjust quantum simulations to the limitations of noisy intermediate-scale quantum devices.

We evaluate all nuclear matrix elements (NMEs) up to next-to-next-to leading order (N$^2$LO) in chiral effective field theory ($\chi$EFT) for the neutrinoless double-beta ($0\nu\beta\beta$) decay of the nuclei most relevant for experiments, including $^{76}$Ge, $^{100}$Mo, and $^{136}$Xe. We use the proton-neutron quasiparticle random-phase approximation (pnQRPA) and the nuclear shell model to calculate the N$^2$LO NMEs from very low-momentum (ultrasoft) neutrinos and from loop diagrams usually neglected in $0\nu\beta\beta$ studies. Our results indicate that, separately, the new N$^2$LO NMEs can have up to a $15\%$ impact on the total NME, but cancellations lead to an overall N$^2$LO contribution centered around $-5\%$ for the shell model and $+10\%$ for the pnQRPA, with sizable uncertainties due to the short-range nature of the loop NMEs. The sign discrepancy between many-body methods is common to all studied nuclei and points to the different behaviour of the intermediate states of the $0\nu\beta\beta$ decay. In addition, our results support the $\chi$EFT prediction that the ultrasoft NME is closely related to contributions usually referred to as ``beyond the closure approximation''., Comment: 10 pages, 4 figures, and 2 tables

We study the shape coexistence in the nucleus $^{28}$Si with the nuclear shell model using numerical diagonalizations complemented with variational calculations based on the projected generator-coordinate method. The theoretical electric quadrupole moments and transitions as well as the collective wavefunctions indicate that the standard USDB interaction in the $sd$ shell describes well the ground-state oblate rotational band, but misses the experimental prolate band. Guided by the quasi-SU(3) model, we show that the prolate band can be reproduced in the $sd$ shell by reducing the energy of the $0d_{3/2}$ orbital. Alternatively, in the extended $sdpf$ configuration space a modification of the SDPF-NR interaction that accommodates cross-shell excitations also reproduces the oblate and prolate bands. Finally, we address the possibility of superdeformation in $^{28}$Si within the $sdpf$ space. Our results indicate that superdeformed structures appear at about $18$-$20$~MeV., Comment: 15 pages, 11 figures. Improved version following the reviewer's comments

The compact muon solenoid (CMS) experiment is a general-purpose detector for high-energy collision at the large hadron collider (LHC) at CERN. It employs an online data quality monitoring (DQM) system to promptly spot and diagnose particle data acquisition problems to avoid data quality loss. In this study, we present semi-supervised spatio-temporal anomaly detection (AD) monitoring for the physics particle reading channels of the hadronic calorimeter (HCAL) of the CMS using three-dimensional digi-occupancy map data of the DQM. We propose the GraphSTAD system, which employs convolutional and graph neural networks to learn local spatial characteristics induced by particles traversing the detector, and global behavior owing to shared backend circuit connections and housing boxes of the channels, respectively. Recurrent neural networks capture the temporal evolution of the extracted spatial features. We have validated the accuracy of the proposed AD system in capturing diverse channel fault types using the LHC Run-2 collision data sets. The GraphSTAD system has achieved production-level accuracy and is being integrated into the CMS core production system--for real-time monitoring of the HCAL. We have also provided a quantitative performance comparison with alternative benchmark models to demonstrate the promising leverage of the presented system., Comment: 23 pages, 15 figures, 3 tables

We investigate neutrinoless double-beta decay ($0\nu\beta\beta$) in the minimal extension of the standard model of particle physics, the $\nu$SM, where gauge-singlet right-handed neutrinos give rise to Dirac and Majorana neutrino mass terms. We focus on the associated sterile neutrinos and argue that the usual evaluation of their contributions to $0\nu\beta\beta$, based on mass-dependent nuclear matrix elements, is missing important contributions from neutrinos with ultrasoft and hard momenta. We identify the hadronic and nuclear matrix elements that enter the new contributions, and calculate all relevant nuclear matrix elements for $^{136}$Xe using the nuclear shell model. Finally, we illustrate the impact on $0\nu\beta\beta$ rates in specific neutrino mass models and show that the new contributions significantly alter the $0\nu\beta\beta$ rate in most parts of the $\nu$SM parameter space., Comment: 8 pages, 3 figures, 2 tables

We explore the relation between the nuclear matrix elements of neutrinoless double-beta ($0\nu\beta\beta$) decay and two other processes: double Gamow-Teller (DGT) and double-magnetic dipole (M1M1) transitions, with focus on medium-mass to heavy nuclei studied with the proton-neutron quasiparticle random-phase approximation (pnQRPA) framework. We explore a wide span of isoscalar proton-neutron pairing strengths covering the typical range of values that describe well $\beta$- and two-neutrino $\beta\beta$-decay data. Our results indicate good linear correlations between $0\nu\beta\beta$ and both DGT and M1M1 matrix elements. Together with future measurements of DGT and M1M1 transitions, these correlations could help constrain the values of the $0\nu\beta\beta$-decay nuclear matrix elements., Comment: 14 pages, 12 figures, 2 tables, corresponds to the published version

We calculate the nuclear matrix element for the two-neutrino $\beta\beta$ decay of $^{136}$Xe into the first excited $0^+$ state of $^{136}$Ba. We use different many-body methods: the quasiparticle random-phase approximation (QRPA) framework, the nuclear shell model, the interacting boson model (IBM-2), and an effective field theory (EFT) for $\beta$ and $\beta\beta$ decays. While the QRPA suggests a decay rate at the edge of current experimental limits, the shell model points to a half-life about two orders of magnitude longer. The predictions of the IBM-2 and the EFT lie in between, and the latter provides systematic uncertainties at leading order. An analysis of the running sum of the nuclear matrix element indicates that subtle cancellations between the contributions of intermediate states can explain the different theoretical predictions. For the EFT, we also present results for two-neutrino $\beta\beta$ decays to the first excited $0^+$ state in other nuclei., Comment: 9 pages, 4 figures, 1 table; corresponds to the published version

We study two-neutrino ($2\nu\beta\beta$) and neutrinoless double-$\beta$ ($0\nu\beta\beta$) decays in the nuclear shell model and proton-neutron quasiparticle random-phase approximation (pnQRPA) frameworks. Calculating the decay half-life of several dozens of nuclei ranging from calcium to xenon with the shell model, and of $\beta\beta$ emitters with a wide range of proton-neutron pairing strengths in the pnQRPA, we observe good linear correlations between $2\nu\beta\beta$- and $0\nu\beta\beta$-decay nuclear matrix elements for both methods. We then combine the correlations with measured $2\nu\beta\beta$-decay half-lives to predict $0\nu\beta\beta$-decay matrix elements with theoretical uncertainties based on our systematic calculations. Our results include two-body currents and the short-range $0\nu\beta\beta$-decay operator., Comment: 12 pages, 8 figures, 4 tables; corresponds to the published version

Lepton-flavor-violating decays of light pseudoscalars, $P=\pi^0,\eta,\eta'\to\mu e$, are stringently suppressed in the Standard Model up to tiny contributions from neutrino oscillations, so that their observation would be a clear indication for physics beyond the Standard Model. However, in effective field theory such decays proceed via axial-vector, pseudoscalar, or gluonic operators, which are, at the same time, probed in spin-dependent $\mu\to e$ conversion in nuclei. We derive master formulae that connect both processes in a model-independent way in terms of Wilson coefficients, and study the implications of current $\mu\to e$ limits in titanium for the $P\to\mu e$ decays. We find that these indirect limits surpass direct ones by many orders of magnitude., Comment: 9 pages, 2 figures; journal version

The discovery of neutrinoless double-beta decay could soon be within reach. This hypothetical ultra-rare nuclear decay offers a privileged portal to physics beyond the Standard Model of particle physics. Its observation would constitute the discovery of a matter-creating process, corroborating leading theories of why the universe contains more matter than antimatter, and how forces unify at high energy scales. It would also prove that neutrinos and anti-neutrinos are not two distinct particles, but can transform into each other, with their mass described by a unique mechanism conceived by Majorana. The recognition that neutrinos are not massless necessitates an explanation and has boosted interest in neutrinoless double-beta decay. The field stands now at a turning point. A new round of experiments is currently being prepared for the next decade to cover an important region of parameter space. In parallel, advances in nuclear theory are laying the groundwork to connect the nuclear decay with the underlying new physics. Meanwhile, the particle theory landscape continues to find new motivations for neutrinos to be their own antiparticle. This review brings together the experimental, nuclear theory, and particle theory aspects connected to neutrinoless double-beta decay, to explore the path toward - and beyond - its discovery., Comment: review article (with illustrations by Laura Manenti)

We study neutrinoless double-beta decay in an effective field theory (EFT) for heavy nuclei, which are treated as a spherical core coupled to additional neutrons and/or protons. Since the low-energy constants of the EFT cannot be fitted to data for this unobserved decay, we follow an alternative strategy to constrain these through a correlation with double Gamow-Teller transitions. This correlation was recently found to hold for shell-model calculations, energy-density functionals, and other nuclear structure models. We therefore first calculate the nuclear matrix elements for double Gamow-Teller transitions in the EFT for heavy nuclei. The combination of the EFT uncertainty with the correlation uncertainty enables predictions of nuclear matrix elements for neutrinoless double-beta decay for a broad range of isotopes with quantified uncertainties. Generally the EFT predicts smaller nuclear matrix elements compared to other approaches, but our EFT results are consistent with recent ab initio calculations., Comment: 13 pages, 5 figures, 4 tables, minor changes, published version

Measurements of the beta-neutrino correlation coefficient (a$_{\beta\nu}$) in nuclear beta decay, together with the Fierz interference term (b$_F$), provide a robust test for the existence of exotic interactions beyond the Standard Model of Particle Physics. The extraction of these quantities from the recoil ion spectra in $\beta$-decay requires accurate knowledge, decay branching ratios, and high-precision calculations of higher order nuclear effects. Here, we report on a new measurement of the $^{23}$Ne $\beta$-decay branching ratio, which allows a reanalysis of existing high-precision measurements. Together with new theoretical calculations of nuclear structure effects, augmented with robust theoretical uncertainty, this measurement improves on the current knowledge of a$_{\beta\nu}$ in $^{23}$Ne by an order of magnitude, and strongly constrains the Fierz term in beta decays, making this one of the first extractions to constrain both terms simultaneously. Together, these results place bounds on the existence of exotic tensor interactions and pave the way for new, even higher precision, experiments.

We evaluate the leading-order short-range nuclear matrix element for the neutrinoless double-beta ($0\nu\beta\beta$) decay of the nuclei most relevant for experiments, including $^{76}$Ge, $^{100}$Mo, $^{130}$Te and $^{136}$Xe. In our calculations, performed with the nuclear shell model and proton-neutron quasiparticle random-phase approximation (pnQRPA) methods, we estimate the coupling of this term by the contact charge-independence-breaking coupling of various nuclear Hamiltonians. Our results suggest a significant impact of the short-range matrix element, which is about $15\%-50\%$ and $30\%-80\%$ of the standard $0\nu\beta\beta$-decay long-range matrix element for the shell model and pnQRPA, respectively. Combining the full matrix elements with the results from current $0\nu\beta\beta$-decay experiments we find that, if both matrix elements carry the same sign, these searches move notably toward probing the inverted mass ordering of neutrino masses., Comment: 9 pages, 3 figures, 2 tables; Fixed sign error in Eq. (5)

We quantify the extent to which future experiments will test the existence of neutrinoless double-beta decay mediated by light neutrinos with inverted-ordered masses. While it remains difficult to compare measurements performed with different isotopes, we find that future searches will fully test the inverted ordering scenario, as a global, multi-isotope endeavor. They will also test other possible mechanisms driving the decay, including a large uncharted region of the allowed parameter space assuming that neutrino masses follow the normal ordering., Comment: 6 pages, 2 figures, version accepted by Phys. Rev. C

Common and community software packages, such as ROOT, Geant4 and event generators have been a key part of the LHC's success so far and continued development and optimisation will be critical in the future. The challenges are driven by an ambitious physics programme, notably the LHC accelerator upgrade to high-luminosity, HL-LHC, and the corresponding detector upgrades of ATLAS and CMS. In this document we address the issues for software that is used in multiple experiments (usually even more widely than ATLAS and CMS) and maintained by teams of developers who are either not linked to a particular experiment or who contribute to common software within the context of their experiment activity. We also give space to general considerations for future software and projects that tackle upcoming challenges, no matter who writes it, which is an area where community convergence on best practice is extremely useful., Comment: 40 pages contribution to Snowmass 2021

The cross section for coherent elastic neutrino-nucleus scattering (CE$\nu$NS) depends on the response of the target nucleus to the external current, in the Standard Model (SM) mediated by the exchange of a $Z$ boson. This is typically subsumed into an object called the weak form factor of the nucleus. Here, we provide results for this form factor calculated using the large-scale nuclear shell model for a wide range of nuclei of relevance for current CE$\nu$NS experiments, including cesium, iodine, argon, fluorine, sodium, germanium, and xenon. In addition, we provide the responses needed to capture the axial-vector part of the cross section, which does not scale coherently with the number of neutrons, but may become relevant for the SM prediction of CE$\nu$NS on target nuclei with nonzero spin. We then generalize the formalism allowing for contributions beyond the SM. In particular, we stress that in this case, even for vector and axial-vector operators, the standard weak form factor does not apply anymore, but needs to be replaced by the appropriate combination of the underlying nuclear structure factors. We provide the corresponding expressions for vector, axial-vector, but also (pseudo-)scalar, tensor, and dipole effective operators, including two-body-current effects as predicted from chiral effective field theory. Finally, we update the spin-dependent structure factors for dark matter scattering off nuclei according to our improved treatment of the axial-vector responses., Comment: 28 pages, 11 figures; journal version

Coherent elastic neutrino-nucleus scattering (CE$\nu$NS), a process recently measured for the first time at ORNL's Spallation Neutron Source, is directly sensitive to the weak form factor of the nucleus. The European Spallation Source (ESS), presently under construction, will generate the most intense pulsed neutrino flux suitable for the detection of CE$\nu$NS. In this paper we quantify its potential to determine the root mean square radius of the point-neutron distribution, for a variety of target nuclei and a suite of detectors. To put our results in context we also derive, for the first time, a constraint on this parameter from the analysis of the energy and timing data of the CsI detector at the COHERENT experiment., Comment: 16 pages, 4 figures. Minor text changes and references added. Matches version accepted by JHEP

We review the main software and computing challenges for the Monte Carlo physics event generators used by the LHC experiments, in view of the High-Luminosity LHC (HL-LHC) physics programme. This paper has been prepared by the HEP Software Foundation (HSF) Physics Event Generator Working Group as an input to the LHCC review of HL-LHC computing, which has started in May 2020., Comment: 20 pages; editors Andrea Valassi, Efe Yazgan and Josh McFayden; addressed additional comments by journal reviewers

Chiral effective field theory allows one to calculate the response of few-nucleon systems to external currents, both for currents that can be probed in the Standard Model and ones that only exist in Standard-Model extensions. In combination with state-of-the-art many-body methods, the constraints from chiral symmetry can then be implemented in nuclear structure factors that describe the response of atomic nuclei in direct-detection searches for dark matter. We review the present status of this approach, including the role of coherently enhanced two-body currents, the discrimination of dark matter candidates based on the nuclear response functions, and limits on Higgs-portal dark matter., Comment: 10 pages, 6 figures, proceedings for The 9th International Workshop on Chiral Dynamics, Durham, North Carolina, USA, September 17-21, 2018

We present nuclear structure factors that describe the generalized spin-independent coupling of weakly interacting massive particles (WIMPs) to nuclei. Our results are based on state-of-the-art nuclear structure calculations using the large-scale nuclear shell model. Starting from quark- and gluon-level operators, we consider all possible coherently enhanced couplings of spin-1/2 and spin-0 WIMPs to one and two nucleons up to third order in chiral effective field theory. This includes a comprehensive discussion of the structure factors corresponding to the leading two-nucleon currents covering, for the first time, the contribution of spin-2 operators. We provide results for the most relevant nuclear targets considered in present and planned dark matter direct detection experiments: fluorine, silicon, argon, and germanium, complementing our previous work on xenon. All results are also publicly available in a Python notebook., Comment: 24 pages, 19 figures, Python notebook available at https://theorie.ikp.physik.tu-darmstadt.de/strongint/ChiralEFT4DM.html; further details on nuclear structure added, journal version

The calculated nuclear matrix elements for the neutrinoless double-beta ($0\nu\beta\beta$) decay suffer from several limitations. Predicted matrix-element values depend on the many-body method used to calculate them and, in addition, they may need to be "quenched", albeit by an unknown amount due to the relatively high momentum transferred in the decay. These uncertainties are hard to figure out from theoretical calculations alone because of the unique character of $0\nu\beta\beta$ decay, not clearly related to other nuclear structure observables. We present recent progress in the determination of the nuclear matrix elements. First, we report improved shell model calculations in an extended configuration space of two harmonic oscillator shells. Second, we note the relation between the $0\nu\beta\beta$ decay and double Gamow-Teller transitions that can in principle be measured in double charge-exchange reactions. These new insights pave the way for a more reliable estimation of the value of the nuclear matrix elements., Comment: 6 pages, 3 figures, presented at the IIRC Symposium "Perspectives of the Physics of Nuclear Structure", Tokyo, Japan, November 1-4, 2017

Neutrinoless $\beta\beta$ decay nuclear matrix elements calculated with the shell model and energy-density functional theory typically disagree by more than a factor of two in the standard scenario of light-neutrino exchange. In contrast, for a decay mediated by sterile heavy neutrinos the deviations are reduced to about 50\%, an uncertainty similar to the one due to short-range effects. We compare matrix elements in the light- and heavy-neutrino-exchange channels, exploring the radial, momentum transfer and angular momentum-parity matrix element distributions, and considering transitions that involve correlated and uncorrelated nuclear states. We argue that the shorter-range heavy-neutrino exchange is less sensitive to collective nuclear correlations, and that discrepancies in matrix elements are mostly due to the treatment of long-range correlations in many-body calculations. Our analysis supports previous studies suggesting that isoscalar pairing correlations, which affect mostly the longer-range part of the neutrinoless $\beta\beta$ decay operator, are partially responsible for the differences between nuclear matrix elements in the standard light-neutrino-exchange mechanism., Comment: 20 pages, 8 figures, published version

Searches of neutrinoless double-beta decay require information on the value of the nuclear matrix elements that rule the process to plan and interpret experiments. At present, however, even the matrix elements obtained with the most reliable many-body approaches do not agree to each other better than a factor two or three. A usual test of the many-body calculations is the comparison to several nuclear observables, but so far no nuclear structure property has been found to show a good correlation to neutrinoless double-beta decay. Here we propose that double charge-exchange experiments can offer a very valuable tool to provide insights on neutrinoless double-beta decay. Double charge-exchange reactions are being currently performed in various laboratories worldwide and aim to find the novel nuclear collectivity given by double Gamow-Teller excitations. Our results suggest a good linear correlation between double Gamow-Teller transitions to the ground state of the final nucleus and neutrinoless double-beta decay nuclear matrix elements. The correlation seems robust across $pf$-shell nuclei., Comment: 6 pages, 2 figures. Proceedings of the Conference on Neutrinos and Nuclear Physics (CNNP2017), Catania, 15-21 October, 2017

We study the double Gamow-Teller (DGT) strength distribution of $^{48}$Ca with state-of-the-art large-scale nuclear shell model calculations. Our analysis shows that the centroid energy of the DGT giant resonance depends mostly on the isovector pairing interaction, while the resonance width is more sensitive to isoscalar pairing. Pairing correlations are also key in neutrinoless $\beta\beta$ ($0\nu\beta\beta$) decay. We find a simple relation between the centroid energy and width of the $^{48}$Ca DGT giant resonance and %the value of the $0\nu\beta\beta$ decay nuclear matrix element. More generally, we observe a very good linear correlation between the DGT transition to the ground state of the final nucleus and the $0\nu\beta\beta$ decay matrix element. The correlation, which originates on the dominant short-range character of both transitions, extends to heavier systems including several $\beta\beta$ emitters, and also holds in energy-density functional results. Our findings suggest that DGT experiments can be a very valuable tool to obtain information on the value of $0\nu\beta\beta$ decay nuclear matrix elements., Comment: 6 pages, 5 figures, published version

We study $\beta$ decays within an effective theory that treats nuclei as a spherical collective core with an even number of neutrons and protons that can couple to an additional neutron and/or proton. First we explore Gamow-Teller $\beta$ decays of parent odd-odd nuclei into low-lying ground-, one-, and two-phonon states of the daughter even-even system. The low-energy constants of the effective theory are adjusted to data on $\beta$ decays to ground states or Gamow-Teller strengths. The corresponding theoretical uncertainty is estimated based on the power counting of the effective theory. For a variety of medium-mass and heavy isotopes the theoretical matrix elements are in good agreement with experiment within the theoretical uncertainties. We then study the two-neutrino double-$\beta$ decay into ground and excited states. The results are remarkably consistent with experiment within theoretical uncertainties, without the necessity to adjust any low-energy constants., Comment: 17 pages, 5 figures, results extended to two-neutrino double beta-minus decays and two-neutrino double electron-capture decays to excited 2+ states, matches published version

Searches for invisible Higgs decays at the Large Hadron Collider constrain dark matter Higgs-portal models, where dark matter interacts with the Standard Model fields via the Higgs boson. While these searches complement dark matter direct-detection experiments, a comparison of the two limits depends on the coupling of the Higgs boson to the nucleons forming the direct-detection nuclear target, typically parameterized in a single quantity $f_N$. We evaluate $f_N$ using recent phenomenological and lattice-QCD calculations, and include for the first time the coupling of the Higgs boson to two nucleons via pion-exchange currents. We observe a partial cancellation for Higgs-portal models that makes the two-nucleon contribution anomalously small. Our results, summarized as $f_N=0.308(18)$, show that the uncertainty of the Higgs-nucleon coupling has been vastly overestimated in the past. The improved limits highlight that state-of-the-art nuclear physics input is key to fully exploiting experimental searches., Comment: 6 pages, 3 figures; version to appear in PRL

Experiments using nuclei to probe new physics beyond the Standard Model, such as neutrinoless $\beta\beta$ decay searches testing whether neutrinos are their own antiparticle, and direct detection experiments aiming to identify the nature of dark matter, require accurate nuclear physics input for optimizing their discovery potential and for a correct interpretation of their results. This demands a detailed knowledge of the nuclear structure relevant for these processes. For instance, neutrinoless $\beta\beta$ decay nuclear matrix elements are very sensitive to the nuclear correlations in the initial and final nuclei, and the spin-dependent nuclear structure factors of dark matter scattering depend on the subtle distribution of the nuclear spin among all nucleons. In addition, nucleons are composite and strongly interacting, which implies that many-nucleon processes are necessary for a correct description of nuclei and their interactions. It is thus crucial that theoretical studies and experimental analyses consider $\beta$ decays and dark matter interactions with a coupling to two nucleons, called two-nucleon currents., Comment: 11 pages, 5 figures, invited parallel talk at the XIIth Quark Confinement & the Hadron Spectrum conference, Thessaloniki, Greece, 2016

The nuclear matrix elements that govern the rate of neutrinoless double beta decay must be accurately calculated if experiments are to reach their full potential. Theorists have been working on the problem for a long time but have recently stepped up their efforts as ton-scale experiments have begun to look feasible. Here we review past and recent work on the matrix elements in a wide variety of nuclear models and discuss work that will be done in the near future. Ab initio nuclear-structure theory, which is developing rapidly, holds out hope of more accurate matrix elements with quantifiable error bars., Comment: 35 pages, 14 figures

We propose a formalism for the analysis of direct-detection dark-matter searches that covers all coherent responses for scalar and vector interactions and incorporates QCD constraints imposed by chiral symmetry, including all one- and two-body WIMP-nucleon interactions up to third order in chiral effective field theory. One of the free parameters in the WIMP-nucleus cross section corresponds to standard spin-independent searches, but in general different combinations of new-physics couplings are probed. We identify the interference with the isovector counterpart of the standard spin-independent response and two-body currents as the dominant corrections to the leading spin-independent structure factor, and discuss the general consequences for the interpretation of direct-detection experiments, including minimal extensions of the standard spin-independent analysis. Fits for all structure factors required for the scattering off xenon targets are provided based on state-of-the-art nuclear shell-model calculations., Comment: 20 pages, 7 figures; v2 includes the two-body current related to the trace anomaly of the QCD energy-momentum tensor