Optimal Pebbling Number of the Square Grid

A pebbling move on a graph removes two pebbles from a vertex and adds one pebble to an adjacent vertex. A vertex is reachable from a pebble distribution if it is possible to move a pebble to that vertex using pebbling moves. The optimal pebbling number p opt is the smallest number m needed to guarantee a pebble distribution of m pebbles from which any vertex is reachable. The optimal pebbling number of the square grid graph P n h P m was investigated in several papers (Bunde et al. in J Graph Theory 57(3):215–238, 2008; Xue and Yerger in Graphs Combin 32(3):1229–1247, 2016; Gy } ori et al. in Period Polytech Electr Eng Comput Sci 61(2):217–223 2017). In this paper, we present a new method using some recent ideas to give a lower bound on p opt . We apply this technique to prove that p opt ð P n h P m Þ (cid:2) 213 nm . Our method also gives a new proof for p opt ð P n Þ ¼ p opt ð C n Þ ¼ 2 n 3 (cid:2) (cid:3) .


Introduction
Graph pebbling is a game on graphs. It was suggested by Saks and Lagarias to solve a number theoretic problem, which was done by Chung [4]. The main framework is the following: a distribution of pebbles is placed on the vertices of a simple graph. A pebbling move removes two pebbles from a vertex and places one pebble on an adjacent vertex. The goal is to reach any specified vertex by a sequence of pebbling moves. This may be viewed as a transportation problem on a graph where the cost of a move is one pebble. We begin with some notation needed to state our results.
Let G be a simple graph. We denote the vertex and edge set of G by V(G) and E(G), respectively. A pebble distribution P is a function from V(G) to the nonnegative integers. We say that G has P(v) pebbles placed at the vertex v under the distribution P. We say that a vertex v is occupied if PðvÞ [ 0 and unoccupied otherwise. The size of a pebble distribution P, denoted P j j, is the total number of pebbles placed on the vertices of G.
Let u be a vertex with at least two pebbles under P, and let v be a neighbor of u. A pebbling move from u to v consists of removing two pebbles from u and adding one pebble to v. That is, a pebbling move yields a new pebble distribution P 0 with P 0 ðuÞ ¼ PðuÞ À 2 and P 0 ðvÞ ¼ PðvÞ þ 1. A sequence of pebbling moves is executable under a pebble distribution if we can apply its moves one after the another so that the number of pebbles is nonnegative after each move on any vertex. We call such a sequence a pebbling sequence. We say that a vertex v is k-reachable under the distribution P if we can obtain, after a pebbling sequence, a distribution with at least k pebbles on v. If k ¼ 1 we say simply that v is reachable under P. More generally, a set of vertices S is k-reachable under the distribution P if, after a pebbling sequence, we can obtain a distribution with at least a total of k pebbles on the vertices in S.
A pebble distribution P on G is solvable if all vertices of G are reachable under P. A pebble distribution on G is optimal if it is solvable and its size is minimal among all of the solvable distributions of G. Note that optimal distributions are usually not unique.
The optimal pebbling number of G, denoted by p opt ðGÞ, is the size of an optimal pebble distribution. In general, the decision problem for this graph parameter is NPcomplete [5].
We denote with P n and C n the path and cycle on n vertices, respectively. The Cartesian product GhH of graphs G and H is defined in the following way: VðGhHÞ ¼ VðGÞ Â VðHÞ and fðg 1 ; h 1 Þ; ðg 2 ; h 2 Þg 2 EðGhHÞ if and only if fg 1 ; g 2 g 2 EðGÞ and h 1 ¼ h 2 or fh 1 ; h 2 g 2 EðHÞ and g 1 ¼ g 2 .
Let u and v be vertices of graph G. The distance between v and u, namely the number of edges contained in the shortest path between u and v, is denoted by d(v, u). The distance k neighborhood of v contains the vertices whose distance from v is exactly k. We denote this set with N k ðvÞ.
The optimal pebbling number is known for several graphs including paths, cycles [1,6,7], caterpillars [8] and m-ary trees [9]. The optimal pebbling number of grids has also been investigated. Exact values were proved for P n hP 2 [1] and P n hP 3 [2].
The question for bigger grids is still open. The best known upper bound for the square grid can be found in [3]. Diagonal induced subgraphs of the square grid was studied in [10].
Instead of the square grid on the plane it is easier to work with the square grid on the torus. As the plane grid is a subgraph of this, any lower bound on the torus grid will also give a lower bound on the plane grid as well. It is well known that the torus grid is a vertex transitive graph, i.e. given any two vertices v 1 and v 2 of G, there is some automorphism f : VðGÞ ! VðGÞ such that f ðv 1 Þ ¼ v 2 . Some of our statements will be stated for all vertex transitive graphs.
In this paper we present a new method giving a lower bound on the optimal pebbling number. We obtain 2 13 VðGÞ as a lower bound for the optimal pebbling number of the square grid, which is better than the previously known bounds.Please provide the complete details for reference [11].
In Sect. 2 we show that the concept of excess-introduced in [2]-can be used to improve the fractional lower bound on the optimal pebbling number. The higher the total excess, the better the obtained bound on the optimal pebbling number is. The problem is that this method is not standalone, because excess can be zero and zero excess does not give us any improvement. Therefore the main objective of the rest of the paper is to give a lower bound on the excess using some other pebbling tools.
In Sect. 3 we study the concept of cooperation. Cooperation is the phenomenon which makes pebbling hard. We show there, that if cooperation can be bounded from above, then we can state a lower bound on the optimal pebbling number. We invent the tool called cooperation excess, which is a mixture of cooperation and excess. In this section we state and prove several small claims which will be required later to prove Lemma 4.1. This lemma is the essence of our work. It shows that if the total excess is small, then there is not much cooperation and if cooperation is huge, then the total excess is also large. Therefore in each case one of our two lower bounds works well.
Unfortunately, the proof of Lemma 4.1 is quite complicated. The third part of Sect. 3 and the whole Sect. 4 contain the parts of this proof. In Sect. 5 we show a general method which can be used to give a lower bound on the optimal pebbling number. This method relies on Lemma 4.1. Using this method we show that p opt ðP n hP m Þ ! 2 13 nm. We also present a new proof for

Improving the Fractional Lower Bound
The optimal pebbling number problem can be formulated as the following integer programming problem [11], where v 1 ; v 2 ; . . .; v n are the vertices of the given graph: Its fractional relaxation can be solved efficiently, and its solution is called the fractional optimal pebbling number, which gives a lower bound on the optimal pebbling number. Originally it was defined in a bit different way, but this is an equivalent definition. You can find the details of fractional pebbling in [11]. Notice that some vertices must be 2-reachable in a solvable distribution if there is an unoccupied vertex. Optimal distributions usually contain many unoccupied and several 2-reachable vertices. However, in some sense, 2, 3, or more reachability wastes the effect of pebbles. Also 3-reachability induces larger waste than 2-reachability. In order to measure this waste we use the notion called excess, which was introduced in [2].
Definition Let Reach ðP; vÞ be the greatest integer k such that v is k-reachable under P. The excess of v under P is Reach ðP; vÞ À 1 if v is reachable and zero otherwise. It is denoted by Exc ðP; vÞ.
We are interested in the total amount of waste, therefore we define the notation of total excess of P, which is TE ðPÞ ¼ P v2 V ðGÞ Exc ðP; vÞ.
Definition An effect of a pebble placed at v is the following: Herscovici et al. proved that the fractional optimal pebbling number of a vertextransitive graph is jVðGÞj= ef ðvÞ, therefore it is a lower bound on the optimal pebbling number. The corollary of the next theorem improves this bound. Naturally, this bound is useless without a proper estimate of total excess. To say something useful about it we look at the optimal pebbling problem from a different angle.

Cooperation Between Distributions
In this section we talk about cooperation, which makes pebbling hard.

Pebbling Cooperation
Definition Let P and Q be pebble distributions on graph G. Now P þ Q is the unique pebble distribution on G which satisfies ðP þ QÞðvÞ ¼ PðvÞ þ QðvÞ. P and Q are disjoint when no vertex has pebbles under both distributions.
Definition The coverage of a distribution P is the set of vertices which are reachable under P. We denote the size of this set with Cov ðPÞ.
A natural idea to find small solvable distributions is finding a distribution with small size and huge coverage and make it solvable by placing some more pebbles.
In the rest of the section we assume that we add disjoint distributions P and Q together. We would like to establish an upper bound using Cov ðPÞ þ Cov ðQÞ on Cov ðP þ QÞ . Similarly, we are interested in some relation between TE ðP þ QÞ and TE ðPÞ þ TE ðQÞ.
Definition A cooperation vertex is neither reachable under P nor Q, but it is reachable under P þ Q. We denote the number of such vertices with Coop ðP; QÞ. A double covered vertex is reachable under both P and Q, we denote the size of their set with DC ðP; QÞ.
The following claim is a trivial consequence of the definitions. Definition We say that a distribution U is a unit, if only one vertex has pebbles under U.
Units are the building blocks of pebble distributions in the following sense: Any distribution P can be written as P ujPðuÞ [ 0 P u , where P u is a unit having P(u) pebbles at u. Units have two main advantages over other distributions. Their coverage and total excess can be easily calculated: Claim 3.2 Let U be a unit distribution which places pebbles at vertex u. Then we have that

Combining Cooperation and Excess
We would like to distinguish the sources of excess. Does it come from P or Q or does it arise from the ''cooperation of P and Q''?
Definition The unit excess of P, denoted by UE ðPÞ, is P ujPðuÞ [ 0 ðTEðP u ÞÞ, where P u is a unit on u containing exactly P(u) pebbles and all of them are placed at u.

Definition
The cooperation excess of a vertex v is Exc ðP þ Q; vÞ À ð Exc ðP; vÞ þ Exc ðQ; vÞÞ. If it is positive, then we say that v has cooperation excess.
Similarly, the cooperation excess between P and Q is the total excess of P þ Q minus the total excesses of P and Q. Denote this with CE(P, Q).
We have mentioned previously, that we can split any pebbling distribution into disjoint unit distributions. If we get t unit distributions, then the application of Claim 3.1 and the definition of cooperation excess gives the following results. Claim 3.3 Let P be a pebble distribution on G and let D be a disjoint decomposition of P to unit distributions. Denote the elements of D with U 1 ; U 2 ; . . .; U t . Now Cov ðU i Þ can be calculated easily. The ''effect'' of cooperation is calculated in the other, more complicated terms. Lemma 4.1 is going to establish a connection between those quantities in a fruitful way.

Connection Between Cooperation and Excess
Now let us consider an arbitrary graph G, and let D be the maximum degree of G. In the rest of the section we assume that Q ¼ U is a unit having pebbles only at vertex u and its size is not zero. Now we state some basic claims about the recently defined objects.

Claim 3.4 Each cooperation vertex c has a neighbor that has cooperation excess.
Proof A cooperation vertex c is not reachable under P or U. Therefore none of its neighbors is 2-reachable under these distributions. On the other hand, c is reachable under P þ U, hence there is a neighbor n of c which is 2-reachable under this distribution. This means that n has cooperation excess.h Definition If a vertex is not a cooperation vertex and it does not have cooperation excess, then we call it cooperation free.
This name is a somewhat misleading, because these vertices can participate in cooperation in a sophisticated way. For an example see Fig. 1.
Definition Let r be a pebbling sequence. P r denotes the pebble distribution which is obtained by the application of r to distribution P. A vertex is utilized by a pebbling sequence if there is a move in the sequence which removes or adds a pebble to the vertex. Let M(v) be the minimal number of cooperation vertices, including v if it is a cooperation vertex, which are utilized by a pebbling sequence r which satisfies that ðP þ UÞ r ðvÞ ! 2. If v is not 2-reachable under P þ U, then we say that MðvÞ ¼ 1.
An example where M() is shown for each vertex can be seen in Fig. 2.

Claim 3.5
If there is an available pebbling move which removes a pebble from a cooperation vertex c, then either two neighbors of c, say e and f, have cooperation excess at least 1 with MðeÞ\MðcÞ and Mðf Þ\MðcÞ or a neighbor d has cooperation excess at least 3 with MðdÞ\MðcÞ.
Proof The condition implies that c can obtain two pebbles by some pebbling moves under P þ U. Consider a pebbling sequence r which does this by utilizing M(c) cooperation vertices. Either r moves the two pebbles to c from two different neighbors e and f, or it can move both pebbles from the same neighbor d. None of the neighbors are 2-reachable under P or U, but e, f and d has to be 2, 2 and 4 reachable under P þ U, respectively. This means that e and f have cooperation excess at least 1 and the cooperation excess of d is at least 3. Furthermore, r moves two pebbles to e and f or to d, then it moves them to c with some more moves. This shows that MðeÞ; Mðf Þ; MðdÞ\MðcÞ.h Claim 3.6 If the cooperation excess of a vertex v is at least 3 and one of its neighbors, say c, is a cooperation vertex, then there is a vertex w that is adjacent to v and MðwÞ MðvÞ. [0] [0] [0] [0] [0] [0] u Proof Vertex v gets a pebble under P þ U, so a neighbor n is 2-reachable under P þ U. If n is not 2-reachable under P or U, then it has cooperation excess.h Remark In fact, a stronger property holds. If a vertex v gains an extra pebble by cooperation, then it can happen in two ways: A neighbor gained extra pebbles and it passes one of them. Or there are two or more neighbors of v such that each of them can give some pebbles to v, but these moves somehow blocks each other. The advantage of the cooperation is that some previously blocked moves can be done simultaneously. This is the way how cooperation free vertices can ''help cooperation''.

Trajectories
Here we introduce a visualization of pebbling sequences, which is slightly different from the signature digraph used in several pebbling papers (i.e. in [5]).
Definition The trajectory of a pebbling sequence r, denoted by TðrÞ, is a digraph on the vertices of G without parallel edges, where (u, v) is a directed edge if and only if a pebbling move u ! v is contained in the sequence.
Definition The size of a pebbling sequence is the total number of moves contained in it. We say that r is a minimal pebbling sequence with property p if its size is minimal among all pebbling sequences having property p.
In the next proof we need a lemma which is frequently used to solve pebbling problems. It is called No-Cycle Lemma and proved in several papers [5,12,13]. We state this lemma in the language of this paper. Lemma 3.8 (No-Cycle [13]) Let P be a pebble distribution on graph G, and r be a pebbling sequence. There is a subsequence d whose trajectory does not contain directed cycles and P r ðvÞ P d ðvÞ for each vertex v.
This implies the following corollary: Corollary 3.9 If r is a minimal pebbling sequence which moves m pebbles to a vertex v, then its trajectory is acyclic. Proof The No-Cycle lemma yields that we can move the maximum possible number of pebbles to u without removing a pebble from u. We can move Reach ðU; uÞ þ 1 pebbles to u, which means that we move here a pebble of P while we keep the pebbles of U, so u is double covered.h The following definition will be crucial in the proof.
Definition We say that a path is a coopexcess path, if each inner vertex of the path has cooperation excess.

Lemma 3.11
Let v be a vertex which is not double covered but it has cooperation excess. There is a coopexcess path between v and a double covered vertex or there are at least two cooperation free vertices such that each of them is connected to v by a coopexcess path. If v is not 2-reachable under both P and U, then these paths does not contain a vertex whose M value is higher than M(v).
An example for the first case is shown in Fig. 2 where y is a double covered vertex and v, x, y is a coopexcess path. The second case can be seen in Fig. 1, where v, x and v, y are coopexcess paths connecting v to cooperation free vertices.
Proof Consider a pebbling sequence r moving Reach ðU; vÞ þ Reach ðP; vÞ þ 1 pebbles to v utilizing M(v) cooperation vertices. Consider some path in the trajectory of r connecting u to v. We can assume that the only sink in the trajectory of r is v. A cooperation vertex without cooperation excess can not be the tail of an arc which is contained in the trajectory, therefore each vertex in the trajectory is either cooperation free or it has cooperation excess.
If there is a path between u and v which is contained in the trajectory such that all vertices of this path have cooperation excess, then according to Claim 3.10 u is double covered and this path is a coopexcess path. If an u, v path which is included in the trajectory contains a vertex d which is double covered and each vertex between d and v has cooperation excess, then it is a coopexcess path which we are looking for. Otherwise, all of the u, v paths which are contained in the trajectory contain cooperation free vertices.
In each such path let w i denote the cooperation free vertex which is the closest vertex to v. If w i 6 ¼ w j exist, then we have found 2 cooperation free vertices such that each of them is connected to v by a coopexcess path.
In the remaining case there is only one such w. Either it is a cut vertex in the trajectory or w ¼ u. Let T be the set of vertices which are included in the trajectory. We divide T to three sets U, V and W in the following way: We remove w from the trajectory obtaining some components, then we place a vertex t of T to U if t is in the component containing u, similarly we place t to V if it is in the component containing v and place the remaining vertices to W. Now we add w to all of these sets. Let r u be the sequence of pebbling moves containing all moves of r which acts only on the vertices of U. We define r w and r v similarly. The sources of the latter two sequences are only w and vertices having pebbles under P.
If there is a cooperation free vertex in V which is not w, then the closest one to v is connected to v by a coopexcess path. Hence, assume that all vertices in the trajectory of d v have cooperation excess.
If w is reachable under U, then r w is empty (w is not double covered) and ðP þ UÞ r u ðwÞ Reach ðU; wÞ. Since w is cooperation free, we can replace r u with a pebbling sequence d which does not use any pebbles of P and ðP þ UÞ r u ðwÞ ¼ ðP þ UÞ d ðwÞ. Therefore dr v is a pebbling sequence under P þ U and ðP þ UÞ r ðvÞ ¼ ðP þ UÞ dr v ðvÞ ¼ Reach ðP; vÞ þ Reach ðU; vÞ þ 1. r v must use a pebble of P to do this, otherwise dr v is executable under U which is a contradiction. The trajectory of r v is connected, therefore there is a vertex which is double covered, furthermore each vertex in this trajectory is connected by a coopexcess path to v, so we are done.
If w is not reachable under U, then ðP þ UÞ r u r w ðwÞ Reach ðP; wÞ. Thus, there is a minimal pebbling sequence d which is executable under P and P d ðwÞ ¼ Reach ðP þ U; wÞ ¼ Reach ðP; wÞ. Clearly dr v is not executable under distribution P or ðP þ UÞ dr v ðvÞ\ðP þ UÞ r ðvÞ. Both cases require that d removes a pebble from a vertex contained in V.
Let X V be the set of vertices from which d removes a pebble. d is executable under P so these vertices are 2-reachable under P. Consider the trajectory of d. If any vertex x from X is connected in the trajectory with a vertex y contained in U without pass-through w, then each vertex in such a connecting path is 2-reachable under P, therefore it is cooperation free or has cooperation excess. So there is either an other cooperation free vertex connected by a coopexcess path to v, or there is a coopexcess path between v and y which is connected to u by a path in the trajectory of r which does not contain w, so that path has to contain a double covered or a cooperation free vertex, which is not w.
The remaining case is when w separates all elements of X from U in the trajectory of d.
Let d uw be a maximal subset of d which is executable without using the pebbles placed at X , and let d v be the remaining subsequence. d uw r v is not executable under P þ U or ðP þ UÞ d uw r v ðvÞ\ðP þ UÞ r ðvÞ ¼ ðP þ UÞ r u r w r v ðvÞ. Therefore r u r w moves more pebbles to w than d uw , but d v is executable under P r u r w , thus P r u r w d v ðwÞ [ P d ðwÞ, therefore w has cooperation excess.
To prove the second claim, consider the paths we have found. If they were part of the trajectory of r, then all of them are 2-reachable under P þ U, so their M value can not be higher than M(v). Otherwise, the path consists of vertices from the trajectory of r and some others whose M value is zero, since they are 2-reachable under P. h The following claim is a trivial consequence of the definitions. Claim 3.12 If u contains at least two pebbles and it is double covered, then one of its neighbors is also double covered.
In the rest of the section we assume, that u contains at least two pebbles, i. e. jUj ! 2. Therefore we can use the previous claim. Proof The previous claim handles the case when d is u, since the neighbor is connected to d ¼ u. So assume that d 6 ¼ u.
Since d is double covered, it is reachable from U, so it is connected to u by a path, whose vertices are 2-reachable under U. Therefore these vertices can not be cooperation vertices. If there is a vertex on this path which does not have cooperation excess, then the vertex closest to d satisfies the conditions of the second type. Otherwise, u has cooperation excess which means that it is double covered.h We are getting closer to establish a connection between the number of cooperation vertices and cooperation excess.
Definition We call a subset Q of V(G) a C-block, if (1) each pair of vertices in Q is connected by a coopexcess path, (2) it contains a vertex having cooperation excess and it is maximal with these properties.
Notice that the intersection of two C-blocks cannot contain a vertex having cooperation excess.

Lemma 3.14 Each C-block either
(3) contains at least two double covered vertices, or (4) contains one double covered vertex and one cooperation free vertex, or (5) contains at least two cooperation free vertices.
Proof Consider an arbitrary element v of Q which having cooperation excess. If the C-block does not have a double covered vertex, then Lemma 3.11 guarantees that two cooperation free vertices are connected to v by a coopexcess path, which means that they are contained in Q, so (5) is satisfied.
Otherwise Q contains a double covered vertex. According to Lemma 3.13, either there is an other double covered vertex in Q, or a cooperation free vertex. Thus either (3)

or (4) is satisfied.h
Later we generalize the notion of C-blocks, so that we keep the properties of 3.14. The following statement will be useful for this.

Lemma 3.15
If a vertex v having cooperation excess is adjacent to a cooperation vertex c such that MðvÞ\MðcÞ , then there are vertices e and f, such that each of them is either double covered or cooperation free and they are connected to v by coopexcess paths containing only vertices whose M values are smaller than M(c).
Proof Vertex v has a cooperation vertex neighbor, therefore v is not 2-reachable under P or U. According to Lemma 3.11 there is a double covered vertex or there are two cooperation free vertices who are connected to v by a coopexcess path containing only vertices whose M values are at most MðvÞ\MðcÞ. In the latter case we are done. Since the double covered vertex is connected to an other double covered or cooperation free by a coopexcess path containing vertices whose M value is zero, according to Lemma 3.13. The concatenation of these two coopexcess paths fulfills the criteria.h

Connection Between Total Cooperation Excess, Number of Cooperation Vertices and Maximum Degree
In this section we prove a crucial lemma. Unfortunately, the proof requires quite a lot of effort, including many small claims. This lemma gives a connection between the total cooperation, the total number of double covered vertices and total cooperation excess. The proof would be relatively easy if the effect of a pebble would appear close to the location of the pebble. The example on Fig. 3. shows, that unfortunately this is not always true.
Another difficulty arises from the fact that a cooperation vertex can have cooperation excess. For such an example see Fig. 4. To prove Lemma 4.1 we get rid of such vertices one by one.
So to prove the lemma we will change the graph in several steps. In the new graph it will be easier to isolate these effects.
We introduce a sequence of auxiliary graphs A 0 ; A 1 ; . . .; A k , whose vertices are labeled with vectors of four coordinates. The first and fourth coordinate is always an integer, while the other coordinates are binary. We denote the vertices of these graphs with underlined letters and the ith coordinate of vertex b with b i . We encode the parameters of the investigated pebbling problem in the auxiliary graph and in the coordinates in the following way: A 0 is isomorphic to G. The first coordinate of each vector is the amount of the cooperation excess of the corresponding vertex. The second coordinate is 1 iff the corresponding vertex is a cooperation vertex. The third coordinate is 1 when the vertex is double covered. Finally, the last coordinate is M(v), i.e. the minimum number of cooperation vertices have to be utilized by a pebbling sequence to obtain 2 pebbles at v, where v is the corresponding vertex. So A 0 is representation of the original configuration, the labels give the values of the various quantities that we are interested in. An example can be seen on Fig. 5.
The other graphs in the sequence A 1 ; . . .; A k will be obtained from A 0 by applying certain operations recursively, until we finally obtain A k with some useful properties. It is important to note that although the labels of A 0 are obtained from the pebble distribution on G, this will not be true any more for the other auxiliary graphs. We are not trying to change the graph and the pebble distribution and then obtain the new labels from these. We just apply the transformation on the abstract, labeled graphs. Now we translate the properties of the pebble distribution to properties of A 0 . Fig. 3 The triangles are cooperation vertices. Notice that they can be far from the added unit if a; b 2 B, then there is an A-path which connects them and B is maximal to these properties.
Note that the concept of A-path and A-block are generalizations of coopexcess path and C-block, respectively. In this language, the statement of Lemma 4.1 can be formulated as: X a2A a 2 À X a2A a 3 X a2A a 1 ðD À 2Þ We state four properties of A 0 which will be inherited to later auxiliary graphs. The significant properties are the first and the last. The other two are technical ones which will help the proof of inheritance stated in Claim 4.5.   Fig. 5 An example graph G with a pebble distribution P, a unit U and the corresponding auxiliary graph A 0 . A 1 is obtained from A 0 by using the first transformation. Note that A 1 does not contain a saturated vertex This claim is equivalent to the following, previously proven, statements with the new notation: Claim 3.5 ! (8), Claim 3.6 ! (9), Lemma 3.15 ! (10) and Lemma 3.14 ! (11).
We will obtain A i from A iÀ1 by applying one of two transformations. Then we repeat this until it is possible to apply at least one of the transformations. Both transformations will preserve P a2A a i , ði 2 f1; 2; 3gÞ, the fourth coordinate of each vertex and D, the maximum degree in the graph. The objective of the transformations is to replace vertices satisfying a 1 a 2 [ 0 (i.e. it has cooperation excess and it is a cooperation vertex) with (one ore more) vertices satisfying b 1 b 2 ¼ 0. From this point, we call these vertices saturated vertices. Both transformations will increase the number of vertices in the auxiliary graph.
Let w be a vertex where w 1 w 2 [ 0 such that its fourth coordinate is maximal among these vertices. By Claim 4.2 (8) there are two cases.
Case 1: If w has a neighbor x such that x 1 ! 3 and w 4 [ x 4 , then we apply the following transformation to A i : Transformation 1 -Choose a neighbor y of x such that its fourth coordinate is minimal among all neighbors of x. -Let R be the set of x's neighbors without y where the product of the first and the second coordinate is positive. -Delete x and add three vertices x 1 , x 2 and x 3 , such that x 1 1 ¼ x 3 1 ¼ 1 and . Connect x 2 with y, x 1 and x 3 . -Delete each element r of R and add two vertices r 1 and r 2 and set the coordinates as: r 1 1 ¼ r 1 3 ¼ 0, r 1 2 ¼ 1, r 1 4 ¼ r 2 4 ¼ r 4 , r 2 1 ¼ r 1 , r 2 2 ¼ 0 and r 2 3 ¼ r 3 . We connect r 1 to x 1 and r 2 to x 3 and to each original neighbor of r.
-We connect the neighbors of x which are not included in R [ y to x 3 . In other words, this transformation replaces each saturated neighbor of x (excluding a chosen y) with two vertices such that one of them is a leaf with zero first coordinate and the other one is act as the original vertex, but its second coordinate is zero. To handle the increased degree of x, we triple it. Also, if x is saturated then we add the additional x 4 vertex. Note that this can be done when D ! 4. If D 3, then we have to handle this case in a slightly different way.
Case 2: If w has two neighbors such that their first coordinates are positive and their fourth coordinates are strictly less than w 4 , then we apply the second transformation: Transformation 2 -We choose neighbors x and y whose fourth coordinate is minimal among all neighbors and x 4 ! y 4 . -We delete w and add vertices w 1 and w 2 . We set the coordinates of these vectors as: -We connect w 1 only with x. In contrast, we connect w 2 with all neighbors of w except y.
Both transformation can be seen on Fig. 6.   Proof We say that a vertex v is created by the ith transformation if v ¼ z j is a vertex of A i and z is a vertex of A iÀ1 . In this situation we say that v is a descendant of z. A vertex is involved in a transformation if either it is created by that or its vector is x y x y w r 1 w 2 r 2 r Fig. 6 Vertices denoted by squares are ''cooperation vertices'' so their second coordinates are one. Edges are shown between vertices contained in the same C-block. The upper transformation is called first, and the lower one is mentioned as the second transformation. Note that in the upper example R ¼ fw; rg changed by it.
Notice that the transformations keep the fourth coordinates of the vertices and if two vertices are descendants of the same vertex, then their fourth coordinates are the same. (8): If c is a saturated vertex in A i , then it is not created by the ith transformation and it was saturated in A iÀ1 also. If none of its neighbors were involved in the last transformation, then the property is clearly holds. Therefore assume the opposite. Assume that c had a neighbor d in A iÀ1 , such that d 1 ! 3 and d 4 \c 4 in A iÀ1 .
If d is contained in A i also, then the ith transformation did not change d 1 . In that case d and c are adjacent in A i and we are done.
Otherwise, the ith transformation created some descendants of d.
If it was Transformation 1 then a descendant of d is connected to c and either its first coordinate equals d 1 or it is d 1 À 2. In the first case we are done and the latter can happen if and only if d acted as x in that transformation. However, this is not possible, because this would mean that c acted as y but y 4 x 4 by (9) and the choice of y in Transformation 1, therefore c 4 ¼ y 4 x 4 ¼ d 4 \c 4 which is a contradiction.
The remaining case is that two descendants of d are created by Transformation 2. Since d 4 \c 4 , vertex c can not be x or y in Transformation 2, therefore it is adjacent to d 2 ¼ w 2 in A i and the first coordinate of this vertex equals d 1 .
Now assume that there are neighbors e and f in A iÀ1 such that their fourth coordinates are smaller than c 4 and e 1 ; f 1 [ 0. We may assume that e 1 ; f 1 \3, otherwise we obtain the previous case. Therefore neither e 1 nor f 1 can act as x in Transformation 1.
If e is contained in A i , then it is still adjacent to c. If e is replaced with some descendants by the ith transformation, then one of its descendants keep its first coordinate and that one is connected to c. Like in the previous case it cannot happen that c ¼ y and e ¼ x in Transformation 1 or c ¼ x and e ¼ w in Transformation 2. We can state the same for f . (9): If a is contained in A iÀ1 then a 1 ! 3 in A iÀ1 also. Therefore according to (9) there is a b which is adjacent to a and b 4 a 4 . Either b 4 or one of its descendants is adjacent to a in A i , therefore we are done.
Otherwise, a is a descendant of a vertex v. v has a neighbor b whose fourth coordinate is at most a 4 . There are several cases: First case: v ¼ r in Transformation 1, where r 2 R. If we remove x from the neighborhood of v and add x 3 we obtains the neighborhood of a. Therefore a has a neighbor whose fourth coordinate is not bigger.
Second case: v ¼ x in Transformation 1. y had the smallest fourth coordinate among the neighbors of x, thus y 4 x 2 4 ¼ a 4 and y and a are adjacent in A i . Third case: v ¼ w in Transformation 2. Since a 1 ! 3, a ¼ w 2 and y 4 w 4 ¼ w 2 4 according to (9).
(10): Transformation 1 keeps the A-paths, because it keeps connectivity and the first coordinate becomes zero only at leaves. Transformation 2 destroys some Apaths but all of them contain the saturated vertex which was handled by the transformation and whose fourth coordinate was at least c 4 .
(11): Transformation 1 does not split an A-block, furthermore it keeps the number of vertices whose third coordinate is one and whose first and second coordinate are both zero in each A-block.
Transformation 2 either does not split an A-block and keeps the investigated quantities, or it splits an A-block to two A-blocks. But (10) guarantees that both blocks contain enough vertices whose third coordinate is one or both first and second coordinates are zero.h   Definition Let B be an A-block in an auxiliary graph. We say that a vertex of B is inner vertex if its first coordinate is positive otherwise, it is called a boundary vertex.

Claim 4.9
Consider an A -block of an auxiliary graph A. Let the number of the boundary and the number of inner vertices denoted by band i, respectively. If a 1 a 2 ¼ 0 holds for each a 2 Athen b ðD À 2Þi þ 2 is satisfied.
Proof Proof by induction: The base case is an A-block with one inner vertex. This A-block is the closed neighborhood of the only inner vertex, therefore the number of boundary vertices is at most D. Now we assume that for any i\k the inequality is true. Let i ¼ k. We take a spanning tree of the inner vertices and consider a leaf vertex l. If we set l 1 to zero, then l becomes a boundary vertex and at most D À 1 boundary vertices, which are neighbors of l, are dropped from the A-block. The number of inner vertices is decreased by one, and the number of boundary vertices is decreased by at most D À 2. Using the induction hypothesis the proof is completed.h Proof of Lemma 4.8: Consider an A-block B and a boundary vertex v of B. Either NðBÞ where NðBÞ is the number of vertices in B whose first two coordinates are zero. It is also clear that P a2B a 1 ! i. Combining these observations and the previous claim: Proof We use Claim 3.3 (3.2) in the following inquality to obtain the second line. To obtain the third line, notice that if Otherwise jU i j ! 2 and we can apply Lemma 4.1.
This result together with the corollary of Theorem 2.1 and Claim 3.3 implies the following: This is a tool that helps to prove lower bounds on optimal pebbling number. Also notice that each element of the formula can be calculated efficiently.

Back to Square Grids
We would like to investigate finite square grids. It is easier to investigate torus graphs instead of square grids, because they are vertex-transitive. Let T m;n be the torus graph which we obtain if we glue together the opposite boundaries of P mþ1 hP nþ1 .
Note that T m;n ffi C m hC n . P m hP n can be obtained from T m;n by deleting some edges. Edge removal can not decrease the optimal pebbling number, therefore p opt ðT m;n Þ p opt ðP m hP n Þ. Therefore we work with T m;n in the rest of the section.
The size of the distance i neighborhood in T m;n is at most 4i. Thus Claim 3.2 gives the following estimates on excess and coverage of any unit placed on T m;n . To obtain these bounds it is enough to check small units and notice that the distance 2 neighborhood of u contains at least 8 5 U j j excess when U j j[ 4.
Claim 5.5 Let v be a vertex of T m;n . Then ef ðvÞ\9.
A similar result is proven in [2] for the square grid. We mimic that calculation.
Now we can obtain our new lower bound on the optimal pebbling number of the square grid: Theorem 5.6 The optimal pebbling number of T m;n is at least 2 13 nm, when m; n ! 5. Proof Let P be an optimal distribution of T m;n and let D be a disjoint decomposition of P to unit distributions. Denote the elements of D with U 1 ; U 2 . . .U t , such that jU i j jU iþ1 j. Let D ! 4 be the subset of D which contains all units whose size is at least four. Furthermore let D 2;3 be the set which contains the units whose size is two or three, and D 1 be the set of units whose size is one. Denote the total number of pebbles which are placed on vertices belonging to D 1 by S 1 . Define S 2;3 and S ! 4 similarly. It is clear that S 1 ¼ jPj À S 2;3 À S ! 4 .
We start with Corollary 5.2 and use the estimates of Claims 5.3, 5.4 and 5.5.
v Fig. 8 A mapping shows that there are at most 4 more vertices in the distance 3 neighborhood than in the distance 2 neighborhood 3k þ r ¼ Cov ðPÞ Cov ðU i Þ 3 2 jPj: So 2k þ 2r 3 jPj. |P| is integer, therefore this is equivalent to 2k þ r jPj.h

Open questions
Question 1 Is there a constant k, which does not depend on n, such that in an optimal distribution of P n hP n no pebble can be moved to the distance-k neighborhood of its initial location?
Question 2 If the answer for Question 1 is yes, then how small can be k?
We think that the answer for Question 1 is yes and we conjecture that k can be 4. Any finite k would improve our lower bound on p opt ðP n hP m Þ.