A combination is presented of the inclusive deep inelastic cross sections measured by the H1 and ZEUS Collaborations in neutral and charged current unpolarised ep scattering at HERA during the period 1994-2000. The data span six orders of magnitude in negative four-momentum-transfer squared, Q^2, and in Bjorken x. The combination method used takes the correlations of systematic uncertainties into account, resulting in an improved accuracy. The combined data are the sole input in a NLO QCD analysis which determines a new set of parton distributions HERAPDF1.0 with small experimental uncertainties. This set includes an estimate of the model and parametrisation uncertainties of the fit result.

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.

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.

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.

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.

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