We numerically investigate the surface states of a strong topological insulator in the presence of strong electron-electron interactions. We choose a spherical topological insulator geometry to make the surface amenable to a finite size analysis. The single-particle problem maps to that of Landau orbitals on the sphere with a magnetic monopole at the center that has unit strength and opposite sign for electrons with opposite spin. Assuming density-density contact interactions, we find superconducting and anomalous (quantum) Hall phases for attractive and repulsive interactions, respectively, as well as chiral fermion and chiral Majorana fermion boundary modes between different phases. Our setup is preeminently adapted to the search for topologically ordered surface terminations that could be microscopically stabilized by tailored surface interaction profiles.Introduction.-Three-dimensional topological insulators (3DTIs) [1][2][3][4][5][6], since their prediction in 2007, realize a quantum state of matter which triggered enormous interest in condensed matter physics, and have been subsequently discovered in various material classes [7][8][9][10][11].When viewed as a symmetry-protected topological phase, 3DTIs exhibit a gapped bulk with twodimensional gapless edge states protected by U(1) electron number conservation and time reversal symmetry (TRS), forbidding any adiabatic deformation into a trivial insulator. The effective electromagnetic field theory characterizing the 3DTI contains a topological axion term, with quantized coefficient θ = π for fermionic 3DTIs as compared to θ = 0 for TRS trivial insulators [12]. On the surface of a 3DTI, where θ changes from 0 to π, the gauge invariance of this topological field theory demands the presence of a half-integer Chern-Simons term at the surface. In the absence of symmetry breaking, this Chern-Simons term is precisely what allows a single Dirac cone surface state to be gauge invariant and hence consistently defined on the 3DTI surface, whereas a single massless Dirac field in odd space-time dimensions would usually exhibit a parity anomaly [13][14][15].When the protecting U(1) particle number symmetry is broken, such as by a superconducting proximity effect, the 3DTI surface yields an unconventional gapped s-wave superconductor with Majorana modes in its vortex cores [16]. Upon breaking TRS, such as by a magnetic coating on the surface, the single surface Dirac cone gaps out, and the Chern-Simons boundary term of the axion bulk action manifests itself as a ν = 1/2 quantum Hall effect [12] without fractionalized excitations. The axion term implies the Witten effect [17] by which a odd-half integer charge binds to magnetic monopoles in the bulk of a 3DTI (see also e.g. Ref. 18).All aforementioned properties of 3DTIs do not involve interactions in the bulk or at the surface. Assuming that the gapped 3DTI bulk is negligibly renormalized by interactions, it remains to be investigated how interactions