Generalized molecular chaos hypothesis and H-theorem: Problem of constraints and amendment of nonextensive statistical mechanics

Quite unexpectedly, kinetic theory is found to specify the correct definition of average value to be employed in nonextensive statistical mechanics. It is shown that the normal average is consistent with the generalized Stosszahlansatz (i.e., molecular chaos hypothesis) and the associated H-theorem, whereas the q-average widely used in the relevant literature is not. In the course of the analysis, the distributions with finite cut-off factors are rigorously treated. Accordingly, the formulation of nonextensive statistical mechanics is amended based on the normal average. In addition, the Shore-Johnson theorem, which supports the use of the q-average, is carefully reexamined, and it is found that one of the axioms may not be appropriate for systems to be treated within the framework of nonextensive statistical mechanics.Comment: 22 pages, no figures. Accepted for publication in Phys. Rev.

Macroscopic proof of the Jarzynski-Wojcik fluctuation theorem for heat exchange

In a recent work, Jarzynski and Wojcik (2004 Phys. Rev. Lett. 92, 230602) have shown by using the properties of Hamiltonian dynamics and a statistical mechanical consideration that, through contact, heat exchange between two systems initially prepared at different temperatures obeys a fluctuation theorem. Here, another proof is presented, in which only macroscopic thermodynamic quantities are employed. The detailed balance condition is found to play an essential role. As a result, the theorem is found to hold under very general conditions.Comment: 9 pages, 0 figure

Liquid-Gas Phase Transition of Supernova Matter and Its Relation to Nucleosynthesis

We investigate the liquid-gas phase transition of dense matter in supernova explosion by the relativistic mean field approach and fragment based statistical model. The boiling temperature is found to be high (T_{boil} >= 0.7 MeV for rho_B >= 10^{-7} fm^{-3}), and adiabatic paths are shown to go across the boundary of coexisting region even with high entropy. This suggests that materials experienced phase transition can be ejected to outside. We calculated fragment mass and isotope distribution around the boiling point. We found that heavy elements at the iron, the first, second, and third peaks of r-process are abundantly formed at rho_B = 10^{-7}, 10^{-5}, 10^{-3} and 10^{-2} fm^{-3}, respectively.Comment: 29 pages, 13 figures. This article is submitted to Nucl. Phys.

Temporal extensivity of Tsallis' entropy and the bound on entropy production rate

The Tsallis entropy, which is a generalization of the Boltzmann-Gibbs entropy, plays a central role in nonextensive statistical mechanics of complex systems. A lot of efforts have recently been made on establishing a dynamical foundation for the Tsallis entropy. They are primarily concerned with nonlinear dynamical systems at the edge of chaos. Here, it is shown by generalizing a formulation of thermostatistics based on time averages recently proposed by Carati [A. Carati, Physica A 348, 110 (2005)] that, whenever relevant, the Tsallis entropy indexed by $q$ is temporally extensive: linear growth in time, i.e., finite entropy production rate. Then, the universal bound on the entropy production rate is shown to be $1/|1-q|$ . The property of the associated probabilistic process, i.e., the sojourn time distribution, determining randomness of motion in phase space is also analyzed.Comment: 25 pages, no figure

Fluctuation theorem for the renormalized entropy change in the strongly nonlinear nonequilibrium regime

Generalizing a recent work [T. Taniguchi and E. G. D. Cohen, J. Stat. Phys. 126, 1 (2006)] that was based on the Onsager-Machlup theory, a nonlinear relaxation process is considered for a macroscopic thermodynamic quantity. It is found that the fluctuation theorem holds in the nonlinear nonequilibrium regime if the change of the entropy characterized by local equilibria is appropriately renormalized. The fluctuation theorem for the ordinary entropy change is recovered in the linear near-equilibrium case. This result suggests a possibility that the the information-theoretic entropy of the Shannon form may be modified in the strongly nonlinear nonequilibrium regime.Comment: 14 pages, no figures. Typos correcte

Stability of Tsallis antropy and instabilities of Renyi and normalized Tsallis entropies: A basis for q-exponential distributions

The q-exponential distributions, which are generalizations of the Zipf-Mandelbrot power-law distribution, are frequently encountered in complex systems at their stationary states. From the viewpoint of the principle of maximum entropy, they can apparently be derived from three different generalized entropies: the Renyi entropy, the Tsallis entropy, and the normalized Tsallis entropy. Accordingly, mere fittings of observed data by the q-exponential distributions do not lead to identification of the correct physical entropy. Here, stabilities of these entropies, i.e., their behaviors under arbitrary small deformation of a distribution, are examined. It is shown that, among the three, the Tsallis entropy is stable and can provide an entropic basis for the q-exponential distributions, whereas the others are unstable and cannot represent any experimentally observable quantities.Comment: 20 pages, no figures, the disappeared "primes" on the distributions are added. Also, Eq. (65) is correcte

Tables of Hyperonic Matter Equation of State for Core-Collapse Supernovae

We present sets of equation of state (EOS) of nuclear matter including hyperons using an SU_f(3) extended relativistic mean field (RMF) model with a wide coverage of density, temperature, and charge fraction for numerical simulations of core collapse supernovae. Coupling constants of Sigma and Xi hyperons with the sigma meson are determined to fit the hyperon potential depths in nuclear matter, U_Sigma(rho_0) ~ +30 MeV and U_Xi(rho_0) ~ -15 MeV, which are suggested from recent analyses of hyperon production reactions. At low densities, the EOS of uniform matter is connected with the EOS by Shen et al., in which formation of finite nuclei is included in the Thomas-Fermi approximation. In the present EOS, the maximum mass of neutron stars decreases from 2.17 M_sun (Ne mu) to 1.63 M_sun (NYe mu) when hyperons are included. In a spherical, adiabatic collapse of a 15$M_\odot$ star by the hydrodynamics without neutrino transfer, hyperon effects are found to be small, since the temperature and density do not reach the region of hyperon mixture, where the hyperon fraction is above 1 % (T > 40 MeV or rho_B > 0.4 fm^{-3}).Comment: 23 pages, 6 figures (Fig.3 and related comments on pion potential are corrected in v3.

Baryon number segregation at the end of the cosmological quark-hadron transition

One of the most interesting questions regarding a possible first order cosmological quark--hadron phase transition concerns the final fate of the baryon number contained within the disconnected quark regions at the end of the transition. We here present a detailed investigation of the hydrodynamical evolution of an evaporating quark drop, using a multi-component fluid description to follow the mechanisms of baryon number segregation. With this approach, we are able to take account of the simultaneous effects of baryon number flux suppression at the phase interface, entropy extraction by means of particles having long mean-free-paths, and baryon number diffusion. A range of computations has been performed to investigate the permitted parameter-space and this has shown that significant baryon number concentrations, perhaps even up to densities above that of nuclear matter, represent an inevitable outcome within this scenario.Comment: 33 pages, Latex file, 6 postscript figures included in the text (psfig.tex). To appear in Phys. Rev. D1

Magnetar Spindown, Hyper-Energetic Supernovae, and Gamma Ray Bursts

The Kelvin-Helmholtz cooling epoch, lasting tens of seconds after the birth of a neutron star in a successful core-collapse supernova, is accompanied by a neutrino-driven wind. For magnetar-strength ($\sim10^{15}$ G) large scale surface magnetic fields, this outflow is magnetically-dominated during the entire cooling epoch.Because the strong magnetic field forces the wind to co-rotate with the protoneutron star,this outflow can significantly effect the neutron star's early angular momentum evolution, as in analogous models of stellar winds (e.g. Weber & Davis 1967). If the rotational energy is large in comparison with the supernova energy and the spindown timescale is short with respect to the time required for the supernova shockwave to traverse the stellar progenitor, the energy extracted may modify the supernova shock dynamics significantly. This effect is capable of producing hyper-energetic supernovae and, in some cases, provides conditions favorable for gamma ray bursts. We estimate spindown timescales for magnetized, rotating protoneutron stars and construct steady-state models of neutrino-magnetocentrifugally driven winds. We find that if magnetars are born rapidly rotating, with initial spin periods ($P$) of $\sim1$ millisecond, that of order $10^{51}-10^{52}$ erg of rotational energy can be extracted in $\sim10$ seconds. If magnetars are born slowly rotating ($P\gtrsim10$ ms) they can spin down to periods of $\sim1$ second on the Kelvin-Helmholtz timescale.Comment: 16 pages, 5 figures, emulateap