Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Spin entanglement in two-proton emission from 6Be

Tomohiro Oishi tomohiro.oishi@ribf.riken.jp RIKEN Nishina Center for Accelerator-Based Science, Wako 351-0198, Japan
Abstract

This paper presents an evaluation of coupled-spin entanglement in the two-proton (2⁢p2𝑝2p2 italic_p) radioactive emission. The three-body model of 6Be 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𝑆Sitalic_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⁢p2𝑝2p2 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⁢p2𝑝2p2 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𝑆Sitalic_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 2H+p⟶2superscript⟶2𝑝absentp\longrightarrow^{2}italic_p ⟶ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPTHe+n𝑛nitalic_n. The spin-correlation function as the CHSH indicator is deduced as Sexpt=2.83±0.24stat±0.07syssubscript𝑆exptplus-or-minus2.83subscript0.24statsubscript0.07sysS_{\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⁢p2𝑝2p2 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⁢p2𝑝2p2 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⁢p2𝑝2p2 italic_p interaction inside finite systems, being in a contrast to the vacuum 2⁢p2𝑝2p2 italic_p interaction supporting no bound state.

In this paper, I evaluate the quantum entaglement in the 2⁢p2𝑝2p2 italic_p emission of 6Be. Main motivation is to utilize this entanglement as a probe into the effective interaction. The 6Be is the lightest 2⁢p2𝑝2p2 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𝑆Sitalic_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⁢p2𝑝2p2 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/2121/21 / 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-productsubscript^𝐴1𝜃1subscript^𝐵1𝜃2tensor-productsubscript^𝐴2𝜃1subscript^𝐵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-productsubscript^𝐴1𝜃1subscript^𝐵2𝜃2tensor-productsubscript^𝐴2𝜃1subscript^𝐵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 θ𝜃\thetaitalic_θ are given as

A^1=σ^z,A^2,θ=σ^z⁢cos⁡2⁢θ+σ^x⁢sin⁡2⁢θ,formulae-sequencesubscript^𝐴1subscript^𝜎𝑧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Ψ12\ket{\Psi(1,2)}| start_ARG roman_Ψ ( 1 , 2 ) end_ARG ⟩, their expectation values are obtained:

⟨Ai⁢Bj⟩=⟨Ψ⁢(1,2)⁢∣A^i,θ⁢(1)⊗B^j,θ⁢(2)∣⁢Ψ⁢(1,2)⟩.expectationsubscript𝐴𝑖subscript𝐵𝑗expectationΨ12delimited-∣∣tensor-productsubscript^𝐴𝑖𝜃1subscript^𝐵𝑗𝜃2Ψ12\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𝑆Sitalic_S is determined as

S=max⁡{|S−⁣+⁣++|,|S+⁣−⁣++|,|S+⁣+⁣−+|,|S+⁣+⁣+−|},𝑆subscript𝑆absentsubscript𝑆absentsubscript𝑆absentsubscript𝑆absentS=\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== −⟨A1⁢B1⟩+⟨A2⁢B1⟩+⟨A1⁢B2⟩+⟨A2⁢B2⟩,expectationsubscript𝐴1subscript𝐵1expectationsubscript𝐴2subscript𝐵1expectationsubscript𝐴1subscript𝐵2expectationsubscript𝐴2subscript𝐵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== ⟨A1⁢B1⟩−⟨A2⁢B1⟩+⟨A1⁢B2⟩+⟨A2⁢B2⟩,expectationsubscript𝐴1subscript𝐵1expectationsubscript𝐴2subscript𝐵1expectationsubscript𝐴1subscript𝐵2expectationsubscript𝐴2subscript𝐵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== ⟨A1⁢B1⟩+⟨A2⁢B1⟩−⟨A1⁢B2⟩+⟨A2⁢B2⟩,expectationsubscript𝐴1subscript𝐵1expectationsubscript𝐴2subscript𝐵1expectationsubscript𝐴1subscript𝐵2expectationsubscript𝐴2subscript𝐵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== ⟨A1⁢B1⟩+⟨A2⁢B1⟩+⟨A1⁢B2⟩−⟨A2⁢B2⟩.expectationsubscript𝐴1subscript𝐵1expectationsubscript𝐴2subscript𝐵1expectationsubscript𝐴1subscript𝐵2expectationsubscript𝐴2subscript𝐵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)⟩=|↑⁣↓⁣−⁣↓⁣↑⟩/2ketΨ12ket↑↓↓↑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𝜃𝜋43𝜋45𝜋4\theta=\pi/4,~{}3\pi/4,~{}5\pi/4italic_θ = italic_π / 4 , 3 italic_π / 4 , 5 italic_π / 4, and 7⁢π/47𝜋47\pi/47 italic_π / 4, this indicator has the maximum value, S=2⁢2𝑆22S=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/4italic_θ = italic_π / 4, except when modified.

Refer to caption
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)⟩=|en⟩ketΨ12ketsubscript𝑒𝑛\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𝑆22S=2\sqrt{2}italic_S = 2 square-root start_ARG 2 end_ARG when θ=π/4𝜃𝜋4\theta=\pi/4italic_θ = italic_π / 4. Here the Bell states read Nielsen and Chuang (2000)

|e1⟩=12⁢|↑⁣↑⁣+⁣↓⁣↓⟩,|e2⟩=i2⁢|↑⁣↑⁣−⁣↓⁣↓⟩,formulae-sequenceketsubscript𝑒112ket↑↑↓↓ketsubscript𝑒2𝑖2ket↑↑↓↓\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 ⟩ ,
|e3⟩=i2⁢|↑⁣↓⁣+⁣↓⁣↑⟩,|e4⟩=12⁢|↑⁣↓⁣−⁣↓⁣↑⟩.formulae-sequenceketsubscript𝑒3𝑖2ket↑↓↓↑ketsubscript𝑒412ket↑↓↓↑\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/2121/21 / 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=14Gi⁢(𝒓1,𝒓2)⁢|ei⟩,ketΨ12superscriptsubscript𝑖14subscript𝐺𝑖subscript𝒓1subscript𝒓2ketsubscript𝑒𝑖\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 Gi⁢(𝒓1,𝒓2)subscript𝐺𝑖subscript𝒓1subscript𝒓2G_{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^12subscript^𝑆12\hat{S}_{12}over^ start_ARG italic_S end_ARG start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT, its eigenvalues are given as S^122⁢|dk⟩=S⁢(S+1)⁢|dk⟩subscriptsuperscript^𝑆212ketsubscript𝑑𝑘𝑆𝑆1ketsubscript𝑑𝑘\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⁢|dk⟩=V⁢|dk⟩subscript^𝑆12𝑧ketsubscript𝑑𝑘𝑉ketsubscript𝑑𝑘\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

|d1⟩ketsubscript𝑑1\displaystyle\ket{d_{1}}| start_ARG italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG ⟩ =\displaystyle== |S=1,V=+1⟩=|↑↑⟩,ketformulae-sequence𝑆1𝑉1ket↑absent↑\displaystyle\ket{S=1,V=+1}=\ket{\uparrow\uparrow},| start_ARG italic_S = 1 , italic_V = + 1 end_ARG ⟩ = | start_ARG ↑ ↑ end_ARG ⟩ ,
|d2⟩ketsubscript𝑑2\displaystyle\ket{d_{2}}| start_ARG italic_d start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG ⟩ =\displaystyle== |S=1,V=−1⟩=|↓↓⟩,ketformulae-sequence𝑆1𝑉1ket↓absent↓\displaystyle\ket{S=1,V=-1}=\ket{\downarrow\downarrow},| start_ARG italic_S = 1 , italic_V = - 1 end_ARG ⟩ = | start_ARG ↓ ↓ end_ARG ⟩ ,
|d3⟩ketsubscript𝑑3\displaystyle\ket{d_{3}}| start_ARG italic_d start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT end_ARG ⟩ =\displaystyle== |S=1,V=0⟩=12⁢|↑⁣↓⁣+⁣↓⁣↑⟩,ketformulae-sequence𝑆1𝑉012ket↑↓↓↑\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 ⟩ ,
|d4⟩ketsubscript𝑑4\displaystyle\ket{d_{4}}| start_ARG italic_d start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT end_ARG ⟩ =\displaystyle== |S=0,V=0⟩=12⁢|↑⁣↓⁣−⁣↓⁣↑⟩.ketformulae-sequence𝑆0𝑉012ket↑↓↓↑\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𝑈Uitalic_U is formulated:

|e1⟩|e2⟩|e3⟩|e4⟩=U⁢|d1⟩|d2⟩|d3⟩|d4⟩=(121200i2−i20000i00001)⁢|d1⟩|d2⟩|d3⟩|d4⟩.ketsubscript𝑒1ketsubscript𝑒2ketsubscript𝑒3ketsubscript𝑒4𝑈ketsubscript𝑑1ketsubscript𝑑2ketsubscript𝑑3ketsubscript𝑑4121200𝑖2𝑖20000𝑖00001ketsubscript𝑑1ketsubscript𝑑2ketsubscript𝑑3ketsubscript𝑑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)⟩=∑iFi⁢(𝒓1,𝒓2)⁢|di⟩ketΨ12subscript𝑖subscript𝐹𝑖subscript𝒓1subscript𝒓2ketsubscript𝑑𝑖\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, Gi⁢(𝒓1,𝒓2)=∑k(U†)i⁢k⁢Fk⁢(𝒓1,𝒓2)subscript𝐺𝑖subscript𝒓1subscript𝒓2subscript𝑘subscriptsuperscript𝑈†𝑖𝑘subscript𝐹𝑘subscript𝒓1subscript𝒓2G_{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 mCsubscript𝑚𝐶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)+vp⁢p⁢(𝒓1,𝒓2)+𝒑1⋅𝒑2mC,subscript^𝐻3𝐵^ℎsubscript𝒓1^ℎsubscript𝒓2subscript𝑣𝑝𝑝subscript𝒓1subscript𝒓2⋅subscript𝒑1subscript𝒑2subscript𝑚𝐶\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^⁢(ri)=−ℏ22⁢μ⁢d2d⁢ri2+ℏ22⁢μ⁢l⁢(l+1)ri2+V⁢(ri),^ℎsubscript𝑟𝑖superscriptPlanck-constant-over-2-pi22𝜇superscript𝑑2𝑑superscriptsubscript𝑟𝑖2superscriptPlanck-constant-over-2-pi22𝜇𝑙𝑙1superscriptsubscript𝑟𝑖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 μ=mp⁢mC/(mp+mC)𝜇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 mp=938.272subscript𝑚𝑝938.272m_{p}=938.272italic_m start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT = 938.272 MeV/c2absentsuperscript𝑐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𝑘kitalic_kth proton-alpha subsystem. For its interaction V⁢(ri)𝑉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

vp⁢p⁢(𝒓1,𝒓2)=vn⁢u⁢c⁢l⁢(d)+e24⁢π⁢ϵ0⁢1d+va⁢d⁢d⁢(𝒓1,𝒓2),subscript𝑣𝑝𝑝subscript𝒓1subscript𝒓2subscript𝑣𝑛𝑢𝑐𝑙𝑑superscript𝑒24𝜋subscriptitalic-ϵ01𝑑subscript𝑣𝑎𝑑𝑑subscript𝒓1subscript𝒓2v_{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𝒓2subscript𝒓1d=\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):

vn⁢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== [VR⁢e−aR⁢d2+VS⁢e−aS⁢d2]⁢P^S=0delimited-[]subscript𝑉𝑅superscript𝑒subscript𝑎𝑅superscript𝑑2subscript𝑉𝑆superscript𝑒subscript𝑎𝑆superscript𝑑2subscript^𝑃𝑆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)
+[VR⁢e−aR⁢d2+VT⁢e−aT⁢d2]⁢P^S=1.delimited-[]subscript𝑉𝑅superscript𝑒subscript𝑎𝑅superscript𝑑2subscript𝑉𝑇superscript𝑒subscript𝑎𝑇superscript𝑑2subscript^𝑃𝑆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=0subscript^𝑃𝑆0\hat{P}_{S=0}over^ start_ARG italic_P end_ARG start_POSTSUBSCRIPT italic_S = 0 end_POSTSUBSCRIPT (P^S=1subscript^𝑃𝑆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 VR=200subscript𝑉𝑅200V_{R}=200italic_V start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT = 200 MeV, VS=−91.85subscript𝑉𝑆91.85V_{S}=-91.85italic_V start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT = - 91.85 MeV, VT=−178subscript𝑉𝑇178V_{T}=-178italic_V start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT = - 178 MeV, aR=1,487subscript𝑎𝑅1487a_{R}=1,487italic_a start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT = 1 , 487 fm-2, aS=0.465subscript𝑎𝑆0.465a_{S}=0.465italic_a start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT = 0.465 fm-2, and aT=0.639subscript𝑎𝑇0.639a_{T}=0.639italic_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 va⁢d⁢dsubscript𝑣𝑎𝑑𝑑v_{add}italic_v start_POSTSUBSCRIPT italic_a italic_d italic_d end_POSTSUBSCRIPT is employed:

va⁢d⁢d⁢(𝒓1,𝒓2)=w0⁢e−(R−R0)2/B02⁢δ⁢(𝒓1−𝒓2),subscript𝑣𝑎𝑑𝑑subscript𝒓1subscript𝒓2subscript𝑤0superscript𝑒superscript𝑅subscript𝑅02superscriptsubscript𝐵02𝛿subscript𝒓1subscript𝒓2v_{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𝒓1subscript𝒓22R=\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 |, R0=1.68subscript𝑅01.68R_{0}=1.68italic_R start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1.68 fm, B0=0.6⁢R0subscript𝐵00.6subscript𝑅0B_{0}=0.6R_{0}italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.6 italic_R start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, and w0=−470subscript𝑤0470w_{0}=-470italic_w start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = - 470 MeV⋅⋅\cdot⋅fm3. This additional term is necessary to reproduce the experimental energy: the 2⁢p2𝑝2p2 italic_p energy and width are obtained as E2⁢p=1.356subscript𝐸2𝑝1.356E_{2p}=1.356italic_E start_POSTSUBSCRIPT 2 italic_p end_POSTSUBSCRIPT = 1.356 MeV and Γ2⁢p=0.055subscriptΓ2𝑝0.055\Gamma_{2p}=0.055roman_Γ start_POSTSUBSCRIPT 2 italic_p end_POSTSUBSCRIPT = 0.055 MeV, respectively, whereas the experimental data read E2⁢p=1.372⁢(5)subscript𝐸2𝑝1.3725E_{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.0926\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.

Refer to caption
FIG. 2: Time-dependent decaying probability of 2⁢p2𝑝2p2 italic_p emission from 6Be. Its spin-singlet (S12=0subscript𝑆120S_{12}=0italic_S start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT = 0) component is also plotted. The factor fp⁢p=1subscript𝑓𝑝𝑝1f_{pp}=1italic_f start_POSTSUBSCRIPT italic_p italic_p end_POSTSUBSCRIPT = 1 indicates that the proton-proton interaction is adjusted to reproduce the experimental 2⁢p2𝑝2p2 italic_p-energy.

Results and Discussions - The 2⁢p2𝑝2p2 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-piketΨ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⁢p2𝑝2p2 italic_p state inside the Coulomb barrier Oishi et al. (2014). The decaying state |Ψd⁢(t)⟩ketsubscriptΨ𝑑𝑡\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)⟩,ketsubscriptΨ𝑑𝑡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:

Pd⁢e⁢c⁢a⁢y⁢(t)=⟨Ψd⁢(t)∣Ψd⁢(t)⟩=1−|β⁢(t)|2.subscript𝑃𝑑𝑒𝑐𝑎𝑦𝑡expectationconditionalsubscriptΨ𝑑𝑡subscriptΨ𝑑𝑡1superscript𝛽𝑡2P_{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, |d4⟩=|e4⟩=|↑⁣↓⁣−⁣↓⁣↑⟩/2ketsubscript𝑑4ketsubscript𝑒4ket↑↓↓↑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,

Ps⁢u⁢r⁢v⁢(t)≡⟨Ψ⁢(0)∣Ψ⁢(t)⟩2=|β⁢(t)|2=1−Pd⁢e⁢c⁢a⁢y⁢(t),subscript𝑃𝑠𝑢𝑟𝑣𝑡superscriptexpectationconditionalΨ0Ψ𝑡2superscript𝛽𝑡21subscript𝑃𝑑𝑒𝑐𝑎𝑦𝑡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, Ps⁢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.055subscriptΓ2𝑝Planck-constant-over-2-pi𝜏0.055\Gamma_{2p}=\hbar/\tau\cong 0.055roman_Γ 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)|2subscript𝜌𝑑𝑡subscript𝒓1subscript𝒓2superscriptsubscriptΨ𝑑𝑡subscript𝒓1subscript𝒓22\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⁢p2𝑝2p2 italic_p emission can be interpreted as the diproton-correlating emission Oishi et al. (2014, 2017).

Refer to caption
FIG. 3: CHSH indicator of the time-dependent decaying state |Ψd⁢(t)⟩ketsubscriptΨ𝑑𝑡\ket{\Psi_{d}(t)}| start_ARG roman_Ψ start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT ( italic_t ) end_ARG ⟩ of the 2⁢p2𝑝2p2 italic_p emission of 6Be ⟶α+p+p⟶absent𝛼𝑝𝑝\longrightarrow\alpha+p+p⟶ italic_α + italic_p + italic_p. Here θ=π/4𝜃𝜋4\theta=\pi/4italic_θ = italic_π / 4. The factor fp⁢psubscript𝑓𝑝𝑝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=14Gi(t,𝒓1,𝒓2))|ei⟩.\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 2222, i.e. the LHV-theory limit, during the time evolution. For c⁢t≥1500𝑐𝑡1500ct\geq 1500italic_c italic_t ≥ 1500 fm, I obtain S≅2.72𝑆2.72S\cong 2.72italic_S ≅ 2.72 with the default setting of vp⁢psubscript𝑣𝑝𝑝v_{pp}italic_v start_POSTSUBSCRIPT italic_p italic_p end_POSTSUBSCRIPT (fp⁢p=1.0subscript𝑓𝑝𝑝1.0f_{pp}=1.0italic_f start_POSTSUBSCRIPT italic_p italic_p end_POSTSUBSCRIPT = 1.0) for reproducing E2⁢p=1.356subscript𝐸2𝑝1.356E_{2p}=1.356italic_E start_POSTSUBSCRIPT 2 italic_p end_POSTSUBSCRIPT = 1.356 MeV of 6Be 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⁢(ri)𝑉subscript𝑟𝑖V(r_{i})italic_V ( italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ). These interactions do not harm the 2⁢p2𝑝2p2 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:

vp⁢p⁢(𝒓1,𝒓2)⟶fp⁢p⁢vp⁢p⁢(𝒓1,𝒓2).⟶subscript𝑣𝑝𝑝subscript𝒓1subscript𝒓2subscript𝑓𝑝𝑝subscript𝑣𝑝𝑝subscript𝒓1subscript𝒓2v_{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 fp⁢p=1subscript𝑓𝑝𝑝1f_{pp}=1italic_f start_POSTSUBSCRIPT italic_p italic_p end_POSTSUBSCRIPT = 1, the vp⁢psubscript𝑣𝑝𝑝v_{pp}italic_v start_POSTSUBSCRIPT italic_p italic_p end_POSTSUBSCRIPT approximates the experimental energy, E2⁢p=1.356subscript𝐸2𝑝1.356E_{2p}=1.356italic_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 6Be, it goes below the LHV-theory limit, S<2𝑆2S<2italic_S < 2, with fp⁢p≤0.2subscript𝑓𝑝𝑝0.2f_{pp}\leq 0.2italic_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 vp⁢psubscript𝑣𝑝𝑝v_{pp}italic_v start_POSTSUBSCRIPT italic_p italic_p end_POSTSUBSCRIPT weakened, the 2⁢p2𝑝2p2 italic_p energy and width also change: E2⁢p=2.864subscript𝐸2𝑝2.864E_{2p}=2.864italic_E start_POSTSUBSCRIPT 2 italic_p end_POSTSUBSCRIPT = 2.864 MeV and Γ2⁢p=0.8579subscriptΓ2𝑝0.8579\Gamma_{2p}=0.8579roman_Γ start_POSTSUBSCRIPT 2 italic_p end_POSTSUBSCRIPT = 0.8579 MeV with fp⁢p=0.5subscript𝑓𝑝𝑝0.5f_{pp}=0.5italic_f start_POSTSUBSCRIPT italic_p italic_p end_POSTSUBSCRIPT = 0.5; E2⁢p=2.925subscript𝐸2𝑝2.925E_{2p}=2.925italic_E start_POSTSUBSCRIPT 2 italic_p end_POSTSUBSCRIPT = 2.925 MeV and Γ2⁢p≅2.9subscriptΓ2𝑝2.9\Gamma_{2p}\cong 2.9roman_Γ start_POSTSUBSCRIPT 2 italic_p end_POSTSUBSCRIPT ≅ 2.9 MeV with fp⁢p=0.2subscript𝑓𝑝𝑝0.2f_{pp}=0.2italic_f start_POSTSUBSCRIPT italic_p italic_p end_POSTSUBSCRIPT = 0.2. Thus, the short-lived (broad-width) 2⁢p2𝑝2p2 italic_p state has the weaker spin entanglement.

Summary - In this paper, the CHSH indicator for the 2⁢p2𝑝2p2 italic_p emission from 6Be is evaluated in order to measure its coupled-spin entanglement. As a product of proton-proton interaction, S≅2.72𝑆2.72S\cong 2.72italic_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 6Be nucleus. Whether the spin entanglement exists commonly in other 2⁢p2𝑝2p2 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⁢p2𝑝2p2 italic_p emitters, including the 6Be, 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⁢p2𝑝2p2 italic_p emitters is rapidly in progress, e.g. that in RIKEN RIBF. Confirmation of spin entanglement in 2⁢p2𝑝2p2 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