during the past ten years , the evolution of structural relaxation in glass - forming liquids has been intensively studied using neutron - scattering spectroscopy , various light - scattering techniques , dielectric - loss spectroscopy , and molecular - dynamics simulation .
results of this work have also been used to test the mode - coupling theory ( mct ) , which interpretes the structural relaxation as precursor of the glass - transition .
originally , the mct was proposed as an approximation approach for the cage effect in liquids @xcite . in its simplest version ,
the mct equations of motion describe an ideal liquid - to - glass transition , i.e. a bifurcation from ergodic to nonergodic dynamics , if control parameters like temperature @xmath3 or packing fraction @xmath4 cross critical values @xmath5 or @xmath6 , respectively .
this bifurcation is connected with the evolution of a two - step relaxation scenario entirely determined by the regularly changing equilibrium structure .
two divergent time scales appear , closely connected to two power - law decay processes .
a detailed description of these results can be found in ref . @xcite and references therein .
comparisons of the theoretical results for simple model systems with experiments done on colloids @xcite , and with computer - simulation studies @xcite demonstrate the validity of the microscopic mct approach .
for the solutions of the mct equations , a variety of results has been derived by asymptotic expansions , using as a small parameter the distance from the critical point , @xmath7 , or @xmath8 , respectively .
the leading - order results of this expansion establish universality features of the mct dynamics .
assessments of the theory have been reached by comparing spectra in the ghz regime or relaxation curves within the pico - second window with the universal results .
the outcome of this work , which is reviewed in ref .
@xcite , leads to the conclusion that mct properly describes some essential features of structural relaxation even for some complicated molecular liquids .
the mct for simple systems has been extended recently to liquids of nonspherical molecules @xcite .
but so far , only the bifurcation equation for the so - called nonergodicity parameters resulting within the new theory could be solved . comparing these results with the findings of molecular - dynamics simulations for a liquid of linear molecules @xcite and for water
@xcite indicates that the mct for molecular liquids is promising .
it was also predicted that there can be two states of nonergodic motion for nonspherical molecules .
these states are connected by a type - a transition if the molecules exhibit a top - down symmetry @xcite . such transitions are generic possibilities in mct , provided there is some symmetry in the problem rendering certain mode - coupling coefficients zero @xcite . at a type - a transition ,
the nonergodicity parameters change continuously , whereas at the conventional mct transition , referred to in this context as a type - b transition , a discontinuity occurs @xcite .    in this paper , correlation functions and susceptibility spectra
shall be discussed , which deal with the glassy dynamics of the orientational degrees of freedom of nonspherical molecules .
the results are obtained as solutions of the equations of motion derived previously @xcite for the dynamics of a linear probe molecule immersed in a simple liquid .
a top - down symmetric dumbbell of two fused hard spheres will be considered as the molecule , and as the solute , a hard - sphere system is chosen .
this model deals with the simplest problem concerning glassy rotational dynamics , namely the influence of the cages formed by the neighbors of the molecule on the molecule s reorientational motion as it is caused by steric hindrance .
the dynamics will be exemplified for two cases : a molecule with a large elongation and a molecule with a small elongation .
it will be shown that large elongations lead to strong coupling of the rotational degrees of freedom to the density fluctuations of the solute , such that the glassy dynamics of the latter enforces the validity of all the universal mct laws for the solvent . moreover ,
the corrections to the leading - order - asymptotic laws show the same qualitative trends as studied for simple liquids @xcite .
a motivation of the present study is the explanation of three general properties of the @xmath0-relaxation in molecular liquids , which are exhibited in fig .
[ ersteabb ] . in this figure ,
experimental susceptibility spectra for the van - der - waals liquid propylene carbonate ( pc ) are reproduced for four temperatures .
one set of data deals with the response for angular - momentum index @xmath1=1 ; it was obtained by dielectric - loss spectroscopy @xcite .
the other set was measured by depolarized - light - scattering spectroscopy @xcite and deals with the ( @xmath1=2 ) reorientational dynamics .
the data show for @xmath9k and @xmath10k @xmath0-relaxation peaks at @xmath11ghz ( @xmath1=1 ) and @xmath12ghz ( @xmath1=2 ) respectively .
these temperatures exceed the melting temperature @xmath13k of pc by more than @xmath14k .
lowering @xmath3 to @xmath15k , the @xmath0-peaks of the spectra are shifted down by about two orders of magnitude .
the shape of the @xmath0-peak is temperature independent , and the ratio of the @xmath0-process - time scales , characterizing the @xmath0-peak - maximum positions for the two values of @xmath1 , is also @xmath3-independent .
these are two features which mct predicts to be universal .
the first nonuniversal feature to be understood is , that the @xmath0-peak intensity , taken relative to that of the band of microscopic excitations at around 1thz , is larger for the ( @xmath1=1 ) response than for the ( @xmath1=2 ) case : the former exceeds the latter by about a factor of @xmath16 .
second , the ( @xmath1=1 ) response is slower than the response for @xmath1=2 : the ratio of the @xmath0-peak positions is about @xmath17 .
third , the @xmath0-peak of the ( @xmath1=1 ) response is less stretched than the peak for @xmath1=2 , i.e. the halfwidth of the ( @xmath1=1 ) peak is smaller than that of the ( @xmath1=2 ) peak .
if one describes these peaks by the spectra of the kohlrausch law , @xmath18 , the stretching exponent @xmath2 for @xmath1=1 , @xmath19 @xcite , is larger than the one for @xmath1=2 , @xmath20 @xcite .
the same three @xmath0-peak features are noticed , if one compares the depolarized - light - scattering spectra of glycerol @xcite with the corresponding dielectric - loss spectra @xcite .
a fourth general feature to be explained is the large ratio of the @xmath0-relaxation - time scale found by depolarized - light - scattering spectroscopy and the one found for the longitudinal elastic modulus by brillouin - scattering spectroscopy . for salol ,
a ratio of about 10 was reported @xcite , while for pc , a factor of about 5 was found @xcite .
the small elongation of concern in this paper is chosen so that it exceeds the critical value for the above - mentioned type - a transition by about @xmath12% . the theory for the corrections to the leading - order asymptotic laws
@xcite implies that these diverge at a type - a transition .
therefore , the range of validity of the universality features of the standard mct bifurcation shrinks upon approaching the type - a transition .
it will be shown that in our example the standard results are not exhibited any more for reasonable choices of the distance parameter @xmath21 .
in particular , it is impossible to identify a two - step scenario for the odd-@xmath22 correlators , nor is @xmath0-relaxation scaling observed .
the paper is organized as follows . in sec .
[ themodelsystem ] , the model is defined , and the mct equations are noted . after an overview of the general scenario for the evolution of the glassy relaxation of the reorientational correlators ( sec .
iiia ) , the differences between the relaxation patterns for the ( @xmath1=1 ) and ( @xmath1=2 ) response are described for strong ( sec .
iiib ) and weak ( sec .
iiic ) steric hindrance . in sec .
iiid it is demonstrated how the @xmath2-relaxation is described by the first scaling law , and in sec .
iiie it is discussed how the @xmath0-relaxation - scaling - law description emerges .
the concluding sec .
iv summarizes the results .
a system of @xmath23 spherical particles shall be considered as the solvent .
the basic variables describing the structure are the density fluctuations for the wave vector @xmath24 : @xmath25 . here
@xmath26 , labels the centers of the particles .
the structure factor @xmath27 provides the simplest information on the equilibrium distribution of the particles ; here , @xmath28 denotes canonical averaging . because of rotational symmetry , @xmath29 depends on the wave - vector modulus @xmath30 only .
the structure factor can be expressed through the direct correlation function @xmath31 via the ornstein - zernicke equation @xmath32 ; where @xmath33 denotes the particle density @xcite .
the simplest quantities , characterizing the structural dynamics in a statistical manner , are the normalized auto - correlation functions for the density fluctuations , called the density correlators @xmath34 .
the evolution with increasing time @xmath35 is given by the canonical equations of motion .
we will also need fourier - laplace transforms for complex frequency @xmath36 , @xmath37 , @xmath38 , using the convention : @xmath39 . for real frequency @xmath40 ,
one gets with @xmath41 : @xmath42 .
the imaginary part @xmath43 is called the fluctuation spectrum , and @xmath44 is the susceptibility spectrum @xcite .
the basic version of mct consists of two equations @xcite .
the first one is exact and derived within the zwanzig - mori formalism :    [ eins ] @xmath45 here , @xmath46 , with @xmath47 denoting the thermal velocity .
the relaxation kernel @xmath48 is a fluctuating - force correlator .
the equation has to be solved with the initial condition @xmath49 @xcite .
equation ( [ einsa ] ) is equivalent to the double fraction @xmath50    the second mct equation is obtained by writing the kernel as a sum of a regular term and a contribution describing the cage effect .
the latter is treated by kawasaki s factorization approximation for the force correlations .
it is found to be a quadratic functional of the density fluctuations : @xmath51 .
for the sake of simplicity , the regular term shall be neglected in the following .
furthermore , the wave - vector modulus will be discretized to @xmath52 values with equal spacing @xmath53 .
thus , @xmath54 , @xmath55 , @xmath56 can be considered as labels running from @xmath57 to @xmath52 . as a result
, the kernel is given as a quadratic mode - coupling polynomial @xmath58 of the @xmath52 correlators @xmath59 , @xmath60 : @xmath61=\sum_{kp } v_{qkp}\,\phi_{k}(t)\,\phi_{p}(t ) .
\label{zwei}\ ] ] the positive coupling coefficients @xmath62 are given by @xmath29 and @xmath31 @xcite .
anticipating these equilibrium quantities to be known , eqs .
( [ einsa ] ) and ( [ zwei ] ) are closed .    equations ( [ eins ] ) and ( [ zwei ] ) exhibit a transition from liquid - state dynamics in the regime @xmath63 or @xmath64 to glass - state dynamics for @xmath65 or @xmath66 . in the former regime
the density fluctuations decay to zero for long times , @xmath67 .
the ideal - glass states exhibit a nontrivial long - time limit , which is called the nonergodicity parameter , @xmath68 .
it is the debye - waller factor of the glass . at the transition ,
this long - time limit is discontinuous , and the jump is called critical nonergodicity parameter or plateau , @xmath69 . at the critical point , the correlators decay algebraically : @xmath70 .
the exponent @xmath71 , @xmath72 , is called the critical exponent , and @xmath73 is denoted as the critical amplitude .
@xmath74 marks the time scale of the transient from the microscopic motion to the relaxation dynamics of the mct .
the mct @xmath0-process is defined as the dynamics for those times , where the correlators of the liquid decay from the plateau @xmath75 to zero .
the mct @xmath2-process deals with the dynamics , where the correlators are near the plateau , i.e. @xmath76 . the first relaxation step of the anomalous dynamics is given by the initial part of the @xmath2-process ; it deals with the decay towards the plateau for times outside the transient : @xmath77 , @xmath78 .
the second step is the @xmath0-decay in the liquid .
its initial part is identical with the final part of the @xmath2-process , and it follows von schweidler s law @xmath79 .
the exponent @xmath80 , @xmath81 , is called the von - schweidler exponent .    in a leading - order expansion in the small parameter @xmath82
one finds the universal results for the @xmath2-process .
there holds the factorization theorem @xmath83 the dependence on time and on control parameters is given by the @xmath54-independent function @xmath84 , which is called the @xmath2-correlator .
it is determined by the equation @xmath85 to be solved with the initial condition @xmath86 .
the number @xmath87 , @xmath88 , is referred to as the exponent parameter .
@xmath89 is a smooth function of the control parameters and is called the separation parameter .
its zero defines the critical point . expanding in leading order in the distance @xmath21
, one can write @xmath90 .    from eq .
( [ viera ] ) , one derives the first scaling law @xmath91 here , @xmath92 denotes the amplitude scale , and @xmath93 abbreviates the first characteristic time scale of the mct - transition scenario :    [ sechs ] @xmath94 the master functions @xmath95 are determined by @xmath87 as solutions of eq .
( [ viera ] ) for @xmath96 , respectively .
they interpolate monotonously between @xmath97 and @xmath98 or @xmath99 .
von schweidler s law is obtained as the long - time limit on the scale @xmath93 in the form : @xmath100 . here
@xmath101 abbreviates the second characteristic scale of the theory : @xmath102    the @xmath0-process obeys for @xmath103 the second scaling - law , called the superposition principle : @xmath104 the control - parameter independent master function @xmath105 exhibits the initial decay @xmath106 .
the parameters @xmath75 , @xmath107 and @xmath87 are determined by @xmath58 from eq .
( [ zwei ] ) for control parameters at the critical point .
the same holds for the function @xmath105 .
the constant @xmath108 is determined by the first taylor coefficient in @xmath21 of the deviations of @xmath58 from its value at the critical point .
formulas for these quantities can be found in ref .
@xcite , where also the original work is cited .
the theory for the leading corrections to the quoted results has been worked out in ref .
@xcite .
the calculations in this paper will be done for the hard - sphere system ( hss ) .
the temperature does not enter the structure , but determines the time scale via the thermal velocity only : @xmath109 .
the relevant control parameter is the packing fraction : @xmath110 , where @xmath111 is the particle diameter .
the structure factor will be calculated within the percus - yevick theory @xcite .
the discretization will be done for @xmath112 wave - vector values with step size @xmath113 . for this model ,
all the mentioned mct quantities have been reported in ref .
in particular it was found : @xmath114 , @xmath115 , @xmath116 the results for the glass - transition of the hss are documented comprehensively in refs .
@xcite , albeit for a brownian microscopic dynamics .
the bifurcation scenario for the model with newtonian dynamics as defined by eqs .
( [ eins ] ) and ( [ zwei ] ) is demonstrated for the wave - vector @xmath117 in ref .
@xcite , where the transient time - scale was determined : @xmath118 . for the presentation of our results in the following figures , the units of length and time will be chosen so that @xmath119 and @xmath120 .
the control parameters @xmath4 shall be cited by the logarithm @xmath121 of the distance parameter @xmath21 : @xmath122 as in the previous work @xcite , the mct equations are solved in the time domain .
the solutions are then laplace - transformed to get @xmath123 .
similarly , the transformed kernel @xmath124 is calculated from @xmath48 in eq .
( [ zwei ] ) .
these results are used to compare the left - hand side of eq .
( [ einsb ] ) with the right - hand side .
thereby a verification of the numerical solutions is obtained .      as a model for a dilute solution of molecules we shall consider a single linear molecule immersed in a simple system .
the position of this molecule is described by the tensor - density fluctuations @xmath125 .
here , @xmath126 denotes the center - of - mass position and @xmath127 abbreviates the axis of the molecule .
the @xmath128 are related to the spherical harmonics by : @xmath129 .
the solute - solvent equilibrium correlations are described by the generalized structure factors @xmath130 , where @xmath131 .
the proper generalization of the density correlators for simple systems are tensor - density correlators for the molecule , @xmath132 .
the mct for these quantities shall be simplified by restricting the correlators to the diagonal elements @xmath133 correlation functions for wave vectors @xmath24 different from @xmath134 can be obtained from the specified ones by elementary transformations @xcite .
the first equation of the mct for the molecule dynamics reads @xcite : @xmath135 here , @xmath136 is the characteristic frequency for the translational motion of a tagged particle .
@xmath137 is the analogue for the rotational dynamics , where @xmath138 plays here and in the following a similar role as @xmath139 for the translational motion .
the frequency @xmath140 denotes the thermal velocity for the rotation .
the relaxation kernels @xmath141 and @xmath142 are approximated along the same lines as indicated above for simple systems .
they are obtained as a functional of the density correlators of the solvent , multiplied by the tensor - density correlators of the solute @xcite .
let us discretize the wave vector to , say , @xmath143 values with equal spacing @xmath144 .
let us also restrict the angular - momentum index by some upper cutoff value @xmath145 .
one obtains the kernels as mode - coupling polynomials @xmath146\\ \label{vierzehn } & = & \sum_{kpl\nu } v_{\alpha qj\mu}(kpl\nu ) ~   \phi(kl\nu , t)~ \phi_{p}(t)\ , , \qquad \alpha={\rm r , t } .\end{aligned}\ ] ] the positive coefficients @xmath147 , @xmath148 are given in ref .
@xcite as specialization of the results in ref .
they are expressed in terms of @xmath29 and @xmath149 for @xmath150 .
anticipating @xmath29 , @xmath151 and @xmath59 as known , eq .
( [ dreizehn ] ) and ( [ vierzehn ] ) are closed equations for the determination of the @xmath152 correlators @xmath153 .
the quantities of main interest for a statistical description of the rotation of the molecule are the reorientational correlators , defined with the legendre polynomials @xmath154 :    [ fuenfzehnneu ] @xmath155 they are the long - wave - length limits of the general correlators ; @xmath156 @xcite .
one gets from eq .
( [ dreizehn ] ) the fraction representation in analogy to eq .
( [ einsb ] ) : @xmath157    here the kernel @xmath158 is the @xmath159 limit of @xmath160 . carrying out the limit in the general formula for @xmath160 @xcite and discretizing the wave - vector integral afterwards , one finds : @xmath161 = \sum_{j~kl\nu}~ v_{kl\nu}^{j}~ \phi(kl\nu , t)\phi_{k}(t).\ ] ] the positive coupling coefficients @xmath162 are listed in ref
after evaluation of @xmath163 and @xmath164 the correlators for the @xmath143 values of @xmath55 , eq .
( [ siebzehn ] ) yields the kernel @xmath165 .
fourier - laplace transformation gives @xmath166 and eq .
( [ sechzehn ] ) provides @xmath167 .
fourier - cosine transformation of the spectrum @xmath168 leads to @xmath169 .
the theory shall be applied for a dumbbell consisting of two equal fused hard spheres of diameter @xmath111 and distance @xmath170 between the centers .
thus , besides the packing fraction @xmath4 , there is the elongation parameter @xmath171 as the second control parameter specifying the structure .
the structure factor @xmath149 and the corresponding pair correlation functions @xmath172 are evaluated within the percus - yevick theory @xcite . figure [ zweiteabb ] exhibits the probability distribution to find a solvent particle in the plane through the symmetry axis of the dumbbell .
the upper panel , calculated for @xmath173 , shows a pronounced quadrupolar pattern extending over several shells .
for the small elongation @xmath174 , the lower panel shows that anisotropy is almost lost from the third shell onwards .
the calculations of the dynamics will be done for such moment of inertia that @xmath175 .
the discretization will be done with @xmath176 wave vectors with spacing @xmath177 .
the cutoff for the angular - momentum index is chosen as @xmath178 for @xmath173 and @xmath179 for @xmath174 .
the equation of motion ( [ dreizehn ] ) is transformed to an integro - differential equation in analogy to eq .
( [ einsa ] ) and then solved by an algorithm similar to that used for the standard mct problem @xcite .
figure [ dritteabb ] demonstrates the transition scenario for the solute correlators for two representative wave vectors @xmath54 and three values of the angular momentum index @xmath1 .
the calculated correlators exhibit a very weak dependence on the helicity index @xmath180 , and therefore only the solutions for @xmath181 are shown .
the wave vector @xmath182 is close to the structure - factor - peak position , and @xmath183 is near the first minimum of @xmath29 .
the correlator for @xmath184 is the probability distribution of the molecule s center - of - mass position , i.e. the analogue of the incoherent - intermediate scattering function for simple liquids : @xmath185 .
results for @xmath1=1 and @xmath1=2 deal with the propagation of the dipole- and quadrupole - density fluctuations , respectively .
the critical - decay curves , i.e. the solutions for @xmath186 , organize the bifurcation pattern .
they deal with the stretched decay towards the plateaus @xmath187 .
if @xmath4 increases above @xmath6 , the long - time limits @xmath188 increase above the plateau because the molecule gets more tightly localized in the frozen solvent .
the @xmath189-versus-@xmath54 curves are bell - shaped , since the molecules are localized with a nearly gaussian probability distribution @xcite . for @xmath64 , the correlators exhibit a long - time decay from the plateau to zero , and this is the @xmath0-process .
the @xmath0-decay time is the larger the smaller the wave vector @xmath54 , while the @xmath0-relaxation stretching increases with increasing @xmath54 .
the @xmath54-dependence of the relaxation features are similar as observed and explained previously for the tagged - particle correlator @xmath190 for simple liquids .
therefore , the following discussion shall be restricted to the ( @xmath191 ) limit , i.e. to the reorientational correlators @xmath169 .
figure [ vierteabb ] exhibits representative decay curves @xmath169 for the liquid state for two separation parameters , and fig .
[ fuenfteabb ] exhibits an extended set of susceptibility spectra @xmath192 .
the plateaus and critical amplitudes shall be denoted by @xmath193 and @xmath194 , respectively .
they have been calculated from the mode - coupling functionals @xcite , and some examples are listed in tab .
[ erstetabelle ] .
these parameters specify the leading - order asymptotic results for the @xmath2-relaxation process as explained in sec .
[ eqnsolv ] for the solvent .
the factorization theorem holds in analogy to eq .
( [ drei ] ) @xcite : @xmath195 the @xmath2-correlator @xmath196 is the same function as explained in connection with eqs .
( [ viera])([sechs ] ) for the solvent .
this implies for the critical correlator the asymptotic law    [ neunzehn ] @xmath197 the nonergodicity parameter of the glass state , @xmath198 , exhibits the @xmath199-singularity @xmath200 and the @xmath0-process initial decay is given by von schweidler s law for @xmath201 and @xmath202 : @xmath203~ , \qquad   \tilde{\tau}_\alpha^j=(f_j^c / h_j)^{1/b } ~ t_\sigma^{\prime } ~ .\ ] ]    let us introduce two ad hoc time scales for the description of the liquid relaxation outside the transient regime .
the center of the @xmath2-relaxation process , @xmath204 , shall be defined as the time , where the correlator has decayed to the plateau : @xmath205 .
the center of the @xmath0-process shall be defined as the time , where the correlator has decayed to 50% of the plateau : @xmath206 .
some values are listed in tab .
[ zweitetabelle ] , and open squares and circles mark these @xmath0- and @xmath2-relaxation times , respectively , in fig .
[ vierteabb ] . the slowing down of the dynamics upon approaching the glass - transition point
is reflected by the increase of the time scales with decreasing distance parameter @xmath207 .
the two step scenario emerges , because the ratio of the scales @xmath208 increases as well .
the @xmath0-decay leads to the @xmath0-peaks of the susceptibility spectrum , which are separated from the microscopic excitation peaks by a susceptibility minimum , as is demonstrated in fig .
[ fuenfteabb ] .
there , the @xmath0-peak - maximum positions , @xmath209 , and the minimum positions , @xmath210 , decrease for @xmath211 . the open squares and circles in fig .
[ fuenfteabb ] demonstrate , that @xmath212 and @xmath213 as @xmath214 .
the two - step scenario implies that the ratio @xmath215 also decreases upon approaching the glass - transition point .
thus , for @xmath211 , the @xmath0-peak gets more and more separated from the rest of the spectrum . in this limit , the plateau height is the relative area under the @xmath216-versus-@xmath217 curve @xcite : @xmath218    figure [ vierteabb ] demonstrates that for @xmath219 the dynamics deals with oscillatory motion , i.e. with rotations and librations which are influenced by steric hindrance affects .
if these effects would lead to some fast decay towards the correlator s long - time limit , one would find a white - noise low - frequency fluctuation spectrum : @xmath220 .
equivalently , one would obtain a regular low - frequency susceptibility spectrum varying linearly with @xmath40 , @xmath221 , as is indicated schematically by the straight dashed - dotted line in the upper left panel of fig .
[ fuenfteabb ] .
a linear susceptibility spectrum is obtained for the glass - spectra for @xmath222 , since the correlators approach the limit @xmath223 exponentially for @xmath224 .
this is shown by the ( @xmath225 ) spectra in fig .
[ fuenfteabb ] .
such regular spectra are also found for the low - frequency wings of the @xmath0-peaks , since the liquid correlators approach zero exponentially for @xmath226 . at the bifurcation point
, however , the critical decay leads to a power law spectrum which , according to eq .
( [ neunzehna ] ) , reads @xmath227 for @xmath228 , the correlators follow the critical decay decay if @xmath229 is small .
therefore the spectra are approaching the asymptotic @xmath230-law for @xmath231 , as is demonstrated for the ( @xmath232 ) results in fig .
[ fuenfteabb ] .
the stretching of the first relaxation step leads to the strong enhancement of the intensity of the spectral minimum @xmath233 relative to any possible estimation of a white - noise - background spectrum .
this enhancement also is exhibited by the experimental data reproduced in fig .
[ ersteabb ] .
let us consider the probability density @xmath234 for the molecule s axis @xmath235 to have the projection @xmath236 on its initial direction @xmath127 : @xmath237 .
since @xmath238 , one gets    [ zweiundzwanzig ] @xmath239 thus , knowledge of the set of @xmath169 , @xmath240 , is equivalent to knowing @xmath241 .
if the summation over @xmath1 is understood with the cutoff @xmath145 , eq .
( [ zweiundzwanziga ] ) describes the evolution of the distribution with the initial value @xmath242    figure [ sechsteabb ] exhibits results for the small distance parameter @xmath243 corresponding to @xmath244 .
the dotted lines exhibit @xmath245 , calculated with @xmath246 for @xmath173 , and @xmath247 for @xmath174 , respectively . within the dynamical window , where the leading - order result for the @xmath2-relaxation , eq .
( [ achtzehn ] ) , applies , one gets    [ dreiundzwanzig ] @xmath248
@xmath249    thus the distribution relaxes towards the distribution @xmath250 , which is frozen for @xmath186 .
the relaxation does not exhibit any correlation between changes in time described by @xmath84 , and variations with angle described by @xmath251 .
this is the scenario expected for relaxation due to dephasing in the random distribution of sizes and shapes of the cages producing steric hindrance for the rotations . for @xmath173 , the @xmath2 regime extends from @xmath252 to about @xmath253 as shown in fig .
[ vierteabb ] , and the upper panel of fig .
[ sechsteabb ] exhibits the described phenomena for @xmath254 and @xmath255 .
the @xmath2-relaxation window is somewhat smaller for @xmath174 , as will be discussed in quantitative detail below in connection with fig .
[ siebteabb ] . the dephasing relaxation for this case is demonstrated in fig . [ sechsteabb ] for @xmath254 and @xmath256 .
the beginning of the @xmath0-relaxation process follows von schweidler s law , eq .
( [ neunzehnc ] ) .
it is identical with the end of the @xmath2-process , and thus it is described within the scenario based on eqs .
( [ dreiundzwanzig ] ) .
the most drastic difference between large- and small - elongation relaxation shows up for the @xmath0-process outside the von - schweidler regime . for @xmath173 ,
the probability decreases monotonically if the angle @xmath257 of the axis increases from its initial value @xmath258 to @xmath259 .
this is shown in the upper panel of fig .
[ sechsteabb ] for @xmath260 . as time increases , the probability for @xmath261 decreases , while it increases for @xmath262 .
thus , the relaxation towards the equilibrium distribution @xmath263 is similar to what one would expect for diffusion on a sphere . for @xmath174 , the correlators for odd @xmath1 decay faster than the corresponding correlators with the even index @xmath264
this is demonstrated in fig .
[ vierteabb ] and by the numbers @xmath265 in tab .
[ zweitetabelle ] .
therefore the @xmath0-process consists of an intermediate time step leading to a probability distribution which is nearly symmetric with respect to the equator @xmath266 . only at later times , the symmetric distribution relaxes to the equilibrium one .
figure [ sechsteabb ] shows that , already for the rather short time @xmath254 , @xmath241 exhibits a minimum . for @xmath267
, there is an overshooting effect of the probability for @xmath268 : @xmath269 ; and this effect increases if the time increases to @xmath270 .
thus , the relaxation pattern is that expected for a random process of large - angle flips of the molecule s axis .
the equations for the nonergodicity parameters @xcite imply that the @xmath223 increase towards unity if the coupling coefficients in eq .
( [ vierzehn ] ) are increased towards infinity . for this strong - coupling limit ,
one derives from eq .
( [ zweiundzwanziga ] ) that @xmath271 .
because of continuity , for strong steric hindrance and for @xmath272 , @xmath241 is a narrowly - peaked distribution centered around @xmath261 .
thus , one expects the expansion coefficients @xmath223 for not too large values of @xmath1 to form a smoothly decreasing sequence of @xmath1 : @xmath273 , @xmath274 .
table [ erstetabelle ] demonstrates this result quantitatively for @xmath173 and @xmath186 :    [ vierundzwanzig ] @xmath275 in particular , the ratio @xmath276 of the relative strengths of the @xmath0-peaks for the dipole relaxation , @xmath277 , and for the quadrupole relaxation , @xmath278 , is larger than unity .
one can not conclude quantitatively from @xmath279 the ratio @xmath280 of the @xmath0-peak heights , since the shapes of the spectra depend on @xmath1
. however , fig .
[ fuenfteabb ] demonstrates that the two ratios are close to each other .
one can also characterize the @xmath0-peak height relative to the microscopic - peak height , @xmath281 , or relative to the minimum intensity , @xmath282 . from fig .
[ fuenfteabb ] , one infers @xmath283 , and @xmath284 , i.e. the ( @xmath1=1)-versus-(@xmath1=2 ) enhancement effect appears even more pronounced .    according to eq .
( [ neunzehnb ] ) , the nonergodicity parameters increase with increasing @xmath285 . on the other hand , @xmath286 .
therefore , @xmath194 must decrease if @xmath193 increases , so that the strongly coupled parameters @xmath223 leave the asymptotic regime for eq .
( [ neunzehnb ] ) for similar magnitudes of @xmath89 .
table [ erstetabelle ] quantifies this result for @xmath173 .
in particular @xmath287    the reasoning assumes @xmath223 to be large , and thus it can not be applied for too large @xmath1 .
there is some @xmath288 , so that @xmath194 decreases with increasing @xmath1 for @xmath289 . within the frequency window , where the leading - order asymptotic law for the critical decay is valid , eq .
( [ einundzwanzig ] ) , one derives an enhancement of the ( @xmath1=2 ) spectrum relative to the ( @xmath1=1 ) spectrum , since : @xmath290 the dotted lines @xmath291 and @xmath292 in the upper panels of fig
. [ fuenfteabb ] demonstrate this result .    for a strongly near-(@xmath293=1)-peaked probability distribution @xmath241 , one can approximately replace averages of functions of @xmath293 by the functions of the average @xmath294 .
thus , lebon et al . concluded : @xmath295 @xcite .
specializing to @xmath186 , one quantifies the sequence of @xmath193 in terms of its first value @xmath277 :    [ sechsundzwnzig ] @xmath296 substituting into eq .
( [ neunzehnb ] ) and specializing to @xmath297 , one can also quantify the sequence of @xmath194 by the first term @xmath298 : @xmath299    from tab .
[ erstetabelle ] one infers , that for @xmath173 the error of eq .
( [ sechsundzwnziga ] ) for @xmath300 is as small as @xmath301 , and eq .
( [ sechsundzwanzigb ] ) is obeyed for @xmath300 within @xmath302 .
the strong nonlinear couplings of the structural - relaxation modes require that all correlators enter the first relaxation step , the second relaxation step , and the equilibrium state nearly at the same respective time .
this is demonstrated in fig .
[ dritteabb ] for @xmath173 .
the most striking manifestation of the coupling effect occurs at the center of the @xmath2-relaxation window for @xmath303 . in this case , the factorization theorem , eq .
( [ achtzehn ] ) , is valid .
all correlators cross their plateau at the same time , say @xmath304 , where @xmath304 is the zero of the @xmath2 correlator @xmath84 . because of the scaling law , eq .
( 5 ) , one gets the result @xmath305 , i.e. @xmath306 here , @xmath307 is the scale from eq .
( [ sechsa ] ) , and @xmath308 is the zero of the master function : @xmath309 . for the hss it reads @xmath310 @xcite .
the open and full circles in fig .
[ vierteabb ] show that the asymptotic eq .
( [ siebenundzwanzig ] ) is obeyed very well for @xmath173 .
since the @xmath0-processes of @xmath311 and @xmath312 start at the same time @xmath313 and reach zero nearly at the same time , one expects from @xmath314 that the decay time for @xmath315 is larger than that for @xmath316 : @xmath317 furthermore , the @xmath318-versus-@xmath319 plot is somewhat steeper than the corresponding graph for @xmath1=2 .
this means that the stretching is larger for the ( @xmath1=2 ) @xmath0 process than for the ( @xmath1=1 ) @xmath0 process .
if one interpolates the decay functions by a kohlrausch law , @xmath320 $ ] , the stretching exponent for @xmath1=1 is larger than that for @xmath1=2 : @xmath321 stretching can also be quantified by the width @xmath322 at half height of the @xmath0-peak of the susceptibility spectrum . for @xmath173 our model yields for @xmath323
@xmath324 decades .
the kohlrausch processes leading to the same @xmath322 require stretching exponents @xmath325 .
the derivation of the inequality for the time scales can be put on a quantitative level by combining eq .
( [ neunzehnc ] ) with the two inequalities in eqs .
( [ vierundzwanzig ] ) .
one gets in analogy to eq .
( [ achtundzwanzig ] ) : @xmath326 .
the @xmath0-relaxation law for the @xmath169 holds in analogy to eq .
( [ sieben ] ) : @xmath327 .
if the shape function @xmath328 would be independent of @xmath1 , the ratio @xmath329 would be equal to the ratio @xmath330^{1/b}$ ] .
but the latter is about @xmath331 times larger than @xmath329 .
there are two universal phenomena which are relevant for a discussion of the dynamics for weak steric hindrance .
the first one concerns the limit @xmath332 of the center - of - mass correlator @xmath333 , which is identical to the tagged - particle - density correlator of the larger of the two spheres forming the dumbbell .
if the radius of this sphere , say @xmath334 , is of the same order or larger than the radius @xmath111 of the solvent spheres , the steric hindrance is very effective . in this case ,
@xmath190 exhibits the canonical bifurcation scenario if @xmath4 crosses @xmath6 , as was discussed comprehensively in ref .
this implies that for @xmath335 the ( @xmath184 ) correlators exhibit only a smooth @xmath171-dependence for @xmath171 decreasing to zero .
a side remark shall be added to this conclusion .
if the ratio of the diameters @xmath336 of a sphere moving in a glass of hard spheres decreases towards zero , there occurs a percolation transition at some critical value @xmath337 .
this is a type - a transition , i.e. a bifurcation where the lamb - mssbauer factor decreases continuously to zero for @xmath338 approaching zero from above @xcite . because of continuity , it is obvious that for a dumbbell built of sufficiently small spheres , @xmath339 , there will be a type - a transition if the elongation @xmath171 decreases to some critical value @xmath340 .
if @xmath171 crosses @xmath341 , the dynamics changes from one dealing with molecules localized in the hard - sphere glass to one dealing with delocalized molecular motion .
this small-@xmath171 phenomenon for small molecules is not considered in this paper .
the second universal phenomenon deals with a type - a transition resulting from the fact , that for top - down symmetrical molecules the mct equations of motion of the even-@xmath1 correlators decouple from the odd-@xmath1 ones @xcite .
the even-@xmath1 correlators couple to the function @xmath190 , and thus the conventional transition scenario of this correlator enforces the same for all other correlators with even @xmath1 .
however , such coupling does not exist for odd @xmath1 . for large @xmath171 , this results in no considerable effect .
but all coupling coefficients in the equations of motion approach zero for odd @xmath1 if @xmath171 tends to zero .
consequently , for all @xmath342 there is some critical elongation @xmath343 for a type - a transition .
for the studied model @xmath344 @xcite .
choosing @xmath171 sufficiently close to @xmath343 , it can happen that for odd @xmath1 : @xmath345 or even @xmath346 @xcite .
the transition at @xmath343 shall not be studied in this paper . for the demonstration of the small-(@xmath347 ) phenomena ,
the value @xmath174 has been chosen so large that the canonical sequence for the plateau values in the ( @xmath191 ) limit , eq .
( [ vierundzwanziga ] ) , is not violated , as is quantified in tab .
[ erstetabelle ] .
but it is chosen so small , that the precursor effects of the type - a transition seriously influence the results for the dumbbell dynamics .
thereby , the results are also representative for such cases , where the type - a transition singularity is avoided @xcite due to a weak breaking of the top - down symmetry of the solvent - solute interaction .
the even-@xmath1 correlators show the conventional behavior .
therefore , the discussion of their trends with decreasing @xmath171 for fixed @xmath1 can be held analogously to that given in sec .
iiib for the trends with increasing @xmath1 for fixed large @xmath171 .
thus one understands that the ( @xmath1=2 ) @xmath0-process for @xmath174 is weaker , faster , and more stretched than that for @xmath173 , as it is demonstrated in figs .
[ vierteabb ] and [ fuenfteabb ] .
the halfwidth of the @xmath0-peak for @xmath348 is @xmath349 decades as for a kohlrausch process with exponent @xmath350 .
notice in particular from fig .
[ vierteabb ] , that the @xmath2-relaxation scale @xmath351 for the @xmath244-result is close to the @xmath171- and @xmath1-independent number @xmath313 from eq .
( [ siebenundzwanzig ] ) . for @xmath352 ,
the asymptotic formula is obeyed reasonably , but the preasymptotic corrections are larger for @xmath174 than for @xmath173 .    the most obvious precursor of the type - a transition is the suppression of the plateau values @xmath223 for odd @xmath1 .
this leads to a violation of the rule @xmath353 , as is quantified in tab .
[ erstetabelle ] . the general qualitative reasoning from sec .
iiib explains , that the suppression of @xmath354 is connected with an enhancement of @xmath298 : @xmath355 .
the amplitude @xmath298 is given by the resolvent of the so - called stability matrix , and at a type - a transition the resolvent exhibits a pole @xcite .
hence @xmath356 , and the regular trend , eq .
( [ vierundzwanzigb ] ) , is reversed : @xmath357 for our example one infers from tab .
[ erstetabelle ] that @xmath358 . according to eq .
( [ fuenfundzwanzig ] ) , the critical spectrum for the dipole relaxation is considerably larger than that for the quadrupole relaxation , as is demonstrated by the dotted lines in the lower two panels of fig .
[ fuenfteabb ] .    combining eq .
( [ einunddreissig ] ) with von schweidler s law , eq .
( [ neunzehnc ] ) , one concludes that the @xmath311-versus-@xmath319 curve crosses its plateau @xmath277 much steeper than the @xmath312-versus-@xmath319 curve .
this is illustrated in the lower panel of fig .
[ vierteabb ] .
hence the @xmath0-relaxation of the ( @xmath1=1 ) response is faster than the one of the ( @xmath1=2 ) response : @xmath359 again the order for large @xmath171 , eq .
( [ achtundzwanzig ] ) , is reversed . from tab .
[ zweitetabelle ] , one infers for @xmath244 : @xmath360 . accordingly , the @xmath0-peak positions for the ( @xmath244 ) spectra for @xmath1=1 and @xmath1=2 in the lower panels of fig . [ fuenfteabb ] differ by about one order of magnitude . for
the ratio of the von schweidler scales in eq .
( [ neunzehnc ] ) , one gets @xmath361^{1/b}\to 0 $ ] for @xmath362 , and this identifies the smallness of the ratio @xmath363 as a precursor of the type - a transition . the preceding discussion is valid more generally and explains that all the odd-@xmath1 correlators decay faster than the nearby even-@xmath1 ones . as a result ,
the probability distribution @xmath241 approaches first an even function of @xmath293 , before the equilibrium value @xmath364 is obtained , as is demonstrated in the lower panel of fig .
[ sechsteabb ] .      the factorization theorem for the @xmath2-relaxation , eq .
( [ achtzehn ] ) , means that the rescaled correlators @xmath365 are given independently from @xmath1 by the @xmath2-correlator @xmath84 of the solvent .
the latter obeys the scaling law , specified by eqs .
( [ viera])-([sechsa ] ) . for fixed rescaled time @xmath366 ,
the cited formulas deal with the results correctly up to order @xmath367 @xcite .
the leading corrections are of order @xmath229 , and they explain the range of validity of the leading results for separations @xmath21 @xcite . figure [ siebteabb ] demonstrates these statements . on a @xmath368-accuracy level
the leading - order results describe @xmath369 of the decay of the correlators around the plateau for @xmath173 , @xmath370 .
these decay intervals are indicated in fig .
[ vierteabb ] by vertical lines . for @xmath371 ,
the corresponding dynamical window extends from about @xmath252 to about @xmath372 , while it extends from about @xmath373 to about @xmath374 for @xmath375 .
this discussion requires a reservation : the corrections to the scaling results can lead to such a violation of eq .
( [ siebenundzwanzig ] ) , which appears as an offset of the plateau @xcite .
this offset can be noted in the lower panel of fig .
[ siebteabb ] for the odd-@xmath1 results . the good description of the @xmath2-decay of the ( @xmath174 ) results for odd @xmath1 holds only after a correction of the offset . for @xmath376 , the correction effects cause the @xmath377
for @xmath173 to differ from @xmath84 ; one infers from fig .
[ siebteabb ] that the @xmath377 increase with increasing @xmath1 .
the general results for the theory of the corrections imply , that then also the @xmath377 increase with @xmath1 for @xmath378 @xcite .
the @xmath377-versus-@xmath379 curves do not intersect for @xmath313 but they touch each other as is demonstrated in the upper panel of fig .
[ siebteabb ] .
corresponding results also hold for @xmath174 after the mentioned offset is eliminated
.    equations ( [ fuenf ] ) and ( [ achtzehn ] ) lead to the scaling law for the susceptibility spectra : @xmath380 .
the master spectra @xmath381 are given by the fourier - cosine transform @xmath382 of the master functions @xmath95 .
the master spectrum for the glass state describes the crossover from a regular spectrum for small rescaled frequencies , @xmath383 , to the critical spectrum at large rescaled frequencies , @xmath384 .
it deals with the knee exhibited by the spectra for @xmath225 and @xmath385 in fig .
[ fuenfteabb ] .
the master spectrum for the liquid describes the crossover from the von - schweidler high - frequency tail of the @xmath0-peak , @xmath386 , to the critical decay for large rescaled frequencies , @xmath387 .
the results describe in the small-@xmath89 limit the @xmath2-relaxation minimum as it can be seen in fig .
[ fuenfteabb ] for the @xmath244 and @xmath232 results .
in particular , the factorization theorem explains why the spectral minima @xmath388 are located at the same position independently of @xmath1 and @xmath171 .
the leading - order formulas imply @xmath389 , where @xmath390 denotes the minimum of the master spectrum @xmath391 . for the hard - sphere system ,
one gets @xmath392 @xcite .    obviously , the @xmath2-relaxation scaling laws can describe the susceptibility minimum only for such small distance parameters , for which @xmath388 is located in that frequency window where the ( @xmath393 ) spectrum exhibits the asymptotic @xmath394-law , eq .
( [ einundzwanzig ] ) .
figure [ fuenfteabb ] shows that for the model under study this window is restricted to @xmath395 .
this means , that @xmath388 has to be located about three decades below the peak of the microscopic susceptibility spectrum . for @xmath396 ,
the critical spectrum is modified by crossover effects to the transient dynamics .
the susceptibility minimum with @xmath397 is due to the crossover of the @xmath0-peak tail to the microscopic excitation spectrum ; it can not be discussed by the universal asymptotic laws for the mct bifurcation .
one concludes from fig .
[ fuenfteabb ] , that @xmath398 needs to be satisfied in order to apply the @xmath2-scaling laws for the model under study .
the @xmath0-relaxation scaling law reads for the reorientational correlators in analogy to eq .
( [ sieben ] ) : @xmath399 the @xmath21-independent master function @xmath400 obeys as initial condition the von schweidler law : @xmath401 .
the superposition principle for the susceptibility spectra reads @xmath402 with @xmath403 denoting the rescaled frequency .
the @xmath21-independent master spectrum is given by the fourier - cosine transform of the master correlators @xmath404 .
consequently , the above defined @xmath0-relaxation time scales @xmath405 and susceptibility - maximum positions @xmath209 read @xmath406 where @xmath407 is defined by @xmath408 and @xmath409 denotes the peak frequency of @xmath410 .
the scaling law implies that a representation of @xmath169 as a function of the rescaled time @xmath411 should superimpose correlators for different distance parameters @xmath21 on the common curve @xmath412 .
asymptotic validity means that the @xmath413 interval , where the scaling law is obeyed , expands to arbitrary size for @xmath103 .
a corresponding statement holds for the representation of the susceptibility peaks as functions of the rescaled frequency .
the corrections to the leading - order asymptotic laws are the larger , the larger the critical amplitude @xmath194 is @xcite .
figures [ fuenfteabb ] and [ achteabb ] demonstrate , that the described scenario for the evolution of the @xmath0-process is valid for @xmath173 , and also for @xmath174 in the case @xmath1=2 . for strong steric hindrance , the @xmath0-scaling law works for larger values of @xmath414 , than the @xmath2-scaling law .
this is so , because the leading corrections to the @xmath0-scaling law are of relative size @xmath415 , while they are of relative size @xmath416 for the @xmath2-scaling law @xcite .
figure [ achteabb ] demonstrates a drastic ( @xmath1=1)-versus-(@xmath1=2 ) effect of the @xmath0-scaling for @xmath174 .
the dipole correlators do not exhibit the superposition principle for @xmath417 , nor do the correlators for the other odd values of @xmath1 .
for @xmath418 the plateau emerges only for extremely small values of the distance parameter @xmath419 .
the @xmath0-peak heights of the dipole spectra decrease with decreasing @xmath207 in fig .
[ fuenfteabb ] in contradiction to the scaling - law prediction .
this anomaly is caused by the large size of the critical amplitude @xmath298 , which was explained in connection with eq .
( [ einunddreissig ] ) .
more precisely , it is caused by the large percentage of the decay of @xmath311 described by the @xmath2-scaling law as is indicated by the vertical lines in fig .
[ vierteabb ] . to formulate this observation quantitatively ,
let us remember that the decay of the correlator near the plateau is described by eqs .
( [ fuenf ] ) , ( [ sechs ] ) and ( [ achtzehn ] ) : @xmath420 . the master function @xmath421 for small positive values and all negative ones is well approximated by :
@xmath422 . here
, @xmath423 is determined by the exponent parameter @xmath87 and for our solvent model reads @xmath424 @xcite .
thus one gets for @xmath425 within the window for the validity of the @xmath2-relaxation scaling law :    [ vierunddreissig ] @xmath426 \}~ .\ ] ] the leading corrections to this formula can explain the possible offset of @xmath193 or , equivalently , of the scales @xmath307 @xcite , which was noticed above in connection with fig .
[ siebteabb ] for @xmath174 .
equation ( [ vierunddreissiga ] ) can be rewritten as @xmath427 + h_j |\sigma| b_1/\tilde{t}^b$ ] . here , the bracket is the @xmath0-scaling - law description of the initial part of the @xmath0-process , and the term proportional to @xmath423 is the leading correction .
the correction term to the @xmath0-scaling law deals with that part of the @xmath2-process below the plateau , which is not given by the von schweidler s large-@xmath428 asymptote .
therefore , one can write the for the @xmath0-process for not too large values of rescaled time @xmath429 : @xmath430    the analogue of this formula was shown in ref .
@xcite to describe the evolution of the @xmath0-process of the density correlators of the hss perfectly for @xmath431 .
it was also shown , that the corresponding spectrum describes the susceptibility peak to increase above the scaling - law constant @xmath432 if the separation @xmath207 increases from @xmath433 to @xmath434 .    in fig .
[ neunteabb ] , the evolution of the ( @xmath1=1 ) @xmath0-process for small steric hindrance is reexamined . instead of rescaling the time with the ad - hoc scale @xmath435
, the theoretically motivated scale @xmath436 is chosen .
one recognizes , that the found scenario does not exhibit any qualitative peculiarity anymore , compared to what is presented in fig .
[ achteabb ] for @xmath174 and @xmath1=2 .
the ( @xmath1=1)-versus-(@xmath1=2 ) anomaly is identified as an anomaly of the size of the corrections only . in the case of the small elongation
, the distance parameter @xmath207 has to be taken almost two orders of magnitude smaller for @xmath1=1 in order to render the corrections to the @xmath0 scaling as small as found for @xmath1=2 . for @xmath437 , even including the leading corrections to von schweidler s law , one can explain the relaxation from the plateau only up to some offset in the time scale .
this is demonstrated in fig .
[ neunteabb ] by the dotted lines for @xmath438 .
two remarks concerning tests of the second scaling law shall be added . the definition of the time scale @xmath265 used in fig .
[ achteabb ] was arbitrary .
let us consider more general definitions to be parameterized by a number @xmath439 and denoted as @xmath440 .
the subscripts @xmath0 and @xmath1 shall be dropped for the sake of simplicity , and the definition shall be : @xmath441 .
if the scaling law is valid , one finds in analogy to eq .
( [ dreiunddreissig ] ) : @xmath442 . here
, @xmath443 is defined by @xmath444 . in this case , the choice of @xmath55 is irrelevant , since the ratio of two different scales is @xmath21-independent : @xmath445 . however ,
if preasymptotic corrections are present , the scales are not equivalent .
the range of validity of the superposition principle expands from large to small rescaled times .
this follows from eq .
( [ vierunddreissigc ] ) and is demonstrated in fig .
[ neunteabb ] .
one gets for @xmath446 : @xmath447 for a detection of the superposition principle for an as large as possible value of @xmath207 , one should therefore choose an as large as possible value of @xmath55 for the rescaling procedure . thereby , the artificial crossing point of the rescaled curves at @xmath448 is suppressed as much as possible .
otherwise , one introduces a time scale @xmath440 for the characterization of a decay process which can not be characterized by a single scale .
the outcome of this ill - defined procedure is demonstrated in the lower left panel in fig .
[ achteabb ] . in this case
, @xmath449 is a parameter extracted from the correlator which , according to fig . [ siebteabb ] , is adequately to be specified by the two scales @xmath450 and @xmath307 of the @xmath2-relaxation scaling law . the dashed line in the inset of fig .
[ neunteabb ] demonstrates explicitly , that the scale @xmath435 does not exhibit the asymptotic behavior for @xmath451 .
the asymptotic law @xmath452 is followed only for @xmath453 .
the second remark concerns the determination of the exponent @xmath454 entering the power - law behavior for the @xmath0-relaxation time scale , as specified by eqs .
( [ sechsb ] ) and ( [ dreiunddreissig ] ) .
these results are based on the validity of the scaling law @xcite .
therefore , one can not appeal to mct if one fits power laws for scaling times for cases where the scaling law is violated .
the dashed line in the inset in fig .
[ neunteabb ] demonstrates , that the scale @xmath435 for @xmath455 can be fitted well by a power law for a two - decade variation of the distance parameter @xmath207 .
the identified effective exponent @xmath456 describes the variation of @xmath457 over three orders of magnitude ; but nevertheless @xmath458 has no well defined meaning for the discussion of our model .
solving the mct equations of motion for the dynamics of a hard - sphere dumbbell moving in a hard - sphere liquid , first - principle results have been obtained for the evolution of the glassy dynamics of the reorientational degrees of freedom of a molecule .
it was found that one has to distinguish between two scenarios , namely between strong steric hindrance as found for large elongations @xmath171 of the dumbbell , and weak steric hindrance as found for small elongations .    for strong steric hindrance , the mode - coupling coefficients for the reorientational degrees of freedom in eq .
( [ vierzehn ] ) are of the same order as the ones entering eq .
( [ zwei ] ) for the description of the translational degrees of freedom of the solvent .
the dependence of the various parameters on the angular - momentum index @xmath1 is similar to the dependence on the wave vector @xmath54 .
one has to view @xmath138 as the analogue of @xmath139 .
while the @xmath54-dependence reflects the decomposition of the direct solute - solvent correlations in plane waves the @xmath1-dependence reflects the decomposition in spherical harmonics .
hence one finds , that all reported results on the @xmath1-dependence of the reorientational correlators @xmath169 are similar  and can be explained in a similar manner  as known from the previous work on the tagged - particle - density correlators @xmath190 in simple liquids @xcite . in particular , it was shown that with increasing @xmath1 the @xmath0-peak - strength parameters @xmath193 , eq .
( [ vierundzwanziga ] ) , the @xmath0-relaxation - time scales @xmath405 , eq .
( [ achtundzwanzig ] ) , and the stretching exponents @xmath459 , eq .
( [ neunundzwanzig ] ) , decrease .
these findings reproduce qualitatively the three general differences between dielectric - loss and depolarized - light - scattering spectra , which were discussed in sec .
i in connection with fig .
[ ersteabb ] . because of eq .
( [ einsa ] ) , the relaxation of the correlator @xmath460 follows that of the kernel @xmath461 .
therefore , the @xmath0-relaxation time scale of the ( @xmath159 ) density fluctuations , say @xmath462 , is larger than the corresponding scale of the longitudinal elastic modulus @xmath463 , say @xmath464 . for strong steric hindrance ,
the decay of the cage is the prerequisite for the reorientation of the molecule , and therefore @xmath465 .
thus one expects the fourth general feature of the @xmath0 relaxation listed in the introduction : @xmath466 . for our model one
gets for @xmath375 : @xmath467 , @xmath468 , @xmath469 .
the ratio @xmath470 is of the same order as cited in sec .
i for pc and salol .
the mode - coupling coefficients in eq .
( [ vierzehn ] ) decrease to zero if @xmath1 tends to infinity .
thus the solutions for large @xmath1 are sums of many small terms , which are not strongly correlated .
each term exhibits the short - time von schweidler - law behavior for the @xmath0-relaxation : @xmath471 .
therefore , one expects for @xmath169 the characteristic function of the stable lvy distribution , @xmath472 $ ] @xcite . for the density correlators of the solvent ,
fuchs has worked out the limit behavior for the @xmath0-relaxation master function for @xmath473 and showed how the kohlrausch law with @xmath474 arises @xcite .
we suspect that a similar derivation can be done for the reorientational correlators .
therefore , we conjecture that the sequence of kohlrausch exponents @xmath475 converges towards the von schweidler exponent @xmath80 . molecular - dynamics - simulation data for a model of water have been interpreted consistently within the standard mct scenario @xcite . in particular ,
the @xmath169 exhibit the conventional behavior @xcite .
one concludes that water exhibits strong - steric - hindrance effects .
therefore it is reinsuring that the sequence of the first five kohlrausch exponents @xmath459 decreases with increasing @xmath1 monotonously towards the von schweidler exponent @xcite .
a further general result , namely the increase in the initial part of the series of critical amplitudes , eq . ( [ vierundzwanzigb ] ) , is also found in the simulation data for @xmath476 @xcite .
figure [ ersteabb ] exhibits as full lines the ( @xmath1=1 ) and ( @xmath1=2 ) spectra calculated for @xmath173 .
the lines for @xmath477 and @xmath352 are the ones discussed in fig .
[ fuenfteabb ] , and the other two refer to @xmath478 and @xmath479 , respectively . in order to transfer the mct results , which are calculated with ad hoc units specified in sec .
ii , to the units used by the experimentalists , one has to introduce three scales .
the first and second scale transfer the calculated dimensionless normalized spectra @xmath480 for @xmath1=1 and @xmath1=2 to the units used by the experimentalists for their dielectric - loss and depolarized - light - scattering spectra , respectively .
the third scale shifts our frequency scale to the ghz scale . in the double - logarithmic representation
, the first two scales define an overall vertical shift of the diagrams in fig .
[ fuenfteabb ] , while the third scale defines a horizontal shift of the figures .
intending to compare data for pc for different temperatures with the mct results for different packing fractions @xmath4 , one gets a mapping of the @xmath3-scale onto the @xmath4-scale via eq .
( [ elf ] ) .
the result is shown as an inset in fig .
[ ersteabb ] .
the inset also includes the point with coordinates of the critical packing fraction of the solvent @xmath6 and the critical temperature @xmath481 .
this value for @xmath5 was determined for pc by analyzing neutron scattering data @xcite , and has recently been corroborated in an mct analysis of various other pc experiments @xcite .
our results in fig .
[ ersteabb ] describe the evolution of the two types of pc spectra semi - quantitatively . in particular , the extrapolation of the @xmath3-@xmath4-relation yields a reasonable estimation of the critical temperature for that system , which is demonstrated through the dashed line in the inset .
there is no obvious reason why the studied dilute solution of hard symmetric dumbbells in a hard - sphere solvent should produce spectra , which are similar to the data for pc .
we consider the found similarities to a large extend as accidental .
the theoretical curves are added in fig .
[ ersteabb ] with the mere intention to justify the conclusion : the model studied in this paper and our choice of parameters are relevant for achieving an understanding of experiments in glass - forming molecular liquids .    in order to further corroborate the preceding conclusion ,
let us consider fig .
[ zehnteabb ] .
the data points exhibit a susceptibility spectrum of pc measured by incoherent - neutron - scattering spectroscopy for the wave vector @xmath482 @xcite .
a remarkable feature of the @xmath0-peak spectrum compared to the spectra shown in fig . [ ersteabb ] is that it is less pronounced relative to the spectrum of the microscopic excitation band and that is is more stretched .
the two dashed lines in fig .
[ zehnteabb ] exhibit the spectra for the center - of - mass correlator @xmath483 for @xmath484 in order to emphasize , that this leading approximation for the scattering function can not easily explain the experimental findings .
the shown @xmath0-peaks of @xmath485 have a half width of @xmath486 decades as produced by a kohlrausch process with exponent @xmath487 .
the scattering function @xmath488 is a sum over the contributions of the molecule s constituents and hence it is a superposition of the density correlators for all angular momentum indices @xmath1 . for the symmetric dumbbell one
gets up to some normalization constant @xcite @xmath489 where @xmath490 are the spherical bessel functions .
the full lines in fig .
10 exhibit the spectra for @xmath488 for @xmath484 .
the @xmath0-peaks have a halfwidth of @xmath491 decades as produced by a kohlrausch law with stretching exponent @xmath492 .
the frequency was rescaled as explained in connection with fig .
[ ersteabb ] and the scale for the spectra was adjusted to meet the one of the data .
comparison of the full line with the dashed one for @xmath352 shows the features distinguishing the @xmath0-processes of @xmath493 from that of @xmath494 .
the result calculated for @xmath477 shows that the finding for our model semiquantitatively accounts for the @xmath0-peak data .    some side remarks considering the comparisons in figs .
[ ersteabb ] and [ zehnteabb ] might be useful .
a schematic - model analysis of the pc data gave the exponent parameter @xmath495 @xcite , in good agreement with the values found from analyses of the susceptibility minima with the @xmath2-relaxation scaling laws @xcite .
the value is close to the result @xmath496 for the hard - sphere system , eq .
( [ neun ] ) .
this accident ensures that the master function for the susceptibility minimum and the values of all anomalous exponents of pc agree within the experimental uncertainties with the corresponding quantities of the model studied in this paper ; and this is a prerequisite of a successful fit .
accidently , the ratio of the @xmath0-relaxation times @xmath497 noted in tab .
[ zweitetabelle ] for @xmath498 is only a bit larger than the ratio of the @xmath0-peak - maximum positions of pc , @xmath499 , exhibited in fig . 1 ; and this is another request for a reasonable fit .
since the ratio decreases with decreasing @xmath171 , some @xmath500 could be chosen to reproduce the specified ( @xmath1=1)-versus-(@xmath1=2 ) effect quantitatively .
nevertheless , it is remarkable that the fit in fig .
[ ersteabb ] reproduces the ratios of @xmath0-peak maximum intensity to @xmath2-minimum intensity @xmath501 reasonably well for both values of @xmath1 .
neither it is trivial , that the model reproduces reasonably the ( @xmath1=1)-versus-(@xmath1=2 ) effect for the stretching .
a new liquid - glass - transition scenario is predicted which is referred to as the regime of weak steric hindrance for reorientational motion .
it is characterized by ( @xmath1=1)-versus-(@xmath1=2 ) effects , more generally by odd-@xmath1-versus - even-@xmath1 effects , which are qualitatively different from the results described above as strong - steric - hindrance results .
the new scenario occurs , if precursor phenomena of a type - a - transition between two nonergodic states strongly influence the asymptotic results for the conventional mct bifurcation .
the scenario appears if the particle interactions deviate not too strongly from spherical symmetry , e.g. , if a linear molecule exhibits only small deviations from a top - down symmetry and if there are not too large elongations .
six features characterize the weak - steric - hindrance scenario .
first ( i ) , the plateaus @xmath193 for the reorientational correlators for odd @xmath1 are suppressed in comparison to what one would expect by interpolating or extrapolating the values for nearby even-@xmath1 plateaus ( tab .
[ erstetabelle ] ) .
most importantly ( ii ) , the critical amplitude @xmath298 is larger than @xmath502 , eq .
( [ einunddreissig ] ) , so that the canonical ordering of the @xmath194 for small @xmath1 , eq .
( [ vierundzwanzigb ] ) , is reversed .
third ( iii ) , the percentage of the decay of the reorientational correlators @xmath169 , which can be explained by the leading- plus next - to - leading - order asymptotic formulas for the @xmath2-relaxation is larger for @xmath1=1 than for @xmath1=2 ; as is indicated by the vertical lines in fig .
the structural relaxation of the reorientations is dominated by large - angle flips ( iv ) , as shown in fig .
[ sechsteabb ] for the dumbbell with elongation @xmath174 .
the @xmath0-relaxation - time scale for @xmath1=1 is smaller than for @xmath1=2 , eq .
( [ zweiunddreissig ] ) , ( v ) so that the canonical order of the @xmath0-relaxation scales , eq . ( [ achtundzwanzig ] ) , is reversed .
this can cause the @xmath169-versus-@xmath319 graphs for @xmath1=1 to cross the graphs for @xmath1=2 , as is shown for the @xmath174 results in figs .
[ dritteabb ] and [ vierteabb ] . finally ( vi ) , for distance parameters @xmath455 , where the conventional @xmath169-versus-@xmath503 plot exhibits the @xmath0-relaxation scaling law for @xmath1=2 , the correlators for @xmath1=1 do not show the validity of the superposition principle , as is demonstrated in fig .
[ achteabb ] . nor does the scale @xmath435 , defined as the time for a 50% decay of the @xmath0-relaxation correlator , exhibit the power law behavior with the correct exponent @xmath454 as is shown in the inset of fig .
[ neunteabb ] .    a side remark concerning a molecular - dynamics study of the evolution of glassy dynamics in a lennard - jones - dumbbell liquid by kmmerer et al .
@xcite might be in order .
it was reported that the correlators dealing with translational degrees of freedom and also for the ones for the reorientational dynamics for angular index @xmath504 could be interpreted qualitatively within the universal asymptotic mct formulas .
however , the evolution of the dipole correlators did not fit into the standard mct pattern .
it was found that @xmath505 and @xmath506 .
a drastic violation of the @xmath0-scaling law was noted quite similar to what is exhibited in the lower left panel of fig .
[ achteabb ] .
the scale @xmath435 exhibited a deviation from the asymptotic law @xmath507 , but a fit by @xmath508 with @xmath458 as discussed in the inset of fig .
[ neunteabb ] was possible .
these simulation results for @xmath1=1 differ from those for water simulations @xcite as well as from the experimental findings for propylene carbonate quoted in fig .
[ ersteabb ] . however , they agree with the features ( ii ) , ( v ) , and ( vi ) specified in the preceding paragraph . moreover , the property ( iv ) concerning the large angle flips is also obtained in ref .
accidently , the @xmath352 results in the lower panel of fig . [ fuenfteabb ] show that the minimum position of the @xmath1=1 spectrum exceeds that of the @xmath1=2 spectrum by nearly one order of magnitude , in agreement with the corresponding finding in ref . @xcite .
furthermore , the @xmath0-peak variation with @xmath121 shown in the lower left panel of fig . [ fuenfteabb ] is in qualitative agreement with that reported in ref . @xcite . in view of these observations
it does not seem impossible , that the scenario studied in refs .
@xcite fits into the framework of the ideal mct .
however , it is not clear , whether or not the results of ref .
@xcite can be explained by our theory for type - a precursors of a dilute solution of molecules .
first , the simulation results for the dumbbell liquid do not exhibit a particular decrease of @xmath277 relative to @xmath278 .
second , the @xmath2-relaxation scaling has not been documented for the dumbbell liquid and so it is unclear , whether or not the feature ( iii ) holds for that case .    summarizing
, it shall be emphasized that all qualitative features for the evolution of the structural relaxation studied in this paper have been explained by means of the formulas for the leading - asymptotic expansions and their leading - order - correction formulas for the bifurcation scenario . in this sense , these asymptotic formulas can be considered as the essence of mct . however , in order to explain the characteristic ( @xmath1=1)-versus-(@xmath1=2 ) differences for the relaxation patterns , it is necessary to also understand the general trends of the nonuniversal parameters with variations of wave - vector @xmath54 and angular - momentum index @xmath1 . and this requires to use mct as a microscopic theory based on the knowledge of the equilibrium structure .
we thank the authors of refs .
@xcite for the permission to use their data files and m. fuchs for many stimulating discussions .
we kindly thank f. sciortino for the permission to cite his unpublished simulation results for the reorientational dynamics of water .
we acknowledge gratefully helpful critique and suggestions for improvements of the manuscript written to us by h.z .
cummins , r. pick , and r. schilling .
this work was supported by verbundprojekt bmbf 03-g05tum .