The reduced cross sections for ep deep inelastic scattering have been measured with the ZEUS detector at HERA at three different centre-of-mass energies, 318, 251 and 225 GeV. From the cross sections, measured double differentially in Bjorken x and the virtuality, Q^2, the proton structure functions FL and F2 have been extracted in the region 5*10^-4 < x <0.007 and 20 < Q^2 < 130 GeV^2.

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.

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.

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.

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.

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.

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.

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.

The reduced cross sections for $e^{+}p$ deep inelastic scattering have been measured with the ZEUS detector at HERA at three different centre-of-mass energies, $318$, $251$ and $225$ GeV. The cross sections, measured double differentially in Bjorken $x$ and the virtuality, $Q^2$, were obtained in the region $0.13\ \leq\ y\ \leq\ 0.75$, where $y$ denotes the inelasticity and $5\ \leq\ Q^2\ \leq\ 110$ GeV$^2$. The proton structure functions $F_2$ and $F_L$ were extracted from the measured cross sections.

Measurements of the cross sections for charged current deep inelastic scattering in e-p collisions with longitudinally polarised electron beams are presented. The measurements are based on a data sample with an integrated luminosity of 175 pb-1 collected with the ZEUS detector at HERA at a centre-of-mass energy of 318 GeV. The total cross section is given for positively and negatively polarised electron beams. The differential cross-sections dsigma/dQ2, dsigma/dx and dsigma/dy are presented for Q2>200 GeV2. The double-differential cross-section d2sigma/dxdQ2 is presented in the kinematic range 280&lt;Q2&lt;30000 GeV2 and 0.015&lt;x&lt;0.65. The measured cross sections are compared with the predictions of the Standard Model.

Measurements of neutral current cross sections for deep inelastic scattering in e+p collisions at HERA with a longitudinally polarised positron beam are presented. The single-differential cross-sections d(sigma)/dQ2, d(sigma)/dx and d(sigma)/dy and the reduced cross-section were measured in the kinematic region Q2 > 185 GeV2 and y < 0.9, where Q2 is the four-momentum transfer squared, x the Bjorken scaling variable, and y the inelasticity of the interaction. The measurements were performed separately for positively and negatively polarised positron beams. The measurements are based on an integrated luminosity of 135.5 pb-1 collected with the ZEUS detector in 2006 and 2007 at a centre-of-mass energy of 318 GeV. The structure functions F3 and F3(gamma)Z were determined by combining the e+p results presented in this paper with previously published e-p neutral current results. The asymmetry parameter A+ is used to demonstrate the parity violation predicted in electroweak interactions. The measurements are well described by the predictions of the Standard Model.

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 beauty and charm quarks in ep interactions has been studied with the ZEUS detector at HERA for exchanged four-momentum squared 5 < Q^2 < 1000 GeV^2 using an integrated luminosity of 354 pb^{-1}. The beauty and charm content in events with at least one jet have been extracted using the invariant mass of charged tracks associated with secondary vertices and the decay-length significance of these vertices. Differential cross sections as a function of Q^2, Bjorken x, jet transverse energy and pseudorapidity were measured and compared with next-to-leading-order QCD calculations. The beauty and charm contributions to the proton structure functions were extracted from the double-differential cross section as a function of x and Q^2. The running beauty-quark mass, m_b at the scale m_b, was determined from a QCD fit at next-to-leading order to HERA data for the first time and found to be 4.07 \pm 0.14 (fit} ^{+0.01}_{-0.07} (mod.) ^{+0.05}_{-0.00} (param.) ^{+0.08}_{-0.05} (theo) GeV.

Mean values and differential distributions of event-shape variables have been studied in neutral current deep inelastic scattering using an integrated {luminosity} of 82.2 pb$^{-1}$ collected with the ZEUS detector at HERA. The kinematic range was $80 &lt; Q^2 &lt; 20 480\gev^2$ and $0.0024 &lt; x &lt; 0.6$, where $Q^2$ is the virtuality of the exchanged boson and $x$ is the Bjorken variable. The data are compared with a model based on a combination of next-to-leading-order QCD calculations with next-to-leading-logarithm corrections and the Dokshitzer-Webber non-perturbative power corrections. The power-correction method provides a reasonable description of the data for all event-shape variables studied. Nevertheless, the lack of consistency of the determination of $\alpha_s$ and of the non-perturbative parameter of the model, $\albar$, suggests the importance of higher-order processes that are not yet included in the model.

Deep inelastic scattering and its diffractive component, $ep \to e^{\prime}\gamma^* p \to e^{\prime}XN$, have been studied at HERA with the ZEUS detector using an integrated luminosity of 52.4 pb$^{-1}$. The $M_X$ method has been used to extract the diffractive contribution. A wide range in the centre-of-mass energy $W$ (37 -- 245 GeV), photon virtuality $Q^2$ (20 -- 450 GeV$^2$) and mass $M_X$ (0.28 -- 35 GeV) is covered. The diffractive cross section for $2 < M_X < 15$ GeV rises strongly with $W$, the rise becoming steeper as $Q^2$ increases. The data are also presented in terms of the diffractive structure function, $F^{\rm D(3)}_2$, of the proton. For fixed $Q^2$ and fixed $M_X$, $\xpom F^{\rm D(3)}_2$ shows a strong rise as $\xpom \to 0$, where $\xpom$ is the fraction of the proton momentum carried by the Pomeron. For Bjorken-$x < 1 \cdot 10^{-3}$, $\xpom F^{\rm D(3)}_2$ shows positive $\log Q^2$ scaling violations, while for $x \ge 5 \cdot 10^{-3}$ negative scaling violations are observed. The diffractive structure function is compatible with being leading twist. The data show that Regge factorisation is broken.

A new method is employed to measure the neutral current cross section up to Bjorken-x values of one with the ZEUS detector at HERA using an integrated luminosity of 65.1 pb-1 for e+p collisions and 16.7 pb-1 for e-p collisions at sqrt{s}=318 GeV and 38.6 pb-1 for e+p collisions at sqrt{s}=300 GeV. Cross sections have been extracted for Q2 >= 648 GeV2 and are compared to predictions using different parton density functions. For the highest x bins, the data have a tendency to lie above the expectations using recent parton density function parametrizations.

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.

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

Exclusive rho^0 electroproduction at HERA has been studied with the ZEUS detector using 120 pb^{-1} of integrated luminosity collected during 1996-2000. The analysis was carried out in the kinematic range of photon virtuality 2 < Q^2 < 160 GeV$^2, and gamma^* p centre-of-mass energy 32 < W < 180 GeV. The results include the Q^2 and W dependence of the gamma^* p --> rho^0 p cross section and the distribution of the squared-four-momentum transfer to the proton. The helicity analysis of the decay-matrix elements of the rho^0 was used to study the ratio of the gamma^* p cross section for longitudinal and transverse photon as a function of Q^2 and W. Finally, an effective Pomeron trajectory was extracted. The results are compared to various theoretical predictions.

Charged particle production has been studied in neutral current deep inelastic ep scattering with the ZEUS detector at HERA using an integrated luminosity of 0.44 fb^-1. Distributions of scaled momenta in the Breit frame are presented for particles in the current fragmentation region. The evolution of these spectra with the photon virtuality, Q^2, is described in the kinematic region 10&lt;Q^2&lt;41000 GeV^2. Next-to-leading-order and modified leading-log-approximation QCD calculations as well as predictions from Monte Carlo models are compared to the data. The results are also compared to e+e- annihilation data. The dependences of the pseudorapidity distribution of the particles on Q^2 and on the energy in the \gamma p system, W, are presented and interpreted in the context of the hypothesis of limiting fragmentation.

A new measurement of inclusive-jet cross sections in the Breit frame in neutral current deep inelastic scattering using the ZEUS detector at the HERA collider is presented. The data were taken in the years 2004 to 2007 at a centre-of-mass energy of $318\,\text{GeV}$ and correspond to an integrated luminosity of $347\,\text{pb}^{-1}$. Massless jets, reconstructed using the $k_t$-algorithm in the Breit reference frame, have been measured as a function of the squared momentum transfer, $Q^2$, and the transverse momentum of the jets in the Breit frame, $p_{\perp,\text{Breit}}$. The measured jet cross sections are compared to previous measurements and to perturbative QCD predictions. The measurement has been used in a next-to-next-to-leading-order QCD analysis to perform a simultaneous determination of parton distribution functions of the proton and the strong coupling, resulting in a value of $\alpha_s(M_Z^2) = 0.1142 \pm 0.0017~\text{(experimental/fit)}$${}^{+0.0006}_{-0.0007}~\text{(model/parameterisation)}$${}^{+0.0006}_{-0.0004}~\text{(scale)}$, whose accuracy is improved compared to similar measurements. In addition, the running of the strong coupling is demonstrated using data obtained at different scales.

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.