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

The inelastic production of J/psi mesons in e p collisions has been studied with the ZEUS detector at HERA using an integrated luminosity of 109 pb-1. The J/psi mesons were identified using the decay channel J/psi -> mu+ mu-. The measurements were performed in the kinematic range 2 < Q^2< 80 Gev^2, 50 < W < 250 Gev, 0.2 < z <0.9 and -1.6 < Ylab < 1.3, where Q^2 is the virtuality of the exchanged photon, W is the photon-proton centre-of-mass energy, z is the fraction of the photon energy carried by the J/psi meson in the proton rest frame and Ylab is the rapidity of the J/psi in the laboratory frame. The measured cross sections are compared to theoretical predictions within the non-relativistic QCD framework including colour-singlet and colour-octet contributions, as well as to predictions based on the kT-factorisation approach. Calculations of the colour-singlet process generally agree with the data, whereas inclusion of colour-octet terms spoils this agreement.

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 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 photoproduction of $D^{*\pm} (2010)$ mesons associated with a leading neutron has been observed with the ZEUS detector in $ep$ collisions at HERA using an integrated luminosity of 80 pb$^{-1}$. The neutron carries a large fraction, {$x_L>0.2$}, of the incoming proton beam energy and is detected at very small production angles, {$\theta_n<0.8$ mrad}, an indication of peripheral scattering. The $D^*$ meson is centrally produced with pseudorapidity {$|\eta|<1.5$}, and has a transverse momentum {$p_{\it T} > 1.9$ GeV}, which is large compared to the average transverse momentum of the neutron of 0.22 GeV. The ratio of neutron-tagged to inclusive $D^*$ production is $8.85\pm 0.93({\rm stat.})^{+0.48}_{-0.61}({\rm syst.})\%$ in the photon-proton center-of-mass energy range {$130 <W<280$ GeV}. The data suggest that the presence of a hard scale enhances the fraction of events with a leading neutron in the final state.

The cross sections for inelastic photoproduction of J/psi and psi' mesons have been measured in ep collisions with the ZEUS detector at HERA, using an integrated luminosity of 38.0 pb-1. The events were required to have 0.1 < z < 0.9 and 50 < W < 180 GeV, where z is the fraction of the incident photon energy carried by the J/psi in the proton rest frame and W is the photon-proton centre-of-mass energy. The psi' to J/psi cross-section ratio was measured in the range 0.55 < z < 0.9. The J/psi data, for various ranges of transverse momentum, are compared to theoretical models incorporating colour-singlet and colour-octet matrix elements. Predictions of a next-to-leading-order colour-singlet model give a good description of the data, although there is a large normalisation uncertainty. The J/psi helicity distribution for z > 0.4 is compared to leading-order QCD predictions.