Standing wave instabilities, breather formation and thermalization in a Hamiltonian anharmonic lattice

Modulational instability of travelling plane waves is often considered as the first step in the formation of intrinsically localized modes (discrete breathers) in anharmonic lattices. Here, we consider an alternative mechanism for breather formation, originating in oscillatory instabilities of spatially periodic or quasiperiodic nonlinear standing waves (SWs). These SWs are constructed for Klein-Gordon or Discrete Nonlinear Schr?dinger lattices as exact time periodic and time reversible multibreather solutions from the limit of uncoupled oscillators, and merge into harmonic SWs in the small-amplitude limit. Approaching the linear limit, all SWs with nontrivial wave vectors (0 < Q < π) become unstable through oscillatory instabilities, persisting for arbitrarily small amplitudes in infinite lattices. The dynamics resulting from these instabilities is found to be qualitatively different for wave vectors smaller than or larger than π/2, respectively. In one regime persisting breathers are found, while in the other regime the system thermalizes. Received 6 October 2001 / Received in final form 1st March 2002 Published online 2 October 2002 RID="a" ID="a"e-mail: mjn@ifm.liu.se     

Birth and long-time stabilization of out-of-equilibrium coherent structures

We study an analytically tractable model with long-range interactions for which an out-of-equilibrium very long-lived coherent structure spontaneously appears. The dynamics of this model is indeed very peculiar: a bicluster forms at low energy and is stable for very long time, contrary to statistical mechanics predictions. We first explain the onset of the structure, by approximating the short time dynamics with a forced Burgers equation. The emergence of the bicluster is the signature of the shock waves present in the associated hydrodynamical equations. The striking quantitative agreement with the dynamics of the particles fully confirms this procedure. We then show that a very fast timescale can be singled out from a slower motion. This enables us to use an adiabatic approximation to derive an effective Hamiltonian that describes very well the long time dynamics. We then get an explanation of the very long time stability of the bicluster: this out-of-equilibrium state corresponds to a statistical equilibrium of an effective mean-field dynamics. Received 28 February 2002 / Received in final form 24 July 2002 Published online 31 October 2002 RID="a" ID="a"e-mail: Thierry.Dauxois@ens-lyon.fr RID="b" ID="b"UMR-CNRS 5672 RID="c" ID="c"UMR 5582     

Dynamics of twists on antiferromagnetic spin chains: Theory

The equation of motion of twists on classical antiferromagnetic Heisenberg spin chains are derived. It is shown that twists interact via position- and velocity-dependent long-range two-body forces. A quiescent regime is identified wherein the system conserves momentum. With increasing kinetic energy the system exits this regime and momentum conservation is violated due to walls annihilation. A bitwist system is shown to be integrable and its exact solution is analysed. Many-twist systems are discussed and novel periodic twist lattice solutions are found on closed chains. The stability of these solutions is discussed. Received 12 June 2002 Published online 2 October 2002 RID="a" ID="a"e-mail: rbbll@phy.cam.ac.uk     

Isomerization and dissociation of small (CH <Emphasis Type="Bold">3</Emphasis>CN) <Emphasis Type="Bold">n</Emphasis> molecular clusters: a computational study

The isomerization and evaporation processes in the neutral homogeneous (CH₃CN)_n molecular clusters (n = 2-7) have been investigated using classical molecular dynamics simulations. The evaporation rate constants and the kinetic energy release in the dissociation have been analysed as a function of the cluster size and as a function of the internal energy in the parent cluster. The competition between monomer and dimer ejections has been also carefully studied. All the dynamical properties in these dissociative processes have been discussed in relation to the static properties of the clusters involved in the dissociation and also in relation to the solid-liquid like transition which appears in these homogeneous molecular clusters. Received 19 November 2002 / Received in final form 5 February 2003 Published online 29 April 2003 RID="a" ID="a"e-mail: pascal.parneix@ppm.u-psud.fr RID="b" ID="b"Laboratoire associé à l'université Paris-Sud.     

DFT studies of ionic vibrations in Na clusters

We use time-dependent density functional theory coupled to molecular dynamics for ionic motion to compute the spectra of ionic vibrations in small Na clusters. Comparison with results from the distance dependent tight-binding approach shows good agreement between these two very different methods. We discuss the evolution of the spectra with cluster size and charge and the impact of ionic vibrations on the optical response. Received 23 July 2001 / Received in final form 5 July 2002 Published online 8 October 2002 RID="a" ID="a"e-mail: suraud@irsamc.ups-tlse.fr     

Resonance interaction between one excited and one ground state atom

The nature of the resonance interaction between two isotropic atoms in an excited configuration is reinvestigated. The currently accepted oscillatory form for the long-range retarded resonance interaction is shown to be a subtle artefact that arises due to too drastic approximations. Formulation of the resonance interaction energy problem in terms of the interacting system leads to a form that it is ∝ r ^-4 in the retarded limit. We also demonstrate that the resonance interaction energy at any finite temperature goes over to purely classical long-range asymptote. This manifestation of the correspondence principle is due to thermal excitation of the electromagnetic field. We finally discuss why the textbook result for the F?rster energy transfer between two atoms is incorrect for the same reasons. Received 31 January 2002 / Received in final form 25 July 2002 Published online 29 October 2002 RID="a" ID="a"e-mail: mtb110@rsphysse.anu.edu.au     

Folding of lattice protein chains with modified Gō potential

We propose a modified Gō model in which the pairwise interaction energies vary as local environment changes. The stability difference between the surface and the core is also well considered in this model. Thermodynamic and kinetic studies suggest that this model has improved folding cooperativity and foldability in contrast with the Gō model. The free energy landscape of this model has broad barriers and narrow denatured states, which is consistent with that of the two-state folding proteins and is lacked for the Gō model. The role of non-native interactions in protein folding is also studied. We find that appropriate consideration of the contribution of the non-native interactions may increase the folding rate around the transition temperature. Our results show that conformation-dependent interaction between the residues is a realistic representation of potential functions in protein folding. Received 10 April 2002 / Received in final form 20 August 2002 Published online 19 December 2002 RID="a" ID="a"e-mail: wangwei@nju.edu.cn     

Quasi-periodic and periodic solutions for dynamical systems related to Korteweg-de Vries equation

We consider quasi-periodic and periodic (cnoidal) wave solutions of a set of n-component dynamical systems related to Korteweg-de Vries equation. Quasi-periodic wave solutions for these systems are expressed in terms of Novikov polynomials. Periodic solutions in terms of Hermite polynomials and generalized Hermite polynomials for dynamical systems related to Korteweg-de Vries equation are found. Received 15 October 2001 / Received in final form 6 March 2002 Published online 2 October 2002 RID="a" ID="a"e-mail: nakostov@ie.bas.bg     

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     

The $-game

We propose a payoff function extending Minority Games (MG) that captures the competition between agents to make money. In contrast with previous MG, the best strategies are not always targeting the minority but are shifting opportunistically between the minority and the majority. The emergent properties of the price dynamics and of the wealth of agents are strikingly different from those found in MG. As the memory of agents is increased, we find a phase transition between a self-sustained speculative phase in which a “stubborn majority” of agents effectively collaborate to arbitrage a market-maker for their mutual benefit and a phase where the market-maker always arbitrages the agents. A subset of agents exhibit a sustained non-equilibrium risk-return profile. Received 5 June 2002 / Received in final form 21 November 2002 Published online 27 January 2003 RID="a" ID="a"e-mail: sornette@unice.fr RID="b" ID="b"CNRS UMR7536 RID="c" ID="c"CNRS UMR6622     

Metastable states of the Ising chain with Kawasaki dynamics

We consider a ferromagnetic Ising chain evolving under Kawasaki dynamics at zero temperature. We investigate the statistics of the blocking time, as well as various characteristics of the metastable configurations reached by the system, including the statistics of the final energy, the spin correlations, and the distribution of domain sizes. Results of extensive numerical simulations are compared with analytical predictions made for the a priori ensemble of all blocked configurations with equal weights. Qualitative differences are found, e.g. in the domain sizes, which are found to be neither statistically independent nor exponentially distributed. Received 24 October 2002 / Received in final form 13 January 2003 Published online 1st April 2003 RID="a" ID="a"e-mail: luck@spht.saclay.cea.fr RID="b" ID="b"URA 2306 of CNRS     

Three-dimensional motion of a vortex filament and its relation to the localized induction hierarchy

Three-dimensional motion of a slender vortex tube, embedded in an inviscid incompressible fluid, is investigated under the localized induction approximation for the Euler equations. Using the method of matched asymptotic expansions in a small parameter ε, the ratio of core radius to curvature radius, the velocity of a vortex filament is derived to O(ε³), whereby the influence of elliptical deformation of the core due to the self-induced strain is taken into account. It is found that there is an integrable line in the core whose evolution obeys a summation of the first and third terms of the localized induction hierarchy. Received 2 October 2001 / Received in final form 10 May 2002 Published online 2 October 2002 RID="a" ID="a"e-mail: yasuhide@math.kyushu-u.ac.jp     

Gamma-ray tracking: Utilizing new concepts in the detection of gamma-radiation

Gamma-ray tracking in a closed array of highly segmented HPGe detectors is a new concept for the detection of γ-radiation. Each of the interacting γ-rays is identified and separated by measuring the energies and positions of individual interactions and by applying tracking algorithms to reconstruct the scattering sequences, even if many γ-rays hit the array at the same time. The three-dimensional position and the energy of interactions are determined by using two-dimensionally segmented Ge detectors along with pulse-shape analysis of the signals. Such a detector will have new and much improved capabilities compared to current γ-ray spectrometer. One implementation of this concept, called GRETA (Gamma-Ray Energy Tracking Array), is currently being under development at LBNL. Received: 21 March 2002 / Accepted: 16 May 2002 / Published online: 31 October 2002 RID="a" ID="a"e-mail: kvetter@lbl.gov     

Vertex cover problem studied by cavity method: Analytics and population dynamics

We study the vertex cover problem on finite connectivity random graphs by zero-temperature cavity method. The minimum vertex cover corresponds to the ground state(s) of a proposed Ising spin model. When the connectivity c > e = 2.718282, there is no state for this system as the reweighting parameter y, which takes a similar role as the inverse temperature β in conventional statistical physics, approaches infinity; consequently the ground state energy is obtained at a finite value of y when the free energy function attains its maximum value. The minimum vertex cover size at given c is estimated using population dynamics and compared with known rigorous bounds and numerical results. The backbone size is also calculated. Received 11 November 2002 Published online 1st April 2003 RID="a" ID="a"e-mail: zhou@mpikg-golm.mpg.de     

K<Superscript>+</Superscript> production in proton-nucleus reactions and the role of momentum-dependent potentials

The production of K⁺-mesons in proton-nucleus collisions from 1.0 to 2.5GeV is analyzed with respect to one-step nucleon-nucleon ( NN → NYK ⁺) and two-step Δ-nucleon ( ΔN → K ⁺ YN) or pion-nucleon ( πN → K ⁺ Y) production channels on the basis of a coupled-channel transport approach (CBUU) including the kaon final-state interactions. The influence of momentum-dependent potentials for the nucleon, hyperon and kaon in the final state are studied as well as the importance of K⁺ elastic rescattering in the target nucleus. The transport calculations are compared to the experimental K⁺ spectra taken at LBL Berkeley, SATURNE, CELSIUS, GSI and COSY-Jülich. It is found that the momentum-dependent baryon potentials affect the excitation function of the K⁺ cross-section; at low bombarding energies of ∼ 1.0GeV the attractive baryon potentials in the final state lead to a relative enhancement of the kaon yield, whereas the net repulsive potential at bombarding energies ∼ 2GeV causes a decrease of the K⁺ cross-section. Furthermore, it is pointed out that especially the K⁺ spectra at low momenta (or kinetic energy T _K) allow to determine the in-medium K⁺ potential almost model independently due to a relative shift of the K⁺ spectra in kinetic energy that arises from the acceleration of the kaons when propagating out of the nuclear medium to free space, i.e. converting the potential energy to the kinetic energy of the free kaon. Received: 28 January 2002 / Accepted: 3 June 2002 / Published online: 19 November 2002 RID="a" ID="a"e-mail: Wolfgang.Cassing@theo.physik.uni-giessen.de Communicated by P. Schuck     

Foam imbibition in microgravity

We report an experimental study of aqueous foam imbibition in microgravity with strict mass conservation. The foam is in a Hele-Shaw cell. The bubble edge width ℓ is measured by image analysis. The penetration of the liquid in the foam, the foam imbibition, the foam inflation, and the rigidity loss are shown all to obey strict diffusion processes. The motion of bubbles needed for the foam inflation is a slow two-dimensional process with respect to the one-dimensional capillary rise of liquid. The foam is found to imbibes faster than it inflates. Received 20 May 2002 / Received in final form 21 January 2003 Published online 23 May 2003 RID="a" ID="a"e-mail: herve.caps@ulg.ac.be     

Selective advantage of topological disorder in biological evolution

We examine a model of biological evolution of Eigen's quasispecies in a so-called holey fitness landscape, where the fitness of a site is either 0 (lethal site) or a uniform positive constant (viable site). The evolution dynamics is therefore determined by the topology of the genome space which is modelled by the random Bethe lattice. We use the effective medium and single-defect approximations to find the criteria under which the localized quasispecies cloud is created. We find that shorter genomes, which are more robust to random mutations than average, represent a selective advantage which we call “topological”. A way of assessing empirically the relative importance of reproductive success and topological advantage is suggested. Received 9 August 2002 / Received in final form 7 November 2002 Published online 14 February 2003 RID="a" ID="a"e-mail: slanina@fzu.cz     

Statics and dynamics of phase segregation in multicomponent fermion gas

We investigate the statics and dynamics of spatial phase segregation process of a mixture of fermion atoms in a harmonic trap using the density functional theory. The kinetic energy of the fermion gas is written in terms of the density and its gradients. Several cases have been studied by neglecting the gradient terms (the Thomas-Fermi limit) which are then compared with the Monte-Carlo results using the full gradient corrected kinetic energy. A linear instability analysis has been performed using the random-phase approximation. Near the onset of instability, the fastest unstable mode for spinodal decomposition is found to occur at q = 0. However, in the strong coupling limit, many more modes with q ≈ K _F decay with comparable time scales. Received 12 September 2001 Published online 24 September 2002 RID="a" ID="a"e-mail: k1@sharif.edu RID="b" ID="b"Permanent address: Dr. Vijay Kumar Foundation, 45 Bazaar Street, K.K. Nagar (West), Chennai 600 078, India.     

Cellular oscillator with a small number of particles

To illustrate complex spatio-temporal effects which are possible in cellular reactions with a small number of particles, we present Monte Carlo simulations of the formation of oscillatory spark-like patterns in a model completely stochastic Ca⁺² oscillator. Our analysis shows that in order to observe such patterns the minimum average number of Ca⁺² ions in the cytosol may be as low as about 50. Received 21 April 2002 / Received in final form 19 June 2002 Published online 14 October 2002 RID="a" ID="a"e-mail: zhdanov@catalysis.nsk.su