Ising transition driven by frustration in a 2D classical model with continuous symmetry

Beyond Kitaev physics in strong spin-orbit coupled magnets

We review the recent advances and current challenges in the field of strong spin-orbit coupled Kitaev materials, with a particular emphasis on the physics beyond the exactly-solvable Kitaev spin liquid point. To this end, we present a comprehensive overview of the key exchange interactions in candidate materials with a specific focus on systems featuring effectiveJeff=1/2magnetic moments. This includes, but not limited to,5d5iridates,4d5ruthenates and3d7cobaltates. Our exploration covers the microscopic origins of these interactions, along with a systematic attempt to map out the most intriguing correlated regimes of the multi-dimensional parameter space. Our approach is guided by robust symmetry and duality transformations as well as insights from a wide spectrum of analytical and numerical studies. We also survey higher spin Kitaev models and recent exciting results on quasi-one-dimensional models and discuss their relevance to higher-dimensional models. Finally, we highlight some of the key questions in the field as well as future directions.

Thermal Ising Transition in the Spin-1/2 J_{1}-J_{2} Heisenberg Model

Using an SU(2) invariant finite-temperature tensor network algorithm, we provide strong numerical evidence in favor of an Ising transition in the collinear phase of the spin-1/2 J_{1}-J_{2} Heisenberg model on the square lattice. In units of J_{2}, the critical temperature reaches a maximal value of T_{c}/J_{2}≃0.18 around J_{2}/J_{1}≃1.0. It is strongly suppressed upon approaching the zero-temperature boundary of the collinear phase J_{2}/J_{1}≃0.6, and it vanishes as 1/log(J_{2}/J_{1}) in the large J_{2}/J_{1} limit, as predicted by Chandra et al., [Phys. Rev. Lett. 64, 88 (1990)PRLTAO0031-900710.1103/PhysRevLett.64.88]. Enforcing the SU(2) symmetry is crucial to avoid the artifact of finite-temperature SU(2) symmetry breaking of U(1) algorithms, opening new perspectives in the investigation of the thermal properties of quantum Heisenberg antiferromagnets.

Locating fixed points in the phase plane

The critical point is a fixed point in finite-size scaling. To quantify the behavior of such a fixed point, we define, at a given temperature and scaling exponent ratio, the width of scaled observables for different sizes. The minimum of the width reveals the position of the fixed point, its corresponding phase transition temperature, and scaling exponent ratio. The value of this ratio tells the nature of the fixed point, which can be a critical point, a point of the first-order phase transition line, or a point of the crossover region. To demonstrate the effectiveness of this method, we apply it to three typical samples produced by the three-dimensional three-state Potts model. Results show the method to be more precise and effective than conventional methods. Finally, we discuss a possible application at the Beam Energy Scan program of the Relativistic Heavy Ion Collider.

Nematic Bond Theory of Heisenberg Helimagnets

We study classical two-dimensional frustrated Heisenberg models with generically incommensurate ground states. A new theory for the lattice-nematic "order by disorder" transition is developed based on the self-consistent determination of the effective exchange coupling bonds. In our approach, fluctuations of the constraint field imposing conservation of the local magnetic moment drive nematicity at low temperatures. The critical temperature is found to be highly sensitive to the peak helimagnetic wave vector, and vanishes continuously when approaching rotation symmetric Lifshitz points. Transitions between symmetry distinct nematic orders may occur by tuning the exchange parameters, leading to lines of bicritical points.

Non-Zhang-Rice Singlet Character of the First Ionization State of T-CuO

We argue that tetragonal CuO (T-CuO) has the potential to finally settle long-standing modeling issues for cuprate physics. We compare the one-hole quasiparticle (qp) dispersion of T-CuO to that of cuprates, in the framework of the strongly correlated (U_{dd}→∞) limit of the three-band Emery model. Unlike in CuO_{2}, magnetic frustration in T-CuO breaks the C_{4} rotational symmetry and leads to strong deviations from the Zhang-Rice singlet picture in parts of the reciprocal space. Our results are consistent with angle-resolved photoemission spectroscopy data but in sharp contradiction to those of a one-band model previously suggested for them. These differences identify T-CuO as an ideal material to test a variety of scenarios proposed for explaining cuprate phenomenology.

Strongly frustrated triangular spin lattice emerging from triplet dimer formation in honeycomb Li2IrO3

Iridium oxides with a honeycomb lattice have been identified as platforms for the much anticipated Kitaev topological spin liquid: the spin-orbit entangled states of Ir(4+) in principle generate precisely the required type of anisotropic exchange. However, other magnetic couplings can drive the system away from the spin-liquid phase. With this in mind, here we disentangle the different magnetic interactions in Li2IrO3, a honeycomb iridate with two crystallographically inequivalent sets of adjacent Ir sites. Our ab initio many-body calculations show that, while both Heisenberg and Kitaev nearest-neighbour couplings are present, on one set of Ir-Ir bonds the former dominates, resulting in the formation of spin-triplet dimers. The triplet dimers frame a strongly frustrated triangular lattice and by exact cluster diagonalization we show that they remain protected in a wide region of the phase diagram.

Emergent Power-Law Phase in the 2D Heisenberg Windmill Antiferromagnet: A Computational Experiment

In an extensive computational experiment, we test Polyakov's conjecture that under certain circumstances an isotropic Heisenberg model can develop algebraic spin correlations. We demonstrate the emergence of a multispin U(1) order parameter in a Heisenberg antiferromagnet on interpenetrating honeycomb and triangular lattices. The correlations of this relative phase angle are observed to decay algebraically at intermediate temperatures in an extended critical phase. Using finite-size scaling we show that both phase transitions are of the Berezinskii-Kosterlitz-Thouless type, and at lower temperatures we find long-range Z(6) order.

Unusual ordered phases of highly frustrated magnets: a review

We review ground states and excitations of a quantum antiferromagnet on triangular and other frustrated lattices. We pay special attention to the combined effects of magnetic field h, spatial anisotropy R and spin magnitude S. The focus of the review is on the novel collinear spin density wave and spin nematic states, which are characterized by fully gapped transverse spin excitations with S(z) = ± 1. We discuss extensively the R - h phase diagram of the antiferromagnet, both in the large-S semiclassical limit and the quantum S = 1/2 limit. When possible, we point out connections with experimental findings.

Angle-resolved photoemission spectroscopy of tetragonal CuO: evidence for intralayer coupling between cupratelike sublattices

We investigate by angle-resolved photoemission the electronic structure of in situ grown tetragonal CuO, a synthetic quasi-two-dimensional edge-sharing cuprate. We show that, in spite of the very different nature of the copper oxide layers, with twice as many Cu in the CuO layers of tetragonal CuO as compared to the CuO(2) layers of the high-T(c) cuprates, the low-energy electronic excitations are surprisingly similar, with a Zhang-Rice singlet dispersing on weakly coupled cupratelike sublattices. This system should thus be considered as a member of the high-T(c) cuprate family, with, however, interesting differences due to the intralayer coupling between the cupratelike sublattices.

Cluster mean-field theory study of J1-J2 Heisenberg model on a square lattice

We study the spin-1/2 J1-J2 Heisenberg model on a square lattice using the cluster mean-field theory. We find a rapid convergence of phase boundaries with increasing cluster size. By extrapolating the cluster size L to infinity, we obtain accurate phase boundaries J(c1)(2) ≈ 0.42 (between the Néel antiferromagnetic phase and non-magnetic phase), and J(c2)(2) ≈ 0.59 (between non-magnetic phase and the collinear antiferromagnetic phase). Our results support the second-order phase transition at J(c1)(2) and the first-order one at J(c2)(2). For the spin-anisotropic J1-J2 model, we present its finite temperature phase diagram and demonstrate that the non-magnetic state is unstable towards the first-order phase transition under intermediate spin anisotropy.

Interlayer-interaction dependence of latent heat in the Heisenberg model on a stacked triangular lattice with competing interactions

We study the phase transition behavior of a frustrated Heisenberg model on a stacked triangular lattice by Monte Carlo simulations. The model has three types of interactions: the ferromagnetic nearest-neighbor interaction J(1) and antiferromagnetic third nearest-neighbor interaction J(3) in each triangular layer and the ferromagnetic interlayer interaction J([perpendicular]). Frustration comes from the intralayer interactions J(1) and J(3). We focus on the case that the order parameter space is SO(3)×C(3). We find that the model exhibits a first-order phase transition with breaking of the SO(3) and C(3) symmetries at finite temperature. We also discover that the transition temperature increases but the latent heat decreases as J([perpendicular])/J(1) increases, which is opposite to the behavior observed in typical unfrustrated three-dimensional systems.

Broken discrete and continuous symmetries in two-dimensional spiral antiferromagnets

We study the occurrence of symmetry breaking, at zero and finite temperatures, in the J1-J3 antiferromagnetic Heisenberg model on the square lattice using Schwinger boson mean field theory. For spin-1/2 the ground state always breaks the SU(2) symmetry with a continuous quasi-critical transition at J3/J1 ∼ 0.38, from Néel to spiral long range order, although local spin fluctuation considerations suggest an intermediate disordered regime around 0.35 ≲ J3/J1 ≲ 0.5, in qualitative agreement with recent numerical results. At low temperatures we find a Z2 broken symmetry region with short range spiral order characterized by an Ising-like nematic order parameter that compares qualitatively well with classical Monte Carlo results. At intermediate temperatures the phase diagram shows regions with collinear short range orders: for J3/J1 < 1 Néel (π,π) correlations and for J3/J1 > 1 a novel phase consisting of four decoupled third neighbour sublattices with Néel (π,π) correlations in each one. We conclude that the effect of quantum and thermal fluctuations is to favour collinear correlations even in the strongly frustrated regime.

Emergent critical phase and Ricci flow in a 2D frustrated Heisenberg model

We introduce a two-dimensional frustrated Heisenberg antiferromagnet on interpenetrating honeycomb and triangular lattices. Classically the two sublattices decouple, and "order from disorder" drives them into a coplanar state. Applying Friedan's geometric approach to nonlinear sigma models, we obtain the scaling of the spin stiffnesses governed by the Ricci flow of a four-dimensional metric tensor. At low temperatures, the relative phase between the spins on the two sublattices is described by a six-state clock model with an emergent critical phase.

Bond orientational order in liquids: Towards a unified description of water-like anomalies, liquid-liquid transition, glass transition, and crystallization: Bond orientational order in liquids

There are at least three fundamental states of matter, depending upon temperature and pressure: gas, liquid, and solid (crystal). These states are separated by first-order phase transitions between them. In both gas and liquid phases a complete translational and rotational symmetry exist, whereas in a solid phase both symmetries are broken. In intermediate phases between liquid and solid, which include liquid crystal and plastic crystal phases, only one of the two symmetries is preserved. Among the fundamental states of matter, the liquid state is the most poorly understood. We argue that it is crucial for a better understanding of liquids to recognize that a liquid generally has the tendency to have a local structural order and its presence is intrinsic and universal to any liquid. Such structural ordering is a consequence of many-body correlations, more specifically, bond angle correlations, which we believe are crucial for the description of the liquid state. We show that this physical picture may naturally explain difficult unsolved problems associated with the liquid state, such as anomalies of water-type liquids (water, Si, Ge, ...), liquid-liquid transition, liquid-glass transition, crystallization and quasicrystal formation, in a unified manner. In other words, we need a new order parameter representing a low local free-energy configuration, which is a bond orientational order parameter in many cases, in addition to a density order parameter for the physical description of these phenomena. Here we review our two-order-parameter model of liquid and consider how transient local structural ordering is linked to all of the above-mentioned phenomena. The relationship between these phenomena is also discussed.

Roles of bond orientational ordering in glass transition and crystallization

It is widely believed that crystallization in three dimensions is primarily controlled by positional ordering, and not by bond orientational ordering. In other words, bond orientational ordering is usually considered to be merely a consequence of positional ordering and thus has often been ignored. This one-order-parameter (density) description may be reasonable when we consider an equilibrium liquid-solid transition, but may not be enough to describe a metastable state and the kinetics of the transition. Here we propose that bond orientational ordering can play a key role in (i) crystallization, (ii) the ordering to quasi-crystal and (iii) vitrification, which occurs under rather weak frustration against crystallization. In a metastable supercooled state before crystallization, a system generally tends to have bond orientational order at least locally as a result of a constraint of dense packing. For a system interacting with hard-core repulsions, the constraint is intrinsically of geometrical origin and thus the basic physics is the same as nematic ordering of rod-like particles upon densification. Furthermore, positional ordering is easily destroyed even by weak frustration such as polydispersity and anisotropic interactions which favour a symmetry not consistent with that of the equilibrium crystal. Thus we may say that vitrification can be achieved by disturbing and prohibiting long-range positional ordering. Even in such a situation, bond orientational ordering still survives, accompanying its critical-like fluctuations, which are the origin of dynamic heterogeneity for this case. This scenario naturally explains both the absence of positional order and the development of bond orientational order upon cooling in a supercooled state. Although our argument is speculative in nature, we emphasize that this physical picture can coherently explain crystallization, vitrification, quasi-crystallization and their relationship in a natural manner. For a strongly frustrated system, even bond orientational order can be destroyed. Even in such a case there may still appear a structural signature of dense packing, which is linked to slow dynamics.

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.

Critical-like behaviour of glass-forming liquids

Recently it has been revealed that when approaching the glass-transition temperature, T(g), the dynamics of a liquid not only drastically slows down, but also becomes progressively more heterogeneous. From our simulations and experiments of six different glass-forming liquids, we find that the heterogeneous dynamics is a result of critical-like fluctuations of static structural order, contrary to a common belief that it is purely of dynamic origin. The static correlation length and susceptibility of a structural order parameter show Ising-like power-law divergence towards the ideal glass-transition point. However, this structural ordering accompanies little density change, which explains why it has not been detected by the static structure factor so far. Our results suggest a far more direct link than thought before between glass transition and critical phenomena. Indeed, the glass transition may be a new type of critical phenomenon where a structural order parameter is directly linked to slowness.

Low-temperature broken-symmetry phases of spiral antiferromagnets

We study Heisenberg antiferromagnets with nearest- (J1) and third- (J3) neighbor exchange on the square lattice. In the limit of spin S-->infinity, there is a zero temperature (T) Lifshitz point at J(3)=1/4J(1), with long-range spiral spin order at T=0 for J3>1/4J(1). We present classical Monte Carlo simulations and a theory for T>0 crossovers near the Lifshitz point: spin rotation symmetry is restored at any T>0, but there is a broken lattice reflection symmetry for 0< or =T<T(c) approximately (J3-1/4J(1))S2. The transition at T=T(c) is consistent with Ising universality. We also discuss the quantum phase diagram for finite S.

Melting of bosonic stripes

We use quantum Monte Carlo simulations to determine the finite temperature phase diagram and to investigate the thermal and quantum melting of stripe phases in a two-dimensional hard-core boson model. At half filling and low temperatures the stripes melt at a first order transition. In the doped system, the melting transitions of the smectic phase at high temperatures and the superfluid smectic (supersolid) phase at low temperatures are either very weakly first order, or of second order with no clear indications for an intermediate nematic phase.

Ising transition in the two-dimensional quantum J1-J2 Heisenberg model

We study the thermodynamics of the spin-S two-dimensional quantum Heisenberg antiferromagnet on the square lattice with nearest (J1) and next-nearest (J2) neighbor couplings in its collinear phase (J(2)/J(1)>0.5), using the pure-quantum self-consistent harmonic approximation. Our results show the persistence of a finite-temperature Ising phase transition for every value of the spin, provided that the ratio J(2)/J(1) is greater than a critical value corresponding to the onset of collinear long-range order at zero temperature. We also calculate the spin and temperature dependence of the collinear susceptibility and correlation length, and we discuss our results in light of the experiments on Li2VOSiO4 and related compounds.