Time scale for energy equipartition in a two-dimensional FPU model

The FPU problem, i.e., the problem of energy equipartition among normal modes in a weakly nonlinear lattice, is here studied in dimension two, more precisely in a model with triangular cell and nearest-neighbors Lennard-Jones interaction. The number n of degrees of freedom ranges from 182 to 6338. Energy is initially equidistributed among a small number n(0) of low frequency modes, with n(0) proportional to n. We study numerically the time evolution of the so-called spectral entropy and the related "effective number" n(eff) of degrees of freedom involved in the dynamics; in this (rather typical) way we can estimate, for each n and each specific energy (energy per degree of freedom) epsilon, the time scale T(n)(epsilon) for energy equipartition. Numerical results indicate that in the thermodynamic limit the equipartition times are short: more precisely, for large n at fixed epsilon we find a limit curve T(infinity)(epsilon), and T(infinity) grows only as epsilon(-1) for small epsilon. Larger equipartition times are obtained by lowering epsilon, at fixed n, below a crossover value epsilon(c)(n). However, epsilon(c) appears to vanish by increasing n (faster than 1n), and the total energy E=nepsilon, rather than epsilon, appears to be the relevant variable when n is large and epsilon相似文献   

Studies of thermal conductivity in Fermi-Pasta-Ulam-like lattices

The pioneering computer simulations of the energy relaxation mechanisms performed by Fermi, Pasta, and Ulam (FPU) can be considered as the first attempt of understanding energy relaxation and thus heat conduction in lattices of nonlinear oscillators. In this paper we describe the most recent achievements about the divergence of heat conductivity with the system size in one-dimensional (1D) and two-dimensional FPU-like lattices. The anomalous behavior is particularly evident at low energies, where it is enhanced by the quasiharmonic character of the lattice dynamics. Remarkably, anomalies persist also in the strongly chaotic region where long-time tails develop in the current autocorrelation function. A modal analysis of the 1D case is also presented in order to gain further insight about the role played by boundary conditions.     

Interaction of interface plasmons and longitudinal optical phonons in a multilayer structure

The interaction between two-dimensional interface plasmons and longitudinal optical phonons in multilayer structures is analysed in this communication. The dispersion relations for the mixed modes are obtained for periodic boundary conditions in the direction normal to the layer planes. The energy gap between the two bands of mixed modes can be seen in the density of states ?_k(ω) at fixed k parallel to the layers. ?_k(ω) also exhibits singularities at the band edges characteristics of the one-dimensional periodic array of parallel layers.     

Quantized spin waves in antiferromagnetic Heisenberg chains

The quantized stationary spin wave modes in one-dimensional antiferromagnetic spin chains with easy axis on-site anisotropy have been studied by means of Landau-Lifshitz-Gilbert spin dynamics. We demonstrate that the confined antiferromagnetic chains show a unique behavior having no equivalent, neither in ferromagnetism nor in acoustics. The discrete energy dispersion is split into two interpenetrating n and n' levels caused by the existence of two sublattices. The oscillations of individual sublattices as well as the standing wave pattern strongly depend on the boundary conditions. Particularly, acoustical and optical antiferromagnetic spin waves in chains with boundaries fixed (pinned) on different sublattices can be found, while an asymmetry of oscillations appears if the two pinned ends belong to the same sublattice.     

The validity of perturbation theory for the O(N) non-linear sigma models

Recently it has been claimed that ordinary perturbation theory (OPT) gives incorrect weak coupling expansions for lattice O(N) non-linear sigma models in the infinite volume limit, and in particular that the two-dimensional non-abelian models are not asymptotically free, contrary to previous findings. Here it is argued that the problem occurs only for one-dimensional infinite lattices, and that in general, OPT gives correct expansions if physical quantities are first computed on a finite lattice, and the infinite volume limit is taken at the end. In one dimension the expansion is sensitive to boundary conditions because of the severe infrared behavior, but this is not expected to happen in higher dimensions. It is concluded that spin configurations which are far from the perturbative vacuum have too small a measure in the path integral to invalidate OPT, even though they are energetically allowed for non-zero values of the coupling.     

Exponential stability of states close to resonance in infinite-dimensional Hamiltonian systems

We develop canonical perturbation theory for a physically interesting class of infinite-dimensional systems. We prove stability up to exponentially large times for dynamical situations characterized by a finite number of frequencies. An application to two model problems is also made. For an arbitrarily large FPU-like system with alternate light and heavy masses we prove that the exchange of energy between the optical and the acoustical modes is frozen up to exponentially large times, provided the total energy is small enough. For an infinite chain of weakly coupled rotators we prove exponential stability for two kinds of initial data: (a) states with a finite number of excited rotators, and (b) states with the left part of the chain uniformly excited and the right part at rest.     

Quantum Communication Through a Two-Dimensional Spin Network

We investigate the state or entanglement transfer through a two-dimensional spin network. We show that for state transfer, better fidelity can be gained along the diagonal direction but for entanglement transfer, when the initial entanglement is created along the boundary, the concurrence is more inclined to propagate along the boundary. This behavior is produced by quantum mechanical interference and the communication quality depends on the precise size of the network. For some number of sites, the fidelity in a two-dimensional channel is higher than one-dimensional case. This is an important result for realizing quantum communication through high dimension spin chain networks.     

Quantum Communication Through a Two-Dimensional Spin Network

We investigate the state or entanglement transfer through a two-dimensional spin network. We show that for state transfer, better fidelity can be gained along the diagonal direction but for entanglement transfer, when the initial entanglement is created along the boundary, the concurrence is more inclined to propagate along the boundary. This behavior is produced by quantum mechanical interference and the communication quality depends on the precise size of the network. For some number of sites, the fidelity in a two-dimensional channel is higher than one-dimensional case. This is an important result for realizing quantum communication through high dimension spin chain networks.     

One-dimensional "turbulence" in a discrete lattice

We study a one-dimensional discrete analog of the von Karman flow, widely investigated in turbulence. A lattice of anharmonic oscillators is excited by both ends in order to create a large scale structure in a highly nonlinear medium, in the presence of a dissipative term proportional to the second order finite difference of the velocities, similar to the viscous term in a fluid. In a first part, the energy density is investigated in real and Fourier space in order to characterize the behavior of the system on a local scale. At low amplitude of excitation the large scale structure persists in the system but all modes are however excited and exchange energy, leading to a power law spectrum for the energy density, which is remarkably stable against changes in the model parameters, amplitude of excitation, or damping. In the spirit of shell models, this regime can be described in terms of interacting scales. At higher amplitude of excitation, the large scale structure is destroyed and the dynamics of the system can be viewed as resulting from the creation, interaction, and decay of localized excitations, the discrete breathers, the one-dimensional equivalents of vortices in a fluid. The spectrum of the energy density is well described by the spectrum of the breathers, and shows an exponential decay with the wave vector. Due to this exponential behavior, the spectrum is dominated by the most intense breathers. In this regime, the probability distribution of the increments of velocity between neighboring points is remarkably similar to the experimental results of turbulence and can be described by distributions deduced from nonextensive thermodynamics as in fluids. In a second part the power dissipated in the whole lattice is studied to characterize the global behavior of the system. Its probability distribution function shows non-Gaussian fluctuations similar to the one exhibited recently in a large class of "inertial systems," i.e., systems that cannot be divided into mesoscopic regions which are independent. The properties of the nonlinear excitations of the lattice provide a partial understanding of this behavior.     

A two-dimensional vibration analysis of piezoelectrically actuated microbeam with nonideal boundary conditions

In this paper, the influences of nonideal boundary conditions (due to flexibility) on the primary resonant behavior of a piezoelectrically actuated microbeam have been studied, for the first time. The structure has been assumed to treat as an Euler–Bernoulli beam, considering the effects of geometric nonlinearity. In this work, the general nonideal supports have been modeled as a the combination of horizontal, vertical and rotational springs, simultaneously. Allocating particular values to the stiffness of these springs provides the mathematical models for the majority of boundary conditions. This consideration leads to use a two-dimensional analysis of the multiple scales method instead of previous works' method (one-dimensional analysis). If one neglects the nonideal effects, then this paper would be an effort to solve the two-dimensional equations of motion without a need of a combination of these equations using the shortening or stretching effect. Letting the nonideal effects equal to zero and comparing their results with the results of previous approaches have been demonstrated the accuracy of the two-dimensional solutions. The results have been identified the unique effects of constraining and stiffening of boundaries in horizontal, vertical and rotational directions. This means that it is inaccurate to suppose the nonideality of supports only in one or two of these directions like as previous works. The findings are of vital importance as a better prediction of the frequency response for the nonideal supports. Furthermore, the main findings of this effort can help to choose appropriate boundary conditions for desired systems.     

Cecotti-Fendley-Intriligator-Vafa index in a box

The Cecotti-Fendley-Intriligator-Vafa (CFIV) index in two-dimensional N=(2,2) models is revisited. We address the problem of “elementary” (nontopological) excitations over a kink solution, in the one-kink sector (using the Wess-Zumino or Landau-Ginzburg models with two vacua as examples). In other words, we limit ourselves to the large-β limit. The excitation spectrum over the BPS-saturated at the classical level kink is discretized through a large box with judiciously chosen boundary conditions. The boundary conditions are designed in such a way that half of supersymmetry is preserved as well as the BPS kink itself, and relevant zero modes. The excitation spectrum acquires a mass gap. All (discretized) excited states enter in four-dimensional multiplets (two bosonic states + two fermionic). Their contribution to indCFIV vanishes level-by-level. The ground state contribution produces |indCFIV|=1.     

A one-to-two-dimensional mapping using a partial Fast Fourier Transform

It will be shown how to map a simple one-dimensional tight binding model with a cosine potential in one dimension exactly to a two-dimensional tight binding model with periodic boundary conditions in the presence of a single flux quantum spread evenly on the torus. The mapping is achieved by a partial sequence of “Fast Fourier Transform” (FFT) steps which if completed would be an exact Fourier transform of the original model. Each step of the FFT recursively maps a tight binding model into two decoupled sublattices of half the lattice length.     

Effective Dynamics for a Kinetic Monte–Carlo Model with Slow and Fast Time Scales

We consider several multiscale-in-time kinetic Monte Carlo models, in which some variables evolve on a fast time scale, while the others evolve on a slow time scale. In the first two models we consider, a particle evolves in a one-dimensional potential energy landscape which has some small and some large barriers, the latter dividing the state space into metastable regions. In the limit of infinitely large barriers, we identify the effective dynamics between these macro-states, and prove the convergence of the process towards a kinetic Monte Carlo model. We next consider a third model, which consists of a system of two particles. The state of each particle evolves on a fast time-scale while conserving their respective energy. In addition, the particles can exchange energy on a slow time scale. Considering the energy of the first particle, we identify its effective dynamics in the limit of asymptotically small ratio between the characteristic times of the fast and the slow dynamics. For all models, our results are illustrated by representative numerical simulations.     

Fluctuations and large deviations in non-equilibrium systems

For systems in contact with two reservoirs at different densities or with two thermostats at different temperatures, the large deviation function of the density gives a possible way of extending the notion of free energy to non-equilibrium systems. This large deviation function of the density can be calculated explicitly for exclusion models in one dimension with open boundary conditions. For these models, one can also obtain the distribution of the current of particles flowing through the system and the results lead to a simple conjecture for the large deviation function of the current of more general diffusive systems.     

磷化物三元混晶中电子-声子相互作用的维度和体积效应

考虑电子-声子耦合强度因维度而异,导出了描述三维、二维和一维混晶中电子-声子相互作用的哈密顿量。考虑构成三元混晶的两种二元晶体的晶格失配会使混晶体积随元素组分比改变,在推导三维、二维和一维三元混晶中极化子自陷能量和重整化有效质量时计入了离子相对位移与二元晶体原胞体积的关系。结果表明:磷化物三元混晶中极化子自陷能量和重整化有效质量随元素组分的变化关系呈明显的非线性特征,对晶格适配明显、电子-声子耦合较强的材料,体积效应不可忽略。维度越低,非线性特征和体积效应越明显。     

Ag/Ni和Cu/Ni界面的分子动力学模拟

用分子动力学方法对金属界面在弯曲状态下的力学行为做了模拟计算．在自行设计的两种弯曲模型中，首先比较了Ａｇ／Ｎｉ在不形成界面、形成界面（错配比约为１５％）以及Ｃｕ／Ｎｉ形成界面（错配比约为３％）时在动态弯曲过程中的势能－应变曲线，应力－应变曲线，模量－应变曲线，通过比较得出的结论是：界面的存在影响很大，失配位错影响界面的性质，并且错配比不同界面的力学性质亦不相同．同时，对计算元胞的尺寸效应做了详细的讨论，给出了用于计算机模拟中比较适宜的计算元胞的尺寸．最后，利用圆弧弯曲模型将静态平移周期性边界条件应用于动态 关键词：     

Finite-size scaling studies of one-dimensional reaction-diffusion systems. Part I. Analytical results

We consider two single-species reaction-diffusion models on one-dimensional lattices of lengthL: the coagulation-decoagulation model and the annihilation model. For the coagulation model the system of differential equations describing the time evolution of the empty interval probabilities is derived for periodic as well as for open boundary conditions. This system of differential equations grows quadratically withL in the latter case. The equations are solved analytically and exact expressions for the concentration are derived. We investigate the finite-size behavior of the concentration and calculate the corresponding scaling functions and the leading corrections for both types of boundary conditions. We show that the scaling functions are independent of the initial conditions but do depend on the boundary conditions. A similarity transformation between the two models is derived and used to connect the corresponding scaling functions.     

Intersubband relaxation of hot carriers in quantum well systems

We use an ensemble Monte Carlo simulation of coupled electrons, holes and nonequilibrium polar optical phonons in multiple quantum well systems to model the intersubband relaxation of hot carriers measured in ultra-fast optical experiments. We have investigated the effect of various models of confined photon modes on the energy relaxation and intersubband transition rate in single quantum well and coupled well systems. In particular, the symmetry of the atomic displacement with respect to the quantum well has a marked effect on the relative intersubband versus intrasubband scattering rates, depending on whether one considers electrostatic boundary conditions(slab modes) or mechanical boundary conditions(guided modes). In single quantum wells systems, the overall intersubband relaxation time is not found to be strongly dependent on the confined mode model used due to competing effects of hot phonons and the relative intrasubband scattering rates. For coupled well systems, the relaxation rate is much more dependent on the exact nature of the phonon amplitude. Large effects are found associated with localized AlAs interface modes which dominate the intersubband relaxation time.     

Interaction properties of the periodic and step-like solutions of the double-Sine-Gordon equation

The periodic and step-like solutions of the double-Sine-Gordon equation are investigated, with different initial conditions and for various values of the potential parameter epsilon. We plot energy and force diagrams, as functions of the inter-soliton distance for such solutions. This allows us to consider our system as an interacting many-body system in 1+1 dimension. We therefore plot state diagrams (pressure vs. average density) for step-like as well as periodic solutions. Step-like solutions are shown to behave similarly to their counterparts in the Sine-Gordon system. However, periodic solutions show a fundamentally different behavior as the parameter epsilon is increased. We show that two distinct phases of periodic solutions exist which exhibit manifestly different behavior. Response functions for these phases are shown to behave differently, joining at an apparent phase transition point.     

Polymer effects on small- and large-scale two-dimensional turbulence

We investigate the effect of dilute polymers on driven two-dimensional turbulence in a soap film. Transitions from strong to weak turbulence are identified by independently varying the polymer concentration phi and the energy injection rate epsilon(inj) . Studies of velocity structures in small scales reveal that strong saddles are suppressed, whereas weak ones become more populated. Interestingly, this redistribution of saddle points in turbulent flows strongly correlates with the quenching of velocity fluctuations on large scales, suggesting that this hydrodynamic structure may play a role in transferring energy from scale to scale.