---
title: Fragile Power Domination in Stochastic PMUs
url: https://www.emergentmind.com/topics/fragile-power-domination
type: topic
---

# Fragile Power Domination in Stochastic PMUs

Fragile power domination is a stochastic extension of graph-theoretic power domination in which each phasor measurement unit (PMU) fails independently with probability \(q\in[0,1]\) before the usual domination and zero-forcing propagation begins. Its central object is the expected number of observed vertices, denoted \(\mathcal E(G;S,q)\), which is a polynomial in \(q\) of degree at most \(|S|\). The framework was introduced as “fragile power domination” in “Power domination with random sensor failure” [2312.12259] and was subsequently developed in “On Fragile Power Domination” [2507.14620], where the expected-coverage polynomial is characterized algebraically and shown to be controllable through graph-gadget constructions.

## 1. Formal model

Let \(G=(V,E)\) be a simple graph and let \(S\subseteq V\) be a set, or more generally a multiset, of PMU locations. In classical power domination, the process begins with the domination step
\[
B:=N[S]=\bigl\{\,v\}\cup\{u:uv\in E,\;v\in S\},
\]
or equivalently \(B_0=\bigcup_{v\in S}N[v]\), and then repeats the zero-forcing step: whenever an observed vertex has exactly one neighbor outside the currently observed set, that unique neighbor is added. At termination, the observed set is \(\Obs(G;S)\), and the power-domination number is
\[
\gamma_P(G)=\min\{\lvert S\rvert:S\subseteq V,\;\Obs(G;S)=V\}.
\]
This is the standard deterministic baseline against which fragile power domination is defined [2312.12259].

Fragile power domination inserts an independent sensor-failure step before observation begins. Given \(S\), each PMU survives with probability \(1-q\), independently, producing a random surviving set \(S^*\subseteq S\). The classical power-domination process is then run from \(S^*\). In the notation of [2312.12259], the terminal random observed set is \(\Obs(G;S,q)\) and the associated random variable is
\[
X_{G,S}(q)=\bigl|\Obs(G;S,q)\bigr|.
\]
In the notation emphasized in [2507.14620], the focus is the expected observed cardinality rather than the full distribution.

The model is motivated by power-grid monitoring, where PMUs may break or otherwise fail. The classical objective of full observation is therefore replaced by a stochastic coverage objective, while preserving the same domination and zero-forcing dynamics after failures occur.

## 2. Expected-coverage polynomial

The expected-coverage polynomial is
\[
\mathcal E(G;S,q)=\mathbb E\bigl[\,|\Obs(G;S)|\bigr]
\]
in the notation of [2507.14620], and equivalently
\[
\mathcal E(G;S,q)=\E\bigl[X_{G,S}(q)\bigr]
\]
in the notation of [2312.12259]. Two equivalent formulas are fundamental:
\[
\mathcal E(G;S,q)
=\sum_{W\subseteq S}|\Obs(G;W)|\,q^{\,|S\setminus W|}(1-q)^{|W|}
\]
and
\[
\mathcal E(G;S,q)=\sum_{v\in V}\Pr\bigl[v\text{ is observed}\bigr].
\]
The first is a subset-sum formula obtained by conditioning on the surviving PMU set; the second is a vertex-sum formula expressing expectation as a sum of observation probabilities [2507.14620].

It follows immediately that \(\mathcal E(G;S,q)\) is a polynomial in \(q\) of degree at most \(|S|\). If one writes
\[
\mathcal E(G;S,q)=\sum_{i=0}^{|S|}C_i\,q^i,
\]
then the coefficient of \(q^{\,|S|-r}\), for \(r=0,1,\dots,|S|\), is
\[
C_{|S|-r} =\sum_{\substack{W\subseteq S\\|W|\ge r}} (-1)^{\,|W|-r}\, \binom{|W|}{r}\; \bigl|\Obs(G;W)\bigr|.
\]
This coefficient formula makes explicit that each polynomial coefficient is an inclusion-exclusion aggregate of classical observed-set sizes over all surviving subsets [2312.12259].

Several boundary cases are immediate. The constant term is
\[
\mathcal E(G;S,0)=|\Obs(G;S)|,
\]
so at \(q=0\) the model collapses to classical power domination. At the opposite extreme, large values of \(q\) emphasize placements whose surviving PMUs still provide local coverage, even when the original set \(S\) was chosen to guarantee full domination only in the no-failure regime [2312.12259].

## 3. Algebraic characterization and robustness criteria

A central structural result in [2507.14620] characterizes when two placements have identical expected-coverage polynomials. If \((G,S)\) and \((G',S')\) satisfy \(|S|=|S'|=s\), define
\[
a_k(G,S):=\sum_{W\subseteq S,\;|W|=k}|\Obs(G;W)|,\qquad 0\le k\le s.
\]
Then
\[
\mathcal E(G;S,q)=\mathcal E(G';S',q)
\]
as polynomials in \(q\) if and only if
\[
a_k(G,S)=a_k(G',S')\quad\text{for all }k=0,1,\dots,s.
\]
The proof expands \(\mathcal E\) in the Bernstein-type basis
\[
\mathcal E(G;S,q)=\sum_{k=0}^s a_k(G,S)\,q^{\,s-k}(1-q)^k,
\]
which is a basis for the space of polynomials of degree \(\le s\) [2507.14620]. The consequence is that equality of expected coverage is determined by the aggregated observed-set sizes over all \(k\)-element surviving subsets, not by graph isomorphism or by the fine-grained chronology of zero-forcing.

The same polynomial also encodes worst-case robustness notions. Let \(B\subseteq V(G)\) with \(|B|=b\). Then \(B\) is \(k\)-PMU-defect-robust if and only if
\[
\mathcal E(G;B,q)=n-q^{\,k+1}h(q)
\]
for some polynomial \(h(q)\). The same factorization characterizes \(k\)-fault-tolerant sets:
\[
\mathcal E(G;S,q)=n-q^{k+1}h(q).
\]
Equivalently, the vanishing of the first \(k\) coefficients of \(\mathcal E(G;S,q)\) is exactly the certificate that \(S\) is resilient to any \(k\) PMU failures [2312.12259].

This relationship prevents a common conflation of models. Fragile power domination is not merely a reformulation of worst-case \(k\)-fault tolerance: it replaces the worst-case assumption by an independent-failure model and studies expected coverage as a \(q\)-dependent polynomial. The factorization above shows precisely how the stochastic model bridges to exactly-\(k\)-fault formulations.

## 4. Full observation, coefficient constraints, and polynomial synthesis

Beyond expectation, [2312.12259] studies the probability of complete observation. Writing
\[
p(G;i)=\text{number of power-dominating sets of size }i
\]
and
\[
\mathcal P(G;x)=\sum_{i=1}^np(G;i)\,x^i,
\]
one has
\[
\Pr\bigl[\Obs(G;V,q)=V\bigr]
=\sum_{i=1}^n p(G;i)\; q^{\,n-i}(1-q)^i
=\mathcal P\bigl(G;1-q\bigr)\,q^n.
\]
If \(S\) is \(k\)-PMU-defect-robust, then
\[
\Pr\bigl[\Obs(G;S,q)=V\bigr]
\ge 1-\sum_{i=0}^{\,|S|-k-1}\binom{|S|}{i}\,q^{\,|S|-i}(1-q)^i,
\]
with equality when \(\lvert S\rvert-k=\gamma_P(G)\) [2312.12259]. These formulas isolate the distinction between expected partial coverage and the event of total observability.

At the same time, [2507.14620] shows that many coverage polynomials can be engineered by local graph modifications. For each integer \(k\ge2\) and each subset \(A\subseteq S\), one can attach a small “star-plus-path” or “bipartite-plus-path” gadget at \(A\). By varying the appended path length \(\ell\), one can control the coefficient of \(q^k\) arbitrarily in the total polynomial without affecting higher powers. This leads to the induced-any-polynomial theorem: if \(|S|=d\), then for any integers \(c_d,c_{d-1},\dots,c_2\), there exists a graph \(G'\supseteq G\) as an induced subgraph such that
\[
\mathcal E(G';S,q)=\sum_{k=0}^d c_k\,q^k
\]
for suitable integers \(c_1,c_0\). In particular, the coefficients of all powers \(q^k\) for \(2\le k\le d\) can be prescribed [2507.14620].

This flexibility is not unconstrained. If
\[
\mathcal E(G;S,q)=\sum_{k=0}^d\alpha_k q^k,
\]
then \(\alpha_0=|\Obs(G;S)|\ge0\), \(\alpha_1\le0\), and all other \(\alpha_k\) with \(k\notin\{\,\text{sum of \(\ge1\) sensor-multiplicities}\}\) vanish [2507.14620]. This suggests a sharp asymmetry: intermediate coefficients are highly tunable, while the constant and linear terms retain strong structural restrictions.

## 5. Complexity and canonical graph families

The algorithmic landscape is explicitly intractable in general. Computing \(\gamma_P(G)\) is NP-complete, and by reduction, checking \(k\)-fault-tolerance or \(k\)-PMU-defect-robustness is NP-hard. A naïve computation of \(\mathcal E(G;S,q)\) via the subset-sum formula requires
\[
O(2^{|S|}\cdot T_{\rm PD})
\]
time, where \(T_{\rm PD}\) is the cost of one power-domination run. For trees, and more generally bounded-treewidth graphs, \(\mathcal E\) can be computed in polynomial time by dynamic programming over a tree decomposition. Even selecting the best set \(S\) of a given size to maximize \(\mathcal E\) is combinatorial and intractable in general [2312.12259].

Several graph families admit explicit formulas. For paths \(P_n\), cycles \(C_n\), and complete graphs \(K_n\), one PMU suffices classically, and a single surviving PMU observes all \(n\) vertices. Consequently,
\[
\mathcal E\bigl(P_n;\{v\},q\bigr)=n(1-q),\qquad
\mathcal E\bigl(C_n;\{v\},q\bigr)=n(1-q),\qquad
\mathcal E\bigl(K_n;\{v\},q\bigr)=n(1-q),
\]
with
\[
\Pr[\text{full obs}]=1-q
\]
in each case [2312.12259]. For paths, the same formula holds whether the PMU is placed at an endpoint or at an interior vertex.

For the grid \(G_{m,n}=P_m\square P_n\), the behavior is more delicate. The source notes that \(\gamma_P(G_{m,n})\) is nontrivial but known to be \(2\) for large grids. If two PMUs are placed in opposite corners, full observation occurs whenever at least one PMU survives, yielding
\[
\mathcal E(G_{m,n};\{v_1,v_2\},q)
=mn(1-q^2)+(mn-1)q^2,
\qquad
\Pr[\text{full obs}]=1-q^2.
\]
More general formulas can be obtained by conditioning on which of the two PMUs survives, but in general computing \(\mathcal E\) on grids is exponential [2312.12259].

## 6. Worked example, interpretive significance, and extensions

A worked example in [2507.14620] illustrates how complicated zero-forcing dynamics can collapse to a very simple expected-coverage polynomial. For the graph in Figure 5 of that paper, with three PMUs at \(u,v,w\),
\[
|\Obs(G;\{u\})|=11,\qquad |\Obs(G;\{v\})|=16,\qquad |\Obs(G;\{w\})|=16,
\]
\[
|\Obs(G;\{u,v\})|=26,\qquad |\Obs(G;\{u,w\})|=26,\qquad |\Obs(G;\{v,w\})|=34,
\]
and
\[
|\Obs(G;\{u,v,w\})|=43.
\]
Using the subset-sum formula,
\[
\mathcal E(G;\{u,v,w\},q)
=\sum_{X\subseteq\{u,v,w\}} |\Obs(G;X)|\;q^{3-|X|}(1-q)^{|X|}
=43(1-q).
\]
The example demonstrates that nontrivial cancellation among subset contributions can compress a high-dimensional family of classical observed-set sizes into a linear expectation polynomial [2507.14620].

In application terms, the fragile model changes the optimization target. In real power-grid monitoring, one places a minimum number of PMUs so that, in the absence of failure, the grid is fully observed, as measured by \(\gamma_P\). Fragile power domination replaces the “worst-case \(k\) faulty sensors” assumption by a stochastic failure probability \(q\) and asks instead for the expected coverage \(\mathcal E(G;S,q)\). One can then optimize sensor placement to maximize expected coverage, or ensure \(\mathcal E(G;S,q)\equiv |V|\) up to a given \(q\). The algebraic characterization further shows that only the multiset of values
\[
\{\,|\Obs(G;W)|:W\subseteq S,\ |W|=k\}
\]
for \(0\le k\le |S|\) matters, which can guide heuristics and integer-programming formulations. Gadget constructions imply that one can inflate or deflate intermediate polynomial coefficients, which is useful in designing hard instances or in reverse-engineering required robustness profiles [2507.14620].

The extensions identified in the sources are broad but technically coherent. They include non-uniform sensor-failure probabilities \(q_v\), time-dependent propagations, continuous-time failure models, and combinations with classical \(k\)-fault-tolerant power domination [2507.14620]. Closely related open directions include weighted or colored graphs where PMU reliability \(q(v)\) depends on the vertex, dynamic failures that occur in stages during the zero-forcing process, sharper formulas for Cartesian-product graphs such as grids and tori, random graph models, structural parameters such as degree sequences and connectivity that control the shape of \(\mathcal E(G;S,q)\), and approximation algorithms for selecting \(S\) of fixed size to maximize expected coverage [2312.12259]. Together these directions indicate that fragile power domination is both a robustness formalism for PMU placement and a polynomial-valued graph invariant with substantial combinatorial structure.

Source: https://www.emergentmind.com/topics/fragile-power-domination