This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Propagation time for zero forcing on a graph

Leslie Hogben Department of Mathematics, Iowa State University, Ames, IA 50011, USA (lhogben@iastate.edu) and American Institute of Mathematics, 360 Portage Ave, Palo Alto, CA 94306, USA (hogben@aimath.org). My Huynh Department of Mathematics, Arizona State University, Tempe, AZ 85287 (mthuynh1@asu.edu). Research supported by DMS 0502354 and DMS 0750986. Nicole Kingsley Department of Mathematics, Iowa State University, Ames, IA 50011, USA (nkingsle@iastate.edu). Sarah Meyer Department of Mathematics, Smith College, Northampton, MA 01063, USA (smeyer@smith.edu). Research supported by DMS 0750986. Shanise Walker Department of Mathematics, University of Georgia, Athens, GA 30602, USA (shanise1@uga.edu). Research supported by DMS 0502354 and DMS 0750986 Michael Young Department of Mathematics, Iowa State University, Ames, IA 50011, USA (myoung@iastate.edu). Research supported by DMS 0946431.
Abstract

Zero forcing (also called graph infection) on a simple, undirected graph GG is based on the color-change rule: If each vertex of GG is colored either white or black, and vertex vv is a black vertex with only one white neighbor ww, then change the color of ww to black. A minimum zero forcing set is a set of black vertices of minimum cardinality that can color the entire graph black using the color change rule. The propagation time of a zero forcing set BB of graph GG is the minimum number of steps that it takes to force all the vertices of GG black, starting with the vertices in BB black and performing independent forces simultaneously. The minimum and maximum propagation times of a graph are taken over all minimum zero forcing sets of the graph. It is shown that a connected graph of order at least two has more than one minimum zero forcing set realizing minimum propagation time. Graphs GG having extreme minimum propagation times |G|1|G|-1, |G|2|G|-2, and 0 are characterized, and results regarding graphs having minimum propagation time 11 are established. It is shown that the diameter is an upper bound for maximum propagation time for a tree, but in general propagation time and diameter of a graph are not comparable.

Keywords zero forcing number, propagation time, graph

AMS subject classification 05C50, 05C12, 05C15, 05C57, 81Q93, 82B20, 82C20

1 Propagation time

All graphs are simple, finite, and undirected. In a graph GG where some vertices are colored black and the remaining vertices are colored white, the color change rule is: If vv is black and ww is the only white neighbor of vv, then change the color of ww to black; if we apply the color change rule to vv to change the color of ww, we say vv forces ww and write vwv\to w (note that there may be a choice involved, since as we record forces, only one vertex actually forces ww, but more than one may be able to). Given an initial set BB of black vertices, the final coloring of BB is the set of black vertices that results from applying the color change rule until no more changes are possible. For a given graph GG and set of vertices BB, the final coloring is unique [1]. A zero forcing set is an initial set BB of vertices such that the final coloring of BB is V(G)V(G). A minimum zero forcing set of a graph GG is a zero forcing set of GG of minimum cardinality, and the zero forcing number, denoted Z(G)\operatorname{Z}(G), is the cardinality of a minimum zero forcing set.

Zero forcing, also known as graph infection or graph propagation, was introduced independently in [1] for study of minimum rank problems in combinatorial matrix theory, and in [3] for study of control of quantum systems. Propagation time of a zero forcing set, which describes the number of steps needed to fully color a graph performing independent forces simultaneously, was implicit in [3] and explicit in [8]. Recently Chilakamarri et al. [4] determined the propagation time, which they call the iteration index, for a number of families of graphs including Cartesian products and various grid graphs. Control of an entire network by sequential operations on a subset of particles is valuable [8] and the number of steps needed to obtain this control (propagation time) is a significant part of the process. In this paper we systematically study propagation time.

Definition 1.1.

Let G=(V,E)G=(V,E) be a graph and BB a zero forcing set of GG. Define B(0)=BB^{(0)}=B, and for t0t\geq 0, B(t+1)B^{(t+1)} is the set of vertices ww for which there exists a vertex bs=0tB(s)b\in\bigcup_{s=0}^{t}B^{(s)} such that ww is the only neighbor of bb not in s=0tB(s)\bigcup_{s=0}^{t}B^{(s)}. The propagation time of BB in GG, denoted pt(G,B)\operatorname{pt}(G,B), is the smallest integer t0t_{0} such that V=t=0t0B(t)V=\bigcup_{t=0}^{t_{0}}B^{(t)}.

Two minimum zero forcing sets of the same graph may have different propagation times, as the following example illustrates.

Refer to caption
Figure 1: The graph GG for Example 1.2
Example 1.2.

Let GG be the graph in Figure 1. Let B1={g,h}B_{1}=\{g,h\} and B2={a,d}B_{2}=\{a,d\}. Then B1(1)={f,c}{B_{1}}^{(1)}=\{f,c\}, B1(2)={e}{B_{1}}^{(2)}=\{e\}, B1(3)={d}{B_{1}}^{(3)}=\{d\}, B1(4)={b}{B_{1}}^{(4)}=\{b\}, and B1(5)={a}{B_{1}}^{(5)}=\{a\}, so pt(G,B1)=5\operatorname{pt}(G,B_{1})=5. However, B2(1)={b}{B_{2}}^{(1)}=\{b\}, B2(2)={c}{B_{2}}^{(2)}=\{c\}, B2(3)={e,h}{B_{2}}^{(3)}=\{e,h\}, and B2(4)={f,g}{B_{2}}^{(4)}=\{f,g\}, so pt(G,B2)=4\operatorname{pt}(G,B_{2})=4.

Definition 1.3.

The minimum propagation time of GG is

pt(G)=min{pt(G,B)|B is a minimum zero forcing set of G}.\operatorname{pt}(G)=\min\{\operatorname{pt}(G,B)\,|\,B\text{ is a minimum zero forcing set of }G\}.
Definition 1.4.

Two minimum zero forcing sets B1B_{1} and B2B_{2} of a graph GG are isomorphic if there is a graph automorphism φ\varphi of GG such that φ(B1)=B2\varphi(B_{1})=B_{2}.

It is obvious that isomorphic zero forcing sets have the same propagation time, but a graph may have non-isomorphic minimum zero forcing sets and have the property that all minimum zero forcing sets have the same propagation time.

Example 1.5.

Up to isomorphism, the minimum zero forcing sets of the dart shown in Figure 2 are {a,c}\{a,c\}, {b,c},\{b,c\}, and {c,d}\{c,d\}. Each of these sets has propagation time 33.

Refer to caption
Figure 2: The dart

The minimum propagation time of a graph GG is not subgraph monotone. For example, it is easy to see that the 4-cycle has Z(C4)=2\operatorname{Z}(C_{4})=2 and pt(C4)=1\operatorname{pt}(C_{4})=1. By deleting one edge of C4C_{4}, it becomes a path P3P_{3}, which has Z(P3)=1\operatorname{Z}(P_{3})=1 and pt(P3)=2\operatorname{pt}(P_{3})=2.

A minimum zero forcing set that achieves minimum propagation time plays a central role in our study, and we name such a set.

Definition 1.6.

A subset BB of vertices of GG is an efficient zero forcing set for GG if BB is a minimum zero forcing set of GG and pt(G,B)=pt(G)\operatorname{pt}(G,B)=\operatorname{pt}(G). Define

Eff(G)={B|B is an efficient zero forcing set of G}.\operatorname{Eff}(G)=\{B\,|\,B\text{ is an efficient zero forcing set of }G\}.

We can also consider maximum propagation time.

Definition 1.7.

The maximum propagation time of GG is defined as

PT(G)=max{pt(G,B)|B is a minimum zero forcing set of G}.\operatorname{PT}(G)=\max\{\operatorname{pt}(G,B)\,|\,B\text{ is a minimum zero forcing set of }G\}.

The bounds in the next remark were also observed in [4].

Remark 1.8.

Let GG be a graph. Then

PT(G)|G|Z(G)\operatorname{PT}(G)\leq|G|-\operatorname{Z}(G)

because at least one force must be performed at each time step, and

|G|Z(G)Z(G)pt(G)\frac{|G|-\operatorname{Z}(G)}{\operatorname{Z}(G)}\leq\operatorname{pt}(G)

because using a given zero forcing set BB, at most |B||B| forces can be performed at any one time step.

Definition 1.9.

The propagation time interval of GG is defined as

[pt(G),PT(G)]={pt(G),pt(G)+1,,PT(G)1,PT(G)}.[\operatorname{pt}(G),\operatorname{PT}(G)]=\{\operatorname{pt}(G),\operatorname{pt}(G)+1,\dots,\operatorname{PT}(G)-1,\operatorname{PT}(G)\}.

The propagation time discrepancy of GG is defined as

pd(G)=PT(G)pt(G).\operatorname{pd}(G)=\operatorname{PT}(G)-\operatorname{pt}(G).

It is not the case that every integer in the propagation time interval is the propagation time of a minimum zero forcing set; this can be seen in the next example. Let S(e1,e2,e3)S(e_{1},e_{2},e_{3}) be the generalized star with three arms having e1,e2,e3e_{1},e_{2},e_{3} vertices with e1e2e3e_{1}\leq e_{2}\leq e_{3}; S(2,5,11)S(2,5,11) is shown in Figure 3.

Refer to caption
Figure 3: The tree S(2,5,11)S(2,5,11)
Example 1.10.

Consider S(e1,e2,e3)S(e_{1},e_{2},e_{3}) with 1<e1<e2<e31<e_{1}<e_{2}<e_{3}. The vertices of degree one are denoted by u1,u2,u3u_{1},u_{2},u_{3}, the vertex of degree three is denoted by vv, and neighbors of vv are denoted by w1,w2,w3w_{1},w_{2},w_{3}. The minimum zero forcing sets and their propagation times are shown in Table 1. Observe that the propagation time interval of S(2,5,11)S(2,5,11) is [12,16][12,16], but there is no minimum zero forcing set with propagation time 1414. The propagation discrepancy is pd(S(2,5,11))=4\operatorname{pd}(S(2,5,11))=4.

Table 1: Minimum zero forcing sets and propagation times of S(e1,e2,e3)S(e_{1},e_{2},e_{3})
BB pt(S(e1,e2,e3),B)\operatorname{pt}(S(e_{1},e_{2},e_{3}),B) pt(S(2,5,11),B)\operatorname{pt}(S(2,5,11),B)
{u1,u2}\{u_{1},u_{2}\} e2+e31e_{2}+e_{3}-1 1515
{u3,w2}\{u_{3},w_{2}\} e2+e31e_{2}+e_{3}-1 1515
{u3,w1}\{u_{3},w_{1}\} e2+e3e_{2}+e_{3} 1616
{u1,u3}\{u_{1},u_{3}\} e2+e31e_{2}+e_{3}-1 1515
{u2,w3}\{u_{2},w_{3}\} e2+e31e_{2}+e_{3}-1 1515
{u2,w1}\{u_{2},w_{1}\} e2+e3e_{2}+e_{3} 1616
{u2,u3}\{u_{2},u_{3}\} e1+e31e_{1}+e_{3}-1 1212
{u1,w3}\{u_{1},w_{3}\} e1+e31e_{1}+e_{3}-1 1212
{u1,w2}\{u_{1},w_{2}\} e1+e3e_{1}+e_{3} 1313

The next remark provides a necessary condition for a graph GG to have pd(G)=0\operatorname{pd}(G)=0.

Remark 1.11.

Let GG be a graph. Then every minimum zero forcing set of GG is an efficient zero forcing set if and only if pd(G)=0\operatorname{pd}(G)=0. In [2], it is proven that the intersection of all minimum zero forcing sets is the empty set. Hence, pd(G)=0\operatorname{pd}(G)=0 implies BEff(G)B=\displaystyle\bigcap_{B\in\operatorname{Eff}(G)}B=\emptyset.

In Section 2 we establish properties of efficient zero forcing sets, including that no connected graph of order more than one has a unique efficient zero forcing set. In Section 3 we characterize graphs having extreme propagation time. In Section 4 we examine the relationship between propagation time and diameter.

2 Efficient zero forcing sets

In [2] it was shown that for a connected graph of order at least two, there must be more than one minimum zero forcing set and furthermore, no vertex is in every minimum zero forcing set. This raises the questions of whether the analogous properties are true for efficient zero forcing sets (Questions 2.1 and 2.14 below).

Question 2.1.

Is there a connected graph of order at least two that has a unique efficient zero forcing set?

We show that the answer to Question 2.1 is negative. First we need some terminology. For a given zero forcing set BB of GG, construct the final coloring, listing the forces in the order in which they were performed. This list is a chronological list of forces of BB [6]. Many definitions and results concerning lists of forces that have appeared in the literature involve chronological (ordered) lists of forces. For the study of propagation time, the order of forces is often dictated by performing a force as soon as possible (propagating). Thus unordered sets of forces are more useful than ordered lists when studying propagation time, and we extend terminology from chronological lists of forces to sets of forces.

Definition 2.2.

Let G=(V,E)G=(V,E) be a graph, BB a zero forcing set of GG. The unordered set of forces in a chronological list of forces of BB is called a set of forces of BB.

Observe that if BB is a zero forcing set and \mathcal{F} is a set of forces of BB, then the cardinality of \mathcal{F} is |G||B||G|-|B|. The ideas of terminus and reverse set of forces, introduced in [2] for a chronological list of forces and defined below for a set of forces, are used to answer Question 2.1 negatively (by constructing the terminus of a set of forces of an efficient zero forcing set).

Definition 2.3.

Let GG be a graph, let BB be a zero forcing set of GG, and let \mathcal{F} be a set of forces of BB. The terminus of \mathcal{F}, denoted Term()\operatorname{Term}(\mathcal{F}), is the set of vertices that do not perform a force in \mathcal{F}. The reverse set of forces of \mathcal{F}, denoted here as Rev()\operatorname{Rev}(\mathcal{F}), is the result of reversing each force in \mathcal{F}. A forcing chain of \mathcal{F} is a sequence of vertices (v1,v2,,vk)(v_{1},v_{2},\dots,v_{k}) such that for i=1,,k1i=1,\dots,k-1, viv_{i} forces vi+1v_{i+1} in \mathcal{F} (k=1k=1 is permitted). A maximal forcing chain is a forcing chain that is not a proper subset of another forcing chain.

The name “terminus” reflects the fact that a vertex does not perform a force in \mathcal{F} if and only if it is the end point of a maximal forcing chain (the latter is the definition used in [2], where such a set is called a reversal of BB). In [2], it is shown that if BB is a zero forcing set of GG and \mathcal{F} is a chronological list of forces, then the terminus of \mathcal{F} is also a zero forcing set of GG, with the reverse chronological list of forces (to construct a reverse chronological list of forces of \mathcal{F}, write the chronological list of forces in reverse order and reverse each force in \mathcal{F}).

Observation 2.4.

Let GG be a graph, BB a minimum zero forcing set of GG, and \mathcal{F} a set of forces of BB. Then Rev()\operatorname{Rev}(\mathcal{F}) is a set of forces of Term()\operatorname{Term}(\mathcal{F}) and B=Term(Rev())B=\operatorname{Term}(\operatorname{Rev}(\mathcal{F})).

When studying propagation time, it is natural to examine sets of forces that achieve minimum propagation time.

Definition 2.5.

Let G=(V,E)G=(V,E) be a graph and BB a zero forcing set of GG. For a set of forces \mathcal{F} of BB, define (0)=B\mathcal{F}^{(0)}=B and for t0t\geq 0, (t+1)\mathcal{F}^{(t+1)} is the set of vertices ww such that the force vwv\to w appears in \mathcal{F}, wi=0t(i)w\notin\bigcup_{i=0}^{t}\mathcal{F}^{(i)}, and ww is the only neighbor of vv not in i=0t(i)\bigcup_{i=0}^{t}\mathcal{F}^{(i)}. The propagation time of \mathcal{F} in GG, denoted pt(G,)\operatorname{pt}(G,\mathcal{F}), is the least t0t_{0} such that V=t=0t0(t)V=\bigcup_{t=0}^{t_{0}}\mathcal{F}^{(t)}.

Let G=(V,E)G=(V,E) be a graph, let BB be a zero forcing set of GG, and let \mathcal{F} be a set of forces of BB. Clearly, i=0t(i)i=0tB(i)\bigcup_{i=0}^{t}\mathcal{F}^{(i)}\subseteq\bigcup_{i=0}^{t}B^{(i)} for all t=0,,pt(G,B)t=0,\dots,\operatorname{pt}(G,B).

Definition 2.6.

Let G=(V,E)G=(V,E) be a graph and let BB be a zero forcing set of GG. A set of forces \mathcal{F} is efficient if pt(G,)=pt(G)\operatorname{pt}(G,\mathcal{F})=\operatorname{pt}(G). Define

eff(G)={| is an efficient set of forces of a minimum zero forcing set B of G}.\mathcal{F}_{eff}(G)=\{\mathcal{F}\,|\,\mathcal{F}\text{ is an efficient set of forces of a minimum zero forcing set $B$ of }G\}.

If \mathcal{F} is an efficient set of forces of a minimum zero forcing set BB of GG, then BB is necessarily an efficient zero forcing set. However, not every efficient set of forces conforms to the propagation process.

Refer to caption
Figure 4: The graph GG for Example 2.7
Example 2.7.

Let GG be the graph in Figure 4. Since every degree one vertex must be an endpoint of a maximal forcing chain and since B={x,z,v}B=\{x,z,v\} is a zero forcing set, Z(G)=3\operatorname{Z}(G)=3. Since (v,c,d,e,f,w)(v,c,d,e,f,w) or (w,f,e,d,c,v)(w,f,e,d,c,v) must be a maximal forcing chain for any set of forces of a minimum zero forcing set, pt(G)=5\operatorname{pt}(G)=5. Then BB is an efficient zero forcing set with efficient set of forces ={vc,za,cd,ab,de,by,ef,fw}.\mathcal{F}=\{v\to c,z\to a,c\to d,a\to b,d\to e,b\to y,e\to f,f\to w\}. Observe that b(2)b\in\mathcal{F}^{(2)} (i.e., bb does not turn black until step t=2t=2 in \mathcal{F}), but bB(1)b\in B^{(1)} (bb can be forced by xx at step t=1t=1).

Definition 2.8.

Let G=(V,E)G=(V,E) be a graph, BB a zero forcing set of GG, and \mathcal{F} a set of forces of BB. Define Q0()=Term()Q_{0}(\mathcal{F})=\operatorname{Term}(\mathcal{F}) and for t=1,,pt(G,)t=1,\ldots,\operatorname{pt}(G,\mathcal{F}), define

Qt()={vV|w(pt(G,)t+1) such that vw}.Q_{t}(\mathcal{F})=\{v\in V\,|\,\exists w\ \in\mathcal{F}^{(\operatorname{pt}(G,\mathcal{F})-t+1)}\mbox{ such that }v\to w\}.

Observe that V=t=0pt(G,)Qt()V=\bigcup_{t=0}^{\operatorname{pt}(G,\mathcal{F})}Q_{t}(\mathcal{F}).

Lemma 2.9.

Let G=(V,E)G=(V,E) be a graph, BB a zero forcing set of GG, and \mathcal{F} a set of forces of BB. Then Qt()i=0tRev()(i)Q_{t}(\mathcal{F})\subseteq\bigcup_{i=0}^{t}{\operatorname{Rev}(\mathcal{F})}^{(i)}.

Proof.

Recall that Rev()\operatorname{Rev}(\mathcal{F}) is a set of forces of Term()\operatorname{Term}(\mathcal{F}). The result is established by induction on tt. Initially, Q0()=Term()=Rev()(0)Q_{0}(\mathcal{F})=\operatorname{Term}(\mathcal{F})=\operatorname{Rev}(\mathcal{F})^{(0)}. Assume that for 0st0\leq s\leq t, Qs()i=0sRev()(i)Q_{s}(\mathcal{F})\subseteq\bigcup_{i=0}^{s}{\operatorname{Rev}(\mathcal{F})}^{(i)}. Let vQt+1()v\in Q_{t+1}(\mathcal{F}). In \mathcal{F}, vuv\to u at time pt(G,)t\operatorname{pt}(G,\mathcal{F})-t. In \mathcal{F}, uu cannot perform a force until time pt(G,B)t+1\operatorname{pt}(G,B)-t+1 or later, so ui=0tQi()i=0tRev()(i)u\in\bigcup_{i=0}^{t}Q_{i}(\mathcal{F})\subseteq\bigcup_{i=0}^{t}{\operatorname{Rev}(\mathcal{F})}^{(i)}. If xN(u){v}x\in N(u)\setminus\{v\} then in \mathcal{F} xx cannot perform a force before time pt(G,)t+1\operatorname{pt}(G,\mathcal{F})-t+1, so xi=0tQi()i=0tRev()(i)x\in\bigcup_{i=0}^{t}Q_{i}(\mathcal{F})\subseteq\bigcup_{i=0}^{t}{\operatorname{Rev}(\mathcal{F})}^{(i)}. So if vi=0tRev()(i)v\notin\bigcup_{i=0}^{t}{\operatorname{Rev}(\mathcal{F})}^{(i)}, then vRev()(t+1)v\in{\operatorname{Rev}(\mathcal{F})}^{(t+1)}. Thus vi=0t+1Rev()(i)v\in\bigcup_{i=0}^{t+1}{\operatorname{Rev}(\mathcal{F})}^{(i)}. ∎

Corollary 2.10.

Let G=(V,E)G=(V,E) be a graph, BB a minimum zero forcing set of GG, and \mathcal{F} a set of forces of BB. Then

pt(G,Rev())pt(G,).\operatorname{pt}(G,\operatorname{Rev}(\mathcal{F}))\leq\operatorname{pt}(G,\mathcal{F}).

The next result follows from Corollary 2.10 and Observation 2.4.

Theorem 2.11.

Let G=(V,E)G=(V,E) be a graph, BB an efficient zero forcing set of GG, and \mathcal{F} an efficient set of forces of BB. Then Rev()\operatorname{Rev}(\mathcal{F}) is an efficient set of forces and Term()\operatorname{Term}(\mathcal{F}) is an efficient zero forcing set. Every efficient zero forcing set is the terminus of an efficient set of forces of an efficient zero forcing set.

The next result answers Question 2.1 negatively.

Theorem 2.12.

Let GG be a connected graph of order greater than one. Then |Eff(G)|2|\operatorname{Eff}(G)|\geq 2.

Proof.

Let BEff(G)B\in\operatorname{Eff}(G) and let \mathcal{F} be an efficient set of forces of BB. By Theorem 2.11, Term(B)Eff(G)\operatorname{Term}(B)\in\operatorname{Eff}(G). Since GG is a connected graph of order greater than one, BTerm()B\neq\operatorname{Term}(\mathcal{F}). ∎

We now consider the intersection of efficient zero forcing sets. The next result is immediate from Theorem 2.11.

Corollary 2.13.

Let GG be a graph. Then BEff(G)B=eff(G)Term()\displaystyle\bigcap_{B\in\operatorname{Eff}(G)}B=\bigcap_{\mathcal{F}\in\mathcal{F}_{eff}(G)}\operatorname{Term}(\mathcal{F}).

Question 2.14.

Is there a connected graph GG of order at least two and a vertex vV(G)v\in V(G) such that vv is in every efficient zero forcing set?

The next example provides an affirmative answer.

Example 2.15.

The wheel W5W_{5} is the graph shown in Figure 5. The efficient zero forcing set {a,b,c}\{a,b,c\} of W5W_{5} shows that pt(W5)=1\operatorname{pt}(W_{5})=1. Up to isomorphism, there are two types of minimum zero forcing sets in W5W_{5}. One set contains aa and two other vertices that are adjacent to each other; the other contains three vertices other than aa. The latter is not an efficient zero forcing set of W5W_{5}, because its propagation time is 22. The possible choices for an efficient zero forcing set are {a,b,c},{a,c,d},{a,d,e},\{a,b,c\},\{a,c,d\},\{a,d,e\}, or {a,b,e}\{a,b,e\}. Therefore, BEff(G)B={a}\displaystyle\bigcap_{B\in\operatorname{Eff}(G)}B=\{a\}.

Refer to caption
Figure 5: The wheel W5W_{5}

We examine the effect of a nonforcing vertex in an efficient zero forcing set. This result will be used in Section 3.2. It was shown in [5] that Z(G)=Z(Gv)+1\operatorname{Z}(G)=\operatorname{Z}(G-v)+1 if and only if there exists a minimum zero forcing set BB containing vv and set of forces \mathcal{F} in which vv does not perform a force. The proof of the next proposition is the same but with consideration restricted to an efficient zero forcing set (the idea is that the same set of forces works for both GG and GvG-v, with vv included in the zero forcing set for GG).

Proposition 2.16.

For a vertex vv of a graph GG, there exists an efficient zero forcing set BB containing vv and an efficient set of forces \mathcal{F} in which vv does not perform a force if and only if pt(Gv)=pt(G)\operatorname{pt}(G-v)=\operatorname{pt}(G) and Z(Gv)=Z(G)1\operatorname{Z}(G-v)=\operatorname{Z}(G)-1.

3 Graphs with extreme minimum propagation time

For any graph GG, it is clear that 0pt(G)PT(G)|G|1.0\leq\operatorname{pt}(G)\leq\operatorname{PT}(G)\leq|G|-1. In this section we consider the extreme values |G|1|G|-1, |G|2|G|-2, i.e., high propagation time, and, 0 and 11, i.e., low propagation time.

3.1 High propagation time

The case of propagation time |G|1|G|-1 is straightforward, using the well known fact [7] that Z(G)=1\operatorname{Z}(G)=1 if and only if GG is a path.

Proposition 3.1.

For a graph GG, the following are equivalent.

  1. 1.

    pt(G)=|G|1\operatorname{pt}(G)=|G|-1.

  2. 2.

    PT(G)=|G|1\operatorname{PT}(G)=|G|-1.

  3. 3.

    Z(G)=1\operatorname{Z}(G)=1.

  4. 4.

    GG is a path.

We now consider graphs GG that have maximum or minimum propagation time equal to |G|2|G|-2.

Observation 3.2.

For a graph GG,

  1. 1.

    pt(G)=|G|2\operatorname{pt}(G)=|G|-2 implies PT(G)=|G|2\operatorname{PT}(G)=|G|-2, but not conversely (see Lemma 3.4 for an example).

  2. 2.

    pt(G)=|G|2\operatorname{pt}(G)=|G|-2 if and only if Z(G)=2\operatorname{Z}(G)=2 and exactly one force is performed at each time for every minimum zero forcing set.

  3. 3.

    PT(G)=|G|2\operatorname{PT}(G)=|G|-2 if and only if Z(G)=2\operatorname{Z}(G)=2 and there exists a minimum zero forcing set such that exactly one force performed at each time.

Lemma 3.3.

Let GG be a disconnected graph. Then the following are equivalent.

  1. 1.

    pt(G)=|G|2\operatorname{pt}(G)=|G|-2.

  2. 2.

    PT(G)=|G|2\operatorname{PT}(G)=|G|-2.

  3. 3.

    G=Pn1˙P1G=P_{n-1}\dot{\cup}P_{1}.

Proof.

Clearly G=Pn1˙P1pt(G)=|G|2PT(G)=|G|2G=P_{n-1}\dot{\cup}P_{1}\Rightarrow\operatorname{pt}(G)=|G|-2\Rightarrow\operatorname{PT}(G)=|G|-2. So assume PT(G)=|G|2\operatorname{PT}(G)=|G|-2. Since Z(G)=2\operatorname{Z}(G)=2, GG has exactly two components. At least one component of GG is an isolated vertex (otherwise, more than one force occurs at time step one), and so G=Pn1˙P1G=P_{n-1}\dot{\cup}P_{1}. ∎

A path cover of a graph GG is a set of vertex disjoint induced paths that cover all the vertices of GG, and the path cover number P(G)\operatorname{P}(G) is the minimum number of paths in a path cover of GG. It is known [6] that for a given zero forcing set and set of forces, the set of maximal zero forcing chains forms a path cover and thus P(G)Z(G)\operatorname{P}(G)\leq\operatorname{Z}(G)

A graph GG is a graph on two parallel paths if V(G)V(G) can be partitioned into disjoint subsets U1U_{1} and U2U_{2} so that the induced subgraphs Pi=G[Ui],i=1,2P_{i}=G[U_{i}],i=1,2 are paths, GG can be drawn in the plane with the paths P1P_{1} and P2P_{2} as parallel line segments, and edges between the two paths (drawn as line segments, not curves) do not cross; such a drawing is called a standard drawing. The paths P1P_{1} and P2P_{2} are called the parallel paths (for this representation of GG as a graph on two parallel paths).

Let GG be a graph on two parallel paths P1P_{1} and P2P_{2}. If vV(G)v\in V(G), then path(v)\operatorname{path}(v) denotes the parallel path that contains vv and path¯(v)\overline{\operatorname{path}}(v) denotes the other of the parallel paths. Fix an ordering of the vertices in each of P1P_{1} and P2P_{2} that is increasing in the same direction for both paths in a standard drawing. With this ordering, let first(Pi)\operatorname{first}(P_{i}) and last(Pi)\operatorname{last}(P_{i}) denote the first and last vertices of Pi,i=1,2P_{i},i=1,2. If v,wV(Pi)v,w\in V(P_{i}), then vwv\prec w means vv precedes ww in the order on PiP_{i}. Furthermore, if vV(Pi)v\in V(P_{i}) and vlast(Pi)v\neq\operatorname{last}(P_{i}), next(v)\operatorname{next}(v) is the neighbor of vv in PiP_{i} such that vnext(v)v\prec\operatorname{next}(v); prev(v)\operatorname{prev}(v) is defined analogously (for vfirst(Pi)v\neq\operatorname{first}(P_{i})).

Row [7] has shown that Z(G)=2Z(G)=2 if and only if GG is a graph on two parallel paths. Observe that for any graph having Z(G)=2\operatorname{Z}(G)=2, a set of forces \mathcal{F} of a minimum zero forcing set naturally produces a representation of GG as a graph on two parallel paths with the parallel paths being the maximal forcing chains. The ordering of the vertices in the parallel paths is the forcing order.

Lemma 3.4.

For a tree GG, PT(G)=|G|2\operatorname{PT}(G)=|G|-2 if and only if G=S(1,1,n3)G=S(1,1,n-3) (sometimes called a T-shaped tree). The graph K1,3K_{1,3} is the only tree for which pt(G)=|G|2\operatorname{pt}(G)=|G|-2.

Proof.

It is clear that PT(S(1,1,n3))=n2\operatorname{PT}(S(1,1,n-3))=n-2, and pt(K1,3)=2\operatorname{pt}(K_{1,3})=2.

Suppose first that GG is a tree such that PT(G)=|G|2\operatorname{PT}(G)=|G|-2. Then GG is a graph on two parallel paths P1P_{1} and P2P_{2}. There is exactly one edge ee between the two paths. Observe that ee must have an endpoint not in {first(Pi),last(Pi),i=1,2}\{\operatorname{first}(P_{i}),\operatorname{last}(P_{i}),i=1,2\}, so without loss of generality first(P1)last(P1)\operatorname{first}(P_{1})\neq\operatorname{last}(P_{1}) and neither first(P1)\operatorname{first}(P_{1}) nor last(P1)\operatorname{last}(P_{1}) is an endpoint of ee. If GG is a graph with multiple vertices in each of P1,P2P_{1},P_{2} (i.e., if first(P2)last(P2)\operatorname{first}(P_{2})\neq\operatorname{last}(P_{2})), then no matter which minimum zero forcing set we choose, more than one force will occur at some time. So assume V(P2)V(P_{2}) consists of a single vertex ww. If the parallel paths were constructed from a minimum zero forcing set BB, then wBw\in B, and without loss of generality B={first(P1),w}B=\{\operatorname{first}(P_{1}),w\}. If N(w)N(first(P1))N(w)\neq N(\operatorname{first}(P_{1})), then at time one, two vertices would be forced. Thus N(w)=N(first(P1))N(w)=N(\operatorname{first}(P_{1})) and G=S(1,1,n3)G=S(1,1,n-3).

Now suppose that GG is a tree such that pt(G)=|G|2\operatorname{pt}(G)=|G|-2. This implies PT(G)=|G|2\operatorname{PT}(G)=|G|-2, so G=S(1,1,n3)G=S(1,1,n-3). Since n3>1n-3>1 implies pt(S(1,1,n3))<n2\operatorname{pt}(S(1,1,n-3))<n-2, G=S(1,1,1)=K1,3G=S(1,1,1)=K_{1,3}. ∎

Observation 3.5.

If GG is one of the graphs shown in Figure 6, then pt(G)<|G|2\operatorname{pt}(G)<|G|-2, because the black vertices are a minimum zero forcing set BB with pt(G,B)<|G|2\operatorname{pt}(G,B)<|G|-2.

Refer to caption
Figure 6: Graphs GG and minimum zero forcing sets BB such that pt(G,B)<|G|2\operatorname{pt}(G,B)<|G|-2 (where light vertices may be absent or repeated and similarly for light edges)

For any graph and vertices x,yx,y, xyx\sim y denotes that xx and yy are adjacent, and xyxy denotes the edge with endpoints xx and yy.

Definition 3.6.

A graph GG on two parallel paths P1P_{1} and P2P_{2} is a zigzag graph if it satisfies the following conditions:

  1. 1.

    There is a path Q=(z1,z2,,z)Q=(z_{1},z_{2},\dots,z_{\ell}) that alternates between the two paths P1P_{1} and P2P_{2} such that:

    1. (a)

      z2i1V(P1)z_{2i-1}\in V(P_{1}) and z2iV(P2)z_{2i}\in V(P_{2}) for i=1,,+12i=1,\dots,\lfloor\frac{\ell+1}{2}\rfloor;

    2. (b)

      zjzj+2z_{j}\prec z_{j+2} for j=1,,2j=1,\dots,\ell-2.

  2. 2.

    Every edge of GG is an edge of one of P1,P2P_{1},P_{2}, or QQ, or is an edge of the form

    zjw where 1<j<,wpath¯(zj), and zj1wzj+1.z_{j}w\mbox{ where }1<j<\ell,w\in\overline{\operatorname{path}}(z_{j}),\mbox{ and }z_{j-1}\prec w\prec z_{j+1}.

The number \ell of vertices in QQ is called the zigzag order.

Examples of zigzag graphs are shown in Figures 6, 7, and 8.

Refer to caption
Figure 7: A zigzag graph (with P1,P2P_{1},P_{2} and QQ in black)
Theorem 3.7.

Let GG be a graph. Then pt(G)=|G|2\operatorname{pt}(G)=|G|-2 if and only if GG is one of the following:

  1. 1.

    Pn1˙P1P_{n-1}\dot{\cup}P_{1}.

  2. 2.

    K1,3K_{1,3}.

  3. 3.

    A zigzag graph of zigzag order \ell such that all of the following conditions are satisfied:

    1. (a)

      GG is not isomorphic to one of the graphs shown in Figure 6.

    2. (b)

      deg(first(P1))>1\deg(\operatorname{first}(P_{1}))>1 or deg(first(P2))>1\deg(\operatorname{first}(P_{2}))>1 (both paths cannot begin with degree-one vertices).

    3. (c)

      deg(last(P1))>1\deg(\operatorname{last}(P_{1}))>1 or deg(last(P2))>1\deg(\operatorname{last}(P_{2}))>1 (both paths cannot end with degree-one vertices).

    4. (d)

      z2first(P2)z_{2}\neq\operatorname{first}(P_{2}) or z2next(z1)z_{2}\sim\operatorname{next}(z_{1})

    5. (e)

      z1last(path(z1))z_{\ell-1}\neq\operatorname{last}(\operatorname{path}(z_{\ell-1})) or z1prev(z)z_{\ell-1}\sim\operatorname{prev}(z_{\ell})

An example of a zigzag graph satisfying conditions (3a) – (3e) is shown in Figure 8.

Proof.

Assume pt(G)=|G|2\operatorname{pt}(G)=|G|-2. If GG is disconnected or a tree, then GG is Pn1˙P1P_{n-1}\dot{\cup}P_{1} or K1,3K_{1,3} by Lemmas 3.3 and 3.4. So assume GG is connected and has a cycle.

First we identify P1,P2P_{1},P_{2} and QQ for GG: By Observation 3.2, there exists a minimum zero forcing set BB of cardinality 2 such that exactly one force is performed at each time for BB. Renumber the vertices of GG as follows: vertices V(G)={1,0,1,2,,n2}V(G)=\{-1,0,1,2,\dots,n-2\}, zero forcing set B={1,0}B=\{-1,0\} with 010\to 1, and vertex tt is forced at time tt. Then GG is a graph on two parallel paths P1P_{1} and P2P_{2}, which are the two maximal forcing chains (with the path order being the forcing order). Observe that deg(0)2\deg(0)\leq 2 and deg(1)2\deg(-1)\geq 2, because 0 can immediately force and 1-1 cannot. If deg(1)=2\deg(-1)=2 and |G|>3|G|>3, then choose P1P_{1} to be path(1)\operatorname{path}(-1), and let z1=1,z2=maxN(1)P2z_{1}=-1,\ z_{2}=\max N(-1)\cap P_{2}. Otherwise, choose P1P_{1} to be path(0)\operatorname{path}(0) and let z1=minN(1),z2=1z_{1}=\min N(-1),\ z_{2}=-1. For j2j\geq 2, define zj+1=maxN(zj)path¯(zj)z_{j+1}=\max N(z_{j})\cap\overline{\operatorname{path}}(z_{j}) until N(zj)path¯(zj)=N(z_{j})\cap\overline{\operatorname{path}}(z_{j})=\emptyset. Define Q=(z1,,z)Q=(z_{1},\dots,z_{\ell}). With this labeling, GG is a zigzag graph.

Now we show that GG satisfies conditions (3a) – (3e). Since pt(G)=|G|2\operatorname{pt}(G)=|G|-2, GG is not isomorphic to one of the graphs shown in Figure 6, i.e., condition (3a) is satisfied. Since 1-1 is the first vertex in one of the paths and deg(1)2\deg(-1)\geq 2, condition (3b) is satisfied. The remaining conditions must be satisfied or there is a different zero forcing set of two vertices with lower propagation time: if (3c) fails, use B={last(P1),last(P2)}B=\{\operatorname{last}(P_{1}),\operatorname{last}(P_{2})\}; if (3d) fails, then z2=first(P2)z_{2}=\operatorname{first}(P_{2}) and z2≁next(z1)z_{2}\not\sim\operatorname{next}(z_{1}), so use B={first(P1),next(z1)}B=\{\operatorname{first}(P_{1}),\operatorname{next}(z_{1})\}; if (3e) fails, this is analogous to (3d) failing, so use B={last(path(z)),prev(z)}B=\{\operatorname{last}(\operatorname{path}(z_{\ell})),\operatorname{prev}(z_{\ell})\}.

For the converse, pt(G)=|G|2\operatorname{pt}(G)=|G|-2 for G=Pn1˙P1G=P_{n-1}\dot{\cup}P_{1} or G=K1,3G=K_{1,3} by Lemmas 3.3 and 3.4. So assume GG is a zigzag graph satisfying conditions (3a) – (3e). The sets B1={first(P1),first(P2)}B_{1}=\{\operatorname{first}(P_{1}),\operatorname{first}(P_{2})\} and B2={last(P1),last(P2)}B_{2}=\{\operatorname{last}(P_{1}),\operatorname{last}(P_{2})\} are minimum zero forcing sets of GG, and pt(G,Bi)=|G|2\operatorname{pt}(G,B_{i})=|G|-2 for i=1,2i=1,2. If z2=first(P2)z_{2}=\operatorname{first}(P_{2}), so first(P2)next(z1)\operatorname{first}(P_{2})\sim\operatorname{next}(z_{1}), then B3={first(P1),next(z1)}B_{3}=\{\operatorname{first}(P_{1}),\operatorname{next}(z_{1})\} is a zero forcing set and pt(G,B3)=|G|2\operatorname{pt}(G,B_{3})=|G|-2. If z1last(path(z1))z_{\ell-1}\neq\operatorname{last}(\operatorname{path}(z_{\ell-1})), so last(path(z1))prev(z)\operatorname{last}(\operatorname{path}(z_{\ell-1}))\sim\operatorname{prev}(z_{\ell}), then B4={last(path(z)),prev(z)}B_{4}=\{\operatorname{last}(\operatorname{path}(z_{\ell})),\operatorname{prev}(z_{\ell})\} is a zero forcing set and pt(G,B4)=|G|2\operatorname{pt}(G,B_{4})=|G|-2. If GG is not isomorphic to one of the graphs shown in Figure 6, these are the only minimum zero forcing sets. ∎

Refer to caption
Figure 8: A zigzag graph GG (with P1,P2P_{1},P_{2} and QQ in black) having pt(G)=|G|2\operatorname{pt}(G)=|G|-2

3.2 Low propagation time

Observation 3.8.

For a graph GG, the following are equivalent.

  1. 1.

    pt(G)=0\operatorname{pt}(G)=0.

  2. 2.

    PT(G)=0\operatorname{PT}(G)=0.

  3. 3.

    Z(G)=|G|\operatorname{Z}(G)=|G|.

  4. 4.

    GG has no edges.

Next we consider pt(G)=1\operatorname{pt}(G)=1. From Remark 1.8 we see that if pt(G)=1\operatorname{pt}(G)=1, then Z(G)|G|2\displaystyle\operatorname{Z}(G)\geq\frac{|G|}{2}. The converse of this statement is false:

Example 3.9.

Let GG be the graph obtained from K4K_{4} by appending a leaf to one vertex. Then Z(G)=3>|G|/2\operatorname{Z}(G)=3>|G|/2 and pt(G)=2\operatorname{pt}(G)=2.

Theorem 3.10.

Let G=(V,E)G=(V,E) be a connected graph such that pt(G)=1\operatorname{pt}(G)=1. For vVv\in V, vBEff(G)B\displaystyle v\in\bigcap_{B\in\operatorname{Eff}(G)}B if and only if for every BEff(G)B\in\operatorname{Eff}(G), |B(1)N(v)|2|B^{(1)}\cap N(v)|\geq 2.

Proof.

If for every BEff(G)B\in\operatorname{Eff}(G), |B(1)N(v)|2|B^{(1)}\cap N(v)|\geq 2, then vv cannot perform a force in an efficient set of forces, so vTerm()v\in\operatorname{Term}(\mathcal{F}) for every eff(G)\mathcal{F}\in\mathcal{F}_{eff}(G). Thus veff(G)Term()=BEff(G)Bv\in\bigcap_{\mathcal{F}\in\mathcal{F}_{eff}(G)}\operatorname{Term}(\mathcal{F})=\bigcap_{B\in\operatorname{Eff}(G)}B.

Now suppose vBEff(G)Bv\in\bigcap_{B\in\operatorname{Eff}(G)}B and let BEff(G)B\in\operatorname{Eff}(G). If vv performs a force in an efficient set of forces \mathcal{F} of BB, then vTerm()v\notin\operatorname{Term}(\mathcal{F}). By Theorem 2.11, Term()Eff(G)\operatorname{Term}(\mathcal{F})\in\operatorname{Eff}(G), so vv cannot perform a force in any such \mathcal{F}. Since vv cannot perform a force, |B(1)N(v)|1|B^{(1)}\cap N(v)|\neq 1. It is shown in [2] that (assuming the graph is connected and of order greater than one) every vertex of a minimum zero forcing set must have a neighbor not in the zero forcing set. Since pt(G)=1\operatorname{pt}(G)=1, |G|>1|G|>1, and so |B(1)N(v)|2|B^{(1)}\cap N(v)|\geq 2. ∎

We now consider the case of a graph GG that has pt(G)=1\operatorname{pt}(G)=1 and Z(G)=12|G|\operatorname{Z}(G)=\frac{1}{2}|G|. Examples of such graphs include the hypercubes QsQ_{s} [1].

Definition 3.11.

Suppose H1=(V1,E1)H_{1}=(V_{1},E_{1}) and H2=(V2,E2)H_{2}=(V_{2},E_{2}) are graphs of equal order and μ:V1V2\mu:V_{1}\to V_{2} is a bijection. Define the matching graph (H1,H2,μ)(H_{1},H_{2},\mu) to be the graph constructed as the disjoint union of H1,H2H_{1},H_{2} and the perfect matching between V1V_{1} and V2V_{2} defined by μ\mu.

Matching graphs play a central role in the study of graphs that have propagation time one.

Proposition 3.12.

Let G=(V,E)G=(V,E) be a graph. Then any two of the following conditions imply the third.

  1. 1.

    |G|=2Z(G)|G|=2\operatorname{Z}(G).

  2. 2.

    pt(G)=1\operatorname{pt}(G)=1.

  3. 3.

    GG is a matching graph

Proof.

(1)&(2)(3)(\ref{cZ})\ \&\ (\ref{cpt})\Rightarrow(\ref{cmatch}): Let BB be an efficient zero forcing set of GG and let B¯=VB\overline{B}=V\setminus B. Since |B|=12|G||B|=\frac{1}{2}|G| and pt(G)=1\operatorname{pt}(G)=1, every element bBb\in B must perform a force at time one. Thus |N(b)B¯|=1|N(b)\cap\overline{B}|=1 and there exists a perfect matching between BB and B¯\overline{B} defined by μ:BB¯\mu:B\to{\overline{B}} where μ(b)N(b)B¯\mu(b)\in N(b)\cap\overline{B}. Then G=(B,B¯,μ)G=(B,{\overline{B}},\mu).

For the remaining two parts, assume G=(H1,H2,μ)G=(H_{1},H_{2},\mu) and n=12|G|n=\frac{1}{2}|G| (=|H1|=|H2|=|H_{1}|=|H_{2}|).

(1)&(3)(2)(\ref{cZ})\ \&\ (\ref{cmatch})\Rightarrow(\ref{cpt}): Since Z(G)=n\operatorname{Z}(G)=n, H1H_{1} is a minimum zero forcing set and pt(G,H1)=1\operatorname{pt}(G,H_{1})=1.

(2)&(3)(1)(\ref{cpt})\ \&\ (\ref{cmatch})\Rightarrow(\ref{cZ}): Since pt(G)=1\operatorname{pt}(G)=1, Z(G)n\operatorname{Z}(G)\geq n, and Z(G)n\operatorname{Z}(G)\leq n because H1H_{1} is a zero forcing set with pt(G,H1)=1\operatorname{pt}(G,H_{1})=1. ∎

We examine conditions that ensure Z((H1,H2,μ))=|Hi|\operatorname{Z}((H_{1},H_{2},\mu))=|H_{i}| and thus pt((H1,H2,μ))=1\operatorname{pt}((H_{1},H_{2},\mu))=1. The choice of matching μ\mu affects the zero forcing number and propagation time, as the next two examples show.

Example 3.13.

The Cartesian product C5P2C_{5}\Box P_{2} is (C5,C5,ι)(C_{5},C_{5},\iota), where ι\iota is the identity mapping. It is known [1] that Z(C5P2)=4\operatorname{Z}(C_{5}\Box P_{2})=4 and thus pt(C5P2)>1\operatorname{pt}(C_{5}\Box P_{2})>1.

Example 3.14.

The Petersen graph PP can be constructed as (C5,C5,μP)(C_{5},C_{5},\mu_{P}) where μP=(1234514253)\mu_{P}=\left(\begin{array}[]{rrrrr}1&2&3&4&5\\ 1&4&2&5&3\end{array}\right). It is known [1] that Z(P)=5\operatorname{Z}(P)=5 and thus pt(P)=1\operatorname{pt}(P)=1.

Let c(G)c(G) denote the number of components of GG.

Theorem 3.15.

Let |H1|=|H2|=n|H_{1}|=|H_{2}|=n and let μ:H1H2\mu:H_{1}\to H_{2} be a bijection. If pt((H1,H2,μ))=1\operatorname{pt}((H_{1},H_{2},\mu))=1, then c(H1)=c(H2)=c((H1,H2,μ))c(H_{1})=c(H_{2})=c((H_{1},H_{2},\mu)).

Proof.

Assume it is not the case that c(H1)=c(H2)=c((H1,H2,μ))c(H_{1})=c(H_{2})=c((H_{1},H_{2},\mu)). This implies μ\mu is not the union of perfect matchings between the components of H1H_{1} and the components of H2H_{2}. Without loss of generality, there is a component H1[C1]H_{1}[C_{1}] of H1H_{1} that is not matched within a single component of H2H_{2}. Then there exist vertices uu and vv in C1C_{1} such that μ(v)Cv\mu(v)\in C_{v}, μ(u)Cu\mu(u)\in C_{u}, and H2[Cv]H_{2}[C_{v}] and H2[Cu]H_{2}[C_{u}] are separate components of H2H_{2}. We show that there is a zero forcing set of size n1n-1 for (H1,H2,μ)(H_{1},H_{2},\mu), and thus pt((H1,H2,μ))>1((H_{1},H_{2},\mu))>1. Let B1=C1\(μ1(Cv){u})B_{1}=C_{1}\backslash(\mu^{-1}(C_{v})\bigcup\{u\}), B2=V2\(μ(B1)μ(u))B_{2}=V_{2}\backslash(\mu(B_{1})\bigcup\mu(u)), and B=B1B2B=B_{1}\bigcup B_{2}, so |B|=n1|B|=n-1. Then 1) xμ1(x)x\rightarrow\mu^{-1}(x) for xCvx\in C_{v}, 2) vuv\rightarrow u, 3) yμ(y)y\rightarrow\mu(y) for yC1\μ1(Cv)y\in C_{1}\backslash\mu^{-1}(C_{v}), and 4) zμ1(z)z\rightarrow\mu^{-1}(z) for all μ1(z)\mu^{-1}(z) in the remaining components of H1H_{1}. Therefore B is a zero forcing set, Z((H1,H2,μ))n1\operatorname{Z}((H_{1},H_{2},\mu))\leq n-1, and thus pt((H1,H2,μ))>1\operatorname{pt}((H_{1},H_{2},\mu))>1. ∎

Theorem 3.16.

Let |H|=n|H|=n and let μ\mu be a bijection of vertices of HH and KnK_{n} (with μ\mu acting on the vertices of HH). Then pt((H,Kn,μ))=1\operatorname{pt}((H,K_{n},\mu))=1 if and only if HH is connected.

Proof.

If HH is not connected, then pt((H,Kn,μ))1\operatorname{pt}((H,K_{n},\mu))\neq 1 by Theorem 3.15. Now assume HH is connected and let G=(H,Kn,μ)G=(H,K_{n},\mu). Let BV(G)B\subseteq V(G) with |B|=n1|B|=n-1. We show BB is not a zero forcing set. This implies Z(G)=n\operatorname{Z}(G)=n and thus pt(G)=1\operatorname{pt}(G)=1. Let X=V(Kn)X=V(K_{n}) and Y=V(H)Y=V(H). For xXx\in X, xx cannot perform a force until at least n1n-1 vertices in XX are black. If |XB|=n1|X\cap B|=n-1 then no force can be performed. So assume |XB|n2|X\cap B|\leq n-2. Until at least n1n-1 vertices in XX are black, all forces must be performed by vertices in YY. We show that no more than n2n-2 vertices in XX can turn black. Perform all forces of the type yyy\to y^{\prime} with y,yYy,y^{\prime}\in Y. For each such force, μ(y)\mu(y) must be black already. Thus at most |XB||X\cap B| such forces within YY can be performed. So there are now at most |YB|+|XB|=n1|Y\cap B|+|X\cap B|=n-1 black vertices in YY. Note first that if at most n2n-2 vertices of YY are black, then after all possible forces from YY to XX are done, no further forces are possible, and at most n2n-2 vertices in XX are black. So assume n1n-1 vertices of YY are black. Let wYw\in Y be white. Since HH is connected, there must be a neighbor uu of ww in YY, and uu is black. Since uN(w)u\in N(w) and ww is white, uu has not performed a force. If μ(u)\mu(u) were black, there would be at most n2n-2 black vertices in YY, so μ(u)\mu(u) is white. After preforming all possible forces from YY to XX, at most n2n-2 vertices in XX are black because all originally black vertices xx have μ1(x)\mu^{-1}(x) black, there are n1n-1 black vertices in YY, and uu cannot perform a force at this time (since both ww and μ(u)\mu(u) are white). Thus not more than n2n-2 vertices of XX can be forced, and BB is not a zero forcing set. ∎

The Cartesian product of GG with P2P_{2} is one way of constructing matching graphs, because GP2=(G,G,ι)G\,\square\,P_{2}=(G,G,\iota). Examples of graphs GG having Z(GP2)=|G|\operatorname{Z}(G\,\square\,P_{2})=|G| include the complete graph KrK_{r} and hypercube QsQ_{s} [1]. Since Z(GP2)2Z(G)\operatorname{Z}(G\,\Box\,P_{2})\leq 2\operatorname{Z}(G) [1], to have Z(GP2)=|G|\operatorname{Z}(G\,\Box\,P_{2})=|G| it is necessary that Z(G)|G|2\operatorname{Z}(G)\geq\frac{|G|}{2}, but that condition is not sufficient.

Example 3.17.

Observe that Z(K1,r)=r112|K1,r|\operatorname{Z}(K_{1,r})=r-1\geq\frac{1}{2}|K_{1,r}| for r3r\geq 3. But Z(K1,rP2)=r<|K1,r|\operatorname{Z}(K_{1,r}\,\square\,P_{2})=r<|K_{1,r}|, so pt(K1,rP2)2\operatorname{pt}(K_{1,r}\,\square\,P_{2})\geq 2.

The next theorem provides conditions that ensure that iterating the Cartesian product with P2P_{2} gives a graph with propagation time one. Recall that one of the original motivations for defining the zero forcing number was to bound maximum nullity, and the interplay between these two parameters is central to the proof of the next theorem. Let G=({v1,,vn},E)G=(\{v_{1},\dots,v_{n}\},E) be a graph. The set of symmetric matrices described by GG is 𝒮(G)={An×n:AT=A and for ij,aij0vivjE}.\mathcal{S}(G)=\{A\in\mathbb{R}^{n\times n}:A^{T}=A\mbox{ and for }i\neq j,a_{ij}\neq 0\Leftrightarrow v_{i}v_{j}\in E\}. The maximum nullity of GG is M(G)=max{nullA:A𝒮(G)}.\operatorname{M}(G)=\max\{{\rm null}\,A:A\in\mathcal{S}(G)\}. It is well known [1] that M(G)Z(G)\operatorname{M}(G)\leq\operatorname{Z}(G). The next theorem provides conditions that are sufficient to iterate the construction of taking the Cartesian product of a graph and P2P_{2} and obtain minimum propagation time equal to one.

Theorem 3.18.

Suppose GG is a graph with |G|=n|G|=n and there exists a matrix L𝒮(G)L\in\mathcal{S}(G) such that L2=InL^{2}=I_{n}. Then

M(GP2)=Z(GP2)=n and pt(GP2)=1.\operatorname{M}(G\,\Box\,P_{2})=\operatorname{Z}(G\,\Box\,P_{2})=n\mbox{ and }\operatorname{pt}(G\,\Box\,P_{2})=1.

Furthermore, for

L^=12[LInInL]\hat{L}=\frac{1}{\sqrt{2}}\begin{bmatrix}L&I_{n}\\ I_{n}&-L\end{bmatrix}

L^𝒮(GP2)\hat{L}\in\mathcal{S}(G\,\Box\,P_{2}) and L^2=I2n{\hat{L}}^{2}=I_{2n}.

Proof.

Given the n×nn\times n matrix LL, define

H=[LInInL].H=\begin{bmatrix}L&I_{n}\\ I_{n}&L\end{bmatrix}.

Then H,L^𝒮(GP2)H,\hat{L}\in\mathcal{S}(G\,\Box\,P_{2}) and L^2=I2n{\hat{L}}^{2}=I_{2n}. Since [In0LIn][LInInL]=[LIn00]\begin{bmatrix}I_{n}&0\\ -L&I_{n}\end{bmatrix}\begin{bmatrix}L&I_{n}\\ I_{n}&L\end{bmatrix}=\begin{bmatrix}L&I_{n}\\ 0&0\end{bmatrix}, null(H)=n\operatorname{null}(H)=n. Therefore, M(GP2)n\operatorname{M}(G\,\Box\,P_{2})\geq n. Then

nM(GP2)Z(GP2)nn\leq\operatorname{M}(G\,\Box\,P_{2})\leq\operatorname{Z}(G\,\Box\,P_{2})\leq n

so we have equality throughout. Since GP2G\,\Box\,P_{2} is a matching graph, pt(GP2)=1\operatorname{pt}(G\,\Box\,P_{2})=1 by Proposition 3.12. ∎

Let G(P2)sG\,(\Box\,P_{2})^{s} denote the graph constructed by starting with GG and performing the Cartesian product with P2P_{2} ss times. For example, the hypercube Qs=P2(P2)s1Q_{s}=P_{2}\,(\Box\,P_{2})^{s-1}, and the proof given in [1] that M(Qs)=Z(Qs)=2s1\operatorname{M}(Q_{s})=\operatorname{Z}(Q_{s})=2^{s-1} is the same as the proof of Theorem 3.18 using the matrix L=[0110]𝒮(P2)L=\begin{bmatrix}0&1\\ 1&0\end{bmatrix}\in\mathcal{S}(P_{2}).

Corollary 3.19.

Suppose GG is a graph such that there exists a matrix L𝒮(G)L\in\mathcal{S}(G) such that L2=I|G|L^{2}=I_{|G|}. Then for s1s\geq 1,

M(G(P2)s)=Z(G(P2)s)=|G|2s1 and pt(G(P2)s)=1.\operatorname{M}(G\,(\Box\,P_{2})^{s})=\operatorname{Z}(G\,(\Box\,P_{2})^{s})=|G|2^{s-1}\mbox{ and }\operatorname{pt}(G\,(\Box\,P_{2})^{s})=1.
Corollary 3.20.

For s1s\geq 1, M(Kn(P2)s)=Z(Kn(P2)s)=n2s1\operatorname{M}(K_{n}\,(\Box\,P_{2})^{s})=\operatorname{Z}(K_{n}\,(\Box\,P_{2})^{s})=n2^{s-1} and pt(Kn(P2)s)=1\operatorname{pt}(K_{n}\,(\Box\,P_{2})^{s})=1.

Proof.

Let L=In2nJnL=I_{n}-\frac{2}{n}J_{n}, where JnJ_{n} is the n×nn\times n matrix having all entries equal to one. Then L𝒮(Kn)L\in\mathcal{S}(K_{n}) and L2=InL^{2}=I_{n}. ∎

Observe that if the matrix LL in the hypothesis of Theorem 3.18 is symmetric (and thus is an orthogonal matrix), then the matrix L^\hat{L} in the conclusion also has these properties. The same is true for Corollary 3.20, where L=In2nJnL=I_{n}-\frac{2}{n}J_{n} is a Householder transformation.

We have established a number of constructions that provide matching graphs having propagation time one. For example, pt(P2(P2)s)=1\operatorname{pt}(P_{2}\,(\Box\,P_{2})^{s})=1, pt(Kn(P2)s)=1\operatorname{pt}(K_{n}\,(\Box\,P_{2})^{s})=1, and if HH is connected, then pt((Kn,H,μ))=1\operatorname{pt}((K_{n},H,\mu))=1 for every matching μ\mu. But the general question remains open.

Question 3.21.

Characterize matching graphs (H1,H2,μ)(H_{1},H_{2},\mu) such that pt((H1,H2,μ))=1\operatorname{pt}((H_{1},H_{2},\mu))=1.

We can investigate when pt(G)=1\operatorname{pt}(G)=1 by deleting vertices that are in an efficient zero forcing set but do not perform a force in an efficient set of forces. The next result is a consequence of Proposition 2.16.

Corollary 3.22.

Let GG be a graph with pt(G)=1\operatorname{pt}(G)=1, BB an efficient zero forcing set of GG containing vv, and \mathcal{F} an efficient set of forces of BB in which vv does not perform a force. Then pt(Gv)=pt(G)=1\operatorname{pt}(G-v)=\operatorname{pt}(G)=1.

Definition 3.23.

Let GG be a graph with pt(G)=1\operatorname{pt}(G)=1, BB an efficient zero forcing set of GG, \mathcal{F} an efficient set of forces of BB, and SS the set of vertices in BB that do not perform a force. Define V=VSV^{\prime}=V\setminus S, G=G[V]=GSG^{\prime}=G[V^{\prime}]=G-S, B=BSB^{\prime}=B\setminus S, and B¯=VB\overline{B^{\prime}}=V^{\prime}\setminus B^{\prime}. The graph GG^{\prime} is called a prime subgraph of GG with associated zero forcing set BB^{\prime}.

Observation 3.24.

Let GG be a graph with pt(G)=1\operatorname{pt}(G)=1. For the prime subgraph GG^{\prime} and associated zero forcing set BB^{\prime} defined from an efficient zero forcing set BB and efficient set of forces \mathcal{F} of BB:

  1. 1.

    B¯=VB\overline{B^{\prime}}=V\setminus B.

  2. 2.

    |B|=|B¯||B|^{\prime}=|\overline{B^{\prime}}| and |G|=2|B||G^{\prime}|=2|B^{\prime}|.

  3. 3.

    GG^{\prime} is the matching graph defined by G[B],G[B¯]G[B^{\prime}],G[\overline{B^{\prime}}] and μ:BB¯\mu:B^{\prime}\to\overline{B^{\prime}} defined by μ(b)(N(b)B¯)\mu(b)\in(N(b)\cap\overline{B^{\prime}}).

  4. 4.

    BB^{\prime} and B¯{\overline{B}^{\prime}} are efficient zero forcing sets of GG^{\prime}.

  5. 5.

    pt(G)=1\operatorname{pt}(G^{\prime})=1.

It is clear that if G=(V,E)G=(V,E) has no isolated vertices, pt(G)=1\operatorname{pt}(G)=1, and if G^\hat{G} is constructed from GG by adjoining a new vertex vv adjacent to every uV(G)u\in V(G), then pt(G^)=1\operatorname{pt}(\hat{G})=1.

We now return to considering  BEff(G)B\bigcap_{B\in\operatorname{Eff}(G)}B, specifically in the case of propagation time one. We have a corollary of Theorem 3.10.

Corollary 3.25.

Let G=(V,E)G=(V,E) be a graph such that pt(G)=1\operatorname{pt}(G)=1. If vBEff(G)B\displaystyle v\in\bigcap_{B\in\operatorname{Eff}(G)}B, then degv4\deg v\geq 4.

Proof.

Let vBEff(G)B\displaystyle v\in\bigcap_{B\in\operatorname{Eff}(G)}B. Since pt(G)=1\operatorname{pt}(G)=1, for any efficient set \mathcal{F} of BB, Term(B)=B(1)\operatorname{Term}(B)=B^{(1)}. By Theorem 3.10, |B(1)N(v)|2|B^{(1)}\cap N(v)|\geq 2, so |Term(B)N(v)|2|\operatorname{Term}(B)\cap N(v)|\geq 2. Since B=Term(Rev())B=\operatorname{Term}(\operatorname{Rev}(\mathcal{F})), |BN(v)|2|{B}\cap N(v)|\geq 2. Thus degv4\deg v\geq 4. ∎

Note that Corollary 3.25 is false without the hypothesis that pt(G)=1\operatorname{pt}(G)=1, as the next example shows.

Example 3.26.

For the graph GG in Figure 9, Z(G)=3\operatorname{Z}(G)=3. Every minimum zero forcing set BB must contain one of {a,c}\{a,c\} and one of {x,z}\{x,z\}; without loss of generality, a,xBa,x\in B. If vBv\in B then pt(G,B)=2\operatorname{pt}(G,B)=2; if not then cc or zz is in BB and pt(G,B)=3\operatorname{pt}(G,B)=3. Thus vv is in every efficient zero forcing set.

Refer to caption
Figure 9: A graph GG with vBEff(G)Bv\in\bigcap_{B\in\operatorname{Eff}(G)}B and degv<4\deg v<4
Proposition 3.27.

Let G=(V,E)G=(V,E) be a graph and vVv\in V. If degv>Z(G)\deg v>\operatorname{Z}(G) and pt(G)=1\operatorname{pt}(G)=1, then vBEff(G)Bv\in\bigcap_{B\in\operatorname{Eff}(G)}B.

Proof.

Suppose vBEff(G)v\notin B\in\operatorname{Eff}(G), and let \mathcal{F} be an efficient set of forces of BB. Then vv performs a force in the efficient set Rev()\operatorname{Rev}(\mathcal{F}) of Term()\operatorname{Term}(\mathcal{F}). Since every force is performed at time 1, degvZ(G)\deg v\leq\operatorname{Z}(G). ∎

The converse of Proposition 3.27 is false, as the next example demonstrates.

Example 3.28.

Let GG be the graph in Figure 10.

Refer to caption
Figure 10: A graph GG with vBEff(G)Bv\in\bigcap_{B\in\operatorname{Eff}(G)}B and degv<Z(G)\deg v<\operatorname{Z}(G)

It can be verified that Z(G)=5\operatorname{Z}(G)=5. Then B1={a,b,c,d,v}B_{1}=\{a,b,c,d,v\} and B2={x,y,z,w,v}B_{2}=\{x,y,z,w,v\} are minimum zero forcing sets and pt(G,B1)=pt(G,B2)=1\operatorname{pt}(G,B_{1})=\operatorname{pt}(G,B_{2})=1, so pt(G)=1\operatorname{pt}(G)=1. Let BB be a zero forcing set of GG not containing vertex vv. In order to have pt(G,B)=1\operatorname{pt}(G,B)=1, some neighbor of vv must be able to force vv immediately. Without loss of generality, this neighbor is dd. Then a,b,c,d,wBa,b,c,d,w\in B. The set {a,b,c,d,w}\{a,b,c,d,w\} is a minimum zero forcing set but has propagation time 22, because no vertex can force zz immediately. Thus BEff(G)B={v}\bigcap_{B\in\operatorname{Eff}(G)}B=\{v\}, and observe that degv=4<5=Z(G)\deg v=4<5=\operatorname{Z}(G).

4 Relationship of propagation time and diameter

In general, the diameter and the propagation time of a graph are not comparable. Let GG be the dart (shown in Figure 2). Then diam(G)=2<pt(G)=3\operatorname{diam}(G)=2<\operatorname{pt}(G)=3. On the other hand, diam(C4)=2>1=pt(C4)\operatorname{diam}(C_{4})=2>1=\operatorname{pt}(C_{4}).

Although it is not possible to a obtain a direct ordering relationship between diameter and propagation time in an arbitrary graph, diameter serves as an upper bound for propagation time in the family of trees. To demonstrate this, we need some definitions. The walk v1v2vpv_{1}v_{2}\dots v_{p} in GG is the subgraph with vertex set {v1,v2,,vp}\{v_{1},v_{2},\dots,v_{p}\} and edge set {v1v2,v2v3,,vp1vp}\{v_{1}v_{2},v_{2}v_{3},\dots,v_{p-1}v_{p}\} (vertices and/or edges may be repeated in these lists but are not repeated in the vertex and edge sets). A trail is a walk with no repeated edges (vertices may be repeated; a path is a trail with no repeated vertices). The length of a trail PP, denoted by len(P)\operatorname{len}(P), is the number of edges in PP. We show in Lemma 4.1 below that for any graph GG and minimum zero forcing set BB, there is a trail of length at least pt(G,B)\operatorname{pt}(G,B). A trail produced by the method in the proof is illustrated in the Example 4.2 below.

Lemma 4.1.

Let GG be a graph and let BB be a minimum zero forcing set of GG. Then there exists a trail PP such that pt(G,B)len(P)\operatorname{pt}(G,B)\leq\operatorname{len}(P).

Proof.

Observe that if u,vV(G)u,v\in V(G) such that uu forces vv at time t>1t>1, then uu cannot force vv at time t1t-1. Thus either uu was forced at time t1t-1 or some neighbor of uu was forced at time t1t-1. So there is a path wuvwuv, where ww forces uu at time t1t-1, or a path wxuvwxuv, where ww forces xx at time t1t-1 and xx is a neighbor of uu.

We construct a trail vpvp+1vp+2v0v_{-p}v_{-p+1}v_{-p+2}\dots v_{0}, such that for each time tt, 1tpt(G,B)1\leq t\leq\operatorname{pt}(G,B), there exists an iti_{t}, pit1-p\leq i_{t}\leq-1, such that vitv_{i_{t}} forces vit+1v_{i_{t}+1} at time tt. Begin with t=pt(G,B)t=\operatorname{pt}(G,B) and work backwards to t=1t=1 to construct the trail, using negative numbering. To start, there is some vertex v0v_{0} that is forced by a vertex v1v_{-1} at time t=pt(G,B)t=\operatorname{pt}(G,B); the trail is now v1v0v_{-1}v_{0}. Assume the trail vjv0v_{-j}\dots v_{0} has been constructed so that for each time t=,,pt(G,B)t=\ell,\dots,\operatorname{pt}(G,B), there exists an iti_{t}, such that vitv_{i_{t}} forces vit+1v_{i_{t}+1} at time tt. If >1\ell>1, then vjvj+1v_{-j}\to v_{-j+1} at t=t=\ell. Thus either vj1vjvj+1v_{-j-1}v_{-j}v_{-j+1} or vj2vj1vjvj+1v_{-j-2}v_{-j-1}v_{-j}v_{-j+1} is a path in GG, and we can extend our trail to vj1vjv0v_{-j-1}v_{-j}\dots v_{0} or vj2vj1vjv0v_{-j-2}v_{-j-1}v_{-j}\dots v_{0}. It should be obvious that no forcing edge will appear in this walk multiple times (by our construction). If uvuv is not a forcing edge, then it can only appear in our walk if uu^{\prime} forced uu and vv forced vv^{\prime}. Let uuvvu^{\prime}uvv^{\prime} be the first occurrence of uvuv in our walk. If uvuv were to occur again, either uu or vv would need to be forced at this later time, but this cannot happen because uu and vv are both black at this point. ∎

Example 4.2.

Let GG be the graph in Figure 11. As shown by the numbering in the figure, pt(G)=4\operatorname{pt}(G)=4, but GG does not contain a path of length 4. The trail produced by the method of proof used in Lemma 4.1 is abcdecfabcdecf and has length 6.

Refer to caption
Figure 11: A graph GG that does not have a path of length pt(G)\operatorname{pt}(G)
Theorem 4.3.

Let TT be a tree and BB be a minimum zero forcing set of TT. Then pt(T,B)diam(T)\operatorname{pt}(T,B)\leq\operatorname{diam}(T). Hence, pt(T)PT(T)diam(T)\operatorname{pt}(T)\leq\operatorname{PT}(T)\leq\operatorname{diam}(T).

Proof.

Choose BB to be a minimum zero forcing set such that pt(T,B)=PT(T)\operatorname{pt}(T,B)=\operatorname{PT}(T). By Lemma 4.1, there exists a trail in TT of length at least pt(T,B)\operatorname{pt}(T,B). Since between any two vertices in a tree there is a unique path, any trail is a path and the diameter of TT must be the length of the longest path in TT. Therefore,

pt(T)PT(T)=pt(T,B)diam(T).\operatorname{pt}(T)\leq\operatorname{PT}(T)=\operatorname{pt}(T,B)\leq\operatorname{diam}(T).\qed

The diameter of a graph GG can get arbitrarily larger than its minimum propagation time. The next example exhibits this result, but first we observe that if GG is a graph having exactly \ell leaves, then Z(G)2\operatorname{Z}(G)\geq\lceil\frac{\ell}{2}\rceil since at most two leaves can be on a maximal forcing chain.

Example 4.4.

To construct a kk-comb, we append a leaf to each vertex of a path on kk vertices, as shown in Figure 12 (our kk-comb is the special case Pk,2P_{k,2} of a more general type of a comb Pk,P_{k,\ell} defined in [8]). Let GG denote a kk-comb where k0mod4k\equiv 0\mod 4. It is clear that diam(G)=k+1\operatorname{diam}(G)=k+1. If we number the leaves in path order starting with one, then the set BB consisting of every leaf whose number is congruent to 2 or 3mod4\mod 4 (shown in black in the Figure 12) is a zero forcing set, and |B|=k2|B|=\frac{k}{2}. Since Z(G)k2\operatorname{Z}(G)\geq\frac{k}{2}, BB is a minimum zero forcing set. Then pt(G)pt(G,B)=3\operatorname{pt}(G)\leq\operatorname{pt}(G,B)=3. Since |G|=2k|G|=2k, Z(G)=k2\operatorname{Z}(G)=\frac{k}{2}, and pt(G)|G|Z(G)Z(G)\operatorname{pt}(G)\geq\frac{|G|-\operatorname{Z}(G)}{\operatorname{Z}(G)}, pt(G)3\operatorname{pt}(G)\geq 3. Therefore, pt(G)=3\operatorname{pt}(G)=3. Thus the diam(G)=k+1\operatorname{diam}(G)=k+1 is arbitrarily larger than pt(G)=3\operatorname{pt}(G)=3.

Refer to caption
Figure 12: A kk-comb

Acknowledgement The authors thank the referees for many helpful suggestions.

References

  • [1] AIM Minimum Rank – Special Graphs Work Group (F. Barioli, W. Barrett, S. Butler, S. M. Cioabă, D. Cvetković, S. M. Fallat, C. Godsil, W. Haemers, L. Hogben, R. Mikkelson, S. Narayan, O. Pryporova, I. Sciriha, W. So, D. Stevanović, H. van der Holst, K. Vander Meulen, A. Wangsness). Zero forcing sets and the minimum rank of graphs. Linear Algebra and its Applications, 428: 1628–1648, 2008.
  • [2] F. Barioli, W. Barrett, S. Fallat, H. T. Hall, L. Hogben, B. Shader, P. van den Driessche, H. van der Holst. Zero forcing parameters and minimum rank problems. Linear Algebra and its Applications 433: 401–411, 2010.
  • [3] D. Burgarth and V. Giovannetti. Full control by locally induced relaxation. Physical Review Letters PRL 99, 100–501, 2007.
  • [4] K. B. Chilakamarri, N. Dean, C. X. Kang, E. Yi. Iteration Index of a Zero Forcing Set in a Graph. Bulletin of the Institute of Combinatorics and its Applications, 64: 57–72, 2012.
  • [5] C. J. Edholm, L. Hogben, M. Huynh, J. Lagrange, D. D. Row. Vertex and edge spread of zero forcing number, maximum nullity, and minimum rank of a graph. Linear Algebra and its Applications, in press.
  • [6] L. Hogben. Minimum rank problems. Linear Algebra and its Applications, 432: 1961–1974, 2010.
  • [7] D. D. Row. A technique for computing the zero forcing number of a graph with a cut-vertex. Linear Algebra and its Applications, in press.
  • [8] S. Severini. Nondiscriminatory propagation on trees. Journal of Physics A, 41: 482002 (Fast Track Communication), 2008.