Quantum Monte Carlo calculations of six-quark states
  The variational Monte Carlo method is used to find the ground state of six quarks confined to a cavity of diameter R_c, interacting via an assumed non-relativistic constituent quark model (CQM) Hamiltonian. We use a flux-tube model augmented with one-gluon and one-pion exchange interactions, which has been successful in describing single hadron spectra. The variational wave function is written as a product of three-quark nucleon states with correlations between quarks in different nucleons. We study the role of quark exchange effects by allowing flux-tube configuration mixing. An accurate six-body variational wave function is obtained. It has only ~13% rms fluctuation in the total energy and yields a standard deviation of ~<.1%; small enough to be useful in discerning nuclear interaction effects from the large rest mass of the two nucleons. Results are presented for three values of the cavity diameter, R_c=2, 4, and 6 fm. They indicate that the flux-tube model Hamiltonian with gluon and pion exchange requires revisions in order to obtain agreement with the energies estimated from realistic two-nucleon interactions. We calculate the two-quark probability distribution functions and show how they may be used to study and adjust the model Hamiltonian. 
Mark W. Paris, Vijay R. Pandharipande