Discrete breathers in a model with Morse potentials

Under harmonic approximation, this paper discusses the linear dispersion relation of the one-dimensional chain. The existence and evolution of discrete breathers in a general one-dimensional chain are analysed for two particular examples of soft (Morse) and hard (quartic) on-site potentials. The existence of discrete breathers in one-dimensional and two-dimensional Morse lattices is proved by using rotating wave approximation, local anharmonic approximation and a numerical method. The localization and amplitude of discrete breathers in the two-dimensional Morse lattice with on-site harmonic potentials correlate closely to the Morse parameter a and the on-site parameter к.     

Influence of momentum-dependent interactions on balance energy and mass dependence⋆

We aim to understand the role of momentum-dependent interactions in transverse flow as well as in its disappearance. For the present study, central collisions involving masses between 24 and 394 are considered. We find that the momentum-dependent interactions have different impact in lighter colliding nuclei compared to heavier colliding nuclei. In lighter nuclei, the contribution of the mean field towards flow is smaller compared to heavier nuclei where binary nucleon-nucleon collisions dominate the scene. The inclusion of momentum-dependent interactions also explains the energy of the vanishing flow in the ¹²C + ¹²C reaction which otherwise was not possible with the static hard equation of state. An excellent agreement of our theoretical attempt is found for balance energy with experimental data throughout the periodic table.     

q-Breathers and the Fermi-Pasta-Ulam problem

The Fermi-Pasta-Ulam (FPU) paradox consists of the non-equipartition of energy among normal modes of a weakly anharmonic atomic chain model. In the harmonic limit each normal mode corresponds to a periodic orbit in phase space and is characterized by its wave number q. We continue normal modes from the harmonic limit into the FPU parameter regime and obtain persistence of these periodic orbits, termed here q-breathers (QB). They are characterized by time periodicity, exponential localization in the q-space of normal modes and linear stability up to a size-dependent threshold amplitude. Trajectories computed in the original FPU setting are perturbations around these exact QB solutions. The QB concept is applicable to other nonlinear lattices as well.     

Thermodynamic Limit for Dipolar Media

We prove existence of a shape- and boundary-condition-independent thermodynamic limit for fluids and solids of identical particles with electric or magnetic dipole moments. Our result applies to fluids of hard-core particles, to dipolar soft spheres and Stockmayer fluids, to disordered solid composites, and to regular crystal lattices. In addition to their permanent dipole moments, particles may further polarize each other. Classical and quantum models are treated. Shape independence depends on the reduction in free energy accomplished by domain formation, so our proof applies only in the case of zero applied field. Existence of a thermodynamic limit implies texture formation in spontaneously magnetized liquids and disordered solids analogous to domain formation in crystalline solids.     

Eigenmodes and thermodynamics of a Coulomb chain in a harmonic potential

The density of ions trapped in a harmonic potential in one dimension is not uniform. Consequently the eigenmodes are not phononlike waves. We calculate the long-wavelength modes in the continuum limit, and evaluate the density of states in the short-wavelength limit for chains of N>1 ions. Remarkably, the results that are found analytically in the thermodynamic limit provide a good estimate of the spectrum of excitations of small chains down to few tens of ions. The spectra are used to compute the thermodynamic functions of the chain. Deviations from the extensivity of the thermodynamic quantities are found. An analytic expression for the critical transverse frequency determining the stability of a linear chain is derived.     

Energy Transport in Weakly Anharmonic Chains

We investigate the energy transport in a one-dimensional lattice of oscillators with a harmonic nearest neighbor coupling and a harmonic plus quartic on-site potential. As numerically observed for particular coupling parameters before, and confirmed by our study, such chains satisfy Fourier’s law: a chain of length N coupled to thermal reservoirs at both ends has an average steady state energy current proportional to 1/N. On the theoretical level we employ the Peierls transport equation for phonons and note that beyond a mere exchange of labels it admits nondegenerate phonon collisions. These collisions are responsible for a finite heat conductivity. The predictions of kinetic theory are compared with molecular dynamics simulations. In the range of weak anharmonicity, respectively low temperatures, reasonable agreement is observed.     

Soft multihadron production from partonic structure and fragmentation functions

We describe a two-chain model for soft multiparticle production in hadronic collisions. The model is formulated in a parton framework and is consistent with the dual topological scheme for the Pomeron. The sole inputs are valence quark structure functions in the colliding hadrons and parton fragmentation functions, both of which are known from “hard” processes. Our model, which contains no adjustable parameters, reproduces the shape, the energy dependence, and the normalization of inclusive spectra both in the central region and in the fragmentation region. The model provides a natural explanation for the ratio of πp topp cross sections.     

Phonon Excitation and Energy Redistribution in Phonon Space for Energy Dissipation and Transport in Lattice Structure with Nonlinear Dispersion

We first propose fundamental solutions of wave propagation in dispersive chain subject to a localized initial perturbation in the displacement. Analytical solutions are obtained for both second order nonlinear dispersive chain and homogenous harmonic chain using stationary phase approximation. Solution is also compared with numerical results from molecular dynamics(MD) simulations. Locally dominant phonon modes(k-space) are introduced based on these solutions. These locally defined spatially and temporally varying phonon modes k(x, t) are critical to the concept of the local thermodynamic equilibrium(LTE). Wave propagation accompanying with the nonequilibrium dynamics leads to the excitation of these locally defined phonon modes. It is found that the system energy is gradually redistributed among these excited phonons modes(k-space). This redistribution process is only possible with nonlinear dispersion and requires a finite amount of time to achieve a steady state distribution. This time scale is dependent on the spatial distribution(or frequency content) of the initial perturbation and the dispersion relation. Sharper and more concentrated perturbation leads to a faster energy redistribution and dissipation. This energy redistribution generates localized phonons with various frequencies that can be important for phonon-phonon interaction and energy dissipation in nonlinear systems.Depending on the initial perturbation and temperature, the time scale associated with this energy distribution can be critical for energy dissipation compared to the Umklapp scattering process. Ballistic type of heat transport along the harmonic chain reveals that at any given position, the lowest mode(k = 0) is excited first and gradually expanding to the highest mode(kmax(x, t)), where kmax(x, t) can only asymptotically approach the maximum mode kBof the first Brillouin zone(kmax(x, t) → kB). No energy distributed into modes with kmax(x, t) k kBdemonstrates that the local thermodynamic equilibrium cannot be established in harmonic chain. Energy is shown to be uniformly distributed in all available phonon modes k ≤ kmax(x, t) at any position with heat transfer along the harmonic chain. The energy flux along the chain is shown to be a constant with time and proportional to the sound speed(ballistic transport).Comparison with the Fourier's law leads to a time-dependent thermal conductivity that diverges with time.     

Rate of equilibration of a one-dimensional Wigner crystal

We consider a system of one-dimensional spinless particles interacting via long-range repulsion. In the limit of strong interactions the system is a Wigner crystal, with excitations analogous to phonons in solids. In a harmonic crystal the phonons do not interact, and the system never reaches thermal equilibrium. We account for the anharmonism of the Wigner crystal and find the rate at which it approaches equilibrium. The full equilibration of the system requires umklapp scattering of phonons, resulting in exponential suppression of the equilibration rate at low temperatures.     

Transverse-momentum spectra in Au+Au and d+Au collisions at sqrt[s(NN)]=200 GeV and the pseudorapidity dependence of high-p(T) suppression

We present spectra of charged hadrons from Au+Au and d+Au collisions at sqrt[s(NN)]=200 GeV measured with the BRAHMS experiment at RHIC. The spectra for different collision centralities are compared to spectra from p+(-)p collisions at the same energy scaled by the number of binary collisions. The resulting ratios (nuclear modification factors) for central Au+Au collisions at eta=0 and eta=2.2 evidence a strong suppression in the high p(T) region (>2 GeV/c). In contrast, the d+Au nuclear modification factor (at eta=0) exhibits an enhancement of the high p(T) yields. These measurements indicate a high energy loss of the high p(T) particles in the medium created in the central Au+Au collisions. The lack of suppression in d+Au collisions makes it unlikely that initial state effects can explain the suppression in the central Au+Au collisions.     

Correlations in ballistic processes

We investigate a class of reaction processes in which particles move ballistically and react upon colliding. We show that correlations between velocities of colliding particles play a crucial role in the long time behavior. In the reaction-controlled limit when particles undergo mostly elastic collisions and therefore are always near equilibrium, the correlations are accounted analytically. For ballistic aggregation, for instance, the density decays as n approximately t(-xi) with xi=2d/(d+3) in the reaction-controlled limit in d dimensions, in contrast with the well-known mean-field prediction xi=2d/(d+2).     

Fluctuations in a one-dimensional mechanical system. I. The Euler limit

We prove the central limit theorem for the density fluctuation field of a one-dimensional mechanical system (hard rods with equal masses and lengths and elastic collisions) in the hydrodynamic limit on the Euler time scale. The limiting process is deterministic and is governed by the linearized Euler equations of the model.     

Ion-induced quark-gluon implosion

We investigate nuclear fragmentation in the central proton-nucleus and nucleus-nucleus collisions at the energies of CERN LHC. Within the semiclassical approximation we argue that because of the fast increase with energy of the cross sections of soft and hard interactions each nucleon is stripped in the average process off "soft" partons and fragments into a collection of leading quarks and gluons with large p(t). Valence quarks and gluons are streaming in the opposite directions when viewed in the c.m. of the produced system. The resulting pattern of the fragmentation of the colliding nuclei leads to an implosion of the quark and gluon constituents of the nuclei. The nonequilibrium state produced at the initial stage in the nucleus fragmentation region is estimated to have densities >/=50 GeV/fm(3) at the LHC energies and probably >/=10 GeV/fm(3) at BNL RHIC.     

Effect of phonon-phonon interactions on localization

We study the heat current J in a classical one-dimensional disordered chain with on-site pinning and with ends connected to stochastic thermal reservoirs at different temperatures. In the absence of anharmonicity all modes are localized and there is a gap in the spectrum. Consequently J decays exponentially with system size N. Using simulations we find that even a small amount of anharmonicity leads to a J approximately 1/N dependence, implying diffusive transport of energy.     

Inelastic scattering and trapping of an atom on a cold,simple-cubic lattice

We report numerical calculations on the inelastic scattering of an incident atom with a simple-cubic harmonic lattice with the collision restricted to the head-on orientation. These calculations represent the first systematic investigation of the behavior of the energy accommodation coefficient and critical trapping energy as the interaction well depth, interaction range, incident energy, and incident mass are varied over significant ranges. These exact results provide a new, substantial “data base” for interpreting experimental results and for testing approximate theories. The incident atom is assumed to collide head-on with only one lattice surface atom. A Morse interaction is used, so that the dynamical effects of both an attractive and a repulsive interaction are included. The lattice is initially cold and classical mechanics is used throughout; the lattice dynamics are treated exactly by means of the response function for the Rosenstock-Newell lattice. We find that long range, soft interactions characteristic of physisorbed species display a broad minimum in energy accommodation at an incident energy near the dissociation energy of the adatom-surface interaction for adatom to lattice atom mass ratios of order one and greater. Stiffer interactions characteristic of chemisorbed species display a pronounced maximum in energy accommodation for mass ratios near 0.7. This is associated with a maximum in the critical trapping energy at the same mass ratio. The results are analyzed to confirm the validity of approximations based on both fast and slow collisions. In the high energy limit, the energy accommodation depends essentially only on the mass ratio, reflecting the pairwise binary collisions. In the slow collision limit, the collision dynamics are described by a newly derived adiabatic approximation for the lattice response.     

Heat flow in an exactly solvable model

A chain of one-dimensional oscillators is considered. They are mechanically uncoupled and interact via a stochastic process which redistributes the energy between nearest neighbors. The total energy is kept constant except for the interactions of the extremal oscillators with reservoirs at different temperatures. The stationary measures are obtained when the chain is finite; the thermodynamic limit is then considered, approach to the Gibbs distribution is proven, and a linear temperature profile is obtained.     

Stationary state of thermostated inelastic hard spheres

We consider a fluid composed of inelastic hard spheres moving in a thermostat modelled by a hard sphere gas. The losses of energy due to inelastic collisions are balanced by the energy transfer via elastic collisions from the thermostat particles. The resulting stationary state is analysed within the Boltzmann kinetic theory. A numerical iterative method permits to study the nature of deviations from the Gaussian state. Some analytic results are obtained for a one-dimensional system.     

High-order harmonic generation spectrum of an excited one-dimensional Coulomb atom in an intense laser field

Response of the wave packet of a one-dimensional Coulomb atom to an intense laser field is calculated using the symmetrized split operator fast Fourier method. The high-order harmonic generation (HHG) of the initial state separately being the ground and excited states is presented. When the hardness parameter \alpha in the soft Coulomb potential V(x)=-1/\sqrt{x^2+\alpha} is chosen to be small enough, the so-called hard Coulomb potential V(x)=-1/|x| can be obtained. It is well known that the hard one-dimensional Coulomb atom has an unstable ground state with an energy eigenvalue of $\sim0.5$ and it has no states corresponding to physical states in the true atoms, and has the first and second excited states being degenerate. The parity effects on the HHG can be seen from the first and second excited states of the hard one-dimensional Coulomb atom. The HHG spectra of the excited states from both the soft and hard Coulomb atom models are shown to have more complex structures and to be much stronger than the corresponding HHG spectrum of the ground state of the soft Coulomb model with $\alpha=2$ in the same laser field. Laser-induced non-resonant one-photon emission is also observed.     

A mechanism for generation of strong energy fluctuations in one-dimensional systems

A mechanism for the generation of strong energy fluctuations in a chain of harmonically and anharmonically coupled atoms in a state of thermodynamic equilibrium has been investigated by a molecular dynamics method. Energy fluctuations multiply exceeding the mean value were localized in a group of three to five atoms with lifetime on the order of an atomic vibrational period. They arise as a consequence of collisions of density waves propagating along the chain. A collision of waves of the same density sign leads to a fluctuation of the potential energy, while for differing signs a fluctuation in the kinetic energy occurs. In particular, collision of dilitatin waves leads to thermal breaking of the chain. The energy concentration at the moment of collision is 30% greater in a harmonic chain than in an anharmonic one.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 5, pp. 47–52, May, 1991.     

Repulsive contact interactions make jammed particulate systems inherently nonharmonic

Many jammed particulate systems, such as granular and colloidal materials, interact via repulsive contact forces. We find that these systems possess no harmonic regime in the large system limit (N→∞) for all compressions Δ? studied, and at jamming onset Δ?→0 for all N. We perform fixed energy simulations following perturbations with amplitude δ along eigendirections of the dynamical matrix. The fluctuations abruptly spread to all modes for δ≈δ(c) (where a single contact breaks) in contrast to linear and weakly nonlinear behavior. For δ > δ(c), all discrete modes disappear into a continuous frequency band. <δ(c)> scales with 1/N and Δ?, which limits harmonic behavior to only overcompressed systems. The density of vibrational modes deviates strongly from that predicted from the dynamical matrix when the system enters the nonharmonic regime, which significantly affects its mechanical and transport properties.