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.

Photoproduction at HERA is studied in $ep$ collisions, with the ZEUS detector, for $\gamma p$ centre-of-mass energies ranging from 130-270 GeV. A sample of events with two high-$p_T$ jets ($p_T > 6$ GeV, $\eta <1.6$) and a third cluster in the approximate direction of the electron beam is isolated using a clustering algorithm. These events are mostly due to resolved photoproduction. The third cluster is identified as the photon remnant. Its properties, such as the transverse and longitudinal energy flows around the axis of the cluster, are consistent with those commonly attributed to jets, and in particular with those found for the two jets in these events. The mean value of the photon remnant $p_T$ with respect to the beam axis is measured to be $2.1 \pm 0.2$ GeV, which demonstrates substantial mean transverse momenta for the photon remnant.

Inclusive photoproduction of D*+- mesons has been measured for photon-proton centre-of-mass energies in the range 130 < W < 280 GeV and a photon virtuality Q^2 < 1 GeV^2. The data sample used corresponds to an integrated luminosity of 37 pb^-1. Total and differential cross sections as functions of the D* transverse momentum and pseudorapidity are presented in restricted kinematical regions and the data are compared with next-to-leading order (NLO) perturbative QCD calculations using the "massive charm" and "massless charm" schemes. The measured cross sections are generally above the NLO calculations, in particular in the forward (proton) direction. The large data sample also allows the study of dijet production associated with charm. A significant resolved as well as a direct photon component contribute to the cross section. Leading order QCD Monte Carlo calculations indicate that the resolved contribution arises from a significant charm component in the photon. A massive charm NLO parton level calculation yields lower cross sections compared to the measured results in a kinematic region where the resolved photon contribution is significant.

Charged particles ($h^\pm$) and \kz mesons have been studied in photoproduced events containing at least one jet of $E_T > 8$ GeV in a pseudorapidity interval (--0.5, 0.5) in the ZEUS laboratory frame. Distributions are presented in terms of transverse momentum, pseudorapidity and distance of the particle from the axis of a jet. The properties of \hpm within the jet are described well using the standard settings of PYTHIA, but the use of the multiparton interaction option improves the description outside the jets. A reasonable overall description of the \kz behaviour is possible with PYTHIA using a reduced value of the strangeness suppression parameter. The numbers of $h^\pm$ and \kz within a jet as defined above are measured to be $3.25\pm0.02\pm0.28$ and $0.431\pm0.013\pm0.088$ respectively. Fragmentation functions are presented for $h^\pm$ and \kz in photoproduced jets; agreement is found with calculations of Binnewies et al. and, at higher momenta, with $p\bar p$ scattering and with standard PYTHIA. Fragmentation functions in direct photoproduced events are extracted, and at higher momenta give good agreement with data from related processes in $e^+e^-$ annihilation and deep inelastic $ep$ scattering.

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.

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

Inclusive jet cross sections in photoproduction for events containing a $D^*$ meson have been measured with the ZEUS detector at HERA using an integrated luminosity of $78.6 {\rm pb}^{-1}$. The events were required to have a virtuality of the incoming photon, $Q^2$, of less than 1 GeV$^2$, and a photon-proton centre-of-mass energy in the range $130&lt;W_{\gamma p}&lt;280 {\rm GeV}$. The measurements are compared with next-to-leading-order (NLO) QCD calculations. Good agreement is found with the NLO calculations over most of the measured kinematic region. Requiring a second jet in the event allowed a more detailed comparison with QCD calculations. The measured dijet cross sections are also compared to Monte Carlo (MC) models which incorporate leading-order matrix elements followed by parton showers and hadronisation. The NLO QCD predictions are in general agreement with the data although differences have been isolated to regions where contributions from higher orders are expected to be significant. The MC models give a better description than the NLO predictions of the shape of the measured cross sections.

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 and semi-leptonic decay of heavy quarks have been studied in the photoproduction process $e^+p -> e^+ + {dijet} + e^- + X with the ZEUS detector at HERA using an integrated luminosity of 38.5 ${\rm pb^{-1}}$. Events with photon-proton centre-of-mass energies, $W_{\gamma p}$, between 134 and 269 GeV and a photon virtuality, Q^2, less than 1 ${\rm GeV^2}$ were selected requiring at least two jets of transverse energy $E_T^{\rm jet1(2)} >7(6)$ GeV and an electron in the final state. The electrons were identified by employing the ionisation energy loss measurement. The contribution of beauty quarks was determined using the transverse momentum of the electron relative to the axis of the closest jet, $p_T^{\rm rel}$. The data, after background subtraction, were fit with a Monte Carlo simulation including beauty and charm decays. The measured beauty cross section was extrapolated to the parton level with the b quark restricted to the region of transverse momentum $p_T^{b} > p_T^{\rm min} =$ 5 GeV and pseudorapidity $|\eta^{b}| <$ 2. The extrapolated cross section is $1.6 \pm 0.4 (stat.)^{+0.3}_{-0.5} (syst.) ^{+0.2}_{-0.4} (ext.) {nb}$. The result is compared to a perturbative QCD calculation performed to next-to-leading order.