the influence of aperiodic ( dis-)order on the thermodynamic properties and the critical behaviour of spin systems has been an active research area , particularly since the discovery of quasicrystals in 1984 . in view of the many articles existing
, we can not give proper credit to all contributors ; instead , we would like to refer the reader to a recent review @xcite which contains a rather complete bibliography on this subject .
for ferromagnetic ising systems , heuristic scaling arguments , put forward by luck @xcite , lead to qualitative predictions on the relevance of disorder implemented through an aperiodic variation of the coupling constants which go beyond the original relevance - irrelevance criterion of harris @xcite .
for the ising quantum chain ( iqc ) , or the equivalent 2@xmath2 classical ising model @xcite with layered disorder , this relevance criterion reads as follows : the critical behaviour is of onsager - type as long as the fluctuations in the sequence of coupling constants are bounded , whereas for unbounded fluctuations the critical singularities resemble those of the randomly disordered case .
this behaviour had first been conjectured by tracy @xcite on ths basis of his results for two different examples of three - letter substitution rules .
of particular interest is the marginal case with logarithmically ( in the system size ) diverging fluctuations . here ,
the critical exponents are predicted to depend continuously on the coupling constants @xcite .
meanwhile , some exact results were obtained for substitution chains with marginal fluctuations of the couplings . for a number of substitution chains ,
the surface magnetisation critical exponent could be calculated @xcite ; and , recently , the correlation length critical exponent was derived for two special substitution rules using an exact renormalisation approach @xcite .
all these results are in accordance with luck s predictions .
the real space renormalisation method we use in this article is based on the decimation process introduced in ref .
our setup is quite different , though ; and in particular we are going to show how _ arbitrary _ substitution rules can be analysed in terms of @xmath3-matrices and their star - product formalism . in this way , the dependence of the critical behaviour on the nature of the fluctuations can be derived explicitly .
the analysis is carried out in detail for two - letter substitution chains , with the result that luck s relevance criterion emerges from the renormalisation relations in a natural way .
an extended renormalisation scheme for the most general case of @xmath4-letter substitution rules is presented in an appendix and discussed for certain examples . for infinite classes of substitution chains with marginal fluctuations , including the two - letter chains , an analytic expression for the correlation length critical exponent
is given .
this article is organised as follows . in sec .  2 , we give a brief outline of the model , including some introductory remarks on substitution rules and on properties of the corresponding substitution sequences .
then , in sec .  3
, we introduce the renormalisation transformation for a restricted class of @xmath4-letter substitution chains , comprising all two - letter substitution chains .
the dependence of the renormalisation flow on the fluctuations of the couplings is shown and the consequences are discussed in detail . subsequently , in sec .  4 , we summarise the implications of our results for the critical properties of the iqc for the three different types of fluctuations in the sequence of the coupling constants . our conclusions and an outlook on future developments are presented in sec .
the paper contains two appendices related to the content of sec .  3 :
in appendix  a , the solution of the eigenvalue problem for the case of a two - letter substitution sequence is derived ; and finally , generalising the discussion of sec .  3 , the renormalisation transformation for arbitrary @xmath4-letter substitution rules is explained in appendix  b.
the iqc in a transverse magnetic field is defined by the hamiltonian @xmath5 acting on the tensor - product space @xmath6 . here ,
the @xmath7 denote the transverse fields at the sites .
the @xmath8 are site - dependent coupling constants , and the operators @xmath9 denote pauli s matrices acting on the @xmath10th site . in what follows ,
we consider closed chains with _ periodic _ boundary conditions which are defined by @xmath11 ; various other kinds of boundary conditions can be chosen without changing the key results .    for a general set of coupling constants @xmath8 and fields @xmath7 ,
the iqc can be written as a free - fermion model via a jordan - wigner transformation and can then be diagonalised by a bogoljubov - valatin transformation , resulting in @xcite @xmath12 where @xmath13 , @xmath14 are @xmath15 fermionic creation and annihilation operators and @xmath16 is some constant .
the dimensionless excitation energies @xmath17 ( which can be chosen to be positive and ordered , i.e. , @xmath18 , while properly adjusting the constant @xmath16 ) satisfy the linear difference equations @xmath19 normally , one proceeds by eliminating either @xmath20 or @xmath21 . in this way
, the problem of diagonalising the @xmath22 matrix @xmath23 can be reduced to the diagonalisation of an @xmath24 matrix . from the set of @xmath15 eigenvalues ( fermion frequencies ) @xmath25 of the latter ,
the total spectrum of the complete hamiltonian @xmath23 is recovered by considering the sums of the elements of any of its @xmath26 subsets .
for periodic boundary conditions , the situation is in fact a little more complicated ; in order to avoid the appearance of a non - local number operator in the fermionic hamiltonian one has to consider so - called _ mixed sector _ hamiltonians .
these are obtained by combining the _ even _ sector ( with respect to the operator @xmath27 that commutes with @xmath28 ) of the hamiltonian ( [ eq : qc ] ) with _ periodic _ boundary conditions with the _ odd _ sector of the same hamiltonian with _ antiperiodic _ boundary conditions ( defined by @xmath29 ) , and vice versa @xcite .
the iqc is equivalent to the 2@xmath2 classical ising model with a layered modulation of the interactions @xcite .
the correspondence shows up through an extremely anisotropic limit ( also called the @xmath30-continuum limit @xcite ) of the coupling constants of the latter , which is believed not to alter the critical behaviour .
explicitly , one establishes the following relations : the _ temperature _ @xmath31 of the classical model corresponds to the _ transverse field _
@xmath32 of the quantum chain , the _ free energy
_ translates into the _ ground - state energy _ and the counterpart of the _ correlation length @xmath33 parallel _ to the layers is the _ inverse mass gap _
@xmath34 of the iqc . in our setting ,
the mass gap is just given by the modulus of the smallest fermionic excitation @xmath35 of ( [ ffh ] ) . according to finite - size scaling
, we thus expect a scaling behaviour of @xmath35 at the critical point as @xmath36 provided the correlation length scales with a power law at all .
hereby , the scaling exponent @xmath37 translates into the correlation length critical exponent @xmath38 parallel to the layers in terminology of the 2@xmath2 statistical system , while @xmath39 retains its unperturbed onsager value @xmath40 .
the hamiltonian ( [ eq : qc ] ) is _ critical _ ( in the sense that the energy gap in the excitation spectrum vanishes , i.e. , @xmath41 ) in the thermodynamic limit if the _ geometric mean _ ( rather than the arithmetic mean ) of the couplings , normalised by the variables of the transverse field , is equal to one , i.e. , @xmath42 compare ref .
the repetitive structure of the substitution sequences we are going to consider guarantees that the mean is well - defined @xcite .
let us now assume that the coupling constants @xmath43 ( @xmath44 ) are chosen from a set of @xmath4 possible values @xmath45 ( @xmath46 ) , the value @xmath45 occuring with a frequency @xmath47 ( where @xmath48 ) . for simplicity ,
we take the field to be constant , i.e. , @xmath49 .
then , condition ( [ eq : critcond ] ) can be fulfilled by choosing @xmath4 positive real numbers @xmath50 with @xmath51 , and setting the couplings to @xmath52 where @xmath53 . unless stated otherwise , this parametrisation of the critical surface will be used in the sequel .
let us now consider coupling constants @xmath8 drawn from a set of values @xmath54 , where the label @xmath55 runs over the letters of a ( finite ) @xmath4-letter alphabet @xmath56 .
we choose the couplings @xmath8 according to sequences generated by iterated application of a substitution rule @xmath57 on @xmath56 . for the introduction of the concept and some notation , consider a simple example on a two - letter alphabet @xmath58 : @xmath59 we start with an initial word , e.g. @xmath60 , and define @xmath61 .
@xmath62 denotes the corresponding substitution matrix ( whose elements count the number of letters @xmath63 and @xmath64 in the words @xmath65 and @xmath66 ) , which has the eigenvalues @xmath67 .
the entries of the , statistically normalised , ( right ) eigenvector @xmath68 , corresponding to the perron - frobenius eigenvalue @xmath69 , give the frequencies of the letters in the limit word @xmath70 . in order to ensure that , starting from any seed , the letter frequencies converge to unique , positive values in the limit chain , we restrict ourselves to _ primitive _ substitution rules @xmath71 throughout this paper .
for the corresponding substitution matrices @xmath72 , this means that the entries of @xmath73 are _ strictly positive _ for any integer @xmath74 ( and a suitably chosen finite @xmath75 ) .
this is obviously true for the examples ( [ eq : sub ] ) .
associated to a word @xmath76 of length ( i.e. , number of letters ) @xmath77 , one obtains finite ( closed ) quantum chains with @xmath15 sites by setting @xmath78 whenever the @xmath10th letter of the word @xmath76 is an @xmath63 , and @xmath79 otherwise , the final letter of @xmath76 determining the coupling between the last and the first spin in the chain . while @xmath80 , the largest eigenvalue of the substitution matrix , describes the asymptotic scaling of the chain length with the number of iterations @xmath81 , @xmath82 , the _ fluctuation _ behaviour of the exchange couplings is connected to the second - largest eigenvalue @xmath83 ( where it is the _ modulus _ of the eigenvalue that matters , hence by `` second - largest eigenvalue '' we always refer to the eigenvalue with the second - largest absolute value ) . in order to quantify the fluctuations ,
let us define @xmath84 where @xmath85 is the average coupling in the limit chain .
fluctuations stay _ bounded _ as long as @xmath86 , but for @xmath87 they grow with iterated substitutions as a power of the chain length , @xmath88 , governed by the so - called _ wandering exponent _
@xmath89 that is given by @xmath90 in the limiting case @xmath91 , @xmath92 stays constant when we prolong the chain by complete iteration steps .
generically , however , finite fluctuations within the substitutes @xmath93 of the letters @xmath94 will increase the total fluctuation in the vicinity of these points by a certain amount @xmath95 : @xmath96 for some word @xmath97 with @xmath98 and fixed @xmath99 . by repeated use of this argument
, we thus expect fluctuations that _ diverge logarithmically _ with the length of the chain ( since the length grows exponentially with the number of iterations ) , see also ref .
@xcite . on the other hand
, the same argument also shows that , in this case , the fluctuation behaviour is not completely determined by the _ substitution matrix _ , but may also depend on the details of the actual _ substitution rule _
@xcite ; and for special choices the fluctuations may still stay bounded even though the second - largest eigenvalue is @xmath91 .
coming back to our examples ( [ eq : sub ] ) , let us concentrate on the three cases @xmath100 ( which yields the famous _ fibonacci sequence _ ) , @xmath101 ( the _ period - doubling sequence _ ) , and @xmath102 ( the _ binary non - pisot sequence _ ) .
as the second - largest eigenvalues of the substitution matrices @xmath62 fulfill @xmath103 , @xmath104 , and @xmath105 , this includes an example of bounded fluctuations ( @xmath100 ) , a marginal case with logarithmically divergent fluctuations ( @xmath101 ) and an example of strong fluctuations ( @xmath102 ) , hence all three interesting cases are illustrated @xcite .
another substitution chain that is referred to several times in the text is the celebrated _ thue - morse sequence _ @xmath106
rather than reducing the set of linear equations ( [ fermeq ] ) to @xmath15 dimensions , we proceed ( as does ref .
@xcite ) with a direct treatment as a @xmath107-dimensional eigenvalue problem of the real symmetric matrix @xmath108{$h_{\scriptscriptstyle n-1}^{}$ } & 0&\makebox[0pt]{$\varepsilon_{\scriptscriptstyle n-1}^{}$}&0\\ 0&0&\cdots&0&0&\makebox[0pt]{$\varepsilon_{\scriptscriptstyle n-1}^{}$ } & 0&h_n^{}\\ \varepsilon_n^{}&0&\cdots&0&0&0&h_n^{}&0 \end{array}\right)\ ; .\ ] ] since @xmath109 is symmetric and bipartite after an appropriate reordering of the basis . ]
, the eigenvalues occur in real pairs of opposite sign . but changing the sign of an eigenvalue simply means interchanging the corresponding creation and annihilation operators in eq .
( [ ffh ] ) , thus we can restrict ourselves to the positive part of the spectrum . in what follows , we concentrate on the case of a constant field @xmath110 , where we may set @xmath111 by an appropriate choice of the energy scale . of course
, one can not possibly expect analytical results for an arbitrary sequence of couplings .
however , for _ substitution _ sequences , more can be said using an exact renormalisation scheme , and , for an interesting class of chains with marginal fluctuations , the scaling exponent @xmath37 ( [ eq : scal ] ) can be determined exactly and given in a closed form .
consider now sequences of couplings chosen according to a substitution rule @xmath112 on an @xmath4-letter alphabet @xmath113 . here
, the @xmath114 are @xmath4 arbitrary words of finite length in the alphabet @xmath56 .
indeed , for the case of _ aperiodic two - letter _ substitution chains , this form imposes no restriction at all , since these can all be generated by substitution rules of the form ( [ subs ] ) . to see this , note that any two - letter substitution rule @xmath115 can be transformed appropriately by inner automorphisms , i.e. , by replacing @xmath116 [ where a natural choice of the word @xmath117 is the common beginning , if any , of the words @xmath118 and @xmath119 , and/or by considering @xmath120 instead .
obviously , this does not alter the limit chain .
note that this is not possible for certain _ periodic _ chains , namely those which are generated by a substitution rule that replaces all letters by the same word ( the period ) , or by a multiple of it .
it is no problem , however , to treat this case separately .
the eigenvalue problem for the matrix @xmath109 ( [ fh ] ) can be reformulated using transfer matrices in the usual way .
our renormalisation transformation ( rt ) is , however , most conveniently described in terms of the star - product formalism , borrowed from scattering theory @xcite .
we therefore introduce new @xmath3-transfer matrices , called @xmath3-matrices from now on , by @xmath121 with @xmath122 where the star - product ( ) of two @xmath123 matrices is defined as @xmath124 in our case , according to eq .
( [ fh ] ) , the elementary @xmath3-matrices are @xmath125 - \varepsilon_{j+1}^{-1 } h & \varepsilon_{j+1}^{-1 } \lambda   \end{array } \right ) \ ; , \ ] ] where @xmath126 denotes an eigenvalue may be _ any _ eigenvalue of @xmath109 . only later ,
after expanding the renormalisation equations in powers of @xmath126 [ eq .  ( [ hren ] ) and below ] , we shall restrict ourselves to eigenvalues that vanish in the limit of an infinite system size .
] of @xmath109 .    in order to establish the renormalisation scheme ,
we introduce a number of additional parameters for these @xmath3-matrices . for each letter @xmath127 , we define two asymmetry - coefficients @xmath128 and a field variable @xmath129 that enter the @xmath130 different elementary @xmath3-matrices through ) and ( [ fh ] ) in that the field variables are now labeled in the same way as the coupling constants ( i.e. , as the bonds ) rather then by the sites of the chain . ]
@xmath131 the rt now just
_ reverses _ the substitution steps . in our framework
, this is simply done by building -products of the @xmath3-transfer matrices that correspond to the words @xmath114 given by the substitution rule ( [ subs ] ) .
schematically , this works in the following fashion @xmath132 \cdots & \mbox{\large { \mbox{\boldmath $ a_i^{}$ } } } & \mbox{\large { \mbox{\boldmath $ a_j^{}$ } } } &            \mbox{\large { \mbox{\boldmath $ a_k^{}$ } } } & \cdots \\ & \multicolumn{3}{c } { \begin{picture}(50,30)(0,0 ) \thicklines \put (   0,20){\oval(100,35)[br ] } \put ( 50,20){\vector(0,1){5 } } \end{picture } \hphantom{\mbox{\large $ a$ } } \begin{picture}(116,30)(0,0 ) \thicklines \put (   0,20){\line(0,1){5 } } \put ( 58,20){\oval(116,35)[b ] } \put(116,20){\vector(0,1){5 } } \put ( 58,15){\makebox[0pt]{\large $ \tilde{s}_{a_i^{}|a_j^{}}$ } } \end{picture } \hphantom{\mbox{\large $ a$ } } \begin{picture}(116,30)(0,0 ) \thicklines \put (   0,20){\line(0,1){5 } } \put ( 58,20){\oval(116,35)[b ] } \put(116,20){\vector(0,1){5 } } \put ( 58,15){\makebox[0pt]{\large $ \tilde{s}_{a_j^{}|a_k^{}}$ } } \end{picture } \hphantom{\mbox{\large $ a$ } } \begin{picture}(50,30)(0,0 ) \thicklines \put (   0,20){\line(0,1){5 } } \put ( 50,20){\oval(100,35)[bl ] } \end{picture } } &   \end{array}\ ] ] where @xmath133 denotes the @xmath134th letter of the word @xmath114 . in this way
, we obtain renormalised @xmath3-matrices @xmath135 - \varepsilon_{a_j}^{-1 } \tilde{h}_{a_i }   & \varepsilon_{a_j}^{-1 } \tilde{\kappa}_{a_i}^-\tilde{\lambda }   \end{array } \right ) \ ; , \ ] ] compare eq .
( [ sj ] ) . as an _ additional condition _ on the rt , and in accordance with their initial values ,
we keep the product of the @xmath136 asymmetry - coefficients @xmath137 fixed : @xmath138 now , the rt for the different parameters can be read off directly from the -product relations .
explicitly , we obtain the following recursion relations @xmath139[r]{$\displaystyle h^{(1)}_{a_i}$ } = \makebox[10ex][l]{$\displaystyle h_{a_i}^ { } \ ; , $ } \qquad \makebox [ 5ex][r]{$\displaystyle h^{(k+1)}_{a_i}$ } = \makebox[25ex][l]{$\displaystyle                      \frac{-\varepsilon_{w_i^k } h_{w_i^k } h^{(k)}_{a_i } }                          { \varepsilon_{w_i^k}^2 -\lambda_{w_i^k+}^ { }                            \lambda_{a_i-}^{(k)}}\ ; , $ } \qquad \makebox [ 5ex][r]{$\displaystyle\tilde{h}_{a_i}^{}$ } = \makebox[15ex][l]{$\displaystyle h^{(|w_i|+1)}_{a_i}\ ; ; $ } &   \\[1 mm ] \label{rg2 } & \makebox [ 5ex][r]{$\displaystyle\lambda_{a_i-}^{(1)}$ } = \makebox[10ex][l]{$\displaystyle\lambda_{a_i-}^{}\ ; , $ } \qquad \makebox [ 5ex][r]{$\displaystyle\lambda_{a_i-}^{(k+1)}$ } = \makebox[25ex][l]{$\displaystyle\lambda_{w_i^k-}^ { } +                      \frac{h_{w_i^k}^2\lambda_{a_i-}^{(k ) } }                          { \varepsilon_{w_i^k}^2 - \lambda_{w_i^k+}^ { }                            \lambda_{a_i-}^{(k)}}\ ; , $ } \qquad \makebox [ 5ex][r]{$\displaystyle\tilde{\lambda}_{a_i-}^{}$ } =   \makebox[15ex][l]{$\displaystyle\lambda_{a_i-}^{(|w_i|+1)}\ ; ; $ } & \\[1 mm ]   \label{rg3 } & \makebox [ 5ex][r]{$\displaystyle\lambda_{a_i+}^{(1)}$ } = \makebox[10ex][l]{$\displaystyle\lambda_{a_i+}^{}\ ; , $ } \qquad \makebox [ 5ex][r]{$\displaystyle\lambda_{a_i+}^{(k+1)}$ } = \makebox[25ex][l]{$\displaystyle\lambda_{a_i+}^{(k ) } +                      \frac{\lambda_{w_i^k+}^ { }                           \left[h^{(k)}_{a_i}\right]^2 }                          { \varepsilon_{w_i^k}^2 - \lambda_{w_i^k+}^ { }                           \lambda_{a_i-}^{(k)}}\ ; , $ } \qquad   \makebox [ 5ex][r]{$\displaystyle\tilde{\lambda}_{a_i+}^{}$ } =   \makebox[15ex][l]{$\displaystyle\lambda_{a_i+}^{(|w_i|+1)}\ ; ; $ } & \end{aligned}\ ] ] for the field variables @xmath140 and for @xmath141 that are defined as @xmath142 using eq .
( [ kappacond ] ) , we obtain the renormalised fermion frequencies @xmath143
as @xmath144 note the following properties of our renormalisation approach :    * the rt is _ exact _ , since no new parameters are introduced in a renormalisation step . altogether , there are @xmath145 parameters for @xmath4 different couplings : effectively @xmath146 asymmetry - parameters @xmath147 , the @xmath4 fields @xmath129 , and the eigenvalue @xmath126 .
note that _ the coupling constants are not renormalised _ in our scheme . * the rescaling factor of the system size
is just given by the leading eigenvalue @xmath80 of the substitution matrix corresponding to the substitution rule under consideration .
* the critical point of the model , defined by the vanishing of the mass gap ( @xmath148 ) , corresponds to a fixed point of the rt for @xmath126 , as can easily be seen from eqs .
( [ rg2 ] ) and ( [ rg3 ] ) . *
the scaling behaviour in the vicinity of the critical point is governed by the leading - order term of the rt as @xmath149 .
we now proceed to establish the rt of the spectrum to leading order for @xmath150 explicitly .
this is most easily done by noting that , in the product ( [ sren ] ) , the common prefactors @xmath151 of the -products ( [ star ] ) gives rise to higher - order corrections only ; and therefore , for the present purpose , the denominator can be neglected completely .
consider first the renormalisation of the field variables .
we obtain @xmath152 instead of working directly with the rt for the fields , it is advantageous to consider new variables @xmath153 to linear order in @xmath126 , the rt of @xmath154 then yields the simple matrix form @xmath155 where @xmath156 is just the _ transpose _ of the substitution matrix of @xmath71 . now , take a look at the initial conditions @xmath157}$ ] for @xmath158 which are @xmath159 } = - 2\log \varepsilon_{a_i}$ ] . as stated above , the entries of the ( statistically normalised ) perron - frobenius eigenvector @xmath68 are the asymptotic letter frequencies @xmath160 .
we thus obtain @xmath161 = 0\ ] ] for an _ arbitrary _ set of couplings fulfilling the criticality condition ( [ eq : critcond ] ) .
hence , @xmath157}$ ] is nothing but an arbitrary ( real ) linear combination of the ( right ) eigenvectors of @xmath156 ( hence of the _ left _ eigenvectors of @xmath162 ) _ excluding _ the perron - frobenius one .
thus , for a generic choice of critical couplings , it scales with the second - largest eigenvalue , @xmath163}\| \sim |\lambda_2|^k$ ] , under the iteration @xmath164}={\mbox{\boldmath $ m$}}_{\varrho}^t { { \mbox{\boldmath $ y$}}}^{[k]}$ ] . translated into the terminology of the renormalisation group , we recognise @xmath158 as a _ scaling field _ transforming according to @xmath165 where the corresponding _ renormalisation group eigenvalue _ is just the wandering exponent @xmath89 defined in eq .  ( [ beta ] ) .
this demonstrates explicitly that power law , logarithmic , or bounded fluctuations lead to relevant , marginal , or irrelevant scaling fields , respectively .
we are going to describe below how this affects the critical behaviour .
we now proceed to the rt for @xmath166 .
defining @xmath167 \prod_{\ell = 1}^{|w_i| } \frac{h_{w_i^\ell}^2}{\varepsilon_{w_i^\ell}^2 } \ ; , \end{aligned}\ ] ] we find , by induction , as the linear part of the rt of @xmath166 @xmath168 setting @xmath169 , this can also be written in a matrix form @xmath170 since all the components of the matrices @xmath171 and of the vectors @xmath172 are positive , the vectors converge to the perron - frobenius eigenvectors of @xmath171 under iteration of the rt .
let @xmath173 denote the perron - frobenius eigenvalues of @xmath171 . from eq .
( [ kappacond2 ] ) , we conclude that @xmath174 for further convenience , we transform @xmath175 ( under conservation of the spectrum ) according to @xmath176 with a diagonal transformation matrix @xmath177 . as new entries ,
we find , dropping the index @xmath31 , @xmath178 this form shows explicitly that , to linear order in @xmath126 , only the _ reduced couplings _
@xmath179 enter in the rt . finally ,
since the rescaling factor of the chain length is the largest eigenvalue @xmath80 of the substitution matrix , we formally obtain the scaling exponent @xmath37 through @xmath180 let us now look at the different cases of fluctuation behaviour in more detail .      because @xmath181 is perpendicular to the perron - frobenius eigenvector of the substitution matrix , and since all other eigenvalues are smaller than one in modulus , we find @xmath163}\|\to 0 $ ] for all substitution chains with bounded fluctuations , hence @xmath182 . evaluating the transformation matrices @xmath171 at this fixed point of the fields , we obtain from eqs .
( [ pminus ] ) and ( [ pplus ] ) @xmath183 and hence @xmath184 thus we have @xmath185 for substitution chains with bounded fluctuations , and the low - energy spectrum at criticality scales as @xmath186 hereby , the @xmath126-transformation induces a finite renormalisation of the _ fermion velocity _ @xmath187 in comparison to the uniform chain @xmath188}/\lambda}\ ] ] where @xmath189}$ ] is the image of @xmath190}=\lambda$ ] after a @xmath134-fold application of the rt .    in general , all renormalisation steps contribute in eq .
( [ vren ] ) . however , if _ all _ eigenvalues of the substitution matrix but the perron - frobenius one are zero , the fixed point of the rt is already reached after the first renormalisation step and @xmath187 is obtained immediately .
this is shown below for periodic chains , but also allows to calculate the fermion velocity of the thue - morse and related chains .
using the equidistribution property of quasiperiodic _ cut - and - project _ chains , like our first example ( [ eq : sub ] ) with @xmath100 ( the fibonacci sequence ) , the limit in eq .
( [ vren ] ) can be carried through explicitly for these cases , too . since this has already been done in ref .
@xcite , we do not repeat the arguments here , and merely cite the result @xmath191 for two couplings parametrised as in eq .
( [ eq : para ] ) .      in this short part
, we briefly show how the known results for the fermion velocity @xmath187 for _ periodic _ quantum chains can be recovered easily within our renormalisation scheme . by a substitution rule of the form @xmath192 _
any _ periodic chain with period @xmath193 can be generated . in this simple situation
, we can take all the fields @xmath194 and asymmetry - parameters @xmath195 to be equal .
by virtue of the criticality condition ( [ eq : critcond ] ) , the reduced couplings transform as @xmath196 [ see eq .
( [ hren ] ) ] , and we obtain the ( fermionic ) hamiltonian of a _ uniform _ chain ( @xmath197 ) after a single renormalisation step . note that the additional asymmetry - parameters result in an irrelevant similarity transformation of the matrix ( [ fh ] ) only .
the renormalisation factor for @xmath126 can be deduced from eqs .
( [ pminus ] ) and ( [ pplus ] ) , explicitly it reads ( using eq .
( [ eq : critcond ] ) , and normalising to @xmath198 ) : @xmath199 \label{v } & = & \sum_{k=1}^{|w|}\sum_{l=1}^{|w| } \varepsilon_k^2 \cdots   \varepsilon_{k+l-1({\rm mod}\,|w|)}^2 \end{aligned}\ ] ] with @xmath200 . the critical scaling of the fermion frequencies now results in @xmath201 where @xmath187 is given in eq .
( [ v ] ) , and where the well - known scaling behaviour of the uniform chain is recovered for @xmath202 .
this result was first derived in the context of the 2@xmath2 classical ising model with layered periodicity in ref .
@xcite .    in general , there is no closed expression for the fermion velocity @xmath187 .
for the special case of a uniform chain of @xmath203-couplings with periodic @xmath204-defects on every @xmath81th site , we derive @xmath205}{\sinh [ \log ( r)]/m } \ ; .\ ] ] note that @xmath206 for @xmath207 ( simple periodic chain ) , but also @xmath208 for @xmath209 since the parametrisation ( [ eq : para ] ) implies that @xmath210 in this limit .
we now consider the case where the second - largest eigenvalue(s ) of the substitution matrix ( and of its transpose ) have @xmath211 .
let @xmath212 denote the corresponding @xmath81-dimensional ( joined ) eigenspace of @xmath213 , where @xmath214 . by eqs .
( [ y1 ] ) and ( [ y2 ] ) , @xmath215}$ ] converges to the projection of @xmath216}$ ] on @xmath212 under iteration of the rt .
the critical surface being @xmath217-dimensional , we can parametrise @xmath218}$ ] in the renormalisation limit as @xmath219 with @xmath81 vectors @xmath220 that span @xmath212 and coefficients @xmath221 , @xmath222 .
this implies @xmath223 and , if the vector components @xmath224 are chosen as integers , the entries of the matrices @xmath171 are polynomials with positive coefficients of @xmath81 parameters @xmath225 .
if all @xmath226 are equal to 1 , we obviously obtain @xmath227 furthermore , if @xmath215}$ ] converges to an eigenvector of @xmath213 to the eigenvalue @xmath228 , @xmath175 and @xmath229 are related by @xmath230 to obtain the scaling exponent @xmath37 from eq .
( [ zz ] ) , we need to know the perron - frobenius eigenvalues @xmath173 of @xmath171 . there is no explicit solution to that problem in the general case . in what follows ,
we present the complete answers for a number of more special , yet infinite , classes of substitution rules .
let @xmath231 be any word from the alphabet , and @xmath232 its @xmath134th cyclic permutation .
consider now a substitution rule @xmath71 of the form ( [ subs ] ) with @xmath233^{l_i } \ ; .\ ] ] one special case is to take equal words @xmath234 for all @xmath235 .
the corresponding substitution matrices indeed have @xmath236 as their second - largest eigenvalue , with degeneracy @xmath237 , since the column vectors of @xmath238 differ by scalar factors only .
thus the fields stay at their original values and are not renormalised at all .
to derive @xmath173 as a function of the parameters @xmath239 , we simply note that the special form of @xmath71 transforms this property to @xmath171 , not just for the point @xmath240 , but _
independently _ of the parameters .
we thus obtain @xmath241 and the scaling exponent is given by @xmath242 where the @xmath243 are given in eqs .
( [ pminus ] ) and ( [ pplus ] ) as functions of the couplings .
since @xmath173 are polynomials of the parameters @xmath244 with only positive coefficients , and @xmath245 , we conclude @xmath246 . as a function of the parameters , @xmath37 is either purely _ convex _ , with unique minimum @xmath247 ( the limit of the uniform chain ) , or it is _ constant _ :
it is interesting to take a closer look at the latter case .
obviously , this means that @xmath249_{ii}$ ] , independently of the parameters .
as can be shown , for substitution rules of the form ( [ subs ] ) with ( [ rhoform2 ] ) , this is possible if and only if @xmath231 is a word ( or a power of a word ) with all letters @xmath250 appearing precisely once and if the last letter of @xmath114 is @xmath250 .
this means , however , that the resulting chain is _ periodic _ and has bounded rather than logarithmically diverging fluctuations .
as stated above , each two - letter substitution chain can be generated by a substitution rule of the form ( [ subs ] ) .
for any marginal substitution rule , we can assume that the second - largest eigenvalue of the substitution matrix is @xmath251 .
if necessary , this can be achieved by going to @xmath252 . as a consequence , @xmath216}$ ] in eq .
( [ y1 ] ) is eigenvector to the eigenvalue @xmath228 of @xmath213 , and as above the fields are already at their fixed points .    in this case , the entries of the transformation matrices @xmath171 are polynomials of a single parameter @xmath253 that determines the critical couplings .
the crucial point now is to show that  like in the case discussed above  the matrices @xmath171 have an eigenvalue @xmath254 independently of @xmath253 and of the detailed form of the substitution rule .
this assertion is proved in appendix  a. then , the perron - frobenius eigenvalues of @xmath171 are given as above ( [ mu ] ) , and we obtain the critical scaling exponent @xmath255 where @xmath80 is the leading eigenvalue of the substitution matrix . as in the case discussed above , @xmath256 is constant ( @xmath257 ) only for the special case of a periodic chain ( with period @xmath258 ) .
otherwise , @xmath37 is a convex function of @xmath253 , with a unique minimum @xmath259 where the ( reduced ) couplings are equal . on the other hand ,
it is easy to see that the periodic chain with period @xmath258 is the only two - letter substitution chain with bounded letter fluctuations that can be generated by a substitution rule of the form ( [ subs ] ) and that has @xmath254 as the second - largest eigenvalue of the substitution matrix .
we thus conclude that @xmath257 _ if and only if _ the fluctuations are bounded , and @xmath260 otherwise .
the _ period - doubling chain _ , given by the substitution rule ( [ eq : sub ] ) with @xmath101 , falls into both classes described above .
we consider double substitution steps , eliminating blocks of the form @xmath261 and obtain the scaling exponent @xmath262 this result was reported before in ref .
if the length of the chain is increased by simple substitution steps , the coefficients @xmath263 in eq .
( [ zz ] ) converge _ separately _ for an even and an odd number of steps .
effectively , this can be seen as a renormalisation of the couplings according to @xmath264 in a single renormalisation step .
hence we obtain @xmath265 .    as an example that is not contained in the above classes ,
consider the three - letter substitution @xmath266 with eigenvalues of the substitution matrix @xmath267 and @xmath236 .
the eigenvector to eigenvalue 1 of the transposed substitution matrix @xmath213 is @xmath268 , and with eq .
( [ y0 ] ) we conclude that @xmath269 and @xmath270 . evaluating the rt of @xmath172 at this fixed point
, we obtain the critical scaling exponent @xmath271 note that , for @xmath272 , the fluctuations of the couplings are indeed bounded , and we obtain @xmath273 as it should be . in our last example , the couplings follow the thue - morse sequence on the _ sites _ of the chain , rather than on the bonds as in all cases discussed so far . the coupling constants for a site - problem
are given as a function of the two adjacent bonds  and thus can take in general @xmath130 different values ( @xmath274 in this case ) .
it can , however , easily be shown that this can be reformulated as a bond - problem with a substitution chain on an @xmath130-letter alphabet @xcite . in this case
, the resulting four - letter substitution rule is @xmath275 the asymptotic letter frequencies are @xmath276 and @xmath277 ; and so we have @xmath278 ( for @xmath279 ) as the criticality condition ( [ eq : critcond ] ) .
a possible parametrisation is @xmath280 , @xmath281 , @xmath282 , and @xmath283 .
the eigenvalues of the substitution matrix are @xmath284 , @xmath285 , @xmath285 , and @xmath286 , we thus expect marginal fluctuations , whereas the thue - morse bond - problem has bounded fluctuations . following the above steps ,
we obtain a scaling exponent @xmath287}{\log ( 2 ) } \ ; .\ ] ] note that for @xmath288 , @xmath272 we obtain another coding of the period - doubling chain , while the special case @xmath289 , @xmath290 leads us back to the thue - morse bond - problem with bounded fluctuations ( with @xmath291 , and thus @xmath273 ) .      for @xmath87 ,
fluctuations diverge with a power law . as a consequence , the vector @xmath215}$ ] finally scales with @xmath292 under the rt .
since the matrix elements of @xmath171 are proportional to @xmath293 , at least a part of them eventually diverge like @xmath294 [ note that the components of @xmath181 have different sign because of eq .
( [ y2 ] ) ] . estimating the product of the perron - frobenius eigenvalues @xmath173 by the product of the corresponding traces , we obtain @xmath295 and thus @xmath296 where the _ wandering exponent _ @xmath89 is defined in eq .
( [ beta ] ) .
let us discuss this in more detail for the case of two - letter substitutions .
we take @xmath297 with corresponding eigenvector @xmath298 of @xmath213 and parametrise the critical couplings as @xmath299 since @xmath300 and @xmath301 have different sign , either one of these two is bigger than one for a fixed value @xmath302 , and without loss of generality we assume @xmath303 .
note that we have the following relation @xmath304 we now show that @xmath305 $ ] . from eq .
( [ pminus ] ) , this is immediately clear for @xmath306 , but it also holds for @xmath307 . to see this ,
consider the products of subsequent couplings of the word @xmath65 .
the product of the entire word is given by eq .
( [ product ] ) .
if the last letter of @xmath65 is an @xmath63 , this enters @xmath307 and we are done . on the other hand ,
if the last letter is a @xmath64 , consider the product of all couplings but the last @xmath308 , which is even larger .
iterating this argument , we finally arrive at a letter @xmath63 which proves our estimate of @xmath307 . using the standard eigenvalue formula
, we obtain the scaling of @xmath173 under the rt , and thereby @xmath309 here , @xmath310 is some constant that may depends on the sign of @xmath311 , but apart from that is independent of @xmath253 and @xmath15 .
this scaling form was also predicted in @xcite .
in this section , we give a short summary of the results of our renormalisation treatment and discuss their consequences for the critical properties of the iqc .      here , the behaviour of the low - energy part of the the spectrum at criticality is essentially the same as for the uniform chain .
the correlation length critical exponents stays at its onsager value @xmath312 . for periodic chains ( or periodic variations of the interactions along the layers of the 2@xmath2 statistical system ) , the results from renormalisation show that the smallest fermion frequencies are only altered by a common factor @xmath187 ( fermion velocity ) , thus also conformal properties persist and the whole model is described by the central charge @xmath313 conformal field theory of a massless free majorana fermion : @xmath314 with a fermion velocity @xmath187 given by ( [ v]).the scaling dimensions @xmath315 are odd integers for periodic boundary conditions ( mixed sector hamiltonian ) and odd half - integers for free boundaries @xcite .
numerical results @xcite confirm this conformal behaviour even for non - periodic chains with bounded fluctuations .
for the critical scaling , we can make the following _ ansatz _ : @xmath316 the index @xmath253 here denotes the dependence on the position on the critical surface .
the predicted scaling behaviour @xcite is confirmed analytically for all marginal two - letter substitution chains and some infinite classes of @xmath4-letter substitution chains .
@xmath317 converges to the scaling coefficients @xmath263 if we increase the length of the chain @xmath15 by ( multiple ) substitution steps , in other words @xmath318 is asymptotically periodic for @xmath319 with period @xmath320 , where @xmath80 is the scaling factor of the system size .
( [ z ] ) implies @xmath321 , and @xmath273 is only obtained in the limit of the uniform chain where the fluctuations actually vanish .
we further note that , in the explicitly solved cases , @xmath37 is invariant under inversion of all couplings . thus this is a symmetry of the lower part of the critical spectrum , while it is easy to see that it is not a symmetry of the entire finite - size spectrum . as a consequence
, @xmath37 quadratically approaches the isotropic value @xmath322 for small deviations of the couplings from the uniform strength @xmath323 , as was also predicted by luck  @xcite .    in the terminology of the 2@xmath2 classical ising model
, we have @xmath324 for any marginal disorder within layers of couplings .
since the correlation length critical exponent perpendicular to the layers remains @xmath325 , marginal fluctuations introduce a relevant anisotropy of the model , as was noticed in ref .
@xcite . due to anisotropic hyperscaling , a _ negative _ specific heat exponent @xmath326
is expected .    within the renormalisation approach
, the length of the chain @xmath15 is increased by ( multiple ) substitution steps .
thus , it is a priori not clear what happens to the scaling dimensions when we increase @xmath15 arbitrarily . indeed ,
numerical studies ( for free boundary conditions ) show that , if we increase @xmath15 by arbitrary amounts , while @xmath37 converges , @xmath327 stays finite but will not tend to a unique limit .
this is in contrast to the bounded fluctuation case . in more detail , the following properties of @xmath327
have been observed in all examples ( of two - letter substitution chains ) studied :    * given the couplings , @xmath327 appears to be bounded for all @xmath15 and @xmath10 . * for @xmath328 ,
the effect of @xmath10 on @xmath329 is a mere shift : @xmath330 .
* as a consequence of the rt , @xmath331 , @xmath318 is asymptotically periodic . *
the amplitude of @xmath329 vanishes in the limit of the uniform chain , @xmath332 , as @xmath329 approaches continuously the constant value @xmath333 .
these points can immediately be translated into properties of the integrated density of states ( idos ) , the inverse of which is given by @xmath334 while the first two points imply the convergence of this quantity , the third point above implies the lower part of the idos to be self - similar with a discrete scaling symmetry @xmath335 where @xmath336 is the asymptotic period of @xmath329 .
mind that in general not the spectrum ( the lowest fermion frequencies ) , but only the idos displays this scaling behaviour for @xmath337 .
this scaling property has consequences for the possible gap structure of the critical spectrum , where the results of the renormalisation approach at least lead to necessary conditions for gap - labelling theorems .
but since within this class more spectral properties can be found analytically using trace maps , we postpone the discussion to a forthcoming publication .      in substitution sequences with strong fluctuations , as in our third example , the binary non - pisot sequence [ @xmath102 in eq .
( [ eq : sub ] ) ] , their divergence is described @xcite by the wandering exponent @xmath89 ( [ beta ] ) which determines the exponential closing of the spectral gaps as @xmath338 here , @xmath339 is the second moment of the distribution of coupling constants along the chain @xcite which for the binary non - pisot sequence is given by @xmath340 and compatibility with the uniform case requires @xmath341 and @xmath342 .    from the results of the renormalisation approach ,
@xmath318 should be asymptotically periodic for @xmath331 .
this is confirmed by numerical investigations .
indeed , the latter suggest the very same properties as given above in the case of marginal fluctuations .
in particular , we observe @xmath343 for @xmath344 .
again , this means that the integrated spectral density converges in the thermodynamic limit and its lower part shows the following scaling property at the critical point : @xmath345
using a renormalisation approach , we obtained analytic results for the critical behaviour of a class of ising models intermediate between uniformly ( or periodically ) ordered and randomly disordered systems .
the systems under consideration are ferromagnetic ising quantum chains whose interaction constants are modulated according to substitution sequences .
luck s criterion @xcite for the relevance of aperiodic disorder on the critical behaviour is fully confirmed for these systems , bringing along increasing evidence that the _ fluctuation of the interactions _ is the basic concept that demarcates the onsager universality class from models with weaker critical singularities . for two - letter substitution rules ,
this analysis was carried out in detail , showing how the renormalisation flow is determined by the nature of the fluctuations .
an exact renormalisation formalism for the most general case of @xmath4-letter substitution rules was also presented and applied to examples .
a number of quantitative results have been displayed , mainly for the correlation length critical exponent ( resp .
the scaling exponent @xmath37 of the mass gap ) , confirming and refining luck s predictions for bounded , marginal , and relevant fluctuations of the coupling constants .
in particular , an exact formula has been obtained for the scaling exponent of arbitrary two - letter substitution chains , containing previous results as special cases . in this light ,
the case of two - letter substitution rules appears to be rather well understood by now .
what remains to be done in the @xmath4-letter case is a thorough discrimination between quantum ( substitution ) chains with bounded and with marginal fluctuations of the couplings .
this is connected with the problem of determining all substitution rules with second - largest eigenvalue of the substitution matrix on the unit circle , which nevertheless result in a chain with bounded rather than marginal fluctuations .
also , the investigation of ordered cases that do not stem from a substitution rule might add some insight as they can possess a higher degree of complexity .
unfortunately , renormalisation techniques will probably be unsuitable here .    of primary interest
are , of course , extensions to higher dimensions .
a natural first step is to study two - dimensional statistical systems with two - dimensional rather than one - dimensional variations of the interactions . as can be seen from some special cases @xcite , analytic results are possible in principle , but they are not generic and solvable cases might actually be misleading .
further analysis is necessary , and some results can be expected from numerical approaches , especially since structures generated by substitution rules are well - suited for ( numerical ) renormalisation .
first results in this direction have been obtained recently @xcite and will be extended soon .
part of this work was done while jh was visiting melbourne university .
it is a great pleasure to thank paul pearce for his kind hospitality . financial support from the german science foundation ( dfg ) and the studienstiftung des deutschen volkes
is gratefully acknowledged .
we consider two - letter chains generated by substitution rules @xmath346 which have @xmath347 as the second - largest eigenvalue of the corresponding substitution matrix . as a consequence ,
one finds the constraints @xmath348 on the words @xmath65 and @xmath66 and the letter frequencies @xmath349 in the limit word , where @xmath350 denotes the number of letters @xmath94 contained in the word @xmath351 . we parametrise the criticality condition ( [ eq : critcond ] ) @xmath352 ( setting @xmath353 ) for the corresponding ising chain by @xmath354 with _ coprime _
integers @xmath134 and @xmath355 , such that @xmath356 .
then @xmath357 divides @xmath358 and @xmath359 .
the scaling exponent of the critical spectrum is a simple expression of the leading eigenvalues of two matrices @xmath360 where the entries @xmath361 are defined in eqs .
( [ pminus ] ) and ( [ pplus ] ) .
they are polynomials in @xmath253 and @xmath362 . in what follows
, we prove that @xmath171 have an eigenvalue @xmath254 , _ independently _ of the value of @xmath253 and of the detailed form of the substitution rule @xmath71 .
since we have @xmath363 , it suffices to concentrate on @xmath175 .
we define two polynomials @xmath364 corresponding to the words @xmath65 and @xmath66 by @xmath365 taking the frequencies @xmath349 of the letters as given , the possible words @xmath65 and @xmath66 are subject to the constraints ( [ eq : wawb ] ) .
we define an _ elementary polynomial _ @xmath366 to be the polynomial @xmath364 corresponding to a word @xmath367 of minimal length fulfilling this condition .
since @xmath134 and @xmath355 are coprime , each @xmath367 consists of @xmath134 @xmath63 s and @xmath355 @xmath64 s and leads to one of @xmath368 different elementary polynomials possible .
let in the following @xmath231 be @xmath65 or @xmath66 .
we now introduce a graphical representation of the word @xmath231 and the accompanying polynomials .
for given @xmath231 , construct iteratively a step function @xmath369 \to \mathbb{z}$ ] as follows : @xmath370 note that the criticality conditions imply @xmath371 .
an example for the case @xmath101 , @xmath372 is shown in figure  [ fig1 ] .    for @xmath373 ( @xmath374 ) , @xmath375 ( resp .
@xmath376 ) is given by the `` sum over the _
upward _ steps '' of @xmath377 ( proceeding in positive @xmath378-direction ) , while @xmath379 ( resp .
@xmath380 ) is given by the `` sum over the _
steps '' , each step contributing a term @xmath381 . in our example
( for @xmath373 ) , @xmath382 and @xmath383 .    in a first step , we now show that @xmath364 can be expressed as sums over elementary polynomials @xmath384 .
_ lemma:_for an arbitrary word @xmath231 fulfilling the criticality condition for given @xmath349 , there are integers @xmath250 and elementary polynomials @xmath385 , such that @xmath386    _ proof:_consider the substrings @xmath387 of @xmath231 of length @xmath357 starting with the letter @xmath388 .
for each @xmath387 containing more than @xmath134 @xmath63 s , there has to be an @xmath389 with less than @xmath134 @xmath63 s , because of the criticality condition for @xmath231 . since the number of @xmath63 s contained in successive substrings @xmath387 and @xmath390 differs at most by one
, there is at least one substring consisting of @xmath134 @xmath63 s and @xmath355 @xmath64 s . in the graphical representation
, this appears as a restriction of @xmath377 to an interval of length @xmath357 with the endpoints of the graph taking the same value .
the corresponding term within @xmath364 is obviously some elementary polynomial times a power of @xmath253 .
we may now eliminate this substring from @xmath231 ( cut out the interval in the graphical picture ) and start the argument again .
we proceed by showing that the elementary polynomials differ only in a ( polynomial ) factor that is independent of @xmath391 .    _
proposition:_given @xmath134 and @xmath355 , there are polynomials @xmath392 and @xmath393 such that for every elementary polynomial @xmath394 @xmath395\ ] ] where @xmath396 is a polynomial that does not depend on @xmath391 .    _
proof:_an elementary polynomial @xmath394 takes the following general form @xmath397 where @xmath398 are integers . defining @xmath399
the proof then follows by direct calculation .
we remark that with this choice of @xmath400 and @xmath401 , @xmath402 can also be shown to be a polynomial in @xmath253 and @xmath362 .
now , we are in the position to complete our argument . from eqs .
( [ eq : pab ] ) , ( [ lem ] ) and ( [ prop ] ) , we conclude that @xmath403 thus it follows that @xmath404 , and hence @xmath236 is an eigenvalue of @xmath175 for an _ arbitrary _ value of @xmath253 .
the leading eigenvalues of @xmath405 are given by eq .
( [ mu ] ) .
in the above discussion , the coupling constants were chosen according to a substitution rule of the special form ( [ subs ] ) . while any ( infinite ) substitution chain with two letters can be generated that way ,
this is no longer the case for three or more different letters ( resp .
couplings ) . in this appendix , we present a generalised version of the renormalisation formalism to deal with the general case . consider a substitution rule of the form @xmath407 the main problem in the general case is to maintain the recursive structure of the renormalisation procedure within the elimination process . by fixing the first letter of the substitutes in eq .
( [ subs ] ) , the appropriate couplings of the renormalised chain were given .
this is no longer possible here ; however , we can proceed as follows .
* fix any letter @xmath250 within the word @xmath408 , thus defining two words @xmath409 and @xmath410 as shown in eq .
( [ subs2 ] ) .
note that the letter @xmath250 appears in the word @xmath114 ( at least for a suitably chosen power @xmath411 of the substitution rule ) for primitive substitutions . *
redefine the substitution rule as a function on _ pairs _ of letters in the following way @xmath412 obviously , ( [ subs3 ] ) leads to the same limit chain as ( [ subs2 ] ) .
the substitution matrix of @xmath115 has a dimension of at most @xmath130 ( pairs @xmath413 that do not appear in the original chain need not be included in @xmath115 ) .
actually , @xmath414 is just the substitution matrix of the corresponding _
site_-problem , where the type of site follows the sequence , and the bond ( resp .
the interaction along it ) is a function of its two endpoints @xcite .
note , however , that the couplings here are still attached to the bonds . * in the renormalisation formalism , now define @xmath3-matrices , fields and asymmetry - parameters corresponding to pairs of letters @xmath415 renormalisation again means inverting the substitution process by integrating out all degrees of freedom attached to the words @xmath409 and @xmath410 .
we obtain renormalised @xmath3-matrices through @xmath416 * the rt ( and their linear orders in @xmath126 ) for the ( at most ) @xmath417 parameters are then obtained analogous to section  3 .    for a classification of the critical behaviour based on the fluctuations we need to know about the spectrum @xmath418 of @xmath414 in dependence of the spectrum @xmath419 of @xmath72 . for an appropriate power of @xmath71 , it is indeed possible to show that @xmath418 contains @xmath419 and that all additional eigenvalues are either 1 or 0 .
we only sketch the proof here . without restriction of the general case
, we can assume that the first and the last letters of all words @xmath114 remain fixed under the substitution ; that is @xmath420 ^ 1 $ ] for the first letter , and equivalently for the last letter .
this is always fulfilled by a finite power @xmath421 of the substitution rule .
we assume the @xmath4 words @xmath114 to have @xmath134 different first letters and @xmath10 different last letters . taking a look at the associated pair substitution matrix @xmath414 ,
it is easy to show that @xmath422 since we always have @xmath423 and the numbers @xmath134 and @xmath10 remain fixed because of the above assumption , the traces of @xmath424 and @xmath425 differ only by the constant @xmath426 for arbitrary @xmath81 .
we conclude that @xmath418 containes @xmath419 and the all additional eigenvalues are fixed points under any power , thus @xmath285 or @xmath286 .
we will , however , not give a detailed discussion of the general case here , but rather illustrate the method by applying it to a special example .
consider the substitution rule of the `` circle sequence '' .
this example of a quasiperiodic chain ( in the sense that it has a pure point fourier spectrum ) has been studied numerically in ref .
@xcite @xmath427 with @xmath252 being of the form ( [ subs3 ] ) .
the eigenvalues of @xmath428 are @xmath429 , @xmath285 , and @xmath430 , where @xmath431 is the golden mean
we thus expect marginal fluctuations . the asymptotic frequencies of the letters are @xmath432 , @xmath433 , and @xmath434 , respectively . we choose @xmath435 , @xmath436 , and @xmath437 as the parametrisation of the critical couplings @xmath438 .
since only the five pairs @xmath258 , @xmath439 , @xmath440 , @xmath441 , and @xmath442 appear in the chain , the dimension of the pair transfer matrix @xmath443 is just five .
the spectrum @xmath444 consists of the spectrum @xmath445 and of two additional eigenvalues 1 . identifying pairs of couplings with the first one , we clearly have @xmath446 , @xmath447 , and @xmath448 .
thus the vector of the logarithms of the reduced couplings , @xmath449 , is sill perpendicular to the perron - frobenius eigenvector of @xmath414 and converges in the renormalisation limit to a linear combination of two eigenvectors corresponding to an eigenvalue 1 of @xmath450 : @xmath451 .
while both parameters , @xmath253 and @xmath452 , enter the transformation matrices @xmath171 of @xmath172 , @xmath452 can actually be eliminated by a similarity transformation .
such a behaviour is always to be expected if a parameter is connected to an eigenvector ( like in this case ) that corresponds to an _ additional _ eigenvalue 1 of @xmath453 , not contained in @xmath445 .
the perron - frobenius eigenvalues of @xmath171 can be determined explicitly , and we obtain for the scaling exponent @xmath454 finally , we like to show an example of a three - letter substitution rule @xmath455 with bounded fluctuations , which nevertheless has @xmath236 as the second - largest eigenvalue of the substitution matrix .
this is the bond substitution that corresponds to the site - problem given by @xmath456 .
fluctuations are indeed bounded here , since we just obtain the fibonacci chain by identifying @xmath63 and @xmath64 , but also by identifying @xmath63 and @xmath310 . calculating the scaling exponent by the above method , we obtain @xmath273 , as was to be expected .    99 h.  au - yang and b.m .
mccoy , theory of layered ising models : thermodynamics , _ phys .
rev .  * b10 * _ ( 1974 ) 88691
. m.  baake , p.  chaselon and m.  schlottmann , the ising quantum chain with defects ( ii ) .
the @xmath457 kac moody spectra , _ nucl .
phys . _ * b314 * ( 1989 ) 62545 .
m.  baake , u.  grimm and r.  baxter , a critical ising model on the labyrinth , _ int .  j.  mod .
* b8 * ( 1994 ) 3579600 .
m.  baake , u.  grimm and d.  joseph , trace maps , invariants and some of their applications , _ int .
j.  mod .
_ * b7 * ( 1993 ) 152750 .
berche and b.  berche , crossover between aperiodic and homogeneous semi - infinite critical behaviors in multilayered two - dimensional ising models , _ phys .
* b * ( to appear ) .
b.  berche , p .- e .
berche , m.  henkel , f.  igli , p.  lajk , s.  morgan and l.  turban , anisotropic scaling in layered aperiodic ising systems , _ j.  phys . _ * a28 * ( 1995 ) l16571 .
b.  berche , p .- e .
berche and l.  turban , marginal anisotropy in layered aperiodic ising systems , _
j.  phys .
i france _ * 6 * ( 1996 ) 62140 .
cardy , effect of boundary conditions on the operator content of two - dimensional conformally invariant theories , _ nucl .
phys . _ * b275 * ( 1986 ) 20018 ; and : conformal invariance , in : _ phase transitions and critical phenomena _ ,
eds .  c.  domb and j.l .
lebowitz , academic press , london ( 1987 ) , pp .
dumont , summation formulae for substitutions on a finite alphabet , in : _ number theory and physics _ ,
luck , p.  moussa , and m.  waldschmidt , springer , berlin ( 1990 ) , pp .
e.  fradkin , l.  susskind , order and disorder in gauge systems and magnets , _ phys .
* d17 * ( 1978 ) 263758 .
u.  grimm , the quantum ising chain with a generalized defect , _ nucl .
_ * b340 * ( 1990 ) 63358 .
u.  grimm and m.  baake , non - periodic ising quantum chains and conformal invariance , _ j.  stat .  phys . _
* 74 * ( 1994 ) 123345 .
u.  grimm and m.  baake , aperiodic ising models , in : _ the mathematics of long - range aperiodic order _ , ed .
.  moody , kluwer , dordrecht ( 1997 ) , pp .
harris , effect of random defects on the critical behaviour of ising models , _
j.  phys .
_ * c7 * ( 1974 ) 167192 .
f.  igli and l.  turban , common trends in the critical behavior of the ising and directed walk models , _ phys .
* 77 * ( 1996 ) 12069 .
lieb , t.  schulz and d.  mattis , two soluble models of an antiferromagnetic chain , _ ann .  phys .
( ny ) _ * 16 * ( 1961 ) 40766 .
j.  kogut , an introduction to lattice gauge theory and spin systems , _ rev .
* 51 * ( 1997 ) 656713 .
luck , critical behavior of the aperiodic quantum ising chain in a transverse magnetic field , _
j.  stat .
* 72 * ( 1993 ) 41758 ; and : a classification of critical phenomena on quasi - crystals and other aperiodic structures , _ europhys .
* 24 * ( 1993 ) 35964 .
f.  lunnon and p.  a.  b.  pleasants , quasicrystallographic tilings , _ j.  math .
pures et appl . _
* 66 * ( 1987 ) 21763 .
pfeuty , an exact result for the 1d random ising model in a transverse field , _ phys .
lett . _ * a72 * ( 1979 ) 24546 m.  queffelec , _ substitution dynamical systems _ , lecture notes in mathematics , vol .
1294 , springer , berlin ( 1987 ) .
r.  redheffer , difference equations and functional equations in transmission - line theory , in : _ modern mathematics for the engineer _
e.  f.  beckenbach , mcgraw hill , new york ( 1961 ) , pp .
. h.  simon , _ ferromagnetische isingmodelle auf aperiodischen strukturen _
, phd thesis , universitt tbingen , dissertations druck darmstadt ( 1997 ) .
tracy , universality classes of some aperiodic ising models , _
j.  phys .
_ * a21 * ( 1988 ) l6035 .
l.  turban , p .- e .
berche and b.  berche , surface magnetisation of aperiodic ising systems : a comparative study of bond and site problems , _ j.  phys . _
* a27 * ( 1994 ) 634966 .
l.  turban , f.  igli and b.  berche , surface magnetization and critical behavior of aperiodic ising quantum chains , _ phys .  rev . _
* 49 * ( 1994 ) 12695702 .