Source author record

Gia-Wei Chern

Gia-Wei Chern appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

57works
16topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

57 published item(s)

preprint2026arXiv

Equivariant Neural Networks for Force-Field Models of Lattice Systems

Machine-learning (ML) force fields enable large-scale simulations with near-first-principles accuracy at substantially reduced computational cost. Recent work has extended ML force-field approaches to adiabatic dynamical simulations of condensed-matter lattice models with coupled electronic and structural or magnetic degrees of freedom. However, most existing formulations rely on hand-crafted, symmetry-aware descriptors, whose construction is often system-specific and can hinder generality and transferability across different lattice Hamiltonians. Here we introduce a symmetry-preserving framework based on equivariant neural networks (ENNs) that provides a general, data-driven mapping from local configurations of dynamical variables to the associated on-site forces in a lattice Hamiltonian. In contrast to ENN architectures developed for molecular systems -- where continuous Euclidean symmetries dominate -- our approach aims to embed the discrete point-group and internal symmetries intrinsic to lattice models directly into the neural-network representation of the force field. As a proof of principle, we construct an ENN-based force-field model for the adiabatic dynamics of the Holstein Hamiltonian on a square lattice, a canonical system for electron-lattice physics. The resulting ML-enabled large-scale dynamical simulations faithfully capture mesoscale evolution of the symmetry-breaking phase, illustrating the utility of lattice-equivariant architectures for linking microscopic electronic processes to emergent dynamical behavior in condensed-matter lattice systems.

preprint2026arXiv

Graph Neural Networks in the Wilson Loop Representation of Abelian Lattice Gauge Theories

Local gauge structures play a central role in a wide range of condensed matter systems and synthetic quantum platforms, where they emerge as effective descriptions of strongly correlated phases and engineered dynamics. We introduce a gauge-invariant graph neural network (GNN) architecture for Abelian lattice gauge models, in which symmetry is enforced explicitly through local gauge-invariant inputs, such as Wilson loops, and preserved throughout message passing, eliminating redundant gauge degrees of freedom while retaining expressive power. We benchmark the approach on both $\mathbb{Z}_2$ and $\mathrm{U}(1)$ lattice gauge models, achieving accurate predictions of global observables and spatially resolved quantities despite the nonlocal correlations induced by gauge-matter coupling. We further demonstrate that the learned model serves as an efficient surrogate for semiclassical dynamics in $\mathrm{U}(1)$ quantum link models, enabling stable and scalable time evolution without repeated fermionic diagonalization, while faithfully reproducing both local dynamics and statistical correlations. These results establish gauge-invariant message passing as a compact and physically grounded framework for learning and simulating Abelian lattice gauge systems.

preprint2026arXiv

Machine learning nonequilibrium phase transitions in charge-density wave insulators

Nonequilibrium electronic forces play a central role in voltage-driven phase transitions but are notoriously expensive to evaluate in dynamical simulations. Here we develop a machine learning framework for adiabatic lattice dynamics coupled to nonequilibrium electrons, and demonstrate it for a gating induced insulator to metal transition out of a charge density wave state in the Holstein model. Although exact electronic forces can be obtained from nonequilibrium Green's function (NEGF) calculations, their high computational cost renders long time dynamical simulations prohibitively expensive. By exploiting the locality of the electronic response, we train a neural network to directly predict instantaneous local electronic forces from the lattice configuration, thereby bypassing repeated NEGF calculations during time evolution. When combined with Brownian dynamics, the resulting machine learning force field quantitatively reproduces domain wall motion and nonequilibrium phase transition dynamics obtained from full NEGF simulations, while achieving orders of magnitude gains in computational efficiency. Our results establish direct force learning as an efficient and accurate approach for simulating nonequilibrium lattice dynamics in driven quantum materials.

preprint2026arXiv

Nonequilibrium DC Current Generation in a Driven Dissipative Haldane Model

The interplay of topology with nonequilibrium driving and dissipation in open quantum systems has recently attracted significant interest in condensed matter physics. In this work, we investigate a driven, dissipative Haldane model using large-scale numerical simulations of Lindblad dynamics. We show that the system evolves into a time-periodic quasi-steady state when subjected to driving and dissipation, with the ground-state topological invariant, the Chern number, no longer being quantized. Nevertheless, remnants of the underlying band topology persist in this state. To quantify this regime, we introduce an occupation-weighted Chern number that captures the topology of this nonequilibrium steady state. We further analyze charge transport in the presence of simultaneous driving and damping and demonstrate that a finite DC bulk current emerges when inversion symmetry is broken by a staggered sublattice potential. The magnitude and direction of this current are controlled by the driving amplitude, revealing a tunable nonequilibrium transport response rooted in broken symmetries and residual topology.

preprint2026arXiv

Phase-space networks and connectivity of the kagome antiferromagnet

We study the coplanar ground-state manifold of the kagome Heisenberg antiferromagnet using a phase-space network representation, in which nodes correspond to coplanar ground states and edges represent transitions generated by weathervane loop rotations. In the coplanar manifold, each configuration can be mapped to a three-coloring problem on the dual honeycomb lattice, where a weathervane mode corresponds to a closed loop of two alternating colors. By comparing networks that include all weathervane loops with networks restricted to elementary six-spin loops, we examine how energetic constraints shape phase-space structure. We find that connectivity distributions are sharply peaked in large systems, while restrictions to short loops reduce typical connectivity. Spectral properties further distinguish the two cases, with short-loop networks exhibiting Gaussian spectra and full networks displaying non-Gaussian features associated with correlated loop updates. Finally, a box-counting analysis reveals distinct fractal properties of the two networks, demonstrating how energetic constraints control the global geometry of configuration space. These results show that the hierarchy of weathervane loop rotations provides a direct link between microscopic constraints and emergent phase-space geometry in a frustrated magnet.

preprint2026arXiv

Pseudospin Formulation of Quench Dynamics in the Semiclassical Holstein Model

We present a pseudospin formulation for the post-quench dynamics of charge-density-wave (CDW) order in the half-filled spinless Holstein model on a square lattice, assuming spatially homogeneous evolution. This Anderson pseudospin description captures the coherent nonequilibrium dynamics of the coupled electron-lattice system. Numerical simulations reveal three distinct dynamical regimes of the CDW order parameter following a quench-locked oscillations, Landau-damped dynamics, and overdamped relaxation-closely paralleling quench dynamics in BCS superconductors and other electronically driven symmetry-breaking phases. Crucially, however, the presence of dynamical lattice degrees of freedom leads to qualitatively different long-time behavior. In particular, while the oscillation amplitude is reduced in the damped regimes, CDW oscillations do not fully decay but instead persist indefinitely due to feedback from the lattice field. We further show that these persistent oscillations are characterized by a nonequilibrium electronic distribution, which provides an intuitive understanding of both their amplitude and the renormalization of the oscillation frequency relative to the bare Holstein phonon frequency. Our results highlight the essential role of lattice dynamics in nonequilibrium ordered phases and establish a clear distinction between electron-lattice-driven CDW dynamics and their purely electronic counterparts.

preprint2026arXiv

Recurrent convolutional neural networks for modeling non-adiabatic dynamics of quantum-classical systems

Recurrent neural networks (RNNs) have recently been extensively applied to model the time-evolution in fluid dynamics, weather predictions, and even chaotic systems thanks to their ability to capture temporal dependencies and sequential patterns in data. Here we present a RNN model based on convolution neural networks for modeling the nonlinear non-adiabatic dynamics of hybrid quantum-classical systems. The dynamical evolution of the hybrid systems is governed by equations of motion for classical degrees of freedom and von Neumann equation for electrons. The physics-aware recurrent convolution (PARC) neural network structure incorporates a differentiator-integrator architecture that inductively models the spatiotemporal dynamics of generic physical systems. We apply our RNN approach to learn the space-time evolution of a one-dimensional semi-classical Holstein model after an interaction quench. For shallow quenches (small changes in electron-lattice coupling), the deterministic dynamics can be accurately captured using a single-CNN-based recurrent network. In contrast, deep quenches induce chaotic evolution, making long-term trajectory prediction significantly more challenging. Nonetheless, we demonstrate that the PARC-CNN architecture can effectively learn the statistical climate of the Holstein model under deep-quench conditions.

preprint2022arXiv

Anisotropic MagnetoMemristance

In the last decade, nanoscale resistive devices with memory have been the subject of intense study because of their possible use in brain-inspired computing. However, operational endurance is one of the limiting factors in the adoption of such technology. For this reason, we discuss the emergence of current-induced memristance in magnetic materials, known for their durability. We show analytically and numerically that a single ferromagnetic layer can possess GHz memristance, due to a combination of two factors: a current-induced transfer of angular momentum (Zhang-Li torque) and the anisotropic magnetoresistance (AMR). We term the resulting effect the anisotropic magneto-memristance (AMM). We connect the AMM to the topology of the magnetization state, within a simple model of a 1-dimensional annulus-shaped magnetic layer, confirming the analytical results with micromagnetic simulations for permalloy. Our results open a new path towards the realization of single-layer magnetic memristive devices operating at GHz frequencies.

preprint2022arXiv

Atomistic simulation of Mott transition in fluid metals: Combining molecular dynamics with dynamical mean-field theory

We present a new quantum molecular dynamics (MD) method where the electronic structure and atomic forces are solved by a real-space dynamical mean-field theory (DMFT). Contrary to most quantum MD methods that are based on effective single-particle wave functions, the DMFT approach is able to describe correlation-induced Mott metal-insulator transitions and the associated incoherent electronic excitations in an atomic liquid. We apply the DMFT-MD method to study Mott transitions in an atomic liquid model which can be viewed as the liquid-state generalization of the Hubbard model. The half-filled Hubbard liquids also provide a minimum model for alkali fluid metals. Our simulations uncover two distinct types of Mott transition depending on the atomic bonding and short-range structures in the electronically delocalized phase. In the first scenario where atoms tend to form dimers, increasing the Hubbard repulsion gives rise to a transition from a molecular insulator to an atomic insulator with a small window of enhanced metallicity in the vicinity of the localization transition. On the other hand, for Hubbard liquids with atoms forming large conducting clusters, the localization of electrons leads to the fragmentation of clusters and is intimately related to the liquid-gas transition of atoms. Implications of our results for metal-insulator transitions in fluid alkali metals are discussed.

preprint2022arXiv

Descriptors for Machine Learning Model of Generalized Force Field in Condensed Matter Systems

We outline the general framework of machine learning (ML) methods for multi-scale dynamical modeling of condensed matter systems, and in particular of strongly correlated electron models. Complex spatial temporal behaviors in these systems often arise from the interplay between quasi-particles and the emergent dynamical classical degrees of freedom, such as local lattice distortions, spins, and order-parameters. Central to the proposed framework is the ML energy model that, by successfully emulating the time-consuming electronic structure calculation, can accurately predict a local energy based on the classical field in the intermediate neighborhood. In order to properly include the symmetry of the electron Hamiltonian, a crucial component of the ML energy model is the descriptor that transforms the neighborhood configuration into invariant feature variables, which are input to the learning model. A general theory of the descriptor for the classical fields is formulated, and two types of models are distinguished depending on the presence or absence of an internal symmetry for the classical field. Several specific approaches to the descriptor of the classical fields are presented. Our focus is on the group-theoretical method that offers a systematic and rigorous approach to compute invariants based on the bispectrum coefficients. We propose an efficient implementation of the bispectrum method based on the concept of reference irreducible representations. Finally, the implementations of the various descriptors are demonstrated on well-known electronic lattice models.

preprint2022arXiv

Enhancement of Atomic Diffusion due to Electron Delocalization in Fluid Metals

We present a general theory of atomic self-diffusion in the vicinity of a Mott metal-insulator transition in fluid metals. Upon decreasing the electron correlation from the Mott insulating phase, the delocalization of electrons gives rise to an increasing attractive interatomic interaction, which is expected to introduce an additional friction, hence reducing the atomic diffusivity. Yet, our quantum molecular dynamics simulations find an intriguing enhancement of the diffusion coefficient induced by the emerging attractive force. We show that this counterintuitive phenomenon results from the reduction of the repulsive core and the suppression of the attractive tail by thermal fluctuations. The proposed scenario is corroborated by the Chapman-Enskog theory and classical molecular dynamics simulations on a standard liquid model based on the Morse potential. Our work not only provides a general mechanism of the attraction-facilitated diffusion enhancement in simple liquids, but also sheds new lights on the nontrivial effects of electron correlation on atomic dynamics.

preprint2022arXiv

Evidence for pressure induced unconventional quantum criticality in the coupled spin ladder antiferromagnet C$_9$H$_{18}$N$_2$CuBr$_4$

Quantum phase transitions in quantum matter occur at zero temperature between distinct ground states by tuning a nonthermal control parameter. Often, they can be accurately described within the Landau theory of phase transitions, similarly to conventional thermal phase transitions. However, this picture can break down under certain circumstances. Here, we present a comprehensive study of the effect of hydrostatic pressure on the magnetic structure and spin dynamics of the spin-1/2 ladder compound C$_9$H$_{18}$N$_2$CuBr$_4$. Single-crystal heat capacity and neutron diffraction measurements reveal that the N$\rm \acute{e}$el-ordered phase breaks down beyond a critical pressure of $P_{\rm c}$$\sim$1.0 GPa through a continuous quantum phase transition. Estimates of the critical exponents suggest that this transition may fall outside the traditional Landau paradigm. The inelastic neutron scattering spectra at 1.3 GPa are characterized by two well-separated gapped modes, including one continuum-like and another resolution-limited excitation in distinct scattering channels, which further indicates an exotic quantum-disordered phase above $P_{\rm c}$.

preprint2022arXiv

Machine learning predictions for local electronic properties of disordered correlated electron systems

We present a scalable machine learning (ML) model to predict local electronic properties such as on-site electron number and double occupation for disordered correlated electron systems. Our approach is based on the locality principle, or the nearsightedness nature, of many-electron systems, which means local electronic properties depend mainly on the immediate environment. A ML model is developed to encode this complex dependence of local quantities on the neighborhood. We demonstrate our approach using the square-lattice Anderson-Hubbard model, which is a paradigmatic system for studying the interplay between Mott transition and Anderson localization. We develop a lattice descriptor based on group-theoretical method to represent the on-site random potentials within a finite region. The resultant feature variables are used as input to a multi-layer fully connected neural network, which is trained from datasets of variational Monte Carlo (VMC) simulations on small systems. We show that the ML predictions agree reasonably well with the VMC data. Our work underscores the promising potential of ML methods for multi-scale modeling of correlated electron systems.

preprint2021arXiv

Distribution of Kinks in an Ising Ferromagnet After Annealing and the Generalized Kibble-Zurek Mechanism

We consider the annealing dynamics of a one-dimensional Ising ferromagnet induced by a temperature quench in finite time. In the limit of slow cooling, the asymptotic two-point correlator is analytically found under Glauber dynamics, and the distribution of the number of kinks in the final state is shown to be consistent with a Poissonian distribution. The mean kink number, the variance, and the third centered moment take the same value and obey a universal power-law scaling with the quench time in which the temperature is varied. The universal power-law scaling of cumulants is corroborated by numerical simulations based on Glauber dynamics for moderate cooling times away from the asymptotic limit, when the kink-number distribution takes a binomial form. We analyze the relation of these results to physics beyond the Kibble-Zurek mechanism for critical dynamics, using the kink number distribution to assess adiabaticity and its breakdown. We consider linear, nonlinear, and exponential cooling schedules, among which the latter provides the most efficient shortcuts to cooling in a given quench time. The non-thermal behavior of the final state is established by considering the trace norm distance to a canonical Gibbs state.

preprint2021arXiv

Kinetics of thermal Mott transitions in the Hubbard model

We present the first-ever multi-scale dynamical simulation of the temperature-controlled Mott metal-insulator transition in the Hubbard model. By integrating advanced electronic structure method and an efficient Gutzwiller/slave-boson solver into molecular dynamics simulations, we demonstrate that the transformation from a correlated metal to the Mott insulating phase proceeds via the nucleation and growth of the Mott droplets. Moreover, the time evolution of the Mott volume fraction is found to follow a universal transformation kinetics. We show that after an initial incubation period, the early stage of the phase transformation is characterized by a constant nucleation rate and an interface-controlled cluster growth mechanism, consistent with the classical theory developed by Kolmogorov, Johnson, Mehl, and Avrami. This is followed by a novel intermediate stage of accelerated phase transformation that is significantly different from the prediction of the classical theory. Moreover, the cluster-growth dynamics in this intermediate stage exhibits an unexpected avalanche behavior, similar to the Barkhausen noise in magnetization dynamics, even in the absence of quenched disorder. Detailed structural characterization further uncovers a universal correlation function for the transient mixed-phase states of the Mott transition. We also discuss implications of our findings for spatially resolved measurements of Mott metal-insulator transition obtained in recent nano-imaging experiments.

preprint2020arXiv

Machine learning dynamics of phase separation in correlated electron magnets

We demonstrate machine-learning enabled large-scale dynamical simulations of electronic phase separation in double-exchange system. This model, also known as the ferromagnetic Kondo lattice model, is believed to be relevant for the colossal magnetoresistance phenomenon. Real-space simulations of such inhomogeneous states with exchange forces computed from the electron Hamiltonian can be prohibitively expensive for large systems. Here we show that linear-scaling exchange field computation can be achieved using neural networks trained by datasets from exact calculation on small lattices. Our Landau-Lifshitz dynamics simulations based on machine-learning potentials nicely reproduce not only the nonequilibrium relaxation process, but also correlation functions that agree quantitatively with exact simulations. Our work paves the way for large-scale dynamical simulations of correlated electron systems using machine-learning models.

preprint2019arXiv

Emergent Snake Magnetic Domains in Canted Kagome Ice

We study the two-dimensional kagome-ice model derived from a pyrochlore lattice with second- and third-neighbor interactions. The canted moments align along the local $\langle 111 \rangle$ axes of the pyrochlore and respond to both in-plane and out-of-plane external fields. We find that the combination of further-neighbor interactions together with the external fields introduces a rich phase diagram with different spin textures. Close to the phase boundaries, metastable $\textit{"snake"}$ domains emerge with extremely long relaxation time. Our kinetic Monte Carlo analysis of the magnetic-field quench process from saturated state shows unusually slow dynamics. Despite that the interior spins are almost frozen in snake domains, the spins on the edge are free to fluctuate locally, leading to frequent creation and annihilation of monopole-anti-monopole bound states. Once the domains are formed, these excitations are localized and can hardly propagate due to the energy barrier of snakes. The emergence of such snake domains may shed light on the experimental observation of dipolar spin ice under tilted fields, and provide a new strategy to manipulate both spin and charge textures in artificial spin ice.

preprint2016arXiv

Magnetic avalanches in granular ferromagnets: Thermal activated collective behavior

We present a numerical study on the thermal activated avalanche dynamics in granular materials composed of ferromagnetic clusters embedded in a non-magnetic matrix. A microscopic dynamical simulation based on the reaction-diffusion process is developed to modeling the magnetization process of such systems. The large-scale simulations presented here explicitly demonstrate inter-granular collective behavior induced by thermal activation of spin tunneling. In particular, we observe an intriguing criticality controlled by the rate of energy dissipation. We show that thermal activated avalanches can be understood in the framework of continuum percolation and the emergent dissipation induced criticality is in the universality class of three-dimensional percolation transition. Implications of these results to the phase-separated states of colossal magnetoresistance materials and other artificial granular magnetic systems are also~discussed.

preprint2016arXiv

Matter-wave propagation in optical lattices: geometrical and flat-band effects

The geometry of optical lattices can be engineered allowing the study of atomic transport along paths arranged in patterns that are otherwise difficult to probe in the solid state. A question readily accessible to atomic systems is related to the speed of propagation of matter-waves as a function of the lattice geometry. To address this issue, we have investigated theoretically the quantum transport of non-interacting and weakly-interacting ultracold fermionic atoms in several 2D optical lattice geometries. We find that the triangular lattice has a higher propagation velocity compared to the square lattice, despite supporting longer paths. The body-centered square lattice has even longer paths, nonetheless the propagation velocity is yet faster. This apparent paradox arises from the mixing of the momentum states which leads to different group velocities in quantum systems. Standard band theory provides an explanation and allows for a systematic way to search and design systems with controllable matter-wave propagation. Moreover, the presence of a flat band such as in a two-leg ladder geometry leads to a dynamical density discontinuity, which contrasts the behavior of mobile and localized atoms in quantum transport. Our predictions are realizable with present experimental capability.

preprint2016arXiv

Resistivity Minimum in Highly Frustrated Itinerant Magnets

We study the transport properties of frustrated itinerant magnets comprising localized {\it classical} moments, which interact via exchange with the conduction electrons. Strong frustration stabilizes a liquidlike spin state, which extends down to temperatures well below the effective Ruderman-Kittel-Kasuya-Yosida interaction scale. The crossover into this state is characterized by spin structure factor enhancement at wave vectors smaller than twice the Fermi wave vector magnitude. The corresponding enhancement of electron scattering generates a resistivity upturn at decreasing temperatures.

preprint2015arXiv

Direct Visualization of Memory Effects in Artificial Spin Ice

We experimentally demonstrate that arrays of interacting nanoscale ferromagnetic islands, known as artificial spin ice, develop reproducible microstates upon cycling an applied magnetic field. The onset of this memory effect is determined by the strength of the applied field relative to the array coercivity. Specifically, when the applied field strength is almost exactly equal to the array coercivity, several training cycles are required before the array achieves a nearly completely repeatable microstate, whereas when the applied field strength is stronger or weaker than the array coercivity, a repeatable microstate is achieved after the first minor loop. We show through experiment and simulation that this memory exhibited by artificial spin ice is due to a ratchet effect on interacting, magnetically-charged defects in the island moment configuration and to the complexity of the network of strings of reversed moments that forms during magnetization reversal.

preprint2015arXiv

PT-symmetric phase in kagome photonic lattices

Kagome lattice is a two-dimensional network of corner-sharing triangles and is often associated with geometrical frustration. In particular, the frustrated coupling between waveguide modes in a kagome array leads to a dispersionless flat band consisting of spatially localized modes. Here we propose a complex photonic lattice by placing $\mathcal{PT}$-symmetric dimers at the kagome lattice points. Each dimer corresponds to a pair of strongly coupled waveguides. With balanced arrangement of gain and loss on individual dimers, the system exhibits a $\mathcal{PT}$-symmetric phase for finite gain/loss parameter up to a critical value. The beam evolution in this complex kagome waveguide array exhibits a novel oscillatory rotation of optical power along the propagation distance. Long-lived local chiral structures originating from the nearly flat bands of the kagome structure are observed when the lattice is subject to a narrow beam excitation.

preprint2015arXiv

Topological defects as relics of emergent continuous symmetry and Higgs condensation of disorder in ferroelectrics

Lars Onsager and Richard Feynman envisioned that the three-dimensional (3D) superfluid-to-normal $λ$ transition in $^{4}$He occurs through the proliferation of vortices. This process should hold for every phase transition in the same universality class. The role of topological defects in symmetry-breaking phase transitions has become a prime topic in cosmology and high-temperature superconductivity, even though direct imaging of these defects is challenging. Here we show that the U(1) continuous symmetry that emerges at the ferroelectric critical point of multiferroic hexagonal manganites leads to a similar proliferation of vortices. Moreover, the disorder field (vortices) is coupled to an emergent U(1) gauge field, which becomes massive by means of the Higgs mechanism when vortices condense (span the whole system) upon heating above the ferroelectric transition temperature. Direct imaging of the vortex network in hexagonal manganites offers unique experimental access to this dual description of the ferroelectric transition, while enabling tests of the Kibble-Zurek mechanism.

preprint2015arXiv

Vortex Crystals with Chiral Stripes in Itinerant Magnets

We study noncoplanar magnetic ordering in frustrated itinerant magnets. For a family of Kondo square lattice models with classical local moments, we find that a double-$Q$ noncoplanar vortex crystal has lower energy than the single-$Q$ helical order expected from the Ruderman-Kittel-Kasuya-Yosida interaction whenever the lattice symmetry dictates four global maxima in the bare magnetic susceptibility. By expanding in the small Kondo exchange and the degree of noncoplanarity, we demonstrate that this noncoplanar state arises from a Fermi surface instability, and it is generic for a wide range of electron filling fractions whenever the two ordering wave vectors connect independent sections of the Fermi surface.

preprint2014arXiv

Avalanches and disorder-induced criticality in artificial spin ices

We show that both square and kagome artificial spin ice systems exhibit disorder-induced nonequilibrium phase transitions, with power law avalanche distributions at the critical disorder level. The different nature of geometrical frustration in the two lattices produces distinct types of critical avalanche behavior. For the square ice, the avalanches involve the propagation of locally stable domain walls separating the two polarized ground states, and the scaling collapse agrees with an interface depinning mechanism. In contrast, avalanches in the fully frustrated kagome ice exhibit pronounced branching behaviors that resemble those found in directed percolation. The kagome ice also shows an interesting crossover in the power-law scaling of the avalanches at low disorder. Our results show that artificial spin ices are ideal systems in which to study nonequilibrium critical point phenomena.

preprint2014arXiv

Dynamically generated flat-band phases in optical kagome lattices

Motivated by recent advances in the realization of complex two-dimensional optical lattices, we investigate theoretically the quantum transport of ultracold fermions in an optical kagome lattice. In particular, we focus on its extensively degenerate localized states (flat band). By loading fermions in a partial region of the lattice and depleting the mobile atoms at the far boundary of the initially unoccupied region, we find a dynamically generated flat-band insulator, which is also a population-inverted state. We further show that inclusion of weak repulsion leads to a dynamical stripe phase for two-component fermions in a similar setup. Finally, by preparing a topological insulating state in a partially occupied kagome lattice, we find that the topological chiral current decays but exhibits an interesting oscillating dynamics during the nonequilibrium transport. Given the broad variety of lattice geometries supporting localized or topological states, our work suggests new possibilities to use geometrical effects and their dynamics in atomtronic devices.

preprint2014arXiv

Novel magnetic orderings in the kagome Kondo-lattice model

We consider the Kondo-lattice model on the kagome lattice and study its weak-coupling instabilities at band filling fractions for which the Fermi surface has singularities. These singularites include Dirac points, quadratic Fermi points in contact with a flat band, and Van Hove saddle points. By combining a controlled analytical approach with large-scale numerical simulations, we demonstrate that the weak-coupling instabilities of the Kondo-lattice model lead to exotic magnetic orderings. In particular, some of these magnetic orderings produce a spontaneous quantum anomalous Hall state.

preprint2014arXiv

Quantum Hall ice

We show that the chiral kagome ice manifold exhibits an anomalous integer quantum Hall effect (IQHE) when coupled to itinerant electrons. Although electron-mediated interactions select a magnetically ordered ground state, the full ice manifold can coexist with the IQHE over a range of finite temperatures. The degenerate ice states provide a natural realization of power-law correlated flux disorder, for which the spectral gap of the system remains robust. The quantized (up to exponentially small finite-temperature corrections) Hall conductance persists over a wide range of electron densities due to the disorder-induced localization of electronic states.

preprint2014arXiv

Raman Scattering Signatures of Kitaev Spin Liquids in A$_2$ IrO$_3$ Iridates

We study theoretically the Raman scattering response $I(ω)$ in the gapless quantum spin liquid phase of the Kitaev-Heisenberg model. The dominant polarization-independent contribution $I_K (ω)$ reflects the density of states of the emergent Majorana fermions in the ground-state flux-sector. The integrability-breaking Heisenberg exchange generates a second contribution, whose dominant part $I_H (ω)$ has the form of a quantum quench corresponding to an abrupt insertion of four $Z_2$ gauge fluxes. This results in a weakly polarization dependent response with a sharp peak at the energy of the flux excitation accompanied by broad features, which can be related to Majorana fermions in the presence of the perturbed gauge field. We discuss the experimental situation and explore more generally the influence of integrability breaking for Kitaev spin liquid response functions.

preprint2014arXiv

Stiffness from Disorder in Triangular-Lattice Ising Thin Films

We study the triangular lattice Ising model with a finite number of vertically stacked layers and demonstrate a low temperature reentrance of two Berezinskii-Kosterlitz-Thouless transitions, which results in an extended disordered regime down to $T=0$. Numerical results are complemented with the derivation of an effective low-temperature dimer theory. Contrary to order by disorder, we present a new scenario for fluctuation-induced ordering in frustrated spin systems. While short-range spin-spin correlations are enhanced by fluctuations, quasi-long-range ordering is precluded at low enough temperatures by proliferation of topological defects.

preprint2013arXiv

Chiral magnetism and helimagnons in a pyrochlore antiferromagnet

Recent neutron scattering measurements on the spinel CdCr2O4 revealed a rare example of helical magnetic order in geometrically frustrated pyrochlore antiferromagnet. The spin spiral characterized by an incommensurate wavevector Q = 2pi (0, delta, 1) with delta ~0.09 is accompanied by a tetragonal distortion. Here we conduct a systematic study on the magnetic ground state resulting from the interplay between the Dzyaloshinskii-Moriya interaction and further neighbor exchange couplings, two of the most important mechanisms for stabilizing incommensurate spin orders. We compute the low-energy spin-wave spectrum based on a microscopic spin Hamiltonian and find a dispersion relation characteristic of the helimagnons. By numerically integrating the Landau-Lifshitz-Gilbert equation with realistic model parameters, an overall agreement between experiment and the numerical spectrum, lending further support to the view that a softened optical phonon triggers the magnetic transition and endows the lattice a chirality.

preprint2013arXiv

Degeneracy and criticality from emergent frustration in artificial spin ice

Although initially introduced to mimic the spin-ice pyrochlores, no artificial spin ice has yet exhibited the expected degenerate ice-phase with critical correlations similar to the celebrated Coulomb phase in the pyrochlore lattice. Here we study a novel artificial spin ice based on a vertex-frustrated rather than pairwise frustrated geometry and show that it exhibits a quasicritical ice phase of extensive residual entropy and, significantly, algebraic correlations. Interesting in its own regard as a novel realization of frustration in a vertex system, our lattice opens new pathways to study defects in a critical manifold and to design degeneracy in artificial magnetic nanoarrays, a task so far elusive.

preprint2013arXiv

Dipolar order by disorder in the classical Heisenberg antiferromagnet on the kagome lattice

Ever since the experiments which founded the field of highly frustrated magnetism, the kagome Heisenberg antiferromagnet has been the archetypical setting for the study of fluctuation induced exotic ordering. To this day the nature of its classical low-temperature state has remained a mystery: the non-linear nature of the fluctuations around the exponentially numerous harmonically degenerate ground states has not permitted a controlled theory, while its complex energy landscape has precluded numerical simulations at low temperature. Here we present an efficient Monte Carlo algorithm which removes the latter obstacle. Our simulations detect a low-temperature regime in which correlations saturate at a remarkably small value. Feeding these results into an effective model and analyzing the results in the framework of an appropriate field theory implies the presence of long-range dipolar spin order with a tripled unit cell.

preprint2013arXiv

Four-coloring model and frustrated superfluidity in the diamond lattice

We propose a novel four-coloring model which describes "frustrated superfluidity" of $p$-band bosons in the diamond optical lattice. The superfluid phases of the condensate wavefunctions on the diamond-lattice bonds are mapped to four distinct colors at low temperatures. The fact that a macroscopic number of states satisfy the constraints that four differently colored bonds meet at the same site leads to an extensive degeneracy in the superfluid ground state at the classical level. We demonstrate that the phase of the superfluid wavefunction as well as the orbital angular momentum correlations exhibit a power-law decay in the degenerate manifold that is described by an emergent magnetostatic theory with three independent flux fields. Our results thus provide a novel example of critical superfluid phase with algebraic order in three dimensions. We further show that quantum fluctuations favor a Néel ordering of orbital angular moments with broken sublattice symmetry through the order-by-disorder mechanism.

preprint2013arXiv

Frustrated colloidal ordering and fully packed loops in arrays of optical traps

We propose that a system of colloidal particles interacting with a honeycomb array of optical traps that each contain three wells can be used to realize a fully packed loop model. One of the phases in this system can be mapped to Baxter's three-coloring problem, offering an easily accessible physical realization of this problem. As a function of temperature and interaction strength, we find a series of phases, including long range ordered loop or stripe states, stripes with sliding symmetries, random packed loop states, and disordered states in which the loops break apart. Our geometry could be constructed using ion trap arrays, BEC vortices in optical traps, or magnetic vortices in nanostructured superconductors.

preprint2013arXiv

Helimagnons in a chiral ground state of the pyrochlore antiferromagnets

The Goldstone mode in a helical magnetic phase, also known as the helimagnon, is a propagating mode with a highly anisotropic dispersion relation. Here we study theoretically the helimagnon excitations in a complex chiral ground state of pyrochlore antiferromagnets such as spinel CdCr2O4 and itinerant magnet YMn2. We show that the effective theory of the soft modes in the helical state possesses a symmetry similar to that of smectic liquid crystals. We compute the low-energy spin-wave spectrum based on a microscopic spin Hamiltonian and find a dispersion relation characteristic of the helimagnons. By performing dynamics simulations with realistic model parameters, we also obtain an overall agreement between experiment and the numerical spin-wave spectrum. Our work thus also clarifies the mechanisms that relive the magnetic frustration in these compounds.

preprint2013arXiv

Magnetic Field Induced Transition in Vanadium Spinels

We study vanadium spinels $A$V$_2$O$_4$ ($A$ = Cd, Mg) in pulsed magnetic fields up to 65 T. A jump in magnetization at $μ_0 H \approx$ 40 T is observed in the single-crystal MgV$_2$O$_4$, indicating a field induced quantum phase transition between two distinct magnetic orders. In the multiferroic CdV$_2$O$_4$, the field-induced transition is accompanied by a suppression of the electric polarization. By modeling the magnetic properties in the presence of strong spin-orbit coupling characteristic of vanadium spinels, we show that both features of the field-induced transition can be successfully explained by including the effects of the local trigonal crystal field.

preprint2013arXiv

Magnetic monopole polarons in spin ice with mixed coordination numbers

Emergent quasiparticles that arise from the fractionalization of the microscopic degrees of freedom have been one of the central themes in modern condensed matter physics. The notion of magnetic monopoles, freely moving quasiparticles fragmented from local dipole excitations, has enjoyed much success in understanding the thermodynamic, static, and transport properties of the so-called spin-ice materials. The artificial version of spin ice, where a lattice of nanoscale magnetic dipoles is sculpted out of a ferromagnetic film, provides a unique opportunity to study these unusual quasiparticles in a material-by-design approach. Here we show that the elementary excitations in the ice phase of a nano-magnetic array arranged in the pentagonal lattice are composite objects comprised of the emergent monopole and a surrounding cloud of opposite uncompensated magnetic charges.

preprint2013arXiv

Realizing three-dimensional artificial spin ice by stacking planar nanoarrays

Artificial spin ice is a frustrated magnetic two-dimensional nano-material, recently employed to study variety of tailor-designed unusual collective behaviours. Recently proposed extensions to three dimensions are based on self-assembly techniques and allow little control over geometry and disorder. We present a viable design for the realization of a three-dimensional artificial spin ice with the same level of precision and control allowed by lithographic nano-fabrication of the popular two-dimensional case. Our geometry is based on layering already available two-dimensional artificial spin ice and leads to an arrangement of ice-rule-frustrated units which is topologically equivalent to that of the tetrahedra in a pyrochlore lattice. Consequently, we show, it exhibits a genuine ice phase and its excitations are, as in natural spin ice materials, magnetic monopoles interacting via Coulomb law.

preprint2013arXiv

Universal Renyi mutual information in classical systems: the case of kagome ice

We study the Renyi mutual information of classical systems characterized by a transfer matrix. We first establish a general relationship between the Renyi mutual information of such classical mixtures of configuration states, and the Renyi entropy of a corresponding Rokhsar-Kivelson--type quantum superposition. We then focus on chiral and nonchiral kagome-ice systems, classical spin liquids on the kagome lattice, which respectively have critical and short-range correlations. Through a mapping of the chiral kagome ice to the quantum Liftshitz critical field theory, we predict a universal subleading term in the Renyi mutual information of this classical spin liquid, which can be realized in the pyrochlore spin ice in a magnetic field. We verify our prediction with direct numerical transfer-matrix computations, and further demonstrate that the nonchiral kagome ice (and the corresponding quantum Rokhsar-Kivelson superposition) is a topologically trivial phase. Finally, we argue that the universal term in the mutual information of the chiral kagome ice is fragile against the presence of defects.

preprint2012arXiv

Broken translational symmetry in an emergent paramagnetic phase of graphene

We show that the spin-density wave state on the partially filled honeycomb and triangular lattices is preempted by a paramagnetic phase that breaks an emergent Z_4 symmetry of the system, associated with the four inequivalent arrangements of spins in the quadrupled unit cell. Unlike other emergent paramagnetic phases in itinerant and localized-spin systems, this state preserves the rotational symmetry of the lattice but breaks its translational symmetry, giving rise to a super-lattice structure that can be detected by scanning tunneling microscopy. This emergent phase also has distinctive signatures in the magnetic spectrum that can be probed experimentally.

preprint2012arXiv

Itinerant half-metal spin-density-wave state on the hexagonal lattice

We consider electrons on a honeycomb or triangular lattice doped to the saddle point of the bandstructure. We assume system parameters are such that spin density wave (SDW) order emerges below a temperature $T_N$ and investigate the nature of the SDW phase. We argue that at $T < T_N$ the system develops a uniaxial SDW phase whose ordering pattern breaks $O(3) \times Z_4$ symmetry and corresponds to an eight site unit cell with non-uniform spin moments on different sites. This state is a half-metal -- it preserves full original Fermi surface, but has gapless charged excitations in one spin branch only. It allows for electrical control of spin currents and is desirable for nano-science.

preprint2012arXiv

Magnetic charge and ordering in kagome spin ice

We present a numerical study of magnetic ordering in spin ice on kagome, a two-dimensional lattice of corner-sharing triangles. The magnet has six ground states and the ordering occurs in two stages, as one might expect for a six-state clock model. In spin ice with short-range interactions up to second neighbors, there is an intermediate critical phase separated from the paramagnetic and ordered phases by Kosterlitz-Thouless transitions. In dipolar spin ice, the intermediate phase has long-range order of staggered magnetic charges. The high and low-temperature phase transitions are of the Ising and 3-state Potts universality classes, respectively. Freeze-out of defects in the charge order produces a very large spin correlation length in the intermediate phase. As a result of that, the lower-temperature transition appears to be of the Kosterlitz-Thouless type.

preprint2012arXiv

Orbital disorder induced by charge fluctuations in vanadium spinels

Motivated by recent experiments on vanadium spinels, $A$V$_2$O$_4$, that show an increasing degree of electronic delocalization for smaller cation sizes, we study the evolution of orbital ordering (OO) between the strong and intermediate-coupling regimes of a multi-orbital Hubbard Hamiltonian. The underlying magnetic ordering of the Mott insulating state leads to a rapid suppression of OO due to enhanced charge fluctuations along ferromagnetic bonds. Orbital double-occupancy is rather low at the transition point indicating that the system is in the crossover region between strong and intermediate-coupling regimes when the orbital degrees of freedom become disordered.

preprint2012arXiv

Spontaneous quantum Hall effect via thermally induced quadratic Fermi point

Gapless electronic systems containing topologically nontrivial Fermi points are sources of various topological insulators. Whereas most of these special band-crossing points are built in the electronic structure of the non-interacting lattice models, we show that a quadratic Fermi point characterized by a non-zero winding number emerges with a collinear triple-Q spin-density-wave state that arises from a perfectly nested but topologically trivial Fermi surface. We obtain a universal low-energy Hamiltonian for the quadratic Fermi point and show that such collinear orderings are unstable against the onset of scalar spin chirality that opens a gap and induces a spontaneous quantum Hall insulator as the temperature tends to zero.

preprint2011arXiv

Orbital ice: an exact Coulomb phase on the diamond lattice

We demonstrate the existence of orbital Coulomb phase as the exact ground state of p-orbital exchange Hamiltonian on the diamond lattice. The Coulomb phase is an emergent state characterized by algebraic dipolar correlations and a gauge structure resulting from local constraints (ice rules) of the underlying lattice models. For most ice models on the pyrochlore lattice, these local constraints are a direct consequence of minimizing the energy of each individual tetrahedron. On the contrary, the orbital ice rules are emergent phenomena resulting from the quantum orbital dynamics. We show that the orbital ice model exhibits an emergent geometrical frustration by mapping the degenerate quantum orbital ground states to the spin-ice states obeying the 2-in-2-out constraints on the pyrochlore lattice. We also discuss possible realization of the orbital ice model in optical lattices with p-band fermionic cold atoms.

preprint2011arXiv

Quantum phase transitions in a strongly entangled spin-orbital chain: A field-theoretical approach

Motivated by recent experiments on quasi-1D vanadium oxides, we study quantum phase transitions in a one-dimensional spin-orbital model describing a Haldane chain and a classical Ising chain locally coupled by the relativistic spin-orbit interaction. By employing a field-theoretical approach, we analyze the topology of the ground-state phase diagram and identify the nature of the phase transitions. In the strong coupling limit, a long-range Néel order of entangled spin and orbital angular momentum appears in the ground state. We find that, depending on the relative scales of the spin and orbital gaps, the linear chain follows two distinct routes to reach the Néel state. First, when the orbital exchange is the dominating energy scale, a two-stage ordering takes place in which the magnetic transition is followed by melting of the orbital Ising order; both transitions belong to the two-dimensional Ising universality class. In the opposite limit, the low-energy orbital modes undergo a continuous reordering transition which represents a line of Gaussian critical points. On this line the orbital degrees of freedom form a Tomonaga-Luttinger liquid. We argue that the emergence of the Gaussian criticality results from merging of the two Ising transitions in the strong hybridization region where the characteristic spin and orbital energy scales become comparable. Finally, we show that, due to the spin-orbit coupling, an external magnetic field acting on the spins can induce an orbital Ising transition.

preprint2011arXiv

Spin superstructure and noncoplanar ordering in metallic pyrochlore magnets with degenerate orbitals

We study double-exchange models with itinerant t2g electrons in spinel and pyrochlore crystals. In both cases the localized spins form a network of corner-sharing tetrahedra. We show that the strong directional dependence of t2g orbitals leads to unusual Fermi surfaces that induce spin superstructures and noncoplanar orderings for a weak coupling between itinerant electrons and localized spins. Implications of our results to ZnV2O4 and Cd2Os2O7 are also discussed.

preprint2011arXiv

Topology of Entanglement Evolution of Two Qubits

The dynamics of a two-qubit system is considered with the aim of a general categorization of the different ways in which entanglement can disappear in the course of the evolution, e.g., entanglement sudden death. The dynamics is described by the function ~n(t), where ~n is the 15-dimensional polarization vector. This representation is particularly useful because the components of ~n are direct physical observables, there is a meaningful notion of orthogonality, and the concurrence C can be computed for any point in the space. We analyze the topology of the space S of separable states (those having C = 0), and the often lower-dimensional linear dynamical subspace D that is characteristic of a specific physical model. This allows us to give a rigorous characterization of the four possible kinds of entanglement evolution. Which evolution is realized depends on the dimensionality of D and of D \cap S, the position of the asymptotic point of the evolution, and whether or not the evolution is "distance-Markovian", a notion we define. We give several examples to illustrate the general principles, and to give a method to compute critical points. We construct a model that shows all four behaviors.

preprint2011arXiv

Two-stage ordering of spins in a dipolar spin ice on the kagome lattice

Spin ice, a peculiar thermal state of a frustrated ferromagnet on the pyrochlore lattice, has a finite entropy density and excitations carrying magnetic charge. By combining analytical arguments and Monte Carlo simulations, we show that spin ice on the two-dimensional kagome lattice orders in two stages. The intermediate phase has ordered magnetic charges and is separated from the paramagnetic phase by an Ising transition. The transition to the low-temperature phase is of the three-state Potts or Kosterlitz-Thouless type, depending on the presence of defects in charge order.

preprint2010arXiv

Halfvortices in flat nanomagnets

We discuss a new type of topological defect in XY systems where the O(2) symmetry is broken in the presence of a boundary. Of particular interest is the appearance of such defects in nanomagnets with a planar geometry. They are manifested as kinks of magnetization along the edge and can be viewed as halfvortices with winding numbers \pm 1/2. We argue that halfvortices play a role equally important to that of ordinary vortices in the statics and dynamics of flat nanomagnets. Domain walls found in experiments and numerical simulations are composite objects containing two or more of these elementary defects. We also discuss a closely related system: the two-dimensional smectic liquid crystal films with planar boundary condition.

preprint2010arXiv

Noncoplanar magnetic ordering driven by itinerant electrons on the pyrochlore lattice

Exchange interaction tends to favor collinear or coplanar magnetic orders in rotationally invariant spin systems. Indeed, such magnetic structures are usually selected by thermal or quantum fluctuations in highly frustrated magnets. Here we show that a complex noncoplanar magnetic order with a quadrupled unit cell is stabilized by itinerant electrons on the pyrochlore lattice. Specifically we consider the Kondo-lattice and Hubbard models at quarter filling. The electron Fermi 'surface' at this filling factor is topologically equivalent to three intersecting Fermi circles. Perfect nesting of the Fermi lines leads to magnetic ordering with multiple wavevectors and a definite handedness. The chiral order might persist without magnetic order in a chiral spin liquid at finite temperatures.

preprint2010arXiv

Pyrochlore antiferromagnet with antisymmetric exchange interactions: critical behavior and order from disorder

We investigate the nature of phase transitions induced by Dzyaloshinskii-Moriya (DM) interactions on a classical pyrochlore antiferromagnet. For both symmetry-allowed antisymmetric exchange interactions, the macroscopic degeneracy due to geometrical frustration is relieved and a long-range magnetic order appears in the ground state. We find an Ising-like phase transition in the case of direct DM interaction which selects a doubly degenerate all-in all-out magnetic order. In the presence of indirect DM coupling, the magnet undergoes an XY-like phase transition into a state with broken Z6 symmetry. We show that the critical behavior in both cases is modified due to the constrained spin fluctuations in the correlated Coulomb phase. We also demonstrate an interesting order from disorder phenomenon where the system switches between two distinct types of broken Z6 symmetry.

preprint2010arXiv

Quantum 120-degree model on pyrochlore lattice: orbital ordering in MnV2O4

We present an analytical model of orbital ordering in vanadium spinel MnV2O4. The model is based on recent first-principles calculation indicating a strong trigonal distortion at the vanadium sites of this compound [Phys. Rev. Lett. 102, 216405 (2009)]. At the single-ion level, the trigonal crystal field leaves a doubly degenerate atomic ground state and breaks the approximate rotational symmetry of t2g orbitals. We find that the effective interaction between the low-energy doublets is described by a quantum antiferromagnetic 120-degree model on the pyrochlore lattice. We obtain the classical ground state and show its stability against quantum fluctuations. The corresponding orbital order consisting of two inequivalent orbital chains is consistent with the experimentally observed tetragonal symmetry. A periodic modulation of electron density function along orbital chains is shown to arise from the staggering of local trigonal axes. In the presence of orbital order, single-ion spin anisotropy arising from relativistic spin-orbit interaction stabilizes the experimentally observed orthogonal magnetic structure.

preprint2010arXiv

Three-dimensional topological phases in a layered honeycomb spin-orbital model

We present an exactly solvable spin-orbital model based on the Gamma-matrix generalization of a Kitaev-type Hamiltonian. In the presence of small magnetic fields, the model exhibits a critical phase with a spectrum characterized by topologically protected Fermi points. Upon increasing the magnetic field, Fermi points carrying opposite topological charges move toward each other and annihilate at a critical field, signaling a phase transition into a gapped phase with trivial topology in three dimensions. On the other hand, by subjecting the system to a staggered magnetic field, an effective time-reversal symmetry essential to the existence of three-dimensional topological insulators is restored in the auxiliary free fermion problem. The nontrivial topology of the gapped ground state is characterized by an integer winding number and manifests itself through the appearance of gapless Majorana fermions confined to the two-dimensional surface of a finite system.

preprint2009arXiv

Large-J approach to strongly coupled spin-orbital systems

We present a novel approach to study the ground state and elementary excitations in compounds where spins and orbitals are entangled by on-site relativistic spin-orbit interaction. The appropriate degrees of freedom are localized states with an effective angular momentum J. We generalize J to arbitrary large values while maintaining the delicate spin-orbital entanglement. After projecting the inter-site exchange interaction to the manifold of effective spins, a systematic 1/J expansion of the effective Hamiltonian is realized using the Holstein-Primakoff transformation. Applications to representative compounds Sr2IrO4 and particularly vanadium spinels AV2O4 are discussed.

preprint2009arXiv

Model for frustrated spin-orbital chains: application to CaV2O4

Motivated by recent interest in quasi-one-dimensional compound CaV2O4, we study the ground states of a spin-orbital chain characterized by an Ising-like orbital Hamiltonian and frustrated interactions between S=1 spins. The on-site spin-orbit interaction and Jahn-Teller effect compete with inter-site superexchange leading to a rich phase diagram in which an antiferro-orbital phase is separated from the orbital paramagnet by a continuous Ising transition. Two distinct spin liquids depending on the underlying orbital order are found in the limit of small spin-orbit coupling. In the opposite limit, the zigzag chain behaves as a spin-2 chain with Ising anisotropy. The implications for CaV2O4 are discussed.