Symmetry analysis and Monte Carlo study of a frustrated antiferromagnetic planar (XY) model in two dimensions

L-based numerical linked cluster expansion for square lattice models

We introduce a numerical linked cluster expansion for square-lattice models whose building block is an L-shape cluster. For the spin-1/2 models studied in this work, we find that this expansion exhibits a similar or better convergence of the bare sums than that of the (larger) square-shaped clusters and can be used with resummation techniques (like the site- and bond-based expansions) to obtain results at even lower temperatures. We compare the performance of weak- and strong-embedding versions of this expansion in various spin-1/2 models and show that the strong-embedding version is preferable because of its convergence properties and lower computational cost. Finally, we show that the expansion based on the L-shape cluster can be naturally used to study properties of lattice models that smoothly connect the square and triangular lattice geometries.

Superlattices of Gadolinium and Bismuth Based Thallium Dichalcogenides as Potential Magnetic Topological Insulators

Using relativistic spin-polarized density functional theory calculations we investigate magnetism, electronic structure and topology of the ternary thallium gadolinium dichalcogenides TlGdZ2 (Z= Se and Te) as well as superlattices on their basis. We find TlGdZ2 to have an antiferromagnetic exchange coupling both within and between the Gd layers, which leads to frustration and a complex magnetic structure. The electronic structure calculations reveal both TlGdSe2 and TlGdTe2 to be topologically trivial semiconductors. However, as we show further, a three-dimensional (3D) magnetic topological insulator (TI) state can potentially be achieved by constructing superlattices of the TlGdZ2/(TlBiZ2)n type, in which structural units of TlGdZ2 are alternated with those of the isomorphic TlBiZ2 compounds, known to be non-magnetic 3D TIs. Our results suggest a new approach for achieving 3D magnetic TI phases in such superlattices which is applicable to a large family of thallium rare-earth dichalcogenides and is expected to yield a fertile and tunable playground for exotic topological physics.

Phase diagrams of the antiferromagnetic XY model on a triangular lattice with higher-order interactions

We study the effects of higher-order antinematic interactions on the critical behavior of the antiferromagnetic (AFM) XY model on a triangular lattice using Monte Carlo simulations. The parameter q of the generalized antinematic (ANq) interaction is found to have a pronounced effect on the phase diagram topology by inducing new quasi-long-range ordered phases due to competition with the conventional AFM interaction as well as geometrical frustration. For values of q divisible by 3, the ground-state competition between the two interactions results in a frustrated canted AFM phase appearing at low temperatures wedged between the AFM and ANq phases. For q nondivisible by 3, with the increase of q one can observe the evolution of the phase diagram topology featuring two (q=2), three (q=4,5), and four (q≥7) ordered phases. In addition to the two phases previously found for q=2, the first new phase with solely AFM ordering arises for q=4 in the limit of strong AFM coupling and higher temperatures by separating from the phase with the coexisting AFM and ANq orderings. For q=7, another phase with AFM ordering but multimodal spin distribution in each sublattice appears at intermediate temperatures. All these algebraic phases also display standard and generalized chiral long-range orderings, which decouple at higher temperatures in the regime of dominant ANq (AFM) interaction for q≥4 (q≥7) preserving only the generalized (standard) chiral ordering.

New ordered phase in geometrically frustrated generalized XY model

Critical properties of a geometrically frustrated generalized XY model with antiferromagnetic (AFM) and third-order antinematic (AN3) couplings on a triangular lattice are studied by Monte Carlo simulation. It is found that such a generalization leads to a phase diagram consisting of three different quasi-long-range ordered (QLRO) phases. Compared to the model with the second-order antinematic (AN2) coupling, besides the AFM and AN3 phases which appear in the limits of relatively strong AFM and AN3 interactions, respectively, it includes an additional complex canted antiferromagnetic (CAFM) phase. It emerges at lower temperatures, wedged between the AFM and AN3 phases as a result of the competition between the AFM and the AN3 couplings, which is absent in the model with the AN2 coupling. The AFM-CAFM and AN3-CAFM phase transitions are concluded to belong to the weak Ising and weak three-state Potts universality classes, respectively. Additionally, all three QLRO phases also feature true LRO of the standard and generalized chiralities, which both vanish simultaneously at second-order phase transitions with non-Ising critical exponents and the critical temperatures slightly higher than the magnetic and nematic order-disorder transition temperatures.

Effects of interlayer and bi-quadratic exchange coupling on layered triangular lattice antiferromagnets

The magnetic field evolution of ground spin states of the stacked planar triangular antiferromagnet with antiferromagnetic interlayer interaction J _c is explored using a minimal 3D classical Heisenberg model. A bi-quadratic coupling is also used to mimic the effect of spin fluctuations (Zhitomirsky 2015 J. Phys.: Conf. Ser. 592 012110) which are known to stabilize the magnetization plateau. A single ion anisotropy is included and states with a magnetic field applied in the ab plane and along the c axis are determined. For [Formula: see text]-plane, an additional new state, in contrast to 2D model (Zhitomirsky 2015 J. Phys.: Conf. Ser. 592 012110), is obtained with weak interlayer interaction, while the magnetization plateau vanishes at large J _c and other new states with z components of spins emerge. For [Formula: see text]-axis, an extra state, compared with 2D model, is obtained with a weak interlayer interaction. When J _c is large enough, only the state corresponding to the Umbrella phase in 2D model exits.

Phase-field crystal for an antiferromagnet with elastic interactions

We introduce a model which contains the essential elements to formulate and study antiferromagnetism, using the phase-field crystal framework. We focus on the question of how magneto-elastic coupling could lift the frustration in the two-dimensional hexagonal antiferromagnetic phase. Using simulations we observe a rich variety of different phases stable in this model. To characterize different phases we calculate the chiral order parameter and identify the scaling behavior of this order parameter. Furthermore, we observe that vortices appear and are stable close to the nonmagnetic defects. Finally, we studied the ferrimagnetic and spin-flop phase transition in the presence of an external magnetic field.

Manifestation of long-range ordered state in layered VX2 (X = Cl, Br, I) systems

We report the first-principles based density functional investigation of correlated magnetic behavior for the layered VX2 (X = Cl, Br, I) system both in bulk and monolayer forms. The bulk system stabilizes in the Néel 120° spin structure with the onset of long-range order at a very low temperature. However, a monolayer of the same system results in a two-dimensional long-range ordered spin-gel configuration in the a-b plane. The occurrence of this two-dimensional non-collinear long-range order for such hexagonal systems with honeycomb topology may provide an additional feature for controlling the use of such materials in magnetic memory devices.

Magnetic phase diagrams of classical triangular and kagome antiferromagnets

We investigate the effect of geometrical frustration on the H-T phase diagrams of the classical Heisenberg antiferromagnets on triangular and kagome lattices. The phase diagrams for the two models are obtained from large-scale Monte Carlo simulations. For the kagome antiferromagnet, thermal fluctuations are unable to lift degeneracy completely and stabilize translationally disordered multipolar phases. We find a substantial difference in the temperature scales of the order by disorder effect related to different degeneracy of the low- and the high-field classical ground states in the kagome antiferromagnet. In the low-field regime, the Kosterlitz-Thouless transition into a spin-nematic phase is produced by unbinding of half-quantum vortices.

Phase and frequency entrainment in locally coupled phase oscillators with repulsive interactions

Recent experiments in one- and two-dimensional microfluidic arrays of droplets containing Belousov-Zhabotinsky reactants show a rich variety of spatial patterns [M. Toiya et al., J. Phys. Chem. Lett. 1, 1241 (2010)]. The dominant coupling between these droplets is inhibitory. Motivated by this experimental system, we study repulsively coupled Kuramoto oscillators with nearest-neighbor interactions, on a linear chain as well as a ring in one dimension, and on a triangular lattice in two dimensions. In one dimension, we show using linear stability analysis as well as numerical study that the stable phase patterns depend on the geometry of the lattice. We show that a transition to the ordered state does not exist in the thermodynamic limit. In two dimensions, we show that the geometry of the lattice constrains the phase difference between two neighboring oscillators to 2π/3. We report the existence of domains with either clockwise or anticlockwise helicity, leading to defects in the lattice. We study the time dependence of these domains and show that at large coupling strengths the domains freeze due to frequency synchronization. Signatures of the above phenomena can be seen in the spatial correlation functions.

String picture for a model of frustrated quantum magnets and dimers

We study the effect of quantum dynamics on geometrically frustrated magnets for a transverse field Ising model at finite temperatures. We develop a microscopic derivation of the Landau-Ginzburg-Wilson action for this model and show that it can be interpreted as the free energy of a 3D elastic lattice of noncrossing strings. As a first application, we quantitatively predict the phase diagram and correlations, confirming excellently a key prediction of recent simulations about the existence of unusual phase transitions and an ordered phase. We discuss the implications of our string picture for the understanding of the effect of quenched disorder in such quantum frustrated systems.

Statistical mechanics of columnar DNA assemblies

Many physical systems can be mapped onto solved or "solvable" models of magnetism. In this work, we have mapped the statistical mechanics of columnar phases of ideally helical rigid DNA--subject to the earlier found unusual, frustrated pair potential (A.A. Kornyshev, S. Leikin, J. Chem. Phys. 107, 3656 (1997))--onto an exotic, unknown variant of the XY model on a fixed or restructurable lattice. Here, the role of the "spin" is played by the azimuthal orientation of the molecules. We have solved this model using a Hartree-Fock approximation, ground-state calculations, and finite-temperature Monte Carlo simulations. We have found peculiar spin order transitions, which may also be accompanied by positional restructuring, from hexagonal to rhombohedric lattices. Some of these have been experimentally observed in dense columnar aggregates. Note that DNA columnar phases are of great interest in biophysical research, not only because they are a useful in vitro tool for the study of DNA condensation, but also since these structures have been detected in living matter. Within the approximations made, our study provides insight into the statistical mechanics of these systems.

Critical behavior of the frustrated antiferromagnetic six-state clock model on a triangular lattice

We study the antiferromagnetic six-state clock model with nearest neighbor interactions on a triangular lattice with extensive Monte Carlo simulations. We find clear indications of two phase transitions at two different temperatures: Below T(I) a chirality order sets in and by a thorough finite-size-scaling analysis of the specific heat and the chirality correlation length we show that this transition is in the Ising universality class (with a nonvanishing chirality order parameter below T(I)). At T(KT) (<T(I)) the spin-spin correlation length as well as the spin susceptibility diverges according to a Kosterlitz-Thouless (KT) form and spin correlations decay algebraically below T(KT). We compare our results to recent x-ray diffraction experiments on the orientational ordering of CF3Br monolayers physisorbed on graphite. We argue that the six-state clock model describes the universal feature of the phase transition in the experimental system and that the orientational ordering belongs to the KT universality class.