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.

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.

Exclusive production of $\rho^0$ and $J/\psi$ mesons in e^+ p collisions has been studied with the ZEUS detector in the kinematic range $0.25 < Q^2 < 50 GeV^2, 20 < W < 167 GeV$ for the $\rho^0$ data and $2 < Q^2 < 40 GeV^2, 50 < W < 150 GeV$ for the $J/\psi$ data. Cross sections for exclusive $\rho^0$ and $J/\psi$ production have been measured as a function of $Q^2, W$ and $t$. The spin-density matrix elements $r^{04}_{00}, r^1_{1-1}$ and $Re r^{5}_{10}$ have been determined for exclusive $\rho^0$ production as well as $r^{04}_{00}$ and $r^{04}_{1-1}$ for exclusive $J/\psi$ production. The results are discussed in the context of theoretical models invoking soft and hard phenomena.

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.

The cross section and the proton structure function F2 for neutral current deep inelastic e+p scattering have been measured with the ZEUS detector at HERA using an integrated luminosity of 30 pb-1. The data were collected in 1996 and 1997 at a centre-of-mass energy of 300 GeV. They cover the kinematic range 2.7 < Q^2 < 30000 GeV2 and 6.10^-5 < x < 0.65. The variation of F2 with x and Q2 is well described by next-to-leading-order perturbative QCD as implemented in the DGLAP evolution equations.

Cross sections for e^-p neutral current deep inelastic scattering have been measured at a centre-of-mass energy of 318 GeV using an integrated luminosity of 15.9 pb^-1 collected with the ZEUS detector at HERA. Results on the double-differential cross-section d^2s/dxdQ^2 in the range 185 < Q^2 < 50000 GeV^2 and 0.0037 < x < 0.75, as well as the single-differential cross-sections ds/dQ^2, ds/dx and ds/dy for Q^2 > 200 GeV^2, are presented. To study the effect of Z-boson exchange, ds/dx has also been measured for Q^2 > 10000 GeV^2. The structure function xF_3 has been extracted by combining the e^-p results presented here with the recent ZEUS measurements of e^+p neutral current deep inelastic scattering. All results agree well with the predictions of the Standard Model.

Cross sections for e^+p neutral current deep inelastic scattering have been measured at a centre-of-mass energy of sqrt{s}=318 GeV with the ZEUS detector at HERA using an integrated luminosity of 63.2 pb^-1. The double-differential cross section, d^2sigma/dxdQ^2, is presented for 200 GeV^2 < Q^2 < 30000 GeV^2 and for 0.005 < x < 0.65. The single-differential cross-sections dsigma/dQ^2, dsigma/dx and dsigma/dy are presented for Q^2 > 200 GeV^2. The effect of Z-boson exchange is seen in dsigma/dx measured for Q^2 > 10000 GeV^2. The data presented here were combined with ZEUS e^+p neutral current data taken at sqrt{s}=300 GeV and the structure function F_2^{em} was extracted. All results agree well with the predictions of the Standard Model.

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.