Inclusive dijet and trijet production in deep inelastic $ep$ scattering has been measured for $10<Q^2<100$ GeV$^2$ and low Bjorken $x$, $10^{-4}<x_{\rm Bj}<10^{-2}$. The data were taken at the HERA $ep$ collider with centre-of-mass energy $\sqrt{s} = 318 \gev$ using the ZEUS detector and correspond to an integrated luminosity of $82 {\rm pb}^{-1}$. Jets were identified in the hadronic centre-of-mass (HCM) frame using the $k_{T}$ cluster algorithm in the longitudinally invariant inclusive mode. Measurements of dijet and trijet differential cross sections are presented as functions of $Q^2$, $x_{\rm Bj}$, jet transverse energy, and jet pseudorapidity. As a further examination of low-$x_{\rm Bj}$ dynamics, multi-differential cross sections as functions of the jet correlations in transverse momenta, azimuthal angles, and pseudorapidity are also presented. Calculations at $\mathcal{O}(\alpha_{s}^3)$ generally describe the trijet data well and improve the description of the dijet data compared to the calculation at $\mathcal{O}(\alpha_{s}^2)$.

Multijet production rates in neutral current deep inelastic scattering have been measured in the range of exchanged boson virtualities 10 < Q2 < 5000 GeV2. The data were taken at the ep collider HERA with centre-of-mass energy sqrt(s) = 318 GeV using the ZEUS detector and correspond to an integrated luminosity of 82.2 pb-1. Jets were identified in the Breit frame using the k_T cluster algorithm in the longitudinally invariant inclusive mode. Measurements of differential dijet and trijet cross sections are presented as functions of jet transverse energy E_{T,B}{jet}, pseudorapidity eta_{LAB}{jet} and Q2 with E_{T,B}{jet} > 5 GeV and -1 < eta_{LAB}{jet} < 2.5. Next-to-leading-order QCD calculations describe the data well. The value of the strong coupling constant alpha_s(M_Z), determined from the ratio of the trijet to dijet cross sections, is alpha_s(M_Z) = 0.1179 pm 0.0013(stat.) {+0.0028}_{-0.0046}(exp.) {+0.0064}_{-0.0046}(th.)

Measurements are reported of the production of dijet events with a leading neutron in ep interactions at HERA. Differential cross sections for photoproduction and deep inelastic scattering are presented as a function of several kinematic variables. Leading order QCD simulation programs are compared with the measurements. Models in which the real or virtual photon interacts with a parton of an exchanged pion are able to describe the data. Next-to-leading order perturbative QCD calculations based on pion exchange are found to be in good agreement with the measured cross sections. The fraction of leading neutron dijet events with respect to all dijet events is also determined. The dijet events with a leading neutron have a lower fraction of resolved photon processes than do the inclusive dijet data.

The dependence of dijet production on the virtuality of the exchanged photon, Q^2, has been studied by measuring dijet cross sections in the range 0 < Q^2 < 2000 GeV^2 with the ZEUS detector at HERA using an integrated luminosity of 38.6 pb^-1. Dijet cross sections were measured for jets with transverse energy E_T^jet > 7.5 and 6.5 GeV and pseudorapidities in the photon-proton centre-of-mass frame in the range -3 < eta^jet <0. The variable xg^obs, a measure of the photon momentum entering the hard process, was used to enhance the sensitivity of the measurement to the photon structure. The Q^2 dependence of the ratio of low- to high-xg^obs events was measured. Next-to-leading-order QCD predictions were found to generally underestimate the low-xg^obs contribution relative to that at high xg^obs. Monte Carlo models based on leading-logarithmic parton-showers, using a partonic structure for the photon which falls smoothly with increasing Q^2, provide a qualitative description of the data.

Dijet production has been studied in neutral current deep inelastic e+p scattering for 470 < Q**2 < 20000 GeV**2 with the ZEUS detector at HERA using an integrated luminosity of 38.4 pb**{-1}. Dijet differential cross sections are presented in a kinematic region where both theoretical and experimental uncertainties are small. Next-to-leading-order (NLO) QCD calculations describe the measured differential cross sections well. A QCD analysis of the measured dijet fraction as a function of Q**2 allows both a precise determination of alpha_s(M_Z) and a test of the energy-scale dependence of the strong coupling constant. A detailed analysis provides an improved estimate of the uncertainties of the NLO QCD cross sections arising from the parton distribution functions of the proton. The value of alpha_s(M_Z), as determined from the QCD fit, is alpha_s(M_Z) = 0.1166 +- 0.0019 (stat.) {+ 0.0024}_{-0.0033} (exp.)} {+ 0.0057}_{- 0.0044} (th.).

Three-jet production is studied for the first time in deep-inelastic positron-proton scattering. The measurement carried out with the H1 detector at HERA covers a large range of four-momentum transfer squared 5 < Q^2 < 5000 GeV^2 and invariant three-jet masses 25 < M_(3jet) < 140 GeV. Jets are defined by the inclusive k_T algorithm in the Breit frame. The size of the three-jet cross section and the ratio of the three-jet to the dijet cross section R_(3/2) are described over the whole phase space by the predictions of perturbative QCD in next-to-leading order. The shapes of angular jet distributions deviate significantly from a uniform population of the available phase space but are well described by the QCD calculation.

A measurement is presented of dijet and 3-jet cross sections in low-|t| diffractive deep-inelastic scattering interactions of the type ep -> eXY, where the system X is separated by a large rapidity gap from a low-mass baryonic system Y. Data taken with the H1 detector at HERA, corresponding to an integrated luminosity of 18.0 pb^(-1), are used to measure hadron level single and double differential cross sections for 4<Q^2<80 GeV^2, x_pom<0.05 and p_(T,jet)>4 GeV. The energy flow not attributed to jets is also investigated. The measurements are consistent with a factorising diffractive exchange with trajectory intercept close to 1.2 and tightly constrain the dominating diffractive gluon distribution. Viewed in terms of the diffractive scattering of partonic fluctuations of the photon, the data require the dominance of qqbarg over qqbar states. Soft colour neutralisation models in their present form cannot simultaneously reproduce the shapes and the normalisations of the differential cross sections. Models based on 2-gluon exchange are able to reproduce the shapes of the cross sections at low x_pom values.

Jet production is studied in the Breit frame in deep-inelastic positron-proton scattering over a large range of four-momentum transfers 5 < Q^2 < 15000 GeV^2 and transverse jet energies 7 < E_T < 60 GeV. The analysis is based on data corresponding to an integrated luminosity of L_int \simeq 33 pb^(-1) taken in the years 1995-1997 with the H1 detector at HERA at a center-of-mass energy sqrt(s)=300 GeV. Dijet and inclusive jet cross sections are measured multi-differentially using k_perp and angular ordered jet algorithms. The results are compared to the predictions of perturbative QCD calculations in next-to-leading order in the strong coupling constant alphas.QCD fits are performed in which alphas and the gluon density in the proton are determined separately. The gluon density is found to be in good agreement with results obtained in other analyses using data from different processes. The strong coupling constant is determined to be alphas(MZ)=0.1186+-0.0059. In addition an analysis of the data in which both alphas and the gluon density are determined simultaneously is presented.

The production of isolated high-energy photons accompanied by jets has been measured in deep inelastic ep scattering with the ZEUS detector at HERA, using an integrated luminosity of 326 pb^{-1}. Measurements were made for exchanged photon virtualities, Q^2, in the range 10 to 350 GeV^2. The photons were measured in the transverse-energy and pseudorapidity ranges 4 < ET^gamma < 15 GeV and -0.7 < eta^gamma < 0.9, and the jets were measured in the transverse-energy and pseudorapidity ranges 2.5 < ET^jet <35 GeV and -1.5 < eta^jet < 1.8. Differential cross sections are presented as functions of these quantities. Perturbative QCD predictions give a reasonable description of the shape of the measured cross sections over most of the kinematic range, but the absolute normalisation is typically in disagreement by 20-30%.

Single- and double-differential inclusive dijet cross sections in neutral current deep inelastic ep scattering have been measured with the ZEUS detector using an integrated luminosity of 374 pb^-1. The measurement was performed at large values of the photon virtuality, Q^2, between 125 and 20000 GeV^2. The jets were reconstructed with the k_T cluster algorithm in the Breit reference frame and selected by requiring their transverse energies in the Breit frame, E_T,B^jet, to be larger than 8 GeV. In addition, the invariant mass of the dijet system, M_jj, was required to be greater than 20 GeV. The cross sections are described by the predictions of next-to-leading-order QCD.