Search

We investigate a diffusion process in heterogeneous media where particles stochastically reset to their initial positions at a constant rate. The heterogeneous media is modeled using a spatial-dependent diffusion coefficient with a power-law dependence on particles' positions. We use the Green function approach to obtain exact solutions for the probability distribution of particles' positions and the mean square displacement. These results are further compared and agree with numerical simulations of a Langevin equation. We also study the first-passage time problem associated with this diffusion process and obtain an exact expression for the mean first-passage time. Our findings show that this system exhibits non-Gaussian distributions, transient anomalous diffusion (sub- or superdiffusion) and stationary states that simultaneously depend on the media heterogeneity and the resetting rate. We further demonstrate that the media heterogeneity non-trivially affect the mean first-passage time, yielding an optimal resetting rate for which this quantity displays a minimum., Comment: 8 pages, 3 figures; accepted for publication Physica A

We establish, in general terms, the conditions to be satisfied by a time-fractional approach formulation of the Poisson-Nernst-Planck model in order to guarantee that the total current across the sample be solenoidal, as required by the Maxwell equation. Only in this case the electric impedance of a cell can be determined as the ratio between the applied difference of potential and the current across the cell. We show that in the case of anomalous diffusion, the model predicts for the electric impedance of the cell a constant phase element behaviour in the low frequency region. In the parametric curve of the reactance versus the resistance, the slope coincides with the order of the fractional time derivative.

The various types of generalized Cattaneo, called also telegrapher's equation, are studied. We find conditions under which solutions of the equations considered so far can be recognized as probability distributions, \textit{i.e.} are normalizable and non-negative on their domains. Analysis of the relevant mean squared displacements enables us to classify diffusion processes described by such obtained solutions and to identify them with either ordinary or anomalous super- or subdiffusion. To complete our study we analyse derivations of just considered examples the generalized Cattaneo equations using the continuous time random walk and the persistent random walk approaches.

Recent experimental findings on anomalous diffusion have demanded novel models that combine annealed (temporal) and quenched (spatial or static) disorder mechanisms. The comb-model is a simplified description of diffusion on percolation clusters, where the comb-like structure mimics quenched disorder mechanisms and yields a subdiffusive regime. Here we extend the comb-model to simultaneously account for quenched and annealed disorder mechanisms. To do so, we replace usual derivatives in the comb diffusion equation by different fractional time-derivative operators and the conventional comb-like structure by a generalized fractal structure. Our hybrid comb-models thus represent a diffusion where different comb-like structures describe different quenched disorder mechanisms, and the fractional operators account for various annealed disorders mechanisms. We find exact solutions for the diffusion propagator and mean square displacement in terms of different memory kernels used for defining the fractional operators. Among other findings, we show that these models describe crossovers from subdiffusion to Brownian or confined diffusions, situations emerging in empirical results. These results reveal the critical role of interactions between geometrical restrictions and memory effects on modeling anomalous diffusion., Comment: 24 pages, 4 figures; accepted for publication in Physical Review E

The spatial dynamics of criminal activities has been recently studied through statistical physics methods; however, models and results have been focused on local scales (city level) and much less is known about these patterns at larger scales such as at a country level. Here we report on a characterization of the spatial dynamics of the homicide crimes along the Brazilian territory using data from all cities (~5000) in a period of more than thirty years. Our results show that the spatial correlation function in the per capita homicides decays exponentially with the distance between cities and that the characteristic correlation length displays an acute increasing trend in the latest years. We also investigate the formation of spatial clusters of cities via a percolation-like analysis, where clustering of cities and a phase transition-like behavior describing the size of the largest cluster as a function of a homicide threshold are observed. This transition-like behavior presents evolutive features characterized by an increasing in the homicide threshold (where the transitions occur) and by a decreasing in the transition magnitudes (length of the jumps in the cluster size). We believe that our work sheds new lights on the spatial patterns of criminal activities at large scales, which may contribute for better political decisions and resources allocation as well as opens new possibilities for modeling criminal activities by setting up fundamental empirical patterns at large scales., Comment: Accepted for publication in EPL

We report on an extensive characterization of the cracking noise produced by charcoal samples when dampened with ethanol. We argue that the evaporation of ethanol causes transient and irregularly distributed internal stresses that promote the fragmentation of the samples and mimic some situations found in mining processes. The results show that, in general, the most fundamental seismic laws ruling earthquakes (Gutenberg-Richter law, unified scaling law for the recurrence times, Omori's law, productivity law and Bath's law) hold under the conditions of the experiment. Some discrepancies were also identified (a smaller exponent in Gutenberg-Richter law, a stationary behavior in the aftershock rates for long times and a double power-law relationship in productivity law) and related to the different loading condition. Our results thus corroborate to elucidate the parallel between seismic laws and fracture experiments caused by a more complex loading condition that also occurs in natural and induced seismicity (such as long-term fluid injection and gas-rock outbursts in mining processes)., Comment: Accepted for publication in PRL

The comb model is a simplified description for anomalous diffusion under geometric constraints. It represents particles spreading out in a two-dimensional space where the motions in the x-direction are allowed only when the y coordinate of the particle is zero. Here, we propose an extension for the comb model via Langevin-like equations driven by fractional Gaussian noises (long-range correlated). By carrying out computer simulations, we show that the correlations in the y-direction affect the diffusive behavior in the x-direction in a non-trivial fashion, resulting in a quite rich diffusive scenario characterized by usual, superdiffusive or subdiffusive scaling of second moment in the x-direction. We further show that the long-range correlations affect the probability distribution of the particle positions in the x-direction, making their tails longer when noise in the y-direction is persistent and shorter for anti-persistent noise. Our model thus combines and allows the study/analysis of the interplay between different mechanisms of anomalous diffusion (geometric constraints and long-range correlations) and may find direct applications for describing diffusion in complex systems such as living cells., Comment: Accepted for publication in New Journal of Physics

We report on a statistical analysis of the lightning activity rates in all Brazilian cities. We find out that the average of lightning activity rates exhibit a dependence on the latitude of the cities, displaying one peak around the Tropic of Capricorn and another one just before the Equator. We verify that the standard deviation of these rates is almost a constant function of the latitude and that the distribution of the fluctuations surrounding the average tendency is quite well described by a Gumbel distribution, which thus connects these rates to extreme processes. We also investigate the behavior of the lightning activity rates versus the longitude of the cities. For this case, the average rates exhibit an approximate plateau for a wide range of longitude values, the standard deviation is an approximate constant function of longitude, and the fluctuations are described by a Laplace distribution. We further characterize the spatial correlation of the lightning activity rates between pairs of cities, where our results show that the spatial correlation function decays very slowly with the distance between the cities and that for intermediate distances the correlation exhibits an approximate logarithmic decay. Finally, we propose to model this last behavior within the framework of the Edwards-Wilkinson equation., Comment: Accepted for publication in EPL

We report on a statistical analysis of the engagement in the electoral processes of all Brazilian cities by considering the number of party memberships and the number of candidates for mayor and councillor. By investigating the relationships between the number of party members and the population of voters, we have found that the functional form of these relationships are well described by sub-linear power laws (allometric scaling) surrounded by a multiplicative log-normal noise. We have observed that this pattern is quite similar to those previously-reported for the relationships between the number candidates (mayor and councillor) and population of voters [EPL 96, 48001 (2011)], suggesting that similar universal laws may be ruling the engagement in the electoral processes. We also note that the power law exponents display a clear hierarchy, where the more influential is the political position the smaller is the value of the exponent. We have also investigated the probability distributions of the number of candidates (mayor and councilor), party memberships and voters. The results indicate that the most influential positions are characterized by distributions with very short-tails, while less influential positions display an intermediate power law decay before showing an exponential-like cutoff. We discuss that, in addition to the political power of the position, limitations in the number of available seats can also be connected with this changing of behavior. We further believe that our empirical findings point out to an underrepresentation effect, where the larger city is, the larger are the obstacles for more individuals to become directly engaged in the electoral process., Comment: Accepted for publication in PRE

We report on the dynamical behavior of defects of strength s = +/- 1/2 in a lyotropic liquid crystal during the annihilation process. By following their positions using time resolved polarizing microscopy technique, we present statistically significant evidence that the relative velocity between defect pairs is Gaussian distributed, anti-persistent and long-range correlated. We further show that simulations of the Lebwohl-Lasher model reproduce quite well our experimental findings., Comment: Accepted for publication in PRE as Brief Report

The complexity of chess matches has attracted broad interest since its invention. This complexity and the availability of large number of recorded matches make chess an ideal model systems for the study of population-level learning of a complex system. We systematically investigate the move-by-move dynamics of the white player's advantage from over seventy thousand high level chess matches spanning over 150 years. We find that the average advantage of the white player is positive and that it has been increasing over time. Currently, the average advantage of the white player is ~0.17 pawns but it is exponentially approaching a value of 0.23 pawns with a characteristic time scale of 67 years. We also study the diffusion of the move dependence of the white player's advantage and find that it is non-Gaussian, has long-ranged anti-correlations and that after an initial period with no diffusion it becomes super-diffusive. We find that the duration of the non-diffusive period, corresponding to the opening stage of a match, is increasing in length and exponentially approaching a value of 15.6 moves with a characteristic time scale of 130 years. We interpret these two trends as a resulting from learning of the features of the game. Additionally, we find that the exponent {\alpha} characterizing the super-diffusive regime is increasing toward a value of 1.9, close to the ballistic regime. We suggest that this trend is due to the increased broadening of the range of abilities of chess players participating in major tournaments., Comment: Accepted for publication in PLoS ONE

Complexity measures are essential to understand complex systems and there are numerous definitions to analyze one-dimensional data. However, extensions of these approaches to two or higher-dimensional data, such as images, are much less common. Here, we reduce this gap by applying the ideas of the permutation entropy combined with a relative entropic index. We build up a numerical procedure that can be easily implemented to evaluate the complexity of two or higher-dimensional patterns. We work out this method in different scenarios where numerical experiments and empirical data were taken into account. Specifically, we have applied the method to i) fractal landscapes generated numerically where we compare our measures with the Hurst exponent; ii) liquid crystal textures where nematic-isotropic-nematic phase transitions were properly identified; iii) 12 characteristic textures of liquid crystals where the different values show that the method can distinguish different phases; iv) and Ising surfaces where our method identified the critical temperature and also proved to be stable., Comment: Accepted for publication in PLoS One

Nowadays we are often faced with huge databases resulting from the rapid growth of data storage technologies. This is particularly true when dealing with music databases. In this context, it is essential to have techniques and tools able to discriminate properties from these massive sets. In this work, we report on a statistical analysis of more than ten thousand songs aiming to obtain a complexity hierarchy. Our approach is based on the estimation of the permutation entropy combined with an intensive complexity measure, building up the complexity-entropy causality plane. The results obtained indicate that this representation space is very promising to discriminate songs as well as to allow a relative quantitative comparison among songs. Additionally, we believe that the here-reported method may be applied in practical situations since it is simple, robust and has a fast numerical implementation., Comment: Accepted for publication in Physica A

We address this work to investigate some statistical properties of symbolic sequences generated by a numerical procedure in which the symbols are repeated following a power law probability density. In this analysis, we consider that the sum of n symbols represents the position of a particle in erratic movement. This approach revealed a rich diffusive scenario characterized by non-Gaussian distributions and, depending on the power law exponent and also on the procedure used to build the walker, we may have superdiffusion, subdiffusion or usual diffusion. Additionally, we use the continuous-time random walk framework to compare with the numerical data, finding a good agreement. Because of its simplicity and flexibility, this model can be a candidate to describe real systems governed by power laws probabilities densities., Comment: Accepted for publication in Physica Scripta

We investigate the dynamics of many interacting bubbles in boiling water by using a laser scattering experiment. Specifically, we analyze the temporal variations of a laser intensity signal which passed through a sample of boiling water. Our empirical results indicate that the return interval distribution of the laser signal does not follow an exponential distribution; contrariwise, a heavy-tailed distribution has been found. Additionally, we compare the experimental results with those obtained from a minimalist phenomenological model, finding a good agreement., Comment: Accepted for publication in Chaos, Solitons & Fractals

We report a statistical analysis about people agglomeration soundscape. Specifically, we investigate the normalized sound amplitudes and intensities that emerge from people collective meetings. Our findings support the existence of nontrivial dynamics characterized by heavy tail distributions in the sound amplitudes, long-range correlations in the sound intensity and non-exponential distributions in the return interval distributions. Additionally, motivated by the time-dependent behavior present in the volatility/variance series, we compare the observational data with those obtained from a minimalist autoregressive stochastic model, a GARCH process, finding a good agreement., Comment: To appear in New Journal of Physics

We report a statistical analysis over more than eight thousand songs. Specifically, we investigate the probability distribution of the normalized sound amplitudes. Our findings seems to suggest a universal form of distribution which presents a good agreement with a one-parameter stretched Gaussian. We also argue that this parameter can give information on music complexity, and consequently it goes towards classifying songs as well as music genres. Additionally, we present statistical evidences that correlation aspects of the songs are directly related with the non-Gaussian nature of their sound amplitude distributions., Comment: Accepted for publication as a Brief Report in Physical Review E

We address this work to investigate symbolic sequences with long-range correlations by using computational simulation. We analyze sequences with two, three and four symbols that could be repeated $l$ times, with the probability distribution $p(l)\propto 1/ l^{\mu}$. For these sequences, we verified that the usual entropy increases more slowly when the symbols are correlated and the Tsallis entropy exhibits, for a suitable choice of $q$, a linear behavior. We also study the chain as a random walk-like process and observe a nonusual diffusive behavior depending on the values of the parameter $\mu$., Comment: Published in the Brazilian Journal of Physics

In this paper a composed A+B magnetic system, with spins J_A=2 and J_B=3/2, is considered within the mean-field approximation, in the framework of Tsallis nonextensive statistics. Our motivation is twofold: (1) to approach the existing experimental data of manganese oxides (manganites), where Mn^{3+} and Mn^{4+} form two magnetic sublattices, and (2) to investigate the structure of nonextensive density matrices of composed systems. By imposing that thermodynamic quantities, such as the magnetization of sublattices A and B, must be invariant weather the calculation is taken over the total Hilbert space or over partial subspaces, we found that the expression for the nonextensive entropy must be adapted. Our argument is supported by calculation of sublattices magnetization M_A and M_B, internal energy, U_A and U_B, and magnetic specific heat, CA and CB. It is shown that only with the modified entropy the two methods of calculation agree to each other. Internal energy and magnetization are additive, but no clear relationship was found between S_A, S_B and the total entropy S_{A+B} for q \neq 1. It is shown that the reason for the failure of the standard way of calculation is the assumption of statistical independence between the two subsystems, which however does not affect the density matrix in the full Hilbert space., Comment: My surname is Navarro and not Navarrro, it is why I replace this paper. 11 pages, 4 figures

We analyse a bidimensional nonlinear Fokker-Planck equation by considering an anisotropic case, whose diffusion coefficients are $D_x \propto |x|^{-\theta}$ and $D_y \propto |y|^{-\gamma}$ with $\theta, \gamma \in {\cal{R}}$. In this context, we also investigate two situations with the drift force $\vec{F}(\vec{r},t)=(-k_{x}x, -k_y y)$. The first one is characterized by $k_x/k_y=(2+\gamma)/(2+\theta)$ and the second is given by $k_{x}=k$ and $k_{y}=0$. In these cases, we can verify an anomalous behavior induced in different directions by the drift force applied. The found results are exact and exhibit, in terms of the $q$-exponentials, functions which emerge from the Tsallis formalism. The generalization for the $D$-dimensional case is discussed., Comment: 6 pages, tex file

In this work we present the logarithmic diffusion equation as a limit case when the index that characterizes a nonlinear Fokker-Planck equation, in its diffusive term, goes to zero. A linear drift and a source term are considered in this equation. Its solution has a lorentzian form, consequently this equation characterizes a super diffusion like a L\'evy kind. In addition is obtained an equation that unifies the porous media and the logarithmic diffusion equations, including a generalized diffusion equation in fractal dimension. This unification is performed in the nonextensive thermostatistics context and increases the possibilities about the description of anomalous diffusive processes., Comment: 5 pages. To appear in Phys. Rev. E

Recently, fractional derivatives have been employed to analyze various systems in engineering, physics, finance and hidrology. For instance, they have been used to investigate anomalous diffusion processes which are present in different physical systems like: amorphous semicondutors, polymers, composite heterogeneous films and porous media. They have also been used to calculate the heat load intensity change in blast furnace walls, to solve problems of control theory \ and dynamic problems of linear and nonlinear hereditary mechanics of solids. In this work, we investigate the scaling properties related to the nonlinear fractional diffusion equations and indicate the possibilities to the applications of these equations to simulate the water transport in unsaturated soils. Usually, the water transport in soils with anomalous diffusion, the dependence of concentration on time and distance may be expressed in term of a single variable given by $\lambda _{q}=x/t^{q}.$ In particular, for $q=1/2$ the systems obey Fick's law and Richards' equation for water transport. We show that a generalization of Richards' equation via fractional approach can incorporate the above property., Comment: 9 pages

We obtain new exact classes of solutions for the nonlinear fractional Fokker-Planck-like equation partial_t rho = partial_x{D(x) partial^{mu -1}_x rho^{nu} - F(x) rho} by considering a diffusion coefficient D = D|x|^{-theta} (theta in R and D>0) and a drift force F = -k_1 x + k-bar_{gamma} x|x|^{gamma-1} (k_1, k-bar_{gamma}, gamma in R). Connection with nonextensive statistical mechanics based on Tsallis entropy is also discussed., Comment: latex, 5 pages. Submitted to Physica A

We analyzed the magnetic behavior of a S-spin system within framework of the Tsallis nonextensive statistics, employing the normalized approach. Unusual properties on magnetization, entropy and susceptibility emerge, as a consequence of nonextensivity. We further show that the nonextensive approach can be relevant to the field of manganites, materials which exhibit long-range interactions and fractality, two basic ingredients for nonextensivity. Our results are in qualitative agreement to experimental data in La$_{0.67}$Ca$_{0.33}$MnO$_3$ and Pr$_{0.5}$Ca$_{0.5}$Mn$_{0.95}$Ga$_{0.05}$O$_3$ manganites., Comment: 5 pages and 6 figures. Submitted to Phys. Rev. B

An $N$-dimensional nonlinear Fokker-Planck equation is investigated here by considering the time dependence of the coefficients, where drift-controlled and source terms are present. We exhibit the exact solution based on the generalized gaussian function related to the Tsallis statistics. Furthermore, we show that a rich class of diffusive processes, including normal and anomalous ones, can be obtained by changing the time dependence of the coefficients., Comment: Latex, 4 pages. To appear in Phys. Rev. E

In this work we incorporate, in a unified way, two anomalous behaviors, the power law and stretched exponential ones, by considering the radial dependence of the $N$-dimensional nonlinear diffusion equation $\partial\rho /\partial{t}={\bf \nabla} \cdot (K{\bf \nabla} \rho^{\nu})-{\bf \nabla}\cdot(\mu{\bf F} \rho)-\alpha \rho ,$ where $K=D r^{-\theta}$, $\nu$, $\theta$, $\mu$ and $D$ are real parameters and $\alpha$ is a time-dependent source. This equation unifies the O'Shaugnessy-Procaccia anomalous diffusion equation on fractals ($\nu =1$) and the spherical anomalous diffusion for porous media ($\theta=0$). An exact spherical symmetric solution of this nonlinear Fokker-Planck equation is obtained, leading to a large class of anomalous behaviors. Stationary solutions for this Fokker-Planck-like equation are also discussed by introducing an effective potential., Comment: Latex, 6 pages. To appear in Phys. Rev. E

In the nonextensive Tsallis scenario, Page's conjecture for the average entropy of a subsystem[Phys. Rev. Lett. {\bf 71}, 1291(1993)] as well as its demonstration are generalized, i.e., when a pure quantum system, whose Hilbert space dimension is $mn$, is considered, the average Tsallis entropy of an $m$-dimensional subsystem is obtained. This demonstration is expected to be useful to study systems where the usual entropy does not give satisfactory results., Comment: Revtex, 6 pages, 2 figures. To appear in Phys. Rev. E

The validity of (1-q) expansion and factorization approximations are analysed in the framework of Tsallis statistics. We employ exact expressions for classical independent systems (harmonic oscillators) by considering the unnormalized and normalized constrainsts. We show that these approxiamtions can not be accurate in the analysis of systems with many degrees of freedom., Comment: Latex, 6 pages, 2 figures

Usually, the study of city population distribution has been reduced to power laws. In such analysis, a common practice is to consider cities with more than one hundred thousand inhabitants. Here, we argue that the distribution of cities for all ranges of populations can be well described by using a $q$-exponential distribution. This function, which reproduces the Zipf-Mandelbrot law, is related to the generalized nonextensive statistical mechanics and satisfies an anomalous decay equation., Comment: Revtex, 4 pages. To appears in Phys. Rev. E

We found, from the analysis of $M$ vs. $T$ curves of some manganese oxides (manganites), that these systems do not follow the traditional Maxwell-Boltzmann statistics, but the Tsallis statistics, within the \QTR{em}{normalized} formalism. Curves were calculated within the mean field approximation, for various ferromagnetic samples and the results were compared to measurements of our own and to various other authors published data, chosen at random from the literature. The agreement between the experimental data and calculated $M_{q}$ vs. $T^{\ast}$ curve, where $T^{\ast}$ is an effective temperature, is excellent for all the compounds. The entropic parameter, $q$, correlates in a simple way with the experimental value of $T_{c}$, irrespect the chemical composition of the compounds, heat treatment or other details on sample preparation. Examples include $q<1$ (superextensivity), $q=1$ (extensivity) and $q>1$ (subextensivity) cases., Comment: 12 pages, 3 figures

We investigate the escape behavior of systems governed by the one-dimensional nonlinear diffusion equation $\partial_t \rho = \partial_x[\partial_x U\rho] + D\partial^2_x \rho^\nu$, where the potential of the drift, $U(x)$, presents a double-well and $D, \nu$ are real parameters. For systems close to the steady state we obtain an analytical expression of the mean first passage time, yielding a generalization of Arrhenius law. Analytical predictions are in very good agreement with numerical experiments performed through integration of the associated Ito-Langevin equation. For $\nu\neq 1$ important anomalies are detected in comparison to the standard Brownian case. These results are compared to those obtained numerically for initial conditions far from the steady state., Comment: to appear in PRE

The nonlinear diffusion equation $\frac{\partial \rho}{\partial t}=D \tilde{\Delta} \rho^\nu$ is analyzed here, where $\tilde{\Delta}\equiv \frac{1}{r^{d-1}}\frac{\partial}{\partial r} r^{d-1-\theta} \frac{\partial}{\partial r}$, and $d$, $\theta$ and $\nu$ are real parameters. This equation unifies the anomalous diffusion equation on fractals ($\nu =1$) and the spherical anomalous diffusion for porous media ($\theta=0$). Exact point-source solution is obtained, enabling us to describe a large class of subdiffusion ($\theta > (1-\nu)d$), normal diffusion ($\theta= (1-\nu)d$) and superdiffusion ($\theta < (1-\nu)d$). Furthermore, a thermostatistical basis for this solution is given from the maximum entropic principle applied to the Tsallis entropy., Comment: 3 pages, 2 eps figures

Feynman's path integral is herein generalized to the nonextensive canonical density matrix based on Tsallis entropy. This generalization is done in two ways by using unnormalized and normalized constraints. Firstly, we consider the path integral formulation with unnormalized constraints, and this generalization is worked out through two different ways, which are shown to be equivalent. These formulations with unnormalized constraints are solutions to two generalized Bloch equations proposed in this work. The first form of the generalized Bloch equation is linear, but with a temperature-dependent effective Hamiltonian; the second form is nonlinear and resembles the anomalous correlated diffusion equation (porous medium equation). Furthermore, we can extend these results to the prescription of field theory using integral representations. The second development is dedicated to analyzing the path integral formulation with normalized constraints. To illustrate the methods introduced here, we analyze the free particle case and a non-interacting scalar field. The results herein obtained are expected to be useful in the discussion of generic nonextensive contexts., Comment: (Univ. Est. de Maringa, PR- Brazil),17 pages, Latex

In this paper, the Green function theory of quantum many-particle systems recently presented is reworked within the framework of nonextensive statistical mechanics with a new normalized $q$-expectation values. This reformulation introduces a renormalization of temperature of the earlier theory and a self-consistency condition. The linear response theory is also presented, along with its two-particle Green function version. Finally, a Boltzmann transport-like equation is also developed here., Comment: RevTeX 16 pages

A unified presentation of the perturbation and variational methods for the generalized statistical mechanics based on Tsallis entropy is given here. In the case of the variational method, the Bogoliubov inequality is generalized in a very natural way following the Feynman proof for the usual statistical mechanics. The inequality turns out to be form-invariant with respect to the entropic index $q$. The method is illustrated with a simple example in classical mechanics. The formalisms developed here are expected to be useful in the discussion of nonextensive systems., Comment: revtex