vibrational excitations in disordered structures is an active field of investigation by different methods .
the first approaches were developed many decades ago for the case of substitutional alloys ( see e.g. ref .
@xcite and references therein ) .
disorder was introduced there by the random positioning of two or more atomic species ( of different masses ) onto lattice sites .
successful analytical methods have also been adapted for such problems @xcite .
the rigidity - percolation problem has been solved to some extent numerically and within the mean - field approximations @xcite .
atomic vibrations in structures with fractal geometries were investigated in refs .
@xcite . some phenomenological @xcite and more realistic @xcite models of atomic vibrations in glasses
have been developed in order to describe the thermodynamical and phonon - transport peculiarities of such materials .
various computer experiments have also been performed on disordered structures ( see e.g. @xcite and references therein ) .
recently , the methods of the mode - coupling theory have been applied to investigate short - time vibrational dynamics , especially in the low - energy regime @xcite .
the low - energy behaviour of the vibrational density of states and dynamical structure factor have been investigated in disordered lattice models within the coherent potential approximation ( cpa ) @xcite , and in positionally disordered models treated perturbatively in the low inverse - particle - density limit @xcite .
the spectral properties of the euclidian random matrices associated with topologically disordered systems @xcite have been investigated within a supersymmetric statistic field theory in ref .
the spectral properties of the laplacian defined on a random graph @xcite and  small - world
lattices @xcite are related to topologically disordered solids and have also been investigated .
previously , vibrational displacements in disordered structures have been treated as scalars , which is a convenient simplification of the vector nature of atomic vibrations ( although vector vibrations were considered in the rigidity - percolation problem @xcite ) . however , a consistent treatment of vector vibrations in disordered lattices has not been performed up to now . the main purpose of this paper is to fill this gap and to give a mean - field description of vector vibrational models , namely : ( i ) to derive the cpa equations for vector vibrations in disordered lattices with force - constant disorder ; ( ii ) to solve them analytically in limiting cases and numerically through the whole energy range ; ( iii ) to apply the solution to investigate the scattering properties of external plane - wave vibrational excitations and , in particular , the ioffe - regel crossover from a weak- to strong - scattering regime ; ( iv ) to demonstrate the importance of sum - rule correlations in the dynamical matrix for scattering properties of vibrational external plane waves .
one of the main results of the paper is the phase diagram for vector vibrations in the representative f.c.c .
lattice showing the ranges of weak- and strong - scattering regimes for different degrees of disorder .
another result concerns a comparative analysis of the vibrational problem with the electronic one , and an explanation of the essential differences in phase diagrams for these problems .
all the major derivations and conclusions are supported by precise numerical solutions for the same model , which demonstrate the very good performance of the mean - field cpa approach , both for the density of states and the spectral density in a plane - wave basis ( dynamical structure factor ) in the whole energy range , except for the high - frequency tail region containing localized states .
the rest of the paper is arranged in the following manner .
the model is defined in sec .
the properties of the effective mean field are analysed in sec .
the spectral densities and the scattering properties of the external plane waves are discussed in sec .
the conclusions are given in sec .
atomic vibrations in the harmonic approximation can be treated in the classical limit @xcite . thus the problem is reduced to the hamiltonian formalism , in which the dynamical operator plays the role of the hamiltonian , and the squared frequency , @xmath0 , is an energy , @xmath1 .
the relevant hamiltonian ( dynamical operator ) for vector vibrations of @xmath2 interacting atoms of unit masses ( @xmath3 ) can be written in the following form : @xmath4 where we have used the site orthonormal basis of vectors @xmath5 , the only non - zero component of which is related to the displacement of atom @xmath6 ( @xmath7 ) along the cartesian direction @xmath8 ( @xmath9 ) .
the ( dynamical ) matrix elements @xmath10 , being the second partial derivatives of the interatomic potential , obey the following sum rules @xcite : @xmath11 these sum - rule correlations between the elements of the @xmath12 diagonal and off - diagonal blocks distinguish the vibrational problems from the electron problems described by similar hamiltonians @xcite ( see sec .
[ s3_e ] ) .
the matrix elements of the hamiltonian ( [ e1_3 ] ) , generally speaking , are functions of all equilibrium atomic positional vectors , @xmath13 .    our aim is to solve the eigenproblem for the hamiltonian ( [ e1_3 ] ) .
this can be easily done for a crystal @xcite but is not so simple for a disordered structure .
the analysis starts with a choice of the vibrational model for a disordered system
. this could be : ( i ) a crystalline lattice of atoms of random masses ( mass disorder ) and/or connected by random springs ( force - constant disorder ) @xcite ; ( ii ) positionally disordered structures created , for example , by random atomic displacements around crystalline sites @xcite , or by other ways , e.g. the bethe lattice or by bond switching @xcite ; ( iii ) atomic structures , created e.g. by molecular dynamics ( classical or _ ab - initio _ ) , in which intrinsic positional disorder defines the force - constant disorder ( see e.g. refs .
the simplest models are from the first class .
analytical approaches can be easily developed for them due to the existence of a well - defined reference crystalline system ( i.e. the same lattice without disorder ) @xcite .
surprisingly at first sight ( see below ) , lattice models with force - constant disorder mimic very well the main features of topologically disordered glasses and even liquids , e.g. the boson peak , localized band tails , main peaks in the density of states and zero - energy singularity in the instantaneous spectrum @xcite .
the known drawback of lattice models is related to the @xmath14-frequency dependence of the width of the dynamical structure factor @xcite as compared to the parabolic behaviour found experimentally ( see e.g. @xcite ) .
this @xmath14-dependence leads to a parabolic frequency scaling of the dynamic structure factor , @xmath15 for @xmath16 , which is in contrast to the frequency - independent @xmath17 found for topologically disordered models @xcite .      in this paper
, we consider as a vibrational model the crystalline f.c.c .
lattice of atoms , in which the nearest neighbours are connected by springs that are unstretched in equilibrium .
the matrix elements of the hamiltonian ( [ e1_3 ] ) in this case have the following form : @xmath18 where @xmath19 is the unit connection vector between atoms @xmath6 and @xmath20 in equilibrium .
the spring constants , @xmath21 , are supposed to be independent of atomic positions and are random values taken from a certain probability distribution , @xmath22 .
different functional forms of @xmath22 have been used for various problems , e.g. a combination of two @xmath23-functions @xcite , ( also used in the rigidity - percolation problem @xcite ) and the normal ( gaussian ) @xcite , box @xcite and hyperbolic @xcite distributions for glassy problems .
below we use the box distribution : @xmath24   \ , \label{e1_2}\ ] ] where the heavyside function @xmath25 and @xmath26 , @xmath27 is the spring constant in the reference crystal and @xmath28 is the full width of the distribution .
the particular form of the peak - shaped distribution is not important for all results presented below from a qualitative point of view , as we have checked by using the normal distribution instead of the box one .
the box distribution also has the advantage of not needing a lower cut - off to ensure positive values of force constants , as in the case for the normal and hyperbolic distributions .    if the distribution given by eq .
( [ e1_2 ] ) is wide enough , i.e. @xmath29 ( @xmath30 below ) , then negative spring constants appear in the system which could result in vibrational instability of the model in equilibrium , i.e. the occurrence of negative eigenvalues .
such a situation is of direct relevance @xcite to an instantaneous mode analysis in liquids @xcite .
the vibrational spectrum of a disordered system is defined by the properties of the hamiltonian ( dynamical ) matrix . in the case of a crystal ,
the hamiltonian matrix is well structured ( e.g. cyclic band diagonal ) . for topologically disordered glasses with conserved local order ,
the hamiltonian matrix is structured similarly to that for the crystalline counterpart ( with the same local order ) but now , due to the presence of topological disorder : ( i ) the positions of elements are distributed about their positions in the crystalline matrix ( locational disorder ) and ( ii ) the values of the elements are also distributed because of different interatomic distances .
the first effect is completely ignored in the disordered - lattice models and the total influence of disorder is attributed entirely to the second effect .
the question is , how good is this approximation ?
one criterion relates to similarities in the probability distributions for the values of the matrix elements for both types of models .
we have compared the distributions of the hamiltonian matrix elements for two vibrational models : ( i ) the f.c.c .
disordered - lattice model discussed above and described by the matrix elements @xmath31 obeying eq .
( [ e1_6 ] ) and ( ii ) a topologically disordered model of a single - component glass with predominantly icosahedral order ( ic - glass ) obtained by classical molecular dynamics ( see refs .
@xcite for more detail ) . in the latter model ,
the hamiltonian matrix elements , @xmath32 , are defined by the pairwise interatomic potential @xcite , @xmath33 , its spatial derivatives , @xmath34 and @xmath35 , and the relative atomic positions , @xmath36 -    \delta_{\alpha,\alpha ' }   \frac{v ' ( |\r_{ij}^{(0)}|)}{|\r_{ij}^{(0)}| } \ .
\label{hij_glass}\ ] ] the hamiltonian matrix contains @xmath12 diagonal and off - diagonal blocks , the diagonal ( @xmath37 and @xmath38 ) and off - diagonal ( @xmath39 and @xmath40 ) elements being distinctively distributed . for the diagonal elements in the off - diagonal blocks , it is convenient to compare the traces taken with opposite signs .
they represent the force constants , @xmath41 and @xmath42 , for interactions beween atoms @xmath6 and @xmath20 in lattice and glassy models , respectively . similarly the traces of diagonal blocks
correspond to the total force constants , @xmath43 and @xmath44 , for atom @xmath6 in lattice and glassy models , respectively , where @xmath45 .
the four distributions for both models presented in fig .
[ distr](a)-(d ) show qualitatively the same features .
the distributions of the diagonal elements in the diagonal blocks , @xmath46 and @xmath47 , are both peak - shaped with remarkably comparable peak positions and peak widths ( see fig .  [ distr](a ) ) .
the distributions of the off - diagonal elements in the diagonal blocks are both centred around zero and have comparable widths ( see fig .
[ distr](b ) ) .
the distributions of the elements for the off - diagonal blocks are much more model dependent . the solid line in fig .  [ distr](c )
represents the distribution @xmath48 for all interacting atomic pairs .
the distribution function has a complicated form , the features of which can be understood in terms of the functional forms of the interatomic potential and the atomic pair correlation function .
for example , if nearest - neighbour interactions only are taken into account ( as in the disordered - lattice model ) , then the distribution function @xmath49 has a single peak ( see the inset in fig .
[ distr](c ) ) , which corresponds to the first peak in the pair correlation function @xcite .
the simplest model distribution which can mimic the function @xmath48 ( at least its scale ) is the box distribution ( the dashed line in fig .
[ distr](c ) ) .
the distribution of the off - diagonal elements in the off - diagonal blocks ( see fig .  [ distr](d ) ) has a maximum at zero and is symmetric about it .
the disordered - lattice models are intrinsically anisotropic , so that the distribution function , @xmath50 , can look quite different to @xmath51 .
for example , for the ordered f.c.c .
lattice ( @xmath52 ) , it consists of two @xmath23-functions symmetrically located around zero . with increasing disorder ,
these @xmath23-functions are transformed into two box distributions which merge at @xmath53 , forming the single box distribution shown in fig .
[ distr](d ) by the dashed line . for larger @xmath54 ,
these two boxes overlap around zero , thus resulting in a box - shaped peak in this region , i.e. the distribution function @xmath50 becomes closer in shape to the peak - shaped function @xmath51 .
the above comparison demonstrates that the very simple disordered - lattice model with a box distribution of force constants is able to reproduce correctly the ranges of the distributions of the hamiltonian matrix elements for a topologically disordered glass with minimal short - range order ( i.e. the ic - glass ) .
this leads to the plausible conjecture that the physical origins of different phenomena , such as the boson peak , ioffe - regel crossover and secondary peak in the dynamic structure factor ( analysed below ) , localization - delocalization transition , etc . , are common for topologically disordered glasses and disordered lattices .    the other important point to emerge concerns the possibility of estimating the degree of disorder in disordered - lattice models at which they are able to mimic glassy models .
for small disorder , all the distributions for the disordered - lattice model presented in fig .
[ distr ] are quite narrow and are not comparable in width with the distributions for the glassy model .
they become similar only for sufficiently large disorder ( @xmath55 ) .
below we use the simplest ( box ) distribution of force constants which allows us to analyse the problem analytically and which mimics some of the features of topologically disordered glasses by disordered lattices . a natural extension of the lattice models
is in the use of more realistic force - constant distributions for crystalline counterparts of topologically disordered glasses and a comparison of the resulting properties of such models with the corresponding glasses .
a possible way to mimic disorder in glasses is to introduce thermal ( positional ) disorder in the crystalline counterpart .
we have done this for the ic - glass and its crystalline counterpart , the frank - kasper @xmath56-phase @xcite .
a @xmath56-phase lattice has been heated by means of molecular dynamics and the dynamical matrix for the instantaneous cofiguration has been calculated .
the elements of such a matrix do not exhibit locational disorder and only their values are distributed . using temperature as a control parameter it is possible to obtain the distributions of matrix elements very similar to those ones for the topologically disordered ic - glass ( see fig .  6 in ref .
the vibrational density of states ( vdos ) for instantaneous configurations of the @xmath56-phase , as a consequence of such a similarity , resembles the vdos for the ic - glass @xcite .
this supports the use of lattice models for mimicking topologically disordered structures .
we should mention , however , that the elements of the dynamical matrix ( belonging to the euclidian ensemble @xcite ) for positionally disordered crystalline counterparts are correlated with each other within a certain row .
these correlations are absent for the simple box distributions used below .
the solution of the eigenproblem for the hamiltonian ( [ e1_3 ] ) is known for a crystal , when @xmath57 in eq .
( [ e1_2 ] ) .
the eigenmodes then are plane waves @xmath58 , characterized by the wavevector @xmath59 , the branch number @xmath60 and the eigenenergies @xmath61 .
the vdos , @xmath62 , forms a band ( see the dashed line in fig .
[ vdos ] ) containing van hove singularities , and behaves as @xmath63 around the band edges @xcite .
disorder ( @xmath64 ) brings new features to the spectrum ( see fig .  [ vdos ] ) : ( i ) the van hove singularities are washed out ; ( ii ) a high - frequency tail containing localized states occurs ; ( iii ) extra states appear in the low - energy regime .
the exact solution of the vibrational problem with random spring constants is not known , but a number of approximate methods , both for electrons and vibrations , have been developed @xcite .
we are interested in the global features of the spectrum , except the tail region where the states are localized .
an adequate approach in this case is a mean - field treatment within the coherent potential approximation ( cpa ) @xcite , which we use below .
the cpa is known to be very successful in a description of the electronic and vibrational properties of substitutional alloys and is expected to work very well in the whole energy range , except the far tails containing strongly localized states , where field - theoretical approaches are much better @xcite .
the hamiltonian ( [ e1_3 ] ) can be rewritten in the bond representation , @xmath65 with @xmath66 , which is a convenient form for a subsequent mean - field treatment .
the vdos of disordered systems , @xmath67 , @xmath68 can be defined via the spectral - density operator @xcite , @xmath69 , where the angular brackets , @xmath70 , denote averaging over random spring constants .
the coefficients , @xmath71 , in eq .
( [ e2_7 ] ) have a simple physical meaning @xcite ; they define the weight of the disordered eigenstate @xmath72 in the crystalline state @xmath58 , which is a plane wave of a certain polarization .
the importance of the spectral - density operator in the plane - wave basis is related to the fact that the matrix elements @xmath73 are proportional to the dynamical structure factor @xcite ( see sec .
[ s3_dd ] ) which can be measured experimentally @xcite .
both the vdos , @xmath67 , and the spectral density , @xmath74 , can be found for vector vibrations within the mean - field approximation ( see appendix a ) : @xmath75   \ ,   \label{e2_17}\ ] ] and @xmath76 ^ 2 +   \left [ z''(\varepsilon)\varepsilon_{\k\beta } \right]^2 } \ ,   \label{e2_18}\ ] ] via the crystalline green s function , @xmath77 , and the complex effective mean field , @xmath78 ( the dimensionless effective force constant ) , which can be found from the solution of the self - consistent eq .
( [ e2_13 ] ) .
as seen from eqs .
( [ e2_17])-([e2_18 ] ) , the spectral properties of a disordered lattice , within the cpa , are characterized by the energy dependence of the effective dimensionless force constant ( effective field ) , @xmath79 . in this section ,
we discuss the properties of the effective field obtained by the solution of the self - consistent equation ( [ e2_13 ] ) ( see appendix a ) for different values of disorder and in different energy ranges . the most important regime from the viewpoint of comparison with experimental data for the dynamical structure factor is the low - energy limit ( @xmath80 ) , analyzed below .    for the particular choice of the box - like force - constant probability distribution given by eq .
( [ e1_2 ] ) , the integration in eq .
( [ e2_13 ] ) can be done analytically , resulting in the following expression : @xmath81 = 0 \ , \label{e22_1}\ ] ] with @xmath82 eq .
( [ e22_1 ] ) should be solved with respect to the complex effective field @xmath79 . without disorder ( @xmath57 ) ,
the solution is trivial : @xmath83 and @xmath84 . for finite disorder
, the solution can be found numerically in the general case and analytically in some limiting cases .
the numerical solution of eq .
( [ e22_1 ] ) for the real and imaginary parts of the effective field for different values of disorder in the whole energy range is presented in fig .
[ z_eff ] . as seen from fig .
[ z_eff](a ) , the real part of the effective spring constant varies around its crystalline value @xmath83 , being less than unity in the lower part of the energy band and greater than unity in the upper part of the band .
this is the expected behaviour , because the value of @xmath85 describes the level - repelling effect for the bare crystalline states when disorder is introduced in the system @xcite ( see sec .
[ s3_1a ] ) .
the deviation of @xmath85 from unity increases with increasing disorder , which reflects the more pronounced degree of level repelling in more highly disordered systems .
the imaginary part of the effective field is non - zero and negative only in the band region .
it basically reproduces the shape of the vdos ( eq .  ( [ e2_17 ] ) ) and
describes the widths of the spectral - density peak ( eq .  ( [ e2_18 ] ) ) which is proportional to the inverse life - time of plane waves in disordered structures ( see sec .  [ s3_b ] ) .
the magnitude of @xmath86 increases with increasing disorder , thus indicating the shortening of the plane - wave life - time in strongly disordered lattices . the imaginary part of the effective field approaches zero in the low - energy limit , @xmath87 , only if the disorder is low enough , for @xmath88 , when the structure in equilibrium is mechanically stable within the cpa ( i.e. there are no negative eigenvalues ) .
using the energy dependence of the effective field in the low - energy regime ( see appendix b ) , the expression ( [ e2_17 ] ) for the vdos can be rewritten as @xmath89 as follows from eq .
( [ e22_7 ] ) , the low - energy vdos in systems with lattice disorder functionally behaves the same as in the reference crystal , i.e. @xmath90 ( see fig .  [ vdos_log ] ) , and differs from the crystalline vdos just by the factor @xmath91 .
in other words , extra states , in addition to the crystalline ones , appear in the low - energy regime due to disorder - induced level - repelling effects @xcite .
the relative density of these extra states is : @xmath92^{-3/2 } - 1   \   .
\label{e22_8}\ ] ] with increasing energy , the disordered vdos increases in energy faster than @xmath93 and this results in the boson peak , the origin of which has been discussed in ref .
@xcite .    in the low - disorder limit
( see appendix b ) , @xmath94 , the extra density of states ( see eq .
( [ e22_8 ] ) ) in the low - energy range scales quadratically with disorder , @xmath95 which has been confirmed numerically ( see fig .
[ vdos_0_vs_delta ] ) .
all the derivations presented above are valid only for sufficiently small disorder , @xmath96 .
this restriction comes from eq .
( [ e22_5 ] ) ( see appendix b ) , which has solutions only for @xmath88 ( see fig .
[ delta_critical ] ) .
the critical value of disorder can be found numerically , giving @xmath97 .
the corresponding value of the real part of the effective field at criticality , @xmath98 , can be found analytically using the fact that @xmath99 ( see eq .
( [ e22_5 ] ) ) , so that @xmath100 eq .  ( [ e22_12 ] )
gives @xmath101 for the f.c.c .
lattice ( with @xmath102 for @xmath103 and @xmath104 ) .    when the disorder approaches the critical value where the vdos behaves as a renormalized debye function ( see eq .
( [ e22_7 ] ) ) and the imaginary part of the effective force constant behaves as @xmath105 ( see eq .  ( [ e22_6 ] ) ) , the low - energy range shrinks to zero ( see fig .
[ vdos_log ] ) . above the critical disorder , both the imaginary part of the effective field and the effective vdos
are finite at zero energy , @xmath106 and @xmath107 , and negative eigenvalues appear within the cpa .
this means that the system becomes vibrationally unstable in equilibrium , even within the mean - field description .
such a situation corresponds to the spectrum of the dynamical matrix for instantaneous metastable configurations in liquid and glassy states .
it turns out that the energy spectrum of the instantaneous dynamical matrix in the unstable regime ( @xmath108 ) exhibits a peculiar ( singular ) universal behaviour around zero energy , which we have also found in numerical experiments on liquid and vitreous sio@xmath109 and the ic - liquid and ic - glass @xcite .
the spectral density in the plane - wave representation , @xmath110 , provides information about the distribution of contributions of the different disordered eigenstates in a particular plane wave from the branch @xmath60 and with wavevector @xmath59 . in a crystal ,
the spectral density is obviously a @xmath23-function , @xmath111 .
disorder introduced in the system transforms @xmath112 into the disordered spectral density which has a shape given by eq .
( [ e2_18 ] ) .
the shape of @xmath113 depends on the energy behaviour of the effective field @xmath79 .
if we assume that the effective field is energy independent , then the spectral density has the shape of a lorentzian peak , the width of which is proportional to the imaginary part , @xmath86 , of the effective field . of course
, the effective field does depend on energy ( see fig .  [ z_eff ] ) , but if the imaginary part of the effective field is small enough , which is true in the whole energy range for sufficiently small disorder , and at least in the low - frequency regime ( @xmath114 , see below ) at any disorder below the critical one ( @xmath115 ) , then the disordered spectral density still has the shape of a narrow peak ( see fig .  [ sd_kpm ] ) .
the position , @xmath116 , and full width at half - maximum , @xmath117 , of such a spectral - density peak can be estimated from the following equations : @xmath118 for a weak energy dependence of the real part of the effective field ( see fig .
[ z_eff](a ) ) , we can approximate eq .
( [ e3_1 ] ) giving : @xmath119 in the low - energy regime ( @xmath120 ) , the imaginary part of the effective field is @xmath121 according to eq .
( [ e22_6 ] ) , so that @xmath122 this equation corresponds to the well - known rayleigh law for the peak width , @xmath123 , in the @xmath124-representation @xcite : @xmath125 or @xmath126 in the three - dimensional case ( see fig .
[ ir_new ] ) .
( [ e3_2 ] ) defines the dispersion law in disordered lattices , @xmath127 ( for @xmath128 ) , which is basically the same as in the reference crystal except for the renormalized sound velocity , @xmath129 , with @xmath130 being the sound velocity for branch @xmath60 in the reference crystal .
bearing in mind that @xmath131 ( fig .
[ z_eff](a ) ) , it is evident that @xmath132 ( see the inset in fig .  [ ir_new ] ) .
we have found a similar relation for the model of the topologically disordered ic - glass , where @xmath133 @xcite .
this value corresponds to @xmath134 ( see the inset in fig .
[ ir_new ] ) , for which the distribution of the diagonal elements in the diagonal blocks for the disordered - lattice model fits very well the similar distribution for the ic - glass ( see fig .
[ distr](a ) ) .      as seen from eq .
( [ e3_2 ] ) , disorder results in a shift of the crystalline spectral density , characterized by the real part of the effective field , @xmath85 .
the bare level moves downwards ( upwards ) if @xmath135 ( @xmath136 ) , i.e. when it is located in the lower ( upper ) part of the band ( see fig .
[ z_eff](a ) ) .
the shift of the bare level , @xmath61 , by disorder is a consequence of level - repelling effects @xcite .
indeed , disorder introduced in the system results in the appearance of non - zero interaction matrix elements between the bare states .
these interactions lead to standard quantum repelling between levels @xcite . if the interaction matrix elements are approximately the same for all levels , then the bare crystalline level , @xmath61 , from the lower part of the energy band is shifted downwards just because the number of the states above the bare level ( and repelling it downwards ) is larger than the number of states below it ( and repelling it upwards ) .
this explains qualitatively why the value of the real part of the effective field is @xmath135 in the lower part of the band and _ vice versa _ , @xmath136 in the upper half of the energy band ( see fig .
[ z_eff](a ) ) .
we have proved quantitatively ( analytically and numerically ) such a picture of level - repelling effects in the vdos for the force - constant - disordered f.c.c .
lattice @xcite and in the topologically disordered ic - glass @xcite .
a similar effect has also been found numerically in models of disordered si @xcite .
an example of the spectral densities for different energies of external plane waves at relatively high disorder ( @xmath137 ) is shown in fig .
[ sd_kpm ] .
their widths at sufficiently large @xmath61 become comparable to their positions and plane waves at such energies are very short - lived quasiparticles @xcite .
the level - repelling effects are clearly seen in pronounced shifts ( to lower energies ) of the peak positions with respect to @xmath61 .
we have also presented in fig .
[ sd_kpm ] the results of precise numerical solutions for the spectral densities obtained by the the kernel polynomial method ( kpm ) @xcite .
the very good agreement between the cpa and kpm results supports the reliability of the mean - field treatment in obtaining the spectral densities .
disorder also broadens a bare level , so that the @xmath23-functions making up the crystalline vdos are broadened into peaks of finite width , related to the imaginary part of the effective field ( see eq .
( [ e3_1 ] ) ) .
this occurs because a plane wave @xmath58 is not an eigenstate in the disordered lattice .
so it decays with time and can be treated as a quasiparticle having a finite life - time .
this life - time is inversely proportional to the width of the peak of the spectral density and this determines the physical meaning of @xmath138 @xcite .
the scattering properties of external plane waves in disordered lattices strongly depend on their bare energy .
if the energy is low enough , the spectral - density peak is narrow compared to its position , and the quasiparticles live for a long time , thus indicating a weak - scattering regime . for higher energies ,
the width of the spectral density becomes comparable to the peak position ( see fig .
[ ir_new ] ) .
this means that the quasiparticles there decay quickly and the bare plane waves are in the region of strong scattering .
the crossover between weak- and strong - scattering regimes is called the ioffe - regel crossover @xcite .
the energy of the bare plane wave , @xmath139 , at which such a crossover occurs is the ioffe - regel crossover energy .
the quantitative condition for the ioffe - regel crossover is @xcite , @xmath140 with @xmath141 standing for the frequency and @xmath142 being the life - time of the quasiparticle .
the quasiparticle life - time can be easily evaluated in the weak - scattering regime from the time dependence of the probability , @xmath143 , to find the system at time @xmath144 in the same plane - wave state @xmath58 as it was at @xmath145 @xcite , i.e. @xmath146 , where @xmath147 the inverse life - time , @xmath148 , of course , coincides with the full width , @xmath149 , of the spectral density in the @xmath124-representation @xcite , @xmath150 , and the ioffe - regel criterion ( [ e3_c_1 ] ) can therefore be rewritten as @xmath151 this equation should be solved with respect to @xmath61 ( or @xmath152 ) and an estimate , @xmath153 ( or @xmath154 ) , for the ioffe - regel crossover energy ( or frequency ) can be found .
such a solution can exist at sufficiently large @xmath155 .
indeed , in the low - frequency limit , the width of the spectral density , @xmath156 , is much smaller than the peak position , @xmath157 , but @xmath158 grows quickly with increasing @xmath159 and can easily reach the peak - position value , if the disorder is large enough ( see fig .
[ ir_new ] ) .    we have solved eq .
( [ e3_c_8 ] ) numerically using the exact energy dependence for the effective field for different values of disorder , @xmath160 ( see fig .  [ ir_new ] ) .
the results for the ioffe - regel crossover energy are presented in the phase diagram for vector vibrations in the f.c.c .
lattice ( see fig .  [ phase_diagram_new ] ) .
there is a minimum value of disorder ( @xmath161 for the disordered f.c.c .
lattice using eq .
( [ e3_c_8 ] ) ) , starting from which the region of the strong - scattering regime appears in the system .
the strong - scattering regime becomes broader with increasing disorder , while the low - energy weak - scattering regime shrinks basically to zero around the critical disorder , @xmath162 .
we should stress that the boundary between weak- and strong - scattering regimes should really be thought of as a relatively wide crossover strip , and eq .
( [ e3_c_8 ] ) serves for an order - of - magnitude estimate of the location of this strip . bearing this in mind , the high - frequency region at sufficiently large disorder
should not properly be regarded as a region of weak - scattering because the width of the spectral - density peak there is of the same order of magnitude as the peak position .
the true weak - scattering regime occurs only either in the low - energy regime at any disorder less than the critical value , @xmath163 , or for all energies at small disorder ( @xmath164 ) .
the existence in these regions of small values of the energy and disorder parameters forces the spectral - density width to be much smaller than the peak position . in order to illustrate
how sensitive the ioffe - regel crossover energy is to the numerical coefficients in eq .
( [ e3_c_8 ] ) , we have calculated @xmath165 from eq .
( [ e3_c_8 ] ) in which the spectral - density peak width , @xmath166 , has been replaced , say , by @xmath167 , thus enhancing the strong scattering .
the results for this ioffe - regel crossover boundary are shown by the dot - dashed curve in fig .
[ phase_diagram_new ] .
it is clearly seen how the region of the strong - scattering regime increases , with just the low - energy and low - disorder regions surviving for the weak - scattering regime .
we have also shown in fig .
[ phase_diagram_new ] the trajectory of the boson peak ( circled solid line ) calculated as in ref .
remarkably , it almost coincides with the trajectory of the ioffe - regel crossover for large enough disorder when the strong - scattering regime appears in the system .
this feature has also been found in real glasses , e.g. v - sio@xmath109 @xcite and in the topologically disordered ic - glass model ( to be published elsewhere ) .
the solid line in the phase diagram shows the upper band edge of the vdos calculated within the cpa .
this appears to be a very good estimate for the trajectory of the threshold for the localization - delocalization transition , which , in fact , we have calculated @xcite more exactly by means of a multi - fractal analysis @xcite with a relatively high precision , @xmath168 to be situated at slightly lower energies than the cpa band edge .
the cpa fails to reproduce the high - energy tail in the vdos ( see fig .
[ vdos ] ) and more accurate methods should be used @xcite .
it should also be mentioned that , for @xmath54 , a low - energy tail of localized states extending into the negative eigenvalue range appears in the true spectrum @xcite .
the cpa is not able to reproduce this effect for @xmath169 but still , for the majority of states with energies far enough away from both localization thresholds , the cpa results are very reliable .      in the previous subsections
, we have analysed the situation at small and intermediate energies of external plane waves .
these energies correspond to a relatively small momentum transfer in scattering experiments ( see e.g. @xcite ) .
however , larger values of @xmath170 are also of interest both from experimental ( see e.g. @xcite ) and theoretical ( see e.g. @xcite ) points of view .
attention has been mainly focused on the behaviour of the dynamic structure factor , @xmath171 , which can be roughly approximated in the high - temperature regime ( @xmath172 ) by the following expression @xcite : @xmath173 where the debye - waller factor has been ignored ( cf .
@xcite ) . in eq .
( [ e3_dd_1 ] ) , we have also neglected the mixture of polarizations for phonons @xmath174 and assume that @xmath60 is related to the longitudinal branch . in this case ,
the dynamical structure factor is obviously ( see eq .
( [ e2_7 ] ) ) proportional to the spectral density , @xmath175 where @xmath176 is the spectral density in the @xmath124-representation .
a peculiar feature has been found in the frequency dependence of the dynamical structure factor in topologically disordered glasses at large values of the momentum transfer .
namely , the appearance of a secondary peak below the brillouin peak in the @xmath124-dependence of @xmath171 has been observed @xcite with increasing value of @xmath170 .
this has been related to the boson peak @xcite , or interpreted as a manifistation of a microscopic relaxation process typical of a topologically disordered glass @xcite .
a natural question is whether or not the secondary peak exists for the disordered - lattice models and , if so , then what is the origin of this peak . in order to answer these questions , first , we plot in fig .
[ sd_1](a ) the evolution of the spectral density @xmath110 for the longitudinal branch with increasing value of the wavevector , taken for definitness to be along the @xmath177 $ ] direction .
it is clearly seen from this figure that the spectral density has a pronounced peak ( brillouin peak ) for all values of @xmath170 taken from the first brillouin zone .
however this peak is no longer lorentzian for large @xmath170 .
it is asymmetric and has a pronounced shoulder in the low - energy range .
a similar picture has been observed in simulations for topologically disordered systems .
the brillouin peak shifts to higher energies , becomes very broad and a smooth shoulder - like feature appears below this peak and may also contain broad peak - shaped features @xcite , so that the spectral density starts to resemble the vdos . in the low - energy regime , the spectral density goes to zero .
the dynamical structure factor differs from the spectral density and contains an additional decaying energy ( frequency ) function ( see eq .
( [ e3_dd_2 ] ) . therefore ,
if we plot the effective dynamic structure factor , @xmath178 , versus frequency ( see fig .  [ sd_1](b ) ) , the low - frequency features are enhanced due to the @xmath179 contribution . the low - energy shoulder in the spectral density of the disordered - lattice models transforms , for large @xmath170 values close to the brillouin - zone boundary , to the secondary peak in the dynamic structure factor ( see the solid line in fig .  [ sd_1](b ) ) .
the intensity of the secondary peak depends on the degree of disorder , so that it becomes very pronounced around the critical disorder ( see fig .  [ sd_delta](b ) ) , although it does not appear for the spectral density in the energy representation ( cf fig .
[ sd_delta](a ) ) .
therefore , the presence of the secondary peak in @xmath171 of disordered - lattice models for large values of wavevector and large degrees of disorder is a distinctive feature of such models .
consequently , the occurrence of a similar peak in topologically disordered structures can hardly be attributed solely to topological disorder ( cf . ref .
@xcite ) .    the origin of the secondary peak in the dynamic structure factor of disordered - lattice models can be understood from the qualitative analysis of eqs .
( [ e2_18 ] ) and ( [ e3_dd_2 ] ) . in the low - frequency limit , @xmath180 ( or @xmath80 ) , the energy dependence of the spectral density in eq .
( [ e2_18 ] ) is dictated by the imaginary part of the effective field , @xmath181 , so that @xmath182 .
the low - energy limiting behaviour of @xmath138 deviates approximately around the energy corresponding to the shifted ( due to the level - repelling effect ) lowest van hove singularity , i.e. the first minimum in the energy dependence of @xmath183 ( cf . figs .
[ z_eff](a ) and ( b ) ) . above this energy ,
the energy dependence of the imaginary part of the effective field is appreciably weaker and basically reproduces the shape of the vdos ( see fig .
[ z_eff](b ) ) .
this results in a weaker dependence of the dynamic structure factor , which shows , in this intermediate frequency range , either the shoulder - like behaviour , if the wavevector magnitude is not large enough and the low - frequency tail of the brillouin peak plays a significant role ( see the long and short dashed lines in fig .  [ sd_1 ] ) , or a secondary peak for large values of @xmath170 ( see the solid line in fig .  [ sd_1 ] ) , when the influence of the brillouin peak is suppressed by the @xmath179 function in eq .
( [ e3_dd_2 ] ) .
therefore , the position of the secondary peak in @xmath171 for the disordered lattice models approximately coincides with the shifted lowest van hove singularity and thus with the position of the boson peak @xcite .
however , we should stress that the presence of the pronounced secondary peak in the dynamic structure factor for large values of @xmath170 does not mean a major contribution of plane waves with such @xmath170 to the boson peak , which is proportional to the sum of all plane waves ( see ref .
@xcite for more detail ) .
the @xmath0-dependence of the dynamic structure factor in the low - frequency regime found for the disordered - lattice models is consistent with a @xmath14-dependence of the spectral density width and is in contrast to the constant ( non - zero ) , limit of @xmath171 found for topologically disordered glasses @xcite , and thus can be considered as a drawback of lattice models .      we have seen in the previous subsection that , independent of the degree of disorder , low - energy external plane waves are weakly scattered in a force - constant disordered lattice .
this is a consequence of the essentially different dependences on the plane - wave energy of the peak width ( @xmath184 ) and peak position ( @xmath185 ) of the spectral density ( see eqs .
( [ e3_2],[e3_3 ] ) ) , which gives rise to the occurrence of weak scattering ( @xmath186 ) for @xmath128 .
the question is how generic is this feature .    in order to answer this question ,
let us consider a more general formulation of the problem using the following hamiltonian in the site representation : @xmath187 the diagonal elements of the hamiltonian have two constituents : random on - site  energies  , @xmath188 , and the correlation contribution proportional to the sum of random off - diagonal elements , @xmath21 , in the same row .
the parameter @xmath189 controls the sum - rule correlations between diagonal and off - diagonal elements . if @xmath190 and @xmath191 , the hamiltonian ( [ e3_e_1 ] ) describes the scalar vibrational problem , while the case @xmath192 corresponds to the well - known electron anderson hamiltonian with on - site and off - diagonal disorder @xcite .
the spectral properties of a similar hamiltonian , but defined on a positionally random set of sites with deterministic transfer integrals , have been studied in ref .
@xcite .
the scattering properties of external plane waves in the disordered structure described by the hamiltonian ( [ e3_e_1 ] ) can be studied using spectral densities for which the mean - field approach works very well ( as shown above for the vibrational problem ) .
it is easy to show that only in the case of @xmath190 and @xmath191 can the hamiltonian ( [ e3_e_1 ] ) be reduced to the form of eq .
( [ e1_8 ] ) which has a multiplicative effective mean field ( the energy of the quasiparticles is the product of the effective field and the bare energy , @xmath193 - see eq .
( [ e2_13_3 ] ) ) . for @xmath194 ,
the situation is quite different , because now the effective field has an additive constituent as well , and this drastically changes the scattering properties of plane waves .
below , we demonstrate this , considering the simplest particular case for the hamiltonian ( [ e3_e_1 ] ) characterized by @xmath192 , @xmath195 and random @xmath188 taken from the box distribution of half - width , @xmath196 , centred around @xmath197 , i.e. the standard anderson hamiltonian with diagonal disorder @xcite , for the simple cubic lattice .
the cpa spectral density , @xmath198 , for the anderson hamiltonian can be expressed via the quasiparticle energy , @xmath199 , which is now the sum of the crystalline energy , @xmath200 , and the effective field , @xmath201 .
the spectral densities @xmath202 are peak shaped and in the low - energy limit , @xmath203 ( the energy @xmath204 is referred to the lower crystalline band edge , @xmath205 , where @xmath206 is the coordination number ) , the position and width of the peak can be estimated from the following relations : @xmath207 and @xmath208 where the peak position , @xmath209 , is measured from the cpa lower band boundary , @xmath210 ( see fig .  [ ir_electrons ] ) .
the last expression is a consequence of the fact that @xmath211 \simeq   -\pi    \chi^{\rm el}\sqrt{\varepsilon -   \varepsilon_{\rm min}^{\rm cryst } } $ ] , if @xmath212 with @xmath213 being the constant coefficient similar to the coefficient @xmath214 in eq .
( [ e22_1a ] ) .    in full analogy to vibrations ( see eq .
( [ e3_c_8 ] ) ) , the ioffe - regel crossover between weak- and strong - scattering regimes for electrons can be defined in terms of the following condition : @xmath215 in contrast to vibrations , the spectral - density peak width for electrons scales more slowly than the peak position as a function of the plane - wave energy ( see eqs .
( [ e3_e_6])-([e3_e_7 ] ) and fig .  [ ir_electrons ] ) .
this means that , in the long - wavelength limit , @xmath216 , the width of the spectral density is always larger than the peak position , and the strong - scattering regime occurs for plane waves having energies around the lower crystalline band edge ( see fig .
[ phase_diagram_electrons ] ) .
the situation is opposite for vibrations , for which @xmath217 ( cf .
( [ e3_3 ] ) and ( [ e3_e_7 ] ) ) . as a consequence ,
the shapes of the regions associated with the weak- and strong - scattering regimes for plane waves in the phase diagram are different for electrons and vibrations ( cf . figs .  [ phase_diagram_new ] and [ phase_diagram_electrons ] ) .
now , for electron plane waves with energies around the lower crystalline band edge , the scattering is always strong and the weak - scattering regime is realized only in the midband range for sufficiently small disorder .
in contrast to vibrations , only a single parameter , the strength of disorder , @xmath196 , rather than two parameters , the degree of disorder @xmath160 and the energy , @xmath61 , for vibrations , controls the appearance of the weak - scattering regime . the boundary between the weak- and strong - scattering regimes ( see the dashed line in fig .  [ phase_diagram_electrons ] ) is again not precise and should be thought of as an estimate for the crossover energy .
it should be mentioned that the electron phase diagram presented in fig .
[ phase_diagram_electrons ] , for the case of very small disorder ( @xmath218 ) , supports a standard concept for electrons @xcite , according to which the ioffe - regel crossover from the weak- to the strong - scattering regime corresponds to the transition from extended states to localized ones .
indeed , the ioffe - regel crossover boundary ( the dashed line in fig .  [ phase_diagram_electrons ] ) almost coincides with the localization - delocalization threshold ( the circles in fig .  [ phase_diagram_electrons ] )
@xcite and the cpa band edge ( solid line ) as @xmath218 . at high disorder ,
the localization - delocalization transition , however , is not associated with the ioffe - regel crossover , because all plane waves for such disorder are strongly scattered .    in the more general case of diagonal and off - diagonal disorder coexisting in the anderson hamiltonian ( but with @xmath192 ) , the most successful homomorphic cluster cpa @xcite results in an effective field which necessarily has an additive constituent and the above conclusions still hold . in the case of @xmath219 , the mean field
can be constructed as well but it also has an additive constituent which disappears only for @xmath190 ( to be published elsewhere ) .
thus , we can conclude that the existence of the low - energy weak - scattering regime is characteristic only of the vibrational problem and it results from the exact sum - rule correlations in the corresponding hamiltonian .
we have analytically investigated classical harmonic atomic vibrations in disordered lattices . the f.c.c .
atomic lattice with force - constant disorder described by a box probability distribution has been chosen for analysis .
the distributions of matrix elements in diagonal and off - diagonal blocks in the dynamical matrix are rather similar to those of topologically disordered models .
the atomic vibrations were treated as vectors and analysed within the framework of the single - bond coherent potential approximation ( cpa ) ; the results are essentially the same as those found from a precise numerical study .
the cpa self - consistency equation has been derived for vector models and solved analytically in the limiting cases of small disorder and low energies .
the cpa solution of the problem has allowed us to investigate the spectral properties of the model in terms of spectral densities .
we have shown that the rayleigh law for the spectral - density width is intrinsic for vector vibrations in disordered lattices .
we have also investigated the ioffe - regel crossover between weak and strong scattering for vector vibrations .
the regime of weak scattering takes place at all energies for small degrees of disorder , but in the low - energy regime occurs only for large disorder .
the existence of a disorder - independent weak - scattering regime for low - energy external plane waves is a consequence of the sum - rule correlations in the dynamical matrix , typical for the vibrational problem only .
we have demonstrated that the ioffe - regel crossover occurs in the boson - peak region and is not related to the localization - delocalization transition for vector vibrations in the disordered lattice studied .
in contrast , the ioffe - regel crossover for electrons ( described by the anderson model with on - site energy disorder in the simple cubic lattice ) exists only at low values of disorder , where it _ is _ associated with the localization - delocalization transition . at higher values of disorder for the anderson model , the transition is between strongly - scattered extended states and localized states .
we have found that a secondary peak appears in the dynamic structure factor below the brillouin peak for sufficiently large values of @xmath170 for the lattice models as previously found in simulations of topologically disordered glasses .
hence this feature can not be due solely to the presence of topological disorder , nor can it be due to relaxations .
here we give the derivation of the self - consistency cpa equation for vector vibrations .
it is known @xcite that the self - consistency equation for the effective force constant , @xmath220 , in the single - bond approximation can be written as follows : @xmath221 where @xmath222 and the effective green s function , @xmath223^{-1 } $ ] , is defined for the effective hamiltonian , @xmath224 .
our aim is to find the matrix element , @xmath225 , for vector vibrations .
this can be rewritten as @xmath226 and this can then be expressed in terms of diagonal elements of the green s function in the site basis for the reference crystal . in order to show this ,
let us first rewrite the matrix elements , @xmath227 in the plane - wave basis , @xmath228 , @xmath229 where @xmath230 is a @xmath231-dimensional unit polarization vector for the branch @xmath60 and the effective energy , @xmath232 , obeys the relation : @xmath233 with @xmath234 being the dimensionless effective force constant and the index @xmath235 running over all nearest neighbours to the site @xmath6 , i.e. @xmath236 , with @xmath206 standing for the number of nearest neighbours ( all the sites in a crystal are equivalent and @xmath61 , in fact , does not depend on @xmath6 ) .    substituting expression ( [ e2_13_5 ] ) back into eq .
( [ e2_13_1 ] ) we obtain : @xmath237 where @xmath238 is the average on - diagonal ( independent of @xmath6 ) element of the effective green s function in the site basis .
the effective green s function , @xmath239 , can be expressed via the crystalline one , @xmath240 , as @xmath241 finally , the expression for @xmath242 reads : @xmath243 \ ,        \label{e2_14}\ ] ] and this solves the problem when substituted back into eq .
( [ e2_13 ] ) .
the dimensionless effective spring constant @xmath244 fully determines the effective mean field .
this is a complex quantity and should be found from the solution of the self - consistent equation ( [ e2_13 ] ) .
the self - consistent equation ( [ e2_13 ] ) for vector vibrations is very similar to that for the scalar models @xcite .
the only difference is the factor @xmath231 in the denominator of eq .
( [ e2_13 ] ) , with @xmath225 given by eq .
( [ e2_14 ] ) . for scalar models ,
this is absent , even for @xmath231-dimensional lattices with @xmath245 .
the crystalline green s function , @xmath246 , with complex argument @xmath247 , can be easily found as an analytical continuation from the real axis of its imaginary part ( vdos ) @xcite .
after that , the self - consistent equation can be solved numerically in the general case , and analytically in some limiting cases ( e.g. @xmath80 and/or @xmath94 ) , and the effective field can be found .
if @xmath79 is known , the task of evaluating the spectral density becomes straightforward .
indeed , bearing in mind that the energy , @xmath248 , of the eigenstates in the effective crystal is connected with the energy , @xmath61 , of the similar eigenstates in the reference crystal according to relation ( [ e2_13_3 ] ) , it is easy to see that the vdos and diagonal elements of the spectral - density operator in the crystalline plane - wave basis satisfy eqs .
( [ e2_17])-([e2_18 ] ) .
here we show how the self - consistent eq .
( [ e22_1 ] ) can be solved analytically in the low - energy limit , @xmath80 , and how the effective field , @xmath79 , can be found in this case . in order to do this
, first we find the expression for the function @xmath249 given by eq .
( [ e22_1aa ] ) via @xmath79 and @xmath250 .
this can be done , bearing in mind the following relations for the crystalline green s function , @xmath251 , in the low - energy limit , when the debye approximation for the crystalline vdos , @xmath252 is valid : @xmath253 \simeq   \mbox{const}\ \ \ \ \mbox{and } \ \ \ \ \mbox{im } \left [        g^{\rm cryst }        \left(\frac{\varepsilon}{z }        \right ) \right ]   \simeq   -\pi\chi_{\rm deb } \left(\frac{\varepsilon}{z ' } \right)^{\frac{d-2}{2 } } \ .
\label{e22_2}\ ] ] substituting eq .
( [ e22_2 ] ) into eqs .
( [ e22_1aa ] ) and ( [ e2_14 ] ) , we obtain the leading terms in @xmath250 for @xmath8 : @xmath254 and @xmath255   \ .
\label{e22_4}\ ] ] these expressions can then be used for an expansion of the logarithm in eq .
( [ e22_1 ] ) in terms of @xmath250 , which , after some algebra , results in the following final relations for @xmath183 and @xmath138 : @xmath256 = 0 \ ,    \label{e22_5}\ ] ] and @xmath257^{-1 } -        \left [            ( \gamma+1)(z')^2       \right]^{-1 }   \right\}^{-1 }   \ ,   \label{e22_6}\ ] ] where @xmath258 .    in fact ,
relation ( [ e22_5 ] ) is the equation to find the real part of the effective force constant , @xmath85 .
the important point is that all the parameters entering this equation are energy independent , so that @xmath259 .
( [ e22_5 ] ) can be solved numerically in the case of arbitrary disorder and analytically in the small - disorder limit .    the imaginary part of the effective force constant , as follows from eq .
( [ e22_6 ] ) , scales with energy as @xmath260 .
such a behaviour is basically a consequence of the debye law for the vdos in the low - energy regime ( see eq .
( [ e22_1a ] ) ) and , in this sense , is general for disordered lattices .
this dependence also determines the scaling behaviour of the spectral - density peak width with energy , the well - known rayleigh law in this case .
( [ e22_5])-([e22_6 ] ) have important consequences for the energy dependence of the vdos in the low - energy limit .
indeed , in the limit @xmath80 , the expression ( [ e2_17 ] ) for the vdos can be rewritten in the form of eq .
( [ e22_7 ] ) .    in the low - disorder limit , @xmath94 , eq .  ( [ e22_5 ] ) can be solved with respect to @xmath85 while eq .
( [ e22_6 ] ) can be simplified , so that @xmath261 @xmath262 these relations are used in deriving eq .
( [ e22_11 ] ) .