Triton-like molecules Li Ma 13/01/2020 Introduction and Motivation - - PowerPoint PPT Presentation

triton like molecules
SMART_READER_LITE
LIVE PREVIEW

Triton-like molecules Li Ma 13/01/2020 Introduction and Motivation - - PowerPoint PPT Presentation

Triton-like molecules Li Ma 13/01/2020 Introduction and Motivation 1 2 Tri-baryon systems , and 3 Application to the NNN system 4 Numerical results for the , and Some other numerical results 5


slide-1
SLIDE 1

Triton-like molecules

Li Ma 13/01/2020

slide-2
SLIDE 2

1

Introduction and Motivation

2

Tri-baryon systems ΛΛΛ, ΞΞΞ and ΣΣΣ

3

Application to the NNN system

4

Numerical results for the ΛΛΛ, ΞΞΞ and ΣΣΣ

5

Some other numerical results

6

Summary

13/01/2020 2 / 20

slide-3
SLIDE 3

Introduction and Motivation

The deuteron-like molecules

π η σ ρω

  • Deuteron

π η σ ρω

  • Di-meson molecule

π η σ ρω

  • Di-baryon molecule

X(3872) → D ¯ D∗, Ds1(2460) → D∗K, Zb(10610) → B ¯ B∗ Y (4260) → D1 ¯ D, Zc(3900) → D ¯ D∗, Zc(4020) → D∗ ¯ D∗ Pc(4380) → Σc ¯ D∗, Pc(4450) → Σ∗

c ¯

D∗

13/01/2020 3 / 20

slide-4
SLIDE 4

Introduction and Motivation

The triton-like molecules

π η σ ρ ω

  • Triton
  • Benzene ring

NDK, ¯ KDN, ND ¯ D, N ¯ KK, DKK, DK ¯ K, BBB∗, DDK, DD∗K, ΩNN and ΩΩN Faddeev equation, or FCA Dimer or isobar formalism Gaussian expansion method (GEM)

13/01/2020 4 / 20

slide-5
SLIDE 5

Introduction and Motivation

Efimov effect

B2 =

2 ma2 (1 + O( r0 a ))

B3(1 + O( r0

a )) =

− 2

ma2 + [e−2πnf (ξ)]

1 s0 κ2 ∗

m

s0 = 1.00624..., ξ is defined by tan ξ = −(mB3)1/2a/, f (ξ)is a universal function with f (−π/2) = 1, K = sgn(E)

  • m|E|.

Figure: Efimov plot

13/01/2020 5 / 20

slide-6
SLIDE 6

Tri-baryon systems ΛΛΛ, ΞΞΞ and ΣΣΣ

Tri-baryon systems ΛΛΛ, ΞΞΞ and ΣΣΣ

The OBE interaction indicates that there is only one virtual boson exchanged by any two constituents as shown in the following.

π Ξ

  • Ξ

Ξ π Λ

  • Λ

Λ π Σ

  • Σ

Σ

Figure: Dynamical illustration of the ΛΛΛ, ΞΞΞ and ΣΣΣ systems with a circle describing the delocalized π bond inside.

Monopole form factor F(q) = Λ2−m2

α

Λ2−q2 , (α = π, η, ρ, ω, σ, φ) with q the

four-momentum of the pion and Λ the cutoff parameter.

13/01/2020 6 / 20

slide-7
SLIDE 7

Tri-baryon systems ΛΛΛ, ΞΞΞ and ΣΣΣ

Born-Oppenheimer potential

Considering that the particle b and c are static with the separation rbc, one can separate the degree of freedom of a from the three-body system. We assume the distance rbc is a parameter. The mesons b and c are static, and have one-pion interactions with meson a, which can be viewed as two static sources. We explore the dynamics for the meson a in the limit rbc → ∞, and subtract the binding energy for the break-up state which is trivial for the three-body bound state.

π π fixed fixed Ξ Ξ

  • Ξ

()

VBO(r bc)

Ξ Ξ

Ξ

  • ()

Figure: Illustration of the BO potential. (a) illustrates the calculation procedure

13/01/2020 7 / 20

slide-8
SLIDE 8

Tri-baryon systems ΛΛΛ, ΞΞΞ and ΣΣΣ

Interpolating wave function of meson a

fixed fixed Ξ Ξ

  • Ξ

ψ ( )

ψ ( )

The zero order of the final wave function for the meson a could be the superposition of these two components ψ( rab, rac) = C[ψ( rab) ± ψ( rac)]|ΞΞΞ Accordingly, one can obtain the energy eigenvalue of the meson a Ea(Λ, rbc) = ψ( rab, rac)|Ha|ψ( rab, rac)

13/01/2020 8 / 20

slide-9
SLIDE 9

Tri-baryon systems ΛΛΛ, ΞΞΞ and ΣΣΣ

BO potential and Its physical meaning (Intensity of ”glue”)

We define the BO potential as VBO(Λ, rbc) = Ea(Λ, rbc) − E2(Λ). The BO potential can describe the contribution for the

  • ne meson on the dynamics of the two remaining
  • mesons. The meson a here works like a mass of ”glue”.

Ξ

  • Ξ

Ξ

0.5 1.0 1.5 2.0 2.5 3.0 r/fm

  • 40
  • 20

20 40 V(r)/MeV

Vϕss Vωss Vρss Vσss Vηss Vπss Vtot

2 4 6 8 10

  • 7
  • 6
  • 5
  • 4
  • 3
  • 2
  • 1

Figure: Here we chose the parameter Λ = 900 MeV. E ΞΞ

I=0 = −3.09 MeV.

ΨT = αΦ( rbc)ψ( rab, rac).

13/01/2020 9 / 20

slide-10
SLIDE 10

Tri-baryon systems ΛΛΛ, ΞΞΞ and ΣΣΣ

The configurations of the three-body systems

Ξ

  • Ξ

Ξ

  • Ξ

Ξ Ξ

  • Ξ

Ξ Ξ

Figure: Every meson can be considered to be a lighter one and separated from the three-body system. Each of them can generate the ”glue” for the remaining mesons.

VBO(r bc)

  • ()

Ξ Ξ Ξ

VBO(r ac)

  • ()

Ξ Ξ Ξ

VBO(r ab)

  • ()

Ξ Ξ Ξ

Figure: (a), (b) and (c) correspond to the wave functions ψ/

a = Φ(

rbc)ψ( rab, rac), ψ/

b = Φ(

rac)ψ( rab, rbc) and ψ/

c = Φ(

rab)ψ( rbc, rac), respectively.

The basis constitute a configuration space {ψ , ψ , ψ }.

13/01/2020 10 / 20

slide-11
SLIDE 11

Tri-baryon systems ΛΛΛ, ΞΞΞ and ΣΣΣ

Interpolating wave functions

The basis constitute a configuration space {ψ/

a, ψ/ b, ψ/ c}.

ΨT = αΦ( rbc)ψ( rab, rac) + βΦ( rac)ψ( rab, rbc) + γΦ( rab)ψ( rbc, rac) = αψ/

a + βψ/ b + γψ/ c =

  α β γ   , Expand Φ( rbc), Φ( rac) and Φ( rab) as a set of Laguerre polynomials χnl(r) =

  • (2λ)2l+3n!

Γ(2l + 3 + n)r le−λrL2l+2

n

(2λr), n = 1, 2, 3... ψ/

a =

  • i

φi( rbc)ψ( rab, rac), ψ/

b =

  • i

φi( rac)ψ( rab, rbc), ψ/

c =

  • i

φi( rab)ψ( rbc, rac).

13/01/2020 11 / 20

slide-12
SLIDE 12

Tri-baryon systems ΛΛΛ, ΞΞΞ and ΣΣΣ

Orthonormalization

We orthonormalize the {ψ/

a, ψ/ b, ψ/ c} into a new

basis { ˜ ψ/

a, ˜

ψ/

b, ˜

ψ/

c}.

˜ ψi

/ a

= 1 Ni

  • (ψi

/ a + ψi / b + ψi / c) −

  • i

xijψj

/ a

  • ,

˜ ψi

/ b

= 1 Ni

  • (ψi

/ a + ψi / b + ψi / c) −

  • i

xijψj

/ b

  • ,

˜ ψi

/ c

= 1 Ni

  • (ψi

/ a + ψi / b + ψi / c) −

  • i

xijψj

/ c

  • ,

ψ

  • ψ
  • ψ
  • ψ

+ ψ + ψ

  • ψ

ψ

ψ

where the xij is a parameter matrix which will be determined later. The Ni are normalization coefficients. Then the eigenvector for the three-body system B(∗)

a B(∗) b B(∗) c

can be written as a vector in the configuration space { ˜ ψ/

a, ˜

ψ/

b, ˜

ψ/

c}. Therefore, we have

ΨT =

  • i

˜ αi ˜ ψi

/ a +

  • j

˜ βj ˜ ψj

/ b +

  • k

˜ γk ˜ ψk

/ c ,

13/01/2020 12 / 20

slide-13
SLIDE 13

Tri-baryon systems ΛΛΛ, ΞΞΞ and ΣΣΣ

First order correction

fixed fixed Ξ Ξ

  • Ξ

ψ *  ( )

ψ *  ( )

The first order of the final wave function for the meson a could be the superposition of these two components ˜ ψ∗( rab, rac) = C ∗[ ˜ ψ∗( rab) ± ˜ ψ∗( rac)]|ΞΞΞ Accordingly, one can obtain the energy eigenvalue of the meson a E ∗

a (Λ,

rbc) = ˜ ψ( rab, rac)|Ha| ˜ ψ( rab, rac) We define the BO potential as V ∗

BO(Λ,

rbc) = E ∗

a (Λ,

rbc) − E2(Λ).

13/01/2020 13 / 20

slide-14
SLIDE 14

Application to the NNN system

Application to the NNN system (Triton or Helium-3 nucleus)

  • =

[]

=/

[] EI=1/2

NNN=1.71 MeV

EI=1/2

NNN=5.38 MeV

π

  • Figure: Dependence of the reduced three-body binding energy on the binding

energy of its two-body subsystem (the deuteron). The result is comparable with the empirical binding energies of the triton (8.48 MeV) and helium-3 (7.80 MeV) nuclei.

13/01/2020 14 / 20

slide-15
SLIDE 15

Application to the NNN system

Numerical results for the NNN system (Triton or Helium-3 nucleus)

Λ(MeV) E2(MeV) E3(MeV) ET (MeV) VBO(0)(MeV) S wave(%) D wave(%) rrms(fm) 830.00

  • 0.18
  • 1.93
  • 2.11
  • 4.54

94.01 5.99 4.21 850.00

  • 0.67
  • 2.71
  • 3.38
  • 5.36

93.36 6.64 4.00 870.00

  • 1.23
  • 3.65
  • 4.88
  • 6.32

92.68 7.32 3.78 890.00

  • 1.88
  • 4.77
  • 6.66
  • 7.42

91.99 8.01 3.54 899.60

  • 2.23
  • 5.38
  • 7.62
  • 8.00

91.66 8.34 3.42 900.00

  • 2.25
  • 5.41
  • 7.66
  • 8.03

91.64 8.36 3.42 920.00

  • 3.05
  • 6.85
  • 9.90
  • 9.35

90.97 9.03 3.18 940.00

  • 3.98
  • 8.51
  • 12.49
  • 10.83

90.35 9.65 2.95 960.00

  • 5.03
  • 10.42
  • 15.45
  • 12.46

89.76 10.24 2.74 980.00

  • 6.21
  • 12.57
  • 18.78
  • 14.23

89.23 10.77 2.54 1000.00

  • 7.55
  • 14.97
  • 22.51
  • 16.14

88.73 11.27 2.37 1020.00

  • 9.04
  • 17.61
  • 26.65
  • 18.19

88.27 11.73 2.23 1040.00

  • 10.69
  • 20.51
  • 31.20
  • 20.37

87.84 12.16 2.10

Table: Bound state solutions for the NNN system with isospin I3 = 1/2. E2 is the energy eigenvalue of its subsystem. E3 is the reduced three-body energy eigenvalue relative to the break-up state of the NNN system. ET is the total three-body energy eigenvalue relative to the NNN threshold.

13/01/2020 15 / 20

slide-16
SLIDE 16

Numerical results for the ΛΛΛ, ΞΞΞ and ΣΣΣ

Numerical results for the ΛΛΛ, ΞΞΞ and ΣΣΣ

  • =

ΞΞ []

=/

ΞΞΞ [] EI=1/2

ΞΞΞ =1.06 MeV

EI=1/2

ΞΞΞ =10.76 MeV

  • =

ΞΞ []

=/

ΞΞΞ [] EI=3/2

ΞΞΞ =1.12 MeV

EI=3/2

ΞΞΞ =11.25 MeV

  • =

ΞΞ []

=/

ΞΞΞ [] EI=1/2

ΞΞΞ =1.06 MeV

  • =

ΞΞ []

=/

ΞΞΞ [] EI=3/2

ΞΞΞ =1.12 MeV

13/01/2020 16 / 20

slide-17
SLIDE 17

Numerical results for the ΛΛΛ, ΞΞΞ and ΣΣΣ

Numerical results for the ΛΛΛ, ΞΞΞ and ΣΣΣ

π

  • =
  • ( )=(/)

Ξ Ξ Ξ

ΞΞΞ

π

  • =
  • ( )=(/)

Ξ Ξ Ξ

ΞΞΞ

π

  • =

( )=()

Σ Σ Σ

ΣΣΣ

uΞΞ

S

uΞΞ

D

uΞΞΞ

S

uΞΞΞ

D
  • []

()=Ψ()* (a)

  • =

ΣΣ []

=

ΣΣΣ[] EI=1

ΣΣΣ=1.07 MeV

EI=1

ΣΣΣ=11.71 MeV
  • =

ΣΣ []

=

ΣΣΣ[] EI=1

ΣΣΣ=1.07 MeV

13/01/2020 17 / 20

slide-18
SLIDE 18

Some other numerical results

Numerical results for the Doubly heavy tri-meson bound states

π (*) (*) (*)

  • π

(*) (*) (*)

  • uDD*
S

uDD

S
  • []

()=Ψ()*

(a)

7.92 MeV

  • =

* []

=/

* []

(b)

1.98 MeV

  • =

*

_

[] =/

*

_

[]

(MeV) DD∗K BB∗ ¯ K D ¯ D∗K B ¯ B∗ ¯ K Binding energy 8.29+4.32

−3.66

41.76+8.84

−8.49

8.29+6.55

−6.13

41.76+9.02

−8.68

Mass 4317.92+3.66

−4.32

11013.65+8.49

−8.84

4317.92+6.13

−6.55

11013.65+8.68

−9.02 13/01/2020 18 / 20

slide-19
SLIDE 19

Some other numerical results

Numerical results for the tri-meson bound state BBB∗

  • =

*[]

=/

*[] EI=3/2

BBB*=0.06 MeV

(a)

  • =

*[]

=/

*[] EI=1/2

BBB*=0.12 MeV

(b)

uBB*

S

uBB*

D

uBBB

* S

uBBB

* D
  • []

()=Ψ()* (a)

uBB*

S

uBB*

D

uBBB

* S

uBBB

* D
  • []

()=Ψ()* (b)

π (*) (*) (*)

  • =

=/

Figure: Here we chose the parameter Λ = 1440 MeV in (a) and Λ = 1107.7 MeV in (b) for a better comparison of all the cases, since they have the same two-body binding energy of 5.08 MeV.

13/01/2020 19 / 20

slide-20
SLIDE 20

Summary

Summary

We predict that the triton-like molecular states for the ΞΞΞ and ΣΣΣ systems are probably existent as long as the molecular states of their two-body subsystems exist. In our calculations, we use the Born-Oppenheimer potential method to construct our interpolating wave functions, which can be regarded as a version of the variational principle which always gives an upper limit of the energy of a system. Other configurations NNΛ, NΛΛ, NNΞ, NΞΞ, NNΣ, NΣΣ will be explored in a future work.

Thank you very much!

13/01/2020 20 / 20