## Clustering structure of nuclei in deep inelastic

## processes

## M. Hirai a , S. Kumano b,c,d , K. Saito d,e , T. Watanabe e

- a Nippon Institute of Technology, Saitama 345-8501, Japan
- b KEK Theory Center, Institute of Particle and Nuclear Studies, KEK and Department of Particle and Nuclear Studies, Graduate University for Advanced Studies 1-1, Ooho, Tsukuba, Ibaraki, 305-0801, Japan
- c Theory Group, Particle and Nuclear Physics Division, J-PARC Center 203-1, Shirakata, Tokai, Ibaraki, 319-1106, Japan
- d J-PARC Branch, KEK Theory Center, Institute of Particle and Nuclear Studies, KEK 203-1, Shirakata, Tokai, Ibaraki, 319-1106, Japan
- e Department of Physics, Faculty of Science and Technology, Tokyo University of Science 2641, Yamazaki, Noda, Chiba, 278-8510, Japan

E-mail: shunzo.kumano@kek.jp

Abstract. A clustering aspect is explained for the 9 Be nucleus in charged-lepton deep inelastic scattering. Nuclear modifications of the structure function F 2 are studied by the ratio R EMC = F A 2 /F D 2 , where A and D are a nucleus and the deuteron, respectively. In a JLab experiment, an unexpectedly large nuclear modification slope | dR EMC /dx | was found for 9 Be, which could be related to its clustering structure. We investigated a mean conventional part of a nuclear structure function F A 2 by a convolution description with nucleon momentum distributions calculated by antisymmetrized (or fermionic) molecular dynamics (AMD) and also by a simple shell model. We found that clustering effects are small in the conventional part, so that the JLab result could be associated with an internal nucleon modification or a short-range nuclear correlation which is caused by high densities due to cluster formation.

## 1. Introduction

Nuclear modifications of the structure function F 2 are known as the EMC effect [1], and such effects are shown by the ratio R EMC ( x ) = F A 2 ( x /F ) D 2 ( x ), where A and D indicate a nucleus and the deuteron, respectively, and x is the Bjorken scaling variable. Various mechanisms contribute to the modifications. At large x &gt; 0 7, nucleon Fermi motion in a nucleus gives rise to increase of . the ratio R EMC ( x ). At medium x (0 3 . &lt; x &lt; 0 7), nuclear binding is the major source of nuclear . effects. In addition, modifications of internal nucleon structure could contribute. At small x &lt; 0 05, . the nuclear shadowing effects become prominent. The nuclear structure functions have been measured from small x to large x for many nuclei, and the nuclear modifications have been determined by global analyses of experimental data [2].

Nuclei are basically described by a shell model, where nucleons are assumed to move in an average central potential created by interactions of all the nucleons in a nucleus. However, cluster structure in a nucleus is not easily described by the shell model with a limited number of shells. We know that some nuclei in the region of the mass number A ∼ 10 have clustering configurations. It is interesting to study whether such clustering features appear in the nuclear structure functions F A 2 .

In 2009, an anomalously large nuclear effect was reported by the JLab measurement on the beryllium-9 structure function F 2 as shown in Fig. 1 [3], where the slope is shown as a functions of the scaled average nuclear density. The scaled means that the factor ( A -1) /A is multiplied for removing the struck nucleon. Although other data are on the smooth curve, the beryllium-9 data is much larger than the one expected from its average nuclear density.

It is known that the beryllium-9 has a clustering configuration with two α nuclei and surrounding neutron clouds. Therefore, such an unexpected nuclear modification could come from the cluster formation in the nucleus. We

Figure 1. Nuclear modification slope | d F ( A 2 /F D 2 ) /dx | .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0660v1-figure-a421ba89-b110-4068-aaee-0e7beeb886e0.jpg)

investigated such a possibility by using a convolution description for the structure function F A 2 with nucleon momentum distributions calculated by antisymmetrized molecular dynamics (AMD) and also by a simple shell model for comparison [4]. Then, possible reasons are suggested for the large slope of 9 Be due to the cluster formation in the nucleus.

## 2. Clustering effects in nuclear structure functions

The standard method for describing the nuclear structure functions F A 2 is to use the convolution description [5] for describing a mean conventional part. The structure function F A 2 is given by the nucleonic one F N 2 convoluted with a nucleon momentum distribution in the nucleus f ( y ) [1, 5] as shown in Fig. 2:

$$\L ^ { \cdot } \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \rho \, \nu \, \nu \, \nu \, \nu \, \nu \, \rho \, \nu \, \nu \, \rho \, \nu \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rp \, \rho \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rrp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \, \rp \,
\text{where $u$ is the momentum fraction $u$ = $M_{\lambda}$ $n_{\lambda}$ $u$} \cdot a / ( M _ { \lambda } n _ { \lambda } \cdot a ) \sim A$$

Figure 2. Convolution for F A 2 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0660v1-figure-ad5d14cc-e47f-4cde-93fe-ba884b58eb57.jpg)

where y is the momentum fraction y = M p A N · q/ M ( N p A · q ) /similarequal Ap + N / p ( + A ) with a light-cone momentum p + .

The momentum density | φ /vector ( p N ) | 2 is calculated by two methods, antisymmetrized molecular dynamics (AMD) [6] or fermionic molecular dynamics (FMD) [7], and a simple shell model for comparison. The AMD and FMD are essentially the same, so that we use the notation AMD hereafter. The AMD is a variational method, and its advantage is that there is no a priori assumption on nuclear structure whether it is a shell or cluster-like configuration. A nuclear wave function is given by the Slater determinant of single-particle wave functions

$$| \Phi ( \vec { r } _ { 1 }, \vec { r } _ { 2 }, \cdot \cdot \cdot, \vec { r } _ { A } ) \rangle = \frac { 1 } { \sqrt { A }! } \det [ \varphi _ { 1 } ( \vec { r } _ { 1 } ), \varphi _ { 2 } ( \vec { r } _ { 2 } ), \cdot \cdot \cdot, \varphi _ { A } ( \vec { r } _ { A } ) ].$$

Here, a nucleon is described by the single-particle wave function ϕ i ( /vector r j ) = φ i ( /vector r j ) χ τ , i i where χ i and τ i indicate spin and isospin states, respectively. The function φ i ( /vector r j ) is the space part of the wave function, and it is assumed to be given by the Gaussian functional form:

$$\phi _ { i } ( \vec { r } _ { j } ) = \left ( \frac { 2 \nu } { \pi } \right ) ^ { 3 / 4 } \exp \left [ - \nu \left ( \vec { r } _ { j } - \frac { \vec { Z } _ { i } } { \sqrt { \nu } } \right ) ^ { 2 } \right ], \\ \text{arameter to express the extent of the wave packet.  The center of the wave packet}$$

where ν is a parameter to express the extent of the wave packet. The center of the wave packet is given by /vector Z / i √ ν , where /vector Z i is a complex variational parameter. The variational parameters are then determined by minimizing the system energy with a frictional-cooling method. We use effective NN interactions which could describe gross properties of nuclei:

$$\Lambda _ { 2 } & = ( 1 - m - m P _ { \sigma } P _ { \tau } ) \left [ v _ { 2 1 } e ^ { - ( r / r _ { 2 1 } ) ^ { 2 } } + v _ { 2 2 } e ^ { - ( r / r _ { 2 2 } ) ^ { 2 } } \right ], \\ 3 \text{-body} \, \ V _ { 3 } & = v _ { 3 } \, \delta ^ { 3 } ( \vec { r } _ { 1 } - \vec { r } _ { 2 } ) \, \delta ^ { 3 } ( \vec { r } _ { 2 } - \vec { r } _ { 3 } ), \\ & \quad \text{LS} \, \ V _ { L S } = v _ { L S } \left [ e ^ { - ( r / r _ { L S 1 } ) ^ { 2 } } - e ^ { - ( r / r _ { L S 2 } ) ^ { 2 } } \right ] P ( ^ { 3 } O ) \, \vec { L } \cdot \vec { S }, \\ \intertext { \ e r m. } \Lambda _ { 2 } \,, \, v _ { 2 1 }, \, v _ { 2 2 }, \, r _ { 2 1 }, \, r _ { 2 2 }, \, v _ { 3 }, \, v _ { L S }, \, r _ { L S 1 }, \, \text{and} \, r _ { L S 2 } \, \text{are constants.} } \Lambda _ { 2 } \,, \, \Lambda _ { 2 } \,, \, \Lambda _ { 2 } \,, \, \Lambda _ { 2 } \,, \, \Lambda _ { 2 } \,, \, \Lambda _ { 2 } \,, \, \Lambda _ { 2 } \,, \, \Lambda _ { 2 } \,, \, \Lambda _ { 2 } \,, \, \Lambda _ { 2 } \,, \, \Lambda _ { 2 } \,, \, \Lambda _ { 2 } \,, \, \Lambda _ { 2 } \,, \, \Lambda _ { 2 } \,, \, \Lambda _ { 2 } \,, \, \Lambda _ { 2 } \,, \, \Lambda }$$

where m v , 21 , v 22 , r 21 , r 22 , v 3 , v LS , r LS 1 , and r LS 2 are constants.

The momentum distributions of Eq. (1) are also calculated in a simple shell model, where the potential is the harmonic-oscillator type ( M ω r / N 2 2 2), for clarifying cluster effects by comparing two results. Its wave function is separated into radial- and angular-dependent parts: ψ n/lscriptm ( r, θ, φ ) = R n/lscript ( r Y ) /lscriptm ( θ, φ ) where r , θ , and φ are spherical coordinates, and n /lscript , , and m are radial, azimuthal, and magnetic quantum numbers, respectively. The function Y /lscriptm ( θ, φ ) is the spherical harmonics, and the radial wave function is given by

$$R _ { n \ell } ( r ) = \sqrt { \frac { 2 \kappa ^ { 2 \ell + 3 } ( n - 1 )! } { [ \Gamma ( n + \ell + 1 / 2 ) ] ^ { 3 } } } r ^ { \ell } e ^ { - \frac { 1 } { 2 } \kappa ^ { 2 } r ^ { 2 } } L _ { n - 1 } ^ { \ell + 1 / 2 } ( \kappa ^ { 2 } r ^ { 2 } ),$$

where L /lscript +1 2 / n -1 ( x ) is the Laguerre polynomial, and κ is defined by κ ≡ √ M ω N .

Figure 3. Density of 9 Be by AMD.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0660v1-figure-9bf68341-ff2e-4158-9fd5-ee72a8fd158c.jpg)

As a result, we obtain the two-dimensional spacial density of the 9 Be nucleus in Fig. 3. It obviously indicates the existence of two α clusters, together with surrounding neutron clouds. In a simple shell model, such clusters do not appear, and the density is a monotonic distribution. If the averages are taken over polar and azimuthal angles θ and φ , the distributions becomes the ones in Fig. 4. After the averages, the cluster structure is no longer apparent. Because the nucleons are confined mainly in two separate regions, the radial distribution is shifted toward larger region for 9 Be in Fig. 4 as shown by the solid curve in comparison with the dashed one.

The momentum distributions in Fig. 5 are calculated from the spacial distributions. Both AMD and shell-model results are shown for 4 He and 9 Be. The momentum distribution of the AMD is shifted toward the high-momentum region in 9 Be, which is caused by the fact that the dense clusters are formed within the 9 Be nucleus. If nucleons are confined in the small space regions of the clusters as indicated by the AMD, it leads to an increase of high momentum components.

Figure 4. Spatial distributions.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0660v1-figure-4f29fec8-94d7-4ed9-b867-9d6aa5444ef0.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0660v1-figure-c236b4a8-a033-4b86-a32f-51d9894f1bed.jpg)

Figure 5. Momentum distributions.

The modifications of the momentum distributions should affect the nuclear structure function in Eq.(1), and the results are shown for the structure function F 2 in Fig. 6. There are two theoretical curves by the AMD and shell model. In general, the theoretical ratios are consistent with the data in the sense that the ratio decreases at medium x and it increases at large x . These decrease and increase are caused by the nuclear binding and the nucleon's Fermi motion, respectively. However, the simple convolution model is not sufficient to explain the experimental data. The differences should be attributed to the effects of shortrange nucleon-nucleon correlations and modifications of

Figure 6. F 9 Be 2 /F D 2 by shell and AMD models.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0660v1-figure-7881d41e-3d43-4be8-bf9c-00a1e368c206.jpg)

internal nucleon structure. The correlation effects change the theoretical ratios toward the experimental data at x = 0 6 . -0 8. . Since we are interested in an order-of-magnitude estimate of clustering effects, we did not step into such details in this work. In any case, the differences of the two curves, namely the clustering effects, are not large and they are of the same order of experimental errors. In Fig. 6, we found that the clustering effects are rather small in the mean conventional part of F A 2 . Therefore, there should be other sources for the large slope of 9 Be. We consider that the nuclear structure functions consist of the mean conventional part and the remaining one depending on the maximum local density:

$$F _ { 2 } ^ { A } = ( \text{mean part} ) + ( \text{part created by large densities due to cluster formation). \quad \text{(6)}$$

The first part was calculated by the convolution model as shown in Fig. 6, and the second part seems to be important for explaining the JLab 9 Be result. The second term is associated with the modification of internal nucleon structure and/or a short-range correlation caused by the cluster formation.

Although we do not step into detailed calculations in this work, we would like to find whether such an idea works. If the physics mechanism behind the JLab 9 Be result is related to the cluster formation, observables could be interpreted by the cluster densities. For such estimates, we calculated maximum local densities by both AMD and shell models, rather than the average nuclear densities. In Fig. 7, the JLab data are plotted by the maximum local densities calculated in the AMD and shell model. The maximum densities are almost the same in the nuclei, 3 He, 4 He, and 12 C, whereas the AMD density is much different from the shell-model one for 9 Be. According to the shell model with a limited number of shells, the cluster formation, namely

Figure 7. Nuclear modification slopes shown by maximum local densities.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0660v1-figure-6ffe14a2-80f8-4182-91e4-f5cd4275e117.jpg)

the large densities, cannot be described, and it leads to a small maximum density of 9 Be in Fig. 7. It is also interesting to find that there is also a small cluster effect in the 12 C. From Fig. 7, we find that all the data are on the smooth curve if the density is calculated by the AMD for 9 Be; however, it is far away from the curve if it is evaluated by the shell model. The figure indicates that the large nuclear modification slope | dR EMC /dx | of 9 Be should be related to the high densities created by the cluster formation.

The JLab measurement could be the first result on a clustering aspect of a nucleus in highenergy nuclear reactions. It could be associated with internal nucleon modifications and/or short-range nuclear correlations, which are caused by the high densities due to the cluster configuration. In future, there will be measurement for other nuclei, such as 6 Li, 7 Li, 10 B, and 11 B [8], so that much detailed information will become available.

## 3. Summary

The large nuclear modification slope of the structure function F 2 was found for 9 Be, which has a typical cluster configuration. We showed that the clustering effects are rather small in the mean conventional part by using the convolution model for the structure functions F A 2 with the momentum distributions calculated in the AMD and also in the shell model. The JLab result should be associated with the high-density formation due to the clusters in the 9 Be nucleus. We showed it by plotting the JLab data as a function of the maximum local density of a nucleus. We have not clarified the physics origin behind its result. However, it is likely to be modifications of internal nucleon structure in the nuclear medium and/or short-range nucleonnucleon correlations in the dense clusters. There is an approved experiment at JLab to measure the EMC effects in the mass region A ∼ 10, so that much details will be investigated in the near future for possible clustering effects in deep inelastic scattering (DIS). The cluster physics has been investigated mainly in low energy nuclear reactions; however, DIS processes could shed light on new aspects in both cluster physics and nuclear structure functions. Cluster physicists have an opportunity to play a major role in high-energy hadron physics.

## Acknowledgements

This work was supported by the MEXT KAKENHI Grant Number 25105010.

## References

- [1] D. F. Geesaman, K. Saito, and A. W. Thomas, Ann. Rev. Nucl. Part. Sci. 45 , 337 (1995).
- [2] M. Hirai, S. Kumano, and M. Miyama, Phys. Rev. D 64 , 034003 (2001); M. Hirai, S. Kumano, and T.-H. Nagai, Phys. Rev. C 70 , 044905 (2004); 76 , 065207 (2007).
- [3] J. Seely et al. , Phys. Rev. Lett. 103 , 202301 (2009).
- [4] M. Hirai, S. Kumano, K. Saito, and T. Watanabe, Phys. Rev. C 83 , 035202 (2011).
- [5] M. Ericson and S. Kumano, Phys. Rev. C 67 , 022201 (2003).
- [6] Y. Kanada-En'yo, M. Kimura, and H. Horiuchi, C. R. Physique 4 , 497 (2003).
- [7] H. Feldmeier and J. Schnack, Rev. Mod. Phys. 72 , 655 (2000).
- [8] Jefferson Lab PAC-35 proposal, PR12-10-008, J. Arrington et al. (2009).
M. Hirai S. Kumano K. Saito T. Watanabe
Clustering structure of nuclei in deep inelastic processes
A clustering aspect is explained for the $^9$Be nucleus in charged-lepton
deep inelastic scattering. Nuclear modifications of the structure function
$F_2$ are studied by the ratio $R_{\rm EMC} = F_2^A /F_2^D$, where $A$ and $D$
are a nucleus and the deuteron, respectively. In a JLab experiment, an
unexpectedly large nuclear modification slope $|dR_{\rm EMC}/dx|$ was found for
$^9$Be, which could be related to its clustering structure. We investigated a
mean conventional part of a nuclear structure function $F_2^A$ by a convolution
description with nucleon momentum distributions calculated by antisymmetrized
(or fermionic) molecular dynamics (AMD) and also by a simple shell model. We
found that clustering effects are small in the conventional part, so that the
JLab result could be associated with an internal nucleon modification or a
short-range nuclear correlation which is caused by high densities due to
cluster formation.