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.

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

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.

Inclusive dijet and trijet production in deep inelastic $ep$ scattering has been measured for $10&lt;Q^2&lt;100$ GeV$^2$ and low Bjorken $x$, $10^{-4}&lt;x_{\rm Bj}&lt;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)$.

The production of the neutral strange hadrons $K^{0}_{S}$, $\Lambda$ and $\bar{\Lambda}$ has been measured in $ep$ collisions at HERA using the ZEUS detector. Cross sections, baryon-to-meson ratios, relative yields of strange and charged light hadrons, $\Lambda$ ($\bar{\Lambda}$) asymmetry and polarization have been measured in three kinematic regions: $Q^2 > 25 \gev^2$: $5 < Q^2 < 25 \gev^2$: and in photoproduction ($Q^2 \simeq 0$). In photoproduction the presence of two hadronic jets, each with at least $5 \gev$ transverse energy, was required. The measurements agree in general with Monte Carlo models and are consistent with measurements made at $e^+ e^-$ colliders, except for an enhancement of baryon relative to meson production in photoproduction.

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.

A new method is employed to measure the neutral current cross section up to Bjorken-x values of one with the ZEUS detector at HERA using an integrated luminosity of 65.1 pb-1 for e+p collisions and 16.7 pb-1 for e-p collisions at sqrt{s}=318 GeV and 38.6 pb-1 for e+p collisions at sqrt{s}=300 GeV. Cross sections have been extracted for Q2 >= 648 GeV2 and are compared to predictions using different parton density functions. For the highest x bins, the data have a tendency to lie above the expectations using recent parton density function parametrizations.

Measurements of neutral current cross sections for deep inelastic scattering in e+p collisions at HERA with a longitudinally polarised positron beam are presented. The single-differential cross-sections d(sigma)/dQ2, d(sigma)/dx and d(sigma)/dy and the reduced cross-section were measured in the kinematic region Q2 > 185 GeV2 and y < 0.9, where Q2 is the four-momentum transfer squared, x the Bjorken scaling variable, and y the inelasticity of the interaction. The measurements were performed separately for positively and negatively polarised positron beams. The measurements are based on an integrated luminosity of 135.5 pb-1 collected with the ZEUS detector in 2006 and 2007 at a centre-of-mass energy of 318 GeV. The structure functions F3 and F3(gamma)Z were determined by combining the e+p results presented in this paper with previously published e-p neutral current results. The asymmetry parameter A+ is used to demonstrate the parity violation predicted in electroweak interactions. The measurements are well described by the predictions of the Standard Model.

The production of D*+-(2010) mesons in deep inelastic scattering has been measured in the ZEUS detector at HERA using an integrated luminosity of 37 pb^-1. The decay channels D*+ -> D0 pi+(+c.c.), with D0 -> K- pi+ or D0 ->K- pi- pi+ pi+, have been used to identify the D mesons. The e+p cross section for inclusive D*+- production with 1<Q^2<600 GeV^2 and 0.02<y<0.7 is 8.31 +- 0.31(stat.) +0.30-0.50(syst.) nb in the kinematic region 1.5< pT(D*+-)<15 GeV and |eta(D*+-)|<1.5. Differential cross sections are consistent with a next-to-leading-order perturbative-QCD calculation when using charm-fragmentation models which take into account the interaction of the charm quark with the proton remnant. The observed cross section is extrapolated to the full kinematic region in pT(D*+-) and eta(D*+-) in order to determine the charm contribution, F^ccbar_2(x,Q^2), to the proton structure function. The ratio F^ccbar_2/F_2 rises from ~10% at Q^2 ~1.8 GeV^2 to ~30% at Q^2 ~130 GeV^2 for x values in the range 10^-4 to 10-3.