Characteristics of the hadronic final state of diffractive deep inelastic scattering events, ep -> eXp, were studied in the kinematic range 4 < M_X < 35 GeV, 4 < Q^2 < 150 GeV^2, 70 < W < 250 GeV and 0.0003 < x_pom < 0.03 with the ZEUS detector at HERA using an integrated luminosity of 13.8 pb^{-1}. The events were tagged by identifying the diffractively scattered proton using the leading proton spectrometer. The properties of the hadronic final state, X, were studied in its center-of-mass frame using thrust, thrust angle, sphericity, energy flow, transverse energy flow and ``seagull'' distributions. As the invariant mass of the system increases, the final state becomes more collimated, more aligned and more asymmetric in the average transverse momentum with respect to the direction of the virtual photon. Comparisons of the properties of the hadronic final state with predictions from various Monte Carlo model generators suggest that the final state is dominated by qqg states at the parton level.

The forward-jet cross section in deep inelastic ep scattering has been measured using the ZEUS detector at HERA with an integrated luminosity of 6.36 pb^-1. The jet cross section is presented as a function of jet transverse energy squared, E(T,jet)^2, and Q^2 in the kinematic ranges 10^-2<E(T,jet)^2/Q^2<10^2 and 2.5 10^-4<x<8.0 10^-2. Since the perturbative QCD predictions for this cross section are sensitive to the treatment of the log(E_T/Q)^2 terms, this measurement provides an important test. The measured cross section is compared to the predictions of a next-to-leading order pQCD calculation as well as to various leading-order Monte Carlo models. Whereas the predictions of all models agree with the measured cross section in the region of small E(T,Jet)^2/Q^2, only one model, which includes a resolved photon component, describes the data over the whole kinematic range.

First inclusive measurements of isolated prompt photons in photoproduction at the HERA ep collider have been made with the ZEUS detector, using an integrated luminosity of 38.4 pb$^{-1}$. Cross sections are given as a function of the pseudorapidity and the transverse energy ($\eta^\gamma$, \eTg) of the photon, for $\eTg > $ 5 GeV in the $\gamma p$ centre-of-mass energy range 134-285 GeV. Comparisons are made with predictions from Monte Carlo models having leading-logarithm parton showers, and with next-to-leading-order QCD calculations, using currently available parameterisations of the photon structure. For forward $\eta^\gamma$ (proton direction) good agreement is found, but in the rear direction all predictions fall below the data.

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.

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.

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.

The exclusive photoproduction of J/psi mesons, gamma p->J/psi p, has been studied in ep collisions with the ZEUS detector at HERA, in the kinematic range 20<W<290 GeV, where W is the photon-proton centre-of-mass energy. The J/psi mesons were reconstructed in the muon and the electron decay channels using integrated luminosities of 38 pb^-1 and 55 pb^-1, respectively. The helicity structure of J/psi production shows that the hypothesis of s-channel helicity conservation is satisfied at the two standard-deviation level. The total cross section and the differential cross-section dsigma/dt, where t is the squared four-momentum transfer at the proton vertex, are presented as a function of W, for |t|<1.8 GeV^2. The t distribution exhibits an exponential shape with a slope parameter increasing logarithmically with W with a value b=4.15 \pm 0.05 (stat.)^{+0.30}_{-0.18} (syst.) GeV^-2 at W=90 GeV. The effective parameters of the Pomeron trajectory are alphapom(0) = 1.200 \pm 0.009(stat.)^{+0.004}_{-0.010}(syst.) and alphappom= 0.115 \pm 0.018(stat.)^{+0.008}_{-0.015}(syst.) GeV^-2.

Dijet production has been studied in neutral current deep inelastic e+p scattering for 470 < Q**2 < 20000 GeV**2 with the ZEUS detector at HERA using an integrated luminosity of 38.4 pb**{-1}. Dijet differential cross sections are presented in a kinematic region where both theoretical and experimental uncertainties are small. Next-to-leading-order (NLO) QCD calculations describe the measured differential cross sections well. A QCD analysis of the measured dijet fraction as a function of Q**2 allows both a precise determination of alpha_s(M_Z) and a test of the energy-scale dependence of the strong coupling constant. A detailed analysis provides an improved estimate of the uncertainties of the NLO QCD cross sections arising from the parton distribution functions of the proton. The value of alpha_s(M_Z), as determined from the QCD fit, is alpha_s(M_Z) = 0.1166 +- 0.0019 (stat.) {+ 0.0024}_{-0.0033} (exp.)} {+ 0.0057}_{- 0.0044} (th.).

Differential cross sections for dijet photoproduction in association with a leading neutron using the reaction e^+ + p --> e^+ + n + jet + jet + X_r have been measured with the ZEUS detector at HERA using an integrated luminosity of 6.4 pb^{-1}. The fraction of dijet events with a leading neutron in the final state was studied as a function of the jet kinematic variables. The cross sections were measured for jet transverse energies E^{jet}_T > 6 GeV, neutron energy E_n > 400 GeV, and neutron production angle theta_n < 0.8 mrad. The data are broadly consistent with factorization of the lepton and hadron vertices and with a simple one-pion-exchange 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.

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 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 cross section for dijet photoproduction at high transverse energies is presented as a function of the transverse energies and the pseudorapidities of the jets. The measurement is performed using a sample of ep-interactions corresponding to an integrated luminosity of 6.3 pb^(-1), recorded by the ZEUS detector.Jets are defined by applying a k_T-clustering algorithm to the hadrons observed in the final state. The measured cross sections are compared to next-to-leading order QCD calculations. In a kinematic regime where theoretical uncertainties are expected to be small, the measured cross sections are higher than these calculations.

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.

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.

Differential dijet cross sections have been measured with the ZEUS detector for photoproduction events in which the hadronic final state containing the jets is separated with respect to the outgoing proton direction by a large rapidity gap. The cross section has been measured as a function of the fraction of the photon (x_gamma^OBS) and pomeron (beta^OBS) momentum participating in the production of the dijet system. The observed x_gamma^OBS dependence shows evidence for the presence of a resolved- as well as a direct-photon component. The measured cross section d(sigma)/d(beta^OBS) increases as beta^OBS increases indicating that there is a sizeable contribution to dijet production from those events in which a large fraction of the pomeron momentum participates in the hard scattering. These cross sections and the ZEUS measurements of the diffractive structure function can be described by calculations based on parton densities in the pomeron which evolve according to the QCD evolution equations and include a substantial hard momentum component of gluons in the pomeron.

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.

The shapes of jets with transverse energies, E_T(jet), up to 45 GeV produced in neutral- and charged-current deep inelastic e+p scattering (DIS) at Q**2 > 100 GeV**2 have been measured with the ZEUS detector at HERA. Jets are identified using a cone algorithm in the eta-phi plane with a cone radius of one unit. The jets become narrower as E_T(jet) increases. The jet shapes in neutral- and charged-current DIS are found to be very similar. The jets in neutral-current DIS are narrower than those in resolved processes in photoproduction and closer to those in direct-photon processes for the same ranges in E_T(jet) and jet pseudorapidity. The jet shapes in DIS are observed to be similar to those in e+e- interactions and narrower than those in pbarp collisions for comparable E_T(jet). Since the jets in e+e- interactions and e+p DIS are predominantly quark initiated in both cases, the similarity in the jet shapes indicates that the pattern of QCD radiation within a quark jet is to a large extent independent of the hard scattering process in these reactions.

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.

The DIS diffractive cross section, $d\sigma^{diff}_{\gamma^* p \to XN}/dM_X$, has been measured in the mass range $M_X < 15$ GeV for $\gamma^*p$ c.m. energies $60 < W < 200$ GeV and photon virtualities $Q^2 = 7$ to 140 GeV$^2$. For fixed $Q^2$ and $M_X$, the diffractive cross section rises rapidly with $W$, $d\sigma^{diff}_{\gamma^*p \to XN}(M_X,W,Q^2)/dM_X \propto W^{a^{diff}}$ with $a^{diff} = 0.507 \pm 0.034 (stat)^{+0.155}_{-0.046}(syst)$ corresponding to a $t$-averaged pomeron trajectory of $\bar{\alphapom} = 1.127 \pm 0.009 (stat)^{+0.039}_{-0.012} (syst)$ which is larger than $\bar{\alphapom}$ observed in hadron-hadron scattering. The $W$ dependence of the diffractive cross section is found to be the same as that of the total cross section for scattering of virtual photons on protons. The data are consistent with the assumption that the diffractive structure function $F^{D(3)}_2$ factorizes according to $\xpom F^{D(3)}_2 (\xpom,\beta,Q^2) = (x_0/ \xpom)^n F^{D(2)}_2(\beta,Q^2)$. They are also consistent with QCD based models which incorporate factorization breaking. The rise of $\xpom F^{D(3)}_2$ with decreasing $\xpom$ and the weak dependence of $F^{D(2)}_2$ on $Q^2$ suggest a substantial contribution from partonic interactions.

Dijet cross sections are presented using photoproduction data obtained with the ZEUS detector during 1994. These measurements represent an extension of previous results, as the higher statistics allow cross sections to be measured at higher jet transverse energy (ETJ). Jets are identified in the hadronic final state using three different algorithms, and the cross sections compared to complete next-to-leading order QCD calculations. Agreement with these calculations is seen for the pseudorapidity dependence of the direct photon events with ETJ > 6 GeV and of the resolved photon events with ETJ > 11 GeV. Calculated cross sections for resolved photon processes with 6 GeV < ETJ < 11 GeV lie below the data.

A measurement of the proton structure function F_2(x,Q^2) is presented in the kinematic range 0.045 GeV^2 < Q^2 < 0.65 GeV^2 and 6*10^{-7} < x < 1*10^{-3}. The results were obtained using a data sample corresponding to an integrated luminosity of 3.9pb^-1 in e^+p reactions recorded with the ZEUS detector at HERA. Information from a silicon-strip tracking detector, installed in front of the small electromagnetic calorimeter used to measure the energy of the final-state positron at small scattering angles, together with an enhanced simulation of the hadronic final state, has permitted the extension of the kinematic range beyond that of previous measurements. The uncertainties in F_2 are typically less than 4%. At the low Q^2 values of the present measurement, the rise of F_2 at low x is slower than observed in HERA data at higher Q^2 and can be described by Regge theory with a constant logarithmic slope. The dependence of F_2 on Q^2 is stronger than at higher Q^2 values, approaching, at the lowest Q^2 values of this measurement, a region where F_2 becomes nearly proportional to Q^2.

Dijet production by almost real photons has been studied at HERA with the ZEUS detector. Jets have been identified using the cone algorithm. A cut on xg, the fraction of the photon energy participating in the production of the two jets of highest transverse energy, is used to define cross sections sensitive to the parton distributions in the proton and in the photon. The dependence of the dijet cross sections on pseudorapidity has been measured for xg $\ge 0.75$ and xg $< 0.75$. The former is sensitive to the gluon momentum density in the proton. The latter is sensitive to the gluon in the photon. The cross sections are corrected for detector acceptance and compared to leading order QCD calculations.

We present measurements of the structure function \Ft\ in $e~+p$ scattering at HERA in the range $3.5\;\Gevsq < \qsd < 5000\;\Gevsq$. A new reconstruction method has allowed a significant improvement in the resolution of the kinematic variables and an extension of the kinematic region covered by the experiment. At $ \qsd < 35 \;\Gevsq$ the range in $x$ now spans $6.3\cdot 10~{-5} < x < 0.08$ providing overlap with measurements from fixed target experiments. At values of $Q~2$ above 1000 GeV$~2$ the $x$ range extends to 0.5. Systematic errors below 5\perc\ have been achieved for most of the kinematic region. The structure function rises as \x\ decreases; the rise becomes more pronounced as \qsd\ increases. The behaviour of the structure function data is well described by next-to-leading order perturbative QCD as implemented in the DGLAP evolution equations.

Measurements of the midrapidity transverse energy distribution, $d\Et/d\eta$, are presented for $p$$+$$p$, $d$$+$Au, and Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV and additionally for Au$+$Au collisions at $\sqrt{s_{_{NN}}}=62.4$ and 130 GeV. The $d\Et/d\eta$ distributions are first compared with the number of nucleon participants $N_{\rm part}$, number of binary collisions $N_{\rm coll}$, and number of constituent-quark participants $N_{qp}$ calculated from a Glauber model based on the nuclear geometry. For Au$+$Au, $\mean{d\Et/d\eta}/N_{\rm part}$ increases with $N_{\rm part}$, while $\mean{d\Et/d\eta}/N_{qp}$ is approximately constant for all three energies. This indicates that the two component ansatz, $dE_{T}/d\eta \propto (1-x) N_{\rm part}/2 + x N_{\rm coll}$, which has been used to represent $E_T$ distributions, is simply a proxy for $N_{qp}$, and that the $N_{\rm coll}$ term does not represent a hard-scattering component in $E_T$ distributions. The $dE_{T}/d\eta$ distributions of Au$+$Au and $d$$+$Au are then calculated from the measured $p$$+$$p$ $E_T$ distribution using two models that both reproduce the Au$+$Au data. However, while the number-of-constituent-quark-participant model agrees well with the $d$$+$Au data, the additive-quark model does not.