The original JETSET program is intimately connected with string fragmentation, in the form of the time-honoured `Lund model'. This is the default for all PYTHIA applications, but independent fragmentation options also exist (although not actively maintained), for applications where one wishes to study the importance of string effects.

All current models are of a probabilistic and iterative nature. This means that the fragmentation process as a whole is described in terms of one or a few simple underlying branchings, of the type jet hadron + remainder-jet, string hadron + remainder-string, and so on. At each branching, probabilistic rules are given for the production of new flavours, and for the sharing of energy and momentum between the products.

To understand fragmentation models, it is useful to start with the simplest possible system, a colour-singlet 2-jet event, as produced in annihilation. Here lattice QCD studies lend support to a linear confinement picture (in the absence of dynamical quarks), i.e. the energy stored in the colour dipole field between a charge and an anticharge increases linearly with the separation between the charges, if the short-distance Coulomb term is neglected. This is quite different from the behaviour in QED, and is related to the presence of a triple-gluon vertex in QCD. The details are not yet well understood, however.

The assumption of linear confinement provides the starting point for the string model. As the and partons move apart from their common production vertex, the physical picture is that of a colour flux tube (or maybe a colour vortex line) being stretched between the and the . The transverse dimensions of the tube are of typical hadronic sizes, roughly 1 fm. If the tube is assumed to be uniform along its length, this automatically leads to a confinement picture with a linearly rising potential. In order to obtain a Lorentz covariant and causal description of the energy flow due to this linear confinement, the most straightforward way is to use the dynamics of the massless relativistic string with no transverse degrees of freedom. The mathematical, one-dimensional string can be thought of as parameterizing the position of the axis of a cylindrically symmetric flux tube. From hadron spectroscopy, the string constant, i.e. the amount of energy per unit length, is deduced to be GeV/fm. The expression `massless' relativistic string is somewhat of a misnomer: effectively corresponds to a `mass density' along the string.

Let us now turn to the fragmentation process. As the and move apart, the potential energy stored in the string increases, and the string may break by the production of a new pair, so that the system splits into two colour-singlet systems and . If the invariant mass of either of these string pieces is large enough, further breaks may occur. In the Lund string model, the string break-up process is assumed to proceed until only on-mass-shell hadrons remain, each hadron corresponding to a small piece of string with a quark in one end and an antiquark in the other.

In order to generate the quark-antiquark pairs which lead to string break-ups, the Lund model invokes the idea of quantum mechanical tunnelling. This leads to a flavour-independent Gaussian spectrum for the of pairs. Since the string is assumed to have no transverse excitations, this is locally compensated between the quark and the antiquark of the pair. The total of a hadron is made up out of the contributions from the quark and antiquark that together form the hadron. Some contribution of very soft perturbative gluon emission may also effectively be included in this description.

The tunnelling picture also implies a suppression of heavy-quark production, . Charm and heavier quarks hence are not expected to be produced in the soft fragmentation, but only in perturbative parton-shower branchings .

When the quark and antiquark from two adjacent string breaks are combined to form a meson, it is necessary to invoke an algorithm to choose between the different allowed possibilities, notably between pseudoscalar and vector mesons. Here the string model is not particularly predictive. Qualitatively one expects a ratio, from counting the number of spin states, multiplied by some wave-function normalization factor, which should disfavour heavier states.

A tunnelling mechanism can also be used to explain the production of baryons. This is still a poorly understood area. In the simplest possible approach, a diquark in a colour antitriplet state is just treated like an ordinary antiquark, such that a string can break either by quark-antiquark or antidiquark-diquark pair production. A more complex scenario is the `popcorn' one, where diquarks as such do not exist, but rather quark-antiquark pairs are produced one after the other. This latter picture gives a less strong correlation in flavour and momentum space between the baryon and the antibaryon of a pair.

In general, the different string breaks are causally disconnected. This means that it is possible to describe the breaks in any convenient order, e.g. from the quark end inwards. One therefore is led to write down an iterative scheme for the fragmentation, as follows. Assume an initial quark moving out along the axis, with the antiquark going out in the opposite direction. By the production of a pair, a meson with flavour content is produced, leaving behind an unpaired quark . A second pair may now be produced, to give a new meson with flavours , etc. At each step the produced hadron takes some fraction of the available energy and momentum. This process may be iterated until all energy is used up, with some modifications close to the end of the string in order to make total energy and momentum come out right.

The choice of starting the fragmentation from the quark end is arbitrary, however. A fragmentation process described in terms of starting at the end of the system and fragmenting towards the end should be equivalent. This `left-right' symmetry constrains the allowed shape of the fragmentation function , where is the fraction of the remaining light-cone momentum (+ for the jet, for the one) taken by each new particle. The resulting `Lund symmetric fragmentation function' has two free parameters, which are determined from data.

If several partons are moving apart from a common origin, the details of the string drawing become more complicated. For a event, a string is stretched from the end via the to the end, i.e. the gluon is a kink on the string, carrying energy and momentum. As a consequence, the gluon has two string pieces attached, and the ratio of gluon to quark string force is 2, a number which can be compared with the ratio of colour charge Casimir operators, . In this, as in other respects, the string model can be viewed as a variant of QCD where the number of colours is not 3 but infinite. Note that the factor 2 above does not depend on the kinematical configuration: a smaller opening angle between two partons corresponds to a smaller string length drawn out per unit time, but also to an increased transverse velocity of the string piece, which gives an exactly compensating boost factor in the energy density per unit string length.

The string will fragment along its length. To first approximation this means that there is one fragmenting string piece between and and a second one between and . One hadron is straddling both string pieces, i.e. sitting around the gluon corner. The rest of the particles are produced as in two simple strings, but strings boosted with respect to the overall c.m. frame. When considered in detail, the string motion and fragmentation is more complicated, with the appearance of additional string regions during the time evolution of the system. These corrections are especially important for soft and collinear gluons, since they provide a smooth transition between events where such radiation took place and events where it did not. Therefore the string fragmentation scheme is `infrared safe' with respect to soft or collinear gluon emission.

Another possible colour topology arises when considering baryon-number-violating processes, or events where more than one valence quark has been knocked out of a beam baryon (as can happen when multiple parton-parton interactions occur). In this case, there will be three (anti-)colour carriers connected antisymmetrically in colour, and of which no two may naturally be considered to form a diquark system. The string topology will thus not be of the simple type, but rather a `Y' shaped topology is spanned between the endpoints. The vertex of the `Y' topology comes to be of special interest in the fragmentation, and will be referred to as a `string junction'. Each of the three string pieces undergo a fragmentation process subject to exactly the same principles as outlined above, only a baryon containing the junction will eventually be formed. The picture is essentially that of three jets going out, with the junction baryon formed `in the middle', hence the junction baryon will tend to have a soft spectrum when the jets are widely separated. Note that, in the limit that two of the endpoints of the 'Y' come close together, the diquark picture for beam remnants mentioned above is effectively recovered, with only minor differences remaining.

For events that involve many partons, there may be several possible topologies for their ordering along the string. An example would be a (the gluon indices are here used to label two different gluon-momentum vectors), where the string can connect the partons in either of the sequences and . The matrix elements that are calculable in perturbation theory contain interference terms between these two possibilities, which means that the colour flow is not always well-defined. Fortunately, the interference terms are down in magnitude by a factor , where is the number of colours, so approximate recipes can be found. In the leading log shower description, on the other hand, the rules for the colour flow are well-defined.

A final comment: in the argumentation for the importance of colour flows there is a tacit assumption that soft-gluon exchanges between partons will not normally mess up the original colour assignment. Colour rearrangement models provide toy scenarios wherein deviations from this rule could be studied. Of particular interest has been the process , where the original singlets and could be rearranged to and . So far, there are no experimental evidence for dramatic effects of this kind, but the more realistic models predict effects sufficiently small that these have not been ruled out. Another example of nontrivial effects is that of Bose-Einstein correlations between identical final-state particles, which reflect the true quantum nature of the hadronization process.