• ### Observation of Anisotropy of TeV Cosmic Rays with Two Years of HAWC(1805.01847)

May 7, 2018 astro-ph.HE
After two years of operation, the High-Altitude Water Cherenkov (HAWC) Observatory has analyzed the TeV cosmic-ray sky over an energy range between $2.0$ and $72.8$ TeV. The HAWC detector is a ground-based air-shower array located at high altitude in the state of Puebla, Mexico. Using 300 light-tight water tanks, it collects the Cherenkov light from the particles of extensive air showers from primary gamma rays and cosmic rays. This detection method allows for uninterrupted observation of the entire overhead sky (2~sr instantaneous, 8.5~sr integrated) in the energy range from a few TeV to hundreds of TeV. Like other detectors in the northern and southern hemisphere, HAWC observes an energy-dependent anisotropy in the arrival direction distribution of cosmic rays. The observed cosmic-ray anisotropy is dominated by a dipole moment with phase $\alpha\approx40^{\circ}$ and amplitude that slowly rises in relative intensity from $8\times10^{-4}$ at 2 TeV to $14\times10^{-4}$ around 30.3 TeV, above which the dipole decreases in strength. A significant large-scale ($>60^{\circ}$ in angular extent) signal is also observed in the quadrupole and octupole moments, and significant small-scale features are also present, with locations and shapes consistent with previous observations. Compared to previous measurements in this energy range, the HAWC cosmic-ray sky maps improve on the energy resolution and fit precision of the anisotropy. These data can be used in an effort to better constrain local cosmic-ray accelerators and the intervening magnetic fields.
• ### Constraining the $\bar{p}/p$ Ratio in TeV Cosmic Rays with Observations of the Moon Shadow by HAWC(1802.08913)

April 22, 2018 astro-ph.IM, astro-ph.HE
An indirect measurement of the antiproton flux in cosmic rays is possible as the particles undergo deflection by the geomagnetic field. This effect can be measured by studying the deficit in the flux, or shadow, created by the Moon as it absorbs cosmic rays that are headed towards the Earth. The shadow is displaced from the actual position of the Moon due to geomagnetic deflection, which is a function of the energy and charge of the cosmic rays. The displacement provides a natural tool for momentum/charge discrimination that can be used to study the composition of cosmic rays. Using 33 months of data comprising more than 80 billion cosmic rays measured by the High Altitude Water Cherenkov (HAWC) observatory, we have analyzed the Moon shadow to search for TeV antiprotons in cosmic rays. We present our first upper limits on the $\bar{p}/p$ fraction, which in the absence of any direct measurements, provide the tightest available constraints of $\sim1\%$ on the antiproton fraction for energies between 1 and 10 TeV.
• ### Search for Dark Matter Gamma-ray Emission from the Andromeda Galaxy with the High-Altitude Water Cherenkov Observatory(1804.00628)

April 3, 2018 astro-ph.HE
The Andromeda Galaxy (M31) is a nearby ($\sim$780 kpc) galaxy similar to our own Milky Way. Observational evidence suggests that it resides in a large halo of dark matter (DM), making it a good target for DM searches. We present a search for gamma rays from M31 using 1017 days of data from the High Altitude Water Cherenkov (HAWC) Observatory. With its wide field of view and constant monitoring, HAWC is well-suited to search for DM in extended targets like M31. No DM annihilation or decay signal was detected for DM masses from 1 to 100 TeV in the $b\bar{b}$, $t\bar{t}$, $\tau^{+}\tau^{-}$, $\mu^{+}\mu^{-}$, and $W^{+}W^{-}$ channels. Therefore we present limits on those processes. Our limits nicely complement the existing body of DM limits from other targets and instruments. Specifically the DM decay limits from our benchmark model are the most constraining for DM masses from 25 TeV to 100 TeV in the $b\bar{b}, t\bar{t}$ and $\mu^{+}\mu{-}$ channels. In addition to DM-specific limits, we also calculate general gamma-ray flux limits for M31 in 5 energy bins from 1 TeV to 100 TeV.
• ### A Search for Dark Matter in the Galactic Halo with HAWC(1710.10288)

Nov. 3, 2017 hep-ph, astro-ph.CO, astro-ph.HE
The High Altitude Water Cherenkov (HAWC) gamma-ray observatory is a wide field-of-view observatory sensitive to 500 GeV - 100 TeV gamma rays and cosmic rays. With its observations over 2/3 of the sky every day, the HAWC observatory is sensitive to a wide variety of astrophysical sources, including possible gamma rays from dark matter. Dark matter annihilation and decay in the Milky Way Galaxy should produce gamma-ray signals across many degrees on the sky. The HAWC instantaneous field-of-view of 2 sr enables observations of extended regions on the sky, such as those from dark matter in the Galactic halo. Here we show limits on the dark matter annihilation cross-section and decay lifetime from HAWC observations of the Galactic halo with 15 months of data. These are some of the most robust limits on TeV and PeV dark matter, largely insensitive to the dark matter morphology. These limits begin to constrain models in which PeV IceCube neutrinos are explained by dark matter which primarily decays into hadrons.
• List of proceedings from the HAWC Collaboration presented at the 35th International Cosmic Ray Conference, 12 July - 20 July 2017, Bexco, Busan, Korea.
• On February 17 2016, the IceCube real-time neutrino search identified, for the first time, three muon neutrino candidates arriving within 100 s of one another, consistent with coming from the same point in the sky. Such a triplet is expected once every 13.7 years as a random coincidence of background events. However, considering the lifetime of the follow-up program the probability of detecting at least one triplet from atmospheric background is 32%. Follow-up observatories were notified in order to search for an electromagnetic counterpart. Observations were obtained by Swift's X-ray telescope, by ASAS-SN, LCO and MASTER at optical wavelengths, and by VERITAS in the very-high-energy gamma-ray regime. Moreover, the Swift BAT serendipitously observed the location 100 s after the first neutrino was detected, and data from the Fermi LAT and HAWC were analyzed. We present details of the neutrino triplet and the follow-up observations. No likely electromagnetic counterpart was detected, and we discuss the implications of these constraints on candidate neutrino sources such as gamma-ray bursts, core-collapse supernovae and active galactic nucleus flares. This study illustrates the potential of and challenges for future follow-up campaigns.
• ### Search for very-high-energy emission from Gamma-ray Bursts using the first 18 months of data from the HAWC Gamma-ray Observatory(1705.01551)

Aug. 4, 2017 astro-ph.HE
The High Altitude Water Cherenkov (HAWC) Gamma-ray Observatory is an extensive air shower detector operating in central Mexico, which has recently completed its first two years of full operations. If for a burst like GRB 130427A at a redshift of 0.34 and a high-energy component following a power law with index -1.66, the high-energy component is extended to higher energies with no cut-off other than from extragalactic background light attenuation, HAWC would observe gamma rays with a peak energy of $\sim$300 GeV. This paper reports the results of HAWC observations of 64 gamma-ray bursts (GRBs) detected by $\mathit{Swift}$ and $\mathit{Fermi}$, including three GRBs that were also detected by the Large Area Telescope ($\mathit{Fermi}$-LAT). An ON/OFF analysis method is employed, searching on the time scale given by the observed light curve at keV-MeV energies and also on extended time scales. For all GRBs and time scales, no statistically significant excess of counts is found and upper limits on the number of gamma rays and the gamma-ray flux are calculated. GRB 170206A, the third brightest short GRB detected by the Gamma-ray Burst Monitor on board the $\mathit{Fermi}$ satellite ($\mathit{Fermi}$-GBM) and also detected by the LAT, occurred very close to zenith. The LAT measurements can neither exclude the presence of a synchrotron self-Compton (SSC) component nor constrain its spectrum. Instead, the HAWC upper limits constrain the expected cut-off in an additional high-energy component to be less than $100~\rm{GeV}$ for reasonable assumptions about the energetics and redshift of the burst.
• ### Alpha-like clustering in $^{20}$Ne from a quartetting wave function approach(1707.04517)

July 14, 2017 nucl-th
Quartetting ($\alpha$-like clustering) occurs in low density matter ($\le 0.03$ fm$^{-3}$) which exists, e.g., at the surface of nuclei. It is of interest for the $\alpha$ preformation to calculate the $\alpha$ decay of heavy nuclei such as $^{212}$Po, but also in light nuclei (e.g., $^{20}$Ne) which shows strong signatures of quartetting. We analyze the intrinsic structure of the $\alpha$-like cluster and the center of mass motion of the quartet, in particular the role of Pauli blocking. The Thomas-Fermi model for the (daughter) core nucleus is improved introducing quasiparticle nucleon states. Calculations performed for harmonic oscillator basis states show that the effective potential for the quartet center of mass motion remains nearly constant within the core nucleus. The relation to the THSR (Tohsaki-Horiuchi-Schuck-R\"opke) approach is discussed.
• ### Dark Matter Limits From Dwarf Spheroidal Galaxies with The HAWC Gamma-Ray Observatory(1706.01277)

June 5, 2017 hep-ph, astro-ph.HE
The High Altitude Water Cherenkov (HAWC) gamma-ray observatory is a wide field of view observatory sensitive to 500 GeV - 100 TeV gamma rays and cosmic rays. It can also perform diverse indirect searches for dark matter (DM) annihilation and decay. Among the most promising targets for the indirect detection of dark matter are dwarf spheroidal galaxies. These objects are expected to have few astrophysical sources of gamma rays but high dark matter content, making them ideal candidates for an indirect dark matter detection with gamma rays. Here we present individual limits on the annihilation cross section and decay lifetime for 15 dwarf spheroidal galaxies within the HAWC field-of-view, as well as their combined limit. These are the first limits on the annihilation cross section and decay lifetime using data collected with HAWC.
• ### The HAWC real-time flare monitor for rapid detection of transient events(1704.07411)

June 1, 2017 astro-ph.HE
We present the development of a real-time flare monitor for the High Altitude Water Cherenkov (HAWC) observatory. The flare monitor has been fully operational since 2017 January and is designed to detect very high energy (VHE; $E\gtrsim100$ GeV) transient events from blazars on time scales lasting from 2 minutes to 10 hours in order to facilitate multiwavelength and multimessenger studies. These flares provide information for investigations into the mechanisms that power the blazars' relativistic jets and accelerate particles within them, and they may also serve as probes of the populations of particles and fields in intergalactic space. To date, the detection of blazar flares in the VHE range has relied primarily on pointed observations by imaging atmospheric Cherenkov telescopes. The recently completed HAWC observatory offers the opportunity to study VHE flares in survey mode, scanning 2/3 of the entire sky every day with a field of view of $\sim$1.8 steradians. In this work, we report on the sensitivity of the HAWC real-time flare monitor and demonstrate its capabilities via the detection of three high-confidence VHE events in the blazars Markarian 421 and Markarian 501.
• ### Search for Very High Energy Gamma Rays from the Northern $\textit{Fermi}$ Bubble Region with HAWC(1703.01344)

May 24, 2017 astro-ph.HE
We present a search of very high energy gamma-ray emission from the Northern $\textit{Fermi}$ Bubble region using data collected with the High Altitude Water Cherenkov (HAWC) gamma-ray observatory. The size of the data set is 290 days. No significant excess is observed in the Northern $\textit{Fermi}$ Bubble region, hence upper limits above $1\,\text{TeV}$ are calculated. The upper limits are between $3\times 10^{-7}\,\text{GeV}\, \text{cm}^{-2}\, \text{s}^{-1}\,\text{sr}^{-1}$ and $4\times 10^{-8}\,\text{GeV}\,\text{cm}^{-2}\,\text{s}^{-1}\,\text{sr}^{-1}$. The upper limits disfavor a proton injection spectrum that extends beyond $100\,\text{TeV}$ without being suppressed. They also disfavor a hadronic injection spectrum derived from neutrino measurements.
• ### Daily monitoring of TeV gamma-ray emission from Mrk 421, Mrk 501, and the Crab Nebula with HAWC(1703.06968)

May 17, 2017 astro-ph.HE
We present results from daily monitoring of gamma rays in the energy range $\sim0.5$ to $\sim100$ TeV with the first 17 months of data from the High Altitude Water Cherenkov (HAWC) Observatory. Its wide field of view of 2 steradians and duty cycle of $>95$% are unique features compared to other TeV observatories that allow us to observe every source that transits over HAWC for up to $\sim6$ hours each sidereal day. This regular sampling yields unprecedented light curves from unbiased measurements that are independent of seasons or weather conditions. For the Crab Nebula as a reference source we find no variability in the TeV band. Our main focus is the study of the TeV blazars Markarian (Mrk) 421 and Mrk 501. A spectral fit for Mrk 421 yields a power law index $\Gamma=2.21 \pm0.14_{\mathrm{stat}}\pm0.20_{\mathrm{sys}}$ and an exponential cut-off $E_0=5.4 \pm 1.1_{\mathrm{stat}}\pm 1.0_{\mathrm{sys}}$ TeV. For Mrk 501, we find an index $\Gamma=1.60\pm 0.30_{\mathrm{stat}} \pm 0.20_{\mathrm{sys}}$ and exponential cut-off $E_0=5.7\pm 1.6_{\mathrm{stat}} \pm 1.0_{\mathrm{sys}}$ TeV. The light curves for both sources show clear variability and a Bayesian analysis is applied to identify changes between flux states. The highest per-transit fluxes observed from Mrk 421 exceed the Crab Nebula flux by a factor of approximately five. For Mrk 501, several transits show fluxes in excess of three times the Crab Nebula flux. In a comparison to lower energy gamma-ray and X-ray monitoring data with comparable sampling we cannot identify clear counterparts for the most significant flaring features observed by HAWC.
• We present the first catalog of TeV gamma-ray sources realized with the recently completed High Altitude Water Cherenkov Observatory (HAWC). It is the most sensitive wide field-of-view TeV telescope currently in operation, with a 1-year survey sensitivity of ~5-10% of the flux of the Crab Nebula. With an instantaneous field of view >1.5 sr and >90% duty cycle, it continuously surveys and monitors the sky for gamma ray energies between hundreds GeV and tens of TeV. HAWC is located in Mexico at a latitude of 19 degree North and was completed in March 2015. Here, we present the 2HWC catalog, which is the result of the first source search realized with the complete HAWC detector. Realized with 507 days of data and represents the most sensitive TeV survey to date for such a large fraction of the sky. A total of 39 sources were detected, with an expected contamination of 0.5 due to background fluctuation. Out of these sources, 16 are more than one degree away from any previously reported TeV source. The source list, including the position measurement, spectrum measurement, and uncertainties, is reported. Seven of the detected sources may be associated with pulsar wind nebulae, two with supernova remnants, two with blazars, and the remaining 23 have no firm identification yet.
• ### Observation of the Crab Nebula with the HAWC Gamma-Ray Observatory(1701.01778)

Jan. 6, 2017 astro-ph.HE
The Crab Nebula is the brightest TeV gamma-ray source in the sky and has been used for the past 25 years as a reference source in TeV astronomy, for calibration and verification of new TeV instruments. The High Altitude Water Cherenkov Observatory (HAWC), completed in early 2015, has been used to observe the Crab Nebula at high significance across nearly the full spectrum of energies to which HAWC is sensitive. HAWC is unique for its wide field-of-view, nearly 2 sr at any instant, and its high-energy reach, up to 100 TeV. HAWC's sensitivity improves with the gamma-ray energy. Above $\sim$1 TeV the sensitivity is driven by the best background rejection and angular resolution ever achieved for a wide-field ground array. We present a time-integrated analysis of the Crab using 507 live days of HAWC data from 2014 November to 2016 June. The spectrum of the Crab is fit to a function of the form $\phi(E) = \phi_0 (E/E_{0})^{-\alpha -\beta\cdot{\rm{ln}}(E/E_{0})}$. The data is well-fit with values of $\alpha=2.63\pm0.03$, $\beta=0.15\pm0.03$, and log$_{10}(\phi_0~{\rm{cm}^2}~{\rm{s}}~{\rm{TeV}})=-12.60\pm0.02$ when $E_{0}$ is fixed at 7 TeV and the fit applies between 1 and 37 TeV. Study of the systematic errors in this HAWC measurement is discussed and estimated to be $\pm$50\% in the photon flux between 1 and 37 TeV. Confirmation of the Crab flux serves to establish the HAWC instrument's sensitivity for surveys of the sky. The HAWC survey will exceed sensitivity of current-generation observatories and open a new view of 2/3 of the sky above 10 TeV.
• ### Isotope effect in superconducting n-doped SrTiO$_3$(1610.10045)

Nov. 2, 2016 cond-mat.supr-con
We report the influence on the superconducting critical temperature $T_c$ in doped SrTiO$_3$ of the substitution of the natural $^{16}$O atoms by the heavier isotope $^{18}$O. We observe that for a wide range of doping this substitution causes a strong ($\sim 50 \%$) enhancement of $T_c$. Also the magnetic critical field $H_{c2}$ is increased by a factor $\sim 2$. Such a strong impact on $T_c$ and $H_{c2}$, with a sign opposite to conventional superconductors, is unprecedented. The observed effect could be the consequence of strong coupling of the doped electrons to lattice vibrations (phonons), a notion which finds support in numerous optical and photo-emission studies. The unusually large size of the observed isotope effect supports a recent model for superconductivity in these materials based on strong coupling to the ferroelectric soft modes of SrTiO$_{3}$.
• ### Scaling behavior of temperature-dependent thermopower in CeAu2Si2 under pressure(1602.02498)

July 28, 2016 cond-mat.str-el
We report a combined study of in-plane resistivity and thermopower of the pressure-induced heavy fermion superconductor CeAu2Si2 up to 27.8 GPa. It is found that thermopower follows a scaling behavior in T/T* almost up to the magnetic critical pressure pc ~ 22 GPa. By comparing with resistivity results, we show that the magnitude and characteristic temperature dependence of thermopower in this pressure range are governed by the Kondo coupling and crystal-field splitting, respectively. Below pc, the superconducting transition is preceded by a large negative thermopower minimum, suggesting a close relationship between the two phenomena. Furthermore, thermopower of a variety of Ce-based Kondo-lattices with different crystal structures follows the same scaling relation up to T/T* ~ 2.
• ### HAWC Contributions to the 34th International Cosmic Ray Conference (ICRC2015)(1508.03327)

Oct. 8, 2015 hep-ex, astro-ph.HE
List of proceedings from the HAWC Collaboration presented at the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands.
• ### Search for TeV Gamma-Ray Emission from Point-like Sources in the Inner Galactic Plane with a Partial Configuration of the HAWC Observatory(1509.05401)

Sept. 17, 2015 astro-ph.HE
A survey of the inner Galaxy region of Galactic longitude l in [+15, +50] degree and latitude b in [-4,+4] degree is performed using one-third of the High Altitude Water Cherenkov (HAWC) Observatory operated during its construction phase. To address the ambiguities arising from unresolved sources in the data, we use a maximum likelihood technique to identify point source candidates. Ten sources and candidate sources are identified in this analysis. Eight of these are associated with known TeV sources but not all have differential fluxes compatible with previous measurements. Three sources are detected with significances $>5\,\sigma$ after accounting for statistical trials, and are associated with known TeV sources.
• ### Determination of Thermal History by Photoluminescence of Core-shelled Quantum Dots Going Through Heating Events(1508.01622)

Aug. 7, 2015 cond-mat.mtrl-sci
A kind of novel thermal history nanosensors were theoretically designed and experimentally demonstrated to permanently record thermal events. The photoluminescence spectrum of core-shelled quantum dots CdSe/ZnS irreversibly shifted with heating histories (temperature and duration) of thermal events. The induced photoluminescence shift of the quantum dots CdSe/ZnS was employed to permanently record thermal histories. We further modeled a kind of thermal history nanosensor based on the thermal induced phenomena of core-shelled quantum dots to permanently record thermal histories at microscale and demonstrated to reconstruct temperature and duration of heating events simultaneously from photoluminescence spectra of the quantum dots. The physical mechanism of the sensors was discussed.
• ### Effect of disorder on the pressure-induced superconducting state of CeAu2Si2(1501.05172)

March 24, 2015 cond-mat.supr-con
CeAu2Si2 is a newly discovered pressure-induced heavy fermion superconductor which shows very unusual interplay between superconductivity and magnetism under pressure. Here we compare the results of high-pressure measurements on single crystalline CeAu2Si2 samples with different levels of disorder. It is found that while the magnetic properties are essentially sample independent, superconductivity is rapidly suppressed when the residual resistivity of the sample increases. We show that the depression of bulk Tc can be well understood in terms of pair breaking by nonmagnetic disorder, which strongly suggests an unconventional pairing state in pressurized CeAu2Si2. Furthermore, increasing the level of disorder leads to the emergence of another phase transition at T* within the magnetic phase, which might be in competition with superconductivity.
• ### Giant overlap between the magnetic and superconducting phases of CeAu2Si2 under pressure(1406.1438)

Sept. 29, 2014 cond-mat.supr-con
High pressure provides a powerful means for exploring unconventional superconductivity which appears mostly on the border of magnetism. Here we report the discovery of pressure-induced heavy fermion superconductivity up to 2.5 K in the antiferromanget CeAu2Si2 (TN ~ 10 K). Remarkably, the magnetic and superconducting phases are found to overlap across an unprecedentedly wide pressure interval from 11.8 to 22.3 GPa. Moreover, both the bulk Tc and TM are strongly enhanced when increasing the pressure from 16.7 to 20.2 GPa. Tc reaches a maximum at a pressure slightly below pc ~ 22.5 GPa, at which magnetic order disappears. Furthermore, the scaling behavior of the resistivity provides evidence for a continuous delocalization of the Ce 4f-electrons associated with a critical endpoint lying just above pc. We show that the maximum Tc of CeAu2Si2 actually occurs at almost the same unit-cell volume as that of CeCu2Si2 and CeCu2Ge2, and when the Kondo and crystal field splitting energies becomes comparable. Dynamical mean-filed theory calculations suggest that the peculiar behavior in pressurized CeAu2Si2 might be related to its Ce 4f-orbital occupancy. Our results not only provide a unique example of the interplay between superconductivity and magnetism, but also underline the role of orbital physics in understanding Ce-based heavy fermion systems.
• 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.
• ### Measurement of the Luminosity in the ZEUS Experiment at HERA II(1306.1391)

June 13, 2013 hep-ex, physics.ins-det
The luminosity in the ZEUS detector was measured using photons from electron bremsstrahlung. In 2001 the HERA collider was upgraded for operation at higher luminosity. At the same time the luminosity-measuring system of the ZEUS experiment was modified to tackle the expected higher photon rate and synchrotron radiation. The existing lead-scintillator calorimeter was equipped with radiation hard scintillator tiles and shielded against synchrotron radiation. In addition, a magnetic spectrometer was installed to measure the luminosity independently using photons converted in the beam-pipe exit window. The redundancy provided a reliable and robust luminosity determination with a systematic uncertainty of 1.7%. The experimental setup, the techniques used for luminosity determination and the estimate of the systematic uncertainty are reported.
• ### Fermiology of Strongly Spin-Orbit Coupled Superconductor Sn1-xInxTe and its Implication to Topological Superconductivity(1212.5886)

We have performed angle-resolved photoemission spectroscopy of the strongly spin-orbit coupled low-carrier density superconductor Sn1-xInxTe (x = 0.045) to elucidate the electronic states relevant to the possible occurrence of topological superconductivity recently reported for this compound from point-contact spectroscopy. The obtained energy-band structure reveals a small holelike Fermi surface centered at the L point of the bulk Brillouin zone, together with a signature of a topological surface state which indicates that this superconductor is essentially a doped topological crystalline insulator characterized by band inversion and mirror symmetry. A comparison of the electronic states with a band-non-inverted superconductor possessing a similar Fermi surface structure, Pb1-xTlxTe, suggests that the anomalous behavior in the superconducting state of Sn1-xInxTe is likely to be related to the peculiar orbital characteristics of the bulk valence band and/or the presence of a topological surface state.
• The production of Z0 bosons in the reaction ep -> eZ0p*, where p* stands for a proton or a low-mass nucleon resonance, has been studied in ep collisions at HERA using the ZEUS detector. The analysis is based on a data sample collected between 1996 and 2007, amounting to 496 pb-1 of integrated luminosity. The Z0 was measured in the hadronic decay mode. The elasticity of the events was ensured by a cut on eta_max < 3.0, where eta_max is the maximum pseudorapidity of energy deposits in the calorimeter defined with respect to the proton beam direction. A signal was observed at the Z0 mass. The cross section of the reaction ep -> eZ0p* was measured to be sigma(ep -> eZ0p*) = 0.13 +/- 0.06 (stat.) +/- 0.01 (syst.) pb, in agreement with the Standard Model prediction of 0.16 pb. This is the first measurement of Z0 production in ep collisions.