Superdiffusion and non-Gaussian statistics in a driven-dissipative 2D dusty plasma

Superstatistics from a dynamical perspective: Entropy and relaxation

Distributions that deviate from equilibrium predictions are commonly observed across a broad spectrum of systems, ranging from laboratory experiments to astronomical phenomena. These distributions are generally regarded as a manifestation of a quasiequilibrium state and can very often be represented as a superposition of statistics, i.e., superstatistics. The underlying idea in this methodology is that the nonequilibrium system consists of a collection of smaller subsystems that remain infinitely close to equilibrium. This procedure has been effectively implemented in a kinetic setting, but thus far, only in the collisionless regime, limiting its scope of application. In this paper, we address the effect of collisions on the relaxation process and time evolution of superstatistical systems. After confronting the superstatistical distributions with experimental and simulation data, relevant to our analysis, we first study the effect of superstatistics on entropy. We explicitly show that, in the absence of long-range interactions, the extensivity of entropy is preserved, albeit influenced by the specific class of temperature fluctuations. Then, we examine how collisions drive the system towards a global equilibrium state, characterized by a maximum entropy, by employing the relaxation time approximation. This allows us to define a dynamical version of superstatistics, characterized by a singular time-varying parameter q(t), which undergoes a continuous evolution towards equilibrium. We show how this approach enables the determination of the evolution of the underlying temperature distribution under the influence of collisions, which act as stochastic forces, gradually narrowing the temperature distribution over time.

Entropy defect in thermodynamics

This paper describes the physical foundations of the newly discovered "entropy defect" as a basic concept of thermodynamics. The entropy defect quantifies the change in entropy caused by the order induced in a system through the additional correlations among its constituents when two or more subsystems are assembled. This defect is closely analogous to the mass defect that arises when nuclear particle systems are assembled. The entropy defect determines how the entropy of the system compares to its constituent's entropies and stands on three fundamental properties: each constituent's entropy must be (i) separable, (ii) symmetric, and (iii) bounded. We show that these properties provide a solid foundation for the entropy defect and for generalizing thermodynamics to describe systems residing out of the classical thermal equilibrium, both in stationary and nonstationary states. In stationary states, the consequent thermodynamics generalizes the classical framework, which was based on the Boltzmann-Gibbs entropy and Maxwell-Boltzmann canonical distribution of particle velocities, into the respective entropy and canonical distribution associated with kappa distributions. In nonstationary states, the entropy defect similarly acts as a negative feedback, or reduction of the increase of entropy, preventing its unbounded growth toward infinity.

Comparison of the static structure factor at long wavelengths for a dusty plasma liquid and other liquids

Especially small values of the static structure factor S(k) at long wavelengths, i.e., small k, were obtained in an analysis of experimental data, for a two-dimensional dusty plasma in its liquid state. For comparison, an analysis of S(k) data was carried out for many previously published experiments with other liquids. The latter analysis indicates that the magnitude of S(k) at small k is typically in a range 0.02-0.13. In contrast, the corresponding value for a dusty plasma liquid was found to be as small as 0.0139. Another basic finding for the dusty plasma liquid is that S(k) at small k generally increases with temperature, with its lowest value, noted above, occurring near the melting point. Simulations were carried out for the dusty plasma liquid, and their results are generally consistent with the experiment. Since a dusty plasma has a soft interparticle interaction, our findings support earlier theoretical suggestions that a useful design strategy for creating materials having exceptionally low values of S(0), so-called hyperuniform materials, is the use of a condensed material composed of particles that interact softly at their periphery.

The Longitudinal Plasma Modes of κ-Deformed Kaniadakis Distributed Plasmas Carrying Orbital Angular Momentum

Based on plasma kinetic theory, the dispersion and Landau damping of Langmuir and ion-acoustic waves carrying finite orbital angular momentum (OAM) were investigated in the κ-deformed Kaniadakis distributed plasma system. The results showed that the peculiarities of the investigated subjects relied on the deformation parameter κ and OAM parameter η. For both Langmuir and ion-acoustic waves, dispersion was enhanced with increased κ, while the Landau damping was suppressed. Conversely, both the dispersion and Landau damping were depressed by OAM. Moreover, the results coincided with the straight propagating plane waves in a Maxwellian plasma system when κ=0 and η→∞. It was expected that the present results would give more insight into the trapping and transportation of plasma particles and energy.

Directional locking in a two-dimensional Yukawa solid modulated by a two-dimensional periodic substrate

Directional depinning dynamics of a two-dimensional (2D) dusty plasma solid modulated by a 2D square periodic substrate are investigated using Langevin dynamical simulations. We observe prominent directional locking effects when the direction of the external driving force is varied relative to the underlying square substrate. These locking steps appear when the direction of the driving force is close to the symmetry direction of the substrate, corresponding to the different dynamical flow patterns and the structures. In the conditions between the adjacent locking steps, moving ordered states are observed. Although the discontinuous transitions often occur between the locking steps and the nonlocking portion, the continuous transitions are also found around the locking step associated with the disordered plastic flow close to its termini. Our results show that directional locking also occurs for underdamped systems, which could be tested experimentally in dusty plasmas modulated by 2D substrates.

Initial measurement of ion nonextensive parameter with geodesic acoustic mode theory

The consideration of nonextensivity effects is crucial to the accurate diagnosis of plasma parameters; common plasma nonextensive parameters include electron nonextensive parameter and ion nonextensive parameter, and the former can be measured, while the ion nonextensive parameter cannot be measured yet. Here we show the measurement of ion nonextensive parameter of plasma based on the theory of nonextensive geodesic acoustic modes. We assume that the plasma to be measured can be described by nonextensive statistical mechanics, and on this basis, the nonextensive geodesic acoustic mode theory is established. Utilizing this theory, we have measured the ion nonextensive parameter [Formula: see text] which cannot be diagnosed even by a nonextensive single electric probe. Our research points out that the proposed measurement method of ion nonextensive parameter may play a role in plasma diagnosis and will help us to grasp the nonextensivity of plasma more precisely. We hope the proposed method of ion nonextensive parameter diagnosis based on the nonextensive geodesic acoustic mode theory can be the starting point of more complex ion nonextensive parameter diagnosis methods. In addition, the measurement of ion nonextensive parameter is closely related to the study of various plasma waves, instabilities, turbulence and abnormal transport, and a defined and quantitative test of nonextensive geodesic acoustic mode theory will bound up deeply with such developments.

Shear softening and hardening of a two-dimensional Yukawa solid

Langevin dynamical simulations are performed to study the elastic behaviors of two-dimensional (2D) solid dusty plasmas under the periodic shear deformation. The frequency- and strain-dependent shear moduli G(ω,γ) of our simulated 2D Yukawa solid are calculated from the ratio of the shear stress to strain in different orientations. The shear-softening and -hardening properties in different lattice orientations are discovered from the obtained G(ω,γ). The component of the elastic constant tensor corresponding to the shear deformation is also calculated, whose variation trend exactly agrees with the discovered shear-softening and -hardening features in different shear directions. It is also found that the shear modulus of the 2D Yukawa solid always increases monotonically with the frequency.

Determination of best particle tracking velocimetry method for two-dimensional dusty plasmas

A practical procedure of determining the best choice of particle tracking velocimetry (PTV) methods for the analysis of two-dimensional (2D) dusty plasmas is provided. Using the measured particle positions with errors, the algorithm of PTV introduces further systematic errors while calculating particle velocities. To better quantify the errors, computer simulations of 2D Yukawa liquids are performed so that the true particle positions and velocities are precisely known. Synthetic data are achieved by varying the time interval between frames and adding the uncertainty with the different magnitude to the particle position data. Various PTV methods are used to calculate the velocity data from the achieved synthetic data, and the obtained velocity data from these PTV methods are compared with the true velocity data using two diagnostics to determine the best PTV method for various conditions. The results of various PTV methods with the simulation data are further confirmed by the demonstration using the 2D dusty plasma experimental data. This practical procedure is also applicable to determine the best PTV method for dusty plasmas in different conditions, such as those in Tokamaks, and for other physical systems, including colloids, granular materials, and fluid mechanics.

Phonon spectra of a two-dimensional solid dusty plasma modified by two-dimensional periodic substrates

Phonon spectra of a two-dimensional (2D) solid dusty plasma modulated by 2D square and triangular periodic substrates are investigated using Langevin dynamical simulations. The commensurability ratio, i.e., the ratio of the number of particles to the number of potential well minima, is set to 1 or 2. The resulting phonon spectra show that propagation of waves is always suppressed due to the confinement of particles by the applied 2D periodic substrates. For a commensurability ratio of 1, the spectra indicate that all particles mainly oscillate at one specific frequency, corresponding to the harmonic oscillation frequency of one single particle inside one potential well. At a commensurability ratio of 2, the substrate allows two particles to sit inside the bottom of each potential well, and the resulting longitudinal and transverse spectra exhibit four branches in total. We find that the two moderate branches come from the harmonic oscillations of one single particle and two combined particles in the potential well. The other two branches correspond to the relative motion of the two-body structure in each potential well in the radial and azimuthal directions. The difference in the spectra between the square and triangular substrates is attributed to the anisotropy of the substrates and the resulting alignment directions of the two-body structure in each potential well.

Polymer pattern-induced self-assembly of inorganic nanoparticles

Functional assemblies of inorganic nanoparticles (NPs) are widely studied owing to their collective electromagnetic properties and various application from nanodrugs and bioimaging. In most cases, the superstructures of NPs are prepared with the assistance of templates or external fields. Therefore, how to prepare the functional assemblies of NPs more simply remains a challenge. Here, a free-template assembly strategy for preparing the superstructures of NPs is proposed in our work. In our strategy, we design poly(glycerol monomethacrylate)-b-poly(2-hydroxypropyl methacrylate) (PGMA-b-PHPMA) coated NPs. Then, using the polymerization-induced self-assembly (PISA), hydrophobic PHPMA blocks resulted in the phase separation to form the orderly patterns, which is expected to induced NPs to self-assemble into the orderly superstructures. By DPD simulations, we find that the disk, ring, composite superstructures can be obtained by regulating the graft density, verifying that our assembly strategy of NPs is feasible. Even more interesting is that NPs are also distributed in an orderly way on the surface of aggregations to form the orderly NP patterns. Besides that, the thermodynamics, dynamics, and structure details in the self-assembly process of HINPs are shown in our work, providing a new idea and elaborate physical picture for the following preparation of the superstructure of NPs.

The initiation of shear band formation in deformed metallic glasses from soft localized domains

It has long been thought that shear band (SB) formation in amorphous solids initiates from relatively "soft" regions in the material in which large-scale non-affine deformations become localized. The test of this hypothesis requires an effective means of identifying "soft" regions and their evolution as the material is deformed to varying degrees, where the metric of "softness" must also account for the effect of temperature on local material stiffness. We show that the mean square atomic displacement on a caging timescale ⟨u²⟩, the "Debye-Waller factor," provides a useful method for estimating the shear modulus of the entire material and, by extension, the material stiffness at an atomic scale. Based on this "softness" metrology, we observe that SB formation indeed occurs through the strain-induced formation of localized soft regions in our deformed metallic glass free-standing films. Unexpectedly, the critical strain condition for SB formation occurs when the softness (⟨u²⟩) distribution within the emerging soft regions approaches that of the interfacial region in its undeformed state, initiating an instability with similarities to the transition to turbulence. Correspondingly, no SBs arise when the material is so thin that the entire material can be approximately described as being "interfacial" in nature. We also quantify relaxation in the glass and the nature and origin of highly non-Gaussian particle displacements in the dynamically heterogeneous SB regions at times longer than the caging time.

Plastic strain rate quantified from dislocation dynamics in dusty plasma shear flows

Dynamics of dislocations and defects are investigated in 2D dusty plasma experiments with two counterpropagating flows. It is experimentally demonstrated that the Orowan equation is able to accurately determine the plastic strain rate from the motion of dislocations, well agreeing with the shear rate defined from the drift velocity gradient. For a higher shear rate, the studied system is in the liquidlike flow state, as a result, the determined shear rate from the Orowan equation deviates from its definition. The obtained probability distribution function of dislocations from the experiments clearly shows that the dislocation motion can be divided into the local and gliding ones. All findings above are further verified by the corresponding Langevin dynamical simulations with various levels of shear rates. The dislocation and defect analysis results from these simulations clearly indicate that the defect and dislocation dynamics in the sheared dusty plasmas clearly exhibit two stages as the shear rate increases.

Fingerprints of nonequilibrium stationary distributions in dispersion relations

Distributions different from those predicted by equilibrium statistical mechanics are commonplace in a number of physical situations, such as plasmas and self-gravitating systems. The best strategy for probing these distributions and unavailing their origins consists in combining theoretical knowledge with experiments, involving both direct and indirect measurements, as those associated with dispersion relations. This paper addresses, in a quite general context, the signature of nonequilibrium distributions in dispersion relations. We consider the very general scenario of distributions corresponding to a superposition of equilibrium distributions, that are well-suited for systems exhibiting only local equilibrium, and discuss the general context of systems obeying the combination of the Schrödinger and Poisson equations, while allowing the Planck's constant to smoothly go to zero, yielding the classical kinetic regime. Examples of media where this approach is applicable are plasmas, gravitational systems, and optical molasses. We analyse in more depth the case of classical dispersion relations for a pair plasma. We also discuss a possible experimental setup, based on spectroscopic methods, to directly observe these classes of distributions.

Heterogeneous diffusion processes and nonergodicity with Gaussian colored noise in layered diffusivity landscapes

Heterogeneous diffusion processes (HDPs) with space-dependent diffusion coefficients D(x) are found in a number of real-world systems, such as for diffusion of macromolecules or submicron tracers in biological cells. Here, we examine HDPs in quenched-disorder systems with Gaussian colored noise (GCN) characterized by a diffusion coefficient with a power-law dependence on the particle position and with a spatially random scaling exponent. Typically, D(x) is considered to be centerd at the origin and the entire x axis is characterized by a single scaling exponent α. In this work we consider a spatially random scenario: in periodic intervals ("layers") in space D(x) is centerd to the midpoint of each interval. In each interval the scaling exponent α is randomly chosen from a Gaussian distribution. The effects of the variation of the scaling exponents, the periodicity of the domains ("layer thickness") of the diffusion coefficient in this stratified system, and the correlation time of the GCN are analyzed numerically in detail. We discuss the regimes of superdiffusion, subdiffusion, and normal diffusion realisable in this system. We observe and quantify the domains where nonergodic and non-Gaussian behaviors emerge in this system. Our results provide new insights into the understanding of weak ergodicity breaking for HDPs driven by colored noise, with potential applications in quenched layered systems, typical model systems for diffusion in biological cells and tissues, as well as for diffusion in geophysical systems.

Determination of viscosity in shear-induced melting two-dimensional dusty plasmas using Green-Kubo relation

Langevin dynamical simulations of shear-induced melting two-dimensional (2D) dusty plasmas are performed to study the determination of the shear viscosity of this system. It is found that the viscosity calculated from the Green-Kubo relation, after removing the drift motion, well agrees with the viscosity definition, i.e., the ratio of the shear stress to the shear rate in the sheared region, even the shear rate is magnified ten times higher than that in experiments. The behaviors of shear stress and its autocorrelation function of shear-induced melting 2D dusty plasmas are compared with those of uniform liquids at the same temperatures, leading to the conclusion that the Green-Kubo relation is still applicable to determine the viscosity for shear-induced melting dusty plasmas.

Entropy formula of N-body system

We prove a proposition that the entropy of the system composed of finite N molecules of ideal gas is the q-entropy or Havrda–Charvát–Tsallis entropy, which is also known as Tsallis entropy, with the entropic index \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q=\frac{D(N-1)-4}{D(N-1)-2}$$\end{document}q=D(N-1)-4D(N-1)-2 in D-dimensional space. The indispensable infinity assumption used by Boltzmann and others in their derivation of entropy formulae is not involved in our derivation, therefore our derived formula is exact. The analogy of the N-body system brings us to obtain the entropic index of a combined system \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q_C$$\end{document}qC formed from subsystems having different entropic indexes \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q_A$$\end{document}qA and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q_B$$\end{document}qB as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{1}{1-q_C}=\frac{1}{1-q_A}+\frac{1}{1-q_B}+\frac{D+2}{2}$$\end{document}11-qC=11-qA+11-qB+D+22. It is possible to use the number N for the physical measure of deviation from Boltzmann entropy.

Dynamical states in two-dimensional charged dust particle clusters in plasma medium

The formation of dynamical states for a collection of dust particles in two dimensions is shown using molecular dynamics simulations. The charged dust particles interact with each other with a Yukawa pair potential mimicking the screening due to plasma. An external radial confining force has also been applied to the dust particles to keep them radially confined. When the particle number is low (say, a few), they get arranged on the radial locations corresponding to multiple rings or shells. For specific numbers, such an arrangement of particles is stationary. However, for several cases, the cluster of dust particles relaxes to a state for which the dust particles on rings display intershell rotation. For a larger number of dust particles (a few hundred, for instance), an equilibrium state with a coherent rigid body displaying angular oscillation of the entire cluster is observed. A detailed characterization of the formation of these states in terms of particle number, coupling parameter, etc., is provided.

Initial measurement of electron nonextensive parameter with electric probe

Theoretical analysis and a large number of experiments have proved that plasma components do not satisfy Boltzmann-Gibbs statistics and can be well described by nonextensive statistical mechanics, while new plasma parameters, electron nonextensive parameters, which are introduced to describe the nonextensive properties of plasma, cannot be diagnosed yet. Here we show measurement of electron nonextensive parameters of plasma with a nonextensive single electric probe. Our results show that nonextensive electric probe may play a role in plasma diagnosis, measuring nonextensivity of plasma and improving diagnostic accuracy of other plasma parameters. We expect the proposed nonextensive single electric probe can be starting point of more complex nonextensive electric probe. In addition, nonextensive electric probe is an important means to study various plasma waves and instability, turbulence, and anomalous transport, and a definite and quantitative test of the theory of nonextensive geodetic acoustic models will be relevant to such development.

Heterogeneity of dynamics in a modulated colloidal liquid

We study the dynamics of a system of two dimensional colloidal particles subjected to a spatially periodic external potential using Brownian dynamics simulations. We characterize the dynamics in the system by the mean square displacements and the self-van Hove function. The static density plots suggest that system gets into modulated liquid phase in presence of the external potential. We find that diffusion coefficients, obtained from long time mean sqaure displacements, decay exponentially with increasing potential strength. The self-van Hove functions computed from the distribution of particle displacemets in a given time interval show non-gaussian behaviour in directions both parallel and transverse to the external modulation. This suggests heterogeneous dynamics and is supported by particle mobilities and residence times.

Entropy Measurement in Strongly Coupled Complex Plasmas

The change in entropy of a system that is transferred between two states at different temperatures is measured in a two-dimensional plasma crystal experiment. One- and especially two-component dust clusters are confined in the plasma sheath and heated to different temperatures using laser manipulation. We find that entropies obtained from the phase space yield consistent results for, i.e., the heat capacity which shows excellent agreement with the Dulong-Petit law. The implications for the validity of basic thermodynamical principles in finite size complex (dusty) plasmas are discussed.

Depinning dynamics of two-dimensional dusty plasmas on a one-dimensional periodic substrate

We investigate the depinning dynamics of two-dimensional dusty plasmas driven over one-dimensional periodic substrates using Langevin dynamical simulations. We find that, for a specific range of substrate strengths, as the external driving force increases from zero, there are three different states, which are the pinned, the disordered plastic flow, and the moving ordered states, respectively. These three states are clearly observed using different diagnostics, including the collective drift velocity, static structural measures, the particle trajectories, the mean-squared displacements, and the kinetic temperature. We compare the observed depinning dynamics here with the depinning dynamics in other systems.

From the nonlinear Fokker-Planck equation to the Vlasov description and back: Confined interacting particles with drag

Nonlinear Fokker-Planck equations endowed with power-law diffusion terms have proven to be valuable tools for the study of diverse complex systems in physics, biology, and other fields. The nonlinearity appearing in these evolution equations can be interpreted as providing an effective description of a system of particles interacting via short-range forces while performing overdamped motion under the effect of an external confining potential. This point of view has been recently applied to the study of thermodynamical features of interacting vortices in type II superconductors. In the present work we explore an embedding of the nonlinear Fokker-Planck equation within a Vlasov equation, thus incorporating inertial effects to the concomitant particle dynamics. Exact time-dependent solutions of the q-Gaussian form (with compact support) are obtained for the Vlasov equation in the case of quadratic confining potentials.

A Generalized Relative (α, β)-Entropy: Geometric Properties and Applications to Robust Statistical Inference

Entropy and relative entropy measures play a crucial role in mathematical information theory. The relative entropies are also widely used in statistics under the name of divergence measures which link these two fields of science through the minimum divergence principle. Divergence measures are popular among statisticians as many of the corresponding minimum divergence methods lead to robust inference in the presence of outliers in the observed data; examples include the ϕ -divergence, the density power divergence, the logarithmic density power divergence and the recently developed family of logarithmic super divergence (LSD). In this paper, we will present an alternative information theoretic formulation of the LSD measures as a two-parameter generalization of the relative α -entropy, which we refer to as the general ( α , β ) -entropy. We explore its relation with various other entropies and divergences, which also generates a two-parameter extension of Renyi entropy measure as a by-product. This paper is primarily focused on the geometric properties of the relative ( α , β ) -entropy or the LSD measures; we prove their continuity and convexity in both the arguments along with an extended Pythagorean relation under a power-transformation of the domain space. We also derive a set of sufficient conditions under which the forward and the reverse projections of the relative ( α , β ) -entropy exist and are unique. Finally, we briefly discuss the potential applications of the relative ( α , β ) -entropy or the LSD measures in statistical inference, in particular, for robust parameter estimation and hypothesis testing. Our results on the reverse projection of the relative ( α , β ) -entropy establish, for the first time, the existence and uniqueness of the minimum LSD estimators. Numerical illustrations are also provided for the problem of estimating the binomial parameter.

Adiabatic Expansion of Electron Gas in a Magnetic Nozzle

A specially constructed experiment shows the near perfect adiabatic expansion of an ideal electron gas resulting in a polytropic index greater than 1.4, approaching the adiabatic value of 5/3, when removing electric fields from the system, while the polytropic index close to unity is observed when the electrons are trapped by the electric fields. The measurements were made on collisionless electrons in an argon plasma expanding in a magnetic nozzle. The collision lengths of all electron collision processes are greater than the scale length of the expansion, meaning the system cannot be in thermodynamic equilibrium, yet thermodynamic concepts can be used, with caution, in explaining the results. In particular, a Lorentz force, created by inhomogeneities in the radial plasma density, does work on the expanding magnetic field, reducing the internal energy of the electron gas that behaves as an adiabatically expanding ideal gas.

Confined sandpile in two dimensions: Percolation and singular diffusion

We investigate the properties of a two-state sandpile model subjected to a confining potential in two dimensions. From the microdynamical description, we derive a diffusion equation, and find a stationary solution for the case of a parabolic confining potential. By studying the systems at different confining conditions, we observe two scale-invariant regimes. At a given confining potential strength, the cluster size distribution takes the form of a power law. This regime corresponds to the situation in which the density at the center of the system approaches the critical percolation threshold. The analysis of the fractal dimension of the largest cluster frontier provides evidence that this regime is reminiscent of gradient percolation. By increasing further the confining potential, most of the particles coalesce in a giant cluster, and we observe a regime where the jump size distribution takes the form of a power law. The onset of this second regime is signaled by a maximum in the fluctuation of energy.

Perturbative thermodynamic geometry of nonextensive ideal classical, Bose, and Fermi gases

We investigate perturbative thermodynamic geometry of nonextensive ideal classical, Bose, and Fermi gases. We show that the intrinsic statistical interaction of nonextensive Bose (Fermi) gas is attractive (repulsive) similar to the extensive case but the value of thermodynamic curvature is changed by a nonextensive parameter. In contrary to the extensive ideal classical gas, the nonextensive one may be divided to two different regimes. According to the deviation parameter of the system to the nonextensive case, one can find a special value of fugacity, z^{*}, where the sign of thermodynamic curvature is changed. Therefore, we argue that the nonextensive parameter induces an attractive (repulsive) statistical interaction for z<z^{*} (z>z^{*}) for an ideal classical gas. Also, according to the singular point of thermodynamic curvature, we consider the condensation of nonextensive Boson gas.

Diffusion and Directionality of Charged Nanoparticles on Lipid Bilayer Membrane

Diffusion dynamics of charged nanoparticles on the lipid membrane is of essential importance to cellular functioning. Yet a fundamental insight into electrostatics-mediated diffusion dynamics of charged nanoparticles on the membrane is lacking and remains to be an urgent issue. Here we present the computational investigation to uncover the pivotal role of electrostatics in the diffusion dynamics of charged nanoparticles on the lipid membrane. Our results demonstrate diffusive behaviors and directional transport of a charged nanoparticle, significantly depending on the sign and spatial distribution of charges on its surface. In contrast to the Fickian diffusion of neutral nanoparticles, randomly charged nanoparticles undergo superdiffusive transport with directionality. However, the dynamics of uniformly charged nanoparticles favors Fickian diffusion that is significantly enhanced. Such observations can be explained in term of electrostatics-induced surface reconstruction and fluctuation of lipid membrane. We finally present an analytical model connecting surface reconstruction and local deformation of the membrane. Our findings bear wide implications for the understanding and control of the transport of charged nanoparticles on the cell membrane.

Effective-medium approximation for lattice random walks with long-range jumps

We consider the random walk on a lattice with random transition rates and arbitrarily long-range jumps. We employ Bruggeman's effective-medium approximation (EMA) to find the disorder-averaged (coarse-grained) dynamics. The EMA procedure replaces the disordered system with a cleverly guessed reference system in a self-consistent manner. We give necessary conditions on the reference system and discuss possible physical mechanisms of anomalous diffusion. In the case of a power-law scaling between transition rates and distance, lattice variants of Lévy-flights emerge as the effective medium, and the problem is solved analytically, bearing the effective anomalous diffusivity. Finally, we discuss several example distributions and demonstrate very good agreement with numerical simulations.

Repulsive particles under a general external potential: Thermodynamics by neglecting thermal noise

A recent proposal of an effective temperature θ, conjugated to a generalized entropy s_{q}, typical of nonextensive statistical mechanics, has led to a consistent thermodynamic framework in the case q=2. The proposal was explored for repulsively interacting vortices, currently used for modeling type-II superconductors. In these systems, the variable θ presents values much higher than those of typical room temperatures T, so that the thermal noise can be neglected (T/θ≃0). The whole procedure was developed for an equilibrium state obtained after a sufficiently long-time evolution, associated with a nonlinear Fokker-Planck equation and approached due to a confining external harmonic potential, ϕ(x)=αx^{2}/2 (α>0). Herein, the thermodynamic framework is extended to a quite general confining potential, namely ϕ(x)=α|x|^{z}/z (z>1). It is shown that the main results of the previous analyses hold for any z>1: (i) The definition of the effective temperature θ conjugated to the entropy s_{2}. (ii) The construction of a Carnot cycle, whose efficiency is shown to be η=1-(θ_{2}/θ_{1}), where θ_{1} and θ_{2} are the effective temperatures associated with two isothermal transformations, with θ_{1}>θ_{2}. The special character of the Carnot cycle is indicated by analyzing another cycle that presents an efficiency depending on z. (iii) Applying Legendre transformations for a distinct pair of variables, different thermodynamic potentials are obtained, and furthermore, Maxwell relations and response functions are derived. The present approach shows a consistent thermodynamic framework, suggesting that these results should hold for a general confining potential ϕ(x), increasing the possibility of experimental verifications.

Role of dimensionality in complex networks

Deep connections are known to exist between scale-free networks and non-Gibbsian statistics. For example, typical degree distributions at the thermodynamical limit are of the form , where the q-exponential form optimizes the nonadditive entropy S_q (which, for q → 1, recovers the Boltzmann-Gibbs entropy). We introduce and study here d-dimensional geographically-located networks which grow with preferential attachment involving Euclidean distances through . Revealing the connection with q-statistics, we numerically verify (for d = 1, 2, 3 and 4) that the q-exponential degree distributions exhibit, for both q and k, universal dependences on the ratio α_A/d. Moreover, the q = 1 limit is rapidly achieved by increasing α_A/d to infinity.

Quasi-two-dimensional complex plasma containing spherical particles and their binary agglomerates

A unique type of quasi-two-dimensional complex plasma system was observed which consisted of monodisperse microspheres and their binary agglomerations (dimers). The particles and their dimers levitated in a plasma sheath at slightly different heights and formed two distinct sublayers. The system did not crystallize and may be characterized as a disordered solid. The dimers were identified based on their characteristic appearance in defocused images, i.e., rotating interference fringe patterns. The in-plane and interplane particle separations exhibit nonmonotonic dependence on the discharge pressure.

Sensitivity to initial conditions of a d-dimensional long-range-interacting quartic Fermi-Pasta-Ulam model: Universal scaling

We introduce a generalized d-dimensional Fermi-Pasta-Ulam model in the presence of long-range interactions, and perform a first-principle study of its chaos for d=1,2,3 through large-scale numerical simulations. The nonlinear interaction is assumed to decay algebraically as d_{ij}^{-α} (α≥0), {d_{ij}} being the distances between N oscillator sites. Starting from random initial conditions we compute the maximal Lyapunov exponent λ_{max} as a function of N. Our N≫1 results strongly indicate that λ_{max} remains constant and positive for α/d>1 (implying strong chaos, mixing, and ergodicity), and that it vanishes like N^{-κ} for 0≤α/d<1 (thus approaching weak chaos and opening the possibility of breakdown of ergodicity). The suitably rescaled exponent κ exhibits universal scaling, namely that (d+2)κ depends only on α/d and, when α/d increases from zero to unity, it monotonically decreases from unity to zero, remaining so for all α/d>1. The value α/d=1 can therefore be seen as a critical point separating the ergodic regime from the anomalous one, κ playing a role analogous to that of an order parameter. This scaling law is consistent with Boltzmann-Gibbs statistics for α/d>1, and possibly with q statistics for 0≤α/d<1.

Dust interferometers in plasmas

An interferometric imaging technique has been proposed to instantly measure the diameter of individual spherical dust particles suspended in a gas discharge plasma. The technique is based on the defocused image analysis of both spherical particles and their binary agglomerates. Above a critical diameter, the defocused images of spherical particles contain stationary interference fringe patterns and the fringe number increases with particle diameters. Below this critical diameter, the particle size has been measured using the rotational interference fringe patterns which appear only on the defocused images of binary agglomerates. In this case, a lower cutoff limit of particle diameter has been predicted, below which no such rotational fringe patterns are observed for the binary agglomerates. The method can be useful as a diagnostics for complex plasma experiments on earth as well as under microgravity conditions.

Phase Separation of Binary Charged Particle Systems with Small Size Disparities using a Dusty Plasma

The phase separation in binary mixtures of charged particles has been investigated in a dusty plasma under microgravity on parabolic flights. A method based on the use of fluorescent dust particles was developed that allows us to distinguish between particles of slightly different size. A clear trend towards phase separation even for smallest size (charge) disparities is observed. The diffusion flux is directly measured from the experiment and uphill diffusion coefficients have been determined.

Validity of the third law of thermodynamics for the Tsallis entropy

Bento et al. [Phys. Rev. E 91, 022105 (2015)] recently stated that the Tsallis entropy violates the third law of thermodynamics for 0<q<1 in the subadditive regime. We first show that the division between the regimes q<1 and q>1 is already inherent in the fundamental incomplete structure of the deformed logarithms and exponentials underlying the Tsallis entropy. Then, we provide the complete deformed functions and show that the Tsallis entropy conforms to the third law of thermodynamics for both superadditive q<1 and subadditive q>1 regimes. Finally, we remark that the Tsallis entropy does not require the use of an escort-averaging scheme once it is expressed in terms of the complete deformed functions.

Replica Symmetry Breaking in Trajectories of a Driven Brownian Particle

We study a Brownian particle passively driven by a field obeying the noisy Burgers' equation. We demonstrate that the system exhibits replica symmetry breaking in the path ensemble with the initial position of the particle being fixed. The key step of the proof is that the path ensemble with a modified boundary condition can be exactly mapped onto the canonical ensemble of directed polymers.

Extended q-Gaussian and q-exponential distributions from gamma random variables

The family of q-Gaussian and q-exponential probability densities fit the statistical behavior of diverse complex self-similar nonequilibrium systems. These distributions, independently of the underlying dynamics, can rigorously be obtained by maximizing Tsallis "nonextensive" entropy under appropriate constraints, as well as from superstatistical models. In this paper we provide an alternative and complementary scheme for deriving these objects. We show that q-Gaussian and q-exponential random variables can always be expressed as a function of two statistically independent gamma random variables with the same scale parameter. Their shape index determines the complexity q parameter. This result also allows us to define an extended family of asymmetric q-Gaussian and modified q-exponential densities, which reduce to the standard ones when the shape parameters are the same. Furthermore, we demonstrate that a simple change of variables always allows relating any of these distributions with a beta stochastic variable. The extended distributions are applied in the statistical description of different complex dynamics such as log-return signals in financial markets and motion of point defects in a fluid flow.

Consistent thermodynamic framework for interacting particles by neglecting thermal noise

An effective temperature θ, conjugated to a generalized entropy s(q), was introduced recently for a system of interacting particles. Since θ presents values much higher than those of typical room temperatures T≪θ, the thermal noise can be neglected (T/θ≃0) in these systems. Moreover, the consistency of this definition, as well as of a form analogous to the first law of thermodynamics, du=θds(q)+δW, were verified lately by means of a Carnot cycle, whose efficiency was shown to present the usual form, η=1-(θ(2)/θ(1)). Herein we explore further the heat contribution δQ=θds(q) by proposing a way for a heat exchange between two such systems, as well as its associated thermal equilibrium. As a consequence, the zeroth principle is also established. Moreover, we consolidate the first-law proposal by following the usual procedure for obtaining different potentials, i.e., applying Legendre transformations for distinct pairs of independent variables. From these potentials we derive the equation of state, Maxwell relations, and define response functions. All results presented are shown to be consistent with those of standard thermodynamics for T>0.

Third law of thermodynamics as a key test of generalized entropies

The laws of thermodynamics constrain the formulation of statistical mechanics at the microscopic level. The third law of thermodynamics states that the entropy must vanish at absolute zero temperature for systems with nondegenerate ground states in equilibrium. Conversely, the entropy can vanish only at absolute zero temperature. Here we ask whether or not generalized entropies satisfy this fundamental property. We propose a direct analytical procedure to test if a generalized entropy satisfies the third law, assuming only very general assumptions for the entropy S and energy U of an arbitrary N-level classical system. Mathematically, the method relies on exact calculation of β=dS/dU in terms of the microstate probabilities p(i). To illustrate this approach, we present exact results for the two best known generalizations of statistical mechanics. Specifically, we study the Kaniadakis entropy S(κ), which is additive, and the Tsallis entropy S(q), which is nonadditive. We show that the Kaniadakis entropy correctly satisfies the third law only for -1<κ<+1, thereby shedding light on why κ is conventionally restricted to this interval. Surprisingly, however, the Tsallis entropy violates the third law for q<1. Finally, we give a concrete example of the power of our proposed method by applying it to a paradigmatic system: the one-dimensional ferromagnetic Ising model with nearest-neighbor interactions.

Superdiffusion of two-dimensional Yukawa liquids due to a perpendicular magnetic field

Stochastic transport of a two-dimensional (2D) dusty plasma liquid with a perpendicular magnetic field is studied. Superdiffusion is found to occur especially at higher magnetic fields with β of order unity. Here, β = ω(c)/ω(pd) is the ratio of the cyclotron and plasma frequencies for dust particles. The mean-square displacement MSD = 4D(α)t(α) is found to have an exponent α > 1, indicating superdiffusion, with α increasing monotonically to 1.1 as β increases to unity. The 2D Langevin molecular dynamics simulation used here also reveals that another indicator of random particle motion, the velocity autocorrelation function, has a dominant peak frequency ω(peak) that empirically obeys ω(peak)(2) = ω(c)(2) + ω(pd)(2)/4.

Perpendicular diffusion of a dilute beam of charged dust particles in a strongly coupled dusty plasma

The diffusion of projectiles drifting through a target of strongly coupled dusty plasma is investigated in a simulation. A projectile's drift is driven by a constant force F. We characterize the random walk of the projectiles in the direction perpendicular to their drift. The perpendicular diffusion coefficient Dp⊥ is obtained from the simulation data. The force dependence of Dp⊥ is found to be a power law in a high force regime, but a constant at low forces. A mean kinetic energy Wp for perpendicular motion is also obtained. The diffusion coefficient is found to increase with Wp with a linear trend at higher energies, but an exponential trend at lower energies.

Planck radiation law and Einstein coefficients reexamined in Kaniadakis κ statistics

Blackbody radiation is reconsidered using the counterpart of the Bose-Einstein distribution in the κ statistics arising from the Kaniadakis entropy. The generalized Planck radiation law is presented and compared to the usual law, to which it reduces in the limiting case κ→0. Effective Einstein's coefficients of emission and absorption are defined in terms of the Kaniadakis parameter κ. It is shown that the Kaniadakis statistics keeps unchanged the first Einstein coefficient A while the second coefficient B admits a generalized form within the present theoretical framework.