We report an attractive critical point occurring in the Anderson localization scaling flow of symplectic models on fractals. The scaling theory of Anderson localization predicts that in disordered symplectic two-dimensional systems weak-antilocalization effects lead to a metal-insulator transition. This transition is characterized by a repulsive critical point above which the system becomes metallic. Fractals possess a noninteger scaling of conductance in the classical limit which can be continuously tuned by changing the fractal structure. We demonstrate that in disordered symplectic Hamiltonians defined on fractals with classical conductance scaling g ∼ L −ε , for 0 < ε < β max ≈ 0.15, the metallic phase is replaced by a critical phase with a scale-invariant conductance dependent on the fractal dimensionality. Our results show that disordered fractals allow an explicit construction and verification of the ε expansion. DOI: 10.1103/PhysRevB.94.161115Introduction. The one-parameter scaling hypothesis [1] is central to the study of disordered electronic systems. The hypothesis states that in disordered noninteracting systems the beta function β = d log g/d log L determining the change of conductance g with the system size L is a universal function of g. Single-parameter scaling is known to be violated in quantum Hall systems [2] or topological insulators [3,4], where the topological invariant is the second scaling variable required to capture the scaling flow, and in systems where disorder itself is an irrelevant scaling variable [5][6][7]. Despite that, the scaling flow of Anderson localization holds in an extremely broad range of systems [8,9].The scaling flow has two universal regimes. In the insulating regime g 1 the exponential localization of the wave functions leads to a further decrease of conductance with the system size leading to β ∝ log g + constant. At high conductance, the beta function recovers the classical Ohm law, lim g→∞ β ≡ β ∞ = d − 2 with d the Euclidean dimension. A successful prediction of the theory was the occurrence of a metal-insulator transition in 3d as the flow passes between these two limits. Later studies have refined the theory in the diffusive regime by taking into account quantum corrections to the Ohmic conductance [8,9]. In time-reversal-invariant systems with spin-orbit interactions, also called symplectic, the corrections to g are positive, yielding weak-antilocalization effects [10]. Consequently, these systems exhibit a metal-insulator transition even in 2d, with logarithmic corrections to conductance g ∝ log L, and a metallic phase at large conductance (see Fig. 1).A successful approach in treating Anderson localization is the ε expansion, which treats d as a continuous variable and constructs a series expansion of the scaling flow [11][12][13][14]. The ε expansion is a mathematical construct which is not expected to have a physical meaning when d is not integer; nevertheless, fractals are examples of systems with noninteger dimensionality. This motivates the main questi...