Abstract. Physics and programming aspects are discussed for a Fortran 77 Monte Carlo program to simulate complete events in deep inelastic leptonnucleon scattering. The parton level interaction is based on the standard model electroweak cross sections, which are fully implemented in leading order for any lepton of arbitrary polarization, and different parametrizations of parton density functions can be used. First order QCD matrix elements for gluon radiation and boson-gluon fusion are implemented and higher order QCD radiation is treated using parton showers. Hadronization is performed using the Lund string model, implemented in Jetset/Pythia. Rapidity gap events are generated through a model based on soft colour interactions.

We introduce soft colour interactions as a novel mechanism to understand the observed events with large rapidity gaps in ep collisions at HERA. Colour exchanges between produced partons and colour-charges in the proton remnant modifies the colour structure for hadronization, such that colour singlet systems may appear well separated in rapidity. Our explicit model show characteristics of diffractive scattering, although no explicit pomeron dynamics have been introduced, and for non-gap events an increased forward energy flow gives agreement with data.

We discuss our model for soft colour interactions and present results on the diffractive structure function F D 2 (β, x IP , Q 2 ) and inclusive transverse energy flows, which agree with available HERA data.

The direct charmonium and bottomonium production rate observed at high p Ќ in p p collisions at the Fermilab Tevatron is factors of 10 larger than predictions based on conventional perturbative QCD. We show that this excess can be accounted for by our model for soft color interactions, previously introduced to describe in a novel way the large rapidity gap events observed at DESY HERA. ͓S0556-2821͑97͒06623-X͔

A simple model is presented for the parton distributions in hadrons. The parton momenta in the hadron rest frame are derived from a spherically symmetric, Gaussian, distribution having a width motivated by the Heisenberg uncertainty relation applied to the hadron size. Valence quarks and gluons originate from the `bare' hadron, while sea partons arise mainly from pions in hadronic fluctuations. Starting from a low Q^2 scale, the distributions are evolved with next-to-leading order DGLAP and give the proton structure function F2(x,Q^2) in good agreement with deep inelastic scattering data.Comment: 12 pages, Latex, 5 figures. Added discussion on the applicability of the model at small Q^2. Figure 3 modified to show smaller Q^