$^{16}$O spectral function from coupled-cluster theory: applications to lepton-nucleus scattering

Joanna E. Sobczyk,Sonia Bacca

2023-09-01

Abstract:We calculate the $^{16}$O spectral function by combining coupled-cluster theory with a Gaussian integral transform and by expanding the integral kernel in terms of Chebyshev polynomials to allow for a quantification of the theoretical uncertainties. We perform an analysis of the spectral function and employ it to predict lepton-nucleus scattering. Our results well describe the $^{16}$O electron scattering data in the quasi-elastic peak for momentum transfers $|\mathbf{q}|\gtrapprox500$ MeV and electron energies up to 1.2 GeV, extending therefore the so-called first principles approach to lepton-nucleus cross sections well into the relativistic regime. To prove the applicability of this method to neutrino-nucleus cross sections, we implement our $^{16}$O spectral functions in the NuWro Monte Carlo event generator and provide a comparison with recently published T2K neutrino data.

Nuclear Theory,High Energy Physics - Phenomenology

What problem does this paper attempt to address?

The problems that this paper attempts to solve mainly focus on the following aspects:

1. **Accurately describe the spectral function of \(^{16}\text{O}\)**: By combining the coupled - cluster theory with Gaussian integral transformation and using Chebyshev polynomial to expand the integral kernel, the authors calculated the spectral function of \(^{16}\text{O}\). This method allows the quantification of theoretical uncertainties.

2. **Predict lepton - nucleus scattering**: Based on the calculated spectral function of \(^{16}\text{O}\), the authors predicted the lepton - nucleus scattering process. Their results well describe the \(^{16}\text{O}\) electron scattering data with momentum transfer \( |q| \geq 500 \, \text{MeV}\) and electron energy up to 1.2 GeV, thus extending the so - called "first - principles method" to the relativistic region.

3. **Apply to neutrino - nucleus scattering**: To prove that this method is suitable for calculating the neutrino - nucleus scattering cross - section, the authors implemented the spectral function of \(^{16}\text{O}\) into the NuWro Monte Carlo event generator and compared it with the recently released T2K neutrino data.

### Main contributions of the paper

- **Theoretical framework**: Proposed a method that combines the coupled - cluster theory and Chebyshev polynomial expansion to calculate the spectral function of \(^{16}\text{O}\) and quantified the theoretical uncertainties.
- **Experimental verification**: Verified the effectiveness of this method by comparing with the electron scattering data of \(^{16}\text{O}\), especially in the quasi - elastic peak region.
- **Neutrino scattering application**: Applied the calculated spectral function to neutrino - nucleus scattering, providing theoretical support for future neutrino experiments (such as T2K and T2HK).

### Formula explanation

- **Definition of spectral function**:
\[
S(p, E) = \sum_{\alpha, \beta} \langle p | \alpha \rangle \langle p | \beta \rangle^\dagger \int d^3 x \, \Phi_{A - 1} \langle 0 | a^\dagger_\beta | \Phi_{A - 1} \rangle \langle \Phi_{A - 1} | a_\alpha | 0 \rangle \delta(E-(E_0 - E_\Phi))
\]

- **Hadronic tensor**:
\[
W^{\mu\nu}(q) = \int \frac{d^3 p}{(2\pi)^3} \frac{dE}{E_p m} \left[ S_n(p, E) w^{\mu\nu}_n(p, q)+S_p(p, E) w^{\mu\nu}_p(p, q) \right] \times \delta(\omega + E - E_{p + q}-E_{\text{kin}}^f)
\]

- **Integral transformation**:
\[
\text{Im}\tilde{G}_h(\alpha, \beta, E) = \int d\omega \, \text{Im}G_h(\alpha, \beta, \omega) K(\omega, E)=-\pi \langle 0 | a^\dagger_\beta | K(H, E - E_0) | a_\alpha | 0 \rangle
\]

- **Chebyshev polynomial expansion**:
\[
K(\omega, E) = \sum_{k = 0}^N c_k(E) T_k(\omega)
\]

- **Uncertainty estimation of spectral function**:
\[
\text{Im}\tilde{G}_h(\Delta - \Lambda)-\Sigma - 2\gamma(\Delta - \Lambda) \leq \text{Im}G_h(\Delta) \leq \text{Im}\tilde{G}