Riemann surfaces are geometric constructions in complex analysis that may represent multi-valued holomorphic functions using multiple sheets of the complex plane. We show that the energy dispersion of surface states in topological semimetals can be represented by Riemann surfaces generated by holomorphic functions in the two-dimensional momentum space, whose constant height contours correspond to Fermi arcs. This correspondence is demonstrated in the recently discovered Weyl semimetals and leads us to predict new types of topological semimetals, whose surface states are represented by double-and quad-helicoid Riemann surfaces. The intersection of multiple helicoids, or the branch cut of the generating function, appears on high-symmetry lines in the surface Brillouin zone, where surface states are guaranteed to be doubly degenerate by a glide reflection symmetry. We predict the heterostructure superlattice [(SrIrO3)2(CaIrO3)2] to be a topological semimetal with double-helicoid Riemann surface states.Introduction The study of topological semimetals [1] has seen rapid progress since the theoretical proposal of a three-dimensional Weyl semimetal in a magnetic phase of pyrochlore iridates [2]. In general, topological semimetals are materials where the conduction and the valence bands cross in the Brillouin zone and the crossing cannot be removed by perturbations preserving certain crystalline symmetry such as the lattice translation. Bloch states in the vicinity of the band crossing possess a nonzero topological index, e.g., the Chern number in case of Weyl semimetals. The nontrivial topology gives rise to anomalous bulk properties of topological semimetals such as the chiral anomaly [3][4][5]. Several classes of topological semimetals have been theoretically proposed so far, including Weyl, [2,[6][7][8][9][10][11][12][13], Dirac [14][15][16][17] and nodal line semimetals [7,[17][18][19][20][21][22][23][24][25][26][27][28][29][30], some among which have been experimentally observed [31][32][33][34][35][36][37][38][39][40][41][42][43][44][45][46][47].Surface states of topological semimetals have attracted much attention. On the surface of a Weyl semimetal, the Fermi surface consists of open arcs connecting the projection of bulk Weyl points onto the surface Brillouin zone [2], instead of closed loops. The presence of Fermi arcs on the surface is a remarkable property that directly reflects the nontrivial topology of the bulk, and plays a key role in the experimental identification of Weyl semimetals [32,33,36]. In contrast, as shown by recent theoretical works [18,21,23,24,28,48,49], existing Dirac and nodal line semimetals do not have robust Fermi arcs that are stable against symmetry-allowed perturbations. Therefore, the general condition for protected Fermi arcs and their existence in topological semimetals beyond Weyl remain open questions.In this work, we report the discovery of a new topological semimetal phase in a wide variety of nonsymmorphic crystal structures with the glide reflection sym-