Recent studies on network geometry, a way of describing network structures as geometrical objects, are revolutionizing our way to understand dynamical processes on networked systems. Here, we cope with the problem of epidemic spreading, using the Susceptible-Infected-Susceptible (SIS) model, in simplicial complexes. In particular, we analyze the dynamics of SIS in complex networks characterized by pairwise interactions (links), and three-body interactions (filled triangles, also known as 2-simplices). This higher-order description of the epidemic spreading is analytically formulated using a microscopic Markov chain approximation. The analysis of the fixed point solutions of the model, reveal an interesting phase transition that becomes abrupt with the infectivity parameter of the 2-simplices. Our results are of outmost importance for network theorists to advance in our physical understanding of epidemic spreading in real scenarios where diseases are transmitted among groups as well as among pairs, and to better understand the behavior of dynamical processes in simplicial complexes.The collective behavior of dynamical systems on networks, has been a major subject of research in the physics community during the last decades [1][2][3][4][5]. In particular, our understanding of both natural and man-made systems has significantly improved by studying how network structures and dynamical processes combined shape the overall systems' behavior. Recently, the network science community has turned its attention to network geometry [6-9] to better represent the kinds of interactions that one can find beyond typical pairwise interactions.These higher-order interactions are encoded in geometrical structures that describe the different kinds of simplex structure present in the network: a filled clique of m + 1 nodes is known as an m-simplex, and together a set of 1-simplexes (links), 2-simplexes (filled triangles), etc., comprise the simplicial complex. While simplicial complexes have been proven to be very useful for the analysis and computation in high dimensional data sets, e.g., using persistent homologies [10][11][12][13][14], little is understood about their role in shaping dynamical processes, save for a handful of examples [15][16][17][18].A more accurate description of dynamical processes on complex systems necessarily requires a new paradigm where the network structure representation helps to include higher-order interactions [19]. Simplicial geometry of complex networks is a natural way to extend manybody interactions in complex systems. The standard approach so far consists in understanding the coexistence of two-body (link) interactions and three-body interactions (filled triangles). Note that this approach is different from considering pairwise interactions among three elements of a triangle, it refers to the interaction of the three elements, in the filled triangle, at unison.Here we present a probabilistic formalization of the higher-order interactions of an epidemic process, represented by the well-known Suscep...