Nuclear Response Functions Research Articles

Dark matter scattering off H2 and He4 nuclei within chiral effective field theory

We study dark matter, assumed to be composed of weak interacting massive particles (WIMPs), scattering off H2 and He4 nuclei. In order to parametrize the WIMP-nucleon interaction, the chiral effective field theory approach is used. Considering only interactions invariant under parity, charge conjugation, and time reversal, we examine five interaction types: scalar, pseudoscalar, vector, axial, and tensor. Scattering amplitudes between two nucleons and a WIMP are determined up to second order of chiral perturbation theory. We apply this program to calculate the interaction rate as a function of the WIMP mass and of the magnitude of the WIMP-quark coupling constants. From our study, we conclude that the scalar nuclear response functions is much greater than the others due to their large combination of low energy constants. We verify that the leading order contributions are dominant in these low energy processes. We also provide an estimate for the background due to atmospheric neutrinos. Published by the American Physical Society 2024

Earth-catalyzed detection of magnetic inelastic dark matter with photons in large underground detectors

Inelastic dark matter with moderate splittings, O\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{O} $$\\end{document}(few to 150) keV, can upscatter to an excited state in the Earth, with the excited state subsequently decaying, leaving a distinctive monoenergetic photon signal in large underground detectors. The photon signal can exhibit sidereal-daily modulation, providing excellent separation from backgrounds. Using a detailed numerical simulation, we examine this process as a search strategy for magnetic inelastic dark matter with the dark matter mass near the weak scale, where the upscatter to the excited state and decay proceed through the same magnetic dipole transition operator. At lower inelastic splittings, the scattering is dominated by moderate mass elements in the Earth with high spin, especially 27Al, while at larger splittings, 56Fe becomes the dominant target. We show that the proposed large volume gaseous detector CYGNUS will have excellent sensitivity to this signal. Xenon detectors also provide excellent sensitivity through the inelastic nuclear recoil signal, and if a future signal is seen, we show that the synergy among both types of detection can provide strong evidence for magnetic inelastic dark matter. In the course we have calculated nuclear response functions for elements relevant for scattering in the Earth, which are publicly available on GitHub.

Impact of shell model interactions on nuclear responses to WIMP elastic scattering

Nuclear recoil from scattering with weakly interacting massive particles (WIMPs) is a signature searched for in direct detection of dark matter. The underlying WIMP-nucleon interactions could be spin and/or orbital angular momentum (in)dependent. Evaluation of nuclear recoil rates through these interactions requires accounting for nuclear structure, e.g., through shell model calculations. We evaluate nuclear response functions induced by these interactions for F19, Na23, Si28,29,30, Ar40, Ge70,72,73,74,76, I127, and Xe128,129,130,131,132,134,136 nuclei that are relevant to current direct detection experiments, and estimate their sensitivity to shell model interactions. Shell model calculations are performed with the NuShellX solver. Nuclear response functions from nonrelativistic effective field theory are evaluated and integrated over transferred momentum for quantitative comparisons. We show that although the standard spin-independent response is barely sensitive to the structure of the nuclei, large variations with the shell model interaction are often observed for the other channels. Significant uncertainties may arise from the nuclear components of WIMP-nucleus scattering amplitudes due to nuclear structure theory and modeling. These uncertainties should be accounted for in analyses of direct detection experiments. Published by the American Physical Society 2024

The Interrelated Roles of Correlations in the Nuclear Equation of State and in Response Functions: Application to a Chiral Confining Theory

We study the role of short-range correlations, as well as pion and rho loops governing long-range RPA correlations, in nuclear matter properties and response functions. We use an adapted formulation of the Brueckner G-matrix approach to generate a pair correlation function satisfying the Beg–Agassi–Gal theorem, providing a natural cutoff to the loop integrals. We present results for the case of a relativistic chiral theory, including the effects of quark confinement and of the chirally broken vacuum in a version where parameters are directly connected to QCD observables or constrained by well-established hadron phenomenology. This provides a unified and coherent view of the nuclear matter equation of state and the effect of correlations on neutrino–nucleus scattering.

A center of mass free equation of motion method and its application to 4He

We report on an upgraded version of the equation of motion phonon method (EMPM) which exploits the singular value decomposition (SVD) of rectangular matrices to decouple completely and exactly the intrinsic from the center of mass (c.m.) translational motion under no restrictions and for an arbitrary single particle (s.p.) basis. Its numerical implementation on 4He illustrates the way such a decoupling is achieved and its important effects on spectra and nuclear response functions.

Dmscatter: A fast program for WIMP-nucleus scattering

Recent work [1,2], using an effective field theory framework, have shown the number of possible couplings between nucleons and the dark-matter-candidate Weakly Interacting Massive Particles (WIMPs) is larger than previously thought. Inspired by an existing Mathematica script that computes the target response [2], we have developed a fast, modern Fortran code, including optional OpenMP parallelization, along with a user-friendly Python wrapper, to swiftly and efficiently explore many scenarios, with output aligned with practices of current dark matter searches. A library of most of the important target nuclides is included; users may also import their own nuclear structure data, in the form of reduced one-body density matrices. The main output is the differential event rate as a function of recoil energy, needed for modeling detector response rates, but intermediate results such as nuclear form factors can be readily accessed. Program summaryProgram Title: dmscatterCPC Library link to program files:https://doi.org/10.17632/6fm7d4m2py.1Developer's repository link:https://github.com/ogorton/dmscatterLicensing provisions: MITProgramming languages: Fortran, PythonNature of problem: Simulating the event rate of nuclear recoils from collisions with dark matter for a variety of different nuclear targets and different nucleon-dark matter couplings.Solution method: The event rate is an integral over the product of the nuclear and dark matter response functions, weighted by the expected dark matter flux. To compute the nuclear response, reduced one-body nuclear density matrices are coupled to multipole expansion of operators computed analytically in a harmonic oscillator basis. The nuclear structure input, in the form of one-body densities matrices, are supplied by outside code; we provide a library of the most prominent targets.

Effective theory of nuclear scattering for a WIMP of arbitrary spin

We introduce a systematic approach to characterize the most general non-relativistic WIMP-nucleus interaction allowed by Galilean invariance for a WIMP of arbitrary spin $j_\chi$ in the approximation of one-nucleon currents. Five nucleon currents arise from the nonrelativistic limit of the free nucleon Dirac bilinears. Our procedure consists in (1) organizing the WIMP currents according to the rank of the $2 j_\chi+1$ irreducible operator products of up to $2 j_\chi$ WIMP spin vectors, and (2) coupling each of the WIMP currents to each of the five nucleon currents. The transferred momentum $q$ appears to a power fixed by rotational invariance. For a WIMP of spin $j_\chi$ we find a basis of 4+20$j_\chi$ independent operators that exhaust all the possible operators that drive elastic WIMP-nucleus scattering in the approximation of one-nucleon currents. By comparing our operator basis, which is complete, to the operators already introduced in the literature we show that some of the latter for $j_\chi=1$ were not independent and some were missing. We provide explicit formulas for the squared scattering amplitudes in terms of the nuclear response functions, which are available in the literature for most of the targets used in WIMP direct detection experiments.

Machine-learning-based inversion of nuclear responses

A microscopic description of the interaction of atomic nuclei with external electroweak probes is required for elucidating aspects of short-range nuclear dynamics and for the correct interpretation of neutrino oscillation experiments. Nuclear quantum Monte Carlo methods infer the nuclear electroweak response functions from their Laplace transforms. Inverting the Laplace transform is a notoriously ill-posed problem; and Bayesian techniques, such as maximum entropy, are typically used to reconstruct the original response functions in the quasielastic region. In this work, we present a physics-informed artificial neural network architecture suitable for approximating the inverse of the Laplace transform. Utilizing simulated, albeit realistic, electromagnetic response functions, we show that this physics-informed artificial neural network outperforms maximum entropy in both the low-energy transfer and the quasielastic regions, thereby allowing for robust calculations of electron scattering and neutrino scattering on nuclei and inclusive muon capture rates.

Ab Initio Study of (νℓ,ℓ−) and (ν¯ℓ,ℓ+) Inclusive Scattering in C12 : Confronting the MiniBooNE and T2K CCQE Data

We carry out an ab initio calculation of the neutrino flux-folded inclusive cross sections, measured on $^{12}$C by the MiniBooNE and T2K collaborations in the charged-current quasielastic (CCQE) regime. The calculation is based on realistic two- and three-nucleon interactions, and on a realistic nuclear electroweak current with one-and two-nucleon terms that are constructed consistently with these interactions and reproduce low-energy electroweak transitions. Numerically exact quantum Monte Carlo methods are utilized to compute the nuclear weak response functions, by fully retaining many-body correlations in the initial and final states and interference effects between one- and two-body current contributions. We employ a nucleon axial form factor of the dipole form with $\Lambda_A = 1.0$ or $1.15$ GeV, the latter more in line with a very recent lattice QCD determination. The calculated cross sections are found to be in good agreement with the neutrino data of MiniBooNE and T2K, and antineutrino MiniBooNE data, yielding a consistent picture of nuclei and their electroweak properties across a wide regime of energy and momenta.

Response functions and cross sections for inclusive neutrino scattering off $^2$H, $^3$H and $^3$He

For several decades there has been an interest in studying neutrino scattering on light nuclei, with the hope of increasing knowledge about the theory of weak interactions and electroweak unification. A study of several neutrino induced reactions on light nuclei (^22H, ^33H and ^33He) is presented here. The cross section for these reactions is evaluated using the nuclear response functions, which holds the information about the nuclear interactions and structure. This approach has been tested for break-up reactions on the deuteron before and is now extended to three-nucleon break-up reactions.

Nuclear response in a finite-temperature relativistic framework

A thermal extension of the relativistic nuclear field theory is formulated for the nuclear response. The Bethe-Salpeter equation (BSE) with the time-dependent kernel for the particle-hole response is treated within the Matsubara Green's function formalism. We show that, with the help of a temperature-dependent projection operator on the subspace of the imaginary time (time blocking), it is possible to reduce the BSE for the nuclear response function to a single frequency variable equation also at finite temperature. The approach is implemented self-consistently in the framework of quantum hadrodynamics based on the meson-nucleon Lagrangian. The method is applied to the monopole, dipole and quadrupole response of $^{48}$Ca and to the dipole response of the tin isotopes $^{100,120,132}$Sn, in particular, to a study of the evolution of nuclear collective oscillations with temperature. The article is dedicated to the memory of Pier Francesco Bortignon and devoted to the developments related to his pioneering ideas.

Nuclear dipole response in the finite-temperature relativistic time-blocking approximation

The radiative neutron capture reaction rates of the r-process nucleosynthesis are immensely affected by the microscopic structure of the low-energy spectra of compound nuclei. The relativistic (quasiparticle) time blocking approximation (R(Q)TBA) has successfully provided a good description of the low-energy strength, in particular, the strength associated with pygmy dipole resonance, describing transitions from and to the nuclear ground state. The finite-temperature generalization of this method is designed for thermally excited compound nuclei and has the potential to enrich the fine structure of the dipole strength, especially in the low-energy region. The finite-temperature RTBA equations are derived using the Matsubara Green's function formalism. We show that with the help of a temperature-dependent projection operator on the subspace of the imaginary time it is possible to reduce the Bethe-Salpeter equation for the nuclear response function to a single frequency variable equation also at finite temperatures. The approach is implemented self-consistently in the framework of quantum hadrodynamics and keeps the ability of connecting the high-energy scale of heavy mesons and the low-energy domain of nuclear medium polarization effects in a parameter-free way. The presented calculations of the dipole response within a self-consistent relativistic framework reveal that, although the Landau damping plays the leading role in the temperature evolution of the strength distribution, (i) at moderate temperatures the PVC effects remain almost as strong as at $T=0$ and (ii) at high temperatures they are tremendously reinforced because of the formation of the new collective low-energy modes. In the dipole channel, the latter effect is responsible for the "disappearance" of the high-frequency GDR or, in other words, brings the GDR to the low-energy domain.

Nuclear response functions with low-momentum interactions

Linear density response functions for nuclear matter are calculated with separable 1S0 and 3S1 2-nucleon interactions obtained from experimental phase-shifts by inverse scattering techniques. It has been shown that results of nuclear matter binding energy Brueckner calculations with these potentials agree closely with those using realistic Bonn interactions. It has also been shown that low-momentum versions give binding energies (nearly) independent of momentum cut-offs $\Lambda \geq \sim 2$ fm-1. Shown here is that the $\Lambda = 2$ rank-one version of the separable interaction (for the S-states) yields a near agreement between a second-order and an all-order (Brueckner) binding energy calculation. Second order calculations of response functions are then made with this version using a method that Kwong and Bonitz used for the Coulomb gas. It is based on time-evolving two-time Green’s functions by Kadanoff-Baym (KB) equations. The effect of correlations, going beyond the conventional HF+RPA method, is included by “dressing” the Green’s function propagators with time-dependent complex self-energies in addition to the real HF-field. It would be of interest to see if the calculated response is also, like the binding energy, independent of cut-offs larger than $\sim 2$ fm-1. That would require using separable potentials of a higher rank, not accommodated by the present computer program. Some results with a 1.5 fm-1 cut-off are however also shown below together with those with a 2.0 fm-1 cut-off. The KB-equations are time-evolved until the nuclear medium is fully correlated. The time-dependent density response to an external perturbation is then registered and Fourier transformed to obtain energy-dependent density response functions. Most previous calculations have been made with in-medium effective (e.g., Gogny or Skyrme) interactions. The medium dependence is in the present calculations contained in the second-order terms.

Molecular dynamics simulations for optical Kerr effect of TIP4P/2005 water in liquid and supercooled states

The optical Kerr effect (OKE) spectroscopy measured with heterodyne detection (HD) is a useful tool to provide information regarding intermolecular vibrations and structural relaxations in liquid water. Recently, the measurements of the OKE spectroscopy have been extended to the supercooled regime of water. Though the measured results can be well described by using a phenomenological model, the time-resolved OKE spectroscopy of liquid and supercooled water still need a comprehensive understanding. In this paper, we investigated the OKE nuclear response functions of this peculiar liquid and their reduced spectral densities by performing molecular dynamics simulations with the TIP4P/2005 water model. The collective polarizability of water was computed via a dipolar induction scheme, which involves the intrinsic polarizability and the first-order hyperpolarizability tensor of water molecule. Our simulation results were qualitatively consistent with the HD-OKE experimental observations for displaying that the polarizability anisotropy relaxation of supercooled water in the high-density liquid phase was fragile-like by following a stretching exponential decay with an exponent βs insensitive to temperature and the temperature dependence of the relaxation time exhibited a power-law divergence at a singular temperature Ts with a critical exponent γs. Indicated by our quantitative results, Ts was predominately determined by the structural arrest, but βs and γs were not only related to the structural relaxation but also influenced by the collective polarizability of the liquid. For all investigations, the effects due to the first-order hyperpolarizability tensor were examined.

A combined analysis of PandaX, LUX, and XENON1T experiments within the framework of dark matter effective theory

Weakly interacting massive particles are a widely well-probed dark matter candidate by the dark matter direct detection experiments. Theoretically, there are a large number of ultraviolet completed models that consist of a weakly interacting massive particle dark matter. The variety of models makes the comparison with the direct detection data complicated and often non-trivial. To overcome this, in the non-relativistic limit, the effective theory was developed in the literature which works very well to significantly reduce the complexity of dark matter-nucleon interactions and to better study the nuclear response functions. In the effective theory framework for a spin-1/2 dark matter, we combine three independent likelihood functions from the latest PandaX, LUX, and XENON1T data, and give a joint limit on each effective coupling. The astrophysical uncertainties of the dark matter distribution are also included in the likelihood. We further discuss the isospin violating cases of the interactions. Finally, for both dimension-five and dimension-six effective theories above the electroweak scale, we give updated limits of the new physics mass scales.

In medium dispersion relation effects in nuclear inclusive reactions at intermediate and low energies

In a well-established many-body framework, successful in modeling a great variety of nuclear processes, we analyze the role of the spectral functions (SFs) accounting for the modifications of the dispersion relation of nucleons embedded in a nuclear medium. We concentrate in processes mostly governed by one-body mechanisms, and study possible approximations to evaluate the particle–hole propagator using SFs. We also investigate how to include together SFs and long-range RPA-correlation corrections in the evaluation of nuclear response functions, discussing the existing interplay between both type of nuclear effects. At low energy transfers (≤50MeV), we compare our predictions for inclusive muon and radiative pion captures in nuclei, and charge–current (CC) neutrino–nucleus cross sections with experimental results. We also present an analysis of intermediate energy quasi-elastic neutrino scattering for various targets and both neutrino and antineutrino CC driven processes. In all cases, we pay special attention to estimate the uncertainties affecting the theoretical predictions. In particular, we show that errors on the σμ∕σe ratio are much smaller than 5%, and also much smaller than the size of the SF+RPA nuclear corrections, which produce significant effects, not only in the individual cross sections, but also in their ratio for neutrino energies below 400 MeV. These latter nuclear corrections, beyond Pauli blocking, turn out to be thus essential to achieve a correct theoretical understanding of this ratio of cross sections of interest for appearance neutrino oscillation experiments. We also briefly compare our SF and RPA results to predictions obtained within other representative approaches.

Density dependence of 2p-2h meson-exchange currents

We analyze the density dependence of the contribution of meson-exchange currents to the lepton-nucleus inclusive cross section in the two-particle two-hole channel. The model is based on the Relativistic Fermi Gas, where each nucleus is characterized by its Fermi momentum $k_F$. We find that the 2p-2h nuclear response functions at their peaks scale as $A k_F^2$ for Fermi momentum going from 200 to 300 MeV/c and momentum transfer $q$ from $2k_F$ to 2 GeV/c. This behavior is different from what is found for the quasielastic response, which scales as $A/k_F$. Additionally, the deep scaling region is also discussed and there the usual scaling behavior is found to be preferable.

Ab initio nuclear response functions for dark matter searches

We study the process of dark matter particles scattering off $^{3,4}$He with nuclear wave functions computed using an ab initio many-body framework. We employ realistic nuclear interactions from chiral effective field theory at next-to-next-to-leading order (NNLO) and develop an ab initio scheme to compute a general set of different nuclear response functions. In particular, we then perform an accompanying uncertainty quantification on these quantities and study error propagation to physical observables. We find a rich structure of allowed nuclear responses with significant uncertainties for certain spin-dependent interactions. The approach and results that are presented in this Paper establish a new framework for nuclear structure calculations and uncertainty quantification in the context of direct and (certain) indirect searches for dark matter.

Chiral effective theory of dark matter direct detection

We present the effective field theory for dark matter interactions with the visible sector that is valid at scales of \U0001d4aa(1 GeV). Starting with an effective theory describing the interactions of fermionic and scalar dark matter with quarks, gluons and photons via higher dimension operators that would arise from dimension-five and dimension-six operators above electroweak scale, we perform a nonperturbative matching onto a heavy baryon chiral perturbation theory that describes dark matter interactions with light mesons and nucleons. This is then used to obtain the coefficients of the nuclear response functions using a chiral effective theory description of nuclear forces. Our results consistently keep the leading contributions in chiral counting for each of the initial Wilson coefficients.

Improvements and testing practical expressions for photon strength functions of E1 gamma-transitions

Analytical expression for the E1 photon strength functions (PSF) is modified to account for the low-energy enhancement due to nuclear structure effects (presence of low-energy state (LES)). A closed-form expression of the E1 PSF function includes response of two nuclear states: LES and giant diplole resonance (GDR). Expression for the nuclear response function on electromagnetic field is based on a model of excitation of two coupled damped states. These approach is tested for different data sets for spherical nuclei. Impact on PSF shape of coupling between LES and GDR excitations is considered.