  Percolation on complex networks is used both as a model for dynamics on
networks, such as network robustness or epidemic spreading, and as a benchmark
for our models of networks, where our ability to predict percolation measures
our ability to describe the networks themselves. In many applications,
correctly identifying the phase transition of percolation on real-world
networks is of critical importance. Unfortunately, this phase transition is
obfuscated by the finite size of real systems, making it hard to distinguish
finite size effects from the inaccuracy of a given approach that fails to
capture important structural features. Here, we borrow the perspective of
smeared phase transitions and argue that many observed discrepancies are due to
the complex structure of real networks rather than to finite size effects only.
In fact, several real networks often used as benchmarks feature a smeared phase
transition where inhomogeneities in the topological distribution of the order
parameter do not vanish in the thermodynamic limit. We find that these smeared
transitions are sometimes better described as sequential phase transitions
within correlated subsystems. Our results shed light not only on the nature of
the percolation transition in complex systems, but also provide two important
insights on the numerical and analytical tools we use to study them. First, we
propose a measure of local susceptibility to better detect both clean and
smeared phase transitions by looking at the topological variability of the
order parameter. Second, we highlight a shortcoming in state-of-the-art
analytical approaches such as message passing, which can detect smeared
transitions but not characterize their nature.
