Neutrino reactions in the resonance region Toru Sato RCNP,Osaka - - PowerPoint PPT Presentation

▶
neutrino reactions in the resonance region
SMART_READER_LITE
LIVE PREVIEW

Neutrino reactions in the resonance region Toru Sato RCNP,Osaka - - PowerPoint PPT Presentation

Dec 07, 2023 144 likes •443 views

Neutrino reactions in the resonance region Toru Sato RCNP,Osaka University J-PARC Branch, KEK Theory Center,KEK T. Sato (Osaka U.) N ( , lM ) N Nov. 12 2018 1 / 29 Motivation CP violation in lepton sector, Mass hierarchy ...

slide-1
SLIDE 1

Neutrino reactions in the resonance region

Toru Sato

RCNP,Osaka University J-PARC Branch, KEK Theory Center,KEK

  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

1 / 29

slide-2
SLIDE 2

Motivation

✞ ✝ ☎ ✆ CP violation in lepton sector, Mass hierarchy ...

CP

δ

3 − 2 − 1 − 1 2 3

ln(L) ∆

  • 2

5 10 15 20 25 30 Normal Inverted T2K Run1-8

T2K Collaboration PRL121(2018) 171802 2σ confidence interval for δCP does not contain δCP = 0, π

  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

2 / 29

slide-3
SLIDE 3

Atmospheric Neutrino:Matter effects depends on MH

5 Eν [GeV^GeV] Eν [GeV^GeV] c

  • s

( z e n i t h ) cos(zenith)

P(νμ→ νe) Vacuum P(νμ→ νe) Matter

Mass hierarchy with atmospheric neutrinos

➢ Order of neutrino mass eigenstates is not fully known ➢ Propagation in matter modifies oscillation probabilities compared to

vacuum, in different ways depending on MH

➢ In particular resonance in muon to electron flavor oscillation

NH: ν only - IH: ν only

  • C. Bronner(Workshop on Shallow and DIS Scattering 2018)
  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

3 / 29

slide-4
SLIDE 4

Neutrino flux of current and future LBL experiments (T. Katori, M. Martini arXiv 1611.0770)

(GeV)

ν

E

1 2 3 4 5 6 7 8

Arbitrary

T2K MiniBooNE/SciBooNE MINOS/MINERvA (LE)

(GeV)

ν

E

1 2 3 4 5 6 7 8

Arbitrary

T2K/Hyper-K MicroBooNE/SBND MINERvA (ME) NOvA DUNE

Neutrino energy of LBL and Atmospheric neutrino experiments ∼ GeV Importance of neutrino reaction in N∗, ∆ resonance region in neutrino detection through ν-nucleus reaction.

  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

4 / 29

slide-5
SLIDE 5

CC Neutrino-nucleon reaction (as building block to describe neutrino-nucleus reaction) ν + N → l + N l + N + π l + N + η l + Y + K l + N + π + π .... l + Y l + N + K .... νµn → µ−X (ANL-Osaka model)

1 2 3 1.2 1.4 1.6 1.8 2 dσ/dW[10-38 cm2/GeV] W[GeV] Eν=1GeV 2GeV 3GeV

  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

5 / 29

slide-6
SLIDE 6

contents Model of neutrino reaction in resonance region Pion production in ∆(1232) region Pion production in N ∗, ∆ region Model of axial vector current Summary

  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

6 / 29

slide-7
SLIDE 7

Models of Neutrino reaction in resonance region

model dependence of single pion production

Neutrino Generator RES+DIS bg, RES=Rein-Sehgal model

Alvarez-Ruso, L. and others, NuSTEC White Paper

  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

7 / 29

slide-8
SLIDE 8

∆(1232) region Beyond ∆(1232) region(W < 2GeV ) Resonane ∆ only several overlaping resonances Non-resonant smaller than Res, Chiral L comparable to Res Channels πN only πN and ππN are comparable ηN, KΛ, KΣ also

(MeV)

Total cross section of γp

  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

8 / 29

slide-9
SLIDE 9

ANL-Osaka DCC model

Model developed for N∗ physics: spectrum of nucleon excited states,transition form factors Fock-Space:isobar(N∗, ∆) , Meson-Baryon (πN, ηN, KΛ, KΣ, ππN(π∆, ρN, σN)) Interaction:isobar excitation and non-resonant meson-baryon interaction Coupled-channel(Lippmann-Schwinger)equation is solved numerically. T = V + V G0T The model is constructed by fitting available data on pion, photon, electron induced meson production reaction(two-body final state). (Recent model: H. Kamano,S.X. Nakamura, T.-S. H. Lee, TS, PRC88,035209(2013) the model is extended for neutron and axial vector current. Axial vector current: gNN∗

A

from gNN∗

π

assuming PCAC and dipole form factor. Neutron: H. Kamano,S.X. Nakamura,T.-S. H. Lee,TS, PRC94 015291 (2016) Neutrino:S. X. Nakamura,H. Kamano, TS,PRD92 07402(2015)

  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

9 / 29

slide-10
SLIDE 10

π0, π+ electroproduction on proton (σT + ϵσL for W = 1.1 − 1.68GeV at Q2 = 0.4(GeV/c)2 electromagnetic NN∗ transition form factors are extracted by analyzing (e, e′π).

  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

10 / 29

slide-11
SLIDE 11

Neutrino induced pion production in ∆(1232) region

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0.5 1 1.5 2 2.5 3 3.5 4 𝜏 [10−38cm2] 𝐹𝜉 [GeV] 𝜉𝜈𝑞 → 𝜈−𝑞𝜌+, 𝑋𝜌𝑂 < 1.4 GeV HNV DCC 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 𝜏 [10−38cm2] 𝐹𝜉 [GeV] 𝜉𝜈𝑜 → 𝜈−𝑜𝜌+, 𝑋𝜌𝑂 < 1.4 GeV HNV DCC 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 𝜏 [10−38cm2] 𝐹𝜉 [GeV] 𝜉𝜈𝑜 → 𝜈−𝑞𝜌0, 𝑋𝜌𝑂 < 1.4 GeV HNV DCC 0.2 0.4 0.6 0.8 1 1.2 1 2 3 4 5 6 𝜏 [10−38cm2] 𝐹𝜉 [GeV] 𝜉𝜈𝑞 → 𝜈−𝑞𝜌+, 𝑋𝜌𝑂 < 2 GeV DCC 0.05 0.1 0.15 0.2 0.25 0.3 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 𝜏 [10−38cm2] 𝐹𝜉 [GeV] 𝜉𝜈𝑜 → 𝜈−𝑜𝜌+, 𝑋𝜌𝑂 < 2 GeV DCC 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 𝜏 [10−38cm2] 𝐹𝜉 [GeV] 𝜉𝜈𝑜 → 𝜈−𝑞𝜌0, 𝑋𝜌𝑂 < 2 GeV DCC

  • J. Sobczyk,E. Hernandez,S.X. Nakamura, J. Nieves,T. Sato PRD98(2018)073001

Re-analyzed ANL/BNL data, C. Wilkinson et al. PRD90 ANL-Osaka DCC,PRD92, Hernandez,Nieves,Valverde PRD76 Recently FSI effects(10 ∼ 30%) found by S.X.Nakamura et al.

  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

11 / 29

slide-12
SLIDE 12

Angular distribution of pion

dσCC dWπNdQ2dΩ∗

π

= G2

F WπN

4πMk2 [A∗ + B∗ cos ϕπ + C∗ cos 2ϕπ + D∗sin ϕπ + E∗sin 2ϕπ]

  • θ

θ*

π

Z* X*

φ*

π

Y *

’ q k k’ k π

PV: sin ϕπ ∼ k × k′ · kπ angular distribution is sensitive to reaction dynamics. (phase, res. and non-res.) parity violating T-odd angular distributions(D∗, E∗): due to final-state-interaction (D∗, E∗) are sensitive to interference among partial waves. In ∆ resonance region, lπ ≤ 1 E∗ ∝ sin2 θπ sin(δP33 − δP31)[|MV

1+||EA 1−| + |MV 1−|(4|MA 1+| + 2|EA 1+|)].

M1+, E1+ : p3/2, M1−, E1− : p1/2

  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

12 / 29

slide-13
SLIDE 13

Pion angular distribution in ∆(1232) region

Flux averaged angular distribution of pion: Comparison with ANL/BNL data

20 25 30 35 40 45 50 55 60 65 −1 −0.5 0.5 1 Nr of events cos 𝜄∗

𝜌

ANL 𝜉𝜈𝑞 → 𝜈−𝑞𝜌+ 30 35 40 45 50 55 60 50 100 150 200 250 300 350 Nr of events 𝜚∗

𝜌

60 80 100 120 140 160 180 −1 −0.5 0.5 1 Nr of events cos 𝜄∗

𝜌

BNL 𝜉𝜈𝑞 → 𝜈−𝑞𝜌+ HNV HNV1 HNV2 DCC 80 90 100 110 120 130 140 50 100 150 200 250 300 350 Nr of events 𝜚∗

𝜌

dσ/d cos θπ dσ/dϕπ

  • J. E. Sobczyk,E. Hernandez,S.X.Nakamura,J. Nieves, T. Sato, PRD98 073001(2018)
  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

13 / 29

slide-14
SLIDE 14

Pion production cross section in N ∗, ∆ region

dσ/dWπN of single pion production Eν = 40GeV

2 4 6 1 1.2 1.4 1.6 1.8 2 Axial Vector 2 4 6 1 1.2 1.4 1.6 1.8 2 0.5 1 1.5 2 1 1.2 1.4 1.6 1.8 2 0.5 1 1.5 2 1 1.2 1.4 1.6 1.8 2 0.5 1 1 1.2 1.4 1.6 1.8 2

+

0.5 1 1 1.2 1.4 1.6 1.8 2

+

P33 S11,D13 D15,P13

Neutrino anti-neutrino BEBC NP343,285(1990) νp: ∆(1232) νn:∆(1232) + higher resonance region

  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

14 / 29

slide-15
SLIDE 15

Pion angular distribution in N ∗, ∆ region

< Ylm >= ∫ dΩπY ∗

lm dσ dW dΩπ

∫ dΩπY ∗

00 dσ dW dΩπ Data:NPB343 (1990)D. Allasia et al.

ANL-Osaka Model (preliminary) ( ¯ νp → µ+pπ−, Eν = 20GeV ) dσ dWdΩπ = σ0 + σc cos ϕπ + σs sin ϕπ + σ2c cos 2ϕπ + σ2s sin 2ϕπ

  • 0.2

0.2 0.4 0.6 0.8 1 1.4 1.8 Y10

  • 0.2

0.2 0.4 0.6 0.8 1 1.4 1.8 Y20

  • 0.2

0.2 1 1.4 1.8 Re(Y11)

  • 0.2

0.2 1 1.4 1.8 Re(Y21)

  • 0.2

0.2 1 1.4 1.8 Im(Y11)

  • 0.2

0.2 1 1.4 1.8 Im(Y21)

  • 0.2

0.2 1 1.4 1.8 Re(Y22)

  • 0.2

0.2 1 1.4 1.8 Im(Y22)

Angular distribution depends on W, Q2. Need to examine models of neutrino event generators.

  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

15 / 29

slide-16
SLIDE 16

Axial Vector current of ANL-Osaka Model

Axial Vector current F CC

2

(total cross section) at Q2 = 0 DCC model : 1π dash, a Total solid πN cross section data: 1π green, a Total brown F CC

2

(Q2 = 0) = 2f2

π

π σ(π + N)

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 F2

CCn(Q2=0)

W (GeV) 1 2 3 1.2 1.4 1.6 1.8 2 F2

CCp(Q2=0)

W (GeV)

DCC(full) DCC(1pi)

Neutron Proton Description of axial vector current at Q2 = 0 is consistent with pion scattering data.

  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

16 / 29

slide-17
SLIDE 17

Adler’s sum rule (Q2 Dependence of axial vector current)preliminary

1 = [gA(Q2)]2 + ∫ ∞

νth

[W A

2,n(ν, Q2) − W A 2,p(ν, Q2)]dν

W A

2,n − W A 2,p = 2(W A 2,I=1/2 − W A 2,I=3/2)/3

gA(Q2): gA = 1.27, MA = 1.1GeV W A

2,n/p: DCC model. ν integration: from threshold up to W = 2GeV

  • 0.5

0.5 1 1.5 0.02 0.04 0.06 0.08 0.1 Q2(GeV)2 N Tot S.L. Adler,PR143(1965)1144,arXiv:0905.2923, E. A. Paschos, D. Schalla, PRD84(2011),013004.

  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

17 / 29

slide-18
SLIDE 18

Resonance to DIS

F2(x, Q2) of electromagnetic current

JJJJJJJJJJ JJJJJJ JJ J J JJJJJJ JJ JJ JJ JJJJJJJJJJJJJJJJ JJJJJJJJ

0.1 0.2 0.3 0.4 0.5 0.2 0.4 0.6 0.8 1

CLAS data

LO Parton Model Our RES Model DIS region W > 2GeV , Q2 > 1GeV 2. Reasonable description of transition region in terms of hadron picture(ANL-Osaka Model).

  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

18 / 29

slide-19
SLIDE 19

F2(x, Q2) of charged current(CC)

0.2 0.4 0.6 0.8 1.2 1.4 1.6 1.8 2 F2

p(Q2=1.5GeV2)

W(GeV)

0.2 0.4 0.6 0.8 1 1.2 1.4 1.2 1.4 1.6 1.8 2 F2

n(Q2=1.5GeV2)

Data: Z. Phys. C28 (1985) Allasia et al.

F p

2 (CC) = F2,I=3/2(CC)

F2,n(CC) = F2,I=3/2(CC) + 2F2,I=1/2(CC) 3 F V

2 (CC) = F A 2 (CC) in LO parton model

F V

2 (CC) >> F A 2 (CC) at large W, Q2 in DCC model

task: construct model of meson production reaction with F V

2 ∼ F A 2 around

W ∼ 2GeV, Q2 ∼ 1 − 2GeV 2

  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

19 / 29

slide-20
SLIDE 20

Summary

ANL-OSAKA DCC model is extended to describe weak meson production reaction up to W < 2GeV . Neutrino induced single pion production in N∗, ∆ resonance region is studied using ANL-Osaka model. Phase among partial waves/ interplay between non-resonant and resonant mechanism are important for angular distribution of pion. Feature of data can be understood within ANL-Osaka model. Comparison with Neutrino event generators(NEUT, GENIE,NuWro,..) and other models will be very useful. Model of axial vector current is examined. At Q2 = 0, DCC model reproduce πN data. Adler sum rule provides a consistency test on description of axial current in resonance region with finite Q2. in progress: improve Q2-dependence of axial current needs new data on νd and νp with improved statistics.

  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

20 / 29

slide-21
SLIDE 21

Back up

  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

21 / 29

slide-22
SLIDE 22

Res Non-res Unit. 1pi 2pi Tot RS

Delta,N*

  • X

O O LPP

Delta,N*

X

X O O HVM

Delta(1232)

chiral

O O

Delta(1232)+N(1440)

chiral

X O O Giessen

Delta, N*

phen.

X O O ANL-Osaka

Delta, N*

O

O O O O Summary of models for neutrino reaction in RES

RS: D. Rein, L. M. Sehgal AP133(81), LPP: O. Lalakulich,E.A. Paschos,G. Piranlshvili,PRD74(2006) HNV: E. Hernandez,J. Nieves,M. Valverde PRD76(2007) Giessen: T. Leitner,O.Buss,L.Alvarez-Ruso,U. Mosel,PRC79(2009) ANL-Osaka DCC:S.X.Nakamura,H. Kamano,TS,PRD92(2015) ,TS,D. Uno,T.-S.H.Lee PRC67(2003)

  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

22 / 29

slide-23
SLIDE 23

single pion production

Figure 12: (left) Comparison of theoretical and event generator calculations available in 2014 with MiniBooNE ν

µCH2

CC π

+ production data [332] (right) Comparison of event generator calculations with MINERν

A ν

µCH CC π + data [289].

Calculations are from NEUT (solid line), GENIE (dotted line), and NuWro (dashed line).

Alvarez-Ruso, L. and others, NuSTEC White Paper

  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

23 / 29

slide-24
SLIDE 24

CCπ0

νµ + CH → µ−π0X Eν ∼ 3GeV

  • O. Altinok et al. PRD96(2018)072003
  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

24 / 29

slide-25
SLIDE 25

Model of neutrino reaction in resonance region

Generator

4

SIS/DIS region in the generators

W 2 GeV/c² Resonances

(1π, 1K, 1η)

+ DIS background

(“Multi-pi” mode)

PYTHIA 5.72 (“DIS” mode) Resonances

+ DIS background (“AGKY model”)

1.7 GeV/c² DIS low W

(“AGKY model”)

2.3 GeV/c² Linear transition to PYTHIA 6 3 GeV/c² PYTHIA 6 W GENIE 1.3 GeV/c² RES Linear transition 1.6 GeV/c² DIS (uses PYTHIA 6 fragmentation routines) W NuWro 1.3 GeV/c² NEUT

  • C. Bronner (Workshop on Shallow and DIS Scattering 2018)

resonance model (Isobar model: Rein-Sehgal Model)

  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

25 / 29

slide-26
SLIDE 26

Isobar model

N N W

Transition form factor pion coupling constant

R

Breit-Wigner formula for partial wave (JπI) gπNR gJNR W − MR + iΓR/2 (Rein-Shegal, O. Lalakulich, E. A. Paschos, Mass(MR), Width(ΓR) of resonance R from PDG. gπNR = √ Γ 2 Bπ, gV NR = √ Γ 2 Bγ gANR = PCAC or Quark model No control of relative phases (non-res vs res, (Jπ, I) channels.)

  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

26 / 29

slide-27
SLIDE 27

Brief explanation of Coupled-channel model( πN, ηN, KΛ, KΣ, ππN ) Physics included inside V

c

  • upled-

c hannels effec t

s-channel u-channel t-channel contact

N*

bare

Differential cross section and polarization data analyzed within the ANL-Osaka Model πp → πN, ηN, KΛ, KΣ γp → πN, ηN, KΛ, KΣ ep → e′πN γn: needed for CC, NC to separate iso-vector from iso-scalar component. Axial vector current: gNN∗

A

from gNN∗

π

assuming PCAC and dipole form factor. (Only in ∆ region, the model can tested by data.)

  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

27 / 29

slide-28
SLIDE 28

High Q2 region (F CC

2

= F CC

2,p + F CC 2,n

2 )

(Left) F em

2

× 18/5 of ANL-Osaka Model. F CC

2

≈ 18 5 F em

2

∼ x(u + ¯ u + d + ¯ d) PDF

0.5 1 1.5 2 0.2 0.4 0.6 0.8 1 Q2=0.5 Q2=2 F

emx 18/5 2

0.5 1 1.5 2 0.2 0.4 0.6 0.8 1 FCC

2

(Right)F CC

2

  • f ANL-Osaka Model

F CC

2

∝ |VIV |2 + |AIV |2 DCC

  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

28 / 29

slide-29
SLIDE 29

Possibility to study ACC

µ

: Parity violating asymmetry N(⃗ e, e′)X

AP V = σ+ − σ− σ+ + σ− = − Q2GF √ 24πα cos2 θ

2 W γZ 2

+ sin2 θ

2 [2W γZ 1

+ (1 − 4 sin2 θW ) Ee+E′

e

MN

W γZ

3

] cos2 θ

2 W em 2

+ sin2 θ

2 W em 1

d(⃗ e, e′)X,Ee = 4.867GeV, θ = 12.9o (PVDIS PRC91 045506 (2015))

1.0 1.2 1.4 1.6 1.8 2.0

)

2

(GeV

2

(ppm) / Q

pv

A

−150 −100 −50 =4.867GeV

b

E

W(GeV)

1.86 1.98 2.10 =6.067GeV

b

E Data I Data III Data IV Data II

  • T. Sato (Osaka U.)

N(ν, lM)N

  • Nov. 12 2018

29 / 29