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

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

Critical behavior of the planar magnet model in three dimensions

We use a hybrid Monte Carlo algorithm in which a single-cluster update is combined with the over-relaxation and Metropolis spin re-orientation algorithm. Periodic boundary conditions were applied in all directions. We have calculated the fourth-order cumulant in finite size lattices using the single-histogram re-weighting method. Using finite-size scaling theory, we obtained the critical temperature which is very different from that of the usual XY model. At the critical temperature, we calculated the susceptibility and the magnetization on lattices of size up to $42^3$. Using finite-size scaling theory we accurately determine the critical exponents of the model and find that $\nu$=0.670(7), $\gamma/\nu$=1.9696(37), and $\beta/\nu$=0.515(2). Thus, we conclude that the model belongs to the same universality class with the XY model, as expected.Comment: 11 pages, 5 figure

Low Temperature Static and Dynamic Behavior of the Two-Dimensional Easy-Axis Heisenberg Model

We apply the self-consistent harmonic approximation (SCHA) to study static and dynamic properties of the two-dimensional classical Heisenberg model with easy-axis anisotropy. The static properties obtained are magnetization and spin wave energy as functions of temperature, and the critical temperature as a function of the easy-axis anisotropy. We also calculate the dynamic correlation functions using the SCHA renormalized spin wave energy. Our analytical results, for both static properties and dynamic correlation functions, are compared to numerical simulation data combining cluster-Monte Carlo algorithms and Spin Dynamics. The comparison allows us to conclude that far below the transition temperature, where the SCHA is valid, spin waves are responsible for all relevant features observed in the numerical simulation data; topological excitations do not seem to contribute appreciably. For temperatures closer to the transition temperature, there are differences between the dynamic correlation functions from SCHA theory and Spin Dynamics; these may be due to the presence of domain walls and solitons.Comment: 12 pages, 14 figure