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.)

Deep inelastic scattering and its diffractive component, ep -> e'gamma*p ->e'XN, have been studied at HERA with the ZEUS detector using an integrated luminosity of 4.2 pb-1. The measurement covers a wide range in the gamma*p c.m. energy W (37 - 245 GeV), photon virtuality Q2 (2.2 - 80 GeV2) and mass Mx. The diffractive cross section for Mx > 2 GeV rises strongly with W: the rise is steeper with increasing Q2. The latter observation excludes the description of diffractive deep inelastic scattering in terms of the exchange of a single Pomeron. The ratio of diffractive to total cross section is constant as a function of W, in contradiction to the expectation of Regge phenomenology combined with a naive extension of the optical theorem to gamma*p scattering. Above Mx of 8 GeV, the ratio is flat with Q2, indicating a leading-twist behaviour of the diffractive cross section. The data are also presented in terms of the diffractive structure function, F2D(3)(beta,xpom,Q2), of the proton. For fixed beta, the Q2 dependence of xpom F2D(3) changes with xpom in violation of Regge factorisation. For fixed xpom, xpom F2D(3) rises as beta -> 0, the rise accelerating with increasing Q2. These positive scaling violations suggest substantial contributions of perturbative effects in the diffractive DIS cross section.

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%.

Differential cross sections for dijet photoproduction and this process in association with a leading neutron, e+ + p -> e+ + jet + jet + X (+ n), have been measured with the ZEUS detector at HERA using an integrated luminosity of 40 pb-1. The fraction of dijet events with a leading neutron was studied as a function of different jet and event variables. Single- and double-differential cross sections are presented as a function of the longitudinal fraction of the proton momentum carried by the leading neutron, xL, and of its transverse momentum squared, pT**2. The dijet data are compared to inclusive DIS and photoproduction results/ they are all consistent with a simple pion-exchange model. The neutron yield as a function of xL was found to depend only on the fraction of the proton beam energy going into the forward region, independent of the hard process. No firm conclusion can be drawn on the presence of rescattering effects.

Inclusive-jet cross sections have been measured in the reaction ep->e+jet+X for photon virtuality Q2 < 1 GeV2 and gamma-p centre-of-mass energies in the region 142 < W(gamma-p) < 293 GeV with the ZEUS detector at HERA using an integrated luminosity of 300 pb-1. Jets were identified using the kT, anti-kT or SIScone jet algorithms in the laboratory frame. Single-differential cross sections are presented as functions of the jet transverse energy, ETjet, and pseudorapidity, etajet, for jets with ETjet > 17 GeV and -1 < etajet < 2.5. In addition, measurements of double-differential inclusive-jet cross sections are presented as functions of ETjet in different regions of etajet. Next-to-leading-order QCD calculations give a good description of the measurements, except for jets with low ETjet and high etajet. The influence of non-perturbative effects not related to hadronisation was studied. Measurements of the ratios of cross sections using different jet algorithms are also presented; the measured ratios are well described by calculations including up to O(alphas2) terms. Values of alphas(Mz) were extracted from the measurements and the energy-scale dependence of the coupling was determined. The value of alphas(Mz) extracted from the measurements based on the kT jet algorithm is alphas(Mz) = 0.1206 +0.0023 -0.0022 (exp.) +0.0042 -0.0035 (th.); the results from the anti-kT and SIScone algorithms are compatible with this value and have a similar precision.

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.

Three- and four-jet final states have been measured in photoproduction at HERA using the ZEUS detector with an integrated luminosity of 121 pb^-1. The results are presented for jets with transverse energy E_T^jet>6 GeV and pseudorapidity |eta^jet|<2.4, in the kinematic region given by the virtuality of the photon Q^2<1 GeV^2 and the inelasticity 0.2<y<0.85 and in two mass regions defined as 25<M_nj<50 GeV and M_nj>50 GeV, where M_nj is the invariant mass of the n-jet system. The four-jet photoproduction cross section has been measured for the first time and represents the highest-order process studied at HERA. Both the three- and four-jet cross sections have been compared with leading-logarithmic parton-shower Monte Carlo models, with and without multi-parton interactions. The three-jet cross sections have been compared to an order(alpha alpha_s^2) perturbative QCD calculation.

Charm production in deep inelastic scattering has been measured with the ZEUS detector at HERA using an integrated luminosity of 82 pb^{-1}. Charm has been tagged by reconstructing D^{*+}, D^0, D^{+} and D_s^+ (+ c.c.) charm mesons. The charm hadrons were measured in the kinematic range p_T(D^{*+},D^0,D^{+}) > 3 GeV, p_T(D_s^+)>2 GeV and |\eta(D)| < 1.6 for 1.5 < Q^2 < 1000 GeV^2 and 0.02 < y < 0.7. The production cross sections were used to extract charm fragmentation ratios and the fraction of c quarks hadronising into a particular charm meson in the kinematic range considered. The cross sections were compared to the predictions of next-to-leading-order QCD, and extrapolated to the full kinematic region in p_T(D) and \eta(D) in order to determine the open-charm contribution, F_2^{c\bar{c}}(x,Q^2), to the proton structure function F_2.

The production of energetic neutrons in $ep$ collisions has been studied with the ZEUS detector at HERA. The neutron energy and $p_T^2$ distributions were measured with a forward neutron calorimeter and tracker in a $40 \pb^{-1}$ sample of inclusive deep inelastic scattering (DIS) data and a $6 \pb^{-1}$ sample of photoproduction data. The neutron yield in photoproduction is suppressed relative to DIS for the lower neutron energies and the neutrons have a steeper $p_T^2$ distribution, consistent with the expectation from absorption models. The distributions are compared to HERA measurements of leading protons. The neutron energy and transverse-momentum distributions in DIS are compared to Monte Carlo simulations and to the predictions of particle exchange models. Models of pion exchange incorporating absorption and additional secondary meson exchanges give a good description of the data.

Exclusive electroproduction of $\phi$ mesons has been studied in $e^\pm p$ collisions at $\sqrt{s}=318 \gev$ with the ZEUS detector at HERA using an integrated luminosity of 65.1 pb$^{-1}$. The $\gamma^*p$ cross section is presented in the kinematic range $2<Q^2<70 \gev^2$, $35<W<145 \gev$ and ${|t|}<0.6 \gev^2$. The cross sections as functions of $Q^2$, $W$, $t$ and helicity angle $\theta_h$ are compared to cross sections for other vector mesons. The ratios $R$ of the cross sections for longitudinally and transversely polarized virtual photons are presented as functions of $Q^2$ and $W$. The data are also compared to predictions from theoretical models.