Constructing a statistical mechanics for Beck-Cohen superstatistics

The basic aspects of both Boltzmann-Gibbs (BG) and nonextensive statistical mechanics can be seen through three different stages. First, the proposal of an entropic functional ($S_{BG} =-k\sum_i p_i \ln p_i$ for the BG formalism) with the appropriate constraints ($\sum_i p_i=1$ and $\sum_i p_i E_i = U$ for the BG canonical ensemble). Second, through optimization, the equilibrium or stationary-state distribution ($p_i = e^{-\beta E_i}/Z_{BG}$ with $Z_{BG}=\sum_j e^{-\beta E_j}$ for BG). Third, the connection to thermodynamics (e.g., $F_{BG}= -\frac{1}{\beta}\ln Z_{BG}$ and $U_{BG}=-\frac{\partial}{\partial \beta} \ln Z_{BG}$). Assuming temperature fluctuations, Beck and Cohen recently proposed a generalized Boltzmann factor $B(E) = \int_0^\infty d\beta f(\beta) e^{-\beta E}$. This corresponds to the second stage above described. In this letter we solve the corresponding first stage, i.e., we present an entropic functional and its associated constraints which lead precisely to $B(E)$. We illustrate with all six admissible examples given by Beck and Cohen.Comment: 3 PS figure

Generalizing the Planck distribution

Along the lines of nonextensive statistical mechanics, based on the entropy $S_q = k(1- \sum_i p_i^q)/(q-1) (S_1=-k \sum_i p_i \ln p_i)$, and Beck-Cohen superstatistics, we heuristically generalize Planck's statistical law for the black-body radiation. The procedure is based on the discussion of the differential equation $dy/dx=-a_{1}y-(a_{q}-a_{1}) y^{q}$ (with $y(0)=1$), whose $q=2$ particular case leads to the celebrated law, as originally shown by Planck himself in his October 1900 paper. Although the present generalization is mathematically simple and elegant, we have unfortunately no physical application of it at the present moment. It opens nevertheless the door to a type of approach that might be of some interest in more complex, possibly out-of-equilibrium, phenomena.Comment: 6 pages, including 2 figures. To appear in {\it Complexity, Metastability and Nonextensivity}, Proc. 31st Workshop of the International School of Solid State Physics (20-26 July 2004, Erice-Italy), eds. C. Beck, A. Rapisarda and C. Tsallis (World Scientific, Singapore, 2005

Two-parameter generalization of the logarithm and exponential functions and Boltzmann-Gibbs-Shannon entropy

The $q$-sum $x \oplus_q y \equiv x+y+(1-q) xy$ ($x \oplus_1 y=x+y$) and the $q$-product $x\otimes_q y \equiv [x^{1-q} +y^{1-q}-1]^{\frac{1}{1-q}}$ ($x\otimes_1 y=x y$) emerge naturally within nonextensive statistical mechanics. We show here how they lead to two-parameter (namely, $q$ and $q^\prime$) generalizations of the logarithmic and exponential functions (noted respectively $\ln_{q,q^\prime}x$ and $e_{q,q^\prime}^{x}$), as well as of the Boltzmann-Gibbs-Shannon entropy $S_{BGS}\equiv -k \sum_{i=1}^Wp_i \ln p_i$ (noted $S_{q,q^\prime}$). The remarkable properties of the $(q,q^\prime)$-generalized logarithmic function make the entropic form $S_{q,q^\prime} \equiv k \sum_{i=1}^W p_i \ln_{q,q^\prime}(1/p_i)$ to satisfy, for large regions of $(q,q^\prime)$, important properties such as {\it expansibility}, {\it concavity} and {\it Lesche-stability}, but not necessarily {\it composability}.Comment: 9 pages, 4 figure

Edge of chaos of the classical kicked top map: Sensitivity to initial conditions

We focus on the frontier between the chaotic and regular regions for the classical version of the quantum kicked top. We show that the sensitivity to the initial conditions is numerically well characterised by $\xi=e_q^{\lambda_q t}$, where $e_{q}^{x}\equiv [ 1+(1-q) x]^{\frac{1}{1-q}} (e_1^x=e^x)$, and $\lambda_q$ is the $q$-generalization of the Lyapunov coefficient, a result that is consistent with nonextensive statistical mechanics, based on the entropy $S_q=(1- \sum_ip_i^q)/(q-1) (S_1 =-\sum_i p_i \ln p_i$). Our analysis shows that $q$ monotonically increases from zero to unity when the kicked-top perturbation parameter $\alpha$ increases from zero (unperturbed top) to $\alpha_c$, where $\alpha_c \simeq 3.2$. The entropic index $q$ remains equal to unity for $\alpha \ge \alpha_c$, parameter values for which the phase space is fully chaotic.Comment: To appear in "Complexity, Metastability and Nonextensivity" (World Scientific, Singapore, 2005), Eds. C. Beck, A. Rapisarda and C. Tsalli

Nonextensive statistical mechanics and central limit theorems II - Convolution of q-independent random variables

In this article we review recent generalisations of the central limit theorem for the sum of specially correlated (or q-independent) variables, focusing on q greater or equal than 1. Specifically, this kind of correlation turns the probability density function known as q-Gaussian, which emerges upon maximisation of the entropy Sq, into an attractor in probability space. Moreover, we also discuss a q-generalisation of a-stable Levy distributions which can as well be stable for this special kind of correlation.Within this context, we verify the emergence of a triplet of entropic indices which relate the form of the attractor, the correlation, and the scaling rate, similar to the q-triplet that connects the entropic indices characterising the sensitivity to initial conditions, the stationary state, and relaxation to the stationary state in anomalous systems.Comment: 14 pages, 4 figures, and 1 table. To appear in the Proceedings of the conference CTNEXT07, Complexity, Metastability and Nonextensivity, Catania, Italy, 1-5 July 2007, Eds. S. Abe, H.J. Herrmann, P. Quarati, A. Rapisarda and C. Tsallis (American Institute of Physics, 2008) in pres

Nonextensive statistical mechanics and central limit theorems I - Convolution of independent random variables and q-product

In this article we review the standard versions of the Central and of the Levy-Gnedenko Limit Theorems, and illustrate their application to the convolution of independent random variables associated with the distribution known as q-Gaussian. This distribution emerges upon extremisation of the nonadditive entropy, basis of nonextensive statistical mechanics. It has a finite variance for q 5/3. We exhibit that, in the case of (standard) independence, the q-Gaussian has either the Gaussian (if q 5/3) as its attractor in probability space. Moreover, we review a generalisation of the product, the q-product, which plays a central role in the approach of the specially correlated variables emerging within the nonextensive theory.Comment: 13 pages, 4 figures. To appear in the Proceedings of the conference CTNEXT07, Complexity, Metastability and Nonextensivity, Catania, Italy, 1-5 July 2007, Eds. S. Abe, H.J. Herrmann, P. Quarati, A. Rapisarda and C. Tsallis (American Institute of Physics, 2008) in pres

An evaluation of biotic ligand models predicting acute copper toxicity to Daphnia magna in wastewater effluent

This is the author's accepted manuscript. The final published article is available from the link below. Copyright @ 2010 SETAC.The toxicity of Cu to Daphnia magna was investigated in a series of 48-h immobilization assays in effluents from four wastewater treatment works. The assay results were compared with median effective concentration (EC50) forecasts produced by the HydroQual biotic ligand model (BLM), the refined D. magna BLM, and a modified BLM that was constructed by integrating the refined D. magna biotic ligand characterization with the Windermere humic aqueous model (WHAM) VI geochemical speciation model, which also accommodated additional effluent characteristics as model inputs. The results demonstrated that all the BLMs were capable of predicting toxicity by within a factor of two, and that the modified BLM produced the most accurate toxicity forecasts. The refined D. magna BLM offered the most robust assessment of toxicity in that it was not reliant on the inclusion of effluent characteristics or optimization of the dissolved organic carbon active fraction to produce forecasts that were accurate by within a factor of two. The results also suggested that the biotic ligand stability constant for Na may be a poor approximation of the mechanisms governing the influence of Na where concentrations exceed the range within which the biotic ligand stability constant value had been determined. These findings support the use of BLMs for the establishment of site-specific water quality standards in waters that contain a substantial amount of wastewater effluent, but reinforces the need for regulators to scrutinize the composition of models, their thermodynamic and biotic ligand parameters, and the limitations of those parameters.EPSRC and Severn Trent Water