recent years have witnessed a surge of interest in understanding the microscopic origin of friction as a result of the increased control in surface preparations and the development of local probes like the atomic force microscopes ( afm ) .
one of the goals of this research is to understand whether extremely low friction can be obtained by an appropriate choice of the sliding conditions .
in particular , commensurability between the sliding lattices is one of the elements that determine friction . for a purely incommensurate infinite contact
, theoretical arguments suggest that static friction should vanish @xcite .
this effect has been called superlubricity @xcite .
commensurability and incommensurability are defined in terms of lattice parameters .
however , specific orientations can also lead to ( mis)matches of lattices in contacts .
very low friction was found in experiments with finite contacts at very low velocities to depend strongly on the orientation  @xcite .
coexisting states of very different friction have also been observed in the sliding of antimony nanoparticles @xcite and have been attributed to contamination or ( in)commensurate interfaces .
meanwhile , recent theoretical studies  @xcite have shown that nanocrystals can slide with constant orientation only for particular orientations .
this paper examines systematically theoretically the ( in)commensurability and friction of sliding nanocrystals with a triangular lattice symmetry [ such as au ( 111 ) ] on a triangular or hexagonal substrate ( such as graphite ) .
gold nanocrystals sliding on graphite are a prototype system for friction that is being investigated both experimentally  @xcite and computationally  @xcite . in this work
, the friction is investigated analytically through the total potential energy of the contact layer on the substrate .
the potential energy corrugation plays an essential role in the survival and appearance of stable ( in)commensurate sliding orientations  @xcite , as well as the order of magnitude of the friction . in the regime of low sliding velocity , and low temperature , which is the typical situation in afm experiments  @xcite ,
the friction is of the order of @xmath4 , with @xmath5 the substrate period and @xmath6 the corrugation of the total potential energy of the interface .
the interaction between the surface and an atom of the contact layer is modelled with a realistic static potential .
the scaling of the friction with contact size is determined , depending on the orientation , and conclusions are drawn from this regarding commensurate and incommensurate orientations .
the calculations are repeated also for contacts between crystals with rectangular lattice geometries .
an illustration of the systems is shown in fig .
[ fig : cartoon ] .
an illustration of the systems under study .
a nanocrystal lies on a triangular or hexagonal substrate with a contact layer that has triangular lattice symmetry .
systems with rectangular lattice symmetry are discussed as an example in the appendix .
, width=325 ]    in sec .
[ sec : geometry ] the geometry of gold on graphite is introduced .
an expression for the potential energy is derived in sec .
[ sec : maths ] . at first ,
rigid crystals are considered , followed by correction terms for crystals with relatively stiff , but deformable latices . in sec .
[ sec : scaling ] , the scaling of the friction is discussed at different ( in)commensurate orientations and orientations with different scaling are identified and classified .
it is found that for perfect crystals there are three types of ( in)commensurability and consequently three possibilities for the scaling of the friction with the size of the contact area @xmath0 , namely @xmath7 , and @xmath0 .
conditions for each are derived .
these three different scaling behaviors for the friction are a direct consequence of the two - dimensional nature of the contact layer , and are thus possible in all geometries , not just triangular . in the appendix ,
the calculations are repeated for a rectangular geometry , to exemplify this .
the scaling of friction for imperfectly and irregularly shaped nanocrystals is discussed in sec .
[ sec : imperfectscaling ] . in sec .
[ sec : discussion ] , the implications for experiments are discussed , and a method is proposed for quantifying the commensurability of finite - size crystal interfaces .
finally , conclusions are presented in sec .
[ sec : conclusions ] .
the interaction between the graphite substrate and the gold crystal are dominated by the bottom , contact layer of the gold crystal and the top layer of the graphite .
commensurability of these two layers is not only controlled by the ratio of the lattice parameters , but also by their relative orientation .
the surface of graphite has a hexagonal lattice of carbon atoms with an inter - atomic distance of @xmath8 .
a two - dimensional hexagonal potential for a gold atom at postition @xmath9 on a graphite substrate can be written as @xmath10 where @xmath11 \vec{e}_x + \sin[(l-1 ) \pi/3 ] \vec{e}_y$ ] is the unit vector in the direction @xmath12 with respect to the @xmath13-axis , and @xmath14 . without any load force ,
the potential corrugation of a single gold atom on graphite is denoted by @xmath15 and is around @xmath16  @xcite .
a load force can be accounted for with a higher corrugation . with a negative @xmath15
, this potential can also be used to describe a substrate with a triangular lattice .
while in general a substrate potential may not be sinusoidal , it is always periodic and hence can be written as a sum of sinusoidal potentials . the analysis described in this work
can also be applied to sums of different sinusoidal potentials .
more details of the geometry are shown in fig .
[ fig : trapezoid ] .    from the triangular shape of the nanocrystal found in experiments  @xcite
, it can be deduced that the contact layer must have similar triangular symmetry .
as bulk gold has an fcc lattice , it is reasonable to conclude that the contact layer must be the ( 111 ) cleavage . in the bulk ,
the lattice parameter is @xmath17   , and consequently the inter - atomic distance in the ( 111 ) layer is @xmath18 .
however , the au(111 ) surface is known to show surface reconstruction under a wide range of conditions  @xcite . in vacuum
, the reconstructed au(111 ) surface looks similar to the unreconstructed bulk ( 111 ) layers , but has a superstructure with a unit cell of size @xmath19 and an average compression in the ( 110 ) direction .
consequently , the rotational symmetry of the contact layer is lost .
in addition , the surface gold atoms do not lie in a plane , but are shifted by up to 0.15   .    as the interaction between the graphite substrate and the au atoms is weak compared to the interaction between the au atoms , the reconstruction of
the contact layer in the au - graphite configuration is likely similar to that of au(111 ) in vacuum . in this work
, we simply assume that the lattice is only slightly distorted , with the atoms arranged in scalene triangles .
as the interaction between the graphite surface and the gold atoms is weak , the shift of 0.15    orthogonal to the surface does not lead to significant changes in the interactions .
as the results derived in this work can be generalized to any regular lattice , they can be applied to any reconstructed surface .    in general , the nanocrystal can be flexible , and may be deformed due to forces exerted by the substrate . in this work ,
the crystal is at first assumed to be rigid , and then corrections are made for the displacements of atoms with respect to their equilibrium position .
elastic deformations can be neglected  @xcite , as the typical coherence length of gold , about 1@xmath20 m , is much larger than the size of a crystal in friction experiments  @xcite .
young s modulus , which for gold is @xmath21  gpa , gives the stress - strain response and can be used as an indication of the displacement @xmath22 of a single atom with respect to the equilibrium lattice positions .
one may consider the force exerted by the substrate on an atom at position @xmath9 , @xmath23 to work on a single unit cell of the contact layer , of diameter @xmath24 .
the displacement of the atom can be estimated as @xmath25 this leads to a typical displacement of a single gold atom in the crystal of about @xmath26   , corresponding to a difference in potential energy two orders of magnitude smaller than the corrugation .
though young s modulus is a bulk property , and we are interested in the displacement of surface atoms , it is sufficient for obtaining the order of magnitude . a similar estimate for the order of magnitude of the displacement of a surface atom can be obtained from the total number of neighboring atoms and parameters used in molecular - dynamics simulations with atomistic force fields , such as that of ref .  @xcite . in the next section , for the total potential energy of the contact layer , the contribution from the displacement of surface atoms are neglected . in sec .
[ sec : systematic_soft ] , corrections that include the displacement of the atoms with respect to their equilibrium positions are also made .
it has been shown that sliding crystals rotate to specific orientations , which are stable and remain ( nearly ) constant for all time  @xcite .
these orientations can be identified with periodic orbits in the dynamics and obtained from the potential energy landscape .
they occur when the potential energy averaged over a scan line and its corrugation simultaneously exhibit an extremum as a function of orientation . however , the incommensurate orientations typically have higher average energy than the commensurate ones , and can be destroyed by sufficiently large thermal fluctuations , leading to commensurate sliding and an increase of friction . in experiments of small graphite flakes on graphite , incommensurate orientations were found to decay  @xcite .
larger crystals were found theoretically to rotate more slowly , so that incommensurate orientations survive .
the sliding gold nanocrystals studied in the experiments by dietzel et al .
@xcite are sufficiently large for this ( with contact areas between @xmath27 and @xmath28 nm@xmath29 ) , though the smaller crystals studied numerically by guerra et al .
@xcite , can rotate easily while sliding . in order to determine the scaling of the friction , and hence the commensurate orientations , as well as
apply the theory of ref .
@xcite , it is necessary to focus on the potential energy of the contact layer on the substrate as a function of position and relative orientation .
diagram of the hexagonal graphite substrate lattice ( left ) and a trapezoidal gold crystal ( right ) with definitions of the size parameters @xmath30 , lattice spacings @xmath31 , shape parameter @xmath32 , substrate lattice vector @xmath33 , orientation @xmath34 , and position @xmath35 .
the dark atoms shown complete the trapezoid crystal into a triangular one .
we use @xmath36 and @xmath37 . for gold , with the reconstruction that has been observed in vacuum  @xcite ,
the parameters for the lattice are @xmath38   , @xmath39   , @xmath40 , @xmath41 . in a perfectly triangular crystal , @xmath42 .
, width=325 ]    let us calculate the total potential energy of a trapezoidal crystal .
the geometry of such as crystal is displayed in fig .
[ fig : trapezoid ] , along with definitions of the diameters @xmath43 and @xmath44 , the inter - atomic distances @xmath45 and @xmath46 , the shape parameter @xmath32 , the orientation @xmath34 , and the position of the center of smallest triangle that contains the trapezoid , @xmath35 .
first , we assume that the crystal is rigid , and then make a systematic correction for displacement of atoms from their equilibrium position .
the total potential energy of the crystal at position @xmath47 and orientation @xmath34 on the substrate is dominated by the interaction of the contact - layer atoms with the periodic substrate potential .
it can therefore be written as a sum over all atoms in the contact layer , @xmath48 where @xmath49 are the positions of the atoms in the contact layer relative to @xmath35 , the center of the contact layer .
if the crystal is rigid , the positions of the atoms can be written as @xmath50 with @xmath51 the @xmath52 matrix that performs a rotation over an angle @xmath34 .    by substituting the substrate potential , eq .
( [ eq : substratepotential ] ) , and rewriting the @xmath53 in terms of imaginary exponentials , we find that for a rigid crystal ,    @xmath54~. \label{eq : s}\end{aligned}\ ] ]    due to the rotational symmetry of the substrate and of a perfectly triangular crystal , @xmath55 consequently , for an irregular trapezoid , it suffices to determine @xmath56 . finally , the sums in eq .
( [ eq : s ] ) can be evaluated by making use of the relation @xmath57}/{\sin\left({c } \right)}~ , \label{eq : sinratio}\end{aligned}\ ] ] with @xmath58 any real number .
one finds that @xmath59\right\}\nonumber\\ & \phantom{\sum_{j=0,1}}\times\frac{\sin\left\ { \frac{\pi b}{\alpha } ( d^\prime+1 ) \left [ \tfrac12 \sqrt{3 } \cos\phi + \left(\gamma + \tfrac12 ( -1)^j \beta \right ) \sin\phi \right]\right\}}{\sin\left\ { \frac{\pi b}{\alpha } \left [ \tfrac12 \sqrt{3 } \cos\phi + \left(\gamma + \tfrac12 ( -1)^j \beta\right ) \sin\phi \right]\right\ } \sin\left(\frac{\pi b}{\alpha } ( -1)^j \beta \sin\phi\right)}~ , \label{eq : pottrap2 } \\ v_{\mathrm{trap . } , d , d^\prime,\alpha,\beta,\gamma}^{(0)}(\vec{r},\phi ) & = \null- \tfrac19 v_0 \sum_{l=1}^3 \sin\left [ \frac{\pi b}{\alpha}\left ( \frac{\vec{r}_\triangle\cdot\vec{e}_l}{b}\right ) \right ] s_{d , d^\prime,\alpha,\beta,\gamma}\left(-\tfrac12 \alpha\vec{e}_l , \phi , l\right ) \nonumber\\ \label{eq : pottrap1 } & \phantom{= \null}{}\null - \tfrac19 v_0 \sum_{l=1}^3 \cos\left [ \frac{\pi b}{\alpha}\left ( \frac{\vec{r}_\triangle\cdot\vec{e}_l}{b}\right ) \right ] s_{d , d^\prime,\alpha,\beta,\gamma}\left(0,\phi , l\right ) ~.\end{aligned}\ ] ] note that the angle - dependent corrugation prefactor contains all the geometric dependence on size and orientation , and does not depend on the position @xmath47 of the center of mass on the substrate .
the total corrugation is @xmath60^{\tfrac{1}{2 } }   \left ( \tfrac32 \cos \frac{\xi}{3 } + \tfrac12 \sqrt{3 } \sin \frac{|\xi|}{3 } \right)~,\\ \xi & = \arctan\left(\frac{s_{d , d^\prime,\alpha,\beta,\gamma } \left ( -\tfrac12 \alpha\vec{e}_1 , \phi,1 \right)}{s_{d , d^\prime,\alpha,\beta,\gamma } \left(0 , \phi,1\right)^2}\right)~,\end{aligned}\ ] ]    where @xmath61 and @xmath62 are used to denote the value of the maxima and minima of the potential respectively .
the average potential energy over any scan line vanishes for all scan lines not exactly orthogonal to one of the @xmath63 . in this case , the conditions for the existence of stable orientations of ref .
@xcite are met . as is also the case for rectangular lattices in ref .
@xcite , the incommensurate stable orientations are distributed evenly over the unit circle and their number grows linearly with the diameter  @xmath43 .
it should be noted that these results are consistent with the numerical results for the distribution of activation energies described in ref .
@xcite for rigid crystals with triangular symmetry on a square lattice .
however , because those results were obtained for relatively small clusters , the authors could not distinguish the commensurate orienations with larger corrugation .      in realistic nanocrystals , atoms of the contact layer
are displaced with respect to their equilibrium position , due to the forces exerted on them by the substrate . from the arguments presented in sec .
[ sec : geometry ] , it is clear that this displacement is small for the system studied in this work .
however , as the contact layer of typical nanocrystals contain many atoms , the energy effect is not necessesarily negligeable compared to the low total potential energy for incommensurate configurations .    to provide some insight into the effects of lattice deformations without making the calculations much more complicated ,
it is assumed here that each atom couples only to its equilibrium position . within this mean - field approximation , small displacements
@xmath64 from the equilibrium position are due to a force equal to @xmath65 where the spring constant @xmath66 can be estimated from eq .
( [ eq : displacementfestimate ] ) to be , @xmath67    because typical sliding velocities in experiments are low , we may limit ourselves to the quasi - static case , where all forces on each atom sum up to zero .
the correction to the potential energy consists of the reduction in potential energy due to the displacement on the lattice , and the additional potential energy stored in the spring . to leading order in the displacement ,
we find a correction to the potential energy of the atom with equilibrium position  @xmath9  of @xmath68 with @xmath69 the force exerted by the substrate on an atom at position @xmath9 .    by substituting eq .
( [ eq : substratepotential ] ) , using the relation @xmath70 , and using the properties of @xmath63 , one obtains    @xmath71 ^ 2\\ & = - \frac{4 \pi^2 v_0 ^ 2}{81\kappa \alpha^2 } \sum_{l=1}^3 \left [ 1 - \cos\left(\frac{4 \pi}{\alpha }
\vec{e}_l\cdot\vec{r}\right ) + \cos\left(\frac{2 \pi}{\alpha } \vec{e}_{l } \cdot\vec{r}\right ) -\cos\left(\frac{2 \sqrt{3}\,\pi}{\alpha } \vec{e}_{l+\tfrac12 } \cdot\vec{r}\right ) \right]\end{aligned}\ ] ]    this expression is of a form very similar to eq .
( [ eq : substratepotential ] ) . the approach used in sec .
[ sec : rigidpotential ] to obtain the total potential energy of a rigid crystal can be applied to this expression as well .
hence , one obtains the first - order correction of @xmath72 } \left .
\null - v_{\mathrm{trap . } , d , d^\prime,\alpha/2,\beta,\gamma}^{(0)}(\vec{r},\phi ) - v_{\mathrm{trap . } , d , d^\prime,\alpha/\sqrt{3},\beta,\gamma}^{(0)}({\mathcal r}(\pi/6 ) \cdot \vec{r},\phi+\pi/6 ) \right ] ~. \label{eq : firstordercorrection}\end{aligned}\ ] ]    though the functional forms of the terms in this correction are very similar to eq .
( [ eq : pottrap1 ] ) , the values of the parameters are different . as is discussed in the next section , this allows the first - order correction to become larger than the leading order in some cases , despite the small prefactor .
higher order corrections and corrections for direct coupling to the neighboring atoms in the contact layer could , in principle , be obtained using the same approach . however
, these corrections would have to include the nonlinear response of the lattice , as well as second - order derivatives of the potential energy .
this would lead to expressions similar in form to eq .
( [ eq : substratepotential ] ) , and therefore to additional terms similar to eq .
( [ eq : firstordercorrection ] ) , but with yet again different parameters and smaller prefacors .
from the structure of the function @xmath73 in eq .
( [ eq : pottrap2 ] ) , one can derive important general geometrical and scaling properties . from eqs .
( [ eq : pottrap2 ] ) ,  ( [ eq : pottrap1 ] ) , and  ( [ eq : firstordercorrection ] ) it is clear that there are a number of parameters which control the potential energy corrugation of the surface which the nanocrystal is subjected to . for a perfectly triangular crystal these are the diameters @xmath43 , the ratio of the lattice constants , @xmath74 , and the orientation @xmath34 .      in eq .
( [ eq : pottrap2 ] ) , the numerator varies rapidly between  -1 and  1 , with @xmath34 , while the denominator provides a size - indepent prefactor that can diverge for particular angles .
this can be seen in figs .
[ fig : numdenom ] and  [ fig : numdenom2 ] , where one of the angle - dependent corrugation prefactors is plotted along with the function    @xmath75\right\ } \sin\left(\frac{\pi b\beta}{\alpha } \sin\phi\right)\right|^{-1 } \nonumber\\ & \phantom{= } \null + \tfrac{1}{9 } v_0 \left|\sin\left\ { \frac{\pi b}{\alpha } \left [ \tfrac12 \sqrt{3 }
\cos\phi + \left(\gamma- \tfrac12 \beta\right ) \sin\phi \right]\right\ } \sin\left(- \frac{\pi b \beta}{\alpha } \sin\phi\right)\right|^{-1}~.\end{aligned}\ ] ]    this function gives an upper bound for the absolute value of the corrugation due to a particular term .
it can be used to estimate the typical corrugation of any triangular or trapezoidal crystal . at particular combinations of orientations @xmath76 and lattice parameter ratio @xmath74 , one or both of the denominators vanishes and @xmath77 diverges .
the numerator in eq .
( [ eq : pottrap2 ] ) vanishes as well , so that @xmath73 itsself does not diverge . as the arguments of the trigonometric functions in the numerator contain factors of @xmath43 and @xmath44 , the resulting contribution to @xmath73
is an order of @xmath43 or @xmath44 higher than at other orientations .
there is some arbitraryness to the choice of numerator and denominator , as both can be multiplied by any function .
this can not however lead to additional factors of @xmath43 . for determining the order of magnitude of the corrugation at any orientation , @xmath77 therefore suffices .     the angle - dependent corrugation prefactor , eq .
( [ eq : pottrapcor ] ) , for a perfectly triangular nanocrystals on a hexagonal substrate with @xmath78 ( solid line ) , along with @xmath79 ( dotted line ) , which does not depend on the diameter @xmath43 .
the ratio @xmath80 .
the denominators control the typical size of the corrugation .
whenever the denominator nearly diverges near the orientations given by eqs .
( [ eq : denominator1 ] ) and  ( [ eq : denominator2 ] ) , a peak appears .
there , the interface is ( partially ) commensurate .
, width=325 ]     a logarithmic plot of a peak in the corrugation for a perfectly triangular rigid nanocrystals on a hexagonal substrate with @xmath81 ( thin solid line ) and @xmath82 ( thick solid line ) , along with @xmath79 ( dotted line ) .
the denominators control the typical size of the corrugation . in this case , the corrugation diverges with @xmath83 , but , very close to @xmath76 , it remains finite , because the numerator vanishes as well . , width=325 ]
the ( in)commensurate configurations can be distinguished and classified by the exponent @xmath84 with which the corrugation scales with the contact diameter .
this order of the commensurability can be obtained from the divergence of the size - independent prefactor in the corrugation , which goes as @xmath85 .
i.e. , it is equal to the number of derivatives ( starting with the zeroth ) with respect to @xmath34 of the denominators in eq .
( [ eq : pottrap2 ] ) that are equal to zero .
the scaling of the friction with the contact area can be directly related to the order of commensurability through @xmath86 for the commensurate orientations generated by the first - order correction in eq .
( [ eq : firstordercorrection ] ) that do not coincide with commensurate orientations of the leading order , the friction is a factor of order @xmath87 smaller .    for most parameter combinations and orientations
the contact is incommensurate with @xmath88 , the denominator does not vanish at all , and both denominator and numerator are of order @xmath89 .
the corrugation of the nanocrystal on the substrate does not scale with the size of the crystal .
the sliding friction consequently is of the order @xmath90 .
it should be noted , however , that without a load force , @xmath15 for gold on graphite is comparable to @xmath91 at room temperature .
temperature effects on the friction can therefore play a role as well in this regime , reducing the friction further .
commensurability order @xmath92 can occur when one of the two @xmath93 functions in the denominator in eq .
( [ eq : pottrap2 ] ) vanishes for a particular orientation @xmath76 . for a rigid crystal , this occurs when one of the following conditions is met : @xmath94 & = p~,\\ \label{eq : denominator2 } \frac{b\beta}{\alpha } \sin\phi_{\mathrm{c . } } & = p'~,\end{aligned}\ ] ] with @xmath95 .
if @xmath74 is sufficiently large , i.  e. large inter - atomic distance in the nanocrystal , solutions exist also for @xmath96 .
this implies that there exist nontrivial commensurate - like orientations with high corrugation . at these orientations ,
there is no obvious lining up of symmetry axes , or other source of commensurability , as there is for the trivial solutions with @xmath97 , but the corrugation and friction nevertheless behave in a similar way .
this is illustrated in fig .
[ fig : nontrivialangles ] , where the corrugation is plotted for several values of @xmath74 . two commensurate combinations of lattices are also shown , one trivial , and one nontrivial . in the nontrivial one , no lining up is immediately visible , yet the potential energy still scales with the diameter , as it does for the trivial commensurate configuration . for bulk
au and graphite , @xmath41 , @xmath98 , @xmath99 , and thus there are commensurate orientations @xmath100 .     the corrugation ( a ) as a function of the orientation for crystals with diameter @xmath101 and various @xmath74 and two commensurate orientations with @xmath92 for triangular crystals with @xmath102 and with ( b ) trivial @xmath103 and ( c ) nontrivial @xmath104 .
though the nontrivial commensurate configuration shown in ( b ) does not have any obvious lining up of symmetry axes , it nevertheless corresponds to a strong peak in the corrugation in ( c ) . for larger @xmath74 , more nontrivial solutions of eqs .
( [ eq : denominator1 ] ) and  ( [ eq : denominator2 ] ) exist and hence more orientations with large corrugations . , title="fig:",width=325 ]   the corrugation ( a ) as a function of the orientation for crystals with diameter @xmath101 and various @xmath74 and two commensurate orientations with @xmath92 for triangular crystals with @xmath102 and with ( b ) trivial @xmath103 and ( c ) nontrivial @xmath104 .
though the nontrivial commensurate configuration shown in ( b ) does not have any obvious lining up of symmetry axes , it nevertheless corresponds to a strong peak in the corrugation in ( c ) . for larger @xmath74 , more nontrivial solutions of eqs .
( [ eq : denominator1 ] ) and  ( [ eq : denominator2 ] ) exist and hence more orientations with large corrugations . , title="fig:",width=309 ]    the first - order correction for lattice deformations [ eq .  ( [ eq : firstordercorrection ] ) ] produces additional commensurate orientations , which satisfy one of the conditions @xmath105 & = \tfrac12 p~,\\ \label{eq : denominator4 } \frac{b\beta}{\alpha } \sin\phi_{\mathrm{c . } } & = \tfrac12 p'~,\\ \frac{b}{\alpha } \left [ \tfrac12 \sqrt{3 } \cos(\phi_{\mathrm{c.}}-\tfrac16 \pi ) + \gamma \sin(\phi_{\mathrm{c.}}-\tfrac16 \pi ) \right.&\nonumber\\\left.\null \label{eq : denominator5 } \pm \tfrac12 \beta \sin(\phi_{\mathrm{c.}}-\tfrac16
\pi ) \right ] & = \tfrac13\sqrt{3 } p~,\\ \frac{b\beta}{\alpha } \sin(\phi_{\mathrm{c.}}-\tfrac16\pi ) & = \tfrac13\sqrt{3 } p'~. \label{eq : denominator6}\end{aligned}\ ] ] these conditions contain additional factors of @xmath106 and @xmath107 when compared to eqs .
( [ eq : denominator1 ] ) and  ( [ eq : denominator2 ] ) , due to the appearance of such factors in eq .
( [ eq : firstordercorrection ] ) .
equations  ( [ eq : denominator3 ] ) through  ( [ eq : denominator6 ] ) yield a set of additional commensurate orientations that are different from the ones produced by a rigid crystal . in fig .
[ fig : corrections ] , the corrugation in the total potential energy , consisting of the rigid - crystal term and the first - order correction for lattice deformations , is shown for a very large perfectly triangular crystal .
peaks corresponding to both the nontrivial commensurate configurations and the additional commensurate configurations associated with lattice deformations are clearly visible .     the the corrugation ( a ) as a function of the orientation for non - rigid triangular crystals with diameter @xmath108 and @xmath109 ( the typical value for au ) , for triangular crystals with @xmath110 , and lattices for commensurate orientations with @xmath92 , for ( b ) a rigid crystal @xmath103 and ( c ) a non - rigid crystal at @xmath111 . the ratio of lattices parameters @xmath110 is incommensurate .
the short lines in ( b ) indicate the direction of the displacement of the atoms as a result of the substrate potential , and their length is proportional to its magnitude . the peaks marked `` rigid crystal '' appear already when the crystal is assumed to be rigid , and satisfy eqs .
( [ eq : denominator1 ] ) and  ( [ eq : denominator2 ] ) . the lower peaks that appear due to lattice deformations
are marked with `` first order '' and are the solutions of eqs .
( [ eq : denominator3 ] )  ( [ eq : denominator6 ] ) . , title="fig:",width=325 ]   the the corrugation ( a ) as a function of the orientation for non - rigid triangular crystals with diameter @xmath108 and @xmath109 ( the typical value for au ) , for triangular crystals with @xmath110 , and lattices for commensurate orientations with @xmath92 , for ( b ) a rigid crystal @xmath103 and ( c ) a non - rigid crystal at @xmath111 .
the ratio of lattices parameters @xmath110 is incommensurate .
the short lines in ( b ) indicate the direction of the displacement of the atoms as a result of the substrate potential , and their length is proportional to its magnitude .
the peaks marked `` rigid crystal '' appear already when the crystal is assumed to be rigid , and satisfy eqs .
( [ eq : denominator1 ] ) and  ( [ eq : denominator2 ] ) .
the lower peaks that appear due to lattice deformations are marked with `` first order '' and are the solutions of eqs .
( [ eq : denominator3 ] )  ( [ eq : denominator6 ] ) . ,
title="fig:",width=309 ]    finally , @xmath112 can occur only if the two @xmath93 functions in the denominators have coinciding zeros .
this happens only for special ratios of the lattice parameters and corresponds to commensurability in two parameters .
the conditions of eqs .
( [ eq : denominator1 ] ) and  ( [ eq : denominator2 ] ) for the orientation can only be satisfied simultaneously if the lattice parameters are related by @xmath113 ^ 2~ , \label{eq : denominatorx1}\end{aligned}\ ] ] with @xmath95 . from the first - order correction in eq .  ( [ eq : firstordercorrection ] ) additional @xmath112 commensurate combinations of parameters are found , @xmath114 ^ 2~,\\ 3 \left(\frac{b}{\alpha}\right)^2 = \left(\frac{p'}{\beta}\right)^2 + \frac43 \left[p - p'\left(\frac\gamma\beta \pm \frac12 \right)\right]^2~. \label{eq : denominatorx3}\end{aligned}\ ] ] equation  ( [ eq : denominatorx2 ] ) corresponds to commensurate configurations at angles satisfying eqs .
( [ eq : denominator3 ] ) and  ( [ eq : denominator4 ] ) , while eq .
( [ eq : denominatorx3 ] ) corresponds to commensurate configurations at angles satisfying eqs .
( [ eq : denominator5 ] ) and  ( [ eq : denominator6 ] ) .
a similar structure appears also in the expressions for rectangular nanocrystals on rectangular lattices shown in the appendix and also treated in ref .
this is due to the fact that any sum over atoms arranged in a regular way in nanocrystal , subjected to a periodic potential , can be written in terms of double sums of the form of eq .
( [ eq : sinratio ] ) .
hence , two @xmath93 functions always appear in the denominator , and similar conditions exist for which the denominator vanishes .
the system studied in ref .
@xcite , w on naf , however does not have a sufficiently high ratio of lattice parameters to produce nontrivial commensurate orientations . consequently , no nontrivial commensurate orientations were observed there .      in experiments of friction of nanocrystals ,
the crystals are either pulled along a particular direction  @xcite , or pushed  @xcite .
when crystals are pulled , the direction of their motion is enforced externally .
when the crystals are pushed , they may move in a direction not precisely along to the force , following some preferred sliding direction .
a comparison between eqs .
( [ eq : pottrap1 ] ) and  ( [ eq : substratepotential ] ) shows that the potential energy of a nanocrystal depends on the position in a similar way to that of a single atom .
consequently , the preferred sliding directions are the same as for a single atom , along vectors othogonal to @xmath63 .
the nontrivial commensurate orienations do not lead to any nontrivial preferred sliding directions .
in addition , the activation energy is always equal to the same fraction of the corrugation .
however , the first - order corrections produce a potential form which is rotated by @xmath115 . for some parameter combinations
, this may lead to additional preferred sliding directions parallel to @xmath63 .
it is interesting to compare the model used here and the results to what is known about the frenkel - kontorova ( fk ) model  @xcite in one and two dimensions .
the fk model consists of a chain or sheet of particles , coupled to one - another and a ( quasi-)periodic substrate .
the coupling between the particles may be weak or strong compared to the substrate , leading to different types of behavior .
when the substrate is weak by comparison , for incommensurate lattice parameters , the interface can slide without friction . in the one - dimensional fk model ,
commensurate configurations occur for rational ratios of the lattice parameters  @xcite .
for the nearly rigid case , in the fk model , there is a cascade of commensurate configurations with decreasing friction .    for weak coupling between monomers in the chain
, there is a breaking of analyticity  @xcite , which results in any ratio of lattice parameters yielding nonvanishing friction for infinite contacts .
other effects also start playing a role that do not affect stiff crystals , for instance nucleation of incommensurate domains within a large sheet  @xcite ) .    in the system
considered here a similar cascade would appear , if more higher order corrections for lattice deformations were included .
more commensurate orientations with increasingly small friction would be found .
each subsequent term in the potential energy would however be multiplied by a small prefactor , and so each additional commensurate orientation would correspond to an increasingly small corrugation .
the one - dimensional fk model has only two length scales , and therefore only one parameter which controls commensurability .
the system studied in this work has two parameters which control commensurability , similar to the frenkel - kontorova ( fk ) model in two dimensions , and the quasi - periodic one - dimensional fk model  @xcite . unlike the model described here
, the quasiperiodic fk model can not have coexistence of commensurate states with different order of commensurability .
the two - dimensional fk model , however , is already too complex to be analytically tractable .
hence , most results dealing with extended two - dimensional deformable interfaces are numerical ( see , for instance , refs  @xcite ) .
recently , though , experiments were performed that model atoms as colloids in suspensions , and reproduce the fk model closely  @xcite .
in real - world , experimental conditions , crystaline particles are almost never perfectly triangular , or even trapezoidal . at the very least , even if the crystal lattice is in tact , some atoms may be missing and the corners may be rounded .
any shape of particle can be written as a sum ( including negative terms , if necessary ) of regularly shaped particles .
the number of terms in this sum can be estimated from the number of step edges @xmath116 on the circumference of the contact layer .
each term contributes an amount to the total potential energy of the order of @xmath117 .
because the total potential energy at a particular orientation varies rapidly with the size of the particle [ see eq .
( [ eq : pottrap2 ] ) ] , one may assume a random phase for every term .
the total corrugation for an irregular particle is therefore of the order of @xmath118 , but can never grow faster than @xmath119 . if the shape of the particle is simply scaled , then the number of step edges increases with @xmath120 .
such crystals therefore have friction scaling similarly to eq .
( [ eq : scaling ] ) with , @xmath121 i.  e.  with exponents @xmath122 , and @xmath89 .
the effects of lattice distortions at the edge of the contact layer can be included by summing up in a similar way over a small number of contact layers with different sizes to account for the different interactions at the edge of the crystal .
therefore , the scaling with the crystal size and the ( in)commensurate configurations is not affected by the edges .
it should also be noted that , regardless of their shape , the friction of amorphous crystals always scales @xmath123  @xcite .
the three different types of ( in)commensurability , quantified with the values of @xmath124 , can occur at different orientations in a single interface already for rigid crystals . in fig .
[ fig : multiplicity ] , such a case is shown .
the trivial commensurate orientations have @xmath92 , corrgation scaling with @xmath43 , and are the result of a single vanishing denominator in eq .
( [ eq : pottrap2 ] ) .
the nontrivial commensurate orientations are the result of both denominators vanishing at the same orientation , and thus correspond @xmath112 and corrugation scaling with @xmath125 .
the corresponding orientations around @xmath126 are particularly interesting , as these are usually assumed to be strongly incommensurate , and are indeed so for most combinations of lattice parameters .     the corrugation ( a ) as a function of the orientation for @xmath127 and @xmath102 for a perfectly triangular nanocrystal with @xmath128 , and two commensurate orientations for @xmath127 , at ( b ) @xmath129 and ( c ) @xmath130 .
the three different friction branches occur in one interface here . for the case of @xmath131 ,
the corrugation is of the order @xmath125 , while for @xmath103 it is of order @xmath43 , and at most other orientations the typical corrugation is not dependent on  @xmath43 .
, title="fig:",width=325 ]   the corrugation ( a ) as a function of the orientation for @xmath127 and @xmath102 for a perfectly triangular nanocrystal with @xmath128 , and two commensurate orientations for @xmath127 , at ( b ) @xmath129 and ( c ) @xmath130 .
the three different friction branches occur in one interface here . for the case of @xmath131 ,
the corrugation is of the order @xmath125 , while for @xmath103 it is of order @xmath43 , and at most other orientations the typical corrugation is not dependent on  @xmath43 .
, title="fig:",width=309 ]    for some lattice parameters , the rigid crystal has only two friction branches , but a third branch appears when lattice deformations are taken into account .
this is shown for a perfectly triangular non - rigid nanocrystal with @xmath132 in fig .
[ fig : correctionsau ] .
the first - order peaks appear at the solutions of eqs .
( [ eq : denominator3 ] )  ( [ eq : denominator6 ] ) , and are generated by both rigid - crystal contributions and first - order corrections . however , for this combination of lattice parameters , eqs .
( [ eq : denominator5 ] ) and  ( [ eq : denominator6 ] ) are satisfied for @xmath133 . as a result , @xmath112 is possible for the first - order correction terms , but not for the zeroth order . consequently , the small prefactor for the stiffness of the crystal is multiplied by a factor that can be as large as @xmath43 .    for sufficiently small crystals with nearly commensurate lattice parameters , the commensurate orientation behaves as if @xmath112 , and its friction grows as @xmath125 . for gold on graphite , with @xmath134 close to @xmath135 , the corrections can easily become as large as the leading order terms .
it should however be noted that , as this effect is very sensitive to the lattice parameters , the surface reconstruction , however minimal , can still enhance or decrease the effect strongly . the surface construction of au(111 ) on graphite
may indeed itsself be affected by the near match in lattice parameters , possibly leading to an even closer match of the reconstructed surface .
the fk model can be used to understand surface recontructions  @xcite , and the calculations of this paper can similarly be applied to get a handle on surface reconstructions in an interface .
the corrugation ( a ) in the potential energy as a function of the orientation and commensurate lattices ( b ) at @xmath129 , ( c ) at @xmath130 for a perfectly triangular nanocrystal with @xmath108 that is not rigid with @xmath109 and @xmath132 , which is a commensurate configuration for the first - order corrections , this plot is similar to fig .
[ fig : corrections ] . the additional peaks that appear because the crystal is not rigid are the solutions of eqs .
( [ eq : denominator3 ] )  ( [ eq : denominator6 ] ) .
the orientation shown in ( c ) is usually strongly incommensurate , even though the symmetry axes area clearly lined up , because the potential energy of the interface is not strongly corrugated . due to the matching lattice parameters , the first - order corrections contribute @xmath136 to the corrugation . for rigid crystals
this relative orientation is incommensurate for almost all values of @xmath74 . even for @xmath135 , which looks commensurate
, a rigid crystal does not have a corrugation that increases with the size .
, title="fig:",width=325 ]   the corrugation ( a ) in the potential energy as a function of the orientation and commensurate lattices ( b ) at @xmath129 , ( c ) at @xmath130 for a perfectly triangular nanocrystal with @xmath108 that is not rigid with @xmath109 and @xmath132 , which is a commensurate configuration for the first - order corrections , this plot is similar to fig .
[ fig : corrections ] .
the additional peaks that appear because the crystal is not rigid are the solutions of eqs .
( [ eq : denominator3 ] )  ( [ eq : denominator6 ] ) .
the orientation shown in ( c ) is usually strongly incommensurate , even though the symmetry axes area clearly lined up , because the potential energy of the interface is not strongly corrugated . due to the matching lattice parameters , the first - order corrections contribute @xmath136 to the corrugation . for rigid crystals this relative orientation is incommensurate for almost all values of @xmath74 . even for @xmath135 , which looks commensurate
, a rigid crystal does not have a corrugation that increases with the size .
, title="fig:",width=309 ]      the order of commensuratebility , as it is defined in this work , is a general quantity that can be used to characterize commensurability of crystalline interfaces .
the @xmath93 functions in the denominator in eq .
( [ eq : pottrap2 ] ) originate from summing over a line of regularly spaced atoms , which must always produce factors of the type shown in eq .
( [ eq : sinratio ] ) .
such sums occur in the potential energy of any rigid regular crystalline contact layer .
thus , the analysis described here in terms of its derivatives can be used to identify and classify commensurate configurations for any regular nanocrystal on any regular substrate . in the appendix ,
the expression for the potential energy of a rectangular crystal on a rectangular substrate is worked out as an example .    for a two - dimensional contact area ,
it is not possible to have @xmath84 larger than 2 , because the total potential energy can never grow faster than linearly with the contact area .
correspondingly , in one dimension , @xmath84 can not exceed unity .
as illustrated by the case of gold on graphite , configurations close to commensurate still show some of the commensurate behavior for finite size crystals .
it is therefore worthwhile to consider a quantity similar to @xmath84 , a finite - size order of commensurability , @xmath137}{\ln { n}}~.\end{aligned}\ ] ] for simplicity , or if @xmath138 is not well - known , a simple periodic sinusoidal function can be used in place of @xmath138 . in the limit of @xmath139 , @xmath140 converges to the order of commensurability @xmath84 .
as it is impossible to define lattice orientation for a single atom , it diverges for @xmath141 . for finite - size crystals , it can be used as a measure of the commensurability of the contact . in ref .  @xcite another measure for commensurability
was defined for determining to what extent weakly - coupled atoms line up with the substrate in an incommensurate contact .
this measure , however , is not suitable for determining the commensurability of two ( nearly ) rigid lattices .
in this work , the ( in)commensurability of finite - size crystalline interfaces was investigated , particularly for crystals with triangular symmetry on triangular or hexagonal substrates such as gold on graphite .
the crystals were assumed to be nearly rigid , with the stiffness being used as an expansion parameter .
a simple method was developed for determining the commensurate configurations , through quantifying the commensurability by the scaling of the potential corrugation of the crystal with the diameter of the contact area at constant orientation and lattice parameters .
this method can be applied not only to the two geometries discussed in this paper ( triangular / hexagonal and rectangular / rectangular ) , but to any regular lattice on a regular substrate .
it was found that , due to the two - dimensional nature of the contact layer , and associated two parameters for commensurability ( ratio of lattice parameters and relative orientation ) , two different types of commensurate contacts can exist for the same materials . if the relative orientation is commensurate , but the lattice parameters are not , the potential corrugation scales with the diameter of the contact area , while if both are commensurate , it scales with the contact area .
consequently , for some combinations of materials three different friction branches can appear for the same crystal shape and materials , scaling with the contact area as @xmath7 , and @xmath142 for perfect crystals [ see eq .
( [ eq : scaling ] ) ] . when the crystals have irregular shapes , the corrugation and friction scale with @xmath143 [ see eq .
( [ eq : scalingimperfect ] ) ] .
nontrivial commensurate configurations were identified in particular for the triangular on triangular on hexagonal geometry .
these commensurate orientations exist if the lattice parameter the nanocrystal is sufficiently high compared to the lattice parameter of the substrate .
they are also a direct result of the two - dimensional nature of the contact layer .
it was found that lattice deformations produce a number of additional commensurate configurations that would be incommensurate for a rigid crystal .
interestingly , the prototype geometry of gold on graphite is very close to the parameters which produce such a friction multiplicity
. it might therefore be possible to detect all three friction branches with different commensurability in friction experiments of gold nanocrystal on graphite , as well as nontrivial commensurate angles .
this work was financially supported by a veni grant of netherlands organisation for scientific research ( nwo ) and by an unga forskare grant from the swedish research council .
the author would also like to thank a.  fasolino , d.  dietzel , and a.  schirmeisen for discussions , and o.  manyuhina for carefully reading the manuscript .
31 natexlab#1#1bibnamefont # 1#1bibfnamefont # 1#1citenamefont # 1#1url # 1`#1`urlprefix[2]#2 [ 2][]#2    , * * , ( ) .
, in _ _ , edited by ( , ) , vol .   of _ _ , pp .
, * * , ( ) .    , , , * * , ( ) .
, , , , , , * * , ( ) .    , , , , , * * , ( ) .    , , , , , , * * , ( ) .    , , , , , , , , , * * , ( ) .
, , , * * , ( ) .    , * * , ( ) .    , , , , ( ) , .
, , , , * * , ( ) .    , * * , ( ) .    , , , * * , ( ) .    , , , , * * , ( ) .    , * * , ( ) .    , * * , ( ) .    , * * , ( ) .    , , , , , * * , ( ) .    , * * , ( ) .    , , , , , * * , ( ) .    , * * , ( ) .    , , , , * * , ( ) .    ,
* * , ( ) .
, , , , * * , ( ) .    , * * , ( ) .    , , , , * * , ( ) .    , ,
, * * , ( ) .    , , , * * , ( ) .
, , , * * , ( ) .    , * * , ( ) .
, lattice parameters @xmath144 , contact layer size @xmath145 , and the potential energy of a contact layer atom on the substrate .
if a contact layer atom lies on top of a substrate atom , its potential energy is @xmath146 . if it lies directly between an atom and its nearest neighbor in the @xmath13 or @xmath147 direction , it has potential energy @xmath148 or @xmath149 respectively .
if it lies in the center of a rectangle , at equal distance from four substrate atoms , without loss of generality , we may set the potential energy to 0 .
the origin of the coordinate system is chosen to lie on top of a substrate atom .
[ fig : rectangular ] , width=294 ]    the identification and classification of commensurate(-like ) configurations described in this paper can easily be repeated for the general rectangular lattices , as described for instance in ref .
a general potential energy of a single atom , @xmath150 is given by @xmath151 where the geometry and associated parameters are defined in fig .
[ fig : rectangular ] .
@xmath152 \nonumber\\ & \phantom{= nnnn } \times \frac { \sin\left[\pi ( m+1 ) \left(\frac{b_1}{a_1}\cos{\phi } + ( -1)^j \frac{b_1}{a_2}\sin{\phi}\right)\right ] \sin\left[\pi ( n+1 ) \left(\frac{b_2}{a_1}\sin{\phi } - ( -1)^j \frac{b_2}{a_2}\cos{\phi}\right)\right ] } { \sin\left[\pi \left(\frac{b_1}{a_1}\cos{\phi } + ( -1)^j \frac{b_1}{a_2}\sin{\phi}\right)\right ] \sin\left[\pi \left(\frac{b_2}{a_1}\sin{\phi } - ( -1)^j \frac{b_2}{a_2}\cos{\phi}\right)\right ] } ~,\end{aligned}\ ] ]    where @xmath153 originate from the second , third , and fourth term on the right - hand side of eq .
( [ eq : va ] ) .
similarly , the higher - order corrections for lattice deformations can also be obtained along similar lines as described in sec .
[ sec : firstordercorrections ] .
this yields terms with higher perdiodicity , i.e. , a term obtained from the derivative of @xmath154 which looks very similar , except that @xmath155 is replaced by @xmath156 .
similarly , there is a term derived from @xmath157 with @xmath158 replaced by @xmath159 , and two terms from @xmath160 with both @xmath155 and @xmath158 replaced .    from these expressions
it is easy to obtain the conditions for vanishing denominators , and thus for @xmath124 .
the three terms @xmath153 each produce a different set of conditions for commensurability at different angles and ratios of lattice parameters . for @xmath154 ,
one finds that @xmath92 commensurability occurs if one of the following conditions is met , @xmath161 for some @xmath95 .
both of these conditions can be met simulataneously , leading to @xmath112 commensurability , if @xmath162          as each of @xmath153 has a different prefactor that depends on the energy barriers of the potential , the different terms do not , in general , cancel each other out .
it is interesting to note that often for square lattices it is assumed that @xmath166 . in this case ,
@xmath160 vanishes and the commensurate configurations associated with it disappear , leading to lower friction .