# Multiple-particle interaction in dimensional lattice model

###### Abstract

Finite volume multiple-particle interaction is studied in a two-dimensional complex lattice model. The existence of analytical solutions to the model in two-dimensional space and time makes it a perfect model for the numerical study of finite volume effects of multi-particle interaction. The spectra from multiple particles are extracted from the Monte Carlo simulation on various lattices in several moving frames. The -matrix of multi-particle scattering in theory is completely determined by two fundamental parameters: single particle mass and the coupling strength of two-to-two particle interaction. These two parameters are fixed by studying single-particle and two-particle spectra. Due to the absence of the diffraction effect in the model, three-particle quantization conditions are given in a simple analytical form. The three-particle spectra from simulation show remarkable agreement with the prediction of exact solutions.

## I Introduction

One of the outstanding but challenging goals in nuclear/hadron physics is to understand the dynamics of particle interaction. Multiple particle interaction is not only important to nuclear/hadron physics, but also plays a crucial role in astrophysics, atomic and condensed matter physics. However, the complication increases dramatically with increasing numbers of dynamical degrees of freedom and poses a significant obstacle in studying and understanding multi-particle interaction. Fortunately, the simplest case of multi-particle interaction turns out to be manageable, three-particle interaction. The dynamics of three-particle interaction were well developed and studied in the past Taylor:1966zza (); Basdevant:1966zzb (); Gross:1982ny (); Faddeev:1960su (); Faddeev:1965 (); Phillips:1966zza (); Gloeckle:1983 (); Fedorov:1993 (); Gloeckle:1995jg (); Khuri:1960zz (); Bronzan:1963xn (); Aitchison:1965kt (); Aitchison:1965zz (); Aitchison:1966kt (); Pasquier:1968zz (); Pasquier:1969dt (); Guo:2014vya (); Guo:2014mpp (); Danilkin:2014cra (); Guo:2015kla (). Recent progress in high statistic experiments, such as GlueX and CLAS programs, have triggered renewed interests in three-body dynamics. One example is the extraction of - and -quark mass difference from decay process Kambor:1995yc (); Anisovich:1996tx (); Colangelo:2009db (); Lanz:2013ku (); Schneider:2010hs (); Kampf:2011wr (); Guo:2015zqa (); Guo:2016wsi (). On the other hand, lattice QCD provides an unprecedented opportunity for the study of multiple particle interaction from the heart of hadrons with quarks and gluons as the fundamental building blocks. Recent advances in lattice computation have made the study of hadron interaction especially possible Aoki:2007rd (); Sasaki:2008sv (); Feng:2010es (); Dudek:2010ew (); Beane:2011sc (); Lang:2011mn (); Aoki:2011yj (); Dudek:2012gj (); Dudek:2012xn (); Wilson:2014cna (); Wilson:2015dqa (); Dudek:2016cru (). Because lattice QCD is formulated in Euclidean space, access to scattering information is not always direct. That adds some additional complication in multi-particle studies in lattice QCD as well as the intense numerical computation and other difficulties. A formalism was proposed nearly 30 years ago by Lüscher Lusher:1991 () to tackle the two-particle elastic scattering problem in finite volume; it is known as the Lüscher formula. Since then, the framework quickly extended to moving frames Gottlieb:1995 (); Lin:2001 (); Christ:2005 (); Bernard:2007 (); Bernard:2008 (), and to coupled-channel scattering Liu:2005 (); Lage:2009 (); Doring:2011 (); Aoki:2011 (); Briceno:2012yi (); Hansen:2012tf (); Guo:2013cp (). In the three-particle sector, many groups have made remarkable progress Kreuzer:2008bi (); Kreuzer:2009jp (); Kreuzer:2012sr (); Polejaeva:2012ut (); Briceno:2012rv (); Hansen:2014eka (); Hansen:2015zga (); Hansen:2016fzj (); Hammer:2017uqm (); Hammer:2017kms (); Guo:2016fgl (); Guo:2017ism (); Meissner:2014dea (); Briceno:2017tce (); Sharpe:2017jej (); Mai:2017bge (); Guo:2017crd () related to the theoretical algorithm of extracting scattering amplitudes from lattice data in recent years.

A three-particle lattice simulation was recently performed based on a complex toy model Romero-Lopez:2018rcb (), the data analysis was carried out by adopting effective theory framework. However, the simulation and analysis are limited solely to ground state energy levels where all three particles are nearly at rest, and the three-particle signals are quite noisy. In present work, we aim to perform a simulation on multiple-particle interaction also using model, and study the finite volume effect on multiple-particle spectra in a better controlled environment and a more systematic way. For this purpose, multiple numbers of multi-particle operators are used in our simulation and variational analysis Michael:1985ne (); Luscher:1990ck (); Blossier:2009kd () is implemented to extract excited state energy levels. The exact scattering solutions of theory in dimensions are known in both free space Thacker:1974kv (); McGuire:1964zt (); Yang:1967bm () and finite volume Guo:2016fgl (). Taking advantage of existing analytic multiple-particle scattering solutions, the simulation is therefore performed in dimensional space and time for various lattice sizes and moving frames. The exact scattering solutions are used in data analysis of multi-particle simulation. In principle, the multiple particle scattering -matrices are completely determined by only two free parameters: single particle mass and coupling strength of two-to-two particle interaction. The single particle mass is obtained from single particle correlation functions, and the coupling strength of pair-wise interaction is extracted by studying two-particle scattering spectra in a lattice. The comparison between three-particle scattering spectra and predicted three-particle energies by using analytic expression of three-particle quantization conditions are presented in the end.

The paper is organized as follows. The exact solutions of theory for two-body and three-body interaction are summarized in Section II. The algorithm of the Hybrid Monte Carlo simulation of lattice model and strategy of data analysis are briefly discussed in Section III. The construction of multi-particle operators, multi-particle spectra in lattice simulation and data analysis are described in Section IV. The summary and outlook are given in Section V.

## Ii Exact solution of model in

In this section, we summarize some results of the two-dimensional model. Classical action of the complex model in two-dimensional Euclidean space is

(1) |

where are temporal and spatial coordinates in two-dimensional Euclidean space, respectively. It is known Thacker:1974kv () that the complex model in Eq.(1) is equivalent to a non-relativistic one-dimensional -body interaction problem of particles interacting with pair-wise -function potentials,

(2) |

where refers to the spatial position of -th particle, and stands for the mass of identical bosons. The coupling strength of -function potential, , differs from renormalized in Eq.(1) by a constant factor. The exact solutions of multi-particle interaction with -function potentials were studied and obtained in both free space Thacker:1974kv (); McGuire:1964zt (); Yang:1967bm () and finite volume Guo:2016fgl (). In fact, the particles interacting with -function potential in is only one of few exactly solvable multi-particle scattering problems. The multi-particle wave function is described completely by the linear superpositions of plane waves with all possible permutation on particle momenta. No new momenta are generated by collision, all the diffraction effects are canceled out as the consequence of Bethe’s hypothesis Bethe:1931hc (); Lieb:1963rt (). The multi-particle -matrix therefore is factorized into the product of a number of two-particle scattering amplitudes, as if the process of multi-particle scattering would be a succession of separated elastic two-particle collisions Guo:2016fgl ().

In finite volume for two-particle scattering, only one quantization condition is required Guo:2016fgl (); Guo:2013vsa ()

(3) |

where and denote center of mass and relative momenta of two particles. The phase shift for -function potential is given by . The stands for the size of the square box in , and center of mass momentum is discretized because of the periodic boundary condition of lattice: .

For three-body scattering in finite volume, three quantization conditions are obtained Guo:2016fgl (). Only two of them are independent,

(4) |

where all the relative momenta are given in terms of two independent particle momenta: and , , , . The momentum of particle-3 is constrained by momentum conservation, . Again, center of mass momentum of three-particle is quantized in the periodic box: .

## Iii The lattice model action

The lattice action is obtained from Eq.(1) by replacing the continuous derivative with discrete difference: , where denotes the unit vector in direction on a periodic square lattice. In addition, by introducing two new parameters: and , and also rescaling the field by , we thus obtain

(5) |

where now refers to discrete coordinates of Euclidean lattice site.

### iii.1 Hybrid Monte Carlo algorithm

The Hybrid Monte Carlo algorithm Duane:1987de (); Duane:1986iw () is adopted in our numerical simulation, the complex model is treated as a coupled two component scalar field model, . In Hybrid Monte Carlo simulation Duane:1987de (); Duane:1986iw (), an auxiliary Hamiltonian is introduced

(6) |

where are fictitious conjugate momenta of field. The auxiliary Hamiltonian in Eq.(6) defines classical evolution of both and fields over a fictitious time within an interval :

(7) |

The trajectory of over time interval is determined by the solutions of motion equations in Eq.(7).

The two pairs of components, and , are updated alternately for each sweep over an entire lattice. Updating each pair is followed with the standard Hybrid Monte Carlo algorithm:

(i) the trajectory begins with choosing a random distribution of fields at initial time . The initial conjugate momenta, , are generated according to the Gaussian probability distribution: .

(ii) solve motion equations in Eq.(7) to evolve over the trajectory up to a time . The motion equations, Eq.(7), are solved numerically by the leapfrog method Duane:1986iw ().

(iii) accept the proposed new fields, , with probability: , where .

The simulations are performed with the choice of the parameters: , and . The temporal extent of the lattice is fixed at , and the spatial extent of lattice, , are from up to . For each set of lattice size and moving frame, one million measurements are generated. The length of trajectory is fixed at , the fields evolve from up to over discrete steps.

### iii.2 Strategy of data analysis

As already mentioned in Section II, the two-dimensional model is exactly solvable, the solutions of the model are given in terms of only two free parameters: particle mass, , and the coupling strength of -function potential, . The mass of identical particles, , can be extracted from one-particle spectra of the lattice simulation. The second parameter, , can be fixed by two-particle spectra from the simulation. Taking advantage of the existence of exact solutions of the two-dimensional model provides an excellent playground and controlled environment for a systematic study of finite volume effects of multi-particle scattering in lattice simulation. In present work, we are not aiming at obtaining any new fundamental information from three-body spectra, such as three-body force effect, etc. Instead, after fixing and from one- and two-particle spectra, we tend to study how well the three-body spectra from simulation match the prediction of exact solutions. In real QCD simulation, the significant difference between simulation results of three-body spectra and prediction based on pair-wise interaction may signal the effect of three-body forces or something more fundamental. The present work serves only as a testbed for more realistic future lattice studies of multi-particle interaction.

To accomplish the goal of this work mentioned above, the following steps are taken in data analysis of the simulation results:

1. measure one-particle spectra for various sizes of lattice and moving frames, and extract continuum limit particle mass, , by using relation Gatteringer:1993 (): .

2. measure two-particle spectra for various sizes of lattice and moving frames, and extract the coupling strength, , from lattice data.

3. three-particle spectra are measured for various sizes of lattice and moving frames as well, three-particle spectra are thus compared with predicted three-particle spectra. The predicted three-particle spectra, , are given in terms of two independent particle momenta, and , by

where and are the solutions of Eq.(4), and .

The particles spectra are extracted by fitting exponential multi-particle correlation functions as a function of time : . See the example of one-, two- and three-particle correlation functions and effective mass, , in Fig.1 and Fig.2, respectively. The construction of multi-particle operators and correlation functions will be explained later in Section IV.

## Iv Particles spectra and data analysis

In this section, we present significant results for multi-particle scattering. Some details on multi-particle operator construction and data analysis are also given.

### iv.1 One particle spectra

The one-particle spectra are extracted from the exponential decay of the correlation functions

(8) |

where the one particle propagator, , is defined by

(9) |

Single particle energy is obtained for multiple lattice sizes, up to . By fitting single particle energies in multiple lattice sizes with relation

(10) |

where and are used as fitting parameters, we thus find the mass of single particle: , see Fig.3. The excited single particle energy levels are used to check the energy-momentum dispersion relations in a finite lattice,

(11) |

The comparison between lattice results and the lattice dispersion relation is presented in Fig.4.

### iv.2 Two particles spectra

In moving frames, the matrix element of the two particle correlation function read

(12) |

where is related to center of mass momentum by , and two-particle operators are constructed by

(13) |

Four two-particle operators are used in our simulation: , so the size of matrix of two-particle correlation functions are , and for .

The spectral decomposition of the correlation function matrices are usually given by

(14) |

where , and labels the -th energy eigenstate . In order to extract excited energy states, a generalized eigenvalue method Luscher:1990ck () is proposed

(15) |

where is a small reference time. Mixing of multi-particle states is protected by the conservation of charge quantum number in the complex model. Also, diagonalized correlation functions barely show the contamination of higher energy states in , see Fig.1 and Fig.2. Therefore, is set to zero and a simple form of is used in the data fitting for . The two particle spectra for various lattice sizes and are presented in Fig.5.

The phase shift of two-body scattering is extracted from two-particle energy levels, , by using relation:

(16) |

where the relative momentum of two particles, , is given by the solutions of two-particle energy momentum dispersion relation

(17) |

where , , see extracted phase shift in Fig.6. The exact expression of phase shift,

(18) |

is used to fit lattice results, , and to fix the coupling strength, . We thus find .

### iv.3 Three particles spectra

For three-particle operators with , four operators are used in present work:

(19) |

and

(20) |

Similar to the two-particle correlation function matrix, the matrix element of the three particle correlation function is given by

(21) |

In the two-particle sector, the generalized eigenvalue method is also applied to extract three-body energy levels,

(22) |

where . An example of the three-particle correlation function, , and effective mass, is given in Fig.1 and Fig.2.

Given the values of particle mass, , and coupling strength, , that we learned from discussion in previous sections, three-particle spectra does not provide any new insight into the fundamental parameters of theory due to the absence of three-body force. However, in general, three-particle spectra are still considered a useful tool to explore and understand the dynamics of three-particle interaction. In reality, it also provides opportunities to investigate the possibility of more fundamental parameters of lattice QCD theory. Nevertheles, since the exact solutions are known, we only tend to demonstrate the consistence of predicted three-particle spectra compared to simulation results. The predicted three-particle spectra are determined by three-body energy-momentum dispersion relations in terms of two independent particle momenta, say and ,

(23) |

where . Two independent particle momenta, and , are the solutions of three-body quantization conditions given in Eq.(4). The main results of three-particle spectra are presented in Fig.7. As we can see, the agreement of predicted spectra against simulation results is quite remarkable.

## V Summary

In summary, the lattice simulation of multi-particle interaction is studied by using a complex lattice model. The simulation is performed in dimensional space and time for various sizes of lattices and multiple moving frames. The two dimensional model is exactly solvable and analytical expressions of multi-particle quantization conditions are known in finite volume Guo:2016fgl (). This feature makes it a perfect testbed for studying multi-particle interaction in a lattice. The typical numbers of multi-particle operators are used in our simulation, and a variational approach is implemented to extract excited state energy levels. Two parameters of theory, single particle mass and coupling strength, are extracted from single particle and two particles spectra, respectively. Then, extracted theory parameters are applied to predict three-particle spectra by using analytical three-body quantization conditions compared with three-particle spectra from simulations. Predicted three-particle spectra and lattice results show quite remarkable agreement.

## Vi Acknowledgments

We acknowledge support from the Department of Physics and Engineering, California State University, Bakersfield, CA. We also thank Vladimir Gasparian and David Gross for numerous fruitful discussions. This research was supported in part by the National Science Foundation under Grant No. NSF PHY-1748958.

## References

- (1) J. G. Taylor, Phys. Rev. 150, 1321 (1966).
- (2) J. -L. Basdevant and R. E. Kreps, Phys. Rev. 141, 1398 (1966).
- (3) F. Gross, Phys. Rev. C 26, 2226 (1982).
- (4) L. D. Faddeev, Zh. Eksp. Teor. Fiz. 39, 1459 (1960) [Sov. Phys.-JETP 12, 1014(1961)].
- (5) L. D. Faddeev, Mathematical Aspects of the Three-Body Problem in the Quantum Scattering Theory, Israel Program for Scientific Translation, Jerusalem, Israel (1965).
- (6) W. Glöckle, The Quantum Mechanical Few-Body Problem, Springer, Berlin, Germany (1983).
- (7) A. C. Phillips, Phys. Rev. 142, 984 (1966).
- (8) D. V. Fedorov and A S. Jensen, Phys. Rev. Lett. 71, 4103 (1993).
- (9) W. Glöckle, H. Witala, D. Hüber, H. Kamada and J. Golak, Phys. Rept. 274, 107 (1996).
- (10) N. N. Khuri and S. B. Treiman, Phys. Rev. 119, 1115 (1960).
- (11) J. B. Bronzan and C. Kacser, Phys. Rev. 132, 2703 (1963).
- (12) I. J. R. Aitchison, II Nuovo Cimento 35, 434 (1965).
- (13) I. J. R. Aitchison, Phys. Rev. 137, B1070 (1965); Phys. Rev. 154, 1622 (1967).
- (14) I. J. R. Aitchison and R. Pasquier, Phys. Rev. 152, 1274 (1966).
- (15) R. Pasquier and J. Y. Pasquier, Phys. Rev. 170, 1294 (1968).
- (16) R. Pasquier and J. Y. Pasquier, Phys. Rev. 177, 2482 (1969).
- (17) P. Guo, I. V. Danilkin and A. P. Szczepaniak, Eur. Phys. J. A 51, 135 (2015).
- (18) P. Guo, Phys. Rev. D 91, 076012 (2015).
- (19) I. V. Danilkin, C. Fernández-Ramírez, P. Guo, V. Mathieu, D. Schott and A. P. Szczepaniak, Phys. Rev. D 91, 094029 (2015).
- (20) P. Guo, Mod. Phys. Lett. A 31, 1650058 (2016).
- (21) J. Kambor, C. Wiesendanger, and D. Wyler, Nucl. Phys. B465, 215 (1996).
- (22) A. V. Anisovich and H. Leutwyler, Phys. Lett. B375, 335 (1996).
- (23) G. Colangelo, S. Lanz, and E. Passemar, PoS CD09, 047 (2009).
- (24) S. Lanz, PoS CD12, 007 (2013).
- (25) S. P. Schneider, B. Kubis, and C. Ditsche, JHEP 1102, 028 (2011).
- (26) K. Kampf, M. Knecht, J. Novotny, and M. Zdrahal, Phys. Rev. D84, 114015 (2011).
- (27) P. Guo, Igor V. Danilkin, D. Schott, C. Fernández-Ramírez, V. Mathieu and A. P. Szczepaniak, Phys. Rev. D92, 054016 (2015).
- (28) P. Guo, Igor V. Danilkin, C. Fernández-Ramírez, V. Mathieu and A. P. Szczepaniak, Phys. Lett. B771, 497 (2017).
- (29) S. Aoki et al. [CP-PACS Collaboration], Phys. Rev. D 76, 094506 (2007)
- (30) K. Sasaki, and N. Ishizuka, Phys. Rev. D 78, 014511 (2008).
- (31) X. Feng, K. Jansen, and D. B. Renner, Phys. Rev. D 83, 094505 (2011).
- (32) J. J. Dudek et al. (Hadron Spectrum Collaboration), Phys. Rev. D 83, 071504 (2011).
- (33) S. R. Beane et al. (NPLQCD Collaboration), Phys. Rev. D 85, 034505 (2012).
- (34) C. B. Lang, D. Mohler, S. Prelovsek and M. Vidmar, Phys. Rev. D 84, 054503 (2011).
- (35) S. Aoki et al. [CS Collaboration], Phys. Rev. D 84, 094505 (2011)
- (36) J. J. Dudek et al. (Hadron Spectrum Collaboration), Phys. Rev. D 86, 034031 (2012).
- (37) J. J. Dudek, R. G. Edwards and C. E. Thomas, Phys. Rev. D 87, 034505 (2013).
- (38) D. J. Wilson, J. J. Dudek, R. G. Edwards and C. E. Thomas, Phys. Rev. D 91, 054008 (2015).
- (39) D. J. Wilson, R. A. Briceno, J. J. Dudek, R. G. Edwards and C. E. Thomas, Phys. Rev. D 92, 094502 (2015).
- (40) J. J. Dudek, et al. (Hadron Spectrum Collaboration), Phys. Rev. D 93, 094506 (2016).
- (41) M. Lüscher, Nucl. Phys. B 354, 531 (1991).
- (42) K. Rummukainen, S. Gottlieb, Nucl. Phys. B 450, 397 (1995).
- (43) C.-J.D. Lin, G. Martinelli, C. T. Sachrajda and M. Testa, Nucl. Phys. B 619, 467 (2001).
- (44) N. H. Christ, C. Kim and T.Yamazaki, Phys. Rev. D 72, 114506 (2005).
- (45) V. Bernard, Ulf-G. Meißner and A. Rusetsky, Nucl. Phys. B 788, 1 (2008).
- (46) V. Bernard, M. Lage, Ulf-G. Meißner and A.Rusetsky, JHEP 0808, 024 (2008).
- (47) S. He, X. Feng, C. Liu, JHEP 0507, 011 (2005).
- (48) M. Lage, Ulf-G. Meißner and A. Rusetsky, Phys. Lett. B 681, 439 (2009)
- (49) M. Döring, Ulf-G. Meißner,E. Oset and A. Rusetsky, Eur. Phys. J. A 47, 139 (2011)
- (50) S. Aoki et al. [HAL QCD Collaboration], Proc. Japan Acad. B 87, 509 (2011)
- (51) R. A. Briceno and Z. Davoudi, arXiv:1204.1110 [hep-lat].
- (52) M. T. Hansen and S. R. Sharpe, Phys. Rev. D 86, 016007 (2012)
- (53) P. Guo, J. Dudek, R. Edwards, A. P. Szczepaniak, [arXiv:1211.0929 [hep-lat]].
- (54) S. Kreuzer and H.-W. Hammer, Phys. Lett. B 673, 260 (2009).
- (55) S. Kreuzer and H.-W. Hammer, Eur. Phys. J. A 43, 229 (2010).
- (56) S. Kreuzer and H.-W. Hammer, Eur. Phys. J. A 48, 93 (2012).
- (57) K. Polejaeva and A. Rusetsky, Eur. Phys. J. A 48, 67 (2012).
- (58) R. A. Briceno and Z. Davoudi, Phys. Rev. D 87, 094507 (2013).
- (59) M. T. Hansen and S. R. Sharpe, Phys. Rev. D 90, 116003 (2014).
- (60) M. T. Hansen and S. R. Sharpe, Phys. Rev. D 92, 114509 (2015).
- (61) M. T. Hansen and S. R. Sharpe, Phys. Rev. D 93, 096006 (2016).
- (62) H. -W. Hammer, J. -Y. Pang and A. Rusetsky, JHEP 1709, 109 (2017).
- (63) H. -W. Hammer, J. -Y. Pang and A. Rusetsky, JHEP 1710, 115 (2017).
- (64) P. Guo, Phys. Rev. D95, 054508 (2017).
- (65) P. Guo and V. Gasparian, Phys. Lett. B774, 441 (2017).
- (66) Ulf-G. Meißner, G. Rios and A. Rusetsky, Phys. Rev. Lett. 114, 091602 (2015). Erratum: Phys. Rev. Lett. 117, 069902 (2016).
- (67) R. A. Briceno, M. T. Hansen and S. R. Sharpe, Phys. Rev. D95, 074510 (2017).
- (68) S. R. Sharpe, Phys. Rev. D96, 054515 (2017).
- (69) M. Mai and M. Döring, Euro. Phys. J A53, 240 (2017).
- (70) P. Guo and V. Gasparian, Phys. Rev. D97, 014504 (2018).
- (71) F. Romero-Löpez, A. Rusetsky and C. Urbach, arXiv:1806.02367 [hep-lat].
- (72) C. Michael, Nucl. Phys. B 259, 58 (1985).
- (73) M. Luscher and U. Wolff, Nucl. Phys. B 339, 222 (1990).
- (74) B. Blossier, M. Della Morte, G. von Hippel, T. Mendes and R. Sommer, JHEP 0904, 094 (2009).
- (75) H. B. Thacker, Phys. Rev. D 11, 838 (1975).
- (76) J. B. McGuire, J. Math. Phys. 5, 622 (1964).
- (77) C. N. Yang, Phys. Rev. Lett. 19, 1312(1967).
- (78) H. A. Bethe, Z. Phys. 71, 205(1931).
- (79) E. H. Lieb and W. Liniger, Phys. Rev. 130, 1605(1963).
- (80) P. Guo, Phys. Rev. D88, 014507 (2013).
- (81) S. Duane, A. D. Kennedy, B. J. Pendleton and D. Roweth, Phys. Lett. B195, 216 (1987).
- (82) S. Duane and J. B. Kogut, Nucl. Phys. B275, 398 (1986).
- (83) C. R. Gattringer and C. B. Lang, Nucl. Phys. B 391, 463 (1993).