Spin entanglement in two-proton emission from 6 Be
Abstract
This paper presents an evaluation of coupled-spin entanglement in the two-proton (2 ⢠p 2 đ 2p 2 italic_p ) radioactive emission.
The three-body model of 6 Be with the proton-proton interaction, which is adjusted to reproduce the experimental energy release, is utilized.
Time-dependent calculation is performed to compute the coupled-spin state of the emitted two protons.
The spin-correlation function S đ S italic_S as Clauser-Horne-Shimony-Holt (CHSH) indicator is obtained as | S | â
2.72 đ 2.72 \left|S\right|\cong 2.72 | italic_S | â
2.72 .
Namely, the 2 ⢠p 2 đ 2p 2 italic_p -spin entanglement beyond the limit of local-hidden-variable (LHV) theory is concluded.
This entanglement is sensitive to the proton-proton interaction.
The short-lived (broad-width) 2 ⢠p 2 đ 2p 2 italic_p  state has the weaker spin entanglement.
In parallel, the core-proton interactions do not harm this entanglement during the time-dependent decaying process.
The CHSH measurement can be a probe into the effective nuclear interaction inside finite systems.
Introduction - Clauser-Horne-Shimony-Holt (CHSH) inequality is one essential property of the quantum-entangled state Clauser et al. (1969 ) .
This is one variant of Bell inequality introduced by John Clauser et. al. for the proof of Bell theorem, which claims that certain consequence of quantum mechanics cannot be reproduced by the local-hidden-variable (LHV) theory Bell (1964 , 2004 ) .
By using the CHSH indicator S đ S italic_S , the limit of LHV theory is symbolically given as | S | ⤠2 đ 2 \left|S\right|\leq 2 | italic_S | ⤠2 .
Throughout the history of Bell-CHSH examinations Clauser and Shimony (1978 ); Aspect et al. (1982 ); Weihs et al. (1998 ); Pan et al. (2000 ); Rowe et al. (2001 ); Scheidl et al. (2010 ); Giustina et al. (2015 ); Shalm et al. (2015 ); Vasilyev et al. (2017 ); Storz et al. (2023 ) , the violation of LHV-theory limit has been confirmed.
In these examinations, the entangled states of photons, electrons, and atoms are populated and measured to satisfy | S | > 2 đ 2 \left|S\right|>2 | italic_S | > 2 .
I emphasize that, for this purpose, many efforts have been devoted to close the loopholes, including those of detection Rowe et al. (2001 ); Giustina et al. (2015 ); Shalm et al. (2015 ) , locality Aspect et al. (1982 ); Weihs et al. (1998 ); Scheidl et al. (2010 ) , and memory Barrett et al. (2002 ) .
In the nuclear physics, the quantum entanglement plays a role in various
scenes Lamehi-Rachti and Mittig (1976 ); Sakai et al. (2006 ); Kanada-Enâyo (2015 ); Bulgac (2023 ); Pazy (2023 ); Gu et al. (2023 ); Hengstenberg et al. (2023 ); PĂŠrez-Obiol et al. (2023 ); Kou et al. (2024 ); Kirchner et al. (2024 ); Robin et al. (2021 ); Kruppa et al. (2021 ); Miller (2023a , b ); Tichai et al. (2023 ); Sun et al. (2023 ) . In Ref. Lamehi-Rachti and Mittig (1976 ) by Lamehi-Rachti and Mittig, the low-energy proton-proton scattering was measured for testing the Bell inequality.
In this pioneer work, some extra assumptions were required to certify the violation of LHV-theory limit.
In Ref. Sakai et al. (2006 ) , Sakai et. al. performed a novel measurement to demonstrate that a strong entanglement is realized between two protons.
This experiment measures the spin-singlet two protons made in the reaction of 2 H+p âś 2 superscript âś 2 đ absent p\longrightarrow^{2} italic_p âś start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT He+n đ n italic_n .
The spin-correlation function as the CHSH indicator is deduced as S expt = 2.83 Âą 0.24 stat Âą 0.07 sys subscript đ expt plus-or-minus 2.83 subscript 0.24 stat subscript 0.07 sys S_{\rm expt}=2.83\pm 0.24_{\rm stat}\pm 0.07_{\rm sys} italic_S start_POSTSUBSCRIPT roman_expt end_POSTSUBSCRIPT = 2.83 Âą 0.24 start_POSTSUBSCRIPT roman_stat end_POSTSUBSCRIPT Âą 0.07 start_POSTSUBSCRIPT roman_sys end_POSTSUBSCRIPT , which is in agreement with the non-local quantum mechanics and beyond the LHV-theory limit.
Recently, several theoretical works have been devoted to compute the entanglement entropy in atomic nuclei Kanada-Enâyo (2015 ); Bulgac (2023 ); Pazy (2023 ); Gu et al. (2023 ); Hengstenberg et al. (2023 ); PĂŠrez-Obiol et al. (2023 ); Kou et al. (2024 ); Kirchner et al. (2024 ) .
In Ref. Bulgac (2023 ) , due to the nuclear short-range correlations, the occupation probabilities of nuclear orbits change to increase the entanglement entropy.
In Ref. Bai (2024 ) , the evaluation of nuclear spin entanglement with the quantum-state tomography is suggested to be feasible.
In Refs. Johnson and Gorton (2023 ); Miller (2023b ) , a finite proton-neutron entanglement is suggested.
Another possible example to observe the nuclear quantum entanglement is the two-proton (2 ⢠p 2 đ 2p 2 italic_p ) radioactive emission Bertulani et al. (2008 ); Grigorenko (2009 ); PfĂźtzner et al. (2012 ); Blank and Ploszajczak (2008 ); Blank and Borge (2008 ); PfĂźtzner et al. (2023 ); Qi et al. (2019 ) .
In this radioactive decay, the parent nuclei spontaneously decay by emitting two protons.
Especially in so-called âpromptâ 2 ⢠p 2 đ 2p 2 italic_p  emission Grigorenko (2009 ); Goldansky (1960 , 1961 ) ,
the two protons are expected to have the diproton-like clustering and/or the dominant spin-singlet configuration.
This is attributable to the effective 2 ⢠p 2 đ 2p 2 italic_p  interaction inside finite systems, being in a contrast to the vacuum 2 ⢠p 2 đ 2p 2 italic_p  interaction supporting no bound state.
In this paper, I evaluate the quantum entaglement in the 2 ⢠p 2 đ 2p 2 italic_p  emission of 6 Be.
Main motivation is to utilize this entanglement as a probe into the effective interaction.
The 6 Be is the lightest 2 ⢠p 2 đ 2p 2 italic_p  emitter and well approximated with the simple three-body model with time dependence Oishi et al. (2014 , 2017 ); Wang and Nazarewicz (2021 ) .
For the two protons spontaneously emitted, the measurement of their coupled-spin correlation is assumed Bertulani et al. (2008 ) .
For quantitizing the entanglement, the CHSH indicator S đ S italic_S ,
which was originally introduced for testing the CHSH inequality Clauser et al. (1969 ); Sakai et al. (2006 ) , is computed.
Because of the three-body problem, the 2 ⢠p 2 đ 2p 2 italic_p  entanglement is under the effect from the third particle,
namely, the daughter alpha nucleus.
Whether this effect destroys or not the entanglement is investigated.
The sensitivity of CHSH indicator to the proton-proton interaction is also discussed.
Formalism and Model - I consider the coupling of two protons, i.e., identical spin-1 / 2 1 2 1/2 1 / 2 fermions.
In such a case, the CHSH indicator is represented as follows.
First the four options of measurement by the two observers, so-called âAliceâ and âBobâ by convenance, are introduced:
A ^ 1 , θ ⢠( 1 ) â B ^ 1 , θ ⢠( 2 ) , A ^ 2 , θ ⢠( 1 ) â B ^ 1 , θ ⢠( 2 ) , tensor-product subscript ^ đ´ 1 đ
1 subscript ^ đľ 1 đ
2 tensor-product subscript ^ đ´ 2 đ
1 subscript ^ đľ 1 đ
2
\displaystyle\hat{A}_{1,\theta}(1)\otimes\hat{B}_{1,\theta}(2),~{}\hat{A}_{2,%
\theta}(1)\otimes\hat{B}_{1,\theta}(2), over^ start_ARG italic_A end_ARG start_POSTSUBSCRIPT 1 , italic_θ end_POSTSUBSCRIPT ( 1 ) â over^ start_ARG italic_B end_ARG start_POSTSUBSCRIPT 1 , italic_θ end_POSTSUBSCRIPT ( 2 ) , over^ start_ARG italic_A end_ARG start_POSTSUBSCRIPT 2 , italic_θ end_POSTSUBSCRIPT ( 1 ) â over^ start_ARG italic_B end_ARG start_POSTSUBSCRIPT 1 , italic_θ end_POSTSUBSCRIPT ( 2 ) ,
A ^ 1 , θ ⢠( 1 ) â B ^ 2 , θ ⢠( 2 ) , A ^ 2 , θ ⢠( 1 ) â B ^ 2 , θ ⢠( 2 ) , tensor-product subscript ^ đ´ 1 đ
1 subscript ^ đľ 2 đ
2 tensor-product subscript ^ đ´ 2 đ
1 subscript ^ đľ 2 đ
2
\displaystyle\hat{A}_{1,\theta}(1)\otimes\hat{B}_{2,\theta}(2),~{}\hat{A}_{2,%
\theta}(1)\otimes\hat{B}_{2,\theta}(2), over^ start_ARG italic_A end_ARG start_POSTSUBSCRIPT 1 , italic_θ end_POSTSUBSCRIPT ( 1 ) â over^ start_ARG italic_B end_ARG start_POSTSUBSCRIPT 2 , italic_θ end_POSTSUBSCRIPT ( 2 ) , over^ start_ARG italic_A end_ARG start_POSTSUBSCRIPT 2 , italic_θ end_POSTSUBSCRIPT ( 1 ) â over^ start_ARG italic_B end_ARG start_POSTSUBSCRIPT 2 , italic_θ end_POSTSUBSCRIPT ( 2 ) ,
(1)
Namely, Alice observes the first fermion with one choosen from the two options, A ^ 1 , θ subscript ^ đ´ 1 đ
\hat{A}_{1,\theta} over^ start_ARG italic_A end_ARG start_POSTSUBSCRIPT 1 , italic_θ end_POSTSUBSCRIPT and A ^ 2 , θ subscript ^ đ´ 2 đ
\hat{A}_{2,\theta} over^ start_ARG italic_A end_ARG start_POSTSUBSCRIPT 2 , italic_θ end_POSTSUBSCRIPT .
Bob does the second fermion with one choosen from B ^ 1 , θ subscript ^ đľ 1 đ
\hat{B}_{1,\theta} over^ start_ARG italic_B end_ARG start_POSTSUBSCRIPT 1 , italic_θ end_POSTSUBSCRIPT and B ^ 2 , θ subscript ^ đľ 2 đ
\hat{B}_{2,\theta} over^ start_ARG italic_B end_ARG start_POSTSUBSCRIPT 2 , italic_θ end_POSTSUBSCRIPT .
Those operators including the parameter angle θ đ \theta italic_θ are given as
A ^ 1 = Ď ^ z , A ^ 2 , θ = Ď ^ z ⢠cos ⥠2 ⢠θ + Ď ^ x ⢠sin ⥠2 ⢠θ , formulae-sequence subscript ^ đ´ 1 subscript ^ đ đ§ subscript ^ đ´ 2 đ
subscript ^ đ đ§ 2 đ subscript ^ đ đĽ 2 đ \hat{A}_{1}=\hat{\sigma}_{z},~{}~{}~{}~{}\hat{A}_{2,\theta}=\hat{\sigma}_{z}%
\cos 2\theta+\hat{\sigma}_{x}\sin 2\theta, over^ start_ARG italic_A end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = over^ start_ARG italic_Ď end_ARG start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT , over^ start_ARG italic_A end_ARG start_POSTSUBSCRIPT 2 , italic_θ end_POSTSUBSCRIPT = over^ start_ARG italic_Ď end_ARG start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT roman_cos 2 italic_θ + over^ start_ARG italic_Ď end_ARG start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT roman_sin 2 italic_θ ,
(2)
for Alice, whereas
B ^ 1 , θ = Ď ^ z ⢠cos ⥠θ + Ď ^ x ⢠sin ⥠θ , B ^ 2 , θ = Ď ^ z ⢠cos ⥠θ â Ď ^ x ⢠sin ⥠θ , subscript ^ đľ 1 đ
subscript ^ đ đ§ đ subscript ^ đ đĽ đ subscript ^ đľ 2 đ
subscript ^ đ đ§ đ subscript ^ đ đĽ đ \begin{array}[]{l}\hat{B}_{1,\theta}=\hat{\sigma}_{z}\cos\theta+\hat{\sigma}_{%
x}\sin\theta,\\
\hat{B}_{2,\theta}=\hat{\sigma}_{z}\cos\theta-\hat{\sigma}_{x}\sin\theta,\end{array} start_ARRAY start_ROW start_CELL over^ start_ARG italic_B end_ARG start_POSTSUBSCRIPT 1 , italic_θ end_POSTSUBSCRIPT = over^ start_ARG italic_Ď end_ARG start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT roman_cos italic_θ + over^ start_ARG italic_Ď end_ARG start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT roman_sin italic_θ , end_CELL end_ROW start_ROW start_CELL over^ start_ARG italic_B end_ARG start_POSTSUBSCRIPT 2 , italic_θ end_POSTSUBSCRIPT = over^ start_ARG italic_Ď end_ARG start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT roman_cos italic_θ - over^ start_ARG italic_Ď end_ARG start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT roman_sin italic_θ , end_CELL end_ROW end_ARRAY
(3)
for Bob.
For an arbitrary two-fermion state | Ψ ⢠( 1 , 2 ) ⊠ket Ψ 1 2 \ket{\Psi(1,2)} | start_ARG roman_Ψ ( 1 , 2 ) end_ARG ⊠, their expectation values are obtained:
⨠A i ⢠B j ⊠= ⨠Ψ ⢠( 1 , 2 ) ⢠⣠A ^ i , θ ⢠( 1 ) â B ^ j , θ ⢠( 2 ) ⣠⢠Ψ ⢠( 1 , 2 ) ⊠. expectation subscript đ´ đ subscript đľ đ expectation Ψ 1 2 delimited-âŁâŁ tensor-product subscript ^ đ´ đ đ
1 subscript ^ đľ đ đ
2 Ψ 1 2 \Braket{A_{i}B_{j}}=\Braket{\Psi(1,2)\mid\hat{A}_{i,\theta}(1)\otimes\hat{B}_{%
j,\theta}(2)\mid\Psi(1,2)}. ⨠start_ARG italic_A start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_B start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_ARG ⊠= ⨠start_ARG roman_Ψ ( 1 , 2 ) ⣠over^ start_ARG italic_A end_ARG start_POSTSUBSCRIPT italic_i , italic_θ end_POSTSUBSCRIPT ( 1 ) â over^ start_ARG italic_B end_ARG start_POSTSUBSCRIPT italic_j , italic_θ end_POSTSUBSCRIPT ( 2 ) ⣠roman_Ψ ( 1 , 2 ) end_ARG ⊠.
(4)
Then the CHSH indicator S đ S italic_S is determined as
S = max ⥠{ | S â ⣠+ ⣠+ + | , | S + ⣠â ⣠+ + | , | S + ⣠+ ⣠â + | , | S + ⣠+ ⣠+ â | } , đ subscript đ absent
subscript đ absent
subscript đ absent
subscript đ absent
S=\max\left\{\left|S_{-+++}\right|,\left|S_{+-++}\right|,\left|S_{++-+}\right|%
,\left|S_{+++-}\right|\right\}, italic_S = roman_max { | italic_S start_POSTSUBSCRIPT - + + + end_POSTSUBSCRIPT | , | italic_S start_POSTSUBSCRIPT + - + + end_POSTSUBSCRIPT | , | italic_S start_POSTSUBSCRIPT + + - + end_POSTSUBSCRIPT | , | italic_S start_POSTSUBSCRIPT + + + - end_POSTSUBSCRIPT | } ,
(5)
where
S â ⣠+ ⣠+ + subscript đ absent
\displaystyle S_{-+++} italic_S start_POSTSUBSCRIPT - + + + end_POSTSUBSCRIPT
= \displaystyle= =
â ⨠A 1 ⢠B 1 ⊠+ ⨠A 2 ⢠B 1 ⊠+ ⨠A 1 ⢠B 2 ⊠+ ⨠A 2 ⢠B 2 ⊠, expectation subscript đ´ 1 subscript đľ 1 expectation subscript đ´ 2 subscript đľ 1 expectation subscript đ´ 1 subscript đľ 2 expectation subscript đ´ 2 subscript đľ 2 \displaystyle-\Braket{A_{1}B_{1}}+\Braket{A_{2}B_{1}}+\Braket{A_{1}B_{2}}+%
\Braket{A_{2}B_{2}}, - ⨠start_ARG italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG ⊠+ ⨠start_ARG italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG ⊠+ ⨠start_ARG italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_B start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG ⊠+ ⨠start_ARG italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_B start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG ⊠,
S + ⣠â ⣠+ + subscript đ absent
\displaystyle S_{+-++} italic_S start_POSTSUBSCRIPT + - + + end_POSTSUBSCRIPT
= \displaystyle= =
⨠A 1 ⢠B 1 ⊠â ⨠A 2 ⢠B 1 ⊠+ ⨠A 1 ⢠B 2 ⊠+ ⨠A 2 ⢠B 2 ⊠, expectation subscript đ´ 1 subscript đľ 1 expectation subscript đ´ 2 subscript đľ 1 expectation subscript đ´ 1 subscript đľ 2 expectation subscript đ´ 2 subscript đľ 2 \displaystyle\Braket{A_{1}B_{1}}-\Braket{A_{2}B_{1}}+\Braket{A_{1}B_{2}}+%
\Braket{A_{2}B_{2}}, ⨠start_ARG italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG ⊠- ⨠start_ARG italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG ⊠+ ⨠start_ARG italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_B start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG ⊠+ ⨠start_ARG italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_B start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG ⊠,
S + ⣠+ ⣠â + subscript đ absent
\displaystyle S_{++-+} italic_S start_POSTSUBSCRIPT + + - + end_POSTSUBSCRIPT
= \displaystyle= =
⨠A 1 ⢠B 1 ⊠+ ⨠A 2 ⢠B 1 ⊠â ⨠A 1 ⢠B 2 ⊠+ ⨠A 2 ⢠B 2 ⊠, expectation subscript đ´ 1 subscript đľ 1 expectation subscript đ´ 2 subscript đľ 1 expectation subscript đ´ 1 subscript đľ 2 expectation subscript đ´ 2 subscript đľ 2 \displaystyle\Braket{A_{1}B_{1}}+\Braket{A_{2}B_{1}}-\Braket{A_{1}B_{2}}+%
\Braket{A_{2}B_{2}}, ⨠start_ARG italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG ⊠+ ⨠start_ARG italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG ⊠- ⨠start_ARG italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_B start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG ⊠+ ⨠start_ARG italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_B start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG ⊠,
S + ⣠+ ⣠+ â subscript đ absent
\displaystyle S_{+++-} italic_S start_POSTSUBSCRIPT + + + - end_POSTSUBSCRIPT
= \displaystyle= =
⨠A 1 ⢠B 1 ⊠+ ⨠A 2 ⢠B 1 ⊠+ ⨠A 1 ⢠B 2 ⊠â ⨠A 2 ⢠B 2 ⊠. expectation subscript đ´ 1 subscript đľ 1 expectation subscript đ´ 2 subscript đľ 1 expectation subscript đ´ 1 subscript đľ 2 expectation subscript đ´ 2 subscript đľ 2 \displaystyle\Braket{A_{1}B_{1}}+\Braket{A_{2}B_{1}}+\Braket{A_{1}B_{2}}-%
\Braket{A_{2}B_{2}}. ⨠start_ARG italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG ⊠+ ⨠start_ARG italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG ⊠+ ⨠start_ARG italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_B start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG ⊠- ⨠start_ARG italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_B start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG ⊠.
(6)
In FIG. 1 , the CHSH indicator for the spin-singlet state,
| Ψ ⢠( 1 , 2 ) ⊠= | â ⣠â ⣠â ⣠â ⣠â ⊠/ 2 ket Ψ 1 2 ket â â â â
2 \ket{\Psi(1,2)}=\ket{\uparrow\downarrow-\downarrow\uparrow}/\sqrt{2} | start_ARG roman_Ψ ( 1 , 2 ) end_ARG ⊠= | start_ARG â â - â â end_ARG ⊠/ square-root start_ARG 2 end_ARG , is presented.
At θ = Ď / 4 , 3 â˘ Ď / 4 , 5 â˘ Ď / 4 đ đ 4 3 đ 4 5 đ 4
\theta=\pi/4,~{}3\pi/4,~{}5\pi/4 italic_θ = italic_Ď / 4 , 3 italic_Ď / 4 , 5 italic_Ď / 4 , and 7 â˘ Ď / 4 7 đ 4 7\pi/4 7 italic_Ď / 4 , this indicator has the maximum value, S = 2 ⢠2 đ 2 2 S=2\sqrt{2} italic_S = 2 square-root start_ARG 2 end_ARG (Tsirelsonâs bound) Cirelâson (1980 ) .
This situation resembles Sakaiâs experiment Sakai et al. (2006 ) .
In the following sections, θ = Ď / 4 đ đ 4 \theta=\pi/4 italic_θ = italic_Ď / 4 , except when modified.
FIG. 1: CHSH indicator of the pure spin-singlet state.
As well known, when the initial state is one of the Bell states, | Ψ ⢠( 1 , 2 ) ⊠= | e n ⊠ket Ψ 1 2 ket subscript đ đ \ket{\Psi(1,2)}=\ket{e_{n}} | start_ARG roman_Ψ ( 1 , 2 ) end_ARG ⊠= | start_ARG italic_e start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT end_ARG ⊠, the CHSH indicator satisfies S = 2 ⢠2 đ 2 2 S=2\sqrt{2} italic_S = 2 square-root start_ARG 2 end_ARG when θ = Ď / 4 đ đ 4 \theta=\pi/4 italic_θ = italic_Ď / 4 .
Here the Bell states read Nielsen and Chuang (2000 )
| e 1 ⊠= 1 2 ⢠| â ⣠â ⣠+ ⣠â ⣠â ⊠, | e 2 ⊠= i 2 ⢠| â ⣠â ⣠â ⣠â ⣠â ⊠, formulae-sequence ket subscript đ 1 1 2 ket â â â â
ket subscript đ 2 đ 2 ket â â â â
\displaystyle\ket{e_{1}}=\frac{1}{\sqrt{2}}\ket{\uparrow\uparrow+\downarrow%
\downarrow},~{}~{}~{}~{}~{}\ket{e_{2}}=\frac{i}{\sqrt{2}}\ket{\uparrow\uparrow%
-\downarrow\downarrow}, | start_ARG italic_e start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG ⊠= divide start_ARG 1 end_ARG start_ARG square-root start_ARG 2 end_ARG end_ARG | start_ARG â â + â â end_ARG ⊠, | start_ARG italic_e start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG ⊠= divide start_ARG italic_i end_ARG start_ARG square-root start_ARG 2 end_ARG end_ARG | start_ARG â â - â â end_ARG ⊠,
| e 3 ⊠= i 2 ⢠| â ⣠â ⣠+ ⣠â ⣠â ⊠, | e 4 ⊠= 1 2 ⢠| â ⣠â ⣠â ⣠â ⣠â ⊠. formulae-sequence ket subscript đ 3 đ 2 ket â â â â
ket subscript đ 4 1 2 ket â â â â
\displaystyle\ket{e_{3}}=\frac{i}{\sqrt{2}}\ket{\uparrow\downarrow+\downarrow%
\uparrow},~{}~{}~{}~{}~{}\ket{e_{4}}=\frac{1}{\sqrt{2}}\ket{\uparrow\downarrow%
-\downarrow\uparrow}. | start_ARG italic_e start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT end_ARG ⊠= divide start_ARG italic_i end_ARG start_ARG square-root start_ARG 2 end_ARG end_ARG | start_ARG â â + â â end_ARG ⊠, | start_ARG italic_e start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT end_ARG ⊠= divide start_ARG 1 end_ARG start_ARG square-root start_ARG 2 end_ARG end_ARG | start_ARG â â - â â end_ARG ⊠.
(7)
These states can be used as basis to represent an arbitrary coupling of two spin-1 / 2 1 2 1/2 1 / 2 fermions.
In addition, there can be also coordinate degrees of freedom.
Thus, an arbitrary two-fermion state is generally expanded as
| Ψ ⢠( 1 , 2 ) ⊠= â i = 1 4 G i ⢠( đ 1 , đ 2 ) ⢠| e i ⊠, ket Ψ 1 2 superscript subscript đ 1 4 subscript đş đ subscript đ 1 subscript đ 2 ket subscript đ đ \ket{\Psi(1,2)}=\sum_{i=1}^{4}G_{i}(\bm{r}_{1},\bm{r}_{2})\ket{e_{i}}, | start_ARG roman_Ψ ( 1 , 2 ) end_ARG ⊠= â start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_G start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( bold_italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) | start_ARG italic_e start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG ⊠,
(8)
where G i ⢠( đ 1 , đ 2 ) subscript đş đ subscript đ 1 subscript đ 2 G_{i}(\bm{r}_{1},\bm{r}_{2}) italic_G start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( bold_italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) is the coordinate part.
For the coupled spin S ^ 12 subscript ^ đ 12 \hat{S}_{12} over^ start_ARG italic_S end_ARG start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT , its eigenvalues are given as S ^ 12 2 ⢠| d k ⊠= S ⢠( S + 1 ) ⢠| d k ⊠subscript superscript ^ đ 2 12 ket subscript đ đ đ đ 1 ket subscript đ đ \hat{S}^{2}_{12}\ket{d_{k}}=S(S+1)\ket{d_{k}} over^ start_ARG italic_S end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT | start_ARG italic_d start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_ARG ⊠= italic_S ( italic_S + 1 ) | start_ARG italic_d start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_ARG âŠ
and S ^ 12 , z ⢠| d k ⊠= V ⢠| d k ⊠subscript ^ đ 12 đ§
ket subscript đ đ đ ket subscript đ đ \hat{S}_{12,z}\ket{d_{k}}=V\ket{d_{k}} over^ start_ARG italic_S end_ARG start_POSTSUBSCRIPT 12 , italic_z end_POSTSUBSCRIPT | start_ARG italic_d start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_ARG ⊠= italic_V | start_ARG italic_d start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_ARG ⊠, where
| d 1 ⊠ket subscript đ 1 \displaystyle\ket{d_{1}} | start_ARG italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG âŠ
= \displaystyle= =
| S = 1 , V = + 1 ⊠= | â â ⊠, ket formulae-sequence đ 1 đ 1 ket â absent â \displaystyle\ket{S=1,V=+1}=\ket{\uparrow\uparrow}, | start_ARG italic_S = 1 , italic_V = + 1 end_ARG ⊠= | start_ARG â â end_ARG ⊠,
| d 2 ⊠ket subscript đ 2 \displaystyle\ket{d_{2}} | start_ARG italic_d start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG âŠ
= \displaystyle= =
| S = 1 , V = â 1 ⊠= | â â ⊠, ket formulae-sequence đ 1 đ 1 ket â absent â \displaystyle\ket{S=1,V=-1}=\ket{\downarrow\downarrow}, | start_ARG italic_S = 1 , italic_V = - 1 end_ARG ⊠= | start_ARG â â end_ARG ⊠,
| d 3 ⊠ket subscript đ 3 \displaystyle\ket{d_{3}} | start_ARG italic_d start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT end_ARG âŠ
= \displaystyle= =
| S = 1 , V = 0 ⊠= 1 2 ⢠| â ⣠â ⣠+ ⣠â ⣠â ⊠, ket formulae-sequence đ 1 đ 0 1 2 ket â â â â
\displaystyle\ket{S=1,V=0}=\frac{1}{\sqrt{2}}\ket{\uparrow\downarrow+%
\downarrow\uparrow}, | start_ARG italic_S = 1 , italic_V = 0 end_ARG ⊠= divide start_ARG 1 end_ARG start_ARG square-root start_ARG 2 end_ARG end_ARG | start_ARG â â + â â end_ARG ⊠,
| d 4 ⊠ket subscript đ 4 \displaystyle\ket{d_{4}} | start_ARG italic_d start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT end_ARG âŠ
= \displaystyle= =
| S = 0 , V = 0 ⊠= 1 2 ⢠| â ⣠â ⣠â ⣠â ⣠â ⊠. ket formulae-sequence đ 0 đ 0 1 2 ket â â â â
\displaystyle\ket{S=0,V=0}=\frac{1}{\sqrt{2}}\ket{\uparrow\downarrow-%
\downarrow\uparrow}. | start_ARG italic_S = 0 , italic_V = 0 end_ARG ⊠= divide start_ARG 1 end_ARG start_ARG square-root start_ARG 2 end_ARG end_ARG | start_ARG â â - â â end_ARG ⊠.
(9)
Then, the unitary transformation U đ U italic_U is formulated:
| e 1 ⊠| e 2 ⊠| e 3 ⊠| e 4 ⊠= U ⢠| d 1 ⊠| d 2 ⊠| d 3 ⊠| d 4 ⊠= ( 1 2 1 2 0 0 i 2 â i 2 0 0 0 0 i 0 0 0 0 1 ) ⢠| d 1 ⊠| d 2 ⊠| d 3 ⊠| d 4 ⊠. ket subscript đ 1 ket subscript đ 2 ket subscript đ 3 ket subscript đ 4 đ ket subscript đ 1 ket subscript đ 2 ket subscript đ 3 ket subscript đ 4 1 2 1 2 0 0 đ 2 đ 2 0 0 0 0 đ 0 0 0 0 1 ket subscript đ 1 ket subscript đ 2 ket subscript đ 3 ket subscript đ 4 \begin{array}[]{l}\ket{e_{1}}\\
\ket{e_{2}}\\
\ket{e_{3}}\\
\ket{e_{4}}\end{array}=U\begin{array}[]{l}\ket{d_{1}}\\
\ket{d_{2}}\\
\ket{d_{3}}\\
\ket{d_{4}}\end{array}=\left(\begin{array}[]{llll}\frac{1}{\sqrt{2}}&\frac{1}{%
\sqrt{2}}&0&0\\
\frac{i}{\sqrt{2}}&\frac{-i}{\sqrt{2}}&0&0\\
0&0&i&0\\
0&0&0&1\end{array}\right)\begin{array}[]{l}\ket{d_{1}}\\
\ket{d_{2}}\\
\ket{d_{3}}\\
\ket{d_{4}}\end{array}. start_ARRAY start_ROW start_CELL | start_ARG italic_e start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG ⊠end_CELL end_ROW start_ROW start_CELL | start_ARG italic_e start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG ⊠end_CELL end_ROW start_ROW start_CELL | start_ARG italic_e start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT end_ARG ⊠end_CELL end_ROW start_ROW start_CELL | start_ARG italic_e start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT end_ARG ⊠end_CELL end_ROW end_ARRAY = italic_U start_ARRAY start_ROW start_CELL | start_ARG italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG ⊠end_CELL end_ROW start_ROW start_CELL | start_ARG italic_d start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG ⊠end_CELL end_ROW start_ROW start_CELL | start_ARG italic_d start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT end_ARG ⊠end_CELL end_ROW start_ROW start_CELL | start_ARG italic_d start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT end_ARG ⊠end_CELL end_ROW end_ARRAY = ( start_ARRAY start_ROW start_CELL divide start_ARG 1 end_ARG start_ARG square-root start_ARG 2 end_ARG end_ARG end_CELL start_CELL divide start_ARG 1 end_ARG start_ARG square-root start_ARG 2 end_ARG end_ARG end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL divide start_ARG italic_i end_ARG start_ARG square-root start_ARG 2 end_ARG end_ARG end_CELL start_CELL divide start_ARG - italic_i end_ARG start_ARG square-root start_ARG 2 end_ARG end_ARG end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL italic_i end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 1 end_CELL end_ROW end_ARRAY ) start_ARRAY start_ROW start_CELL | start_ARG italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG ⊠end_CELL end_ROW start_ROW start_CELL | start_ARG italic_d start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG ⊠end_CELL end_ROW start_ROW start_CELL | start_ARG italic_d start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT end_ARG ⊠end_CELL end_ROW start_ROW start_CELL | start_ARG italic_d start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT end_ARG ⊠end_CELL end_ROW end_ARRAY .
(10)
In the following sections, the total-spin basis is also utilized:
| Ψ ⢠( 1 , 2 ) ⊠= â i F i ⢠( đ 1 , đ 2 ) ⢠| d i ⊠ket Ψ 1 2 subscript đ subscript đš đ subscript đ 1 subscript đ 2 ket subscript đ đ \ket{\Psi(1,2)}=\sum_{i}F_{i}(\bm{r}_{1},\bm{r}_{2})\ket{d_{i}} | start_ARG roman_Ψ ( 1 , 2 ) end_ARG ⊠= â start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_F start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( bold_italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) | start_ARG italic_d start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG ⊠.
Thus, G i ⢠( đ 1 , đ 2 ) = â k ( U â ) i ⢠k ⢠F k ⢠( đ 1 , đ 2 ) subscript đş đ subscript đ 1 subscript đ 2 subscript đ subscript superscript đ â đ đ subscript đš đ subscript đ 1 subscript đ 2 G_{i}(\bm{r}_{1},\bm{r}_{2})=\sum_{k}\left(U^{\dagger}\right)_{ik}F_{k}(\bm{r}%
_{1},\bm{r}_{2}) italic_G start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( bold_italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) = â start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ( italic_U start_POSTSUPERSCRIPT â end_POSTSUPERSCRIPT ) start_POSTSUBSCRIPT italic_i italic_k end_POSTSUBSCRIPT italic_F start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ( bold_italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) to convert to Bell basis.
Notice that, for computing the expectation values in Eq. (6 ),
these coordinate parts can be simply integrated.
I employ the three-body model, which has been developed and utilized in Refs. Suzuki and Ikeda (1988 ); Bertsch and Esbensen (1991 ); Esbensen et al. (1997 ); Hagino and Sagawa (2005 ); Oishi et al. (2014 , 2017 ) .
The system contains an alpha particle as the rigid core with mass m C subscript đ đś m_{C} italic_m start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT and two valence protons.
The two valence protons feel the spherical mean field V ⢠( đ ) đ đ V(\bm{r}) italic_V ( bold_italic_r ) generated by the alpha core.
Thus, the three-body Hamiltonian reads
H ^ 3 ⢠B = h ^ ⢠( đ 1 ) + h ^ ⢠( đ 2 ) + v p ⢠p ⢠( đ 1 , đ 2 ) + đ 1 â
đ 2 m C , subscript ^ đť 3 đľ ^ â subscript đ 1 ^ â subscript đ 2 subscript đŁ đ đ subscript đ 1 subscript đ 2 â
subscript đ 1 subscript đ 2 subscript đ đś \displaystyle\hat{H}_{3B}=\hat{h}(\bm{r}_{1})+\hat{h}(\bm{r}_{2})+v_{pp}(\bm{r%
}_{1},\bm{r}_{2})+\frac{\bm{p}_{1}\cdot\bm{p}_{2}}{m_{C}}, over^ start_ARG italic_H end_ARG start_POSTSUBSCRIPT 3 italic_B end_POSTSUBSCRIPT = over^ start_ARG italic_h end_ARG ( bold_italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) + over^ start_ARG italic_h end_ARG ( bold_italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) + italic_v start_POSTSUBSCRIPT italic_p italic_p end_POSTSUBSCRIPT ( bold_italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) + divide start_ARG bold_italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT â
bold_italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT end_ARG ,
(11)
h ^ ⢠( r i ) = â â 2 2 ⢠Ο ⢠d 2 d ⢠r i 2 + â 2 2 ⢠Ο ⢠l ⢠( l + 1 ) r i 2 + V ⢠( r i ) , ^ â subscript đ đ superscript Planck-constant-over-2-pi 2 2 đ superscript đ 2 đ superscript subscript đ đ 2 superscript Planck-constant-over-2-pi 2 2 đ đ đ 1 superscript subscript đ đ 2 đ subscript đ đ \displaystyle\hat{h}(r_{i})=-\frac{\hbar^{2}}{2\mu}\frac{d^{2}}{dr_{i}^{2}}+%
\frac{\hbar^{2}}{2\mu}\frac{l(l+1)}{r_{i}^{2}}+V(r_{i}), over^ start_ARG italic_h end_ARG ( italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) = - divide start_ARG roman_â start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_Îź end_ARG divide start_ARG italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_d italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + divide start_ARG roman_â start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_Îź end_ARG divide start_ARG italic_l ( italic_l + 1 ) end_ARG start_ARG italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + italic_V ( italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) ,
(12)
where Îź = m p ⢠m C / ( m p + m C ) đ subscript đ đ subscript đ đś subscript đ đ subscript đ đś \mu=m_{p}m_{C}/(m_{p}+m_{C}) italic_Îź = italic_m start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT italic_m start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT / ( italic_m start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT + italic_m start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT ) and m p = 938.272 subscript đ đ 938.272 m_{p}=938.272 italic_m start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT = 938.272 MeV/ c 2 absent superscript đ 2 /c^{2} / italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT for protons.
Here h ^ ⢠( đ k ) ^ â subscript đ đ \hat{h}(\bm{r}_{k}) over^ start_ARG italic_h end_ARG ( bold_italic_r start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) is of the k đ k italic_k th proton-alpha subsystem.
For its interaction V ⢠( r i ) đ subscript đ đ V(r_{i}) italic_V ( italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) , the same Woods-Saxon and Coulomb potentials in the previous work Oishi et al. (2017 ) are employed.
The proton-proton interaction reads
v p ⢠p ⢠( đ 1 , đ 2 ) = v n ⢠u ⢠c ⢠l ⢠( d ) + e 2 4 â˘ Ď â˘ Ďľ 0 ⢠1 d + v a ⢠d ⢠d ⢠( đ 1 , đ 2 ) , subscript đŁ đ đ subscript đ 1 subscript đ 2 subscript đŁ đ đ˘ đ đ đ superscript đ 2 4 đ subscript italic-Ďľ 0 1 đ subscript đŁ đ đ đ subscript đ 1 subscript đ 2 v_{pp}(\bm{r}_{1},\bm{r}_{2})=v_{nucl}(d)+\frac{e^{2}}{4\pi\epsilon_{0}}\frac{%
1}{d}+v_{add}(\bm{r}_{1},\bm{r}_{2}), italic_v start_POSTSUBSCRIPT italic_p italic_p end_POSTSUBSCRIPT ( bold_italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) = italic_v start_POSTSUBSCRIPT italic_n italic_u italic_c italic_l end_POSTSUBSCRIPT ( italic_d ) + divide start_ARG italic_e start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 4 italic_Ď italic_Ďľ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG divide start_ARG 1 end_ARG start_ARG italic_d end_ARG + italic_v start_POSTSUBSCRIPT italic_a italic_d italic_d end_POSTSUBSCRIPT ( bold_italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) ,
(13)
where d = | đ 2 â đ 1 | đ subscript đ 2 subscript đ 1 d=\left|\bm{r}_{2}-\bm{r}_{1}\right| italic_d = | bold_italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - bold_italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT | .
The nuclear-force term is described by the spin-dependent Gaussian potential Thompson et al. (1977 ) :
v n ⢠u ⢠c ⢠l ⢠( d ) subscript đŁ đ đ˘ đ đ đ \displaystyle v_{nucl}(d) italic_v start_POSTSUBSCRIPT italic_n italic_u italic_c italic_l end_POSTSUBSCRIPT ( italic_d )
= \displaystyle= =
[ V R ⢠e â a R ⢠d 2 + V S ⢠e â a S ⢠d 2 ] ⢠P ^ S = 0 delimited-[] subscript đ đ
superscript đ subscript đ đ
superscript đ 2 subscript đ đ superscript đ subscript đ đ superscript đ 2 subscript ^ đ đ 0 \displaystyle\left[V_{R}e^{-a_{R}d^{2}}+V_{S}e^{-a_{S}d^{2}}\right]\hat{P}_{S=0} [ italic_V start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT - italic_a start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT + italic_V start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT - italic_a start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT ] over^ start_ARG italic_P end_ARG start_POSTSUBSCRIPT italic_S = 0 end_POSTSUBSCRIPT
(14)
+ [ V R ⢠e â a R ⢠d 2 + V T ⢠e â a T ⢠d 2 ] ⢠P ^ S = 1 . delimited-[] subscript đ đ
superscript đ subscript đ đ
superscript đ 2 subscript đ đ superscript đ subscript đ đ superscript đ 2 subscript ^ đ đ 1 \displaystyle+\left[V_{R}e^{-a_{R}d^{2}}+V_{T}e^{-a_{T}d^{2}}\right]\hat{P}_{S%
=1}. + [ italic_V start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT - italic_a start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT + italic_V start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT - italic_a start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT ] over^ start_ARG italic_P end_ARG start_POSTSUBSCRIPT italic_S = 1 end_POSTSUBSCRIPT .
The operator P ^ S = 0 subscript ^ đ đ 0 \hat{P}_{S=0} over^ start_ARG italic_P end_ARG start_POSTSUBSCRIPT italic_S = 0 end_POSTSUBSCRIPT (P ^ S = 1 subscript ^ đ đ 1 \hat{P}_{S=1} over^ start_ARG italic_P end_ARG start_POSTSUBSCRIPT italic_S = 1 end_POSTSUBSCRIPT ) indicates the projection into the
spin-singlet (spin-triplet) channel of the proton-proton subsystem.
Parameters are given as
V R = 200 subscript đ đ
200 V_{R}=200 italic_V start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT = 200 MeV,
V S = â 91.85 subscript đ đ 91.85 V_{S}=-91.85 italic_V start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT = - 91.85 MeV,
V T = â 178 subscript đ đ 178 V_{T}=-178 italic_V start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT = - 178 MeV,
a R = 1 , 487 subscript đ đ
1 487
a_{R}=1,487 italic_a start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT = 1 , 487 fm-2 ,
a S = 0.465 subscript đ đ 0.465 a_{S}=0.465 italic_a start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT = 0.465 fm-2 , and
a T = 0.639 subscript đ đ 0.639 a_{T}=0.639 italic_a start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT = 0.639 fm-2 .
These parameters correctly reproduce the experimental vacuum-scattering length of two protons Thompson et al. (1977 ) .
In addition, the surface-dependent term v a ⢠d ⢠d subscript đŁ đ đ đ v_{add} italic_v start_POSTSUBSCRIPT italic_a italic_d italic_d end_POSTSUBSCRIPT is employed:
v a ⢠d ⢠d ⢠( đ 1 , đ 2 ) = w 0 ⢠e â ( R â R 0 ) 2 / B 0 2 ⢠δ ⢠( đ 1 â đ 2 ) , subscript đŁ đ đ đ subscript đ 1 subscript đ 2 subscript đ¤ 0 superscript đ superscript đ
subscript đ
0 2 superscript subscript đľ 0 2 đż subscript đ 1 subscript đ 2 v_{add}(\bm{r}_{1},\bm{r}_{2})=w_{0}e^{-(R-R_{0})^{2}/B_{0}^{2}}\delta(\bm{r}_%
{1}-\bm{r}_{2}), italic_v start_POSTSUBSCRIPT italic_a italic_d italic_d end_POSTSUBSCRIPT ( bold_italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) = italic_w start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT - ( italic_R - italic_R start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT italic_δ ( bold_italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - bold_italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) ,
(15)
where R = | ( đ 1 + đ 2 ) / 2 | đ
subscript đ 1 subscript đ 2 2 R=\left|(\bm{r}_{1}+\bm{r}_{2})/2\right| italic_R = | ( bold_italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + bold_italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) / 2 | ,
R 0 = 1.68 subscript đ
0 1.68 R_{0}=1.68 italic_R start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1.68 fm,
B 0 = 0.6 ⢠R 0 subscript đľ 0 0.6 subscript đ
0 B_{0}=0.6R_{0} italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.6 italic_R start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , and
w 0 = â 470 subscript đ¤ 0 470 w_{0}=-470 italic_w start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = - 470 MeVâ
â
\cdot â
fm3 .
This additional term is necessary to reproduce the experimental energy:
the 2 ⢠p 2 đ 2p 2 italic_p  energy and width are obtained as E 2 ⢠p = 1.356 subscript đ¸ 2 đ 1.356 E_{2p}=1.356 italic_E start_POSTSUBSCRIPT 2 italic_p end_POSTSUBSCRIPT = 1.356 MeV and Î 2 ⢠p = 0.055 subscript Î 2 đ 0.055 \Gamma_{2p}=0.055 roman_Î start_POSTSUBSCRIPT 2 italic_p end_POSTSUBSCRIPT = 0.055 MeV, respectively,
whereas the experimental data read E 2 ⢠p = 1.372 ⢠( 5 ) subscript đ¸ 2 đ 1.372 5 E_{2p}=1.372(5) italic_E start_POSTSUBSCRIPT 2 italic_p end_POSTSUBSCRIPT = 1.372 ( 5 ) MeV and Î 2 ⢠p = 0.092 ⢠( 6 ) subscript Î 2 đ 0.092 6 \Gamma_{2p}=0.092(6) roman_Î start_POSTSUBSCRIPT 2 italic_p end_POSTSUBSCRIPT = 0.092 ( 6 ) MeV Ajzenberg-Selove (1988 , 1991 ) .
Notice that this additional potential vanishes when one of the three particles is infinitely separated.
Thus, the vacuum properties of two-body subsystems can be conserved in the time-development calculations.
I especially notify that the emitted two protons should be unbound in vacuum.
FIG. 2:
Time-dependent decaying probability of 2 ⢠p 2 đ 2p 2 italic_p  emission from 6 Be.
Its spin-singlet (S 12 = 0 subscript đ 12 0 S_{12}=0 italic_S start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT = 0 ) component is also plotted.
The factor f p ⢠p = 1 subscript đ đ đ 1 f_{pp}=1 italic_f start_POSTSUBSCRIPT italic_p italic_p end_POSTSUBSCRIPT = 1 indicates that the proton-proton interaction is adjusted to reproduce the experimental 2 ⢠p 2 đ 2p 2 italic_p -energy.
Results and Discussions - The 2 ⢠p 2 đ 2p 2 italic_p -emitting process is simulated with the time-dependent method Oishi et al. (2014 , 2017 ) :
| Ψ ⢠( t ) ⊠= exp ⥠[ â i ⢠t ⢠H ^ 3 ⢠B â ] ⢠| Ψ ⢠( 0 ) ⊠, ket Ψ đĄ đ đĄ subscript ^ đť 3 đľ Planck-constant-over-2-pi ket Ψ 0 \ket{\Psi(t)}=\exp\left[-it\frac{\hat{H}_{3B}}{\hbar}\right]\ket{\Psi(0)}, | start_ARG roman_Ψ ( italic_t ) end_ARG ⊠= roman_exp [ - italic_i italic_t divide start_ARG over^ start_ARG italic_H end_ARG start_POSTSUBSCRIPT 3 italic_B end_POSTSUBSCRIPT end_ARG start_ARG roman_â end_ARG ] | start_ARG roman_Ψ ( 0 ) end_ARG ⊠,
(16)
where the initial state is solved as the confined 2 ⢠p 2 đ 2p 2 italic_p  state inside the Coulomb barrier Oishi et al. (2014 ) .
The decaying state | Ψ d ⢠( t ) ⊠ket subscript Ψ đ đĄ \ket{\Psi_{d}(t)} | start_ARG roman_Ψ start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT ( italic_t ) end_ARG ⊠, which describes the emitted component outside the barrier,
is determined as
| Ψ d ⢠( t ) ⊠= | Ψ ⢠( t ) ⊠â β ⢠( t ) ⢠| Ψ ⢠( 0 ) ⊠, ket subscript Ψ đ đĄ ket Ψ đĄ đ˝ đĄ ket Ψ 0 \ket{\Psi_{d}(t)}=\ket{\Psi(t)}-\beta(t)\ket{\Psi(0)}, | start_ARG roman_Ψ start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT ( italic_t ) end_ARG ⊠= | start_ARG roman_Ψ ( italic_t ) end_ARG ⊠- italic_β ( italic_t ) | start_ARG roman_Ψ ( 0 ) end_ARG ⊠,
(17)
where β ⢠( t ) đ˝ đĄ \beta(t) italic_β ( italic_t ) is the survival coefficient,
β ⢠( t ) = ⨠Ψ ⢠( 0 ) | Ψ ⢠( t ) ⊠đ˝ đĄ inner-product Ψ 0 Ψ đĄ \beta(t)=\Braket{\Psi(0)}{\Psi(t)} italic_β ( italic_t ) = ⨠start_ARG roman_Ψ ( 0 ) end_ARG | start_ARG roman_Ψ ( italic_t ) end_ARG ⊠.
In FIG 2 , the time-dependent decaying probability is displayed:
P d ⢠e ⢠c ⢠a ⢠y ⢠( t ) = ⨠Ψ d ⢠( t ) ⣠Ψ d ⢠( t ) ⊠= 1 â | β ⢠( t ) | 2 . subscript đ đ đ đ đ đŚ đĄ expectation conditional subscript Ψ đ đĄ subscript Ψ đ đĄ 1 superscript đ˝ đĄ 2 P_{decay}(t)=\Braket{\Psi_{d}(t)\mid\Psi_{d}(t)}=1-\left|\beta(t)\right|^{2}. italic_P start_POSTSUBSCRIPT italic_d italic_e italic_c italic_a italic_y end_POSTSUBSCRIPT ( italic_t ) = ⨠start_ARG roman_Ψ start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT ( italic_t ) ⣠roman_Ψ start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT ( italic_t ) end_ARG ⊠= 1 - | italic_β ( italic_t ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT .
(18)
One clearly finds that the decaying probability increases in time development.
The spin-singlet state, | d 4 ⊠= | e 4 ⊠= | â ⣠â ⣠â ⣠â ⣠â ⊠/ 2 ket subscript đ 4 ket subscript đ 4 ket â â â â
2 \ket{d_{4}}=\ket{e_{4}}=\ket{\uparrow\downarrow-\downarrow\uparrow}/\sqrt{2} | start_ARG italic_d start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT end_ARG ⊠= | start_ARG italic_e start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT end_ARG ⊠= | start_ARG â â - â â end_ARG ⊠/ square-root start_ARG 2 end_ARG , is always dominant.
I confirmed that the survival probability,
P s ⢠u ⢠r ⢠v ⢠( t ) ⥠⨠Ψ ⢠( 0 ) ⣠Ψ ⢠( t ) ⊠2 = | β ⢠( t ) | 2 = 1 â P d ⢠e ⢠c ⢠a ⢠y ⢠( t ) , subscript đ đ đ˘ đ đŁ đĄ superscript expectation conditional Ψ 0 Ψ đĄ 2 superscript đ˝ đĄ 2 1 subscript đ đ đ đ đ đŚ đĄ P_{surv}(t)\equiv\Braket{\Psi(0)\mid\Psi(t)}^{2}=\left|\beta(t)\right|^{2}=1-P%
_{decay}(t), italic_P start_POSTSUBSCRIPT italic_s italic_u italic_r italic_v end_POSTSUBSCRIPT ( italic_t ) ⥠⨠start_ARG roman_Ψ ( 0 ) ⣠roman_Ψ ( italic_t ) end_ARG ⊠start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = | italic_β ( italic_t ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = 1 - italic_P start_POSTSUBSCRIPT italic_d italic_e italic_c italic_a italic_y end_POSTSUBSCRIPT ( italic_t ) ,
(19)
is well approximated by the exponential damping.
Namely, P s ⢠u ⢠r ⢠v ⢠( t ) â
e â t / Ď subscript đ đ đ˘ đ đŁ đĄ superscript đ đĄ đ P_{surv}(t)\cong e^{-t/\tau} italic_P start_POSTSUBSCRIPT italic_s italic_u italic_r italic_v end_POSTSUBSCRIPT ( italic_t ) â
italic_e start_POSTSUPERSCRIPT - italic_t / italic_Ď end_POSTSUPERSCRIPT , where the width can be evaluated as Î 2 ⢠p = â / Ď â
0.055 subscript Î 2 đ Planck-constant-over-2-pi đ 0.055 \Gamma_{2p}=\hbar/\tau\cong 0.055 roman_Î start_POSTSUBSCRIPT 2 italic_p end_POSTSUBSCRIPT = roman_â / italic_Ď â
0.055 MeV by numerical fitting.
Also, by observing the decaying density distribution, Ď d ⢠( t , đ 1 , đ 2 ) = | Ψ d ⢠( t , đ 1 , đ 2 ) | 2 subscript đ đ đĄ subscript đ 1 subscript đ 2 superscript subscript Ψ đ đĄ subscript đ 1 subscript đ 2 2 \rho_{d}(t,\bm{r}_{1},\bm{r}_{2})=\left|\Psi_{d}(t,\bm{r}_{1},\bm{r}_{2})%
\right|^{2} italic_Ď start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT ( italic_t , bold_italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) = | roman_Ψ start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT ( italic_t , bold_italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ,
I confirmed that the present 2 ⢠p 2 đ 2p 2 italic_p  emission can be interpreted
as the diproton-correlating emission Oishi et al. (2014 , 2017 ) .
FIG. 3:
CHSH indicator of the time-dependent decaying state | Ψ d ⢠( t ) ⊠ket subscript Ψ đ đĄ \ket{\Psi_{d}(t)} | start_ARG roman_Ψ start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT ( italic_t ) end_ARG ⊠of the 2 ⢠p 2 đ 2p 2 italic_p  emission of 6 Be ✠ι + p + p âś absent đź đ đ \longrightarrow\alpha+p+p âś italic_Îą + italic_p + italic_p .
Here θ = Ď / 4 đ đ 4 \theta=\pi/4 italic_θ = italic_Ď / 4 .
The factor f p ⢠p subscript đ đ đ f_{pp} italic_f start_POSTSUBSCRIPT italic_p italic_p end_POSTSUBSCRIPT indicates the strength of the proton-proton interaction.
In Fig. 3 , the CHSH indicator is evaluated for the decaying state by using the expansion,
| Ψ d ⢠( t , đ 1 , đ 2 ) ⊠= â i = 1 4 G i ( t , đ 1 , đ 2 ) ) | e i ⊠. \ket{\Psi_{d}(t,\bm{r}_{1},\bm{r}_{2})}=\sum_{i=1}^{4}G_{i}(t,\bm{r}_{1},\bm{r%
}_{2}))\ket{e_{i}}. | start_ARG roman_Ψ start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT ( italic_t , bold_italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) end_ARG ⊠= â start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_G start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_t , bold_italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) ) | start_ARG italic_e start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG ⊠.
(20)
One clearly finds that the CHSH indicator becomes larger than 2 2 2 2 , i.e. the LHV-theory limit, during the time evolution.
For c ⢠t ⼠1500 đ đĄ 1500 ct\geq 1500 italic_c italic_t ⼠1500 fm, I obtain S â
2.72 đ 2.72 S\cong 2.72 italic_S â
2.72 with the default setting of v p ⢠p subscript đŁ đ đ v_{pp} italic_v start_POSTSUBSCRIPT italic_p italic_p end_POSTSUBSCRIPT (f p ⢠p = 1.0 subscript đ đ đ 1.0 f_{pp}=1.0 italic_f start_POSTSUBSCRIPT italic_p italic_p end_POSTSUBSCRIPT = 1.0 ) for reproducing E 2 ⢠p = 1.356 subscript đ¸ 2 đ 1.356 E_{2p}=1.356 italic_E start_POSTSUBSCRIPT 2 italic_p end_POSTSUBSCRIPT = 1.356 MeV of 6 Be Ajzenberg-Selove (1988 , 1991 ) .
This coupled-spin entanglement is of course attributable to the dominant spin-singlet component in FIG 2 .
Notice that the present time evolution in Eq. (16 ) also contains the effect of core-proton interactions V ⢠( r i ) đ subscript đ đ V(r_{i}) italic_V ( italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) .
These interactions do not harm the 2 ⢠p 2 đ 2p 2 italic_p  entanglement.
For deeper investigation, the sensitivity of CHSH indicator to the proton-proton interaction is studied.
For this purpose, an effective tuning factor is employed:
v p ⢠p ⢠( đ 1 , đ 2 ) âś f p ⢠p ⢠v p ⢠p ⢠( đ 1 , đ 2 ) . âś subscript đŁ đ đ subscript đ 1 subscript đ 2 subscript đ đ đ subscript đŁ đ đ subscript đ 1 subscript đ 2 v_{pp}(\bm{r}_{1},\bm{r}_{2})\longrightarrow f_{pp}v_{pp}(\bm{r}_{1},\bm{r}_{2%
}). italic_v start_POSTSUBSCRIPT italic_p italic_p end_POSTSUBSCRIPT ( bold_italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) âś italic_f start_POSTSUBSCRIPT italic_p italic_p end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT italic_p italic_p end_POSTSUBSCRIPT ( bold_italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) .
(21)
With f p ⢠p = 1 subscript đ đ đ 1 f_{pp}=1 italic_f start_POSTSUBSCRIPT italic_p italic_p end_POSTSUBSCRIPT = 1 , the v p ⢠p subscript đŁ đ đ v_{pp} italic_v start_POSTSUBSCRIPT italic_p italic_p end_POSTSUBSCRIPT approximates the experimental energy, E 2 ⢠p = 1.356 subscript đ¸ 2 đ 1.356 E_{2p}=1.356 italic_E start_POSTSUBSCRIPT 2 italic_p end_POSTSUBSCRIPT = 1.356 MeV Ajzenberg-Selove (1988 , 1991 ) .
From FIG. 3 , one can read that this CHSH indicator becomes smaller when the interaction is weakened.
In the present case of 6 Be, it goes below the LHV-theory limit, S < 2 đ 2 S<2 italic_S < 2 ,
with f p ⢠p ⤠0.2 subscript đ đ đ 0.2 f_{pp}\leq 0.2 italic_f start_POSTSUBSCRIPT italic_p italic_p end_POSTSUBSCRIPT ⤠0.2 .
Consequently, the present entanglement is a product of the proton-proton interaction.
With the v p ⢠p subscript đŁ đ đ v_{pp} italic_v start_POSTSUBSCRIPT italic_p italic_p end_POSTSUBSCRIPT weakened, the 2 ⢠p 2 đ 2p 2 italic_p  energy and width also change:
E 2 ⢠p = 2.864 subscript đ¸ 2 đ 2.864 E_{2p}=2.864 italic_E start_POSTSUBSCRIPT 2 italic_p end_POSTSUBSCRIPT = 2.864 MeV and Î 2 ⢠p = 0.8579 subscript Î 2 đ 0.8579 \Gamma_{2p}=0.8579 roman_Î start_POSTSUBSCRIPT 2 italic_p end_POSTSUBSCRIPT = 0.8579 MeV with f p ⢠p = 0.5 subscript đ đ đ 0.5 f_{pp}=0.5 italic_f start_POSTSUBSCRIPT italic_p italic_p end_POSTSUBSCRIPT = 0.5 ;
E 2 ⢠p = 2.925 subscript đ¸ 2 đ 2.925 E_{2p}=2.925 italic_E start_POSTSUBSCRIPT 2 italic_p end_POSTSUBSCRIPT = 2.925 MeV and Î 2 ⢠p â
2.9 subscript Î 2 đ 2.9 \Gamma_{2p}\cong 2.9 roman_Î start_POSTSUBSCRIPT 2 italic_p end_POSTSUBSCRIPT â
2.9 MeV with f p ⢠p = 0.2 subscript đ đ đ 0.2 f_{pp}=0.2 italic_f start_POSTSUBSCRIPT italic_p italic_p end_POSTSUBSCRIPT = 0.2 .
Thus, the short-lived (broad-width) 2 ⢠p 2 đ 2p 2 italic_p  state has the weaker spin entanglement.
Summary - In this paper, the CHSH indicator for the 2 ⢠p 2 đ 2p 2 italic_p  emission from 6 Be
is evaluated in order to measure its coupled-spin entanglement.
As a product of proton-proton interaction, S â
2.72 đ 2.72 S\cong 2.72 italic_S â
2.72 is obtained in the time-dependent decaying state: an entanglement beyond the limit of LHV theory is suggested.
This entanglement is not harmed by the core-proton interactions.
The CHSH measurement can be a probe into the effective nuclear interaction inside finite systems.
This work is limited to the lightest 6 Be nucleus.
Whether the spin entanglement exists commonly in other 2 ⢠p 2 đ 2p 2 italic_p  emitters or not is an open question.
The dominance of spin-singlet Bell state is not trivial in other systems.
For evaluating the quantum entanglement, the von-Neumann entropy has been considered as an essential
quantity Kanada-Enâyo (2015 ); Bulgac (2023 ); Pazy (2023 ); Gu et al. (2023 ); Hengstenberg et al. (2023 ); PĂŠrez-Obiol et al. (2023 ); Kou et al. (2024 ); Kirchner et al. (2024 ) ,
whereas this work focuses only on the CHSH indicator.
Evaluation and discussion of the von-Neumann entropy for the time-dependent state is in progress.
Whether the entanglement is apparent or not may depend on the observable quantity of interest, e.g. the coupled spin, momenta, and energy distribution.
These topics will be addressed in forthcoming studies.
In the experimental side, mass production of 2 ⢠p 2 đ 2p 2 italic_p  emitters, including the 6 Be, for a sufficient statistics is still challenging.
For measuring the two protons by independent detectors after decay, an advanced design of experiment will be necessary.
Closing all the loopholes of
detection Rowe et al. (2001 ); Giustina et al. (2015 ); Shalm et al. (2015 ) ,
locality Aspect et al. (1982 ); Weihs et al. (1998 ); Scheidl et al. (2010 ) , and
memory Barrett et al. (2002 )
should require another lot of efforts.
On the other side, the experimental survey of 2 ⢠p 2 đ 2p 2 italic_p  emitters is rapidly in progress, e.g. that in RIKEN RIBF.
Confirmation of spin entanglement in 2 ⢠p 2 đ 2p 2 italic_p  radioactivity can be one landmark in the nuclear physics as well as quantum many-body science.
Acknowledgments - I especially thank Yutaka Shikano, Masaki Sasano, Tomoya Naito,
Tokuro Fukui, and Masaaki Kimura for fruitful discussions.
Numerical calculations are supported by the cooperative project of supercomputer Yukawa-21 in Yukawa Institute for Theoretical Physics, Kyoto University.
I appreciate the Multi-disciplinary Cooperative Research Program (MCRP) 2023-2024 by Center for Computational Sciences, University of Tsukuba (project ID wo23i034), allocating computational resources of supercomputer Wisteria/BDEC-01 (Odyssey) in Information Technology Center, University of Tokyo.
References
Clauser et al. (1969)
J. F. Clauser, M. A. Horne,
A. Shimony,  and R. A. Holt, Phys. Rev. Lett. 23 , 880 (1969) .
Bell (1964)
J. Bell, Physics 1 , 195
(1964).
Bell (2004)
J. S. Bell, Speakable and Unspeakable
in Quantum Mechanics , 2nd ed., Collected Papers on
Quantum Philosophy (Cambridge University Press, Cambridge, UK, 2004).
Clauser and Shimony (1978)
J. F. Clauser and A. Shimony, Reports on Progress in Physics 41 , 1881 (1978) .
Aspect et al. (1982)
A. Aspect, P. Grangier,  and G. Roger, Phys. Rev. Lett. 49 , 91 (1982) .
Weihs et al. (1998)
G. Weihs, T. Jennewein,
C. Simon, H. Weinfurter,  and A. Zeilinger, Phys. Rev. Lett. 81 , 5039 (1998) .
Pan et al. (2000)
J.-W. Pan, D. Bouwmeester,
M. Daniell, H. Weinfurter,  and A. Zeilinger, Nature 403 , 515 (2000) .
Rowe et al. (2001)
M. A. Rowe, D. Kielpinski,
V. Meyer, C. A. Sackett, W. M. Itano, C. Monroe,  and D. J. Wineland, Nature 409 , 791 (2001) .
Scheidl et al. (2010)
T. Scheidl, R. Ursin,
J. Kofler, S. Ramelow, X.-S. Ma, T. Herbst, L. Ratschbacher, A. Fedrizzi, N. K. Langford, T. Jennewein,  and A. Zeilinger, Proceedings of the National Academy of
Sciences 107 , 19708
(2010) , https://www.pnas.org/doi/pdf/10.1073/pnas.1002780107 .
Giustina et al. (2015)
M. Giustina, M. A. M. Versteegh, S. Wengerowsky, J. Handsteiner, A. Hochrainer, K. Phelan,
F. Steinlechner, J. Kofler, J.-A. Larsson, C. Abellån, W. Amaya, V. Pruneri, M. W. Mitchell, J. Beyer, T. Gerrits,
A. E. Lita, L. K. Shalm, S. W. Nam, T. Scheidl, R. Ursin, B. Wittmann,  and A. Zeilinger, Phys. Rev. Lett. 115 , 250401 (2015) .
Shalm et al. (2015)
L. K. Shalm, E. Meyer-Scott,
B. G. Christensen,
P. Bierhorst, M. A. Wayne, M. J. Stevens, T. Gerrits, S. Glancy, D. R. Hamel, M. S. Allman, K. J. Coakley,
S. D. Dyer, C. Hodge, A. E. Lita, V. B. Verma, C. Lambrocco, E. Tortorici, A. L. Migdall, Y. Zhang, D. R. Kumor, W. H. Farr, F. Marsili, M. D. Shaw,
J. A. Stern, C. Abellån, W. Amaya, V. Pruneri, T. Jennewein, M. W. Mitchell, P. G. Kwiat, J. C. Bienfang, R. P. Mirin,
E. Knill,  and S. W. Nam, Phys. Rev. Lett. 115 , 250402 (2015) .
Vasilyev et al. (2017)
D. Vasilyev, F. O. Schumann, F. Giebels,
H. Gollisch, J. Kirschner,  and R. Feder, Phys.
Rev. B 95 , 115134
(2017) .
Storz et al. (2023)
S. Storz, J. Schär,
A. Kulikov, P. Magnard, P. Kurpiers, J. Lßtolf, T. Walter, A. Copetudo, K. Reuer, A. Akin, J.-C. Besse, M. Gabureac, G. J. Norris, A. Rosario,
F. Martin, J. Martinez, W. Amaya, M. W. Mitchell, C. Abellan, J.-D. Bancal, N. Sangouard, B. Royer,
A. Blais,  and A. Wallraff, Nature 617 , 265
(2023) .
Barrett et al. (2002)
J. Barrett, D. Collins,
L. Hardy, A. Kent,  and S. Popescu, Phys.
Rev. A 66 , 042111
(2002) .
Lamehi-Rachti and Mittig (1976)
M. Lamehi-Rachti and W. Mittig, Phys. Rev. D 14 , 2543 (1976) .
Sakai et al. (2006)
H. Sakai, T. Saito,
T. Ikeda, K. Itoh, T. Kawabata, H. Kuboki, Y. Maeda, N. Matsui, C. Rangacharyulu, M. Sasano, Y. Satou, K. Sekiguchi, K. Suda, A. Tamii, T. Uesaka,  and K. Yako, Phys. Rev. Lett. 97 , 150405 (2006) .
Kanada-Enâyo (2015)
Y. Kanada-Enâyo, Progress of Theoretical and Experimental
Physics 2015 , 043D04
(2015) .
Bulgac (2023)
A. Bulgac, Phys. Rev. C 107 , L061602 (2023) .
Pazy (2023)
E. Pazy, Phys. Rev. C 107 , 054308 (2023) .
Gu et al. (2023)
C. Gu, Z. H. Sun,
G. Hagen,  and T. Papenbrock, Phys. Rev. C 108 , 054309 (2023) .
Hengstenberg et al. (2023)
S. M. Hengstenberg, C. E. P. Robin,  and M. J. Savage, The European Physical Journal A 59 , 231 (2023) .
PÊrez-Obiol et al. (2023)
A. PÊrez-Obiol, S. Masot-Llima, A. M. Romero, J. MenÊndez, A. Rios,
A. Garcia-Såez,  and B. Juliå-Diaz, The European Physical Journal A 59 , 240 (2023) .
Kou et al. (2024)
W. Kou, J. Chen,  and X. Chen, Physics Letters B 849 , 138453 (2024) .
Kirchner et al. (2024)
T. Kirchner, W. Elkamhawy,
 and H.-W. Hammer, Few-Body Systems 65 , 29 (2024) .
Robin et al. (2021)
C. Robin, M. J. Savage,  and N. Pillet, Phys. Rev. C 103 , 034325 (2021) .
Kruppa et al. (2021)
A. T. Kruppa, J. Kovåcs,
P. Salamon,  and O. Legeza, Journal of Physics G: Nuclear and Particle Physics 48 , 025107 (2021) .
Miller (2023a)
G. A. Miller, Phys. Rev. C 108 , L041601 (2023a) .
Miller (2023b)
G. A. Miller, Phys. Rev. C 108 , L031002 (2023b) .
Tichai et al. (2023)
A. Tichai, S. Knecht,
A. Kruppa, O. Legeza, C. Moca, A. Schwenk, M. Werner,  and G. Zarand, Physics Letters B 845 , 138139 (2023) .
Sun et al. (2023)
Z. H. Sun, G. Hagen,  and T. Papenbrock, Phys. Rev. C 108 , 014307 (2023) .
Bai (2024)
D. Bai, Phys. Rev. C 109 , 034001 (2024) .
Johnson and Gorton (2023)
C. W. Johnson and O. C. Gorton, Journal of Physics G: Nuclear and Particle Physics 50 , 045110 (2023) .
Bertulani et al. (2008)
C. Bertulani, M. Hussein,
 and G. Verde, Physics Letters B 666 , 86 (2008) .
Grigorenko (2009)
L. V. Grigorenko, Physics of Particles and Nuclei 40 , 674 (2009) .
Pfßtzner et al. (2012)
M. Pfßtzner, M. Karny,
L. V. Grigorenko,  and K. Riisager, Rev. Mod. Phys. 84 , 567 (2012) .
Blank and Ploszajczak (2008)
B. Blank and M. Ploszajczak, Reports on Progress in Physics 71 , 046301 (2008) .
Blank and Borge (2008)
B. Blank and M. Borge, Progress in Particle and Nuclear Physics 60 , 403 (2008) .
Pfßtzner et al. (2023)
M. Pfßtzner, I. Mukha,
 and S. Wang, Progress in Particle and Nuclear Physics 132 , 104050 (2023) .
Qi et al. (2019)
C. Qi, R. Liotta,  and R. Wyss, Progress in Particle and Nuclear Physics 105 , 214 (2019) .
Goldansky (1960)
V. I. Goldansky, Nucl. Phys. 19 , 482
(1960).
Goldansky (1961)
V. I. Goldansky, Nucl. Phys. 27 , 648
(1961).
Oishi et al. (2014)
T. Oishi, K. Hagino,  and H. Sagawa, Phys. Rev. C 90 , 034303 (2014) .
Oishi et al. (2017)
T. Oishi, M. Kortelainen,
 and A. Pastore, Phys. Rev. C 96 , 044327 (2017) .
Wang and Nazarewicz (2021)
S. M. Wang and W. Nazarewicz, Phys. Rev. Lett. 126 , 142501 (2021) .
Cirelâson (1980)
B. S. Cirelâson, Letters in Mathematical Physics 4 , 93 (1980) .
Nielsen and Chuang (2000)
M. Nielsen and I. Chuang, Quantum Computation and
Quantum Information  (Cambridge University Press, Cambridge, UK, 2000).
Suzuki and Ikeda (1988)
Y. Suzuki and K. Ikeda, Phys. Rev. C 38 , 410 (1988) .
Bertsch and Esbensen (1991)
G. Bertsch and H. Esbensen, Annals of Physics 209 , 327 (1991) .
Esbensen et al. (1997)
H. Esbensen, G. F. Bertsch,  and K. Hencken, Phys. Rev. C 56 , 3054 (1997) .
Hagino and Sagawa (2005)
K. Hagino and H. Sagawa, Phys. Rev. C 72 , 044321 (2005) .
Thompson et al. (1977)
D. Thompson, M. Lemere,  and Y. Tang, Nuclear Physics A 286 , 53 (1977) .
Ajzenberg-Selove (1988)
F. Ajzenberg-Selove, Nuclear Physics A 490 , 1 (1988) , note: several versions with the same title has been published.
Ajzenberg-Selove (1991)
F. Ajzenberg-Selove, Nuclear Physics A 523 , 1 (1991) , with revised version at http://www.tunl.duke.edu/nucldata/.