Thermodynamics and transport in an active Morse ring chain

We investigate the stochastic dynamics of an one-dimensional ring with N self-driven Brownian particles. In this model neighboring particles interact via conservative Morse potentials. The influence of the surrounding heat bath is modeled by Langevin-forces (white noise) and a constant viscous friction coefficient γ. The Brownian particles are provided with internal energy depots which may lead to active motions of the particles. The depots are realized by an additional nonlinearly velocity-dependent friction coefficient γ ₁(v) in the equations of motions. In the first part of the paper we study the partition functions of time averages and thermodynamical quantities (e.g. pressure) characterizing the stationary physical system. Numerically calculated non-equilibrium phase diagrams are represented. The last part is dedicated to transport phenomena by including a homogeneous external force field that breaks the symmetry of the model. Here we find enhanced mobility of the particles at low temperatures. Received 21 July 2001     

Short-range plasma model for intermediate spectral statistics

We propose a plasma model for spectral statistics displaying level repulsion without long-range spectral rigidity, i.e. statistics intermediate between random matrix and Poisson statistics similar to the ones found numerically at the critical point of the Anderson metal-insulator transition in disordered systems and in certain dynamical systems. The model emerges from Dysons one-dimensional gas corresponding to the eigenvalue distribution of the classical random matrix ensembles by restricting the logarithmic pair interaction to a finite number k of nearest neighbors. We calculate analytically the spacing distributions and the two-level statistics. In particular we show that the number variance has the asymptotic form Σ²(L) ∼χL for large L and the nearest-neighbor distribution decreases exponentially when s→∞, P(s) ∼ exp(- Λs) with Λ = 1/χ = kβ + 1, where β is the inverse temperature of the gas (β = 1, 2 and 4 for the orthogonal, unitary and symplectic symmetry class respectively). In the simplest case of k = β = 1, the model leads to the so-called Semi-Poisson statistics characterized by particular simple correlation functions e.g. P(s) = 4s exp(- 2s). Furthermore we investigate the spectral statistics of several pseudointegrable quantum billiards numerically and compare them to the Semi-Poisson statistics. Received 13 September 2000     

Dynamical effects of a one-dimensional multibarrier potential of finite range

We discuss the properties of a large number N of one-dimensional (bounded) locally periodic potential barriers in a finite interval. We show that the transmission coefficient, the scattering cross section σ, and the resonances of σ depend sensitively upon the ratio of the total spacing to the total barrier width. We also show that a time dependent wave packet passing through the system of potential barriers rapidly spreads and deforms, a criterion suggested by Zaslavsky for chaotic behaviour. Computing the spectrum by imposing (large) periodic boundary conditions we find a Wigner type distribution. We investigate also the S-matrix poles; many resonances occur for certain values of the relative spacing between the barriers in the potential. Received 1st August 2001 and Received in final form 18 November 2001     

The glassy phase of Gallager codes

Gallager codes are the best error-correcting codes to date. In this paper we study them by using the tools of statistical mechanics. The corresponding statistical mechanics model is a spin model on a sparse random graph. The model can be solved by elementary methods (i.e. without replicas) in a large connectivity limit. For low enough temperatures it presents a completely frozen glassy phase (q _EA = 1). The same scenario is shown to hold for finite connectivities. In this case we adopt the replica approach and exhibit a one-step replica symmetry breaking order parameter. We argue that our ansatz yields the exact solution of the model. This allows us to determine the whole phase diagram and to understand the performances of Gallager codes. Received 9 April 2001     

Numerical results for ground states of spin glasses on Bethe lattices

The average ground state energy and entropy for ±J spin glasses on Bethe lattices of connectivities k + 1 = 3..., 26 at T = 0 are approximated numerically. To obtain sufficient accuracy for large system sizes (up to n = 2¹²), the Extremal Optimization heuristic is employed which provides high-quality results not only for the ground state energies per spin e_k+1 but also for their entropies s_k+1. The results indicate sizable differences between lattices of even and odd connectivities. The extrapolated ground state energies compare very well with recent one-step replica symmetry breaking calculations. These energies can be scaled for all even connectivities k + 1 to within a fraction of a percent onto a simple functional form, e _{k + 1} = E _SK - (2E _SK + )/, where E _SK = - 0.7633 is the ground state energy for the broken replica symmetry in the Sherrington-Kirkpatrick model. But this form is in conflict with perturbative calculations at large k + 1, which do not distinguish between even and odd connectivities. We also find non-zero entropies per spin s_k+1 at small connectivities. While s_k+1 seems to vanish asymptotically with 1/(k + 1) for even connectivities, it is numerically indistinguishable from zero already for odd k + 1 ≥ 9. Received 9 August 2002 Published online 27 January 2003 RID="a" ID="a"e-mail: sboettc@emory.edu www.physics.emory.edu/faculty/boettcher     

Observation of Δ+→ pπ0 decay in heavy-ion collisions

Proton-π⁰ coincidences have been measured at the beam energy of 180A MeV in the reaction Ar+Ca studied by TAPS at SIS/GSI. In the proton-π⁰ invariant mass spectrum we observe a significant excess of strength above the background obtained by event mixing. We attribute this signal to the strength distribution N _Δ of the Δ baryonic resonance. No correlation is observed in the case of deuteron-π⁰ coincidences. Assuming isotropic emission of π⁰ and Δ⁺ from a midrapidity thermal source and isospin symmetry, we determined the global N _Δ/N _πratio of 0.79 ±0.30(stat) ±0.2(syst). This value indicates that most pions produced at subthreshold energy in heavy-ion reaction are mediated by the Δ-resonance. Received: 7 April 2000 / Accepted: 29 August 2000     

Spin glasses on Bethe lattices for large coordination number

We study spin glasses on random lattices with finite connectivity. In the infinite connectivity limit they reduce to the Sherrington Kirkpatrick model. In this paper we investigate the expansion around the high connectivity limit. Within the replica symmetry breaking scheme at two steps, we compute the free energy at the first order in the expansion in inverse powers of the average connectivity (z), both for the fixed connectivity and for the fluctuating connectivity random lattices. It is well known that the coefficient of the 1/z correction for the free energy is divergent at low temperatures if computed in the one step approximation. We find that this annoying divergence becomes much smaller if computed in the framework of the more accurate two steps breaking. Comparing the temperature dependance of the coefficients of this divergence in the replica symmetric, one step and two steps replica symmetry breaking, we conclude that this divergence is an artefact due to the use of a finite number of steps of replica symmetry breaking. The 1/z expansion is well defined also in the zero temperature limit. Received 15 July 2002 Published online 31 December 2002     

SU(3) coupling schemes for odd-odd nuclei in the interacting boson – fermion – fermion model with both odd proton and odd neutron in natural parity orbits

The SU^BFF(3) dynamical symmetry limits of interacting boson – fermion – fermion model are identified and they are appropriate for heavy deformed odd – odd nuclei for configurations with both the odd proton and odd neutron occupying all the natural parity orbits in the corresponding valence shells. There are three symmetry limits and their correspondence with two quasi-particle (proton-neutron) Nilsson configurations is established; one of the limits mixes both Nilsson n_z's and Λ's and other two limits mix only Nilsson Λ's. The ¹⁹¹Ir (d,t) ¹⁹⁰Ir single nucleon transfer spectroscopic strengths are well described by one of the symmetry limits that mixes only Nilsson Λ's. Received: 22 June 1998     

Phase equilibria in random multiblock copolymers

A mean-field theory for domain structures in random multiblock copolymer melts is developed. We focus on the finite molecular weight effects resulting in a competition between macroscopic phase separation and microdomain formation in the system. We identify an essential parameter N ε controlling the phase behavior of the system, where N is the number of blocks per chain and ε is the composition asymmetry parameter (= the difference between the mean copolymer composition and its critical value). The phase diagram involving N ε and the reduced temperature as variables is obtained. The regions of coexistence of two or more phases are identified. We show that a superstructure formation on cooling is always pre-empted by a macroscopic phase separation of the macroscopically homogeneous (disordered) system yielding two homogeneous phases: H ₀↦H ₁ + H ₂. The third (lamellar) phase separates on further cooling. Then hexagonal and body-centred-cubic phases take over if N ε 1. As the Flory interaction parameter χ increases further, the standard transitions BCC↦HEX↦LAM take place. Received 13 July 2001     

Dynamical properties of the slithering-snake algorithm: A numerical test of the activated-reptation hypothesis

Correlations in the motion of reptating polymers in a melt are investigated by means of Monte Carlo simulations of the three-dimensional slithering-snake version of the bond-fluctuation model. Surprisingly, the slithering-snake dynamics becomes inconsistent with classical reptation predictions at high chain overlap (created either by chain length N or by the volume fraction φ of occupied lattice sites), where the relaxation times increase much faster than expected. This is due to the anomalous curvilinear diffusion in a finite time window whose upper bound (N) is set by the density of chain ends φ/N. Density fluctuations created by passing chain ends allow a reference polymer to break out of the local cage of immobile obstacles created by neighboring chains. The dynamics of dense solutions of “snakes” at t ≪ is identical to that of a benchmark system where all chains but one are frozen. We demonstrate that the subdiffusive dynamical regime is caused by the slow creeping of a chain out of its correlation hole. Our results are in good qualitative agreement with the activated-reptation scheme proposed recently by Semenov and Rubinstein (Eur. Phys. J. B, 1 (1998) 87). Additionally, we briefly comment on the relevance of local relaxation pathways within a slithering-snake scheme. Our preliminary results suggest that a judicious choice of the ratio of local to slithering-snake moves is crucial to equilibrate a melt of long chains efficiently. Received: 18 December 2002 / Accepted: 3 April 2003 / Published online: 12 May 2003 RID="a" ID="a"e-mail: jwittmer@dpm.univ-lyon1.fr RID="b" ID="b"Current address: University of Illinois at Urbana-Champaign.     

Order-parameter fluctuations (OPF) in spin glasses: Monte Carlo simulations and exact results for small sizes

The use of parameters measuring order-parameter fluctuations (OPF) has been encouraged by the recent results reported in referenece [2,3] which show that two of these parameters, G and G _c, take universal values in the . In this paper we present a detailed study of parameters measuring OPF for two mean-field models with and without time-reversal symmetry which exhibit different patterns of replica symmetry breaking below the transition: the Sherrington-Kirkpatrick model with and without a field and the Ising p-spin glass (p = 3). We give numerical results and analyze the consequences which replica equivalence imposes on these models in the infinite volume limit. We give evidence for the transition in each system and discuss the character of finite-size effects. Furthermore, a comparative study between this new family of parameters and the usual Binder cumulant analysis shows what kind of new information can be extracted from the finite T behavior of these quantities. The two main outcomes of this work are: 1) Parameters measuring OPF give better estimates than the Binder cumulant for T _c and even for very small systems they give evidence for the transition. 2) For systems with no time-reversal symmetry, parameters defined in terms of connected quantities are the proper ones to look at. Received 20 September 2000 and Received in final form 10 January 2001     

πNN coupling determined beyond the chiral limit

Within the conventional QCD sum rules, we calculate the πNN coupling constant, g _πN, beyond the chiral limit using two-point correlation function with a pion. For this purpose, we consider the Dirac structure, iγ₅, at m _π ² order in the expansion of the correlator in terms of the pion momentum. For a consistent treatment of the sum rule, we include the linear terms in quark mass as they constitute the same chiral order as m _π ². In this sum rule, we obtain g _πN= 13.3 ± 1.2, which is very close to the empirical πNN coupling. This demonstrates that going beyond the chiral limit is crucial in determining the coupling. Received: 8 July 1999 / Revised version: 20 August 1999     

Statistical physics analysis of the computational complexity of solving random satisfiability problems using backtrack algorithms

The computational complexity of solving random 3-Satisfiability (3-SAT) problems is investigated using statistical physics concepts and techniques related to phase transitions, growth processes and (real-space) renormalization flows. 3-SAT is a representative example of hard computational tasks; it consists in knowing whether a set of αN randomly drawn logical constraints involving N Boolean variables can be satisfied altogether or not. Widely used solving procedures, as the Davis-Putnam-Loveland-Logemann (DPLL) algorithm, perform a systematic search for a solution, through a sequence of trials and errors represented by a search tree. The size of the search tree accounts for the computational complexity, i.e. the amount of computational efforts, required to achieve resolution. In the present study, we identify, using theory and numerical experiments, easy (size of the search tree scaling polynomially with N) and hard (exponential scaling) regimes as a function of the ratio α of constraints per variable. The typical complexity is explicitly calculated in the different regimes, in very good agreement with numerical simulations. Our theoretical approach is based on the analysis of the growth of the branches in the search tree under the operation of DPLL. On each branch, the initial 3-SAT problem is dynamically turned into a more generic 2+p-SAT problem, where p and 1 - p are the fractions of constraints involving three and two variables respectively. The growth of each branch is monitored by the dynamical evolution of α and p and is represented by a trajectory in the static phase diagram of the random 2+p-SAT problem. Depending on whether or not the trajectories cross the boundary between satisfiable and unsatisfiable phases, single branches or full trees are generated by DPLL, resulting in easy or hard resolutions. Our picture for the origin of complexity can be applied to other computational problems solved by branch and bound algorithms. Received 10 March 2001     

Polymer chain in a quenched random medium: slow dynamics and ergodicity breaking

The Langevin dynamics of a self-interacting chain embedded in a quenched random medium is investigated by making use of the generating functional method and one-loop (Hartree) approximation. We have shown how this intrinsic disorder causes different dynamical regimes. Namely, within the Rouse characteristic time interval the anomalous diffusion shows up. The corresponding subdiffusional dynamical exponents have been explicitly calculated and thoroughly discussed. For the larger time interval the disorder drives the center of mass of the chain to a trap or frozen state provided that the Harris parameter, (Δ/b ^d)N ^{2 - νd}≥1, where Δ is a disorder strength, b is a Kuhnian segment length, N is a chain length and ν is the Flory exponent. We have derived the general equation for the non-ergodicity function f (p) which characterizes the amplitude of frozen Rouse modes with an index p = 2πj/N. The numerical solution of this equation has been implemented and shown that the different Rouse modes freeze up at the same critical disorder strength Δ _c ∼ N ^{- γ} where the exponent γ ≈ 0.25 and does not depend from the solvent quality. Received 17 December 2002 Published online 23 May 2003 RID="a" ID="a"e-mail: vilgis@mpip-mainz.mpg.de     

ωN final state interactions and ω-meson production from heavy-ion collisions

We calculate the elastic and inelastic ωN→ωN, →πN, →ρN, →ρπN, →ππN, →σN reactions within a boson exchange approximation where the ωρπ coupling constant and form factor are fixed by the reaction πN→ωN in comparison to the experimental data. We find rather large ωN cross sections at low relative momenta of the ω-meson which leads to a substantial broadening of the ω-meson width in nuclear matter. The implications of the ωN final state interactions are studied for ω production in ¹²C +¹²C, ⁴⁰Ca +⁴⁰Ca and ⁵⁸Ni +⁵⁸Ni reactions at about 2 · A GeV within the HSD transport approach; the drastic changes of the transverse mass spectra relative to a general m _T-scaling (for π⁰ and η mesons) might be controlled experimentally by the TAPS Collaboration. Received: 28 April 1999 / Revised version: 7 June 1999     

The Nambu–Jona–Lasinio model of light nuclei

The Nambu–Jona–Lasinio model of the deuteron suggested by Nambu and Jona–Lasinio (Phys. Rev. 124 (1961) 246) is formulated from the first principles of QCD. The deuteron appears as a neutron–proton collective excitation, i.e. a Cooper np–pair, induced by a phenomenological local four–nucleon interaction in the nuclear phase of QCD. The model describes the deuteron coupled to itself, nucleons and other particles through one–nucleon loop exchanges providing a minimal transfer of nucleon flavours from initial to final nuclear states and accounting for contributions of nucleon–loop anomalies which are completely determined by one–nucleon loop diagrams. The dominance of contributions of nucleon–loop anomalies to effective Lagrangians of low–energy nuclear interactions is justified in the large N _C expansion, where N _C is the number of quark colours. Received: 10 March 2000     

Mixing and sorting of bidisperse two-dimensional bubbles

We have examined a number of candidates for the minimum-surface-energy arrangement of two-dimensional clusters composed of N bubbles of area 1 and N bubbles of area λ ( λ≤1). These include hexagonal bubbles sorted into two monodisperse honeycomb tilings, and various mixed periodic tilings with at most four bubbles per unit cell. We identify, as a function of λ, the minimal configuration for N → ∞. For finite N, the energy of the external (i.e., cluster-gas) boundary and that of the interface between honeycombs in “phase-separated” clusters have to be taken into account. We estimate these contributions and find the lowest total energy configuration for each pair (N,λ). As λ is varied, this alternates between a circular cluster of one of the mixed tilings, and “partial wetting” of the monodisperse honeycomb of bubble area 1 by the monodisperse honeycomb of bubble area λ. Received 1 August 2002 RID="a" ID="a"e-mail: paulo@ist.utl.pt     

On N-wave type systems and their gauge equivalent

The class of nonlinear evolution equations (NLEE) - gauge equivalent to the N-wave equations related to the simple Lie algebra are derived and analyzed. They are written in terms of (x, t) ∈ satisfying r = rank nonlinear constraints. The corresponding Lax pairs and the time evolution of the scattering data are found. The Zakharov-Shabat dressing method is appropriately modified to construct their soliton solutions. Received 20 October 2001 / Received in final form 30 April 2002 Published online 2 October 2002 RID="a" ID="a"e-mail: gerjikov@inrne.bas.bg     

Microscopic three-body force for asymmetric nuclear matter

Brueckner calculations including a microscopic three-body force have been extended to isospin-asymmetric nuclear matter. The effects of the three-body force on the equation of state and on the single-particle properties of nuclear matter are discussed with a view to possible applications in nuclear physics and astrophysics. It is shown that, even in the presence of the three-body force, the empirical parabolic law of the energy per nucleon vs. isospin asymmetry β = (N - Z)/A is fulfilled in the whole asymmetry range 0≤β≤1 up to high densities. The three-body force provides a strong enhancement of the symmetry energy which increases with density in good agreement with the predictions of relativistic approaches. The Lane's assumption that proton and neutron mean fields linearly vary vs. the isospin parameter is violated at high density due to the three-body force, while the momentum dependence of the mean fields turns out to be only weakly affected. Consequently, a linear isospin split of the neutron and proton effective masses is found for both cases with and without the three-body force. The isospin effects on multifragmentation events and collective flows in heavy-ion collisions are briefly discussed along with the conditions for direct URCA processes to occur in the neutron star cooling. Received: 18 February 2002 / Accepted: 16 May 2002     

Study of N*(1440) from J/Ψ decays

For J/Ψ↦ pπ⁰ and pπ⁺π^-, the π⁰ p and pπ⁺π^- systems are limited to be pure isospin-(1/2) due to isospin conservation. This is a big advantage in studying N^* resonances from J/Ψ decays, compared with πN and γN experiments. The process J/Ψ↦ N ^* or p provides a new way to probe the internal structure of the N^* resonances. Here we report a quark model calculation for J/Ψ↦ p, N ^*(1440) and N ^*. The implication for the internal structure of N ^*(1440) is discussed. Received: 1 June 2001 / Accepted: 20 June 2001