in the classical graph coloring problem we assign a color to each vertex of a given graph such that the graph is properly colored , i.e. , no two adjacent vertices have the same color . the chromatic number @xmath2 of a graph @xmath0 is the smallest @xmath3 such that @xmath0 can be colored with @xmath3 distinct colors .
deciding whether the chromatic number of a graph is at most @xmath3 is well known to be np - complete , even in the case with three colors .
the online variant of graph coloring can be defined as follows : the vertices of @xmath0 arrive one by one , and an online algorithm must color vertices as they arrive so that the revealed graph is properly colored at all times .
when a vertex arrives , the algorithm sees edges to previously colored vertices . to give the online algorithm slightly more power ,
a copy of @xmath0 is sent at the start of the input . however , the exact correspondence between the incoming vertices and the vertices of the copy of @xmath0 is not known to the algorithm .
this problem is called online graph coloring .    in this paper
we focus on a graph parameter called the _ online chromatic number _
@xmath4 of a graph @xmath0 .
this parameter is analogous to the standard chromatic number of a graph : it denotes the smallest number @xmath3 such that there exists a deterministic online algorithm which is able to color the specified graph @xmath0 using @xmath3 colors .
the online chromatic number has been studied since 1990 @xcite .
one of the main open problems in the area is the computational complexity of deciding whether @xmath5 for a specified simple graph @xmath0 , given @xmath0 and @xmath3 on input ; see e.g. kudahl @xcite .
we denote this decision problem as online chromatic number . in this paper
, we fully resolve this problem :    [ thm : main ] the decision problem online chromatic numberis pspace - complete .    as is usual in the online computation model
, we can view online graph coloringas a game between two players , which we call (representing the online algorithm ) and ( often called adversary in the online algorithm literature ) . in each round
chooses an uncolored vertex @xmath1 from @xmath0 and sends it to without telling him to which vertex of @xmath0 it corresponds , only revealing the edges to the previously sent vertices
. then must properly color ( `` paint '' ) @xmath1 , i.e. , cannot use a color of a neighbor of @xmath1 .
we stress that in this paper is restricted to be deterministic .
the game continues with the next round until all vertices of @xmath0 are colored .
deciding the outcome of many two - player games is pspace - complete ; among those are amazons , checkers or hex , to name a few . however , in most of these games both players have roughly the same power .
this does not hold for online graph coloringwhich is highly asymmetric , since has perfect information ( knows which vertices are sent and how they are colored ) , but does not .
may only guess to which part of the graph does the colored subgraph really belong .
this is the main difficulty in proving pspace - hardness .
* examples . *
consider a path @xmath6 on four vertices .
initially , sends two nonadjacent vertices .
if assigns different colors to them , then these are the first and the third vertex of @xmath6 , thus the second vertex must get a third color ; otherwise they obtained the same color @xmath7 and they are the endpoints of @xmath6 , therefore the second and the third vertex get different colors which are not equal to @xmath7 . in both cases , there are three colors on @xmath6 and thus @xmath8 , while @xmath9 .
note also that we may think of deciding where an incoming vertex belongs at some time _
after _ it is colored provided that his choice still allows for at least one isomorphism to the original @xmath0 .
this is possible only for a deterministic .    a particularly interesting class of graphs in terms of @xmath10 is the class of binomial trees .
a binomial tree of order @xmath3 is defined inductively : the binomial tree of order 0 is a single vertex ( the root ) and the binomial tree of order @xmath3 is created by taking two disjoint copies of binomial trees of order @xmath11 , adding an edge between their roots and choosing one of their roots as the root for the resulting tree .
thus @xmath6 is a binomial tree of order 2 with root on the second vertex of @xmath6 .
it is not hard to generalize the example of @xmath6 and show that the online chromatic number of the binomial tree of order @xmath3 is @xmath12  @xcite .
this shows that the ratio between @xmath10 and @xmath2 may be arbitrarily large even for the class of trees .    * history and related work .
* the online problem online graph coloringhas been known since 1976 @xcite , originally studied in the variant where the algorithm has no extra information at the start of the input .
bean  @xcite showed that no online algorithm that is compared to an offline algorithm can perform well under this metric .
the notion of online chromatic number was first defined in 1990 by @xcite .    for the online problem , lovsz ,
saks and trotter @xcite show an algorithm with a _ competitive ratio _
@xmath13 , where the competitive ratio is a ratio of the number of colors used by the online algorithm to the ( standard ) chromatic number .
this was later improved to @xmath14 by kierstad @xcite using a deterministic algorithm .
there is a better @xmath15-competitive randomized algorithm against an oblivious adversary by halldrsson  @xcite .
a lower bound on the competitive ratio of @xmath16 was shown by halldrsson and szegedy @xcite .
our variant of online graph coloring , where the algorithm receives a copy of the graph at the start , was suggested by halldrsson @xcite , where it is shown that the lower bound @xmath16 also holds in this model .
( note that the previously mentioned algorithmic results are valid for this model also . )
kudahl  @xcite recently studied online chromatic numberas a complexity problem .
the paper shows that the problem is conp - hard and lies in pspace .
later  @xcite he proved that if some part of the graph is precolored , i.e. , some vertices are assigned some colors prior to the coloring game between and and also reveals edges to the precolored vertices for each incoming vertex , then deciding whether @xmath17 is pspace - complete .
we call this decision problem online chromatic number with precoloring .
the paper  @xcite conjectures that online chromatic number(with no precolored part ) is pspace - complete too .
interestingly , it is possible to decide @xmath18 in polynomial time  @xcite .
keep in mind that while online graph coloringis an online problem , online chromatic numberis an ( offline ) decision problem of checking whether @xmath19
.    * proof outline .
* the fact that online chromatic numberbelongs to pspace is not hard to see : the online coloring is represented by a game tree which is evaluated using the minimax algorithm .
this can be done in polynomial space , since the number of rounds in the game is bounded by @xmath20 , i.e. , the number of vertices , and possible moves of each player can be enumerated in polynomial space : has at most @xmath20 possible moves , because it either uses a color already used for a vertex , or it chooses a new color , and has @xmath21 moves where @xmath22 is the number of colored vertices , since it chooses which colored vertices shall be adjacent to the incoming vertex . must ensure that sent vertices form an induced subgraph of @xmath0 , but this can be checked in polynomial space .    inspired by @xcite , we prove the pspace - hardness of online chromatic numberby a reduction to q3dnf - sat , i.e. , the satisfiability of a fully quantified formula in the 3-disjunctive normal form ( 3-dnf ) .
an example of such a formula is    @xmath23    the similar problem of satisfiability of a fully quantified formula in the 3-conjunctive normal form is well known to be pspace - complete .
since pspace is closed under complement , q3dnf - sat is pspace - complete as well .
note that by an easy polynomial reduction , we can assume that each 3-dnf clause contains exactly three literals .
we show the hardness in several iterative steps .
first , in section  [ sec : largeprecol ] , we present a new , simplified proof of the pspace - hardness of online chromatic number with precoloringin which the sizes of both precolored and non - precolored parts of our construction are linear in the size of the formula .
then , in section  [ sec : smallprecol ] , we strengthen the result by reducing the size of the precolored part to be logarithmic in the size of the formula .
this is achieved by adding linearly many vertices to our construction .
finally , in section  [ sec : noprecol ] , we show how to remove one precolored vertex and replace it by a non - precolored part , while keeping the pspace - hardness proof valid .
the cost for removing one vertex is that the size of the graph is multiplied by a constant , but since we apply it only logarithmically many times , we obtain a graph of polynomial size and with no precolored vertex .
this will complete the proof of theorem  [ thm : main ] .
we remark that removing the last precolored vertex is the most difficult part of proving pspace - hardness of online chromatic number .
still , our technique for removing a precolored vertex can be used for any graph satisfying a few assumptions . in our analysis
, often uses the natural greedy algorithm , which is ubiquitous in the literature ( see @xcite , @xcite ) :    the online algorithm colors an incoming vertex @xmath24 using the smallest color not present among colored vertices adjacent to @xmath24 .
our first construction will reduce the pspace - complete problem q3dnf - sat to online coloring with precoloring with a large precolored part . given a fully quantified formula @xmath25 in the 3-disjunctive normal form
, we will create a graph @xmath26 that will simulate this formula .
we assume that the formula contains @xmath20 variables @xmath27 and @xmath28 clauses @xmath29 , @xmath30 , and that variables are indexed in the same order as they are quantified .
our main resource will be a large precolored clique @xmath31 on @xmath3 vertices and naturally using @xmath3 colors ; the number @xmath3 will be specified later .
using such a precolored clique , we can restrict the allowed colors on a given uncolored vertex @xmath1 by connecting it with the appropriate vertices in @xmath31 , i.e. , we connect @xmath1 to all vertices in @xmath31 which do not have a color allowed for @xmath1 .    for simplicity
we use the precoloring in the strong sense , i.e. , is able to recognize which vertex in @xmath31 is which .
we use this to easily recognize colors .
however , it is straightforward to avoid the strong precoloring by modifying the precolored part ; for example by creating @xmath32 independent and identical copies of the @xmath32-th vertex in @xmath31 , each having the same color and the same edges to other vertices in @xmath31 and the rest of the graph .
with such a modification , would able to recognize each color by the number of its vertices in @xmath31 .
each vertex in @xmath31 thus corresponds to a color .
colors used by are naturally denoted by numbers 1 , 2 , 3 , ...
@xmath3 , but we shall also assign meaningful names to them .
we want to construct a graph @xmath26 that has the online chromatic number @xmath3 if and only if the quantified 3-dnf formula can be satisfied .
see figure  [ fig : gadgetexample ] for an example of @xmath26 and an overview of our construction .
we use the following gadgets for variables and clauses :    for a variable @xmath33 which is quantified universally , we will create a gadget consisting of _ universal _ vertices @xmath34 and @xmath35 , connected by an edge .
the vertex @xmath34 represents the positive literal @xmath33 , while @xmath35 represents the negative literal @xmath36 .
both vertices have exactly two allowed colors : @xmath37 and @xmath38 .
if @xmath34 is assigned the color @xmath39 , it corresponds to setting the variable @xmath33 to @xmath40 , and vice versa .    note that if presents a vertex @xmath41 to
, is able to recognize that it is a vertex corresponding to the variable @xmath42 , but it is not able to recognize whether it is the vertex @xmath41 or @xmath43 .    for a variable @xmath42 which is quantified existentially , we will create a gadget consisting of three _ existential _
vertices @xmath41 ( for the positive literal @xmath42 ) , @xmath43 ( for the literal @xmath44 ) and @xmath45 ( the helper vertex ) , connected as a triangle .
coloring of the first two variables also corresponds to setting the variable @xmath42 to true or false , but in a different way : @xmath41 has allowed colors @xmath46 and @xmath47 , while @xmath43 has allowed colors @xmath48 and @xmath47 .
we want to avoid both @xmath41 and @xmath43 to have the color of type @xmath49 , and so the `` helper '' vertex @xmath45 can be colored only by @xmath46 or @xmath48 .
note that the color choice for the vertices of @xmath42 means that if is presented any vertex of this variable , can recognize it and decide whether to set @xmath42 to @xmath40 ( and color accordingly ) or to @xmath50 .
we call existential and universal vertices together _ variable vertices_.    for each clause @xmath29 , we will add four new vertices .
first , we create a vertex @xmath51 for each literal in the clause , which is connected to one of the vertices @xmath34 and @xmath35 corresponding to the sign of the literal .
for example if @xmath52 , then @xmath51 is connected to @xmath34 , @xmath53 is connected to @xmath43 and @xmath54 to @xmath55 .
the allowed colors on a vertex @xmath51 are @xmath56 .    finally , we add a fourth vertex @xmath57 connected to the three vertices @xmath58 .
this vertex can be colored only using the color @xmath59 or the color @xmath60 .
the color @xmath60 is used to signal that this particular clause is evaluated to be @xmath50 .
if the color @xmath59 is used for the vertex @xmath57 , this means that the clause is evaluated to @xmath40 , because @xmath59 is not present on any of @xmath58 , thus they have colors of type @xmath61 and their neighbors corresponding to literals have colors of type @xmath49 .
the last vertex we add to the construction will be @xmath62 , a final vertex .
the vertex @xmath62 is connected to all the vertices @xmath57 corresponding to the clauses .
the allowed colors of the vertex @xmath62 are @xmath63 .
this final vertex corresponds to the final evaluation of the formula .
if all clauses are evaluated to @xmath50 , the vertex @xmath62 has no available color left and must use a new color .
we have listed all the vertices and colors in our graph @xmath26 and the functioning of our gadgets , but we will need slightly more edges .
the reasoning for the edges is as follows : if presents any vertex of the type @xmath64 or @xmath62 before presenting the variable vertices , or in the case when the variable vertices are presented out of the quantifier order , we want to give an advantage to so it can finalize the coloring .    this will be achieved by allowing to treat all remaining universal vertices as existential vertices , i.e. , can recognize which of the two universal vertices @xmath65 corresponds to setting @xmath42 to @xmath40 .    to be precise , we add the following edges to @xmath26 :    every existential vertex @xmath66 is connected to all previous universal vertices @xmath34 , that is to all such @xmath34 for which @xmath67 .
every universal vertex @xmath65 is connected to all previous universal vertices @xmath34 such that @xmath67 .
every vertex of type @xmath51 is connected to all the universal vertices @xmath68 for @xmath69 . note that @xmath51 is connected either to @xmath34 , or to @xmath35 ; we do not add an edge to such vertices .
every vertex of type @xmath57 is connected to all the universal vertices @xmath34 for all @xmath32 .
the vertex @xmath62 is connected to all the universal vertices @xmath34 for all @xmath32 .
we call all non - precolored vertices the _ gadgets _ for variables and clauses
.        the full analysis can be found in appendix  [ app : section2 ] .
we now make a step to the general case without precoloring by reducing the size of the precolored part so that it has only logarithmic size .
our construction is based on the one with a large precolored part ; namely , all the vertices @xmath70 ( the gadgets for variables and clauses ) and the whole color clique @xmath31 will be connected the same way .
let @xmath26 denote the gadgets for variables and clauses and @xmath31 .
since @xmath31 is now not precolored and may send it after the gadgets , we help by a structure for recognizing vertices in @xmath26 or for saving colors .
we remark that there is also a simpler construction with a logarithmic number of precolored vertices .
if we just add precolored vertices to recognize vertices in @xmath26 , the following proof would work and be easier .
however , when we replace a precolored vertex @xmath1 by some non - precolored graph in section  [ sec : noprecol ] , we will use some conditions on the graph @xmath71 that this construction would not satisfy .
our structure will consist of many small _ nodes _ , all of them have the same internal structure , only their adjacencies with other vertices vary .
each node consists of three vertices and a single edge ; vertices are denoted by @xmath72 and the edge leads between @xmath73 and @xmath74 .
we call the vertices @xmath75 and @xmath73 the _ lower partite set _ of the node , @xmath74 form the _ upper partite set_. see figure  [ fig : node ] for an illustration of a node .
r2.8 cm     clearly , the online chromatic number of a node is two .
the intuition behind the nodes is as follows :    if presents vertices of a node in the correct way , needs to use two colors in the lower partite set of the node .
no color can be used in two different nodes .
each vertex @xmath76 ( in the gadgets and in @xmath31 ) has its own node @xmath77 that identifies @xmath1 .
if the vertex @xmath74 from @xmath77 does not arrive before @xmath1 is sent , can color @xmath74 and @xmath1 with the same color , thus save a color .
otherwise , can use the node to recognize @xmath1 .
universal vertices @xmath78 for each universally quantified variable @xmath33 should be distinguishable only by the same vertices as in the previous section .
therefore they are identified by the same two nodes .
let @xmath79 be the number of vertices in @xmath26 .
we create @xmath79 nodes , denoted by @xmath80 , one for each vertex in @xmath26 . for any two distinct nodes @xmath81 and @xmath82 ( @xmath83 ) ,
there is an edge between each vertex in @xmath81 and each vertex in @xmath82 .
therefore , no color can be used in two nodes .
we now define which nodes identifies which vertices and add edges between the original construction @xmath26 and nodes .
let @xmath84 be the vertices in @xmath26 ( in an arbitrary order ) .
then @xmath81 identifies the vertex @xmath85 .
moreover , if @xmath85 is a vertex @xmath55 or @xmath86 for a universally quantified variable @xmath87 and @xmath88 is the other vertex , then @xmath82 also identifies @xmath85 and @xmath81 also identifies @xmath88 .
thus each node identifies one or two vertices and each vertex is identified by one or two nodes .    for a vertex @xmath76 and for a node @xmath77 ,
if @xmath77 identifies @xmath1 , then we connect only the whole lower partite set of @xmath77 to @xmath1 , i.e. , we add two edges from @xmath1 to both @xmath75 and @xmath73 of @xmath77 .
otherwise , we add three edges  one between @xmath1 and every vertex in @xmath77 .
the only precolored part @xmath89 of the graph is intended for distinguishing nodes . since there are @xmath79 nodes in total , we have @xmath90 precolored vertices @xmath91 with no edges among them .
precolored vertices have a color that may be used later for coloring @xmath26 ( the gadgets and @xmath31 ) . for simplicity
, we again use the precoloring in the strong sense , i.e. , is able to recognize which precolored vertex is which .
we connect all vertices in the node @xmath81 to @xmath92 if the @xmath93-th bit in the binary notation of @xmath32 is 1 ; otherwise @xmath92 is not adjacent to any vertex in @xmath81 .
clearly , the node to which an incoming vertex belongs can be recognized by its adjacency to the precolored vertices .
note that a vertex from nodes is connected to at least one precolored vertex and there is no edge between @xmath26 and precolored vertices .
so far , we have introduced all vertices and edges in our construction of the graph @xmath71 .
the full analysis can be found in appendix  [ app : section3 ] .
in this section we show how to replace one precolored vertex by a large nonprecolored graph whose size is a constant factor of the size of the original graph , while keeping  s winning strategy in the case of a satisfiable formula .
s winning strategy in the other case is of course preserved also and easier to see .
we prove the following lemma which holds for all graphs with precolored vertices satisfying a few assumptions .
[ l : removeprecolvertex ] let @xmath0 be a graph with precolored subgraph @xmath94 created from a fully quantified formula @xmath95 , and let @xmath96 be its precolored vertex .
let @xmath97 be the induced subgraph with all nonprecolored vertices that are _ not _ connected to @xmath98 and let @xmath99 be the induced subgraph with all nonprecolored vertices that are connected to @xmath98 .
let @xmath3 be an integer .
assume that the following holds :    @xmath5 if and only if @xmath95 is satisfiable ,    in the winning strategy of in the case if @xmath95 is satisfiable , can color @xmath99 using before two nonadjacent vertices from @xmath97 arrive .
moreover , in this case if assigns the same color to a vertex in @xmath97 and to a vertex in @xmath99 before two nonadjacent vertices from @xmath97 arrive , can still color @xmath0 using @xmath3 colors .
then there exists an integer @xmath100 and a graph @xmath101 with the following properties :    @xmath101 has only @xmath102 precolored vertices , and @xmath103 ,    @xmath101 can be constructed from @xmath0 in polynomial time with respect to the size of @xmath0 ,    it holds that @xmath104 if and only if @xmath95 is satisfiable .
the proof of theorem  [ thm : main ] using lemma  [ l : removeprecolvertex ] is in appendix  [ app : proofoftheorem ] . *
construction of @xmath101 .
* let @xmath79 be the total number of vertices in @xmath97 and @xmath99 and let @xmath105 .
our graph @xmath101 consists of precolored part @xmath106 , graphs @xmath97 and @xmath99 and three huge cliques @xmath107 and @xmath108 of size @xmath109 ; cliques @xmath107 and @xmath108 together form a _ supernode_. we keep the edges inside and between @xmath97 and @xmath99 and the edges between @xmath110 and @xmath111 as they are in @xmath0 .
we add a complete bipartite graph between cliques @xmath112 and @xmath108 , i.e. , @xmath113 forms a clique of size @xmath114 .
no vertex in @xmath77 is connected to @xmath112 or @xmath108 . in other words ,
the supernode is created from a node by replacing each vertex by a clique of size @xmath109 and the only edge in the node by a complete bipartite graph .
there are no edges between the supernode ( cliques @xmath77 and @xmath113 ) and a precolored vertex in @xmath110 .
it remains to add edges between the supernode and @xmath115 .
there is an edge between each vertex in @xmath99 and each vertex in the supernode , while every vertex in @xmath97 is connected only to the whole @xmath77 and @xmath112 , but not to any vertex in @xmath108 .
the fact that @xmath97 and @xmath108 are not adjacent at all is essential in our analysis .
our construction is depicted in figure  [ fig : supernode ] .     where parts @xmath107 and @xmath108 form a supernode and parts @xmath97 and @xmath99 form the original graph @xmath0 ( the remaining precolored vertices are not shown )
. a thick line connecting two parts of @xmath101 corresponds to a complete bipartite graph between them .
the edges between @xmath97 and @xmath99 remain unchanged from the original graph @xmath0 .
if there is no line between two parts , then there is no edge between them in @xmath101.,width=151 ]    let @xmath101 be the graph defined as above . note that the number of vertices in @xmath101 is at most @xmath116 , @xmath101 can be constructed from @xmath0 in polynomial time and @xmath101 has only @xmath102 precolored vertices .
therefore , it remains to prove @xmath104 if and only if @xmath95 is satisfiable for some @xmath100 .
we set @xmath100 to @xmath117 , since there will be at most @xmath114 colors used in the supernode .    assuming that @xmath95 is not satisfiable
, it is straightforward to design a winning strategy for on @xmath101 ; we only need to adapt the approach of sections  [ sec : largeprecol ] and  [ sec : smallprecol ] .
see appendix  [ app : noprecolwinstratd ] for a full description of the strategy .    in the rest of this section we focus on the opposite direction : assuming that @xmath95 is satisfiable , we show that can color @xmath101 with @xmath100 colors regardless of the strategy of .    in the following , when we refer to the colored part of @xmath101 , we do not take precolored vertices into account .
actually does not look at precolored vertices unless it uses its winning strategy for coloring @xmath0 , which exists by the assumptions .
* intuition .
* at the beginning has too little data to infer anything about the vertices .
therefore , shall wait for two nonadjacent vertices from @xmath97 and for two large cliques ( larger than @xmath118 ) with a small intersection .
before such vertices arrive , it will color greedily .
note that the greedy coloring algorithm eventually stops before everything is colored .
having two large cliques , one mostly from @xmath77 and the other mostly from @xmath113 , and two nonadjacent vertices from @xmath97 , is able to recognize where an incoming vertex belongs .
therefore , can use the supernode like a precolored vertex and colors the remaining vertices from @xmath97 and @xmath99 by its original winning strategy on @xmath0 .
this approach may fail if a part of @xmath97 is already colored by
s application of the greedy strategy . to remedy this
, we prove that colors used on @xmath97 so far are also used in @xmath108 or @xmath99 , or will be used on @xmath108 later .
the other obstacle is that might not be able to distinguish between one clique from @xmath97 and vertices in @xmath77 if nothing from @xmath112 arrives .
nevertheless , each vertex @xmath24 in such a `` hidden '' clique is connected to all other colored vertices in @xmath97 and to the whole colored part of @xmath99 , otherwise it would be distinguishable from vertices in @xmath77 .
hence , it does not matter on the color of @xmath24 .    in summation , the sheer size of the supernode should allow the player to be able to use it as if it would be precolored .
still , we need to allow for some small margin of error .
this leads us to the following definition :    [ dfn : almost ] let @xmath79 be the number of vertices of @xmath111 as in the construction of @xmath101 . for subgraphs @xmath119
, we say that @xmath120 is _ practically a subgraph _ of @xmath121 if @xmath122 , and @xmath120 is _ practically disjoint _ with @xmath121 if @xmath123 .
we also say that a vertex @xmath1 is _ practically universal _ to a subgraph @xmath124 if it is adjacent to all vertices in @xmath120 except at most @xmath79 of them .
similarly , we say a vertex @xmath1 is _ practically independent _ of a subgraph @xmath124 if @xmath1 has at most @xmath79 neighbors in @xmath120 .    at first
, the player uses the following algorithm for coloring incoming vertices , which may stop when it detects two useful vertices @xmath125 and @xmath126 :    while the algorithm may seem to use a huge amount of computation for one step , we should realize that we are not concerned with time complexity when designing the strategy for . in fact , even a non - constructive proof of existence of a winning strategy would be enough to imply existence of a pspace algorithm for finding it  we have observed already in section [ sec : intro ] that online chromatic numberlies in pspace .
let @xmath1 be the incoming vertex @xmath24 when stops ; note that @xmath1 is not colored by the algorithm and @xmath1 can be from any part of @xmath101 .
one of the cliques @xmath127 and @xmath128 is practically a subgraph of @xmath129 and we denote this clique by @xmath130 .
the other clique must be practically a subgraph of @xmath77 and we denote it by @xmath131 .
( keep in mind that both cliques may contain up to @xmath79 vertices from @xmath115 . )
we remark that some vertices from @xmath108 must have arrived , as @xmath77 and @xmath112 alone are indistinguishable by step  [ wfd : if ] of . by the same argument
, the player knows whether @xmath132 or @xmath133 .
let @xmath125 and @xmath126 be the nonadjacent vertices that caused the algorithm to stop .
we observe that @xmath134 by eliminating all other possibilities :    neither of @xmath125 and @xmath126 can be from @xmath99 , since any vertex of @xmath99 is practically universal to both cliques .    both @xmath125 and @xmath126 can not be from @xmath113 or both from @xmath77 , as they would be adjacent .
if @xmath125 is in @xmath113 and @xmath126 in @xmath77 , then we have a contradiction with the fact that @xmath125 and @xmath126 are not practically universal to the same clique .
if @xmath135 and @xmath126 would be from @xmath77 or @xmath112 , then @xmath125 and @xmath126 are adjacent .
finally , if @xmath135 and @xmath136 , then the clique to which they are not practically universal can not be the same for both , since @xmath125 is universal to the whole @xmath77 and @xmath126 to the whole @xmath113 .    having cliques @xmath137 and @xmath130 and vertices @xmath134 , uses the following rules to recognize where an incoming or an already colored vertex @xmath24 belongs :    if @xmath24 is practically universal to both @xmath130 and @xmath137 , then @xmath138
.    if @xmath24 is practically universal to @xmath130 and practically independent of @xmath137 and @xmath24 is adjacent to @xmath125 , then @xmath139 .    if @xmath24 is practically universal to @xmath130 and practically independent of @xmath137 , but there is no edge between @xmath125 and @xmath24 , then @xmath140 .
if @xmath24 is not practically universal to @xmath130 , but it is practically universal to @xmath137 , then @xmath141 or @xmath142 .    among such vertices , if there is a vertex not adjacent to @xmath24 or @xmath24 is not adjacent to a vertex in @xmath99 or @xmath24 is adjacent to a vertex in @xmath112 , then @xmath142 ; we say that such @xmath24 is _ surely in _
@xmath97 .
otherwise , the player cannot yet recognize whether @xmath141 or @xmath142 .
the reader should take a moment to verify that indeed , the set of rules covers all possible cases for @xmath24 .
let @xmath143 be the colored parts of @xmath101 when stops .
we observe that in the last case of the recognition the vertices from @xmath144 which are indistinguishable from @xmath77 form a clique ; we denote it by @xmath145 .
note that all vertices in @xmath145 are connected to all vertices surely from @xmath97 that arrived and @xmath145 contains all vertices in @xmath144 that are not surely in @xmath97 .
we stress that does not know @xmath145 or even its size .    * intuition for the next step . *
since can now recognize the parts of the construction ( with an exception of @xmath145 ) , we may basically use the winning strategy for on @xmath0 and on the rest .
more precisely , creates a virtual copy of @xmath0 , adds vertices into it and simulates the winning strategy on this virtual graph .
our main problem is that some part of @xmath97 ( namely @xmath144 ) is already colored .
we shall prove that if @xmath144 is not a clique , can ignore colors used in @xmath144 ( but not the colors that it will use on @xmath97 ) , since they are already present in @xmath108 or @xmath99 or they may be used later in @xmath108 .
if @xmath144 is a clique , it may be the case that @xmath108 and @xmath77 arrived completely and have the same colors , thus cannot ignore colors used on @xmath144 .
another obstacle in the simulation is @xmath145 , the hidden part of @xmath97 . to overcome this , tries to detect vertices in @xmath145 and reclassify them as surely in @xmath97 .
shall keep that all vertices in @xmath145 are connected to all currently colored vertices in @xmath97 and @xmath99 , therefore it does not matter much on colors in @xmath145 .    when discovers a vertex from @xmath145 , it adds the vertex immediately to its simulation of @xmath0 . on the other hand ,
the size of @xmath145 increases when sends a vertex from @xmath97 which is indistinguishable from @xmath77 .
the precise algorithm used by to finish the coloring of @xmath101 using @xmath100 colors and its analysis can be found in appendix  [ app : noprecolwinstratp ] .
* acknowledgments . *
the authors thank christian kudahl and their supervisor ji sgall for useful discussions on the problem .
_ effective coloration_. the journal of symbolic logic , vol .
469 - 480 ( 1976 ) .    .
_ first fit and on - line chromatic number of families of graphs_. ars combinatorica 29c , 168176 ( 1990 ) .    .
_ on - line graph coloring and finite basis problems_. combinatorics : paul erdos is eighty volume 1 . , 207214 ( 1993 ) .
j. algorithms , 23 , 265280 ( 1997 ) .
the electronic journal of combinatorics , 7(1 ) , r7 ( 2000 ) .
theoretical computer science 130.1 : 163 - 174 ( 1994 ) .    israel journal of mathematics 105 , 93 - 104 ( 1998 ) .    .
_ on - line graph coloring .
_ master s thesis , university of southern denmark ( 2013 ) .    .
to appear in proceedings of _ 9th international conference on algorithms and complexity _
( ciac 2015 ) . also arxiv 1406.1623 .
annals of discrete mathematics 43 : 319 - 325 ( 1989 ) .
the number of colors allowed for (the same as the size of @xmath31 ) is @xmath146 where @xmath28 is the number of clauses , @xmath147 the number of universally quantified variables and @xmath148 the number of existentially quantified variables .
the next two lemmas contain the analysis of our construction .
if anytime during the game uses a color not present in the color clique @xmath31 , the lemma is proven .
we therefore assume this does not happen and show that the vertex @xmath62 can not be colored with any color present in the precolored clique @xmath31 .
s strategy is to first present the vertices @xmath149 in the same order of the variables as in the formula . whenever sends a vertex @xmath34 or @xmath35 for a universally quantified variable @xmath33 , is
not able to detect whether it is setting the value of @xmath33 to @xmath40 or @xmath50 .
as the formula is not satisfiable , can therefore present these vertices in a sequence such that chooses the value of @xmath33 so that the final evaluation is false .
we now know that all clauses are evaluated to @xmath50 .
this means that at least one of the vertices @xmath58 of each clause @xmath7 is assigned the color @xmath59 .
it follows that the vertex @xmath57 of this clause must be assigned the color @xmath60 .
this holds for all clauses , thus the vertex @xmath62 has a neighbor of the color @xmath60 for each clause @xmath7 .
hence , @xmath62 can not get any of the @xmath3 colors allowed and needs to be assigned a new color , which completes the proof .        if presents a vertex , is able to recognize it by the set of edges to the clique @xmath31 , i.e. , by the set of allowed colors , with the exception of vertices @xmath78 for a universally quantified variable @xmath33 .
if presents a vertex of type @xmath64 or @xmath62 before any of the universal vertices @xmath78 is presented for a universally quantified variable @xmath33 , will be able to distinguish the vertex @xmath34 from the vertex @xmath35 and therefore choose the assignment of @xmath33 .
furthermore , for universal vertices @xmath78 and for a vertex @xmath1 for a variable @xmath42 with @xmath150 , if presents @xmath1 before any of these universal vertices is sent , is then able to distinguish the universal vertices @xmath34 and @xmath35 .
the second observation is easy to see by noting that presenting @xmath1 means that it has an edge with @xmath34 , but not with @xmath35 ( or , in the case of @xmath51 , it may have an edge with @xmath35 , but not with @xmath34 , but this is symmetric ) .
this means that we can distinguish @xmath34 from @xmath35 by the presence or absence of the edge .
we continue with the proof of the lemma .
we say that a variable @xmath33 is _ set _ , if at least one of its vertices was colored or if has assigned a value to it ( e.g. , if a vertex of a variable @xmath42 was sent before a vertex of @xmath33 for @xmath151 ) ; otherwise @xmath33 is _
unset_. until presents a vertex of type @xmath152 , or @xmath62 , colors an incoming vertex @xmath1 that corresponds to a variable @xmath33 using the following strategy :    let @xmath153 be the set of all unset variables @xmath42 with @xmath154 ,
i.e. , all previous unset variables .
chooses an assignment for all variables in @xmath153 according to satisfiability of the formula .
then can color vertices for variables in @xmath153 according to this assignment , since it is now able to distinguish vertices of each universally quantified variable in @xmath153 .
therefore every variable in @xmath153 becomes set .
if @xmath33 is not set and it is quantified existentially , knows how to set this variable to satisfy the formula , thus colors @xmath1 according to the value of @xmath33 , since @xmath1 can be recognized .      if @xmath33 is set and cannot recognize @xmath1 , then @xmath33 must be quantified universally and the other vertex for @xmath33 is colored , therefore there is a single color left for @xmath1 .
( note that in this case no vertex for a variable @xmath42 with @xmath155 arrived . )
when the vertex @xmath24 of type @xmath64 , or @xmath62 is sent , let @xmath153 be the set of all unset variables .
chooses an assignment for all variables in @xmath153 according to satisfiability of the formula .
then can color vertices for variables in @xmath153 according to this assignment , since it is now able to distinguish vertices of a universally quantified variable in @xmath153 . then all variables are set and decides how to color all remaining vertices in the graph using the following rules :        a vertex of type @xmath57 obtains color @xmath59 if all of its adjacent vertices of type @xmath51 have colors of type @xmath156 , i.e. , the clause @xmath7 evaluates to @xmath40 ; otherwise it obtains color @xmath60 .
[ o : eachnode2colors ] if reveals all the nodes @xmath81 before it sends any vertex from @xmath26 ( the gadgets and @xmath31 ) , can force to use two colors in the lower partite set of each node and these colors are different for each node .
thus uses the following strategy for each node @xmath77 independently : first presents the vertex @xmath75 from @xmath77 ; let @xmath157 be its color
. then sends a vertex @xmath158 that is one of vertices @xmath73 and @xmath74 .
cannot deduce which of them is @xmath158 , because both are not connected to @xmath75 ( and both are connected to all vertices in other nodes ) . if assigns the color @xmath157 to @xmath158 , then @xmath159 and sends @xmath73 which must get another color than @xmath157 . otherwise if @xmath158 obtains another color than @xmath157 , then @xmath160 .
[ o : nodestwocoloredbygreedy ] suppose that is using , i.e. , always assigns the smallest color not present among colored vertices adjacent to the incoming vertex .
then there are at most two colors used on each node .
@xmath161 : let the color of @xmath73 be @xmath157 and the color of @xmath74 be @xmath162 ( clearly @xmath163 ) .
since the edges from @xmath75 and @xmath73 to other vertices in the construction except @xmath74 are exactly the same , can use @xmath157 for @xmath75 , but he can not use any color forbidden for @xmath73 when @xmath73 was colored with the possible exception of @xmath162 . hence @xmath75 obtains @xmath157 or @xmath162 , but in both cases the node @xmath77 has two colors .
@xmath164 : let the color of @xmath75 be @xmath157 and the color of @xmath74 be @xmath162 .
if @xmath163 , uses @xmath157 for @xmath73 for the same reason as in the previous case .
otherwise @xmath165 and must color @xmath73 with another color , but there are two colors on @xmath77 .
@xmath166 : we show that @xmath75 and @xmath73 have the same color . without loss of generality , @xmath75 arrives first and obtains a color @xmath157 .
when @xmath73 arrives , it can be colored by @xmath157 , but by no other color @xmath167 , since @xmath168 would be available for @xmath75 also when @xmath75 was colored , because the edges from @xmath75 and @xmath73 to other vertices except @xmath74 are the same .
hence uses another color for @xmath74 and @xmath77 is 2-colored .
has @xmath169 colors for the @xmath79 nodes and @xmath3 colors for @xmath26 ( with the same names as in the previous section ) , thus overall is allowed to use @xmath170 colors .
the precolored vertices have a color that may be used later in @xmath26 .
if at any point of the game there are @xmath171 colors in the graph , the lemma is proven .
s strategy is to first present all nodes ( in any order ) and force to use two colors in the lower partite set of each node .
moreover , has to use different colors in distinct nodes .
forcing such a coloring is possible by observation  [ o : eachnode2colors ] .    then sends @xmath31 and the situation is similar to the one with the precolored @xmath31 ,
since none of two colors used in the lower partite set of a node is allowed for a vertex in @xmath26 , thus must use the @xmath3 colors in @xmath31 for coloring the gadgets . is
able to recognize vertices in @xmath31 by nodes , but nodes does not give an additional knowledge compared to @xmath31 .
in particular , is not able to distinguish vertices @xmath78 for a universally quantified variable @xmath33 when one of them arrives .
hence we conclude the proof by applying lemma  [ l : precolnotsat ] .    in the case of a satisfiable formula , nodes and
precolored vertices become important .
 s strategy to win is basically the following : color vertices in nodes greedily until a vertex @xmath24 from @xmath26 arrives .    at this time ,
some ( but maybe not all ) vertices in the gadgets can be recognized by their nodes .
uses the winning strategy from lemma  [ l : precolsat ] on the gadgets for vertices that it can recognize , even if @xmath31 has arrived only partially . for each vertex
@xmath1 that can not be recognized , is able to color ( or already colored ) the lower partite set of the @xmath1 s node @xmath77 with only one color , therefore it can use the same color for @xmath1 and the vertex @xmath172 , since they are not adjacent .    in the following proof , waits for two nonadjacent vertices in @xmath26 , even although one vertex from @xmath26 suffices .
the reason is the we shall need such a condition in section  [ sec : noprecol ] when we remove precolored vertices .
let @xmath173 be the set of @xmath169 colors for nodes and let @xmath174 be the set of @xmath3 colors for @xmath26 ( the gadgets and @xmath31 ) including the color used on precolored vertices ; both @xmath173 and @xmath174 are ordered arbitrarily and @xmath175 .
let @xmath176 and @xmath177 be the two nonadjacent vertices from the stopping condition of .
observe that the nodes identifying @xmath176 and @xmath177 are different , since these nodes can be the same only for both universal vertices for one variable , but these vertices are adjacent .      if @xmath24 is from @xmath26 and can recognize it by a nonadjacent vertex @xmath1 from nodes ( line  [ alg - ln : win - precol - g1vertrecognized ] ) then coloring it with the winning strategy is correct , since vertices that have a color from @xmath174 in @xmath71 have the same color in @xmath178 .    otherwise ,
if cannot recognize @xmath24 from @xmath26 ( line  [ alg - ln : win - precol - g1vertnotrecognized ] ) , the vertex @xmath74 from the node @xmath77 that serves for identifying @xmath24 is not yet sent .
it follows that there is at most one color used in the lower partite set of @xmath77 by observation  [ o : nodestwocoloredbygreedy ] .
if @xmath24 is from a node @xmath77 , a vertex @xmath1 from @xmath26 identified by @xmath77 is colored by @xmath157 , there is no edge between @xmath1 and @xmath24 , @xmath179 and @xmath157 is not present among nodes ( line  [ alg - ln : win - precol - nodessavecolor ] ) , then @xmath24 is the vertex @xmath74 of the node @xmath77 .
there is at most one color used in the lower partite set of @xmath77 before coloring @xmath24 by observation  [ o : nodestwocoloredbygreedy ] and that @xmath157 is not used anywhere except @xmath1 .
an unrecognized vertex @xmath1 from @xmath26 which is not a universal vertex has the same color as the vertex @xmath74 from the single node identifying @xmath1 , since @xmath77 identifies only @xmath1 and other vertices are adjacent to @xmath74 .
an unrecognized universal vertex @xmath1 identified by two nodes @xmath77 and @xmath180 obtains some color @xmath179 .
let @xmath162 be the color of the other vertex @xmath181 for the same universal variable .
vertices @xmath74 from @xmath77 or @xmath180 arrive after @xmath1 and one of them obtains @xmath157 , since the algorithm prefers to use colors of @xmath1 and @xmath181 if their colors are in @xmath173 .
recognized vertices from @xmath26 are colored by the wining strategy on @xmath178 using @xmath3 colors from @xmath174 .
as we have shown , each unrecognized vertex @xmath1 has the same color as a vertex @xmath74 in one of the nodes identifying @xmath1 .
there are at most @xmath169 colors used in nodes by observation  [ o : nodestwocoloredbygreedy ] and all of them are from @xmath173 .
thus uses at most @xmath170 colors .
suppose that @xmath95 is not satisfiable and consider @xmath101 constructed from @xmath95 .
starts by sending @xmath77 and then @xmath113 such that no color is both in @xmath77 and in @xmath112 ; this is possible , since cliques @xmath77 and @xmath108 have the same size .
more precisely , if would assign a color from @xmath77 to a vertex from @xmath113 , decides that the vertex is in @xmath108 ; otherwise the vertex is in @xmath112 .
this forces to use @xmath114 different colors on @xmath77 and @xmath112 and no such color can be used for @xmath97 and @xmath99 .
now the colored part of @xmath101 ( precolored vertices and the supernode ) may be viewed as a precolored part in @xmath0 except that may not be able to use the supernode in the same way as a precolored vertex .
since this only helps , we use the `` only if '' direction of the first assumption , i.e. , sends vertices according to its winning strategy for @xmath0 .
this proves the `` only if '' direction of the first proposition of lemma  [ l : removeprecolvertex ] .
after the algorithm finishes , the player applies an algorithm which simulates its winning strategy on @xmath0 and maintains disjoint sets of colors @xmath173 ( for the supernode and @xmath97 ) and @xmath182 ( mainly for @xmath99 ) .
recall that is not able to distinguish between @xmath77 and @xmath145 , thus vertices in @xmath145 are treated in the same way as vertices in @xmath77 .
more precisely , we define three color sets @xmath173 , @xmath183 and @xmath182 as follows :    if colored vertices surely in @xmath97 form a clique , let @xmath183 be the set of colors that are used on a colored vertex surely in @xmath97 and that are not used in @xmath184 .
otherwise , if there are two nonadjacent colored vertices surely in @xmath97 , let @xmath183 be an empty set .
let @xmath173 be the set of colors present currently in the supernode or in @xmath145 and among vertices surely in @xmath97 ( i.e. , in @xmath185 , or @xmath144 ) except colors from @xmath183 and colors present also in @xmath186 .
we will update the color sets when we apply our algorithms and as more vertices arrive on input ; the precise updating procedure is specified throughout the algorithms . the winning algorithm for is split into two parts , initialization and coloring :
when adding a colored vertex @xmath24 that is surely in @xmath97 to the virtual graph @xmath178 , if @xmath187 , coloring @xmath24 with @xmath157 in @xmath178 does not harm the winning strategy by the assumptions of lemma  [ l : removeprecolvertex ]  if @xmath183 is nonempty , colored vertices surely in @xmath97 form a clique and may use any proper coloring for them in the winning strategy . moreover
, some colors may be both in @xmath144 and in @xmath186 , but the winning strategy still works by the assumptions . otherwise ,
if @xmath188 , @xmath24 may obtain a different color in @xmath178 than in the actual @xmath101 .
we remark that the first vertex @xmath24 colored by the algorithm is @xmath1 , the vertex on which stops . in a sense , if @xmath144 , the part of @xmath97 colored by , is not a clique , pretends that in @xmath144 there are colors from the simulation , not the colors assigned by . otherwise , if @xmath144 is a clique , uses colors from @xmath101 in @xmath178 which may result in renaming the colors in the winning strategy on @xmath178 .
observe that any color used in @xmath101 is in @xmath173 or @xmath182 and that these sets of colors are disjoint .
thus the colors used for @xmath115 are different from those in the supernode except colors used in @xmath144 or colors of vertices that were in @xmath145 when they arrived .
observe that if @xmath24 is from @xmath99 or surely from @xmath97 , coloring with @xmath157 is sound ; this follows from the winning strategy on @xmath0 and from the fact that colors in @xmath182 are not used in the supernode ( recall that colors in @xmath183 are not used in @xmath184 ) . it remains to show the following claim :    [ clm : colorsfromdinc ] for each vertex @xmath24 from @xmath144 ( either surely from @xmath97 or from @xmath145 when starts ) colored with @xmath157 in @xmath101 , the color @xmath157 is present in @xmath178 , or in @xmath108 when the whole @xmath101 is colored .
we assume that @xmath157 is not used in @xmath178 , otherwise the claim holds for @xmath24 . colors in @xmath144 and not in @xmath178 ( thus also not in @xmath183 ) are in @xmath173 . observe that @xmath157 can be used only for a vertex in @xmath108 .
let @xmath189 be the time of coloring @xmath24 .
if at time @xmath189 there is some color in @xmath108 not used in @xmath77 , let @xmath168 be the smallest such color in the ordering of colors .
we observe that uses the color @xmath168 for @xmath24 , since other colors allowed for coloring @xmath24 are some colors used in @xmath99 at time @xmath189 ( which we assume that @xmath24 does not get ) , or a new color not yet used anywhere , but the new color is greater than @xmath168 .
hence the greedy algorithm prefers @xmath168 and @xmath190 .    otherwise , when colors @xmath24 , all colors used in @xmath108 are also present in @xmath77 which also means that there are at least as many colored vertices in @xmath77 as in @xmath108 at time @xmath189 . in this case
we show that @xmath157 will be used in @xmath108 after time @xmath189 .
suppose that at least @xmath79 vertices from @xmath108 are colored by after @xmath24 arrives .
then , must assign @xmath157 to a vertex in @xmath108 at some point , since it prefers colors used in @xmath144 before colors in @xmath77 or new colors , plus we know at least @xmath191 vertices from @xmath108 will arrive after @xmath24 .    assuming the opposite , more than @xmath192 vertices arrived from @xmath108 before time @xmath189 , which means that also more than @xmath192 vertices arrived from @xmath77 . since @xmath24 is colored by
, there can not be two nonadjacent vertices from @xmath97 at time @xmath189 , because otherwise the algorithm would stop .
observe that in this case stops immediately after there are two nonadjacent vertices in @xmath97 , thus @xmath144 is a clique and in particular , @xmath157 is in @xmath178 or @xmath157 was used in @xmath184 which concludes the proof .
it follows that all colors in @xmath173 are used in the supernode .
there are always @xmath114 colors used in the supernode , since otherwise there must be a color @xmath157 in @xmath77 that is not in @xmath113 , but both and prefer coloring @xmath112 using colors in @xmath77 and coloring @xmath77 using colors from @xmath113 . thus @xmath193 .
since all colors in @xmath182 are used in the virtual graph @xmath178 ( this includes also colors in @xmath183 ) and uses the winning strategy on @xmath178 , we know that @xmath194 .
overall , uses at most @xmath195 colors and we proved @xmath104 for @xmath101 constructed from a satisfiable formula @xmath95 .
this concludes the proof of lemma  [ l : removeprecolvertex ] .
let @xmath95 be a formula of size @xmath20 .
we first construct a graph @xmath71 with @xmath196 precolored vertices @xmath91 as described in section  [ sec : largeprecol ] and  [ sec : smallprecol ] . then for each precolored vertex , we apply lemma  [ l : removeprecolvertex ] iteratively until we obtain a graph @xmath197 with no precolored vertex .
we check that the graph satisfies both assumptions of lemma  [ l : removeprecolvertex ] in each iteration @xmath32 . the first assumption , @xmath5 if and only if @xmath95 is satisfiable , holds for the first iteration by lemmas  [ l : smallprecolsat ] and  [ l : smallprecolnotsat ] , and for every other iteration by the previous use of lemma  [ l : removeprecolvertex ] .
recall the second assumption : in the winning strategy of in the case if @xmath95 is satisfiable , can color @xmath99 using before two nonadjacent vertices from @xmath97 arrive .
moreover , if assigns the same color to a vertex in @xmath97 and a vertex in @xmath99 before two nonadjacent vertices from @xmath97 arrive , can still color @xmath0 using @xmath3 colors .
let @xmath198 be @xmath71 and for each iteration @xmath32 let @xmath199 be the graph before removing the @xmath32-th precolored vertex @xmath200 .
let the nonprecolored vertices of @xmath199 be partitioned into two disjoint induced subgraphs @xmath201 and @xmath202 such that all vertices from @xmath202 are connected to @xmath200 and no vertex in @xmath201 is connected to @xmath200 .
let @xmath203 be the graph created from @xmath199 by the construction in lemma  [ l : removeprecolvertex ] ; we denote by @xmath204 and @xmath205 the cliques in the supernode .
thus @xmath203 consists of parts @xmath206 , and @xmath202 ( and possibly some precolored vertices which we do not take into account ) .
note that @xmath202 contains only some nodes and @xmath201 contains the whole @xmath26 ( i.e. , @xmath31 and the gadgets for variables and clauses ) together with supernodes from previous iterations and nodes that are not in @xmath202 .
the first part of the second assumption holds for the first iteration , since the algorithm from the winning strategy from lemma  [ l : smallprecolsat ] colors the nodes by before two nonadjacent vertices from @xmath26 arrive and @xmath207 contains @xmath26 .    for
the `` moreover '' part in the first iteration , if assigns the same color @xmath157 to a vertex @xmath162 in @xmath207 and a vertex @xmath208 in @xmath209 before two nonadjacent vertices from @xmath207 arrive , then @xmath208 must be from a node that identifies @xmath162 and in this case can just pretend that @xmath162 is an unrecognized vertex .
in other words , the color @xmath157 is intended for nodes and not used in @xmath26 .    for an iteration @xmath210 , recall that @xmath203 consists of parts @xmath206 , and @xmath202 ( except precolored vertices ) and that @xmath201 and @xmath202 form the graph @xmath199 from the previous iteration .
similarly , @xmath199 consists of @xmath211 , and @xmath212 .
note that @xmath213 , and @xmath214 are in @xmath201 and that there may be some nodes both in @xmath201 and in @xmath212 and also some nodes both in @xmath202 and in @xmath215 .    the first part of the second assumption holds for the iteration @xmath32 , because in the winning strategy on @xmath199 colors @xmath199 using before two large and practically disjoint cliques from the supernode arrive , one mostly from @xmath216 and the other mostly from @xmath217 .
since there are many pairs of nonadjacent vertices between these cliques and these cliques are contained in @xmath201 , can color @xmath202 using before two nonadjacent vertices from @xmath201 arrives .
for the `` moreover '' part , if assigns the same color @xmath157 to a vertex @xmath162 in @xmath201 and a vertex @xmath208 in @xmath202 before two nonadjacent vertices from @xmath201 arrive , then @xmath208 is from nodes ( as @xmath202 contains only nodes ) and @xmath162 is either from @xmath26 as in the first iteration , or from the supernode created in the previous iteration ; in the latter case , @xmath208 is in @xmath215 and @xmath162 is in @xmath214 , since otherwise they would be connected .
( note that @xmath162 can not be from nodes , since there is complete bipartite graph between every two nodes and a node is never split between @xmath99 and @xmath97 . )    in both cases , this does not harm the winning strategy of on @xmath199 , because the algorithm assumes that may use the same colors for more vertices in @xmath218 or in @xmath219 .
this concludes that the assumptions of lemma  [ l : removeprecolvertex ] are satisfied in each iteration .
since @xmath220 iff @xmath95 is satisfiable , it holds that @xmath221 iff @xmath95 is satisfiable by lemma  [ l : removeprecolvertex ] ( for some @xmath100 ) .
the number of vertices in @xmath197 is polynomial in @xmath20 , because @xmath71 has linearly many vertices and the number of vertices is multiplied by a constant with each of logarithmically many applications of lemma  [ l : removeprecolvertex ] .
the constructions in sections  [ sec : largeprecol ] and  [ sec : smallprecol ] and lemma  [ l : removeprecolvertex ] yield a polynomial - time algorithm for computing @xmath101 from @xmath95 .


in recent years , the study of entanglement in quantum many - body systems has become a common issue in condensed matter physics , statistical mechanics , and quantum information theory .
entanglement is a purely quantum mechanical phenomenon in which quantum systems are linked together even when they are spatially separated so that one system can not be correctly described without full mention of its counterpart .
this concept was first introduced by schrdinger in  @xcite . from the viewpoint of quantum information science
, entanglement is a fundamental measure of how much quantum effects we can observe and use to control one quantum system by another , and it is the primary resource for quantum information processing and communication  @xcite . on the other hand , in condensed matter physics and statistical mechanics , the concept of quantum entanglement has recently been used to investigate quantum phase transitions  @xcite ( see  @xcite for a review ) , topological order  @xcite , and macroscopic properties of solids  @xcite .
a fundamental and practical problem common in all the areas is how to detect entanglement  @xcite and quantify the degree of it .
a general approach is to study the reduced density matrix of a certain subsystem in an entangled state .
the concept of the reduced density matrix was first introduced by dirac in  @xcite .
there are many kinds of measures of entanglement related to the reduced density matrix ( see  @xcite and references therein ) . among them , one of the most popular characteristic functions is the von neumann entropy ( entanglement entropy ) of a reduced density matrix  @xcite .
it serves as the measure of entanglement for a pure state .    generally speaking , to elucidate quantum entanglement in many - body systems is a formidable task since we need a precise description of many - body quantum states .
so far , much insight has been obtained by studying exactly solvable models in which there is a possibility to obtain the reduced density matrix and calculate the von neumann entropy exactly or approximately .
representative models are the harmonic models  @xcite , the xy spin chain  @xcite , the xxz spin chain  @xcite , the calogero - sutherland model  @xcite , one - dimensional critical models  @xcite , conformal field theories  @xcite , the affleck - kennedy - lieb - tasaki ( aklt ) model  @xcite , the kitaev model  @xcite , and the quantum dimer ( rokhsar - kivelson ) model  @xcite . among them
, the aklt model  @xcite has some important meanings since this model is directly related to one of the schemes of quantum computation , measurement based quantum computation  @xcite .
historically , the aklt model and its exact ground state known as the valence - bond - solid ( vbs ) state were proposed to understand ground - state properties of the haldane gap systems  @xcite .
the authors ( aklt ) rigorously proved that the vbs state is a unique ground state and there is an energy gap immediately above the ground state .
exponential decay of correlation functions in the vbs ground state has also been shown .
there is an intimate relation between the vbs states and fractional quantum hall states which has been revealed by using a spin coherent state representation  @xcite .
a hidden topological order in the vbs state is characterized by string order parameters  @xcite and existence of edge states which are degenerate ground states in a chain with open boundary conditions  @xcite .
the vbs picture and the edge states were also experimentally supported by electron spin resonance ( esr ) in the @xmath1 haldane chain ni(c@xmath2h@xmath3n@xmath2)@xmath2no@xmath2(clo@xmath4 ) ( nenp )  @xcite .
although most of the exactly solvable models are limited to one dimension , the aklt model and vbs ground state can be generalized to higher dimensions and/or inhomogeneous ( non - translational invariant ) lattices  @xcite .
this feature is especially important in connection with the measurement based quantum computation as we discussed later . for studying nature of the vbs states ,
the non - trivial observation has been done by klmper _
et al_. , through the @xmath5-deformation of the aklt model  @xcite .
they found that the one - dimensional ( 1d ) spin-1 vbs state can be expressed in a form of the matrix product state ( mps ) .
this representation gives an efficient method to calculate correlation functions . mathematical theory and generalization of the vbs states ,
i.e. , finitely correlated states ( fcs ) , were essentially developed by fannes _
et al_. , in  @xcite .
another interesting generalization of the vbs state is to replace su(2 ) spins with other lie algebra generators . in this direction , there have been constructed many versions of the aklt model such as su(3 )  @xcite and su(@xmath6 ) modes  @xcite , su(@xmath7)-extended model  @xcite , symplectic model  @xcite , so(5 ) model  @xcite , so(@xmath8 ) model  @xcite , supersymmetric model  @xcite , and several @xmath5-deformations  @xcite .
the vbs state has recently been attracting renewed interest from the viewpoint of quantum information theory .
in particular , it has been proposed that the vbs state can be used as a resource state in measurement - based quantum computation ( mqc ) .
mqc is one of the schemes of quantum computation in which we can perform universal quantum computation using only measurements as computational steps ( see  @xcite for an introduction ) .
there are two kinds of measurement based models .
one is the teleportation quantum computation ( tqc )  @xcite .
the other is the so - called one - way quantum computation ( 1wqc ) or cluster state computation  @xcite . in this model ,
any quantum gate array can be implemented as a pattern of single - qubit measurement on a highly entangled cluster state .
note that 1wqc with one - dimensional cluster states is insufficient for universal quantum computation and thus two - dimensional cluster states are inevitably needed .
it has then been proved in  @xcite by verstraete and cirac that tqc and 1wqc are equivalent using the vbs picture .
gross , eisert and collaborators have introduced novel schemes for mqc based on finitely correlated or tensor network state beyond the cluster state  @xcite .
another recent important development by brennen and miyake is that the vbs state ( more precisely dynamically coupled aklt chains ) can be used as a resource state in mqc instead of the cluster state  @xcite .
the vbs state is more advantageous to 1wqc than the cluster state since it is robust against noise due to the energy gap .
implementations of the aklt hamiltonian using bosonic atoms  @xcite or spin-1 polar molecules  @xcite in optical lattices have also been proposed .    in this paper , motivated by considerable current interest from several viewpoints ,
we study entanglement properties of vbs states on two - dimensional ( 2d ) graphs which is relevant to mqc . before discussing the 2d vbs states ,
let us briefly summarize the obtained results in 1d vbs states .
the reduced density matrix and entanglement properties in the 1d spin-1 vbs state was first studied in  @xcite .
the authors found that the von neumann entropy of the block with the rest of the chain approaches a constant value exponentially fast , as the size of the block increases .
the constant value is given by @xmath9 .
then , the 1d integer-@xmath10 vbs states were studied in  @xcite . in this case
, the von neumann entropy of the large block with the rest is given by @xmath11 and this value is interpreted in terms of the edge states , the degenerate ground states of the open aklt chain .
this picture leads to an idea that the reduced density matrix of the block in the aklt chain is _ exactly _ spanned by the degenerate ground states of the block hamiltonian ( see eq .
( [ block_ham ] ) for the definition ) .
this idea has been confirmed in the series of papers  @xcite and was summarized in the review  @xcite .
another interesting quantities such as the entanglement length  @xcite , the geometric entanglement  @xcite , the single - copy entanglement  @xcite , and the boundary effect  @xcite in the 1d spin-1 vbs state have also been studied .
more generally , the entanglement of formation in the 1d fcs was estimated in  @xcite and an area law in 1d gapped quantum system was proven in  @xcite . compared to the 1d vbs states , few results have been obtained in 2d vbs states so far .
entanglement properties of the vbs state on the cayley tree  @xcite was studied in  @xcite .
the authors showed that the von neumann entropy of the block does not depend on the whole size of the system and its asymptotic value is linearly proportional to the number of valence bonds crossing the boundary .
this gives an explicit proof of the area law in this system ( see  @xcite for a review of area laws and the entanglement entropy ) .
general entanglement properties of the vbs state on an arbitrary graph  @xcite have been studied in  @xcite .
the authors showed that the eigenspace of the reduced density matrix of the block is spanned by the degenerate ground states of the block hamiltonian .
the authors have shown that the von neumann entropy of a large block with the rest approaches a value less than @xmath12 , where @xmath13 is the ground - state degeneracy of the block hamiltonian .
although some general results for 2d vbs states have been obtained so far , a quantitative analysis of entanglement was missing since a powerful analytical method such as a transfer matrix technique does not work sufficiently in the analysis of 2d systems .
therefore , entanglement properties of the 2d vbs states must be analyzed with the aid of numerics .    in this paper , we study entanglement properties of the 2d vbs states on symmetric graphs .
symmetries such as reflection and inversion enable us to develop an efficient method to study the reduced density matrix .
due to the symmetry , one can find the relation between the eigenvalues of the reduced density matrix and those of the _ overlap _ matrix ( see eqs .
( [ overlap1]-[overlap2 ] ) for the definition ) .
the overlap matrix is a real symmetric and its elements are given by the inner product of the degenerate ground states of the block hamiltonian which are linearly independent but may not be orthonormal .
the inner product can be , for example , calculated by monte carlo integration using the schwinger boson representation of the vbs state . even without numerical integration , one can formally write the overlap matrix as a matrix - valued correlation function in the vbs state restricted to one of the subsystems . from this , we conjecture that the von neumann entropy per valence bond ( boundary site ) is strictly less than @xmath0 even in the 2d infinite system .
holographic _ interpretation of the deviation from @xmath0 is given in terms of the hidden spin chain along the boundary between two subsystems ( reflection axis ) . to confirm this conjecture , we numerically study the eigenvalues of the overlap matrix by exact diagonalization . since the dimension of the overlap matrix does not depend on the size of the whole system but is given by @xmath13 , we can study relatively large systems .
from the eigenvalues of the reduced density matrix , we calculate the von neumann entropy .
we find that the von neumann entropy per valence bond on the boundary is strictly less than @xmath0 for finite size systems .
we also find a function which fits the obtained numerical results well .
the long - distance behavior of this function strongly supports our conjecture .
these numerical results are analytically or algebraically confirmed for the graphs with a ladder geometry , where transfer matrix technique can be applied .
the organization of this work is as follows . in section 2
, we review the construction of the aklt model on an arbitrary graph and introduce the schwinger boson representation of the vbs state . in section 3
, we provide a general argument on the schmidt decomposition and discuss a relation between the reduced density matrix and the overlap matrix .
we show how to apply this method to study the von neumann entropy of the vbs state .
then we make a conjecture on the von neumann entropy in the 2d infinite system . in section 4 ,
numerical results of the vbs states on the square and hexagonal lattices are shown . in section 5 , the obtained numerical results
are partly confirmed analytically for graphs in a ladder geometry .
conclusions are given in the last section . in appendix , we provide a graphical interpretation of the overlap matrix in terms of closed loops and open strands .
let us first define the _ basic _ aklt model on a connected graph . a graph consists of two types of elements , namely vertices ( sites ) and edges ( bonds ) . as shown in fig .
[ fig : vbs_basic_model ] , a vertex is drawn as a ( large ) circle and an edge is as a solid line .
has the spin @xmath14 in this graph . ]
a graph is called _ connected _ if every pair of distinct vertices in the graph can be connected through some path .
let @xmath15 be a connected graph , where @xmath16 and @xmath17 are the sets of vertices and edges , respectively . henceforth , we assume @xmath18 , otherwise there is no interaction between spins
. the spin operator @xmath19 is located at vertex @xmath20 and its spin value is denoted by @xmath21 . in the basic model , we require that @xmath22 , where @xmath23 is the coordination number , i.e. , the number of incident edges connected to @xmath20 . to guarantee uniqueness of the ground state , this relation must be true for any vertex @xmath20 , including boundaries .
let us now define the spin hamiltonian for the basic aklt model .
the hamiltonian is written as a sum of interactions on all the edges : h= _ k , l e a(k , l ) _ s_k+s_l ( k , l ) .
[ basic_aklt_ham ] here , @xmath24 is an arbitrary positive real coefficient ( it may depend on the edge @xmath25 ) and the operator @xmath26 projects the total spin on the edge @xmath25 , @xmath27 , on the subspace with the highest possible spin value @xmath28 . since the hamiltonian in eq .
( [ basic_aklt_ham ] ) is a linear combination of projection operators with positive coefficients , @xmath29 is positive semi - definite .
the explicit form of the projector @xmath30 in terms of @xmath31 and @xmath32 is given by _
s_k+s_l(k , l)=_|s_k - s_l| j s_k+s_l-1 .
[ projector_pi ] note that one can expand @xmath33 and hence the projector defined in eq .
( [ projector_pi ] ) is a polynomial in the scalar product @xmath34 of the degree @xmath35 , where @xmath36 is the minimum of @xmath21 and @xmath37 . as an example , consider a case where @xmath38 , which corresponds to the @xmath1 homogeneous aklt chain .
the projector @xmath39 is written as _ 2(k , l)= ( s_k s_l)^2 + ( s_k s_l)+ .
the projectors @xmath40 and @xmath41 are related to the hamiltonians on the hexagonal and the square lattices , respectively : @xmath42 in general , explicit expressions in terms of @xmath34 become more complicated .
let us now discuss the ground state of @xmath29 in eq .
( [ basic_aklt_ham ] ) with condition s_k = z_k .
the ground state is unique and known as the valence - bond - solid ( vbs ) state .
although there are several possible representations for the vbs state ( see  @xcite ) , the most convenient one for our purposes is the schwinger boson representation  @xcite . in this representation ,
we introduce a pair of bosonic creation and annihilation operators at each vertex to realize su(2 ) lie algebra .
we define a pair of bosonic operators @xmath43 and @xmath44 for each vertex @xmath20 with the canonical commutation relations : = [ b_k , b^_l ] = _ kl with all other commutators vanishing : = [ b_k , b_l]=[a_k , b_l]=[a_k , b^_l]=0 ,  ^(k , l ) .
spin operators are represented by the schwinger bosons as s^+_k = a^_k b_k ,  s^-_k = b^_k a_k ,   s^z_k = ( a^_k a_k - b^_k b_k ) . to reproduce the dimension of the spin-@xmath21 hilbert space at each vertex @xmath20 , the following constraint on the total boson occupation number
is required : a^_k a_k + b^_k b_k = 2 s_k . with this constraint in mind
, the vbs ground state is written as |vbs = _ k , le ( a^_k b^_l - b^_k a^_l ) |vac , [ vbs_gs ] where the product runs over all edges and the vacuum @xmath45 is defined as the direct product of vacuums of each vertex , i.e. , |vac= _ k |vac^[k ] .
note that @xmath46}\ra$ ] is destroyed by any annihilation operators @xmath43 or @xmath44 . using the schwinger boson representation
, one can generalize the vbs state which is the ground state of the _ generalized _ aklt model . by associating a positive integer @xmath47 to each edge @xmath48 of @xmath49
, the generalized vbs state is written as |vbs_gen = _ k , le ( a^_k b^_l - b^_k a^_l)^m_k , l |vac .
for the construction of the generalized aklt hamiltonian and the condition of the uniqueness of the ground state , please see original references  @xcite . in this paper , we shall focus on the basic model , i.e. , @xmath50 for any @xmath25 .
let us now prove that the vbs state in eq .
( [ vbs_gs ] ) is the zero - energy groud state of the basic aklt hamiltonian @xmath29 . to prove this , we have only to show for any vertex @xmath20 and edge @xmath25 ( i ) the total power of @xmath51 and @xmath52 is @xmath53 so that the spin value at the vertex @xmath20 @xmath21 , and ( ii ) the maximum value of the total spin of the edge @xmath25 is @xmath54 .
( i ) can be shown by expanding r.h.s .
( [ vbs_gs ] ) and finding that the total power of @xmath51 and @xmath52 is @xmath55 .
let us prove ( ii ) . by counting the power of @xmath56 s minus @xmath57 s on the bond @xmath25
, we find that the maximal eigenvalue of @xmath58 is @xmath59 , i.e. , @xmath60 .
since the state @xmath61 is invariant under global rotations , the maximum value of the total spin of the edge @xmath25 is @xmath54 .
therefore , from ( i ) and ( ii ) , the vbs state @xmath62 in eq .
( [ vbs_gs ] ) is a zero - energy ground state of the basic hamiltonian @xmath29 .
it was shown in @xcite that this ground state is unique .
then , eq . ( [ pre_schmidt ] ) can be written as @xmath80}_\tau d^{[b]}_\eta } ( { \tilde x}^{\rm t}{\tilde y})_{\tau \eta } |e_\tau\ra \otimes |f_\eta\ra , \label{pre_schmidt2}\ ] ] where @xmath81 denotes matrix transpose .
so far , we have considered a general case .
let us now focus on the special situation where @xmath69}$ ] and @xmath70}$ ] are the same , i.e. , @xmath69}=m^{[b]}=m$ ] .
furthermore , if @xmath82 is a real symmetric matrix , it can be diagonalized by some orthogonal matrix @xmath83 as @xmath84 .
in such a case , eq .
( [ pre_schmidt2 ] ) becomes a simpler form : @xmath85 note that @xmath86}_\tau = d^{[b]}_\tau$ ] .
one may think that the above situation is not generic .
however , as we will see later , it is enough for our purpose to restrict ourselves to such a special case .
let us now consider the density matrix for the total system consisting of @xmath65 and @xmath66 .
it is defined by @xmath87 . then the reduced density matrix for @xmath65
is given by @xmath88 . in our case , it is given by @xmath89 therefore , the von neumann entropy in this bipartite partitioning is given by @xmath90 therefore , we can obtain the von neumann entropy in terms of the overlap matrix @xmath82 . to obtain @xmath91 , we need all the eigenvalues of @xmath82 .
let us now apply the schmidt decomposition method to study entanglement in vbs states on graphs with symmetries such as reflection and/or inversion . for simplicity , we shall focus on the graphs with reflection symmetry .
however , one can easily generalize our argument to graphs with inversion symmetry .
we define reflection symmetry such that any site ( vertex ) in a subgraph @xmath65 has a reflection partner in @xmath66 and vice versa .
a graphical example is shown in fig .
[ fig : ref_sym ] .     and
@xmath92 denote the sets of sites ( vertices ) on the boundaries of subgraphs @xmath65 and @xmath66 .
any site @xmath93 has a unique reflection symmetric partner @xmath94 . ]    in the schwinger boson language , the vbs state on a reflection symmetric graph can be written as @xmath95 where @xmath96 and @xmath97 are sets of edges ( bonds ) in subsystems @xmath65 and @xmath66 while @xmath98 and @xmath92 are sets of vertices ( sites ) on the boundaries of @xmath65 and @xmath66 .
we note that we only consider the _
basic _ vbs state which is the ground state of the basic aklt model . in this state
, there is one valence bond on any edge @xmath99 .
the generalization of our argument to generic vbs states is straightforward . to apply the schmidt decomposition , we slightly change the above expression by a local gauge transformation : @xmath100 and @xmath101 for @xmath102 .
( [ ref_vbs ] ) becomes @xmath103 since the graph is reflection symmetric , @xmath93 necessarily has a reflection symmetric partner @xmath104 ( see fig . [
fig : ref_sym ] ) .
for simplicity , we henceforth assume that any @xmath105 is nearest neighbor to @xmath106 in the graph @xmath49 .
therefore , the above expression can be written as follows : @xmath107}_{\{\al \ } } \ra \otimes |\phi^{[b]}_{\ { \al \}}\ra,\end{aligned}\ ] ] where @xmath108 @xmath109 corresponds to the spin state of the boundary state , i.e. , degenerate ground state of the _ block hamiltonian _ defined within the subsystem @xmath65 or @xmath66 .
the block hamiltonians are explicitly defined as @xmath110 where @xmath111 is defined in eq .
( [ projector_pi ] ) .
the states @xmath112}_{\{\al \ } } \ra$ ] and @xmath113}_{\ { \al \}}\ra$ ] are given by @xmath114 respectively , where @xmath115}\ra$ ] ( @xmath116}\ra$ ] ) is the vacuum for bosons in @xmath65 ( @xmath66 ) and @xmath117}\ra \otimes |{\rm vac}^{[b]}\ra$ ] .
equations in ( [ phi_a_b ] ) immediately yield that overlap matrices @xmath69}$ ] and @xmath70}$ ] are the same .
the matrix element of @xmath118}=m^{[b]})$ ] is given by @xmath119}|\prod_{\la i , j \ra\in{\cal b}_a}(a_i b_j - b_i a_j )   \prod_{k \in \lambda_a } ( a_k)^{1/2+\al_k } ( b_k)^{1/2-\al_k } \nonumber \\ & \times \prod_{k \in \lambda_a } ( a\dag_k)^{1/2+\be_k } ( b\dag_k)^{1/2-\be_k } \prod_{\la i , j \ra\in{\cal b}_a}(a\dag_i b\dag_j - b\dag_i a\dag_j)|{\rm vac}^{[a]}\ra .
\label{integral0}\end{aligned}\ ] ] the rank of @xmath82 is @xmath120 , where we denote the number of elements in a set @xmath10 by @xmath121 .
the matrix element of @xmath82 is real , i.e. , @xmath122 , because the commutations of bosons do not produce any complex phases .
therefore , one can easily show the symmetric property @xmath123 .
we now rewrite the above expression in terms of the spherical angles .
this can be done by introducing spin coherent state representation .
let us first introduce spinor coordinates @xmath124 at each vertex ( site ) @xmath20 : @xmath125 where @xmath126 and @xmath127 .
then for a point @xmath128 on the unit sphere , the spin-@xmath21 coherent state is defined as follows : @xmath129}\ra.\ ] ] here we have fixed the overall phase ( u@xmath130 gauge degree of freedom ) since it has no physical content .
the coherent state is not orthogonal but complete and the resolution of identity is given by @xmath131 where @xmath132 denotes the simultaneous eigenstate of @xmath133 and @xmath134 , and @xmath135 is the @xmath136-dimensional identity matrix . by inserting this resolution of the identity , eq .
( [ integral0 ] ) can be recast in an integral form : @xmath137 where we have used the fact @xmath138 .
therefore , if we can obtain the matrix element of @xmath82 by some method such as monte carlo integration , the only thing to do is to diagonalize the matrix @xmath82 .
once the eigenvalues of @xmath82 , i.e. , @xmath139 ( @xmath140 ) , are obtained , the von neumann entropy can be calculated through eq .
( [ ee_in_terms_of_m ] ) .
let us now estimate the upper bound on the von neumann entropy .
suppose that all the eigenvalues of @xmath82 are the same .
this means that the subsystems @xmath65 and @xmath66 are _ maximally entangled_. then , the von neumann entropy is given by @xmath141 .
this value gives an upper bound on the von neumann entropy , i.e. , @xmath142 .
however , as we will see in the next section , numerical results indicate a deviation from this naive estimate and the von neumann entropy @xmath91 is strictly less than @xmath141 in 2d vbs .
it is interesting to note that there is a _ holographic _ interpretation of the overlap matrix @xmath82 .
the matrix @xmath82 can be regarded as an operator acting on the auxiliary one dimensional chain which is along the boundary of @xmath65 , i.e. , @xmath98 . from the fact that @xmath143 , one can rewrite @xmath82 as @xmath144 here
we have used the fact that @xmath145 . note that the pauli matrices @xmath146 ( @xmath147 ) act on the auxiliary space , @xmath148}\ra , |\downarrow^{[i]}\ra \}$ ] . then using the following relation : @xmath149 the overlap matrix can be expressed as a matrix - valued correlation function : @xmath150}|\prod_{i \in \lambda_a } \left ( s_i+\frac{1}{2}+{\vec s}_i \cdot { \vec \sigma}_i \right)|{\rm vbs}^{[a]}\ra , \label{overlap_and_cor}\ ] ] where @xmath151}\ra \equiv \prod_{\la i , j \ra \in { \cal b}_a}(a\dag_i b\dag_j - b\dag_i a\dag_j ) |{\rm vac}^{[a ] } \ra \label{half_vbs}\ ] ] note that the matrix rank of @xmath152 for @xmath93 is @xmath153 while that for @xmath154 is @xmath155
. one can interpret the overlap matrix as a hamiltonian acting on @xmath156-spin space .
the matrix element of this hamiltonian is governed by the boundary correlation in the half - vbs state given in eq .
( [ half_vbs ] ) .
the naive upper bound value of the von neumann entropy , @xmath141 , corresponds to the strongly disordered limit where any boundary correlation among @xmath10-spins vanishes .
however , as shown numerically in @xcite , the spin - spin correlation functions at small distances are nonzero in the vbs - type state even though they decay exponentially fast .
therefore , the leading term in eq .
( [ overlap_and_cor ] ) is given by a constant plus a term proportional to the spin-@xmath157 heisenberg hamiltonian whose coupling @xmath158 is given by the nearest neighbor @xmath10-spin correlation function .
the eigenvalues of this matrix can be different and hence the entanglement spectrum may not be flat .
therefore , it is plausible that the von neumann entropy is strictly less than @xmath159 even in the 2d infinite system . to support this conjecture ,
we perform numerical and analytical calculations for the vbs states on square and hexagonal lattices in the following sections .
in the previous section , we derived the integral formula for the overlap matrix ( eq .
( [ integral1 ] ) ) which helps us to reduce computational cost .
let us explain this point in more detail .
suppose that the graph @xmath49 has @xmath160 vertices and every vertex has the same spin @xmath10 .
then , the needed dimension for representing the vbs state on @xmath49 is @xmath161 without consideration of symmetry .
however , if we focus on the von neumann entropy of reflection symmetric vbs state , the needed dimension is greatly reduced from @xmath161 to @xmath120 , where @xmath162 denotes the number of sites on the boundary between @xmath65 and @xmath66 .
therefore , we can study the von neumann entropy of 2d vbs state on a relatively large system by combining monte carlo integration and exact diagonalization . after obtaining the overlap matrix and all its eigenvalues ,
the von neumann entropy can be calculated according to eq .
( [ ee_in_terms_of_m ] ) . in this section ,
we calculate von neumann entropies of vbs states on square and hexagonal lattices with open boundary conditions using monte carlo integration .
we first study the von neumann entropy of the vbs state on square lattices with open boundary conditions .
the reflection axis of the square lattice is taken to be the boundary between @xmath65 and @xmath66 as shown in fig .
[ fig : squarelattice ] .
@xmath163 in fig .
[ fig : squarelattice ] is the number of sites along @xmath164-axis while @xmath165 is that along @xmath166-axis . for square lattice , the number of boundary sites @xmath162 is the same as @xmath165 , i.e. , @xmath167 .     and @xmath168 denote the numbers of sites along @xmath164- and @xmath166-axes in @xmath65 , respectively . ]
we first study the @xmath163-dependence of the von neumann entropy .
figure [ graph : ee - square - nx ] shows the von neumann entropy per boundary site , i.e. , valence bond , for various @xmath165 as a function of @xmath163 .
the von neumann entropy per valence bond , @xmath169 , rapidly converges to certain values depending on @xmath165 .
the case of @xmath170 reduces to the linear aklt chain along @xmath164-axis . in this case , @xmath171 monotonically decreases with increasing @xmath163 and approaches to @xmath172 exponentially fast .
this correctly reproduces the exact result in @xcite .
the obtained results show that @xmath169 for @xmath173 and @xmath174 approach to @xmath175 and @xmath176 , respectively .
these values are consistent with the analytical results shown in the next section .     for the square vbs states with various @xmath165 .
the dotted lines indicate the exact results in @xmath177 limit . ]
next we consider the @xmath165-dependence of the von neumann entropy .
figure [ graph : ee - square - ny ] shows the von neumann entropy per boundary site as a function of @xmath165 .
we find that @xmath169 is well fitted by the following function : @xmath178 where @xmath179 , @xmath180 and @xmath181 are fitting parameters .
the constant term @xmath181 means the von neumann entropy per valence bond in the limit of @xmath182 .
note that @xmath179 , @xmath180 and @xmath181 depend on @xmath163 .
the blue curves in fig .
[ graph : ee - square - ny ] represent the fitting curves .
the obtained numerical data show good agreement with the power law behavior assumed in eq .
( [ eq : fitting ] ) .     for the square lattice vbs states with various @xmath163 .
blue curves represent the fitting curves ( eq .
( [ eq : fitting ] ) ) . ]
.the coefficients of the fitting function eq .
( [ eq : fitting ] ) for square lattices . [ cols="^,^,^,^",options="header " , ]      next we shall study the vbs states on the 3-leg and 4-leg hexagonal ladders . as an illustration , in fig .
[ fig : hexagonal_vbs_ladder2 ] , we show an example of 3-leg hexagonal ladders .
we first study the 3-leg ladders . in this case
, the @xmath183 matrix is defined similarly to eq .
( [ z_n_3leg ] ) : @xmath184 with @xmath185 and @xmath186 .
the recursion relation becomes slightly complicated compared to the square cases due to the wiggly nature of the hexagonal ladders .
it is composed of 2 steps as follows : @xmath187 where two matrices @xmath188 and @xmath189 are given by @xmath190 respectively , with @xmath191 . using the above recursion relation ,
one can obtain the coefficients @xmath192 . by diagonalizing the matrix in eq .
( [ z_n_3leg_h ] ) with these coefficients , we can calculate the von neumann entropy along the same lines as sec .
[ sec : 3-leg square ladder ] .
the obtained results are summarized in table [ tab : hexagonal ladders ] , which are compared with numerical results obtained in sec .
[ sec : mc hexagonal ] .
we now consider the large - block - size limit .
let @xmath193 be the perron - frobenius vector of @xmath194 .
this vector is numerically given by @xmath195 using the above values , we obtain the von neumann entropy in the large - block - size limit as @xmath196 and hence the von neumann entropy per valence bond is given by @xmath197 , which is strictly less than @xmath0 .    as is obvious , the argument so far can be straightforwardly generalized to @xmath198-leg ladders with @xmath199  @xcite .
in table [ tab : hexagonal ladders ] , we show the results for the 4-leg hexagonal ladders .
all the analytically obtained results indicate that the boundary correlation does not vanish even in the large - block - size limit and support the conjecture in sec .
[ sec : schmidt ] , i.e. , the eigenvalues of the reduced density matrix can be different . a numerical implementation of the transfer matrix technique for the @xmath198-leg ladders with @xmath200 will be discussed elsewhere  @xcite .
so far , we have studied square or hexagonal ladders in a horizontal geometry . in this section , we shall study the opposite limit , i.e. , the ladders in a vertical geometry .
let us first consider the square ladder with @xmath201 and the height @xmath165 shown in fig .
[ fig : vertical_ladder](a ) .
+ then the @xmath183 matrix is given by @xmath202 similarly to the horizontal ladder cases , we can use the integral formula eq .
( [ 3.3 ] ) and rewrite eq .
( [ n mat 1 ] ) as @xmath203 where @xmath204 .
here we have abbreviated @xmath205 .
the above formula can also be obtained using loops and strands shown in [ sec : loop model ] .
note that @xmath206 is @xmath207 identity matrix .
we now look for a recursion relation between @xmath208 and @xmath209 .
we first decompose @xmath209 into four matrices as @xmath210 where @xmath211 trivially acts on the @xmath165-th vector space @xmath212 while @xmath213 acts on it as @xmath214 .
this decomposition is graphically shown in fig .
[ fig : decomp ] .    .
]    using the graphical representation , we find the following set of recursion relations : @xmath215 \otimes \sigma^\alpha , \\ \fl~~   z(n_y+1 ) & = & z^0 ( n_y+1 ) + \sum^3_{\alpha=1 } n^\alpha(n_y+1),\end{aligned}\ ] ] where ` . '
denotes matrix multiplication while ` @xmath216 ' denotes tensor product , and ` @xmath217 ' is the @xmath207 identity matrix . solving the above recursions by _ mathematica _ , we obtain the overlap matrix @xmath209 exactly up to @xmath218
. then we numerically diagonalize @xmath209 and calculate the von neumann entropy as a function of @xmath165 .
figure [ fig : algebraic_nx1](a ) shows the result .
the obtained values and also fitting parameters in eq .
( [ eq : fitting ] ) show good agreement with the monte carlo results .
@xmath219{ee_l_s_b_o_nx1_exact.eps } & \hspace{5 mm } & \includegraphics[width=0.45\columnwidth]{ee_l_h_b_o_nx1_exact.eps } \\ ( { \rm a } )   & \hspace{5 mm } & ( { \rm b } ) \end{array}$ ]    let us next consider the hexagonal ladder with @xmath201 and the height @xmath165 shown in fig .
[ fig : vertical_ladder](b ) .
it is convenient to set @xmath220 since @xmath165 is always odd in this case . in this case , from eq .
( [ eq : def of z ] ) , the @xmath183 matrix is given by   z(m ; _ 1 , ... , _ m)= ( ^n_y_i=1 ) ^m_j=1(1+_2j-1_j ) ^n_y-1_k=1(1-_k _ k+1 ) . along the same lines as the square ladder case
, we obtain the set of recursion relations for the @xmath183 matrix : @xmath221 \otimes \sigma^\alpha , \\ \fl~~ z(m+1)=z^0(m+1)+\sum^3_{\alpha=1}z^\alpha ( m+1).\end{aligned}\ ] ] here we have abbreviated @xmath222 . again , solving the above recursions , we find the @xmath183 matrix up to @xmath223 ( @xmath224 ) . the numerically calculated von neumann entropy by this method is shown in fig . [
fig : algebraic_nx1](b ) .
the obtained values and the fitting parameters also show good agreement with the monte carlo results .
the obtained results for square and hexagonal ladders with @xmath201 strongly support that the assumed fitting function eq .
( [ eq : fitting ] ) correctly describes the behavior of large-@xmath165 systems .
in conclusion , we have studied entanglement properties of the vbs states on two - dimensional graphs with reflection symmetry . we have shown that the reflection symmetry permits us to develop an efficient method to study the reduced density matrix of the subsystem which is a mirror image of the other one .
we found the relation between the reduced density matrix and the overlap matrix whose element is given by the inner product between degenerate ground states of the block hamiltonian defined in the subsystem .
this relation enables us to conjecture that the von neumann entropy per valence bond on the boundary between the subsystems is strictly less than @xmath0 in 2d systems .
this conjecture means that the entanglement spectrum is not flat in this system and is in contrast to the case of 1d vbs where the reduced density matrix is proportional to the identity matrix in the large - block limit .
we have given a _
holographic _ interpretation of the reduced density matrix in terms of the spin chain along the boundary between the subsystems , which describes a hidden correlation among the degenerate ground states of the block hamiltonian . to confirm the conjecture ,
we have numerically studied the eigenvalues of the reduced density matrix for finite square and hexagonal graphs by combining monte carlo integration with exact diagonalization . from the eigenvalues of the reduced density matrix
, we have calculated the von neumann entropy and have found that the von neumann entropy per valence bond is well fitted by a constant term with a function decaying algebraically in the length of the boundary .
the constant is strictly less than @xmath0 and supports our conjecture .
we have also analytically and algebraically studied the quasi-1d cases where graphs are on ladders .
the analytical results are totally consistent with the numerical results and strongly support our conjecture .
although our conjecture is very plausible , we still lack a rigorous proof of it in the infinite 2d system .
therefore , it would of course be interesting to study the problem by a completely different approach based on inequalities and perhaps reflection positivity .
the authors are grateful to zhengcheng gu , akimasa miyake , masaki oshikawa , yutaka shikano , xiaoliang qi , and ying xu for their valuable comments and discussions .
the creative atmosphere of erwin schrdinger institute was also important for advancing the work .
the present work was supported in part by nsf grant dms 0905744 , by mext grant - in - aid for scientific research ( b ) ( 19340109 ) and for scientific research on priority areas `` novel states of matter induced by frustration '' ( 19052004 ) , and by next generation supercomputing project , nanoscience program , mext , japan .
is partly supported by the jsps postdoctral fellow for research abroad . s.t .
is partly supported by grant - in - aid for young scientists start - up ( 21840021 ) from the jsps .
the computation is executed on computers at the supercomputer center , issp .
in this appendix , we illustrate a graphical interpretation of the overlap matrix @xmath183 ( see eq .
( [ eq : integral=2leg_square ] ) ) in terms of loops and strands . as we have already mentioned in sec .
[ sec : square ladders ] , the terms containing odd number of @xmath225 s in @xmath226 vanish after integration over @xmath225 since the measure of integration is invariant under the local change of variable , @xmath227 .
this can be graphically explained as follows .
suppose that one vertex @xmath228 is connected to two vertices .
they are labeled 1 and 2 ( see fig . [
fig : klt](a ) ) .
then we consider the following integral : @xmath229.\end{aligned}\ ] ] from the lemma 3.3 in ref .
@xcite , we can show @xmath230 and hence we obtain @xmath231 where @xmath191 . this procedure
is graphically shown in fig .
[ fig : klt](a ) .
next , suppose that one vertex @xmath106 is surrounded by three vertices labeled by 1 , 2 , and 3
. then we encounter with the following integral : @xmath232 using eq .
( [ eq : integral_formula_1 ] ) , again , we obtain @xmath233 \nonumber \\ & = & 1 + q\ , \oh_1 \cdot \oh_2 + q\ , \oh_2 \cdot \oh_3 + q\ , \oh_3 \cdot \oh_1 . \label{loop2}\end{aligned}\ ] ] this procedure can also be graphically interpreted as shown in fig .
[ fig : klt](b ) .    ) .
the mark @xmath216 denotes the integration over the spherical angle @xmath225 .
the broken and colored ( thick ) lines have weight @xmath234 and @xmath191 , respectively .
( b ) graphical representation for eq . ( [ loop2 ] ) .
notations are the same as ( a ) . ]
therefore , as we have seen , the integrand containing odd number of @xmath225 s will vanish . in terms of the graphical representation ,
the matrix @xmath235 for the 2-leg square ladder can be written as @xmath236 where @xmath237 and @xmath238 denote the numbers of colored edges ( bonds ) and loops in the configuration @xmath239 , respectively .
the minus sign stems from the product like @xmath240 . in the first sum over @xmath241 ,
only the closed loops are allowed . on the other hand , in the second sum over @xmath242 , open ends from @xmath234 and @xmath243
are allowed ( see fig .
[ fig : config ] ) .     and open strand @xmath244 configurations for the length @xmath245 , and @xmath246 .
the weight of each configuration is assigned at the bottom right of each figure . ]    therefore , we can obtain the coefficients @xmath247 and @xmath248 in eq .
( [ eq : zmatrix_2leg_square ] ) by a completely graphical or combinatorial way . from fig .
[ fig : config ] , one can read @xmath249 , @xmath250 , @xmath251 and @xmath252 , @xmath253 , @xmath254 .
although we only show the case of 2-leg square ladders for simplicity , one can easily generalize this approach to any aklt ladders .
10 url # 1#1urlprefix[2][]#2 schrdinger e 1935 _ naturwissenschaften _ * * 23 * * 807