Atomistic simulation of kink structure on edge dislocation in bcc iron

Based on the general theory of dislocation and kink, we have constructed the three kink models corresponding to the 1/2 （111）{011} and 1/2 （111）{112} edge dislocations （EDs） in bcc Fe using the molecular dynamics method. We found that the geometric structure of a kink depends on the type of ED and the structural energies of the atom sites in the dislocation core region, as well as the geometric symmetry of the dislocation core and the characteristic of the stacking sequence of atomic plane along the dislocation line. The formation energies and widths of the kinks on the 1/2 （111）{011} and 1/2 （111）{112} EDs are calculated, the formation energies are 0.05eV and 0.04eV, and widths are 6.02b and 6.51b, respectively （b is the magnitude of the Burgers vector）. The small formation energies indicate that the formation of kink in the edge dislocation is very easy in bcc Fe.     

Oscillating Propagation of Kink in Nondissipative Frenkel-Kontorova Chain Due to External DC Force

We report the oscillating propagation of kink in a nondissipative Frenkel-Kontorova （FK） chain driven by external DC force, which is different from the usual propagation of localized modes with equal speed. When the kink moves in the opposite direction of the external DC force, the kink will be accelerated and the potential of the FK chain in the external force field is transformed to be the kinetic energy of the kink. If the kink reaches the boundary of the FK chain, the kink will be bounced back and moves in the opposite direction, then the kink will be decelerated gradually and the kinetic energy of the kink ＆ transformed to be the potential of the FK chain in the external force field. If the speed of the kink reaches zero, the kink will move in the opposite direction again driven by the external DC force, and a new oscillating cycle begins. Simulation result demonstrates exactly the transformation between the kinetic energy of the kink and the potential of the FK chain in the external force field. The interesting energy exchange is induced by the special topology of kinks, and other localized modes, such as breathers and envelope solitons, have no the interesting phenomenon.     

Dynamic Properties of Two-Dimensional Polydisperse Granular Gases

We propose a two-dimensional model of polydisperse granular mixtures with a power-law size distribution in the presence of stochastic driving. A fractal dimension D is introduced as a measurement of the inhomogeneity of the size distribution of particles. We define the global and partial granular temperatures of the multi-component mixture. By direct simulation Monte Carlo, we investigate how the inhomogeneity of the size distribution influences the dynamic properties of the mixture, focusing on the granular temperature, dissipated energy, velocity distribution, spatial clusterization, and collision time. We get the following results： a single granular temperature does not characterize a multi-component mixture and each species attains its own ＂granular temperature＂; The velocity deviation from Gaussian distribution becomes more and more pronounced and the partial density of the assembly is more inhomogeneous with the increasing value of the fractal dimension D; The global granular temperature decreases and average dissipated energy per particle increases as the value olD augments.     

Pattern Formation in a Vibrated Granular Layer on an Inclined Base

We carry out the simulations of pattern formation in a two-dimensional vibrated granular layer on an inclined base by molecular dynamics. It is found that the maximum amplitude of the pattern is greater at the lower part than at the higher part of the base, and is proportional to the thickness of the layer. Meanwhile, the wavelength varies non-monotoniclly as the inclined angle of the base is increased.     

Propagation of kink–antikink pair along microtubules as a control mechanism for polymerization and depolymerization processes

Among many types of proteinaceous filaments, microtubules(MTs) constitute the most rigid components of the cellular cytoskeleton. Microtubule dynamics is essential for many vital cellular processes such as intracellular transport,metabolism, and cell division. We investigate the nonlinear dynamics of inhomogeneous microtubulin systems and the MT dynamics is found to be governed by a perturbed sine-Gordon equation. In the presence of various competing nonlinear inhomogeneities, it is shown that this nonlinear model can lead to the existence of kink and antikink solitons moving along MTs. We demonstrate kink–antikink pair collision in the framework of Hirota’s bilinearization method. We conjecture that the collisions of the quanta of energy propagating in the form of kinks and antikinks may offer a new view of the mechanism of the retrograde and anterograde transport direction regulation of motor proteins in microtubulin systems.     

MOVING NONLINEAR LOCALIZED MODES FOR ONE-DIMENSIONAL KLEIN－GORDON DIATOMIC LATTICE

We study analytically the moving nonlinear localized vibrational modes (discrete breathers) for a one-dimensional Klein－Gordon diatomic lattice in the whole ω(q) plane of the system by means of a semi- discrete approximation, in which the carrier wave of the modes is treated explicitly while the envelope is described in the continuum approximation. We find that both pulse and kink envelope moving modes for this lattice system can occur with certain carrier wave vectors and vibrational frequencies in separate regions of the ω(q) plane. However, the kink envelope moving modes have not been reported previously for this lattice system.     

Pion transverse-momentum spectrum,elliptic flow,and Hanbury-Brown-Twiss interferometry in a viscous granular source model

We examine the evolution of quark-gluon plasma(QGP) droplets with viscous hydrodynamics and analyze the pion transverse-momentum spectrum, elliptic flow, and Hanbury-Brown-Twiss(HBT) interferometry in a granular source model consisting of viscous QGP droplets. The shear viscosity of the QGP droplet speeds up and slows down the droplet evolution in the central and peripheral regions of the droplet, respectively. The effect of the bulk viscosity on the evolution is negligible. Although there are viscous effects on the droplet evolution, the pion momentum spectrum and elliptic flow change little for granular sources with and without viscosity. On the other hand,the influence of viscosity on HBT radius Rout is considerable. It makes Rout decrease in the granular source model.We determine the model parameters of granular sources using the experimental data of pion transverse-momentum spectrum, elliptic flow, and HBT radii together, and investigate the effects of viscosity on the model parameters. The results indicate that the granular source model may reproduce the experimental data of pion transverse-momentum spectrum, elliptic flow, and HBT radii in heavy-ion collisions of Au-Au at√s_(NN)=200 GeV and Pb-Pb at√s_(NN) =2.76 Te V in different centrality intervals. The viscosity of the droplet leads to an increase in the initial droplet radius and a decrease of the source shell parameter in the granular source model.     

Binary collision approximation for solitary waves in a Y-shaped granular chain

We implement a binary collision approximation to study solitary wave propagation in a two-dimensional double Y-shaped granular chain. The solitary wave was transmitted and reflected when it met the interface of the bifurcated branches of the Y-shaped granular chains. We obtain the analytic results of the ratios of the transmitted and reflected speeds to the incident speed of the solitary wave, the maximum force between the two neighbor beads in a solitary wave, and the total time taken by the pulse to pass through each branch. All of the analytic results are in good agreement with the experimental observations from Daraio et al. [Phys. Rev. E 82 036603 (2010)]. Moreover, we also discuss the delay effects on the arrival of split pulses, and predict the recombination of the split waves traveling in branches in the final stem of asymmetric systems. The prediction of pulse recombination is verified by our numerical results.     

Vortices in dipolar Bose–Einstein condensates in synthetic magnetic field

We study the formation of vortices in a dipolar Bose–Einstein condensate in a synthetic magnetic field by numerically solving the Gross–Pitaevskii equation. The formation process depends on the dipole strength, the rotating frequency, the potential geometry, and the orientation of the dipoles. We make an extensive comparison with vortices created by a rotating trap, especially focusing on the issues of the critical rotating frequency and the vortex number as a function of the rotating frequency. We observe that a higher rotating frequency is needed to generate a large number of vortices and the anisotropic interaction manifests itself as a perceptible difference in the vortex formation. Furthermore, a large dipole strength or aspect ratio also can increase the number of vortices effectively. In particular, we discuss the validity of the Feynman rule.     

Wavelength Selection of Patterns in a Vibrated Granular Layer

We present a numerical study on the pattern formation in a two-dimensional vibrated granular layer by a molecular dynamics algorithm.Through analysing the granular density distribution,we can explore the inner movement process of particles.It is confirmed that there are a dispersive regime and a saturation regime for frequency dependence,between which a critical frequency exists,It is found that there is another saturation regime for thickness dependence,The wavelength increases with increasing layer thickness,but there is a critical thickness after which the wavelength is unchanged.