Spin-dynamics simulations of the triangular antiferromagnetic XY model

Using Monte Carlo and spin-dynamics methods, we have investigated the dynamic behavior of the classical, antiferromagnetic XY model on a triangular lattice with linear sizes $L \leq 300$. The temporal evolutions of spin configurations were obtained by solving numerically the coupled equations of motion for each spin using fourth-order Suzuki-Trotter decompositions of exponential operators. From space- and time-displaced spin-spin correlation functions and their space-time Fourier transforms we obtained the dynamic structure factor $S({\bf q},w)$ for momentum ${\bf q}$ and frequency $\omega$. Below $T_{KT}$(Kosterlitz-Thouless transition), both the in-plane ($S^{xx}$) and the out-of-plane ($S^{zz}$) components of $S({\bf q},\omega)$ exhibit very strong and sharp spin-wave peaks. Well above $T_{KT}$, $S^{xx}$ and $S^{zz}$ apparently display a central peak, and spin-wave signatures are still seen in $S^{zz}$. In addition, we also observed an almost dispersionless domain-wall peak at high $\omega$ below $T_{c}$(Ising transition), where long-range order appears in the staggered chirality. Above $T_{c}$, the domain-wall peak disappears for all $q$. The lineshape of these peaks is captured reasonably well by a Lorentzian form. Using a dynamic finite-size scaling theory, we determined the dynamic critical exponent $z$ = 1.002(3). We found that our results demonstrate the consistency of the dynamic finite-size scaling theory for the characteristic frequeny $\omega_{m}$ and the dynamic structure factor $S({\bf q},\omega)$ itself.Comment: 8 pages, RevTex, 10 figures, submitted to PR

Monte Carlo study of the critical temperature for the planar rotator model with nonmagnetic impurities

We performed Monte Carlo simulations to calculate the Berezinskii-Kosterlitz-Thouless (BKT) temperature $T_{BKT}$ for the two-dimensional planar rotator model in the presence of nonmagnetic impurity concentration $(\rho)$. As expected, our calculation shows that the BKT temperature decreases as the spin vacancies increase. There is a critical dilution $\rho_c \approx 0.3$ at which $T_{BKT} =0$. The effective interaction between a vortex-antivortex pair and a static nonmagnetic impurity is studied analytically. A simple phenomenological argument based on the pair-impurity interaction is proposed to justify the simulations.Comment: 5 pages, 5 figures, Revetex fil

Stratified spatiotemporal chaos in anisotropic reaction-diffusion systems

Numerical simulations of two dimensional pattern formation in an anisotropic bistable reaction-diffusion medium reveal a new dynamical state, stratified spatiotemporal chaos, characterized by strong correlations along one of the principal axes. Equations that describe the dependence of front motion on the angle illustrate the mechanism leading to stratified chaos

On the existence of internal modes of sine-Gordon kinks

We study whether or not sine-Gordon kinks exhibit internal modes or ``quasimodes.'' By considering the response of the kinks to ac forces and initial distortions, we show that neither intrinsic internal modes nor ``quasimodes'' exist in contrast to previous reports. However, we do identify a different kind of internal mode bifurcating from the bottom edge of the phonon band which arises from the discretization of the system in the numerical simulations, thus confirming recent predictions.Comment: 4 pages, 2 figures, REVTeX, to appear as a Rapid Communication in Phys Rev E (July 1st

Anomalous resonance phenomena of solitary waves with internal modes

We investigate the non-parametric, pure ac driven dynamics of nonlinear Klein-Gordon solitary waves having an internal mode of frequency $\Omega_i$. We show that the strongest resonance arises when the driving frequency $\delta=\Omega_i/2$, whereas when $\delta=\Omega_i$ the resonance is weaker, disappearing for nonzero damping. At resonance, the dynamics of the kink center of mass becomes chaotic. As we identify the resonance mechanism as an {\em indirect} coupling to the internal mode due to its symmetry, we expect similar results for other systems.Comment: 4 pages, 4 figures, to appear in Phys Rev Let

Internal Modes and Magnon Scattering on Topological Solitons in 2d Easy-Axis Ferromagnets

We study the magnon modes in the presence of a topological soliton in a 2d Heisenberg easy-axis ferromagnet. The problem of magnon scattering on the soliton with arbitrary relation between the soliton radius R and the "magnetic length" Delta_0 is investigated for partial modes with different values of the azimuthal quantum numbers m. Truly local modes are shown to be present for all values of m, when the soliton radius is enough large. The eigenfrequencies of such internal modes are calculated analytically on limiting case of a large soliton radius and numerically for arbitrary soliton radius. It is demonstrated that the model of an isotropic magnet, which admits an exact analytical investigation, is not adequate even for the limit of small radius solitons, R<<Delta_0: there exists a local mode with nonzero frequency. We use the data about local modes to derive the effective equation of soliton motion; this equation has the usual Newtonian form in contrast to the case of the easy-plane ferromagnet. The effective mass of the soliton is found.Comment: 33 pages (REVTeX), 12 figures (EPS

Full Genome Characterization of the Culicoides-Borne Marsupial Orbiviruses: Wallal Virus, Mudjinbarry Virus and Warrego Viruses

Viruses belonging to the species Wallal virus and Warrego virus of the genus Orbivirus were identified as causative agents of blindness in marsupials in Australia during 1994/5. Recent comparisons of nucleotide (nt) and amino acid (aa) sequences have provided a basis for the grouping and classification of orbivirus isolates. However, full-genome sequence data are not available for representatives of all Orbivirus species. We report full-genome sequence data for three additional orbiviruses: Wallal virus (WALV); Mudjinabarry virus (MUDV) and Warrego virus (WARV). Comparisons of conserved polymerase (Pol), sub-core-shell 'T2' and core-surface 'T13' proteins show that these viruses group with other Culicoides borne orbiviruses, clustering with Eubenangee virus (EUBV), another orbivirus infecting marsupials. WARV shares <70% aa identity in all three conserved proteins (Pol, T2 and T13) with other orbiviruses, consistent with its classification within a distinct Orbivirus species. Although WALV and MUDV share <72.86%/67.93% aa/nt identity with other orbiviruses in Pol, T2 and T13, they share >99%/90% aa/nt identities with each other (consistent with membership of the same virus species - Wallal virus). However, WALV and MUDV share <68% aa identity in their larger outer capsid protein VP2(OC1), consistent with membership of different serotypes within the species - WALV-1 and WALV-2 respectively

Magnetic vortex as a ground state for micron-scale antiferromagnetic samples

Here we consider micron-sized samples with any axisymmetric body shape and made with a canted antiferromagnet, like hematite or iron borate. We find that its ground state can be a magnetic vortex with a topologically non-trivial distribution of the sublattice magnetization $\vec{l}$ and planar coreless vortex-like structure for the net magnetization $\vec{M}$. For antiferromagnetic samples in the vortex state, in addition to low-frequency modes, we find high-frequency modes with frequencies over the range of hundreds of gigahertz, including a mode localized in a region of radius $\sim$ 30--40 nm near the vortex core.Comment: 20 pages, 1 figur

Phase transition and landscape statistics of the number partitioning problem

The phase transition in the number partitioning problem (NPP), i.e., the transition from a region in the space of control parameters in which almost all instances have many solutions to a region in which almost all instances have no solution, is investigated by examining the energy landscape of this classic optimization problem. This is achieved by coding the information about the minimum energy paths connecting pairs of minima into a tree structure, termed a barrier tree, the leaves and internal nodes of which represent, respectively, the minima and the lowest energy saddles connecting those minima. Here we apply several measures of shape (balance and symmetry) as well as of branch lengths (barrier heights) to the barrier trees that result from the landscape of the NPP, aiming at identifying traces of the easy/hard transition. We find that it is not possible to tell the easy regime from the hard one by visual inspection of the trees or by measuring the barrier heights. Only the {\it difficulty} measure, given by the maximum value of the ratio between the barrier height and the energy surplus of local minima, succeeded in detecting traces of the phase transition in the tree. In adddition, we show that the barrier trees associated with the NPP are very similar to random trees, contrasting dramatically with trees associated with the $p$ spin-glass and random energy models. We also examine critically a recent conjecture on the equivalence between the NPP and a truncated random energy model

Vortex behavior near a spin vacancy in 2D XY-magnets

The dynamical behavior of anisotropic two dimensional Heisenberg models is still a matter of controversy. The existence of a central peak at all temperatures and a rich structure of magnon peaks are not yet understood. It seems that the central peaks are related, in some way, to structures like vortices. In order to contribute to the discussion of the dynamical behavior of the model we use Monte Carlo and spin dynamics simulations as well analytical calculations to study the behavior of vortices in the presence of nonmagnetic impurities. Our simulations show that vortices are attracted and trapped by the impurities. Using this result we show that if we suppose that vortices are not very much disturbed by the presence of the impurities, then they work as an attractive potential to the vortices explaining the observed behavior in our simulations.Comment: 4 pages, 6 figure