Monte Carlo calculational methods for the generation of events with Bose-Einstein correlations

The momentum-space correlations for n bosons emitted by a source distributed in space and time are discussed in the context of intensity interferometry for identical pions. The Metropolis algorithm is used to generate events containing such correlations to all orders via a Monte Carlo technique. Direct calculation of the probability for multiparticle correlations is practicable for n-body states up to n∼20. Beyond this, a method which samples the symmetrized n-body probability must be used. It is observed that the traditional pair-correlation function is distorted for events with high phase-space density in a fashion consistent with the results of a simple model calculation.