The production of neutrons carrying at least 20% of the proton beam energy ($\xl > 0.2$) in $e^+p$ collisions has been studied with the ZEUS detector at HERA for a wide range of $Q^2$, the photon virtuality, from photoproduction to deep inelastic scattering. The neutron-tagged cross section, $e p\to e' X n$, is measured relative to the inclusive cross section, $e p\to e' X$, thereby reducing the systematic uncertainties. For $\xl >$ 0.3, the rate of neutrons in photoproduction is about half of that measured in hadroproduction, which constitutes a clear breaking of factorisation. There is about a 20% rise in the neutron rate between photoproduction and deep inelastic scattering, which may be attributed to absorptive rescattering in the $\gamma p$ system. For $0.64 < \xl < 0.82$, the rate of neutrons is almost independent of the Bjorken scaling variable $x$ and $Q^2$. However, at lower and higher $\xl$ values, there is a clear but weak dependence on these variables, thus demonstrating the breaking of limiting fragmentation. The neutron-tagged structure function, ${{F}^{\rm\tiny LN(3)}_2}(x,Q^2,\xl)$, rises at low values of $x$ in a way similar to that of the inclusive \ff of the proton. The total $\gamma \pi$ cross section and the structure function of the pion, $F^{\pi}_2(x_\pi,Q^2)$ where $x_\pi = x/(1-\xl)$, have been determined using a one-pion-exchange model, up to uncertainties in the normalisation due to the poorly understood pion flux. At fixed $Q^2$, $F^{\pi}_2$ has approximately the same $x$ dependence as $F_2$ of the proton.

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.

Inclusive photoproduction of $\dspm$ in ep collisions at HERA has been measured with the ZEUS detector for photon-proton centre of mass energies in the range \linebreak \wrang and photon virtuality Q~2 < 4 \g2. The cross section $\sigma_{ep \to \ds X} $ integrated over the kinematic region \ptrangand \etarang is {\xsecs}. Differential cross sections as functions of $p_{\perp}~{\ds}$, $\eta~{\ds}$ and W are given. The data are compared with two next-to-leading order perturbative QCD predictions. For a calculation using a massive charm scheme the predicted cross sections are smaller than the measured ones. A recent calculation using a massless charm scheme is in agreement with the data.

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.

The photoproduction of beauty quarks in events with two jets and a muon has been measured with the ZEUS detector at HERA using an integrated luminosity of 110 pb$^{- 1}$. The fraction of jets containing b quarks was extracted from the transverse momentum distribution of the muon relative to the closest jet. Differential cross sections for beauty production as a function of the transverse momentum and pseudorapidity of the muon, of the associated jet and of $x_{\gamma}^{jets}$, the fraction of the photon's momentum participating in the hard process, are compared with MC models and QCD predictions made at next-to-leading order. The latter give a good description of the data.

The production of $D^{*\pm}$ mesons in deep inelastic $ep$ scattering has been measured for exchanged photon virtualities $ 5<Q^2<1000 \gev^2 $, using an integrated luminosity of 363 pb$^{-1}$ with the ZEUS detector at HERA. Differential cross sections have been measured and compared to next-to-leading-order QCD calculations. The cross-sections are used to extract the charm contribution to the proton structure functions, expressed in terms of the reduced charm cross section, $\sigma_{\rm red}^{c\bar{c}}$. Theoretical calculations based on fits to inclusive HERA data are compared to the results.

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.

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

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