animal packing in social aggregations is of fundamental interest in ecology , and their conformations have been extensively studied @xcite .
pioneer tank observations of pelagic fishes ( herring , sprat and mackerel ) were conducted by @xcite , and the `` persistently uniform density '' of a school was noted . from experiments on schools of saithe _ pollachius virens _ , herring _ clupea harengus _ and cod _ gadus morhua _ cruising in a ten meters circular gantry tank , @xcite validated that all the fish in a school occupy a volume proportional to @xmath2 , where @xmath3 is the number of fish and @xmath4 is the mean nearest - neighbor distance ( approximate to the mean fish - body length ) in the school .
@xcite reported from field observations that the number densities of herring schools are almost independent of the dimensional size of school but are an order of magnitude lower than the packing densities observed when schooling in small tanks .
tauti and colleagues @xcite had experimentally shown that a fish school can be viewed as a group of inorganic particles and treated with the methods of physics .
@xcite explicitly applied physical equations to such fish schools . recently
, a number of theoretical and numerical models of schooling and flocking behavior have been studied by biologists , mathematicians and physicists @xcite . as for the animal group geometry , many models for social aggregations , however , predict that densities increase as the group size ( in number ) increases [ overviews of such models are given in refs .
@xcite mathematically tackled this problem of constant density and revealed the condition for a well - spaced group , i.e. what class of mutual interactions results in a relatively constant individual distance in the interior of the aggregate .
more recently , by means of underwater acoustics , the school sizes ( in number or biomass of fish ) have been quantitatively measured for different values of the dimensional size of schools in the wild .
precise data on conformations of large - sized schools of pelegic fishes became available rapidly , which were extremely helpful in elucidating a certain geometric law , i.e. the relation between dimensional and biomass sizes of pelagic fish schools , bringing about some important changes in our viewpoints .
misund and colleagues @xcite found that the power - law scaling in dimension - to - biomass relationship exists , and is robust across a broad range of pelagic species as well as across diverse environments .
they demonstrated that if the biomass @xmath3 in a school is , say , doubled , the cross - sectional area of a school is increased by a factor @xmath5 , i.e. @xmath6 and that the exponent @xmath7 looks universal , reading @xmath8 ( from the field data @xmath7 ranges from @xmath9 to @xmath10 ) .
the geometric law they found implies that the mean density of a school scales as @xmath11 in three dimensions of space and the conformation of social aggregations swells ( i.e. @xmath12 ) .
such a relationship as eq.([eqn : misund - school - geometry ] ) with @xmath13 has been utilized for the dimensions to biomass conversion ( e.g. transforming the school diameter @xmath14 to the number @xmath3 of fish ) in analyzing school - size distributions @xcite .
it might come as a surprise that packing densities decrease as the group size increases , contrary to previous observations and predictions .
laboratory observations for school geometry and internal structure have been made exclusively in small tanks , which generally show the constant density of fish that the school volume is proportional to the number of individuals . in all set - ups , the factors resulting in homogeneous ,
cohesive school may be especially pronounced , and it therefore seems that the quantified structure is skewed .
@xcite claimed that homogeneity in schools has been over - emphasized . _ in - situ _ observations of herring and sprat schools with a high resolution sonar revealed that the packing structure within the schools is rather heterogeneous @xcite .
this has been confirmed by measurements of free - swimming schools using photography and high - resolution echo integration , which showed that the packing density distribution in capelin _ mallotus villosus _ and clupeoid schools varies considerably @xcite .
regions of high density are usually found within the schools , and even empty vacuoles have been recorded .
@xcite observed by repeated echo integration that there were large variations in internal packing density of herring schools ( i.e. high - density regions or empty lacunas within a school ) , and that the packing density structure was quasi - stationary . besides artificial environments in small tanks , the discrepancy in former observations could have been caused by too small numbers of fish in the schools .
for instance , in @xcite , @xmath3 takes a few tens of fish .
since such a geometric law above is always defined only in a certain limit @xcite , the scaling in the dimension - to - biomass relationship is expected to hold for large - sized schools of pelagics .    in this paper , the exponent @xmath7
is estimated according to the established universal scaling law in the school - size distribution of pelagic fishes @xcite : choose the suitable value of @xmath7 to achieve the best data - collapse on the size distributions in terms of the school dimension .
notice that the dimension - to - biomass relationship is a property of the single school and the scaling exponent @xmath7 is determined by the behavioral algorithm of fish schooling at individual level , while the scaling in the school - size distribution emerges from the inter - school interactions at population level ( i.e. a global property of the interacting school system ) .
animal group size is a focal issue in ecology that , in contrast to scaling , has introduced a single preferred size ( i.e. optimal or compromise size ) for any organism living in groups @xcite .
figure  [ fig:1 ] shows an example of the histogram of school dimensions of japanese sardine _
sardinops melanostictus_. a peak frequency and a right skew are typical of pelagic fishes @xcite
. their linear dimension , e.g. the vertical thickness , of the school falls into a certain range below a few tens of meters . one fish may not be the right atomic unit in schooling , since field observations suggest that no school exists under a certain minimal size .
this may cause binned data of school sizes to exhibit the fake peak frequency . @xcite and @xcite have ignored slowly decaying ( or fat - tailed ) distributions , including scaling laws , in such data of school sizes short ranged with a fake peak frequency ( the possibility of power laws was already presented in their models but not exploited ) . a traditional , widely used gauss statistics says that , for the data from @xcite shown in fig.[fig:1 ] , finding sardine schools ranging from 18 to 20 meters in vertical thickness should only occur about once every @xmath15 detections of schools [ for detail , consult @xcite ] .
in other words , it is not the real world !
aquatic observations actually say that finding such schools occurs about once every 500 detections .
the probability that such schools are found is @xmath16 times as large !    ) may be stray ones . ,
width=245 ]    the possibility of scaling in such distributions was found by @xcite .
fat - tailed distributions have been found to quite generally characterize the size heterogeneity of pelagic fish schools @xcite .
only lately , @xcite showed that the school - size distribution of pelagic fishes has the universal functional form of a power - law decay and a crossover towards the exponential distribution .
scaling laws have been found to characterize many patterns in ecological systems @xcite .
there is an empirical rule about the relationship between the areas ( @xmath17 ) of island and the numbers ( @xmath18 ) of species on them : whether counting birds , butterflies , plants or copepods in ponds , a tenfold smaller area contains roughly half the species .
this can be fitted by a power - law function @xmath19 with the exponent @xmath20 @xcite .
the species - area relationship is widely used in rough assessments of likely future rates of species extinction ( because observed rates of tropical deforestation are equated to loss of habitat area , we can assess the annual production of species committed to extinction ) .
so , a more secure understanding of such relationships has great practical importance @xcite .
@xcite applied finite - size scaling @xcite to test the hypothesis of scaling invariance in the species - area relationship resting on a model proposed by @xcite for the species abundance distribution .
following the conventional notation in the physics literature , the system sizes are hereafter denoted by @xmath17 .
the finite - size scaling ( fss ) hypothesis assumes that the fraction of objects with size @xmath21 for a finite system of size @xmath17 is written , with a scaling exponent @xmath22 , as @xmath23 where @xmath24 is a universal scaling function .
@xmath25 reads , for instance , the species - abundance distribution , which defines the probability that any given species on a census patch of area @xmath17 has @xmath21 individuals , where @xmath26 is the expected number of species with @xmath21 individuals .
@xmath27 is a measure of the width of the probability distribution @xmath25 , e.g. the mean or the standard deviation of the distribution .
equation  ( [ eqn : fss - hypothesis ] ) expresses the principle that the behavior of the finite system , i.e. @xmath25 , is controlled by the ratio @xmath28 .
we test the fss hypothesis by data collapsing : when @xmath29 is plotted versus @xmath30 , if the distributions across different @xmath17 s fall on a single curve ( e.g. collapsing distinct abundance distributions measured in different areas and with different total numbers of individuals onto one scaling curve ) , then we should be able to identify a scaling function @xmath31 . in order to determine the exponent @xmath22
, we try a best - fitting procedure such that the family of curves @xmath25 collapses onto a single curve as well as possible .
we then obtain other scaling exponents ( e.g. @xmath32 of the species - area relationship ) resting on the scaling relation between scaling exponents @xcite . by noting that the total number of individuals of all species is equal to @xmath33 , the species - area relationship leads to a scaling relation @xmath34 assuming that the total number of individuals of all species is proportional to area @xmath17 @xcite .
the universal scaling is fundamental to data collapse .
the power - law scaling emerges as ubiquitous properties in ecology . in statistical analysis in ecology ,
data collapsing across different environments or species has been , however , observed only lately in dynamics of breeding bird populations @xcite , food web structure @xcite , microbial body - mass spectra in marine ecosystems @xcite , school - size distributions of pelagic fishes @xcite , and in the context of ecological economics , capture fisheries productions of countries @xcite .
in order to understand the geometric properties of school configurations , the school biomass @xmath3 must be measured for different values of dimensional size @xmath14 , and we must compare them .
misund and colleagues extensively performed the simultaneous observations of the values of @xmath14 and @xmath3 for pelagics by the following two methods : ( i )  two - dimensional ( cross - sectional ) acoustic measurements ( unit in square meters ) and subsequent purse seine capture ( unit in tonnes ) of schools @xcite , and ( ii )  acoustic measurements of the three - dimensional structures and backscattered echo energy @xcite .
the surveys were conducted on stocks of anchovy _ engraulis capensis _
, herring _ clupea harengus _ , horse mackerel _
trachurus trachurus capensis _ ,
mackerel _ scomber scombrus _ , pilchard _ sardinops ocellatus _ , round herring _ etrumeus whiteheadi _ , saithe _ pollachius virens _ , sardine _ sardinops sagax _ , and sprat _ sprattus sprattus _ , in different seasons and geographic regions ( the barents sea , the north sea , the norwegian sea , the northeastern atlantic and off namibia , and off the coast of south africa ) .
they found that there is some biological universality in the dimensions - to - biomass relationships , eq.([eqn : misund - school - geometry ] ) .
their finding implies that the radius @xmath14 of the school scales as @xmath35 in a statistical sense .
@xmath3 , denoting the school biomass , is hereafter defined by the number of fish in a school .
now i propose another way to establish the scaling law for school dimensions , eq.([eqn : scaling - school - geometry ] ) , by applying universal scaling law in fish school - biomass distributions @xcite . the biomass distributions
@xmath36 follow a power law with exponent @xmath37 up to a cut - off size @xmath38 , @xmath39 where @xmath40 is a crossover scaling function with a strong drop for @xmath41 , and the cut - off size ( crossover size from power - law to exponential decay ) is calculated from the biomass histogram data @xmath42 , @xmath43 where histogram bins are chosen with width @xmath44 .
the dimension data of fish schools are binned with width @xmath45 , giving the set of frequencies @xmath46 . from eqs.([eqn : scaling - school - geometry ] ) and ( [ eqn : biomass - distribution ] ) , the distribution of geometric dimensions of fish schools is represented as @xmath47 where @xmath48 and @xmath49 therefore , the school - dimension distribution follows a power law with the same exponent `` @xmath50 '' as the school - biomass distribution .
the power - law distribution @xmath51 is truncated at a cut - off size @xmath52 , which is also calculated from histogram data of school geometric dimensions , @xmath53^{\nu}.    \label{eqn : depend - cutoff - nu}\ ] ] the following normalizations are adopted for the scaling function @xmath54 and the histogram data of geometric dimensions of fish schools , because the cut - off size @xmath38 is proportional to the total number of fish in the school system [ @xcite ; see also eq.([eqn : cutoff - population ] ) ] : @xmath55 and @xmath56 respectively .
fat - tailed school - size distributions are necessarily truncated because the population is finite . since the size @xmath57 depends on the exponent @xmath7 following eq.([eqn : depend - cutoff - nu ] ) , so that we can determine the value of @xmath7 by evaluating the location of the cut - off in the power - law distribution @xmath51 .
the ordinary least squares regressions might , however , lead to a `` wrong '' estimation of the exponent @xcite .
i make use of the data collapse to extract the `` right '' exponent . from eqs.([eqn : dimension - distribution ] ) and ( [ eqn : dimension - distribution - normalization ] ) , when @xmath58 is plotted against @xmath59 with correct parameter @xmath7 , all the empirical data should collapse onto each other .
the power - law exponent of school - dimension distributions , @xmath7 , is then evaluated through data collapse .
let us search for the value of @xmath7 that places all the points most accurately on a single curve .
we use a set of histogram data of vertical dimension of japanese sardine _ sardinops melanostictus _ schools , from 22 acoustic surveys by @xcite off southeastern hokkaido for seven years , 19811987 .
@xcite reported that japanese sardine migrate as a huge - sized school in number from a few hundreds of thousands to a few million of fish . to obtain the best data collapse , the @xmath60-axis
is divided into bins ( fig.[fig:2]a ) , and for each bin two - dimensional variance @xmath61 is calculated , where @xmath62 and @xmath63 denote the standard deviation of the mean @xmath64 and @xmath65 , respectively .
the parameter @xmath7 is then estimated at value that minimize the mean of two - dimensional variance for the bins ( fig.[fig:2]b ) .
the mean of two - dimensional variance , @xmath66 , is a measure to determine the goodness of collapse @xcite .
a good data collapse can be obtained by using the value @xmath67 .
the resulting plot of empirical school data is shown in fig.[fig:2]a .
experimentally fitting the parameter @xmath7 to achieve a good data collapse , `` @xmath0 '' , is reminiscent of the flory value of the exponent in a power - law dependence of the coil radius of a polymer chain ( in three - dimensional solutions ) on the degree of polymerization @xcite .    .
( a )  scaled distribution of geometric dimensions of sardine schools .
@xmath68 is plotted versus @xmath69 with @xmath70 on double - logarithmic scale .
the bins are chosen equally spaced on a logarithmic scale as @xmath71 with @xmath72 .
the rectangle in gray reads the interval @xmath73 , i.e. the error @xmath63 on the mean @xmath65 ( indicated by the slit ) for each bin .
the solid line is a prediction of the mean - field theory @xcite .
( b )  the mean of two - dimensional variance , @xmath66 , versus the power - law exponent @xmath7 .
although @xmath66 shows noisy fluctuations , it takes a minimum around @xmath74 .
data from @xcite are analyzed .
, title="fig:",width=245 ] .
( a )  scaled distribution of geometric dimensions of sardine schools .
@xmath68 is plotted versus @xmath69 with @xmath70 on double - logarithmic scale .
the bins are chosen equally spaced on a logarithmic scale as @xmath71 with @xmath72 .
the rectangle in gray reads the interval @xmath73 , i.e. the error @xmath63 on the mean @xmath65 ( indicated by the slit ) for each bin . the solid line is a prediction of the mean - field theory @xcite .
( b )  the mean of two - dimensional variance , @xmath66 , versus the power - law exponent @xmath7 .
although @xmath66 shows noisy fluctuations , it takes a minimum around @xmath74 .
data from @xcite are analyzed .
, title="fig:",width=245 ]    here we see that the power - law regime of the distribution is too short , which misled @xcite and @xcite into overlooking the power - law distributions of school sizes .
notice that many power laws that appeared in the ecology literature span less than two orders of magnitude of scale @xcite .
the power law range of too few scales is not unique to ecology ; the largest numbers of power laws reported in the physical science are for small ranges @xcite .
the acoustic - survey data are converted into a school - biomass histogram as follows @xmath75 we now crosscheck the value of @xmath7 through finite - size scaling ( fss ) analysis of school - biomass distribution @xcite . since the finite population size causes the truncation of power - law distribution @xmath76 ,
there is a well - defined quantity @xmath77 which depends on the system population size . in order to characterize the finite size effects ,
fss hypothesis is used : the distribution function depends on @xmath3 only through the ratio @xmath78 , @xmath79 where @xmath31 is a universal function independent of fish population size .
the prefactor @xmath80 is required to ensure the normalization @xmath81 where @xmath82 $ ] now represents the probability distribution of school - biomass sizes . from the fss hypothesis ,
it is expected that when @xmath83 is plotted against @xmath78 with correct parameters @xmath22 and @xmath84 all the data collapse onto a single curve .
the power - law exponent of biomass distributions , @xmath85 , is then evaluated through fss analysis .
the value of @xmath86 is the estimate of the power - law exponent @xmath87 let us analyze a set of 22 acoustic - survey data of sardine schools ( same as fig.[fig:2 ] ) converted into biomass histograms by using eq.([eqn : biomass - dimension - convert ] ) with @xmath88 . in a simultaneous best - fitting procedure ( fig.[fig:3 ] ) , a good data collapse can be obtained by using the values @xmath89 and @xmath90 .
the power - law exponent derived from the fss collapse is @xmath91 .
the resulting plot is shown in fig.[fig:3]a .
the school - biomass distribution follows a power - law decay with exponent @xmath50 , and is truncated at the cut - off size of eq.([eqn : cut - off - size ] ) .
the fss collapse confirms the scaling laws for school sizes , eq.([eqn : scaling - school - geometry ] ) with @xmath88 and eq.([eqn : biomass - distribution ] ) with @xmath37 .    )
are converted by eq.([eqn : biomass - dimension - convert ] ) with @xmath88 . here
@xmath92 is plotted versus @xmath93 with @xmath94 .
two - dimensional variances [ same as eq.([eqn:2d - var ] ) ] are calculated for bins chosen equally spaced on a logarithmic scale as @xmath95 with @xmath96 .
the rectangle in gray is same as fig.[fig:2 ] .
the solid line is a prediction of the mean - field theory @xcite .
( b )  the region of the @xmath97-plane in which the minimum of the mean of two - dimensional variance exists .
a measure of data collapse for scaling , i.e. the mean of two - dimensional variance , @xmath66 , takes a minimum @xmath98 for the right choice of @xmath99 .
the minimum is found with the precision i.e. width of the minimum , @xmath100 in black region ( @xmath101 ) .
the values of the parameters lie in the intervals @xmath102 and @xmath103 , and therefore @xmath104 .
, title="fig:",width=245 ] ) are converted by eq.([eqn : biomass - dimension - convert ] ) with @xmath88 . here
@xmath92 is plotted versus @xmath93 with @xmath94 .
two - dimensional variances [ same as eq.([eqn:2d - var ] ) ] are calculated for bins chosen equally spaced on a logarithmic scale as @xmath95 with @xmath96 . the rectangle in gray is same as fig.[fig:2 ] .
the solid line is a prediction of the mean - field theory @xcite .
( b )  the region of the @xmath97-plane in which the minimum of the mean of two - dimensional variance exists .
a measure of data collapse for scaling , i.e. the mean of two - dimensional variance , @xmath66 , takes a minimum @xmath98 for the right choice of @xmath99 .
the minimum is found with the precision i.e. width of the minimum , @xmath100 in black region ( @xmath101 ) .
the values of the parameters lie in the intervals @xmath102 and @xmath103 , and therefore @xmath104 .
, title="fig:",width=245 ]    finite - size scaling techniques have been applied to ecology and shown to predict scaling relations between scaling exponents in relative abundance of species [ eq.([eqn : scaling - ecology ] ) ] @xcite , in dynamics of breeding bird populations @xcite , in biomass - size distributions of seston @xcite , and in exploitation of fish stocks @xcite . in school - size statistics
, we expect to have a scaling relation .
we now choose the normalization @xmath105 where @xmath106 denotes the total fish population in the school system ( @xmath107 gives the total number of schools ) .
since eq.([eqn : sys - size - school ] ) implies @xmath108 with a scaling exponent @xmath109 , the fss relation for the school - biomass distribution is written as @xmath110 accordingly , the normalization of eq.([eqn : norm - phi ] ) yields the scaling relation @xmath111 the best - fitting value in data collapse , @xmath112 , gives @xmath113 . as a consequence , the location of the cut - off in the power - law distribution of school sizes simply reads @xmath114 which is verified by means of extensive numerical simulations @xcite .
let us now investigate cohesive motion of schools in a three - dimensional space from the viewpoint of the behavioral algorithms which govern their formation and dynamics : attraction between neighbors maintains cohesion of the school ; a tendency to align with neighbors produces collective motion of the school .
the minimal model of cohesion is a linkage of neighbors consisting of harmonic spring , because the inter - fish distance follows a gaussian distribution @xcite .
using the relative coordinates to the center of school , the equation of motion of fish ( as noisy self - propelled particles ) in a large school of size @xmath3 is written by a one - body approximation as the following langevin equation @xcite : @xmath115 providing that individuals are sufficiently sensitive to behavior of their neighbors , where @xmath116 is the cohesive force ; @xmath117 denotes the strength of alignment ; and @xmath118 is a random perturbation of the velocity with strength @xmath119 and a @xmath120-correlated time dependence .
the attractive interaction force acting on one body due to a system of bodies is the neighbor - joining harmonic spring .
one essential approximation is to replace the many - fish problem by the problem of solving the motion of one fish in a certain self - consistent field .
the cohesive force is then written in the following form by the one - body approximation : @xmath121 where @xmath122 is the effective spring constant for the springs strung out from the center of the school to a fish : the number of consecutive springs joining the fish to the center of the school via other companions is proportional to @xmath123 .
the total `` elastic energy '' of inter - fish bonds in the school depends linearly on number @xmath3 of individuals in the school . the overall elastic energy after integration over a sphere of radius @xmath14 ( denoting the average radius of the school ) results in @xmath124 , which is derived from the ideal harmonic spring , eq.([eqn : cohesion ] ) .
therefore , the exponent @xmath125 will be determined self - consistently , @xmath126 the solution of the langevin equation  ( [ eqn : langevin ] ) with eq.([eqn : cohesion ] ) takes the following asymptotic forms @xcite : @xmath127
\approx     { 3\epsilon_v}/j ,    \label{eqn : fluctuation - dissipation}\ ] ] @xmath128 where the root - mean - square @xmath129 gives the average radius of the school , @xmath14 . from eq.([eqn : self - consistent - spring ] ) together with eq.([eqn : mean - square - radius ] ) , the self - consistent value of the exponent @xmath125 is obtained : @xmath130 as a consequence , the self - consistent cohesive force reads @xmath131 where @xmath4 denotes the effective inter - fish distance ( a constant independent of @xmath3 ) : @xmath132 . the average radius of the school is then given by @xmath133 .
note that eq.([eqn : fluctuation - dissipation ] ) is an example of a more general principle called the fluctuation - dissipation theorem @xcite .
notice that the above model has been developed in the absence of the excluded volume interactions .
social cohesion by harmonic spring interactions between neighboring fish reminds us of the classical picture of a polymer chain based on the bead - spring model @xcite .
the exponent @xmath7 close to the value @xmath0 may then be understood by taking the excluded volume effect @xcite .
such a school without excluded volume effect may be called the `` gaussian '' school in line with polymer physics .      from eq.([eqn :
langevin ] ) , the probability of the school radius being between @xmath134 and @xmath135 is given by the following @xcite : @xmath136 in a stationary state in three dimensions ( i.e. the position vector @xmath137 follows a gaussian distribution ) .
the gaussian school model considered above permits fish to occupy the same region in space . of course
this is a physical impossibility since each fish possesses its own finite volume .
therefore , in the `` excluded volume '' school , there are a number of gaussian school configurations which are disallowed due to the steric effect .
let @xmath138 be the probability that a gaussian school configuration , as counted in eq.([eqn : prob - gaussian - radius ] ) , is also allowable under the excluded volume condition .
we now calculate the probability that no overlaps occur when we place @xmath3 fish within a region of volume ( @xmath139 ) , which will lead to an estimation for @xmath138 .
the approach is due to @xcite .
letting @xmath140 be the volume which is effectively excluded to one fish by the presence of another ( @xmath141 ) , the probability that one particular fish will not overlap with another is then given by @xmath142 .
since there are @xmath143 possible combinations of pairs , the probability that no overlap occurs in all of these combinations is given by @xmath144 where @xmath145 and @xmath146 .
therefore , the probability distribution of the school radius @xmath134 can then be written as @xmath147 for the excluded volume school .
both @xmath148 and @xmath149 have a maximum at certain values of @xmath134 .
let us estimate the average size of the school radius in each model by calculating the positions of these maxima .
the maximum of @xmath148 occurs at @xmath150 .
the maximum of @xmath149 occurs at @xmath14 , which satisfies the following equation obtained by differentiating the logarithm of eq.([eqn : prob - excluded - radius ] ) : @xmath151 if @xmath152 , the second term on the left - hand side of eq.([eqn : diff - log ] ) can be neglected and hence @xmath153 the exponent @xmath0 is the very value extracted through the data collapse .
the scaling in the relationship between geometric dimensions and biomass of pelagic fish schools is analogous to that developed in polymer physics @xcite .
the result suggests that the dynamics are unexpectedly `` simple '' and depend primarily on common cohesive motion in these animate and inanimate systems .
though the above is a very rough theory of school conformation with the excluded volume effect , the overall statistical properties do not depend on the details of the model , which is a consequence of universality @xcite .      as indicated in the introduction , large internal variations in packing density
occur within a school .
@xcite pointed out that a source of substantial variation in internal school structure is the formation of subgroups .
such subgroups have been observed in saithe schools @xcite , and in schools of minnow _ phoxinus phoxinus _ @xcite and herring @xcite .
relatively independent movements of such clusters of individuals can open up empty spaces and cause large variation in school volume .
the polymer - chain analogue of the subgroup in school conformation is the `` blob '' @xcite : the polymer chain behaves as a series of blobs . based on the blob concept , for small value of the excluded volume @xmath140 , a subgroup within a school , with a number of @xmath154 of fish ,
must be nearly gaussian .
we see this from eq.([eqn : prob - excluded - radius ] ) when we find no effect of @xmath140 if @xmath155 .
there is a certain value of @xmath156 $ ] , beyond which excluded volume effects become important .
a single school will appear gaussian at scales @xmath157 where @xmath158 while at scales @xmath159 it will show excluded volume effects . according to the `` blob '' approach
, an @xmath3-sized school can be described as a cluster of @xmath160 subgroups . inside the subgroup the core repulsion ( by excluded volume effects ) is a weak perturbation leading to a gaussian state .
the school radius is then written as @xmath161 in the native state .
notice the plasticity of school geometry and internal structure .
the compact , dense packing , as observed in artificial environments like small tanks , may be caused by the very strong stress .
when the school is confronted by danger such as predators , the inter - fish distance decreases rapidly and all the vacuoles within the school collapse quickly @xcite .
this change may be closely related to the coil - globule transition in a polymer solution @xcite .    in real schools the nature of the short range interaction
is quite complicated like van der waals - type inter - molecular forces .
accordingly , all fish within a school interact via the two - body potential deformed from a parabolic potential of harmonic spring [ grouping forces have been expressed as the gradient of a potential function in previous model studies @xcite ] .
the potential change , denoted by @xmath162 , is given as the correction term to cohesive interaction potential between the @xmath163-th and the @xmath164-th fish : @xmath162 will include steric effects and also may involve weak attractions .
the correction factor to the probability distribution @xmath149 of the school radius @xmath134 , in place of eq.([eqn : prob - no - overlap-1 ] ) , is then given by the configuration integral @xmath165 where @xmath166 denotes the position vector of the @xmath163-th fish , and @xmath167 is integrated over the school configuration space . we can estimate @xmath138 by means of the cluster expansion @xcite . introducing a function of separation between pairs , @xmath168 , @xmath169 we integrate eq.([eqn : prob - no - overlap ] ) in the mayer cluster expansion over diagrams up to the second irreducible cluster , yielding @xmath170 with the excluded volume parameter @xmath171 and the other parameter defined by @xmath172 the actual configuration is decided by maximizing the distribution function @xmath173 , and this maximization condition leads to the equation determining the school radius @xmath14 , @xmath174 notice that , when the first irreducible cluster integral in mayer s expansion is only taken into account , the configuration integral , eq.([eqn : prob - no - overlap ] ) , is evaluated by eq.([eqn : prob - no - overlap-1 ] ) . in case that @xmath162 consists of a hard core repulsion and a short - range weak attraction ( i.e. at the distance of near - collision ( @xmath175 ) the attraction gradient between fish becomes steeper ) , eq.([eqn : second - virial ] ) is estimated as @xmath176 , expressing the dependence of @xmath140 on the velocity - fluctuation @xmath119 of fish , where @xmath177 and @xmath178 are constants @xcite . during predation threats or laboratory observations ,
the value of excluded volume parameter @xmath140 will change sign from positive to negative , provided that the velocity fluctuation @xmath119 is suppressed below a certain value @xmath178 .
decreasing @xmath119 below @xmath178 , the dimensional size of the school becomes much smaller than that of a gaussian school , as depicted in fig.[fig:4 ] . according to eq.([eqn : diff - log-2 ] ) the so called `` expansion factor '' @xmath179 is determined not by @xmath140 but by @xmath180 , and so if @xmath3 is large only a small change in @xmath119 will cause a big change in dimensional size .
for example , for school of @xmath16 fish , a small variation in @xmath119 will induce a dramatic change in the school radius .
equation  ( [ eqn : diff - log-2 ] ) states that when @xmath181 , @xmath179 is proportional to @xmath182 and hence the solution is as follows : @xmath183 which shows the closest packing of fish in a school , i.e. the constant density independent of the number of individuals .
this change may be called the swollen - dense packing transition of schools .    ) ] . here
@xmath184 is plotted as a function of a combined variable @xmath185 .
shown is the transition region @xmath186 . at @xmath187
the excluded volume parameter @xmath140 equals @xmath188 .
the repulsive excluded volume effect then balances the attractive forces between fish , and the school behaves as a gaussian school .
the property of the swollen - dense packing transition is quite different in the two regions separated by the critical value @xmath189 .
the transition between the the swollen conformation and the compact packing is gradual when @xmath190 ( dotted line ) , while it becomes extremely sharp and the transition is discontinuous when @xmath191 ( gray line ) .
the solid line shows the marginal transition behavior , and the crossover takes place when @xmath192 ( indicated by the arrow ) .
notice that the solution of eq.([eqn : diff - log-2 ] ) is multivalued [ depicted by a sigmoid curve ( white cutout ) ] , showing a first - order discrete phase change @xcite , when @xmath191
. , width=245 ]      the approach of treating fish schools as interacting particle systems naturally leads to the idea of applying successful methods of statistical physics to the description of moving together without a leader @xcite . in statistical physics
, the presence of non - trivial scaling is usually taken to mean that the dynamics are largely governed by simple geometric properties of the system and do not depend strongly on detailed properties of the system components @xcite .
as for the school geometry , the basic units are the effective inter - fish length @xmath4 and the number of fish , @xmath3 .
we assemble @xmath193 neighboring fish in a school into @xmath194 groups , e.g. in the sense of the voronoi - dirichlet diagram @xcite .
the length between centers of neighboring subgroups is defined by @xmath195 with an exponent @xmath196 .
after such a scale change as @xmath197 the macroscopic quantities which determine the overall properties of the school satisfy @xmath198 where @xmath199 is an exponent which depends on the quantity under consideration .
for this to hold true for arbitrary @xmath193 , the function @xmath200 must take the scaling form @xmath201 where a new function @xmath202 is introduced .
this is the widom - kadanoff scaling law in critical phenomena @xcite .
for example , the average radius of the school , @xmath14 , should be unaltered under the transformation of eq.([eqn : scale - change ] ) : @xmath203 because of eq.([eqn : scaling - school - geometry ] ) , @xmath204 and we must have @xmath205 where @xmath206 in three dimensions .    as an example of the application of the above scaling law ,
let us consider the pair correlation function inside a school , where fish movement take place in a three - dimensional space .
a pair correlation function @xmath207 is defined as follows .
we pick one fish at random in the school , and we place it at origin
. then we ask , what is the number density of other fish at a ( directional ) distance @xmath137 from the first , and we average the result over all choices of the first fish ( and of directions ) .
the function @xmath207 has an integral which is just the total number of fish per school , @xmath208 .
hence , from dimensional analysis , we can write @xmath209 under the transformation that neighboring @xmath193 fish are grouped to form one subunit in a school , @xmath207 will be reduced by @xmath210 , since the pair correlation function is proportional to the number density of fish .
therefore @xmath211 then , the function @xmath207 obeys a simple scaling rule @xmath212 where @xmath213 is a dimensionless universal function , and @xmath214 is the dimensional factor . focusing on the limit @xmath215 , we can reach the form of @xmath207 by a simple argument . in a sphere of radius @xmath134
we have a certain number of fish , @xmath216 , related to @xmath134 by the excluded volume exponent : @xmath217 . the function @xmath207 scales like the density of fish in the sphere , @xmath218 which gives @xmath219 .
the possible observable property , eq.([eqn : pair - correlation - scaling ] ) , tells us that if we were to measure @xmath207 for fish schools of different size @xmath3 , there would be superposition of the curves obtained by plotting @xmath220 against @xmath221 .
thus this kind of scaling relation may be verified experimentally on large schools of pelagic fish .
the power - law scaling generally exists in dimension - to - biomass relationship of pelagic fish schools in nature . here
i have estimated the power - law exponent @xmath7 of the geometric relation , based on the dependence of the power - law regime of the school - dimension distribution @xmath51 on @xmath7 , i.e. eq.([eqn : depend - cutoff - nu ] ) .
we have tested whether the distributions @xmath51 are self - similar ( i.e. exhibit scaling ) , relying on both the power law for the dimension - to - biomass relationship of eq.([eqn : scaling - school - geometry ] ) and the fss relation in the power - law school - biomass distribution of eq.([eqn : biomass - distribution ] ) .
plotting the scaled histogram - data from the 22 _ in - situ _ observations , we have found that the 22 curves do indeed collapse onto each other ( fig.[fig:2]a ) , suggesting that @xmath51 follows a universal functional form [ eq.([eqn : dimension - distribution ] ) ] .
we have extracted the power - law exponent of dimension - to - biomass scaling relationship via a minimization of a measure to quantify the nature of scaling collapse .
next i have explained the value of the exponent @xmath7 , proposing the gaussian school model for the fish with excluded volume @xmath140 . by exploiting the analogy between fish - school and polymer conformations
, we have examined the behavioral algorithm governing the swollen conformation of large - sized schools .
we have seen that the exponent is modified strongly by the steric effect . what is universal in eq.([eqn : scaling - school - geometry ] )
is the exponent @xmath222 : it is independent of species as well as environmental conditions , and the same for all schools .
the constant that multiplies @xmath223 in eq.([eqn : scaling - school - geometry ] ) is non - universal [ @xmath224 with effective inter - fish distance @xmath4 , predicted from eq.([eqn : rouse - excluded - volume - effect ] ) ] , and depends on the details of interactions between fish which may vary with species and environmental conditions @xcite .    to understand the scaling law for school geometry , eq.([eqn : scaling - school - geometry ] ) with @xmath88 , it is essential to see what the value of @xmath7 represents .
obviously , for a regular object embedded into a @xmath225-dimensional euclidean space , eq.([eqn : scaling - school - geometry ] ) would have the form @xmath226 expressing the fact that the volume of a @xmath225-dimensional object grows with its linear size @xmath14 as @xmath227 .
contrary to an integer dimensionality , it has been shown that the packing structure within the schools is characterized by a non - integer ( i.e. fractal ) dimensionality of @xmath228 : the number of fish in a school of radius @xmath14 scales as @xmath229 . during the last decades of the twentieth century it has widely been recognized by researchers working in diverse areas of science
that many of the structures commonly observed possess a rather special kind of geometrical complexity ; the name `` fractal '' was coined by benoit mandelbrot for these complex shapes @xcite .
objects of biological origin are many times fractal - like @xcite .
fss plots of data lead to possibly important and practical insights in ecology . in the context of the fisheries mission ,
the demonstrated geometric relation between dimensions and biomass of pelagic fish schools , when applied to mapping pelagic schooling fish , will largely improve the precision in the fish stock assessment .
because of the linear dependence of the fish population on @xmath38 , as given by eq.([eqn : cutoff - population ] ) , the fish stock abundance can be inferred from an index ( i.e. cut - off size ) that can be determined directly from observations [ see eq.([eqn : cut - off - size ] ) ] .
99 aji , v. , goldenfeld , n. , 2001 .
fluctuations in finite critical and turbulent systems .
86 , 10071010 .
, doi:10.1103/physrevlett.86.1007 .
anderson , j. j. , 1981 . a stochastic model for the size of fish schools .
fish . bull .
79 , 315323 .
aoki , i. , 1980 .
an analysis of the schooling behavior of fish : internal organization and communication process .
ocean res .
tokyo 12 , 165 .
aoki , i. , inagaki , t. , 1988 .
photographic observations on the behavior of japanese anchovy _ engraulis japonica _ at night in the sea .
43 , 213221 .
azovsky , a. i. , 2000 .
concept of scale in marine ecology : linking the words or the worlds ?
web ecol . 1 , 2834 .
banavar , j. r. , green , j. l. , harte , j. , maritan , a. , 1999 .
finite size scaling in ecology .
83 , 42124214 .
beecham , j. a. , farnsworth , k. d. , 1999 .
animal group forces resulting from predator avoidance and competition minimization .
198 , 533548 .
binder , k. , heermann , d. w. , 1988 .
monte carlo simulation in statistical physics .
springer - verlag , berlin .
bonabeau , e. , dagorn , l. , 1995 . possible universality in the size distribution of fish schools .
e 51 , r5220r5223 .
bonabeau , e. , dagorn , l. , fron , p. , 1998 .
space dimension and scaling in fish school - size distributions .
a : math . gen .
31 , l731l736 .
bonabeau , e. , dagorn , l. , fron , p. , 1999 .
scaling in animal group - size distribution .
usa 96 , 44724477 .
breder , jr .
, c. m. , 1954 .
equation descriptive of fish schools and other animal aggregations .
ecology 35 , 361370 .
camacho , j. , guimer , r. , amaral , l. a. n. , 2002 .
robust patterns in food web structure .
88 , 228102 .
, doi:10.1103/physrevlett.88.228102 .
camazine , s. , deneubourg , j .- l . ,
franks , n. r. , sneyd , j. , theraulaz , g. , bonabeau , e. , 2001 .
self - organization in biological systems .
princeton university press , princeton .
chave , j. , levin , s. , 2003 . scale and scaling in ecological and economic systems .
resource econ .
26 , 527557 .
coetzee , j. , 2000 .
use of a shoal analysis and patch estimation system ( shapes ) to characterise sardine schools .
living resour .
13 , 110 .
cushing , d. h. , 1977 .
observations of fish schools with the arl scanner . rapp .
p .- v run
mer 170 , 1520 . de gennes , p .-
g . , 1975 .
collapse of a polymer chain in poor solvents .
j. de physique 36 , l55l57 . de gennes , p .- g . ,
1979 . scaling concepts in polymer physics .
cornell university press , ithaca , new york .
dill , l. m. , dunbrack , r. l. , major , p. f. , 1981 . a new stereophotographic technique for analyzing the three - dimensional structure of fish schools .
biol . fish .
dirichlet , g. l. , 1850 .
ber die reduktion der positiven quadratischen formen mit drei unbestimmten ganzen zahlen .
j. reine angew .
40 , 209227 .
doi , m. , 1996 .
introduction to polymer physics .
oxford university press , new york .
doi , m. , edwards , s. f. , 1986 .
the theory of polymer dynamics .
oxford university press , new york .
fisher , m. e. , barber , m. n. , 1972 . scaling theory for finite - size effects in the critical region .
28 , 15161519 .
flory , p. , 1953
principles of polymer chemistry .
cornell university press , ithaca , new york .
fron , p. , gerlotto , f. , soria , m. , 1992 .
changes in school structure according to external stimuli : description and influence on acoustic assessment . fish .
15 , 4566 .
fron , p. , misund , o. a. , 1999 .
dynamics of pelagic fish distribution and behaviour : effects on fisheries and stock assessment .
fishing news books , blackwell , osney mead , oxford .
goldenfeld , n. , 1992 .
lectures on phase transitions and the renormalization group .
vol.85 of frontiers in physics .
addison - wesley , reading , massachusetts .
graves , j. , 1976 .
photographic method for measuring spacing and density within pelagic fish schools at sea . u.s .
75 , 230234 . hamburger , d. a. , biham , o. , avnir , d. , 1996 .
apparent fractality emerging from models of random distributions .
e 53 , 33423358 .
halley , j. m. , hartley , s. , kallimanis , a. s. , kunin , w. e. , lennon , j. j. , sgardelis , s. p. , 2004 .
uses and abuses of fractal methodology in ecology .
. lett . 7 , 254271 . , doi:10.1111/j.1461 - 0248.2004.00568.x .
hara , i. , 1984 .
distribution and school size of japanese sardine in the waters off the southeastern coast of hokkaido on the basis of echo sounder surveys .
tokai reg . fish .
res . lab .
113 , 6778 .
hara , i. , 1986 .
stock assessment of japanese sardine in the waters off the southeast coast of hokkaido using line transect method .
sci . fish .
52 , 6973 .
hara , i. , 1990 .
stock assessment and migration of japanese sardine schools in the waters off the southeastern coast of hokkaido . ph.d .
thesis , kyoto university , kyoto , japan .
harte , j. , kinzig , a. , green , j. , 1999 .
self - similarity in the distribution and abundance of species .
science 284 , 334336 .
higashi , m. , yamamura , n. , 1993 .
what determines animal group size ?
insider - outsider conflict and its resolution .
142 , 553563 .
hoare , d. j. , couzin , i. d. , godin , j .-
j. , krause , j. , 2004 .
context - dependent group size choice in fish .
67 , 155164 .
, doi:10.1016/j.anbehav.2003.04.004 .
hori , j. , 1977 . langevin equation .
iwanami shoten , publ . ,
keitt , t. h. , stanley , h. e. , 1998 .
dynamics of north american breeding bird populations .
nature 393 , 257260 .
keitt , t. h. , amaral , l. a. n. , buldyrev , s. v. , stanley , h. e. , 2002 .
scaling in the growth of geographically subdivided populations : invariant patterns from a continent - wide biological survey .
b 357 , 627633 . , doi:10.1098/rstb.2001.1013 .
kubo , r. , 1966 . the fluctuation - dissipation theorem .
29 ( 1 ) , 255284 .
lillo , f. , farmer , j. d. , mantegna , r. n. , 2002 .
single curve collapse of the price impact function for the new york stock exchange .
preprint cond - mat/0207428 at @xmath230 lillo , f. , farmer , j. d. , mantegna , r. n. , 2003 .
master curve for price - impact function .
nature 421 , 129130 .
macarthur , r. h. , wilson , e. o. , 1967 .
the theory of island biogeography .
princeton university press , princeton .
mandelbrot , b. b. , 1982 . the fractal geometry of nature .
w. h. freeman and co. , new york .
may , r. m. , 1999 .
unanswered questions in ecology .
b 354 , 19511959 .
mayer , m. g. , mayer , j. e. , 1940 . statistical mechanics .
john wiley & sons , new york .
misund , o. a. , 1990 . sonar observations of schooling herring : school dimensions , swimming behaviour , and avoidance of vessel and purse seine .
p .- v run .
int . explor .
mer 189 , 135146 .
misund , o. a. , 1993a .
dynamics of moving masses : variability in packing density shape and size among herring sprat and saithe schools .
ices j. mar .
50 , 145160 .
misund , o. a. , 1993b .
abundance estimation of fish schools based on a relationship between school area and school biomass
living resour .
6 , 235241 .
misund , o. a. , aglen , a. , beltestad , a. k. , dalen , j. , 1992 .
relationship between the geometric dimensions and biomass of schools .
ices j. mar .
49 , 305315 .
misund , o. a. , aglen , a. , frns , e. , 1995 . mapping the shape , size , and density of fish schools by echo integration and a high - resolution sonar .
ices j. mar .
52 , 1120 .
misund , o. a. , aglen , a. , hamre , j. , ona , e. , rttingen , i. , skagen , d. , valdemarsen , j. w. , 1996 .
improved mapping of schooling fish near the surface : comparison of abundance estimates obtained by sonar and echo integration .
ices j. mar .
53 , 383388 .
misund , o. a. , coetzee , j. , 2000 .
recording fish schools by multi - beam sonar : potential for validating and supplementing echo integration recordings of schooling fish .
47 , 149159 .
, doi:10.1016/s0165 - 7836(00)00166 - 1 .
misund , o. a. , coetzee , j. , fron , p. , gardener , m. , olsen , k. , svellingen , i. , hampton , i. , 2003 .
schooling behaviour of sardine _ sardinops sagax _ in false bay , south africa .
25 , 185193 .
misund , o. a. , floen , s. , 1993 . packing density structure of herring schools .
196 , 2629 .
mogilner , a. , edelstein - keshet , l. , bent , l. , spiros , a. , 2003 .
mutual interactions , potentials , and individual distance in a social aggregation .
47 , 353389 .
, doi:10.1007/s00285 - 003 - 0209 - 7 .
morgan , m. j. , 1988 .
the effect of hunger , shoal size and the presence of a predator on shoal cohesiveness in bluntnose minnows , _ pimephales notatus _ rafinesque .
32 , 963971 .
niwa , h .- s . , 1994
. self - organizing dynamic model of fish schooling .
171 , 123136 .
niwa , h .- s . , 1996a .
newtonian dynamical approach to fish schooling .
181 , 4763 .
niwa , h .- s . , 1996b . mathematical model for the size distribution of fish schools .
computers math .
32 ( 11 ) , 7988 .
niwa , h .- s . , 1998 .
school size statistics of fish .
195 , 351361 .
niwa , h .- s . , 2003 .
power - law versus exponential distributions of animal group sizes .
224 , 451457 .
, doi:10.1016/s0022 - 5193(03)00192 - 9 .
niwa , h .- s .
space - irrelevant scaling law for fish school sizes .
228 , 347357 .
, doi:10.1016/j.jtbi.2004.01.011 .
niwa , h .- s . , 2004b .
exploitation dynamics of fish stocks .
meeting abstracts of the physical society of japan 59 ( 2 ) , part 2 , p.201 .
okubo , a. , chiang , h. c. , 1974 .
an analysis of the kinematics of swarming of _ anarete pritchardi _
kim ( diptera : cecidomyiidae ) .
16 , 142 .
okubo , a. , levin , s. a. , 2001 . diffusion and ecological problem : modern perspectives .
springer - verlag , new york .
parr , a. e. , 1927 . a contribution to the theoretical analysis of the schooling behavior of fishes .
bulletin of the bingham oceanographic collection 1 ( 1 ) , 132 .
parrish , j. k. , hamner , w. m. ( eds . ) , 1997 . animal groups in three dimensions .
cambridge university press , cambridge .
partridge , b. l. , 1980 .
the effect of school size or the structure and dynamics of minnow schools .
28 , 6877 .
partridge , b. l. , 1981 .
internal dynamics and interrelations of fish in schools .
144 , 313325 .
partridge , b. l. , 1982 .
the structure and function of fish schools .
246 , 9099 .
pitcher , t. j. , 1973 .
the three - dimensional structure of fish schools in the minnow , _
phcxinus phcxinus_. anim .
21 , 673686 .
pitcher , t. j. , parrish , j. k. , 1993 .
functions of shoaling behaviour in teleosts , in : pitcher , t. j. ( ed . ) , behaviour of teleost fishes .
chapman & hall , london , ch.12 , pp.363439 .
pitcher , t. j. , partridge , b. l. , 1979 . fish school density and volume .
54 , 383394 .
preston , f. w. , 1962 .
the canonical distribution of commonness and rarity : part i , and part ii .
ecology 43 , 185215 , and 410432 .
ptitsyn , o. b. , kron , a. k. , eizner , y. y. , 1968 .
the models of the denaturation of globular proteins .
i. theory of globula - coil transitions in macromolecules .
c 16 , 35093517 .
pulliam , h. r. , caraco , t. , 1984 . living in groups
: is there an optimal group size ?
, in : krebs , j.r
. , davies , n.b .
( eds . ) , behavioural ecology . an evolutionary approach .
blackwell science , oxford , pp . 122147 .
rinaldo , a. , maritan , a. , cavender - bares , k. k. , chisholm , s. w. , 2002 .
cross - scale ecological dynamics and microbial size spectra in marine ecosystems .
b 269 , 20512059 . , doi : 10.1098/rspb.2002.2102 .
rouse , p. e. , 1953 . a theory of the linear viscoelastic properties of dilute solutions of coiling polymers .
21 , 12721280 .
stanley , h. e. , 1995 .
power laws and universality .
nature 378 , 554 .
symons , p. e. k. , 1971 .
spacing and density in schooling threespine sticklebacks ( _ gasterosteus aculeatus _ ) and mummichog ( _ fundulus heteroclitus _ ) .
28 , 9991004 .
tauti , m. , miyosi , k. , 1929 . on the movement of a fish - group .
i. on the direction of a fish - group on passing the end of a wall along which it has moved .
j. imperial fish .
24 , 7376 .
tauti , m. , hudino , t. , 1929 . on the movement of a fish - group .
the case when fishes were barred from their progress by a wall
. j. imperial fish .
inst . jpn .
24 , 7679 .
tauti , m. , yasuda , h. , 1929 . on the movement of a fish - group .
the case when fishes were barred from their progress by a net . j. imperial fish .
24 , 8086 .
tauti , m. , yasuda , h. , 1930 . on the movement of a fish - group .
the shape of the entrance of fishing net and the fish - group moving in and out . j. imperial fish .
25 , 4150 .
tauti , m. , yasuda , h. , 1933a .
the direction in which a fish - group goes on after having met with a net .
2 , 1719 .
tauti , m. , yasuda , h. , 1933b .
how will a fish - group be attracted ?
2 , 124127 .
uhlenbeck , g. e. , ford , g. w. , 1962 .
the theory of linear graphs with applications to the theory of the virial development of the properties of gases , in : de boer , j. , uhlenbeck , g. e. ( eds . ) , studies in statistical mechanics .
north - holland , amsterdam , vol.1 , ch.2 , pp.119211 .
vicsek , t. , 2001 . a question of scale .
nature 411 , 421 .
vicsek , t. ( ed . ) , 2001 .
fluctuations and scaling in biology .
oxford university press , oxford .
voronoi , g. m. , 1908 .
nouvelles applications des paramtres continus  la thorie des formes quadratiques deuxime mmoire : rechrches sur les parallllodres primitifs .
j. reine angew .
134 , 198287 .
voronoi , g. m. , 1909 .
deuxime mmoire : recherches sur les parallloedres primitifs
. j. reine angew .
136 , 67181 .
warburton , k. , lazarus , j. , 1991 .
tendency - distance models of social cohesion in animal groups .
150 , 473488 .
wilson , k. g. , 1983 .
the renormalization group and the critical phenomena .
55 , 583 - 600 .


over the last decades , efforts have been made to get a better understanding of non - gaussian time series in the time and frequency domains .
in particular , time series whose marginal distributions exhibit power law behavior have attracted a lot of attention .
the need for such models arises from applications in areas as diverse as insurance , geophysics , finance and telecommunications .
infinite fourth moments are not untypical for series of daily log - returns from exchange rates , stock indices , and other speculative prices , whereas infinite second moments can be observed in time series from insurance such as for windstorm , industrial fire and earthquake insurance @xcite .
infinite first moments are typical for the marginal distribution of the magnitudes of earthquakes @xcite .
infinite variances are observed for the sizes of teletraffic data in the world wide web @xcite ; see also the recent books @xcite .
classical time series analysis deals with the second ( or higher ) order moment structure of a stationary sequence .
heavy - tailed modeling requires , in addition , that one takes into account the interplay between the dependence structure and the tails of the series .
an important task is to understand the classical statistical estimators and test procedures when big shocks to the underlying system are present . when the marginal distributions have infinite variance , the notions of autocovariance , autocorrelation and spectral distribution lose their meaning .
however , various studies over the last twenty years have shown that the analysis of linear processes @xmath3 , @xmath4 , with heavy - tailed i.i.d .
_ innovations _ @xmath5 and constant coefficients @xmath6 is very similar to classical ( finite variance ) time series analysis , where notions such as autocovariances and spectral density are defined only in terms of the @xmath7 s and the innovation variance @xmath8 .
most estimators and test statistics from classical time series analysis can be modified insofar that one considers self - normalized ( or studentized ) versions of them and for these versions , an asymptotic theory exists which parallels the classical theory with gaussian limit processes .
in contrast to the latter theory , the limits involve infinite variance stable distributions and processes @xcite .
one of the main goals of classical ( finite variance ) time series analysis is to study the spectral properties of the linear process @xmath9 . in this context , the _ periodogram _
@xmath10 , \ ] ] plays a prominent role as an estimator of the spectral density . numerous estimation and test procedures
are based on this statistic and integrated versions of the form @xmath11 for appropriate classes of real - valued functions @xmath12 on @xmath13 $ ] . in applications ,
the class @xmath14}\dvt x\in[0,\curpi]\}$ ] is most important .
the resulting integrated periodogram is a process indexed by @xmath15 $ ] . under general conditions , @xmath16})$ ] converges uniformly with probability 1 to the function @xmath17 , where @xmath18 , \end{aligned}\ ] ] and @xmath19 is the corresponding _ power transfer function_.
the latter is one of the essential building blocks of the _ spectral density _ of the stationary process @xmath9 :
@xmath20 .\ ] ] this is the fourier series based on the _ autocovariance function _ @xmath21 since @xmath16})$ ] estimates the spectral distribution function of the stationary process @xmath9 , it has been used for a long time as the _ empirical spectral distribution function _ , both as an estimator and a basic tool for constructing goodness - of - fit tests for the underlying spectral distribution function @xcite .
since the limit process of the properly centered and normalized process @xmath16})$ ] depends on the ( in general unknown ) spectral density @xmath22 , bartlett @xcite proposed to consider the integrated periodogram  based on @xmath23}/f_x(x)\dvt x\in[0,\curpi]\}$ ] . under general conditions , this process converges uniformly , with probability 1 , to the function @xmath24 and the limit process can be shown to be independent of the coefficients of the linear process , but depends on the fourth moment of @xmath25 .
more generally , weighted integrated periodograms based on the classes @xmath26 } g(x)\dvt x \in[0,\curpi]\}$ ] for suitable functions @xmath27 are used to estimate the spectral density or to perform various tests concerning the spectrum of the underlying stationary sequence @xcite .
the weighted integrated periodogram is also the basis of one of the classical estimators for fitting arma and fractional arima models .
this method goes back to early work by whittle @xcite . in this context , one considers the functional @xmath28 , @xmath29 , where @xmath30 is a class of spectral densities indexed by a parameter @xmath31 .
the _ whittle estimator _
@xmath32 of the true parameter @xmath33 is the minimizer of @xmath28 over @xmath34 .
this estimation technique is one of the backbones of quasi - maximum likelihood estimation in parametric time series modeling .
the whittle estimator is known to be asymptotically equivalent to the corresponding least - squares and gaussian quasi - maximum likelihood estimators @xcite .
when proving the asymptotic normality and consistency of @xmath32 , one has to study the properties of the sequence @xmath35 which , again , can be considered as a weighted integrated periodogram indexed by a class of functions .
the above examples have in common that one always considers a stochastic process @xmath36 for some class of functions . in all cases ,
one is interested in the asymptotic behavior of the process @xmath37 , uniformly over the class @xmath38 .
this is analogous to the case of the empirical distribution function indexed by classes of functions .
general references in this context are the monographs @xcite . early on
, this analogy was discovered by dahlhaus @xcite , who gave the uniform convergence theory for @xmath37 under entropy and exponential moment conditions .
the almost sure and weak convergence theory under entropy and power moment conditions was given in @xcite . a recent survey of nonparametric statistical methods related to the empirical spectral distribution indexed by classes of functions is @xcite .
it is the aim of this paper to develop an analogous weak convergence theory for heavy - tailed stationary processes .
we will understand ` heavy - tailedness ' in the sense of infinite variance of the marginal distributions .
our focus will be on linear processes @xmath9 with i.i.d .
symmetric @xmath0-stable ( s@xmath0s ) innovations @xmath39 for some @xmath40 .
recall that a random variabl @xmath41 has a _
symmetric stable distribution _
@xmath42 if there are parameters @xmath43 , and @xmath44 such that its characteristic function has the form latexmath:[$e \mathrm{e}^{\mathrm{i } t y_{\alpha } } = \mathrm{e}^{-\sigma ^\alpha    of @xmath25 . for @xmath46 ,
the random variable @xmath25 is known to have infinite variance @xcite .
much of the theory given below depends on tail estimates for random quadratic forms in i.i.d .  infinite variance random variables .
such results are available for i.i.d
.  stable sequences .
although it seems feasible that the theory can be extended to the more general class of processes whose innovations have regularly varying  tails , we do not attempt to achieve this goal
. the price would be more technicalities and the gain would be negligible .
we intend to show how the classical ( finite variance ) tools and methods have to be modified in the infinite variance stable situation which can be considered as a boundary case of the classical one when some of the innovations assume extremely large values . by now , there exists quite a clear picture concerning the asymptotic theory of the sample autocovariances , the periodogram and its integrated versions when the innovation sequencein a linear process has infinite variance ; see @xcite , chapter 7 .
in addition to the latter reference , goodness - of - fit tests for heavy - tailed processes ( corresponding to the class @xmath47 ) were considered for short- and long - memory linear processes @xcite , and whittle estimation for infinite variance arma and farima processes was also studied @xcite .
the paper is organized as follows . in section [ sec:1 ]
, we introduce some useful notation for the integrated periodogram .
our main goal is to prove the weak convergence of the integrated periodograms indexed by suitable classes of functions .
we achieve this goal for an i.i.d .  sequence in section  [ sec : iid ] , first by showing the _ convergence of the finite - dimensional distributions _
( section  [ s1 ] ) , then the _ tightness_. the conditions and methods are rather different in the cases @xmath1 ( section  [ s2 ] ) and @xmath2 ( section  [ s3 ] ) .
the case @xmath1 is treated in the more general context of random quadratic forms with toeplitz coefficient matrices satisfying some summability condition .
the case @xmath2 requires entropy conditions and the corresponding techniques . in section  [ sextra ]
, we extend the limit theory for the integrated periodograms from an i.i.d .
sequence to linear processes .
the contains some auxiliary results concerning tail estimates of random quadratic forms in stable random variables .
the weak convergence results of this paper might also be of separate interest in the context of infinite variance random quadratic forms .
the theory for such quadratic forms is not well studied .
we also refer to an extended version of this paper @xcite which covers the class of stochastic volatility processes with regularly varying marginal distributions .
the following decomposition of the periodogram is fundamental : @xmath48 where @xmath49 , @xmath50 , denotes the _ sample autocovariance function _ of the sample @xmath51 .
note that the definition of @xmath52 deviates slightly from the usual one where the @xmath53 s are centered by the sample mean .
however , for the theory given below , this centering is not essential . centering with the sample
mean @xmath54 is not the most natural choice when dealing with infinite variance processes . in what follows
, we will frequently make use of the _ self - normalized periodogram _
@xmath55 where @xmath56 , @xmath50 , denotes the _ sample autocorrelation function _ of @xmath51 . in view of ( [ decomp ] ) ,
we can rewrite @xmath57 as follows : @xmath58 where @xmath59 are the _ fourier coefficients _ of @xmath60 .
we also introduce the self - normalized version of @xmath61 : @xmath62
in this section , we study the limit behavior of the integrated periodograms @xmath61 indexed by classes of functions for an i.i.d .  s@xmath0s sequence with @xmath40 . in section  [ s1 ]
, we consider the convergence of the finite - dimensional distributions . in sections  [ s2 ] and
[ s3 ] , we prove the tightness of the processes in the cases @xmath1 and @xmath2 , respectively . in the case
@xmath63 we solve a more general weak convergence problem for random quadratic forms in the i.i.d .
sequence @xmath39 ; the convergence of the integrated periodograms indexed by classes of functions is only a special case
. the case @xmath2 is more involved . among others ,
entropy conditions will be needed and we only prove results on the weak convergence of the empirical spectral distribution , that is , we focus on random quadratic forms with toeplitz coefficient matrices given by the fourier coefficients @xmath64 defined in .
a glance at decomposition ( [ decom2 ] ) is enough to see that the convergence of the finite - dimensional distributions of @xmath61 is essentially determined by the weak limit behavior of the sample autocovariances @xmath65 .
for this reason , we recall a well - known result due to davis and resnick @xcite ; see also @xcite , section  13.3 .
[ le1 ] for every @xmath66 , @xmath67 where @xmath68 denotes weak convergence , the @xmath69 s are independent , @xmath70 is @xmath71 and @xmath72 are i.i.d .
@xmath73 for some @xmath74 , @xmath75 .
in particular , @xmath76    the latter result is an immediate consequence of and the continuous mapping theorem .
lemma  [ le1 ] yields the weak convergence for any finite linear combination of the sample autocovariances and autocorrelations .
it also suggests that the weak limit of the standardized process @xmath57 will be determined by the infinite series @xmath77 .
but this also means that we need to require additional assumptions on the sequence @xmath78 .
we will treat this problem in a more general context .
consider a sequence @xmath79 that is , @xmath80 satisfies the summability condition @xmath81 .
for such an @xmath80 , we define the sequences of processes @xmath82 , & \quad$y(\mathbf{a})= \displaystyle\sum_{k=1}^{\infty } a_k   y_k , $ \vspace*{1pt}\cr \displaystyle\widetilde x_n(\mathbf{a})= ( n/\log n)^{1/\alpha } \sum_{k=1}^{n-1}a_k \rho_{n,\varepsilon}(k ) , & \quad$\widetilde y(\mathbf { a})=y(\mathbf{a})/y_0 .$ } \ ] ] here , @xmath83 are independent stable random variables , as described in lemma  [ le1 ] .
the 3-series theorem @xcite implies that @xmath84 is equivalent to the a.s .
convergence of the infinite series @xmath85 in ( [ eqrto ] ) .
however , for the weak convergence of @xmath86 and @xmath87 , we need a slightly stronger assumption : @xmath88 this assumption ensures the weak convergence of the random quadratic forms in ( [ eqrto ] ) ; see the proof of theorem  [ pr1 ] below .
assumptions of this type frequently occur in the literature on infinite variance quadratic forms ( e.g. , @xcite ) .
they appear in a natural way in tail estimates for quadratic forms in i.i.d .
stable random variables ; see the .
we can now formulate our result concerning the convergence of the finite - dimensional distributions .
[ pr1 ] for any @xmath40 , @xmath89    using a cramr  wold argument , it suffices to prove the convergence of the one - dimensional distributions . from ( [ aqua ] ) and the continuous mapping theorem
, it immediately follows that for every @xmath90 , @xmath91 \longrightarrow y_m(\mathbf{a})=\sum_{k=1}^{m } a_k   y_k , \ ] ] where @xmath68 denotes weak convergence . since @xmath92 , as @xmath93 follows from the 3-series theorem .
according to @xcite , theorem 4.2 , it remains to show that @xmath94\biggr|>\epsilon\biggr)=0\end{aligned}\ ] ] for every @xmath95 and @xmath96 .
we write @xmath97 for the above probabilities . applying lemma  [ le2 ] in the and the fact that @xmath98
, we conclude that @xmath99 ( the constant on the right - hand side depends on @xmath100 . )
this proves the theorem for @xmath86 ; the convergence of @xmath87 can be shown analogously by utilizing ( [ aqua1 ] ) .
as an immediate corollary of theorem  [ pr1 ] , we obtain the following result which solves the problem of finding the limits of the finite - dimensional distributions for the integrated periodogram @xmath61 in ( [ decom2 ] ) and its self - normalized version @xmath101 in ( [ decom2a ] ) .
[ co1 ] let @xmath40 and @xmath102\dvt   \mathbf { a}(f)=(a_1(f),a_2(f),\ldots)\in\ell^{\alpha } \log\ell\ } , \ ] ] where @xmath103 is specified in @xmath104 . we then have @xmath105_{f\in\mathcal{f}}&\stackrel{\mathit{fidi}}{\longrightarrow } & 2 [ y(\mathbf{a}(f))]_{f\in\mathcal{f } } , \\
( n/\log n)^{1/\alpha}[\widetilde j_{n,\varepsilon } ( f)-a_0(f)]_{f\in\mathcal{f}}&\stackrel { \mathit{fidi}}{\longrightarrow } & 2 [ \widetilde y(\mathbf{a}(f))]_{f\in\mathcal{f } } .\end{aligned}\ ] ]    [ rem:1 ] the condition @xmath106 is , in general , not easily verified . however , if @xmath60 represents the spectral density of a stationary process @xmath86 with absolutely summable autocovariance function @xmath107 , then , up to a constant multiple , @xmath60 is represented by the fourier series of  @xmath107 , and the rate of decay of @xmath108 as @xmath109 is well known for numerous time series models . for example , if @xmath60 is the spectral density of an arma process , @xmath108 at an exponential rate and then @xmath110 is satisfied for every @xmath111 .    conditions ensuring that @xmath112 can be found in the literature on fourier series , for example , in @xcite .
theorem ( 3.10 ) on page 243 in volume i of that reference yields , for lipschitz continuous functions @xmath60 with exponent @xmath113 $ ] , that @xmath112 for @xmath114 , but not necessarily for @xmath115 .
this means , in particular , that lipschitz continuous functions do not necessarily satisfy @xmath112 for small values @xmath116 .
zygmund s theorem ( 3.13 ) , @xcite , page 243 , volume i , states that @xmath117 if @xmath60 is of bounded variation and lipschitz continuous with exponent @xmath113 $ ] such that @xmath118 , but this statement is not necessarily valid for @xmath119 .
we also note that @xmath120 for @xmath121}$ ] , @xmath122 $ ] and @xmath123 .
indeed , then @xmath124 , @xmath125 and @xmath126 diverges a.s .  by the 3-series theorem and the 01 law
hence , corollary  [ co1 ] does not apply to the important class of indicator functions when @xmath116 .
moreover , @xmath127 is not tight .
indeed , it follows from the argument above and from @xcite , theorem 6.2.1 that for some @xmath128 , for every @xmath129 , @xmath130\biggr|>k\biggr ) \\ &
\le & \mathit{const}\ , \liminf_{n\to\infty } p\bigl(n ( n \log n)^{-1/\alpha } |j_{n,\varepsilon } ( f)-a_0(f ) \gamma_{n,\varepsilon}(0)|>k\bigr ) .\end{aligned}\ ] ]      in order to derive a full weak convergence counterpart of the convergence in terms of the finite - dimensional distributions in corollary [ co1 ] , it remains to establish tightness of the corresponding family of laws .
we start , once again , in the more general context of random fields indexed by sequences in @xmath132 .
since we are dealing with the weak convergence of infinite - dimensional objects , we may expect difficulties which are due to the geometric properties of the underlying path spaces .
it is also not completely surprising that the case @xmath1 is the ` better one ' in comparison with @xmath2 ; see , for example , the results on boundedness , continuity and oscillations of @xmath0-stable processes in @xcite , chapter 10 .
note , however , that the constraint @xmath133 is harder to satisfy for smaller @xmath0 than for larger @xmath0 ; see also remark  [ rem:1 ] .    in the present case @xmath1
, we introduce the function @xmath134 where @xmath135 is chosen so large that @xmath136 is concave on @xmath137 .
note that @xmath138 and this set is a linear metric space when endowed with the metric @xmath139    assume that @xmath140 is a compact set of @xmath141 with the additional property that @xmath142 observe that @xmath140 is then also a compact subset of @xmath143 and @xmath144 is sample - continuous as a random element with values in @xmath145 , the space of continuous functions on @xmath140 equipped with the uniform topology ; see @xcite , section  10.4 .
the following is our main result on the weak convergence of the sequences @xmath86 and @xmath87 of infinite variance random quadratic forms in the case @xmath1 .
[ th1 ] assume @xmath1 .
for a compact set @xmath140 of @xmath132 satisfying @xmath146 the following weak convergence result holds in @xmath145 : @xmath147 where @xmath148 , @xmath149 , @xmath150 and @xmath151 are defined in @xmath152 and the processes @xmath150 and @xmath151 are sample - continuous .
we restrict ourselves to showing that @xmath153 . in view of theorem
[ pr1 ] , it suffices to prove the tightness of the processes @xmath148 in @xmath154 , where @xmath155 is the restriction of @xmath156 to @xmath140 .
we have , for positive @xmath100 and @xmath157 , @xmath158\\[-8pt ] & & \quad \le p\biggl(\sum_{k=1}^{n-1 } \sup_{d_{\mathcal{a}}(\mathbf { a},\mathbf{b})<\delta}|a_k - b_k| \bigl[n   \gamma_{n,|\varepsilon
= p_n(\epsilon,\delta ) .\nonumber\end{aligned}\ ] ] we want to show that @xmath159 can be made arbitrarily small for all @xmath160 , provided that @xmath157 is small .
we solve this problem in a modified form : let @xmath161 be a sequence of i.i.d .
@xmath162 random variables , independent of @xmath39 , and @xmath163 a double array of real numbers .
we then have @xmath164 by virtue of this argument , it suffices to replace the products @xmath165 in the quadratic form in ( [ eq4a ] ) with the products @xmath166 .
this means that it suffices to show that @xmath167 can be made arbitrarily small for all @xmath160 , provided that @xmath157 is small , where @xmath168 now , first apply lemma  [ le3 ] to the @xmath169 s and then condition ( [ eq2 ] ) : @xmath170    theorem  [ th1 ] provides the limit process for a very general class of random quadratic forms with infinite first moments .
the coefficient matrices of these quadratic forms are given by toeplitz matrices .
the conditions on the parameter set @xmath140 are nothing but restrictions on the infinite toeplitz matrices @xmath171 .
when specified to the particular case of fourier coefficients , as in ( [ eq45 ] ) , theorem  [ th1 ] yields the following .
[ co2 ] assume that @xmath1 and let @xmath172\dvt   \mathbf{a}(f)=(a_1(f),a_2(f),\ldots ) \in\mathcal{a}\ } , \ ] ] where @xmath140 is a compact set of @xmath132 satisfying @xmath173 and @xmath103 is specified in @xmath104 .
we then have @xmath174_{f\in\mathcal{f}}\longrightarrow 2 [ y(\mathbf{a}(f))]_{f\in\mathcal{f } } , \vspace*{3pt}\cr ( n/\log n)^{1/\alpha}[\widetilde j_{n,\varepsilon } ( f)-a_0(f)]_{f\in\mathcal{f}}\longrightarrow 2 [ \widetilde y(\mathbf{a}(f))]_{f\in\mathcal{f } } , } \ ] ] where the convergence holds in @xmath175 .
let @xmath176 be defined by @xmath177 .
we claim that @xmath178 is closed , hence compact .
indeed , if @xmath179 is such that @xmath180 converges in @xmath132 to some point @xmath181 , then ( as @xmath182 ) , the sequence @xmath183 , \ ] ] @xmath184 converges in @xmath185 $ ] to some function @xmath60 that must necessarily be in @xmath38 .
therefore , @xmath186 , and the latter set is compact .
the above argument shows that @xmath187 $ ] convergence in @xmath38 is equivalent to @xmath132 convergence in @xmath188 .
since theorem [ th1 ] implies weak convergence of the left - hand side of to its right - hand side in @xmath145 ( when each function @xmath12 is identified with @xmath189 ) , we conclude that weak convergence in also holds in  @xmath175.=1      establishing full weak convergence in the case @xmath2 is more difficult than in the case @xmath1 .
indeed , for @xmath1 , we were allowed to switch from the random variables @xmath191 to their absolute values , due to the specific geometry of the spaces @xmath192 and , in particular , @xmath141 .
the spaces @xmath192 , @xmath2 , have a much more complicated structure and , therefore , the particular geometry of these spaces will be need to be invoked in proving tightness for the random quadratic forms @xmath148 and @xmath149 .
the requirements prescribed by the geometry are usually given by entropy conditions ; see @xcite for a general treatment of random elements with values in banach spaces .
entropy conditions are typically needed when @xmath0-stable processes with @xmath2 appear ; see the discussion in @xcite , chapter 12 .    in this section ,
we only consider vectors @xmath193 of the form ( [ eq45 ] ) , that is , they are the fourier coefficients of some functions @xmath60 .
corollary  [ co1 ] determines the structures of the limit processes of the quadratic forms @xmath61 via the convergence of their finite - dimensional distributions .
it hence suffices to show the tightness in @xmath175 for suitable classes @xmath38 .
@xcite considered the special case of the one - dimensional class @xmath194 of indicator functions on @xmath13 $ ] .
we extend their approach to more general classes of functions , using an entropy condition .    for @xmath195 ,
let @xmath196 each @xmath197 defines a pseudo - metric on @xmath38 .
let @xmath198 recall that the @xmath100-covering number @xmath199 of @xmath200 is the minimal integer @xmath201 for which we can find functions @xmath202 such that @xmath203    [ thm3 ] assume that @xmath2 , define @xmath103 as in @xmath204 and let @xmath38 be a subset of @xmath187 $ ] satisfying :    @xmath106 for all @xmath205 ;    @xmath206 such that @xmath207 , \qquad \epsilon > 0 ,    k \ge0 .\end{aligned}\ ] ]    the weak convergence result @xmath208 then holds in @xmath209 .
in contrast to the finite variance case @xcite , the entropy condition is a rather strong one .
indeed , in the papers mentioned , integrability of some power of @xmath210 in a neighborhood of the origin suffices .
however , conditions such as are common in problems of continuity and boundedness for stable processes ; see @xcite , chapter 10 .
proof of theorem [ thm3 ] the convergence of the finite - dimensional distributions follows from theorem  [ pr1 ] .
we restrict ourselves to proving tightness for @xmath61 , which follows by proving that @xmath211 where @xmath212 and @xmath213 , \qquad j=1,2,\ldots .\ ] ] as in @xcite , ( 6.4 ) , one can argue that it suffices in ( [ eq56 ] ) to consider @xmath201 and @xmath160 of some specific form .
let @xmath214 be two positive integers , @xmath215 and @xmath216 , and consider numbers @xmath217 , @xmath218 with @xmath219 .
for @xmath220 sufficiently small and @xmath221 large enough , we have @xmath222 consider an array @xmath223 of positive numbers such that @xmath224 as @xmath225 for each @xmath226 . then , @xmath227    where @xmath228 \widehat\gamma_{n,\varepsilon}(j ) \biggr| >
2^{-(l+1)}\epsilon_k \biggr ) .\end{aligned}\ ] ] by virtue of lemma  [ le2 ] , we have , for all @xmath195 , @xmath229 \widehat\gamma_{n,\varepsilon}(j ) \biggr|>2^{-(l+1)}\epsilon_k \biggr ) \le \mathit{const}\ , b_{k , l } , \end{aligned}\ ] ] where @xmath230 b_{k , l}= \epsilon_k^{-\alpha}2^{\alpha l } \sum_{j=2^k}^{2^{k+1}-1}|a_j(f)-a_j(g)|^{\alpha } \bigl[1 + \log^+\bigl(1/|a_j(f)-a_j(g)|\bigr)\bigr ] .\end{aligned}\ ] ] assuming that @xmath231 , we have @xmath232 \\ % a similar argument yields % b_{k , l } & \le & \mathit{const}\ , \epsilon_k^{-\alpha } 2^{\alpha l } \epsilon_{k , l-1}^{\alpha } 2^{-k(\alpha-1 ) } [ 1+k \log^+\epsilon_{k , l-1}^{-1 } ] .\end{aligned}\ ] ] hence , we are left to consider @xmath233\nonumber \\ & & \quad= \sum_{k = a}^b 2^{-k(\alpha-1-\alpha\theta ) } \sum_{l=1}^{\infty } n(\epsilon_{k , l},\mathcal{f},\rho_k)\epsilon_{k , l-1}^{\alpha } [ 1+k \log^+\epsilon_{k , l-1}^{-1 } ] 2^{\alpha l } \\ & & \quad\le \mathit{const}\ , \sum_{k = a}^b 2^{-k(\alpha-1-\alpha\theta ) } \sum _ { l=1}^\infty\biggl[1 + \biggl ( \frac{2^k}{\epsilon_{k , l } } \biggr)^\beta\biggr]\epsilon_{k , l-1}^\alpha[1 + \log^+\epsilon _ { k , l-1}^{-1 } ] 2^{\alpha l } .\nonumber
% 2^{-k(\alpha-1-\alpha\theta ) } % [ 1+k \log^+\epsilon_{k , l-1}^{-1 } ] 2^{\alpha l}\end{aligned}\ ] ] assume that @xmath220 is so small that ( [ neu ] ) holds .
define the numbers @xmath234 , @xmath235 with @xmath236 such that @xmath237 and @xmath238 . for these parameter choices ,
it is not difficult to see that ( [ nanuchen ] ) converges to zero by first letting @xmath239 ( i.e. , @xmath240 ) and then @xmath93 ( i.e. , @xmath241 ) .
this proves ( [ eq56 ] ) , hence the tightness of the processes considered in @xmath175 .    in what follows ,
we give examples of function spaces @xmath38 satisfying condition ( ii ) of theorem  [ thm3 ] .
consider a space of indexed functions @xmath242 that are defined on @xmath243 $ ] such that @xmath244 is a compact metric space , the mapping @xmath245 is hlder continuous with exponent @xmath246 and constant @xmath129 , that is , @xmath247 and the number of balls ( in metric @xmath248 ) of radius at most @xmath249 necessary to cover @xmath34 is of the order @xmath250 for some @xmath251 .
then , @xmath252 satisfies @xmath253 with @xmath254 .
indeed , let @xmath255 .
we can find @xmath256 balls of radius at most @xmath257 covering @xmath34 . call them @xmath258 , with centers @xmath259 .
now , given @xmath260 , we have @xmath261 for some @xmath262 and @xmath263 the desired bound now follows since @xmath264    consider a vapnik  ervonenkis ( vc ) class @xmath265 of functions defined on @xmath13 $ ] with vc index @xmath266 ; see @xcite , section 2.6.2 for more information on vc classes of functions .
given @xmath267 and @xmath268 , we can find @xmath269 balls of radius at most @xmath270 that cover @xmath265 in the norm @xmath271 ; see , for example , @xcite , theorem 2.6.7 .
therefore , there exist @xmath272 such that for any @xmath273 , @xmath274 we then have @xmath275 it follows that @xmath276
it is the aim of this section to show that the results for the case of an i.i.d .
sequence @xmath39 translate to the linear process case
. the following decomposition will be crucial : @xmath277 this decomposition is analogous to the decomposition @xmath278 of the spectral density @xmath22 of a linear process .
we will show that the normalized integrated remainder term @xmath279 is negligible , uniformly over the class of functions @xmath280 , in comparison to the normalized main part @xmath281 which can be treated by the methods of the previous section .
note that , for a given sequence of coefficients @xmath282 , the functions @xmath283 constitute just another class of functions on @xmath13 $ ] , @xmath284 say , and therefore we will study the process @xmath285 , @xmath286 , for suitable classes @xmath287 .
[ remainder ] let @xmath288 be the remainder term appearing in the decomposition of the periodogram @xmath289 .
suppose that the linear filter @xmath290 of the process @xmath291 satisfies @xmath292 and @xmath38 is a collection of real - valued functions defined on @xmath243 $ ] such that @xmath293 @xmath294 .
we then have @xmath295    from @xcite , proposition 5.1 , substituting @xmath296 for @xmath297 , we have the following decomposition for @xmath288 : @xmath298 where @xmath299 is the transfer function as defined before and @xmath300 we first show that @xmath301 note that @xmath302 & \le & % |u_{nj}(x)|^2 dx\\ % & = & \mathit{const}\ , \biggl ( \sum_{j=-\infty}^{-1 } + \sum_{j=1}^\infty\biggr )    the convergence in ( [ izmir ] ) will follow if we can show that the suprema over @xmath303 of the two infinite sums in the last expression are bounded in probability as @xmath304 .
we will prove this for the second sum ; the first one can be handled analogously .
we have , by definition of the terms @xmath305 , the cauchy
schwarz inequality and the fact that , by assumption , @xmath306 , @xmath307 & & \quad\le \sup_{f\in{\mathcal{f}}}\sum_{j=1}^n |\psi_j| \int_0^\curpi|f(x)| \biggl|\sum_{t=1-j}^0 \varepsilon_t \mathrm{e}^{-\mathrm{i}xt } - \sum_{t = n - j+1}^n \varepsilon_t \mathrm{e}^{-\mathrm{i}xt}\biggr|^2\ , \mathrm{d}x \\[-2pt ] & & \qquad { } + \sup_{f\in{\mathcal{f}}}\sum_{j = n+1}^\infty|\psi_j| \int _ 0^\curpi    \sum _ { t=1}^n \varepsilon_t \mathrm{e}^{-\mathrm{i}xt}\biggr|^2\ , \mathrm{d}x \\[-2pt ] & & \quad\le c   [ i_1(n ) + i_2(n ) + i_3(n ) + i_4(n ) ] , \end{aligned}\ ] ] where @xmath308 i_2(n)&= & \sum_{j=1}^n |\psi_j| \biggl(\int_0^\curpi\biggl|\sum _ { t = n - j+1}^n \varepsilon_t \mathrm{e}^{-\mathrm{i}xt}\biggr|^4\ , \mathrm{d}x \biggr)^{1/2 } , \\[-2pt ] i_3(n)&= & \sum_{j = n+1}^\infty|\psi_j| \biggl(\int_0^\curpi \biggl|\sum_{t=1-j}^{n - j } \varepsilon_t \mathrm{e}^{-\mathrm{i}xt}\biggr|^4 \,\mathrm{d}x \biggr)^{1/2 } , \\[-2pt ] i_4(n)&= & \sum_{j = n+1}^\infty|\psi_j| \biggl(\int_0^\curpi\biggl| \sum_{t=1}^n \varepsilon_t \mathrm{e}^{-\mathrm{i}xt}\biggr|^4 \,\mathrm{d}x \biggr)^{1/2 } .\end{aligned}\ ] ] it remains to show that each sequence @xmath309 , is tight .
now , @xmath310 let @xmath267 .
choose @xmath311 so large that the following holds , for @xmath312 : @xmath313 write @xmath314 then , for @xmath315 , we have , for @xmath157 chosen as above , @xmath316^{1/2 } \\ & & \quad \le \frac{c}{k } \sum_{j=1}^n |\psi_j|\biggl [ \biggl(\sum_{m=1}^j \bigl ( m^{1/\alpha}(1 + \log m)^{{1/\alpha } + \delta } \bigr)^{4 - \alpha } \biggr)^{1/2 } \\ & & { } \qquad + \biggl ( \biggl[\sum_{m=1}^j \bigl ( m^{1/\alpha } ( 1 + \log m)^{{1/\alpha } + \delta } \bigr)^{2 - \alpha } \biggr]^2 \biggr)^{1/2}\biggr ] \\ % & \le&\frac{c}{k } \sum_{j=1}^n |\psi_j| [ ( \sum_{m=1}^j m^ { % + \sum_{m=1}^j m^{\frac{2}{\alpha}-1}(1 + \log m)^{(2-\alpha)(\frac{1 } { % & \le & \frac{c}{k } \sum_{j=1}^n |\psi_j| j^{2/\alpha } [ ( 1 + \log % j)^{\frac{1}{2}(4-\alpha)(\frac{1}{\alpha } + \delta ) } % + ( 1 + \log j)^{(2-\alpha)(\frac{1}{\alpha } + \delta)}]\\ % & \le & \frac{c}{k } \sum_{j=1}^n |\psi_j| j^{2/\alpha } ( 1 + \log j)^ { & & \quad \le \frac{c}{k } \sum_{j=1}^\infty|\psi_j|   j^{2/\alpha }   ( 1 + \log j)^{{(4-\alpha)/(2\alpha ) } + \tau } .\end{aligned}\ ] ] by virtue of , the last expression can be made smaller than @xmath317 by choosing @xmath318 large enough , which proves the tightness of @xmath319 .
similar arguments show that @xmath320 , @xmath321 are tight sequences as well .
the convergence in ( [ izmir ] ) follows . by the decomposition ( [ vienna ] )
, the proof will be finished if we can also establish that @xmath322 we have , by the cauchy
schwarz inequality and the identity @xmath323 , @xmath324 similar arguments as for ( [ izmir ] ) ensure the tightness of the sequence @xmath325 .
the tightness of the term @xmath326 follows from @xcite theorem 2.1 ( for @xmath327 ) and proposition 3.1 ( for @xmath328 ) .
thus , we conclude that ( [ prague ] ) holds , and lemma [ remainder ] is proved .
a referee kindly pointed out that lemma  [ remainder ] remains valid under the following condition , which is weaker than : assume that there exists a sequence @xmath329 of positive numbers such that @xmath330 condition follows by taking @xmath331 for some positive @xmath157 .    by ( [ munich ] ) , we may write , for each @xmath60 , @xmath332 where @xmath333 stands for @xmath334 . combining this decomposition with lemma [ remainder ] , we can now state the following analogs to corollary [ co2 ] and theorem [ thm3 ] .
assume that @xmath1 or @xmath2 and let @xmath335 be as defined as in corollary [ co2 ] or theorem  [ thm3 ] , respectively .
suppose that the set @xmath336
\rightarrow\mathbb{r}\dvt f|\psi|^2 \in\mathcal{f}\ } = \mathcal{f}_\psi$ ] satisfies @xmath337 and holds for some @xmath338 .
we then have @xmath339_{f\in\mathcal{f}_\psi}\longrightarrow 2 [ y(\mathbf{a}(f|\psi|^2))]_{f\in\mathcal{f}_\psi } , \vspace*{3pt}\cr ( n/\log n)^{1/\alpha}[\widetilde j_{n , x}(f)-a_0(f|\psi    2 [ \widetilde y(\mathbf{a}(f|\psi|^2))]_{f\in\mathcal{f}_\psi } , } \vspace*{2pt}\ ] ] where the convergence holds in @xmath340 .
for an array @xmath341 of real numbers , define the quadratic forms @xmath342 and @xmath343 the following lemma is a consequence of @xcite , theorem  3.1 ; see also @xcite .
[ le2 ] for @xmath40 , there exists a positive constant @xmath344 such that for all @xmath345 , @xmath346    now , let @xmath161 be a sequence of i.i.d .
@xmath162 random variables , independent of @xmath39 , and let @xmath347 be as above .
the following lemma is a consequence of lemma  [ le2].=1    [ le3 ] for @xmath1 , there exists a positive constant @xmath348 such that for all @xmath345 ,    @xmath349    apply lemma  [ le2 ] to @xmath350 , conditionally on @xmath351 : @xmath352 because @xmath1 , we also have , for @xmath345 , @xmath353 and combining ( [ eq8 ] ) and ( [ eq9 ] ) , we thus obtain ( [ eq18 ] ) .
thomas mikosch s research was supported in part by the danish research council ( fnu ) grant 272 - 06 - 0442 .
gennady samorodnitsky s research was supported in part by the aro grant w911nf-07 - 1 - 0078 at cornell university and a villum kann rasmussen visiting professor grant at the university of copenhagen .
the paper was written when sami umut can and gennady samorodnitsky visited the department of mathematical sciences at the university of copenhagen in 2008/2009 .
they take pleasure in thanking the department for hospitality and financial support .
the authors also thank the anonymous referees for suggestions that helped to improve the paper .
can , s.u . , mikosch , t. and samorodnitsky , g. ( 2009 ) .
weak convergence of the function - indexed integrated periodogram for infinite variuance processes . extended technical report .
available at http://www.math.ku.dk/\textasciitilde mikosch[www.math.ku.dk/~mikosch ] .
crovella , m. and bestavros , a. ( 1996 ) .
self - similarity in world wide web traffic : evidence and possible causes . in _ proceedings of the 1996 acm sigmetrics international conference on measurement and modeling of computer systems _
* 24 * 160169 .
dahlhaus , r. and polonik , w. ( 2002 ) .
empirical processes and nonparametric maximum likelihood estimation for time series . in _ empirical process techniques for dependent data
dehling , t. mikosch and m. srensen , eds . )
boston : birkhuser .
davis , r.a . and mikosch , t. ( 2009 ) .
extreme value theory for garch processes . in
_ the handbook of financial time series _ ( t.g .
andersen , r.a .
kreiss and t. mikosch , eds . ) 187200 .
heidelberg : springer .
davis , r.a . and mikosch , t. ( 2009 ) . extremes of stochastic volatility models . in _ the handbook of financial time series _
andersen , r.a .
davis , j .-
kreiss and t. mikosch , eds . ) 355364 .
heidelberg : springer .
klppelberg , c. and mikosch , t. ( 1996 ) . self - normalized and randomly centred spectral estimates . in _ proceedings of the athens international conference on applied probability and time series , vol
. 2 : time series _ ( c.c .
heyde , yu.v .
prokhorov , r. pyke and s.t .
rachev , eds . ) 259271 .
berlin : springer .
mikosch , t. ( 1998 ) .
periodogram estimates from heavy - tailed data . in _ a practical guide to heavy tails : statistical techniques for analysing heavy - tailed distributions _ ( r.a .
adler , r. feldman and m.s .
taqqu , eds . ) 241258 .
boston : birkhuser .
mikosch , t. ( 2003 ) . modelling dependence and tails of financial time series . in _
extreme values in finance , telecommunications and the environment _ ( b. finkenstdt and h. rootzn , eds . ) 185286 .
london : chapman & hall .
willinger , w. , taqqu , m.s . ,
sherman , r. and wilson , d. ( 1995 ) .
self - similarity through high variability : statistical analysis of ethernet lan traffic at the source level .
proceedings of the acm / sigcomm95 , cambridge , ma . _ computer communications review _ * 25 * 100113 .