transition metal dichalcogenides ( tmds ) such as molybdenum disulfide ( mos@xmath0 ) have layered structures , where the atoms within a tmd monolayer form a covalently bonded planar network , and the interaction between these layers is a weak , van der waals interaction.@xcite a monolayer of mos@xmath0  consists of a layer of molybdenum atoms sandwiched between two layers of sulfur atoms .
monolayers can be exfoliated through micro - mechanical cleavage , similar to graphene or boron nitride.@xcite unlike graphene ( a metal ) , or boron nitride ( an insulator ) , mos@xmath0  is a semiconductor .
moreover , whereas bulk mos@xmath0  has an indirect band gap ( 1.2 ev ) , monolayer mos@xmath0  has a direct band gap ( @xmath21.8 - 1.9 ev ) , and shows a strong optical absorption and luminescence.@xcite at present mos@xmath0 , and tmds in general , are vehemently pursued as promising materials for applications in electronics and optoelectronics.@xcite    contacting mos@xmath0  to metal electrodes proves to be a problem ; it tends to produce unexpectedly high interface resistances , indicative of a high schottky barrier at the interface.@xcite a high barrier could be caused by strong interface bonding creating interface states that pin the fermi level,@xcite or by weak bonding creating a potential step due to pauli repulsion at the interface.@xcite the nature of the interaction at the mos@xmath0/metal interface is far from trivial .
on the one hand , one could argue that , as mos@xmath0  has no dangling bonds at its surface , its interaction with metal substrates should be weak and van - der - waals - like . on the other hand , many metal species form ( di)chalcogenide compounds,@xcite and when adsorbing mos@xmath0  onto a metal substrate , there could be a competition between the metal surface and the mo atoms for interacting with the sulfur atoms at the interface . in that case ,
not only the mos@xmath0/metal bonding would be a much stronger chemical bonding , but also the structure and electronic structure of the mos@xmath0  adsorbate could be significantly perturbed .    in this paper
we explore the adsorption of mos@xmath0  on a variety of metal substrates by means of first - principles density functional theory ( dft ) calculations , following up on work briefly reported in a short paper.@xcite previous dft studies have concentrated foremost on the schottky barrier formed at mos@xmath0/metal interfaces using the local density approximation ( lda).@xcite lda gives a reasonable description of the adsorption of graphene and _
h_-bn  on metal surfaces , but such results can not be generalized to other systems , as it is known that lda often leads to an unrealistic overbinding.@xcite other studies have used a generalized gradient approximation ( gga ) functional , such as pbe,@xcite which apparently works well for tmds adsorbed on metals,@xcite although it generally gives bad results for weakly bonded systems.@xcite    here we focus on the interface interaction and its implications for the structure and electronic structure of the mos@xmath0  adsorbate and the schottky barrier .
we choose a wide range of metal substrates : the ( 111 ) surfaces of al , ni , cu , pd , ag , pt and au , and the ( 0001 ) surfaces of mg , ti , and co , which are expected to have a wide range of interaction strength with the adsorbate .
as the interface interaction can vary from weak ( physisorption ) to strong ( chemisorption ) , it is a priori not clear which dft functional describes such bonding .
we test and compare results obtained with a van der waals functional , designed to describe weak , van der waals , interactions,@xcite to results obtained with gga and lda functionals , which are conventionally used to describe chemical bonding .
we assess the importance of van der waals interactions for the interface interaction , and evaluate its effect on the structure and electronic structure of the mos@xmath0  adsorbant .
we consider the situation where a mos@xmath0  layer is adsorbed as a whole on a metal substrate , making it more likely that the integrity of the mos@xmath0  layer is preserved in the adsorption process .
if the mos@xmath0/metal interaction is not too strong , and the mos@xmath0  and metal surface lattices are not matched , the interface structure is likely to be incommensurable . in a supercell calculation one
is forced to approximate such a structure by a commensurable one .
previous calculations have used small supercells , where in some cases appreciable artificial strain is generated because of the mismatch between the mos@xmath0  and the metal surface lattices.@xcite we apply a strategy for choosing supercells such that the artificial strain is minimal , and test the influence of strain on the electronic properties of the interface .
this paper is organized as follows .
section  [ sec : calculations ] describes the dft calculations , comparing different functionals in sec .
[ sec : functionals ] and discussing the effect of lattice mismatch in sec .
[ sec : mismatch ] .
results are discussed in sec .
[ sec : results ] , with the metal / mos@xmath0  interaction in sec .
[ sec : bonding ] and its effects on the interface potential step and the schottky barrier in sec .
[ sec : sb ] .
strong chemisorption is discussed in more detail in sec .
[ sec : ti ] , and a summary and the conclusions are presented in sec .
[ sec : conclusions ] .
we calculate ground - state energies and optimize geometries at the density functional theory ( dft ) level , using projector - augmented waves ( paws ) as implemented in the vasp code.@xcite the plane - wave kinetic - energy cutoff is set at 400 ev .
the surface brillouin zone is integrated with the methfessel - paxton technique using a smearing parameter of 0.05 ev,@xcite and a @xmath3-point sampling grid with a spacing of 0.01 @xmath4 .
the mos@xmath0/metal interface is modeled as a slab of 4 - 6 layers of metal atoms with one or two layers of mos@xmath0  adsorbed on one side and a vacuum region of @xmath212  .
the in - plane supercell is chosen such as to minimize the mismatch between the mos@xmath0  and metal lattices , which is discussed in more detail in sec .
[ sec : mismatch ] .
a dipole correction is applied to avoid spurious interactions between periodic images of the slab.@xcite we allow the positions of the atoms to relax until the force on each atom is smaller than 0.01 ev@xmath4 , except for the bottom layer of metal atoms , whose positions are kept fixed .
the electronic self - consistency criterion is set to 10@xmath5 ev .
it is well known that commonly used dft exchange - correlation functionals , based upon lda@xcite or gga,@xcite give decent descriptions of covalent and ionic bonding , but they may fail for weakly bonded systems , as such functionals do not contain a description of van der waals interactions . for example , gga functionals such as pw91 or pbe,@xcite do not capture the bonding between _
h_-bn  or graphene layers , nor that between _
h_-bn  or graphene and transition metal(111 ) surfaces.@xcite a priori
we do nt know how important van der waals interactions are in the bonding between mos@xmath0  and a metal surface . in sec .
[ sec : functionals ] we compare results obtained using a van der waals density functional ( vdw - df),@xcite with results obtained with gga and lda functionals .    structure with corresponding plane - averaged electrostatic potential @xmath6 .
@xmath7 is defined as the difference between the work function on the metal side and on the mos@xmath0  side . ]
one way of visualizing bonding at a mos@xmath0/metal interface is by the electron density difference @xmath8 where @xmath9 , @xmath10 , and @xmath11 are the electron densities of mos@xmath0  adsorbed on the metal , of the metal surface and of the free standing mos@xmath0 , respectively .
the system as a whole is neutral , and @xmath12 is localized around the metal / mos@xmath0  interface , i.e. @xmath13 for @xmath14 sufficiently far from the interface .
solving the poisson equation with @xmath15 as source then gives a potential step across the interface @xmath16 here @xmath17 is the direction normal to the interface , @xmath18 is the interface area , and @xmath19 is the difference between the asymptotic values of the potential left and right of the interface .
figure  [ fig : deltav ] also illustrates an alternative expression for the potential step @xmath20 where @xmath21 , @xmath22 are the work functions of the clean metal surface , and of the metal surface covered by mos@xmath0 , respectively .
a practical way of obtaining work functions from dft calculations is to track the plane - averaged electrostatic ( hartree ) potential @xmath6 into the vacuum , see fig .
[ fig : deltav ] , where typically the asymptotic value is reached with a few   from the surface . in converged calculations the expressions of eqs .
[ eq : sb2 ] and [ eq : sb3 ] give results that are with a few mev of one another .
the interface potential step plays an important role in determining the schottky barrier , i.e. , the potential barrier for charge carrier transport from a metal into a semiconductor .
the schottky barrier for electrons can be written as @xmath23 with @xmath24 the fermi level and @xmath25 the electron affinity of mos@xmath0 .
there are several ways to extract the schottky barrier height from mos@xmath0/metal slab calculations .
one could determine @xmath26 by measuring @xmath27 in the band structure or in the density of states of the mos@xmath0/metal slab , as in refs .   and .
in order to identify the mos@xmath0  related states , one needs to calculate the amplitudes of the projections of the wave functions of the slab on the mos@xmath0  layer .
there is always some arbitrariness involved in such a projection if the adsorbate and the substrate are in close connection . in addition ,
identification of states belonging to the adsorbate is possible only if its electronic structure is not significantly perturbed in the adsorption process , which is only the case if the adsorbate is ( weakly ) physisorbed on the substrate.@xcite in practice we find that this procedure for obtaining the schottky barrier height at mos@xmath0/metal interfaces is not sufficiently accurate when applied to the ( projected ) density of states , and of practical use only when applied to the ( projected ) band structure of a small supercell .
in contrast , starting from the right - hand side of eq .
[ eq : sb ] , @xmath21 , @xmath25 and @xmath7 are easily obtained in separate calculations on the clean metal surface , the free - standing mos@xmath0  layer , and the mos@xmath0/metal slab , respectively .
of course , if the mos@xmath0  electronic structure is very strongly perturbed by adsorption , one has to reconsider the definition of the schottky barrier , see sec .
[ sec : ti ] . in the following the potential step @xmath19
is used to characterize the mos@xmath0/metal interface , along with the binding energy and the structure .
an alternative way of locating the conduction band edge @xmath25 in a mos@xmath0/metal slab calculation , without having to resort to wave function projections , is by aligning the core levels of the mo or s atoms in the slab with the corresponding core levels in free - standing mos@xmath0 .
it allows us to compare the densities of states of free - standing and adsorbed mos@xmath0 , see sec .
[ sec : sb ] . in principle
, this procedure also allows for calculating the schottky barrier height . as the method described in the previous paragraph , this only makes sense if the mos@xmath0  electronic structure is not perturbed too strongly by the adsorption .
materials such as graphite , _ h_-bn , and mos@xmath0have a layered structure , where the atoms within one layer form strong covalent bonds , but the interaction between the layers consists of weak , van der waals , forces .
common gga functionals , such as pbe,@xcite lack a description of van der waals interactions , which results in a severe underestimation of the interlayer binding energy in graphite and _ h_-bn , and an overestimation of the interlayer bonding distance.@xcite similar problems are encountered when graphene or _
adsorbed on a metal substrate.@xcite the lda functional also lacks a description of van der waals interactions , but it , somewhat fortuitously , gives reasonable binding energies and geometries for graphite , _
h_-bn , and for the adsorption of these materials on metals.@xcite in general however , the lda functional tends to overestimate binding energies , which is regularly accompanied by an underestimation of the bonding distance .
[ sec : functionals ]    ccccccccccccc    & & & & +   + & @xmath28( ) & @xmath7(ev ) & @xmath29(ev ) & @xmath28( ) & @xmath7(ev ) & @xmath29(ev ) & @xmath28( ) & @xmath7(ev ) & @xmath29(ev ) & @xmath28( ) & @xmath7(ev ) & @xmath29(ev ) + lda & 2.6 & 0.54 & @xmath300.27 & 2.5 & 0.10 & @xmath300.33 & 2.2 & 0.50 & @xmath300.69 & 2.4 & 0.85 & @xmath300.43 + pbe & 3.3 & 0.38 & @xmath300.02 & 2.8 & 0.10 & @xmath300.08 & 2.3 & 0.34 & @xmath300.25 & 2.6 & 0.66 & @xmath300.11 + vdw - df & 2.9 & 0.41 & @xmath300.33 & 2.8 & 0.11 & @xmath300.35 & 2.4 & 0.30 & @xmath300.54 & 2.6 & 0.71 & @xmath300.43 +    [ tab : functionals ]    many of these problems are mitigated when using vdw - dfs,@xcite which , for instance , describe the bonding in graphite very well.@xcite the exchange - correlation energy in a vdw - df takes the form @xmath31 where @xmath32 , @xmath33 and @xmath34 describe the exchange part , and the local and nonlocal electron - electron correlations , respectively .
for @xmath34 we use the vdw kernel developed by dion _
_ @xcite and for @xmath33 the correlation part of the lda functional . for the exchange part
@xmath32 we use the optb88 functional.@xcite the opt88-vdw - df has given good results for binding energies and geometries of graphite , _
h_-bn , and the adsorption of these materials on metals.@xcite    in the following we test the gga / pbe , lda , and opt88-vdw - df functionals for the adsorption of mos@xmath0  on metals . as test cases we use the 4d and 5d metals ag , au , pd , and pt .
we place a mos@xmath0  monolayer on top of the ( 111 ) surface of these metals , choosing a @xmath35 in - plane mos@xmath0  unit cell on top of a @xmath36 ( 111 ) surface cell .
the in - plane mos@xmath0  lattice parameters are kept at their optimized values for a free - standing layer , and the in - plane metal lattice parameter is adapted accordingly .
the size of the adaption is maximal for au , where it results in a compression of the in - plane au lattice by 4.2% .
the effects of this artificial strain are discussed in the next section .
figure  [ fig : binding](a ) shows the binding curves of mos@xmath0on au and pt(111 ) for the three functionals .
the binding energy is defined as the total energy per mos@xmath0formula unit of the metal / mos@xmath0slab minus the total energies of the clean metal slab and the free - standing mos@xmath0layer , as a function of the distance @xmath37 between the top layer of metal atoms and the bottom layer of sulfur atoms . for mos@xmath0  on au(111 ) , pbe gives virtually no bonding , and opt88-vdw - df gives a sizable binding energy .
the opt88-vdw - df result suggests that mos@xmath0  is physisorbed on au(111 ) , with van der waals interactions playing the decisive role in the bonding .
pbe does not capture this at all .
lda gives a equilibrium binding distance that is 0.3   smaller , and an equilibrium binding energy that is 32% larger .    for mos@xmath0  on pt(111 ) all three functionals give equilibrium bonding distances that are shorter than for mos@xmath0  on au(111 ) , and a bonding that is stronger , which suggests that mos@xmath0  may be weakly chemisorbed on pt(111 ) . pbe and opt88-vdw - df give a similar equilibrium distance , although pbe captures only 26% of the binding energy , indicating that van der waals interactions still play an important role here .
lda gives a similar binding energy as opt88-vdw - df , but an equilibrium binding distance that is 0.2   smaller .
table [ tab : functionals ] shows the equilibrium binding distances and energies obtained with the three functionals for mos@xmath0on au , ag , pd , and pt(111 ) . treating the results for opt88-vdw - df as a benchmark
, pbe is seen to severely underestimate binding energies , whereas lda gives quite reasonable binding energies .
lda however gives binding distances that are up to 0.3 shorter than those obtained with opt88-vdw - df , in particular for cases where the bonding is weak , such as au and ag .
in contrast , pbe gives binding distances that are similar to those obtained with opt88-vdw - df , except for au , where pbe essentially fails to give any significant bonding .
potential steps @xmath7 as a function of the distance @xmath37 between the top layer of metal atoms and the bottom layer of sulfur atoms , calculated according to eq .
[ eq : sb3 ] , are shown in fig .  [
fig : binding](b ) for au and pt . the curves for the pbe and the opt88-vdw - df functionals are within 0.05 ev of one another in the range @xmath38 - 3  , whereas lda gives a potential step that is 0.10 - 0.15 ev higher . in view of the considerable differences in the binding curves for these three functionals ,
the differences in the potential steps are remarkably small .
this is true for all metal substrates listed in table  [ tab : functionals ] .    in ref .
the main contribution to the potential step in the adsorption of _
h_-bn  on metal substrates was attributed to pauli repulsion .
this can be modeled by an electron density that is obtained by anti - symmetrizing the product of the metal and the adsorbate wave functions . as long
as these wave functions do not strongly depend on the functional , the electron density and the potential step are also relatively insensitive to the functional used .
this is unlike the total energy , which for a given electron density is very dependent on the functional .
for the potential step to be accurate it is however important to obtain the proper equilibrium binding distance.@xcite      the absolute values of the binding energies given in table  [ tab : functionals ] are much smaller than what one expects to find for true chemical bonding .
the differences between the values obtained with pbe and opt88-vdw - df indicate that van der waals interactions play a significant role in the bonding .
with such a weak metal / adsorbate bonding it is unlikely that the metal substrate can enforce its lattice periodicity onto the mos@xmath0  overlayer .
therefore , a metal / mos@xmath0  interface very likely becomes incommensurable if the metal / mos@xmath0  lattice mismatch is substantial .
in electronic structure calculations one is forced to use commensurable structures to model incommensurable systems .
obviously care must be taken to ensure that the artificial strain introduced this way , does not alter the electronic structure in an unrealistic way .
based upon previous experience , we expect that modifying the in - plane lattice constant of a close - packed metal surface by a few percent only affects its electronic properties mildly.@xcite in contrast , changing the lattice parameter of mos@xmath0  by just one percent already alters the band gap by @xmath20.1 ev , and changes it from direct to indirect .
a larger change in the lattice parameter has an even more dramatic effect .
applying a tensile strain of @xmath39% to mos@xmath0  reduces the band gap by @xmath40 ev.@xcite    monolayer in a @xmath41@xmath42@xmath41 cell , where the direct band gap appears at @xmath43 ; ( b ) band structure of mos@xmath0/au(111 ) with the in - plane au lattice compressed by 4.2% to match the mos@xmath0  lattice ; the blue color indicates the weight of a projection of the wave functions on the mos@xmath0  sites ; ( c ) as ( b ) but with the mos@xmath0  lattice stretched by 4.2% to match the au(111 ) lattice.,scaledwidth=50.0% ]    as an example , the pbe optimized in - plane lattice parameters of mos@xmath0  and au(111 ) are 3.19 and 2.88  .
placing a @xmath44 mos@xmath0  cell on top of 2@xmath422 au(111 ) surface cell then leads to a lattice mismatch of 4.2% .
figure  [ ref : strain](b ) shows the electronic band structure of mos@xmath0/au(111 ) where the in - plane au(111 ) is compressed by 4.2% to match the lattice parameter of mos@xmath0 . as the interaction between
mos@xmath0  and the au surface is relatively small , it is not surprising to see that the band structure of adsorbed mos@xmath0  resembles that of free - standing mos@xmath0 , shown in fig .
[ ref : strain](a ) .
note that in the @xmath45 mos@xmath0  cell the bands are folded such that the direct band gap appears at the @xmath43 point .
the work function of clean au(111 ) is changed by only 0.08 ev by the 4.2% compression of its lattice .    for comparison ,
[ ref : strain](c ) shows the band structure of mos@xmath0/au(111 ) when mos@xmath0  is stretched by 4.2% to match the au(111 ) lattice . clearly the band structure of mos@xmath0  is now changed significantly .
it no longer shows a direct band gap at @xmath43 , but an indirect band gap , and the size of the band gap is reduced to @xmath21 ev , which is consistent with previous studies.@xcite the schottky barrier for electrons ( the energy difference between the bottom of the conduction band and the fermi level ) , which is a sizable 0.7 ev in fig .
[ ref : strain](b ) , is reduced to zero in fig .
[ ref : strain](c ) as in ref .  .
the latter is clearly unphysical : one would not expect a high work - function metal such as au to form a barrierless contact for electrons . indeed experimentally au is found to form a substantial schottky barrier with mos@xmath0.@xcite    in the following we base the in - plane lattice constant of the mos@xmath0/metal slab on the optimized values of free - standing mos@xmath0 , which are 3.13 , 3.18 , and 3.19 for the lda , optb88-vdw - df , and pbe functionals , respectively . experimentally reported bulk mos@xmath0  lattice constants are in the range 3.13 - 3.16 ,@xcite suggesting that the lda result may be more accurate and both pbe and the vdw functional are overestimating the lattice constant somewhat .
in making a commensurable structure we adapt the metal to the mos@xmath0  lattice . to minimize the artificial strain that is introduced by this adaptation
, we construct in - plane supercells following the procedure of ref .  .
we denote a basis vector of a mos@xmath0  supercell by @xmath46 , with @xmath47 the basis vectors of the primitive cell , and @xmath48 integers .
similarly , @xmath49 is a basis vector of a metal surface supercell , with @xmath50 the basis vectors of the primitive cell , and @xmath51 integers .
we search for a set of values for @xmath52,@xmath53,@xmath54 , and @xmath55 , such that the difference between the mos@xmath0  and the metal supercell basis vectors is less than a margin @xmath56 , @xmath57 we then rotate the mos@xmath0  lattice by an angle @xmath58 such , that the directions of the @xmath59 and @xmath60 vectors coincide . because of the symmetry of the lattice the second basis vector of the supercell
is easily obtained by a 120@xmath61 rotation , @xmath62 .
the commonly used surface science notation of this supercell is a @xmath63 mos@xmath0  lattice on top of a @xmath64 metal lattice , where @xmath65 and @xmath66 .    .
in - plane supercell defined by the mos@xmath0  lattice vector @xmath67 , where @xmath46 and the metal lattice vector @xmath49 .
@xmath56 gives the mismatch between the mos@xmath0  and metal lattices , eq .
[ eq : sb6 ] ( pbe values ) . [ cols="<,^,^,^,^",options="header " , ]     [ tab : structures ]    the displacements are quite large for the mos@xmath0/ti(0001 ) structure , indicating that there is a significant distortion of the geometries of both the mos@xmath0  adsorbate and the ti surface , which strongly suggests that mos@xmath0  is chemisorbed on ti . at the opposite end of the scale we find mos@xmath0/au(111 ) , where the atomic displacements are small , indicating that here we are in the physisorption regime .
the behavior of the other metal substrates is in between these two extreme cases but more to the physisorption side .
the 3d transition metals co , ni , cu and the low work function simple metal mg show somewhat larger distortions than the 4d and 5d metals pd , ag , pt and the simple metal al .      table  [ tab : deltav ] gives the interface potential steps @xmath19 created by the adsorption of mos@xmath0  on a metal substrate .
this potential step strongly influences the schottky barrier at metal / mos@xmath0  contacts , see eq .
[ eq : sb ] , and as such it plays an important role in the physics of mos@xmath0  semiconductor devices .
the potential steps can be divided into two groups , i.e. , positive @xmath19 for metals with a high work function , and negative @xmath19 for low work function metals .
a positive @xmath19 means that adsorption of mos@xmath0  effectively lowers the work function of the substrate .
the mos@xmath0  layer has no intrinsic dipole moment perpendicular to the layer that could create such a potential step .
so the work function lowering is a purely electronic effect that results from the displacement of surface electron density into the metal by physisorption of the adsorbate .
lccccccc & & +   + & @xmath68 & @xmath69 & @xmath70 & @xmath71 & @xmath69 & @xmath70 + mg & @xmath72 & @xmath73 & @xmath74 & 3.96 & @xmath75 & @xmath76 + al & @xmath77 & @xmath78 & @xmath79 & 4.20 & @xmath80 & @xmath81 + ag & @xmath82 & @xmath83 & @xmath84 & 4.82 & @xmath85 & @xmath86 + ti & @xmath87 & @xmath88 & 0.53@xmath89 & 4.80 & @xmath90 & 0.54@xmath89 + cu & @xmath91 & @xmath92 & @xmath93 & 5.10 & 0.39 & @xmath86 + au & @xmath94 & @xmath95 & @xmath96 & 5.58 & @xmath97 & @xmath98 + pd & @xmath99 & @xmath92 & @xmath100 & 5.48 & @xmath101 & @xmath102 + pt & @xmath103 & @xmath104 & @xmath105 & 5.96 & @xmath96 & @xmath106 + co & @xmath107 & @xmath108 & @xmath109 & 5.42 & @xmath110 & @xmath111 + ni & @xmath112 & @xmath113 & @xmath114 & 5.40 & @xmath115 & @xmath116 +    @xmath89 see sec .  [
sec : ti ] .
[ tab : deltav ]    this effect is known as the push - back effect or the pillow effect , which is a general phenomenon observed in the physisorption of closed - shell atoms , molecules , and layers on metal substrates . in ref .
we have developed a quantitative model for this effect , based upon an anti - symmetrization of the product of the metal and adsorbate wave functions .
when an adsorbate approaches a metal surface , the wave functions of the two systems overlap .
pauli exchange repulsion between these states leads to a spatial redistribution of the electron density , in particular to a decrease of the density in the overlap region .
since the metal wave functions are usually more extended and more easily deformable than those of the adsorbate , the net result of this redistribution is that electrons are pushed back into the metal , which effectively lowers the work function .    in the adsorption of graphene and _ h_-bnon high work function metals we found potentials steps of up to 1 - 2 ev .
the potential steps for mos@xmath0  adsorbed on the same metals are generally smaller , and more typically around 0.3 - 0.4 ev .
the wave functions of first - row elements ( b , c , n ) are compact and not easily deformable , as compared to the wave functions of the metal substrate .
the effect of pauli repulsion in the metal / adsorbate overlap region is then very asymmetric .
it is foremost the metal electron density that is deformed , i.e. pushed back , which gives a large work function lowering . if the adsorbate contains heavier elements , such as mos@xmath0 , the effect of pauli repulsion is more symmetric , i.e. , both the metal and the adsorbate electrons are pushed out of the overlap region in a more symmetric way .
this gives a smaller effect on the work function .
note that if the effect of pauli repulsion would be completely symmetric , the work function would be unchanged .
/al(111 ) slab ; ( red ) the sum of the densities of states of free - standing mos@xmath0  and of the al(111 ) slab .
the densities of states are aligned by aligning the mo 4s core levels , and the al 2p core levels . ]
low work function metals experience an increase of the work function upon adsorption of mos@xmath0 , i.e. a negative @xmath19 , which indicates a net transfer of electrons from the metal to the mos@xmath0  adsorbate .
as mos@xmath0  is a semiconductor it can only receive electrons in its conduction band . therefore , for low work function metal substrates one expects the fermi level to be in the conduction band of mos@xmath0 .
analysis of the electronic structure of the mos@xmath0/metal slab however shows that this is not the case .
the interaction between mos@xmath0  and the metal at the interface leads to interface states with energies in the mos@xmath0  band gap .
that seems obvious if mos@xmath0  is chemisorbed onto the substrate , as in the case of mos@xmath0/ti(0001 ) , which we will discuss in the next section .     of the mos@xmath0/al(111 ) interface at the equilibrium distance @xmath117  , and ( b ) at a distance @xmath118  ; ( c ,
d ) the corresponding band structures ; the blue color indicates the weight of a projection of the wave functions on the mos@xmath0  sites . ]
somewhat surprisingly , a significant density of interface states also forms if the interaction between mos@xmath0  and the metal substrate is relatively weak .
for example , as discussed in the previous section , the interaction between mos@xmath0  and al(111 ) can be classified as physisorption , see fig .
[ fig : physchem ] .
nevertheless , states with energies inside the mos@xmath0  band gap are formed at the mos@xmath0/al(111 ) interface , as is immediately obvious when comparing the density of states of the interface with that of free - standing mos@xmath0 , see fig .
[ fig : dos_al ] .
the density of these interface states is not extremely high , yet sufficiently high to pin the fermi level below the mos@xmath0  conduction band , as demonstrated by fig .
[ fig : dos_al ] .
the density of interface states increases with increasing mos@xmath0/metal interaction , but even for physisorption it seems sufficiently high to prevent the fermi level from reaching the mos@xmath0  conduction band .
one can prove that these interface states are indeed responsible for pinning the fermi level by artificially enlarging the distance between the mos@xmath0  layer and the metal surface .
this breaks the direct mos@xmath0/metal contact that is responsible for the formation of interface states . in the absence of interface states ,
the fermi level is at the bottom of the mos@xmath0  conduction band , see fig .
[ fig : alfinal](d ) , which is what one would expect if the work function of the metal @xmath21 is smaller than the electron affinity of mos@xmath0  @xmath119 , see eq .
[ eq : sb ] .
a transfer of electrons between the metal and the mos@xmath0  overlayer then yields a charge distribution that can be associated with a simple interface dipole , see fig .
[ fig : alfinal](b ) .
in contrast , if interface states are formed , the band structure of adsorbed mos@xmath0  is perturbed , see fig .
[ fig : alfinal](c ) , and the pattern of the charge distribution at the interface is much more complicated , as shown in fig .
[ fig : alfinal](a ) . in that case
, the interface states pin the fermi level below the bottom of the mos@xmath0  conduction band , see fig .
[ fig : dos_al ] .    schottky barrier heights ( sbhs ) for electrons , calculated according to eq .
[ eq : sb ] , are also listed in table  [ tab : deltav ] .
the functional causes some uncertainty , as the work functions of the clean metal surfaces obtained with the opt88-vdw - df functional tend to be somewhat higher than those obtained with the pbe functional .
lda in general gives even higher work functions , so opt88-vdw - df gives work functions that are in between those of pbe and lda.@xcite note that the interface potential steps @xmath19 do not depend strongly on the functional . as the opt88-vdw - df functional also gives a larger electron affinity for mos@xmath0 , the schottky barrier @xmath26 according to eq .
[ eq : sb ] , also does not depend strongly on the functional .
we see that the sbh decreases with decreasing metal work function but does not go to zero .
instead it goes through a minimum for cu and ag , and then increases again for low work function metals like al and mg .
as discussed above , this phenomenon is caused by interface states .
the only way to get rid of such states is to break the direct interaction between mos@xmath0  and the metal substrate .
.   discusses a practical way of doing this by inserting an atomic layer between the metal surface and the mos@xmath0  layer .
if this intermediate layer is purely van der waals - bonded to mos@xmath0 , no gap states are formed at its interface with mos@xmath0 .
in addition the intermediate layer should be transparent to electrons , such that the interface resistance is not dramatically increased .
a monolayer of _
h_-bn  or graphene satisfy these criteria.@xcite      adsorbed on ti(0001 ) .
]    as discussed in sec .
[ sec : bonding ] , mos@xmath0  is chemisorbed in ti(0001 ) .
the binding energy and the equilibrium bonding distance of mos@xmath0  on ti(0001 ) does not seem to be qualitatively different from those for mos@xmath0  on metal substrates such as co(0001 ) or ni(111 ) , see table  [ tab : energies ] and fig .
[ fig : physchem ] .
however , the structural deformation of the mos@xmath0  layer adsorbed on ti(0001 ) is much larger than that of mos@xmath0  on other metals , see table  [ tab : structures ] .
this deformation is illustrated in fig .
[ fig : ti_structure ] . atoms of the bottom sulfur layer make a bond with ti atoms of the top layer of the substrate , where several of these metal atoms are pulled up from the substrate .
the mos@xmath0  and the ti(0001 ) lattices do not fit very well ; one needs a @xmath120 mos@xmath0  supercell on top of a @xmath121 ti(0001 ) supercell to get a mismatch below 1% , see table  [ tab : supercell ] .
the result therefore is a mos@xmath0/ti(0001 ) interface that contains a substantial local strain , which explains why the binding energy is not very large , despite the bonding being chemisorption .
the potential step @xmath19 at the mos@xmath0/ti(0001 ) interface is negative , in contrast to the potential step at the mos@xmath0/ag(111 ) interface for instance , which is positive , despite the fact that the work functions of ti and ag are very similar , see table  [ tab : deltav ] .
we argued that physisorption should lead to a positive potential step because of the pauli repulsion effect , and indeed mos@xmath0  is physisorbed on ag(111 ) .
chemisorption , as in the case of mos@xmath0  on ti(0001 ) , leads to a more drastic reorganization of the charge distribution at the interface , because of the formation of new chemical bonds . upon the formation of these bonds there
is apparently a net displacement of electronic density towards the sulfur atoms , which is not unreasonable as sulfur is more electronegative than ti .
this displacement results in an increase of the work function , i.e. , a negative @xmath19 .
one expects that chemisorption also leads to a strong perturbation of the electronic structure of the adsorbate .
figure  [ fig : dos_ti ] shows the density of states ( dos ) of a mos@xmath0  bilayer projected on the individual mos@xmath0  layers .
the dos of the first ( chemisorbed ) layer is indeed strongly perturbed as compared to the dos of a free - standing mos@xmath0  layer .
the mos@xmath0  wave functions strongly hybridize with those of the underlying ti substrate , and the resulting hybridized states give a non - zero dos for energies all through the mos@xmath0  band gap .
it is sometimes argued that such interface states promote having a good ( ohmic ) mos@xmath0/metal contact.@xcite one could however also argue that chemisorption is harmful to obtaining a good contact , because it damages the integrity of the mos@xmath0  layer . in ref .
it is found that mos@xmath0/ag gives a better contact than mos@xmath0/ti , due to a much smoother interface in the former case , suggesting to prefer physisorption over chemisorption .
layer of mos@xmath0  adsorbed on ti(0001 ) and on the second mos@xmath0  layer , respectively . ]
it is not possible to define a schottky barrier for a single mos@xmath0  layer adsorbed on ti(0001 ) .
chemisorption affects the electronic structure of mos@xmath0  to such an extend that its semiconducting character is lost .
it is however possible to define a schottky barrier for a second mos@xmath0  layer that is adsorbed on the first layer . the first and second mos@xmath0  layer are bonded by a van der waals interaction , which does not perturb the electronic structure of the second layer significantly . indeed
the dos of the second mos@xmath0  layer is quite similar to that of a free - standing mos@xmath0  layer with a clear band gap , see fig .
[ fig : dos_ti ] .
the height of the schottky barrier to the second layer is a sizable 0.54 ev .
this mainly results from the fact that the first adsorbed mos@xmath0  layer effectively increases the work function of the ti substrate , see table  [ tab : deltav ] .
the size of the schottky barrier indicates that it is not possible to make an ohmic contact to pristine ( undoped ) mos@xmath0  with ti.@xcite
in this paper , we explore the adsorption of mos@xmath0  on a range of metal substrates by means of first - principles dft calculations .
the ( 111 ) surfaces of al , ni , cu , pd , ag , pt and au , and the ( 0001 ) surfaces of mg , co and ti cover a range of metals with different interaction strengths , allowing for a systematic study of the metal - mos@xmath0  interface .
we compare the results obtained with different dft functionals . in many cases
the gga / pbe density functional only captures a small part of the binding energy of mos@xmath0  on a metal substrate , as compared to the opt88-vdw - df van der waals density functional , which indicates the importance of van der waals interactions in the interface bonding . nevertheless , the equilibrium binding distances obtained with both functionals are generally very similar , and so are the interface potential steps and schottky barrier heights .
exceptions are cases for which van der waals interactions essentially describe the whole bonding , where pbe fails completely .
lda tends to overbind , leading to shorter binding distances and larger interface potential steps .
the interface structure that results from adsorbing an mos@xmath0  layer on a metal surface will be incommensurable in most cases , as the two lattices have a mismatch .
we investigate the effects of the artificial strain introduced by approximating the lattice using a commensurable supercell .
we conclude that these effects are moderate provided the mos@xmath0  lattice parameter is kept at its optimized value , and the metal lattice is strained .
large lattice mismatches should however be avoided , and straining the mos@xmath0  lattice can lead to very unphysical results.@xcite    of the metal substrates studied ti is the one on which mos@xmath0  is clearly chemisorbed .
adsorption of mos@xmath0  on ti(0001 ) is accompanied by a clear structural deformation of the ti surface and of the mos@xmath0  overlayer , due to the formation of bonds between the surface ti atoms and the sulfur atoms at the interface .
formation of these interface bonds significantly alters the electronic structure of the mos@xmath0  adsorbate .
in particular , the interface states fill up the band gap of mos@xmath0 , which makes defining a schottky barrier for this layer meaningless . however , for a second , unperturbed , adsorbed mos@xmath0  layer a schottky barrier of 0.54 ev can be extracted .
mos@xmath0  is physisorbed on au(111 ) , where the bonding is almost completely due to van der waals interactions , and the structure and electronic structure of mos@xmath0  are hardly perturbed by the adsorption .
the properties of mos@xmath0  adsorbed on other metal substrates fall in the range between the two extreme cases ( ti and au ) , without the possibility of drawing a clear dividing line , as has been done for the adsorption of graphene or _
h_-bn  on metal substrates.@xcite    experiments have focused foremost on schottky barrier heights .
transport measurements on multilayer mos@xmath0  devices generally yield small numbers for the schottky barrier heights , i.e. 0.03 - 0.2 ev , for different metals,@xcite whereas photoemission , photoconduction , and scanning tunneling spectroscopy give higher values 0.2 - 0.9 ev.@xcite it has been suggested that the mos@xmath0  samples used in devices is quite defective and inhomogeneous , such that the position of the fermi level does not reflect an intrinsic property of mos@xmath0  or of the mos@xmath0/metal contact,@xcite which obstructs a comparison to calculated results .    our results for the 4d and 5d metals ag , au , pd , and
pt agree qualitatively with those reported in previous calculations,@xcite provided the mos@xmath0  lattice is not stretched.@xcite quantitatively , the reported schottky barrier heights for these metals are @xmath122 ev larger than our results .
these calculations were based upon the lda functional , which tends to overbind , and to overestimate the metal work functions.@xcite compressing the metal lattice , which is sometimes required to accommodate a lattice mismatch in a small supercell , does not help either , as that gives an even higher work function.@xcite the same is likely true for simple metals such as al and in.@xcite the strong interaction we find for ti is also found in lda calculations.@xcite in those calculations the lattice mismatch used was large , however , which can alter the interface interactions .
the overall picture emerging from these calculations is that mos@xmath0  interacts strongly with the early transition metals , where it is clearly chemisorbed .
the interaction with the late transition metals is much weaker , where the 3d metals interact stronger than the 4d and 5d metals .
interacts rather weakly with the simple metals , but the interaction increases for very low work function metals . in all but the strongly chemisorbed case ,
van der waals forces play an important role in the interface interactions .    in case
the interface interaction is weak ( physisorption ) the interface potential step can be understood as resulting from pauli repulsion , which effectively decreases the substrate work function .
the schottky barrier is then simply calculated from the modified work function .
strong interaction ( chemisorption ) leads to the formation of bonds between the substrate metal atoms and the adsorbate sulfur atoms .
it increases the substrate work function if the electronegativity of the adsorbate is higher than that of the metal .
if the mos@xmath0  layer is chemisorbed , its electronic structure is perturbed to an extend that a schottky barrier can not be defined . however
, a second adsorbed mos@xmath0  layer then shows the characteristics of a single unperturbed layer .
we acknowledge taher amlaki and deniz ak@xmath123r for fruitful discussions
. this work is part of the research program of the foundation for fundamental research on matter ( fom ) , which is part of the netherlands organization for scientific research ( nwo ) .
the use of supercomputer facilities was sponsored by the physical sciences division ( ew ) of nwo .
67ifxundefined [ 1 ] ifx#1 ifnum [ 1 ] # 1firstoftwo secondoftwo ifx [ 1 ] # 1firstoftwo secondoftwo `` `` # 1''''@noop [ 0]secondoftwosanitize@url [ 0 ]
+ 12$12  & 12#1212_12%12@startlink[1]@endlink[0]@bib@innerbibempty link:\doibase 10.1038/nchem.1589 [ * * ,   ( ) ] link:\doibase    10.1021/cr300263a [ * * ,   ( ) ] link:\doibase 10.1038/nature12385 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.105.136805 [ * * ,   ( ) ] link:\doibase 10.1038/nnano.2012.193 [ * * ,   ( ) ] link:\doibase 10.1021/ar500274q [ * * ,   ( ) ] link:\doibase 10.1021/nn303513c [ * * ,   ( ) ] link:\doibase 10.1021/nl303583v [ * * ,   ( ) ] link:\doibase    10.1021/nl4010157 [ * * ,   ( ) ] link:\doibase 10.1038/srep01634 [ * * ,   ( ) ] \doibase    http://dx.doi.org/10.1063/1.4895767 [ * * ,   ( ) ] \doibase http://dx.doi.org/10.1063/1.4866340 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.91.161304 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.90.085415 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.90.201411 [ * * ,   ( ) ] link:\doibase 10.1021/jp212558p [ * * ,   ( ) ] link:\doibase dx.doi.org/10.1063/1.4893875 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.108.156802 [ * * ,   ( ) ] link:\doibase    10.1021/nl303909f [ * * ,   ( ) ] link:\doibase    10.1103/physrevx.4.031005 [ * * ,   ( ) ] link:\doibase    10.1021/nl403465v [ * * ,   ( ) ] link:\doibase    10.1039/c4tc01794c [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.89.245403 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.92.246401 [ * * ,   ( ) ] link:\doibase    10.1103/physrevb.76.125112 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.83.195131 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.47.558 [ * * , ( ) ] link:\doibase 10.1103/physrevb.49.16223 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.54.11169 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.59.1758 [ * * , ( ) ] link:\doibase 10.1103/physrevb.40.3616 [ * * , ( ) ] link:\doibase 10.1103/physrevb.46.16067 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.23.5048 [ * * , ( ) ] link:\doibase 10.1103/physrevlett.77.3865 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.76.073103 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.90.155448 [ * * ,   ( ) ] link:\doibase    10.1021/acsami.5b02782 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.84.195414 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.78.045409 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.106.186102 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.107.156401 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.86.085405 [ * * ,   ( ) ] \doibase    http://dx.doi.org/10.1063/1.4807855 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.101.026803 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.79.195425 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.108.235502 [ * * , ( ) ] link:\doibase 10.1103/physrevb.89.121103 [ * * ,   ( ) ] link:\doibase    10.1103/physrevb.87.155304 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.85.235407 [ * * , ( ) ] link:\doibase    10.1103/physrevb.85.033305 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.88.075409 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.88.245309 [ * * ,   ( ) ] link:\doibase    10.1021/nl4014748 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.108.156802 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.36.1647 [ * * ,   ( ) ] \doibase    http://dx.doi.org/10.1016/j.susc.2005.10.042 [ * * ,   ( ) ] http://stacks.iop.org/0022-3727/1/i=7/a=416 [ * * ,   ( ) ] \doibase http://dx.doi.org/10.1016/0022-0248(72)90129-7 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.64.235305 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.88.085318 [ * * ,   ( ) ] \doibase    http://dx.doi.org/10.1063/1.3524232 [ * * ,   ( ) ] link:\doibase    10.1021/nn506567r [ * * ,   ( ) ] link:\doibase 10.1109/led.2014.2313340 [ * * ,   ( ) ] link:\doibase 10.1021/am506921y [ * * ,   ( ) ] \doibase    http://dx.doi.org/10.1063/1.3696045 [ * * ,   ( ) ] link:\doibase    10.1021/nn500044q [ * * ,   ( ) ] link:\doibase 10.1021/nl5047736 [ * * ,   ( ) ]


entanglement is the property of states of multipartite quantum systems that is the most important resource for quantum information processing @xcite .
for this reason , one of the most important tasks of quantum information science is to establish entanglement at long distances in quantum networks , and to optimize final entanglement and probability of success .
quantum networks @xcite play a key role in quantum information processing . here
we limit our attention to those networks in which quantum states can be prepared initially and shared .
that is , entanglement can be generated between neighboring or , at least , not - too - widely - separated nodes ( or stations ) .
there are two instances in which the above - mentioned tasks become obviously relevant . on one hand
, one can consider macroscopic quantum communication networks , such as cryptographic networks , or more generally quantum communication nets @xcite , or distributed quantum computation @xcite involving arbitrary nodes of the network .
the second instance concerns microscopic or mesoscopic networks that could constitute architectures of quantum computers ( cf . ref .  ) .    despite enormous progress in experimental techniques ( cf .
ref .   and references therein ) , it is in principle a very hard task to establish entanglement at large distances due to decoherence and attenuation effects .
two remedies for this problem have been proposed :    * _ quantum repeaters_. this concept has been developed for 1d quantum communication chains @xcite .
although the simple entanglement swapping @xcite procedure can lead to quantum communication at large distances ( see fig .  [ fig : swapping ] ) , for imperfect resources , the performance of such communication chains decays exponentially with the distance ( i.e. the number of repeaters ) .
however , one can use more sophisticated quantum repeater protocols , which use purification and swapping methods that lead to polynomial decay only .
* _ entanglement percolation_. recently , our collaborators , together with two of us , proposed using networks in which properties of the connectivity of the network enable the establishment of , and determine the probability of , entanglement on large distances . in refs .   and  we considered in particular _ pure - state networks _ on regular lattices , where the nearest - neighbor ( nn ) nodes share a non - maximally entangled pair of qubits , or more generally qudits ( an entangled bond ) .
we searched for local operations and classical communication ( locc ) protocols that lead to establishment of entanglement between remote nodes of the network .
entanglement swapping and 1d repeater .
circles are qubits .
heavy lines represent pure non - maximally - entangled states .
loops represent the entanglement swapping measurement on a pair of qubits .
a ) two states @xmath0 and @xmath1 of the form specified in ( [ stateab ] ) . after the operation ,
qubits @xmath2 and @xmath3 may be in an entangled , mixed state .
b ) swapping as the first step in a repeater .
c ) a 1d chain of repeaters . ]    at present , only a handful of studies that explicitly apply percolation theory to problems in quantum information have appeared in the literature . these articles , comprising refs .   and  as well as applications to cluster states ,
are reviewed in ref .  .
our results leading to the present work can be summarized as follows :    * for 1d chains we proved that even optimal locc strategies only allow the establishment of entanglement between distant nodes that decays exponentially with both the distance , as well as the quality of entanglement of nn bonds . * in 2d and higher dimensions , the possibilities for protocols are greatly expanded .
the most straightforward , naive protocol the one we use as a baseline to evaluate other protocols is the one we term _ classical entanglement percolation _ ( cep ) , although it does involve obviously some quantum operations .
this protocol begins with converting ( using locc ) each of the entangled bonds into a singlet ( _ i.e. _ maximally entangled state ) with probability @xmath4 ( the so - called singlet conversion probability ( scp ) @xcite ) .
after the conversion , each of the parties ( nodes ) knows obviously which of the bonds are now perfectly entangled . using classical communications ,
the parties establish whether there exists an infinite percolating cluster ( or one spanning the lattice ) and who belongs to it .
then , by performing a series of entanglement swappings , it is possible to propagate entanglement between any two ( widely ) separated nodes that belong to the percolating cluster .
we call this scheme `` classical '' because it essentially maps the problem onto a classical bond percolation problem @xcite , and its success or failure is equivalent to the success or failure of bond percolation in the same lattice .
namely , if @xmath5 , where @xmath6 is the lattice - dependent critical percolation threshold ( or , in other words , critical open - bond density ) , then entanglement between any two remote nodes can be established with probability @xmath7 , which asymptotically does not depend on the distance . * at the same time that we introduced cep , we presented several schemes that went beyond the simple application of singlet conversion everywhere followed by entanglement swapping along a path .
we call these _ quantum entanglement percolation _ ( qep ) protocols because they use some kind of quantum pre - processing for instance quantum measurements to transform one percolation problem to another one . in these protocols cep
is also used , but is preceded by application of certain locc , which remove and replace bonds in the network lattice , resulting in a new lattice geometry with fundamentally different long - range properties . the pre - processing in qep
may greatly improve the possibility of entanglement percolation , either by reducing @xmath6 , or by increasing @xmath8 .
in this paper we present a more systematic and thorough study of qeps based on lattice transformations .
in particular , we explore the entire parameter space with the addition of monte carlo and series - expansion methods .
although general principles remain to be found , we do discover certain rules and regularities governing such strategies ( see sec .
[ sec : examples ] ) .
the paper is organized as follows . in sec .
[ sec : models ] we formulate and describe the models and fix the notation .
section [ sec : nets_and_perc ] discusses quantum networks and percolation , whereas sec .
[ sec : transformations ] deals with loccs used for lattice transformations .
our main results are presented in sec .
[ sec : examples ] , which contains descriptions of five examples of transformations enhancing cep : i ) the transformation of the kagom lattice to the square lattice , ii ) the transformation of the double - bond honeycomb ( hexagonal ) lattice to the triangular lattice , iii ) the transformation of the square lattice to two decoupled copies of the square lattice , iv ) the transformation of the bowtie lattice to the decoupled triangular and square lattices , and v ) the transformation of a triangular lattice with two different degrees of entanglement into the decoupled square and triangular lattices .
the cases ii ) , iii ) , and v ) were discussed already in ref .  , but we present here more general and stronger results .
we conclude in sec .
[ sec : conclusions ] .    over the past few decades
, classical percolation theory has seen the development of a number of quite sophisticated and powerful methods , which we have adapted to the questions at hand . in order to demonstrate the supremacy of qep over cep in each of the cases
i)-v ) we employ these methods combined with the methods of quantum information theory .
an overview of these methods is presented in the appendices : appendix  [ sec : majorization ] presents the necessary facts from majorization theory and singlet conversion protocols @xcite , appendix  [ sec : montecarlo ] explains some details of our monte carlo ( mc ) simulations , and finally appendix  [ sec : series ] deals with series expansions .
in this section , we describe the classes of quantum networks and communication protocols that we investigated .
the networks we consider consist , prior to application of communication protocols , of a collection of qubits partitioned into pairs , each pair being prepared in an identical pure state @xmath9 . with appropriate choice of bases
, any such state can be written @xmath10 where the schmidt coefficients @xmath11 satisfy @xmath12 and @xmath13 .
we identify these pairs with bonds on a two - dimensional lattice or edges on a graph , with two spatially - separated qubits , one occupying each end of the bond . at regular positions , a small set of these qubits
are arranged near enough to one another to allow measurements on any subset .
such a set of qubits then constitutes a vertex ( or node or site ) which is incident to each edge that contributes a qubit to the vertex .
this defines a lattice to which we apply methods of statistical physics .
in particular , many results on the possibility of long - range entanglement are described using percolation theory @xcite , and depend only on the graph structure of the system .      given a lattice prepared as described above , we search for the protocols consisting of local operations and classical communication ( locc ) that yield the maximum probability of achieving entanglement between two nodes separated by an arbitrarily large distance .
the answer depends on the single parameter @xmath14 characterizing the state @xmath15 .
even with a small palette of possible operations , finding the globally optimal solution is a difficult task .
instead , we search for promising classes of protocols .
the simplest protocol consists of attempting to convert the state associated with each bond to a singlet via the `` procrustean method '' of entanglement concentration , which is the optimal strategy at the level of a single bond @xcite .
this conversion succeeds with probability @xmath16 ( see appendix  [ sec : majorization ] . ) , while failure leaves the pair in a state with no entanglement . in this way the system is described exactly by a bond percolation process with open - bond density @xmath4 . if there is a path of open ( maximally entangled ) bonds connecting two nodes , a sequence of entanglement swapping measurements , one at each intermediate node ,
is then applied in order to entangle the first and last node .
this protocol is the simplest example of classical entanglement percolation .
here we review a few fundamental ideas in percolation theory necessary to analyze cep .
the nodes in the lattice can be partitioned into sets such that each node within a set is connected to each other node in the set via a path of open bonds .
such a set of nodes is called an open cluster , or sometimes simply a cluster .
the central fact of percolation theory is that percolation processes on most commonly - studied lattices in dimension @xmath17 and higher exhibit a continuous phase transition as the bond density passes through a critical value @xmath6
. for @xmath5 there exists with probability one a unique ( for the lattices we study here ) cluster of infinite mass ( number of nodes ) , while for @xmath18 all clusters are finite with probability one .
it follows that improving an locc protocol to obtain a small change in @xmath4 can have a dramatic effect on the probability of long - range entanglement . in the supercritical phase ,
long - range entanglement is possible , while in the subcritical phase it is not possible .
serving as the order parameter is the density of the infinite cluster @xmath19 which we define via @xmath20,\ ] ] the probability that a fixed node @xmath21 ( say , the node at the origin ) is in the the infinite cluster @xmath22 .
when referring to @xmath23 on a specific lattice , we use a symbol representing the lattice as a subscript .
the probability that two selected nodes are members of the same cluster decays roughly exponentially in their separation distance to an asymptotic value of @xmath24 .
( the length scale of decay is the correlation length @xmath25 . )
this means , for the problem at hand , that in cep the probability that information can be propagated between two nodes is asymptotically @xmath24 .
it was proven in refs .   and  that cep is not the optimal strategy for establishing long - distance entanglement .
the demonstration is based on applying certain locc prior to cep .
all of these pre - processing locc act on pairs of qubits , and they either transform the state on a given bond , or they replace two adjacent pure - state bonds by one , in general , mixed state bond .
there are essentially three types of generalized measurement used :    * _ singlet conversion _ the optimal locc singlet protocol @xcite is used in three situations . in cep with single - bond lattices we apply it directly to the state ( 1 ) , which results in the singlet conversion probability @xmath16 . if the protocol is successful , the bond that is converted to the singlet can be used for entanglement propagation ( swapping ) , otherwise it is useless . in cep with a double - bond lattice we apply it directly to the two copies of the state ( 1 ) , which lives in @xmath26 , and has schmidt coefficients @xmath27 , _ ergo _ the singlet conversion probability is @xmath28 finally , in all of the qep strategies that we study , we also apply singlet conversion to any remaining untouched bonds after the lattice transformation . *
_ entanglement swapping _ this protocol  @xcite , illustrated in fig .
[ fig : swapping]a , consists of performing the so  called bell measurement on a pair of qubits ( b and c ) in a node , i.e. a von neumann measurement in a basis of 4 maximally entangled orthonormal states ( in the computational basis ) .
it allows conversion of a pair of adjacent singlets into a singlet connecting end points ( a and d ) with probability 1 , i.e. allows for perfect entanglement propagation in a connected cluster of singlets . at the same time ,
when applied to a pair of imperfect states ( 1 ) , it produces a mixed state , which , amazingly , has the average singlet conversion probability equal to @xmath16 .
unfortunately , this effect can not be iterated . when applied to the mixed states ,
entanglement swapping reduces the singlet conversion probability .
this last point places a significant constraint on our choice of lattice transformations .
+ we use entanglement swapping in both cep and qep , but the two uses have very different effects . in the the case of cep , entanglement swapping is used to locally move entanglement between neighboring nodes , an operation that is repeated in hopes of transporting entanglement over long distances ( this is , roughly speaking , a brute - force method ) . but with qep , before attempting to transfer entanglement , we search for a way to selectively apply entanglement swappings to alter the geometry of the lattice and hence its long - range connectivity properties as given by percolation theory . the goal of this paper and future work is to enumerate the rich possibilities and point a way towards a general description of qep .
* _ worst case entanglement _ finally , in refs .   and   the worst case protocol was used , which maximizes minimal entanglement over all measurement outcomes .
this protocol consists also of bell measurement , but the basis is computational for one qubit , and corresponding to eigenstates of @xmath29 for the second . when applied to qubits , it produces a mixed state with the property that for all measurement outcomes , the resulting pure states have the same singlet conversion probability .
we will not use this protocol here .    in order to describe the lattice transformations involved in qep
, we need a more general formulation of the lattice than the one given in sec .
[ networksbipartite ] for cep .
it is useful to define the percolation processes a bit more precisely .
we begin with a graph , that is , a set of edges ( bonds ) @xmath30 and a set of vertices ( nodes ) @xmath31 .
we consider embedding the graph in @xmath32 in order to treat geometric properties .
in fact , the important properties do nt change if we force the vertices to occupy points in @xmath33 .
the configuration of open and closed bonds can be described by a probability sample space @xmath34 , with points @xmath35 , where @xmath36 takes the values @xmath37 and @xmath38 .
we allow each bond to be open with a different probability , _
i.e. _ @xmath39 .
the appropriate measure is a product measure on @xmath40 with marginal probabilities defined by @xmath41 , @xmath42 . to allow all transformations possible via locc , we take the graph to be complete that is , all possible edges @xmath43 with @xmath44 are present ( in fact , sometimes we need double bonds , as well . ) for percolation properties , any graph with fewer edges can be identified with this complete graph by setting the appropriate @xmath45 to zero .
( color online ) a ) mapping a quantum network to a percolation problem .
qubits are represented by small circles .
each vertex ( node ) contains six qubits .
b ) using entanglement swapping to transform the lattice structure . the two red ( grey ) bonds in a ) are replaced by the red bond in b ) .
this process can be continued to produce a double - bond hexagonal lattice . ]
now we describe the correspondence of entanglement on the physical system of qubits to this percolation formulation .
we denote the set of the indices of all qubits by @xmath21 .
each qubit @xmath46 is assigned to a vertex @xmath47 , with , in general , multiple qubits assigned to each vertex . for example , a portion of a network that can be described by a triangular lattice is shown in figure fig .
[ fig : show_qep]a , where vertex @xmath48 contains the six circled qubits . for every pair of qubits @xmath49 with @xmath50 ,
we denote the reduced state by @xmath51 , and by @xmath52 the singlet conversion probability ( scp ) , that is , the probability of conversion of @xmath51 to a singlet , maximized over all possible measurements
. then we assign @xmath53 for every @xmath54 .
initially , all reduced bipartite states have either @xmath55 for separable states , or @xmath56 for the prepared , pure , partially entangled states . in fig .
[ fig : show_qep]a the edges with @xmath57 are shown in black , while edges with @xmath55 are simply absent in the diagram .
two bonds are colored red ( grey ) only to show that they will be replaced by another bond ( in the sense of altering the scp ) via entanglement swapping .
for qep , a successful preprocessing , entanglement swapping operation such as that depicted in fig .
[ fig : swapping ] , alters the percolation process by setting @xmath58 and @xmath59 and @xmath60 .
likewise , in the example in fig .
[ fig : show_qep ] , an entanglement swapping measurement on qubits @xmath61 and @xmath62 sets @xmath63 , @xmath64 , and @xmath65 . each of the examples discussed below conforms to this description .
we consider five examples of lattice transformations in this paper , each exhibiting a different combination of features , with implications for the analysis of protocols .
in particular , two transformations convert one lattice to another , while the other three convert a lattice into two decoupled lattices .
four transformations involve lattices with single bonds , while the fifth involves double bonds .
quantities arising in the analysis , such as @xmath19 , are analytic about @xmath66 .
in fact , the lowest - order term in the expansion about @xmath66 is typically much larger than the remaining terms even relatively far from @xmath66 , with the result that most of the interesting crossover behavior in comparing protocols occurs for smaller values of @xmath4 .
but this behavior does not appear when using the distillation procedure for double bonds that leads to ( [ fourqubits ] ) .
the reason is that distillation produces a saturation point for @xmath19 smaller than @xmath66 .
this results in a much stronger difference between classical and quantum protocols in the high - density regime than does singlet conversion on single bonds ( two - qubit pure states ) .
four of the transformations produce a smaller critical density on at least one of the resulting lattices , which gives the most pronounced advantage in the regime near the critical density .
the fifth example shows that on some lattices where a particular qep strategy is not advantageous , allowing bonds of different strengths @xmath4 and @xmath67 can produce regions in the phase space @xmath68 where qep is indeed advantageous .
finally , consider comparing cep in which singlet conversion is applied to each bond separately , with qep that results in a single transformed lattice . in every such case
we find that qep is better than cep over the entire ( non - trivial ) range of @xmath4 .
it is an open question whether this is generic behavior .
we analyze the results of the transformation of each lattice in three regimes : near the critical density ( or densities ) ; near @xmath66 ; and between these two regimes .
arguments near the critical densities typically rely simply on the fact that long - range entanglement is impossible on a lattice with density below the critical density .
these results are the most insensitive to details of the definitions of entanglement and connectivity .
near @xmath66 we compute high - density expansions for @xmath19 and related quantities ( see sec .  [
sec : expansions ] . )
often the difference between the performance of cep and qep in this regime is marginal .
it is important nonetheless to carry out the analysis if we hope to make general statements about transformations that hold for all @xmath4 . between the critical regime and high - density regime
there are some techniques widely used in percolation theory that may be useful , such as russo s formula , which is used to prove inequalities in the rate of change of @xmath19 .
however , we leave these techniques for future work and use monte carlo computations in the present paper .      in our first example
we compare cep on the kagom lattice to quantum entanglement percolation consisting of transformation of the kagom lattice to the square lattice via entanglement swapping at nodes specified in fig .
[ fig : kagfive ] .
although we do not treat them here , there are at least two more ways to transform the kagom lattice to the square lattice using the same kind of entanglement swapping .
rigorous bounds have been obtained for @xmath69 @xcite ( as well as a high - precision monte carlo estimate @xcite ) while @xmath70 is known exactly , proving that the transformation gives an advantage for @xmath4 lying between @xmath70 and the lower bound for @xmath69 .
as shown in table  [ tab : series ] , the series for @xmath71 and @xmath72 are the same to the first non - trivial order in @xmath73 .
but the next term shows that @xmath74 for @xmath4 close enough to @xmath38 .
the mc data provides strong evidence that @xmath74 everywhere except near @xmath66 , where the statistical error is too large to distinguish the curves . but as is evident from the lower plot in fig .
[ fig : kagfive ] , the terms in the expansions that we computed dominate @xmath75 already at values of @xmath4 for which the mc data is still accurate .
thus we find that qep is advantageous over the entire range of @xmath4 .
one of three transformations of the kagom lattice to the square lattice .
pairs of qubits that are subjected to entanglement swapping are marked with loops . ]
_ @xmath16 .
the solid line is on the square lattice .
the dashed line is on the kagom lattice .
lower plot : normalized difference between @xmath71 and @xmath72 .
the solid line is computed from mc data .
the dashed line is computed from expansions from table  [ tab : series ] .
, title="fig : " ]   _ v.s .
_ @xmath16 .
the solid line is on the square lattice .
the dashed line is on the kagom lattice .
lower plot : normalized difference between @xmath71 and @xmath72 .
the solid line is computed from mc data .
the dashed line is computed from expansions from table  [ tab : series ] .
, title="fig : " ]      the hexagonal lattice with double bonds , each in the state specified by ( [ stateab ] ) ( see fig .
[ fig : hontrtwo]a ) , was discussed in refs .   and ,
in which it was shown that a transformation of the lattice to a triangular lattice offers an advantage over cep protocols for values of @xmath4 between @xmath76 and another critical value defined below .
( a ) double - bond hexagonal lattice .
( b ) transformed into the triangular lattice ( b ) .
pairs of qubits that are subjected to entanglement swapping are marked with loops . ]    here we extend the analysis to determine over the entire range of @xmath4 which of three protocols gives the highest probability of long - range entanglement .
the first of the three protocols , which we call cep i , consists of performing an optimum singlet conversion to each bond separately so that the probability of getting at least one singlet connecting two nodes is @xmath77 .
communication on the resulting lattice is then determined by bond percolation on the hexagonal lattice with bond density @xmath67 .
in particular , the critical density ( see table  [ tab : pcs ] ) is @xmath78 . in the second protocol
, we perform a more efficient conversion , namely distillation , on all four qubits in a double bond .
this protocol succeeds in producing a maximally entangled pair with probability @xmath79 ( see ( [ fourqubits ] ) and appendix  [ sec : majorization ] . )
this results in bond percolation on the hexagonal lattice with density @xmath80\right\}.\ ] ] we refer to this method as cep ii , with critical density @xmath81 . in the third method ,
the entanglement swapping that maximizes scp is applied at every other node ( see fig .  [
fig : hontrtwo ] ) , converting the double - bond hexagonal lattice to the triangular lattice with bond density @xmath4 , with @xmath82 .     of @xmath83 _ v.s .
_ @xmath16 for various strategies on the double - bond hexagonal lattice .
( a ) solid lines are monte carlo data .
dashed lines are series expansions .
( b ) solid lines are mc .
series expansions are indistinguishable from mc over the entire plot .
( c ) monte carlo data . [
fig : trihexone],title="fig : " ]   of @xmath83 _ v.s .
_ @xmath16 for various strategies on the double - bond hexagonal lattice .
( a ) solid lines are monte carlo data .
dashed lines are series expansions .
( b ) solid lines are mc .
series expansions are indistinguishable from mc over the entire plot .
( c ) monte carlo data .
[ fig : trihexone],title="fig : " ]   of @xmath83 _ v.s .
_ @xmath16 for various strategies on the double - bond hexagonal lattice .
( a ) solid lines are monte carlo data .
dashed lines are series expansions .
( b ) solid lines are mc .
series expansions are indistinguishable from mc over the entire plot .
( c ) monte carlo data .
[ fig : trihexone],title="fig : " ]    the three protocols are compared in fig .
[ fig : trihexone ] .
we first observe that each method fails below its corresponding critical density , and these are known exactly .
to compare the methods away from the critical points , we computed series expansions of @xmath19 about @xmath84 with results listed in table  [ tab : series ] .
defining @xmath85 we see that @xmath86 comparing ( [ thetacepi ] ) to @xmath87 from table  [ tab : series ] we see that for @xmath88 the conversion to the triangular lattice is better than cep  i. as is evident from figs .
[ fig : trihexone]a and [ fig : trihexone]b , the two curves have the same leading behavior as @xmath89 , with the result that the the monte carlo can not distinguish them in this region .
however , the series expansion lies well within the statistical error of the monte carlo points , even relatively far from @xmath66 , so we can be confident that @xmath90 even in the region where the mc can not distinguish the curves .
finally , the mc clearly shows that @xmath90 for smaller @xmath4 , where the expansion fails .
turning now to cep ii , we know that @xmath91 for @xmath92 and from ( [ distillp ] ) that @xmath93 for @xmath94 .
it follows that there must be a crossover point .
the mc data suggests that this occurs for @xmath95 . in summary
, we see that qep gives an advantage over cep i for all @xmath4 , but that cep ii , which is more efficient in its use of double bonds , is better than qep at high densities .
this transformation replaces the square lattice with two decoupled copies of the square lattice . to effect the transformation , at selected nodes the two horizontally opposing
bonds are joined into one bond and likewise with the vertically opposing bonds .
this procedure is applied at every other node , while staggering by one node when shifting rows as shown in fig .
[ fig : doubsq ] .
_ @xmath4 .
upper points : @xmath96 _ v.s .
_ @xmath4 .
solid lines are series expansions . ]
_ @xmath4 .
lower points : @xmath97 _ v.s .
_ @xmath4 . ]
because the transformation splits the original lattice into two disjoint lattices , each taking half the surviving nodes , we can not compare connection between a single pair of nodes before and after the transformation .
rather , we consider on the original lattice connection between , on the one hand , a pair of nodes @xmath98 separated by a small distance , and on the other hand , the same pair translated a distance much larger than the correlation length to a pair @xmath99 . we choose
@xmath21 and @xmath100 so that each one goes to a different lattice in the transformation .
after the transformation , there is the possibility of connection between @xmath21 and @xmath101 on one lattice and between @xmath100 and @xmath102 on the other .
we choose the nodes as shown in fig .
[ fig : doubsq ] .
in addition to the complication of splitting , this example is subtle because the critical density is the same on all three lattices . in ref .
it was shown that just above @xmath103 , the transformation gives an advantage over cep as measured by the choice of nodes @xmath98 .    here
we examine the connectivity for all @xmath104 .
the probability that at least one of the pairs on the decoupled lattices is connected is @xmath105 we compare @xmath106 to the probability on the original lattice that at least one of @xmath21 and @xmath100 is connected to at least one of @xmath101 and @xmath102 .
this probability is @xmath107 where @xmath108.\ ] ] ( in this paper ` or ' in ` @xmath109 or @xmath110 ' does _ not _ mean exclusive or . )
we prove that the doubling is advantageous in the limit @xmath88 using the power series about @xmath84 for @xmath71 and @xmath111 given in table  [ tab : series ] .
this follows from comparing @xmath112 and @xmath113 moreover , monte carlo evidence suggests that the transformation improves on cep for all @xmath114 in fig .
[ fig : pi_pdoub_asymp ] one sees from the accuracy of the series expansion in the region of high - quality mc data that that the advantage is maintained even when the monte carlo becomes noisy near @xmath66 .
the same mc data is not useful near @xmath70 for two reasons .
firstly , both @xmath115 and @xmath107 vanish at the same critical point with infinite slope .
secondly , the systematic error due to finite lattice size and large fluctuations in cluster statistics further complicate distinguishing the curves .
we instead make use of an alternate expression for @xmath116 , that is @xmath117 where @xmath118 $ ] .
the condition for advantage over cep then becomes @xmath119 the mc data for @xmath120 was generated by considering the largest cluster in the finite lattice to represent the infinite cluster @xmath22 even if it is not a spanning cluster .
the mc data supports ( [ alternate ] ) , and furthermore , shows no evidence of non - analyticity at @xmath70 , as is evident in fig .
[ fig : theta_omega_nolog ] .
the foregoing analysis of the doubled square lattice was based on the choice of @xmath98 shown in fig .
[ fig : doubsq ] . but
this is not the only reasonable choice .
in fact the question of whether the doubling transformation is better than cep depends somewhat on the details of how the nodes are chosen .
although this is an extreme example , similar questions arise in analyzing other transformations .
thus , the ambiguity in the measures comparing the various protocols must be addressed in future work .
transformation of the bowtie lattice to decoupled square and triangular lattices .
loops marking pairs of qubits represent swapping measurements . ]
this is our second example that transforms a lattice into two decoupled lattices .
figure [ fig : bowtie ] shows the measurements that decouple the bowtie lattice into the square lattice and the triangular lattice . for @xmath4 satisfying @xmath121 ( the exact critical values are listed in table  [ tab : pcs ] ) , the transformation is obviously advantageous .    in carrying out further analysis , the bowtie lattice presents a complication not present in other lattices studied in this paper : while the other lattices are regular in the sense that each vertex has the same environment up to rotations and reflections , the bowtie lattice has two kinds of vertices in this sense .
notice first that all the surviving nodes on the square and triangular lattices shown in fig .
[ fig : bowtie ] were generated from nodes of coordination number @xmath122 on the bowtie lattice , with all nodes of @xmath123 disappearing .
each node of @xmath122 on the bowtie lattice has four nearest neighbors of @xmath122 , discounting nodes of @xmath123 ( for example node @xmath21 in fig .
[ fig : bowtie ] . )
furthermore , two of these nearest neighbors are connected by a diagonal bond ( node @xmath124 ) and two are not ( node @xmath100 . )    in treating the example of doubling the square lattice in sec .
[ sec : doubsquare ] , the fact that the lattice decouples forced us to consider connections between two widely separated pairs of nodes .
we treat the present example in the same way , except that the more complicated local structure of the bowtie lattice forces us to consider a cluster of three nodes @xmath125 rather than a pair . in fig .
[ fig : bowtie]b we see that @xmath126 are sent to the square lattice and @xmath21 is sent to the triangular lattice .
we consider a distant cluster of nodes @xmath127 related to @xmath125 by translation , and examine the probability that a connection exists between at least one of @xmath125 and at least one of @xmath127 . on the bowtie lattice , this probability is @xmath128 where @xmath129.\ ] ] if instead we decouple the lattices via the lattice transformation , this probability is given by @xmath130 where @xmath131.\ ] ] a comparison of monte carlo estimates and high - density expansions of these quantities is shown in fig .
[ fig : bt_sq_tr_log ] .
monte carlo estimates of @xmath132 and @xmath133 .
the inset is the same data on a linear scale .
solid lines are high - density expansions . ]
the mc data shows that for small densities , the lattice transformation gives an advantage , while for higher densities cep is the better protocol .
the cross - over occurs at @xmath134 . as with the other examples , the curves representing the series expansions suggests that no cross - over occurs in the high - density region .
.[tab : pcs ] @xmath6 for bond percolation on some lattices .
@xmath135 is the unique root of @xmath136 .
all critical densities are exact@xcite except for @xmath69@xcite . [ cols="<,<",options="header " , ]       a ) the asymmetric triangular lattice .
solid bonds represent the state @xmath15 ( density @xmath4 ) .
dashed bonds represent the state @xmath137 ( density @xmath67 ) .
b ) entanglement swapping is performed on pairs of solid bonds resulting in a triangular lattice with bond density @xmath4 .
the remaining dashed bonds form a kagom lattice which is then transformed into a square lattice with bond density @xmath67 as in sec .
[ sec : kag ] . ]    with this example we demonstrate that qep can succeed when the initial bonds are not all in the same state , but fail if they are in the same state . consider the lattice in fig .
[ fig : asymtr]a composed of two different kinds of bonds , each of the form given by ( [ stateab ] ) .
the solid and dashed bonds represent the states @xmath15 ( with @xmath16 ) and @xmath137 ( with @xmath138 ) respectively . in general ,
we consider two entanglement distribution protocols : a ) classical entanglement percolation on the original lattice and b ) qep consisting of creating two decoupled lattices ( the square and triangular lattices ) via entanglement swapping , followed by cep on each of the resulting lattices . the transformation is described in fig .  [ fig : asymtr ] .
note that , for @xmath140 , the initial lattice is a triangular lattice with serial double bonds and therefore the critical point @xmath141 . on the other hand for @xmath142 and @xmath143
the initial lattice is the kagom lattice .
we first examine the most robust measure of a protocol s effectiveness the binary measure that tells whether long - range entanglement is possible or not .
the phase diagrams for this example before and after the transformation are shown in fig .
[ fig : asymtr_phase ] .
we see that there are regions in the phase space for which qep is better than cep , and vice versa .    in the region in which both protocols allow long - range entanglement , a more detailed measure similar to those in previous sections
is necessary .
for example , for @xmath144 , we performed an analysis similar to the one used for doubling the square lattice .
the two critical boundaries intersect at @xmath145 , but comparing connectivity just above the critical point shows that the transformation does not improve the probability of long - range entanglement .
likewise for @xmath144 and @xmath4 near @xmath38 , series expansions showed that the transformation is not an improvement on cep .
we have not yet determined whether qep in this scenario succeeds for some other @xmath144 , but this seems unlikely .
phase diagram for the asymmetric triangular lattice .
the heavy solid curve ( obtained by monte carlo ) separates the supercritical and subcritical regions on the original asymmetric triangular lattice ( cep ) .
the dashed curve composed of two line segments separates the supercritical and subcritical regions after the transformation ( qep ) .
after transformation , the supercritical region is defined to be the region for which at least one of the disjoint lattices is supercritical . in the light gray regions , long - range entanglement is possible with qep , but not with cep . in the dark gray region
it is possible with cep , but not with qep . ]
in this paper we have considered the problem of establishing long - distance entanglement in quantum pure - state networks on regular 2d lattices .
we have discussed in detail several examples of quantum entanglement percolation strategies that are better than the corresponding classical strategies , i.e. those consisting of a direct attempt to convert bonds into singlets .
our results illustrate nicely the interplay between quantum information theory and classical percolation theory . despite the fact that we do find certain rules and regularities governing qep strategies , many questions remain open .
for instance , we can not say anything about optimality of our qep protocols and most probably they are not optimal .
all of our protocols involve locc acting on pairs of qubits only ; can multipartite locc , which leads inevitably to creation of multipartite entanglement , help ? we know that for sufficiently large initial entanglement , perfect entanglement between remote nodes may be established with distance - independent probability @xmath7 .
is it also possible for arbitrarily small initial entanglement ?
we thank a. acn , j.i .
cirac , and d. cavalcanti for discussions .
we thank s. perseguers for kindly supplying figures .
we acknowledge support of the eu ip programme `` scala '' , and spanish mec grants ( fis 2005 - 04627 , consolider ingenio 2010 `` qoit '' ) .
j.w . thanks lluis torner and icfo for hospitality .
he was partially supported by nsf grant dms 0623941 .
m. l. thanks also the alexander von humboldt foundation for support .
majorization theory has been applied to questions of transforming one bipartite pure state to another by locc . in particular , a theorem due to vidal@xcite
gives the probability that such a transformation can be achieved via an optimal protocol ( without specifying that protocol . )
we state the result and apply it to the distillation procedure used in this paper .
we must first introduce a certain partial order on vectors .
consider two real , @xmath3-dimensional vectors @xmath146 and @xmath147 .
we define the vector @xmath148 by reordering the elements of @xmath146 into non - decreasing order , and likewise with @xmath147 .
we say that @xmath146 is submajorized by @xmath147 , denoted by @xmath149 , if @xmath150 for all @xmath151 . denoting the vector of schmidt coefficients of @xmath152 by @xmath153 ,
the theorem states that @xmath154 can be transformed into @xmath155 with probability @xmath4 , where @xmath4 is the largest number on @xmath156 $ ] such that @xmath157 .
a simple , relevant application is the computation of the probability that optimal conversion of a state @xmath158 to a singlet will succeed .
application of ( [ submaj ] ) gives @xmath159 .
now we consider distillation , defined here as the optimal protocol for converting @xmath160 pure states @xmath161 to @xmath162 product states and one maximally entangled state in the schmidt bases .
this operation is useful , for instance , in attempting to get a single , maximally - entangled bond from @xmath160 bonds connecting two nodes .
explicitly , @xmath163 with @xmath164 and elements of the ordered set @xmath165 non - increasing . as usual , @xmath166 is shorthand for @xmath167 . for @xmath168
define the bijective numeration @xmath169 and @xmath170 .
then we can write @xmath171 in order to apply the theorem , we need to show that this state can be written as a bipartite state with schmidt coefficients that can be chosen to satisfy our distillation condition . to this end
collect all the first qubits of the bipartite states and all the second qubits , defining @xmath172 so that ( [ refactor ] ) becomes @xmath173 note that @xmath174 so that ( [ bigbipartite ] ) defines a state in @xmath175 already in schmidt form .
the submajorization condition is @xmath176 for which the only nontrivial inequality is the penultimate one @xmath177 .
thus the maximum distillation probability is @xmath178 note that this result agrees with the special case in ( [ fourqubits ] ) .
we computed monte carlo estimates of @xmath19 and @xmath179 using the hoshen - kopelmann @xcite algorithm with modifications for efficiency @xcite and the mersenne twister @xcite random number generator .
together with series expansions and exactly known critical densities , the quality of the data we obtained is more than sufficient to determine which transformations are advantageous .
we estimated @xmath19 by computing the mean density of the largest cluster on an @xmath180 lattice .
near the critical density we typically used lattices of size @xmath181@xmath182 .
we took the value of @xmath183 at the inflection point of plots of @xmath19 _ v.s .
_ @xmath183 at fixed @xmath4 as estimates of the correlation length @xmath25 .
we typically found that @xmath184 for @xmath185 . we only expect significant systematic error in this region , but this poses no problem because the curves are never close to one another in these regions .
( we reformulated the problem when this is the case . )
we computed statistical error , but the error bars are at most barely visible on the plots , so we omitted them .
exceptions are very near @xmath186 where fluctuations in the size of the largest cluster become large , and near @xmath66 , where finite clusters are rare , so collecting sufficient samples to distinguish curves is too expensive .
we discuss the effects of these errors in the main body of the paper .
we wrote a single computer code to study all the lattices .
the code supports lattices with vertices that occupy points on @xmath33 with bonds connecting each pair of nearest neighbors as well as a diagonal bond from @xmath187 to @xmath188 .
because the connectivity properties that we calculated depend only on the graph structure of the lattice , we embedded the graph of each lattice in the square lattice plus diagonals described above .
depending on the lattice being modeled , some of the bonds in the underlying lattice are closed with probability one , and some are open with probability one , with the vertices identified in the embedded lattice .
the graph structures of all lattices appearing in this paper were modeled in this way .
we follow the ideas of the perimeter method @xcite to compute high - density series expansions of @xmath19 and @xmath179 listed in table  [ tab : series ] . here
we discuss the method for computing the series for @xmath19 , but our method for computing @xmath179 is similar . here , a cluster is any connected subgraph that contains at least one vertex .
we denote by @xmath189 the collection of all finite clusters that include the vertex at the origin , and by @xmath190 the partition of @xmath189 induced by equivalence under translation , in other words the collection of free clusters . choosing an enumeration @xmath191 of clusters in @xmath190 , the probability that a randomly selected site is in the infinite open cluster
is easily seen to be @xmath192 where @xmath193 is the number of sites and @xmath194 the number of bonds in cluster @xmath191 , and @xmath195 is the number of perimeter bonds ( bonds adjacent to cluster @xmath191 ) .
clearly , we can find the series expansion in @xmath196 by enumerating the clusters in an order that is non - decreasing in @xmath195 .
for instance , on the hexagonal lattice , only the cluster consisting of an isolated site has a perimeter @xmath197 of size less than or equal to @xmath198 ( see fig .  [
fig : theta_hex_b ] ) , so that to lowest non - trivial order @xmath199 .
percolation theory is a difficult subject precisely because the full enumeration is difficult .
the series for @xmath200 ( for bond percolation ) was calculated by machine to high order in ref .  .
although tables of cluster numbers have been published , all the others that we are aware of are either for site percolation or for bond percolation with cluster size measured by the number of bonds rather than sites ( neither of which can be mapped to our problem . )    for our results we counted a few small clusters by hand , which is not so difficult . in practice
however , we find that it is also essential to categorize the clusters by symmetry . as an example , the clusters contributing to @xmath23 on the hexagonal lattice to sixth order in @xmath196 are shown in fig .
[ fig : theta_hex_b ] .    10 pt   clusters on the hexagonal lattice contributing to @xmath23 to sixth order in @xmath196 .
symmetries and multiplicities @xmath201 are : none@xmath202 , one axis@xmath203 , rotation@xmath204 , two axes@xmath198 , rotation and axis@xmath17 , all@xmath38 .
the number of contributing clusters per figure is given by @xmath205 , the factor of @xmath206 accounting for the fact that each site only supports half the rotations .
the contribution to @xmath23 is then @xmath207 .
, title="fig : " ]