The Hikami-Larkin-Nagaoka (HLN) formula [Prog. Theor. Phys. 63, 707 (1980)] describes the quantum corrections to the magnetoconductivity of a quasi-2D electron gas (quasi-2DEG) with parabolic dispersion. It predicts a crossover from weak localization to antilocalization as a function of the strength of scattering off spin-orbit impurities. Here, we derive the conductivity correction for massless Dirac fermions in 3D topological insulators (3DTIs) in the presence of spin-orbit impurities. We show that this correction is always positive and therefore we predict weak antilocalization for every value of the spin-orbit disorder. Furthermore, the correction to the diffusion constant is surprisingly linear in the strength of the impurity spin-orbit. Our results call for a reinterpretation of experimental fits for the magnetoconductivity of 3D TIs which have so far used the standard HLN formula.PACS numbers: 73.20.Fz, 73.43.Qt, 73.25.+i Introduction. The problem of the diffusion of the surface states of 3DTIs is a complex one due to a variety of competing phenomena. The most striking is the fact that these surface states are described by a Dirac Hamiltonian [1,2], which gives rise to weak antilocalization (WAL) in the presence of scalar disorder [3][4][5]. The WAL correction can be affected by the interaction of the surface states with the residual bulk states [6], the thickness of the film [7], or electron-electron interactions [8,9] changing its sign and turning it into weak localization (WL). At the same time, since 3DTIs have strong spin-orbit coupling, one expects spin-orbit coupled impurities to have a strong effect on transport, different however than in graphene where two valleys are present [10,11]. Surprisingly this problem has received virtually no attention [12] and so far the Dirac nature of the states, that manifest itself in the angular dependence of the Green functions, has not been taken into account in studies of spin-orbit impurities [7]. The problem is even more complicated in a transverse magnetic field.The formula commonly used to fit the magnetoconductance experiments on 3DTIs [13][14][15][16][17][18][19][20] was derived by Hikami, Larkin and Nagaoka (HLN) [21]. The HLN formula, however, lacks important features relevant to 3DTIs: it is derived for quasi-2DEGs with a parabolic electron dispersion (impurities are treated as threedimensional objects), so it accounts neither for the Dirac nature of the surface states nor for their strictly twodimensional character. Fig. 1 is the main result of this paper, and shows the different effects of spin-orbit scatter- * The two first authors contributed equally.† Author to whom correspondence should be addressed : hankiewicz@physik.uni-wuerzburg.de ing for the HLN formula, namely WL to WAL crossover with the strength of spin-orbit impurity scattering, and for Dirac fermions, where WAL appears regardless of the strength of spin-orbit impurities. We only observe a convergence of the two formulas for large values of the strength of the impurity spin-orbit sc...