Researcher profile

Hae-Young Kee

Hae-Young Kee contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
21works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

21 published item(s)

preprint2026arXiv

Spin-orbit coupling controlled two-dimensional magnetism in chromium trihalides

CrX$_3$ (X = Cl, Br, I) have the same crystal structure and Hamiltonian but different ligand spin-orbit coupling (SOC) constant $λ_X$, providing excellent material platform exploring for exotic two-dimensional (2D) spin orders. Their microscopic mechanism underlying 2D spin physics remain unestablished, along with experimental corroboration of Kitaev exchange interaction, central to realizing topological quantum spin liquids. Finding direct evidence for Kitaev interaction and determining its value has been an essential but formidable challenge in Kitaev physics. Here we report the direct Kitaev interaction signature in magnetic anisotropy measured by ferromagnetic resonance (FMR) spectroscopy. We present measured values of Heisenberg $J$, Kitaev $K$, and off-diagonal symmetric $Γ$ exchange interactions in CrX$_3$ determined using FMR and exact diagonalization. $K$ and $Γ$ exhibit dominant dependencies on $λ_X$, indicating its central role in 2D magnetism. Our study provides a foundation for designing 2D magnetic materials exhibiting novel behaviors by tuning intrinsic material parameters such as SOC.

preprint2023arXiv

Reconciling $π$ phase shift in Josephson junction experiments with even-parity superconductivity in Sr$_2$RuO$_4$

The superconducting state of Sr$_2$RuO$_4$ was once thought to be a leading candidate for $p$-wave superconductivity. A constant Knight shift below the transition temperature provided evidence for spin-triplet pairing, and a $π$ phase shift observed in Josephson junction tunneling experiments suggested odd-parity pairing, both of which are described by $p$-wave states. However, with recent experiments observing a significant decrease in the Knight shift below the transition temperature, signifying a spin-singlet state, the odd-parity results are left to be reconciled. In this work, we show that an even-parity pseudospin-singlet state originating from interorbital pairing via spin-orbit coupling can explain what has been assumed to be evidence for an odd-parity state. In the presence of small mirror symmetry breaking, interorbital pairing is uniquely capable of displaying odd-parity characteristics required to explain these experimental results. Further, we discuss how these experiments may be used to differentiate the proposed pairing states of Sr$_2$RuO$_4$.

preprint2022arXiv

Competing multipolar orders in a face-centered cubic lattice: Application to the osmium double perovskites

In 5$d^2$ Mott insulators with strong spin-orbit coupling, the lowest pseudospin states form a non-Kramers doublet, which carries quadrupolar and octupolar moments. A family of double perovskites where magnetic ions form a face-centered cubic (fcc) lattice was suggested to unveil an octupolar order offering a rare example in $d$-orbital systems. The proposed order requires a ferromagnetic (FM) octupolar interaction, since the antiferromagnetic (AFM) Ising model is highly frustrated on the fcc lattice. A microscopic model was recently derived for various lattices: for an edge-sharing octahedra geometry, AFM Ising octupolar and bond-dependent quadrupolar interactions were found when only dominant inter- and intraorbital hopping integrals were taken into account. Here we investigate all possible intra- and interorbital exchange processes and report that interference of two intraorbital exchanges generates a FM octupolar interaction. Applying the strong-coupling expansion results together with tight-binding parameters obtained by density functional theory, we estimate the exchange interactions for the osmium double perovskites, Ba$_2$BOsO$_6$ (B=Mg,Cd, and Ca). Using classical Monte Carlo simulations, we find that these systems are close to the phase boundary between AFM type-I quadrupole and FM octupole orders. We also find that exchange processes beyond second-order perturbation theory including virtual processes via pseudospin-triplet states may stabilize an octupolar order.

preprint2022arXiv

Counter-rotating spiral, zigzag, and 120$^\circ$ orders from coupled-chain analysis of Kitaev-Gamma-Heisenberg model, and relations to honeycomb iridates

We study the nearest neighboring spin-1/2 Kitaev-Heisenberg-Gamma ($KJΓ$) model on the honeycomb lattice in the parameter region of ferromagnetic (FM) Kitaev and antiferromagnetic (AFM) Heisenberg couplings relevant for honeycomb iridates, using a coupled-chain analysis. Starting from the gapless Luttinger liquid phase of a decoupled $KJΓ$ chain, the inter-chain interactions in the two-dimensional model is treated within a self-consistent mean field approach based on the Luttinger liquid theory. In the FM Gamma region, our analysis recovers the reported 120$^\circ$ magnetic order, previously obtained by classical analysis and exact diagonalization method. On the other hand, new physics is revealed in the AFM Gamma region, where three magnetic orders are found, including 120$^\circ$, commensurate counter-rotating spiral, and zigzag orders. Interestingly, the two first order phase transition lines separating these three magnetic orders merge at a single point at $K = -2 Γ$ and $J=0$, which is predicted to be a quantum critical point. The current theory captures the experimentally observed counter-rotating spiral order in $α$-Li$_22$IrO$_3$ and the zigzag order in Na$_2$IrO$_3$, thereby indicating that the spin-1/2 $KJΓ$ model may serve as a minimal model for honeycomb iridates. Limitations of the mean field theory presented in this work and the $J \rightarrow 0$ regime are also discussed.

preprint2022arXiv

Evolution of Interorbital Superconductor to Intraorbital Spin-Density Wave in Layered Ruthenates

The ruthenate family of layered perovskites has been a topic of intense interest with much work dedicated to understanding the superconducting state of the single layer, Sr$_2$RuO$_4$. Another longstanding puzzle is the lack of superconductivity in its sister compound, Sr$_3$Ru$_2$O$_7$, which constrains the possible superconducting mechanisms of Sr$_2$RuO$_4$. Here we address a microscopic mechanism that unifies superconducting and spin-density wave order in this family of materials. Beginning from a model of Sr$_2$RuO$_4$ featuring intraband pseudospin-singlet superconductivity originating from interorbital spin-triplet pairing via Hund's and spin-orbit couplings, the addition of bilayer coupling and staggered rotations of the oxygen octahedra are investigated. We find that the bilayer coupling alone enhances superconductivity, while staggered rotations destroy interorbital superconductivity, leaving a paramagnetic metal. The magnetic field then shifts van Hove singularities near the Fermi level, allowing intraorbital SDW order to form in the bilayer. Our theory predicts that bilayer Sr$_3$Ru$_2$O$_7$ without staggered rotations exhibits interorbital superconductivity with a possibly higher transition temperature.

preprint2022arXiv

Interplay of magnetic field and trigonal distortion in the honeycomb $Γ$ model: Occurrence of a spin-flop phase

In candidate Kitaev materials, the off-diagonal $Γ$ and $Γ'$ interactions are identified to come from the spin-orbit coupling and trigonal distortion, respectively. They have generated intense research efforts because of their intimate relation to the field-induced magnetically disordered state reported in $α$-RuCl$_3$. Theoretically, while a plethora of field-induced phases has been proposed in the honeycomb lattice, a stable intermediate phase that can survive in a wide parameter region regardless of the underlying phases is still lacking. Here we focus on the interplay of an out-of-plane magnetic field and a symmetry-allowed $Γ'$ term due to trigonal distortion in the dominant antiferromagnetic $Γ$ region. By using multifaceted approaches ranging from classical Monte Carlo and semiclassical spin-wave theory to density-matrix renormalization group, we identify an intriguing spin-flop phase in the presence of magnetic field and antiferromagnetic $Γ'$ interaction, before it eventually enters into a fully polarized state. As the $Γ'$ interaction approaches the size of $Γ$ one, the $Γ$-$Γ'$ model maps to the easy-axis XXZ antiferromagnet, where the spin-flop phase can be understood as a superfluid phase in the extended Bose-Hubbard model. Our finding thus demonstrates an exciting path from the honeycomb $Γ$ model towards a $U(1)$-symmetric XXZ antiferromagnet in a magnetic field.

preprint2022arXiv

Spontaneous Chiral-Spin Ordering in Spin-Orbit Coupled Honeycomb Magnets

Frustrated magnets with highly degenerate ground states are at the heart of hunting exotic states of matter. Recent studies in spin-orbit coupled honeycomb magnets have generated immense interest in bond-dependent interactions, appreciating a symmetric off-diagonal $Γ$ interaction which exhibits a macroscopic degeneracy in the classical limit. Here, we study a generic spin model and discover a novel chiral-spin ordering with spontaneously broken time-reversal symmetry near the dominant $Γ$ region. The chiral-spin phase is demonstrated to possess a staggered chirality relation in different sublattices, and it exhibits gapless excitations as revealed by the vanishing energy gap and the finite central charge on cylinders. Although there is a vestige of a tiny peak in the corner of the second Brillouin zone, the magnetic order is likely to vanish as the system size increases. Finally, we also attempt to gain insight into the possible topological signature of the chiral-spin phase by calculating the dynamic structure factor and the modular $\mathcal{S}$ matrix.

preprint2022arXiv

Strategy to Extract Kitaev Interaction using Symmetry in Honeycomb Mott Insulators

The Kitaev spin liquid, a ground state of the bond-dependent Kitaev model in a honeycomb lattice has been a centre of attraction, since a microscopic theory to realize such an interaction in solid-state materials was discovered. A challenge in real materials though is the presence of the Heisenberg and another bond-dependent Gamma interactions detrimental to the Kitaev spin liquid, and there have been many debates on their relative strengths. Here we offer a strategy to extract the Kitaev interaction out of a full microscopic model by utilizing the symmetries of the Hamiltonian. Two tilted magnetic field directions related by a two-fold rotational symmetry generate distinct spin excitations originating from a specific combination of the Kitaev and Gamma interactions. Together with the in- and out-of-plane magnetic anisotropy, one can determine the Kitaev and Gamma interactions separately. Dynamic spin structure factors are presented to motivate future experiments. The proposed setups will advance the search for Kitaev materials.

preprint2022arXiv

Theoretical analysis of single-ion anisotropy in $d^3$ Mott insulators

An effective spin model for Mott insulators is determined by the symmetries involved among magnetic sites, electron fillings, and their interactions. Such a spin Hamiltonian offers insight to mechanisms of magnetic orders and magnetic anisotropy beyond the Heisenberg model. For a spin moment S bigger than 1/2, single-ion anisotropy is in principle allowed. However, for $d^3$ Mott insulators with large cubic crystal field splitting, the single-ion anisotropy is absent within the LS coupling, despite S = 3/2 local moment. On the other hand, preferred magnetic moment directions in $d^3$ materials have been reported, which calls for a further theoretical investigation. Here we derive the single-ion anisotropy interaction using the strong-coupling perturbation theory. The cubic crystal field splitting including $e_g$ orbitals, trigonal distortions, Hund's coupling, and spin-orbit coupling beyond the LS scheme are taken into account. For compressed distortion, the spin-orbit coupling at magnetic sites can favor either the easy-axis or the easy-plane while that of anions leads to easy-axis anisotropy. We apply the theory on $\rm{CrX}_3$ with X = Cl and I, and show the dependence of the single-ion anisotropy on the strength of the spin-orbit couplings of both magnetic and anion sites. Significance of the single-ion anisotropy in ideal two-dimensional magnets is also discussed.

preprint2021arXiv

Characterizing spin-one Kitaev quantum spin liquids

Material realizations of the bond-dependent Kitaev interactions with $S$=1/2 local moments have vitalized the research in quantum spin liquids. Recently, it has been proposed that higher-spin analogues of the Kitaev interactions may also occur in a number of materials with strong spin-orbit coupling. In contrast to the celebrated $S$=1/2 Kitaev model on the honeycomb lattice, the higher-spin Kitaev models are not exactly solvable. Hence, the existence of quantum spin liquids in these systems remains an outstanding question. In this work, we use the density matrix renormalization group (DMRG) methods to numerically investigate the $S$=1 Kitaev model with both ferromagnetic (FM) and antiferromagnetic (AFM) interactions. Using results on a cylindrical geometry with various circumferences, we conclude that the ground state of the $S$=1 Kitaev model is a quantum spin liquid with a $\mathbb{Z}_2$ gauge structure. We also put a bound on the excitation gap, which turns out to be quite small. The magnetic field responses for the FM and AFM models are similar to those of the $S$=1/2 counterparts. In particular, in the AFM $S$=1 model, a gapless quantum liquid state emerges in an intermediate window of magnetic field strength, before the system enters a trivial polarized state.

preprint2021arXiv

Classical spin order near antiferromagnetic Kitaev point in the spin-1/2 Kitaev-Gamma chain

A minimal Kitaev-Gamma model has been recently investigated to understand various Kitaev systems. In the one-dimensional Kitaev-Gamma chain, an emergent SU(2)$_1$ phase and a rank-1 spin ordered phase with $O_h\rightarrow D_4$ symmetry breaking were identified using non-Abelian bosonization and numerical techniques. However, puzzles near the antiferromagnetic Kitaev region with finite Gamma interaction remained unresolved. Here we focus on this parameter region and find that there are two new phases, namely, a rank-1 ordered phase with an $O_h\rightarrow D_3$ symmetry breaking, and a peculiar Kitaev phase. Remarkably, the $O_h\rightarrow D_3$ symmetry breaking corresponds to the classical magnetic order, but appears in a region very close to the antiferromagnetic Kitaev point where the quantum fluctuations are presumably very strong. In addition, a two-step symmetry breaking $O_h\rightarrow D_{3d}\rightarrow D_3$ is numerically observed as the length scale is increased: At short and intermediate length scales, the system behaves as having a rank-2 spin nematic order with $O_h\rightarrow D_{3d}$ symmetry breaking; and at long distances, time reversal symmetry is further broken leading to the $O_h\rightarrow D_3$ symmetry breaking. Finally, there is no numerical signature of spin orderings nor Luttinger liquid behaviors in the Kitaev phase whose nature is worth further studies.

preprint2021arXiv

Insights into the Anisotropic Spin-$S$ Kitaev Chain

Recently there has been a renewed interest in properties of the higher-spin Kitaev models, especially their low-dimensional analogues with additional interactions. These quasi-1D systems exhibit rich phase diagrams with symmetry-protected topological phases, Luttinger liquids, hidden order, and higher-rank magnetism. However, the nature of the pure spin-$S$ Kitaev chains are not yet fully understood. Earlier works found a unique ground state with short-ranged correlations for $S = 1$, and an intriguing double-peak structure in the heat capacity associated with an entropy plateau. To understand the low-energy excitations and thermodynamics for general $S$ we study the anisotropic spin-$S$ Kitaev chain. Starting from the dimerized limit we derive an effective low-energy Hamiltonian at finite anisotropy. For half-integer spins we find a trivial effective model, reflecting a non-local symmetry protecting the degeneracy, while for integer $S$ we find interactions among the flux degrees of freedom that select a unique ground state. The effective model for integer spins is used to predict the low-energy excitations and thermodynamics, and we make a comparison with the semiclassical limit through linear spin wave theory. Finally, we speculate on the nature of the isotropic limit.

preprint2021arXiv

Optical magnons with dominant bond-directional exchange interactions in a honeycomb lattice iridate $α$-Li$_{2}$IrO$_{3}$

We have used resonant inelastic x-ray scattering to reveal optical magnons in a honeycomb lattice iridate $α$-Li$_{2}$IrO$_{3}$. The spectrum in the energy region 20-25 meV exhibits momentum dependence, of which energy is highest at the location of the magnetic Bragg peak, ($\textit{h}, \textit{k}$) = ($\pm$0.32, 0), and lowered toward (0, 0) and ($\pm$1, 0). We compare our data with a linear spin-wave theory based on a generic nearest-neighbor spin model. We find that a dominant bond-directional Kitaev interaction of order 20 meV is required to explain the energy scale observed in our study. The observed excitations are understood as stemming from optical magnon modes whose intensity is modulated by a structure factor, resulting in the apparent momentum dependence. We also observed diffuse magnetic scattering arising from the short-range magnetic correlation well above $\textit{T}_{N}$. In contrast to Na$_{2}$IrO$_{3}$, this diffuse scattering lacks the $C_3$ rotational symmetry of the honeycomb lattice, suggesting that the bond anisotropy is far from negligible in $α$-Li$_{2}$IrO$_{3}$.

preprint2021arXiv

Possible Realization of Optical Quadratic and Dirac Points in Woodpile Photonic Crystals

The simulation of fermionic relativistic physics, e.g., Dirac and Weyl physics, has led to the discovery of many unprecedented phenomena in photonics, of which the optical-frequency realization is, however, still challenging. Here, surprisingly, we discover that the woodpile photonic crystals commonly used for optical frequency applications host exotic fermion-like relativistic degeneracies: a Dirac nodal line and a fourfold quadratic point, as protected by the nonsymmorphic crystalline symmetry. Deforming the woodpile photonic crystal leads to the emergence of type-II Dirac points from the fourfold quadratic point. Such type-II Dirac points can be detected by its anomalous refraction property which is manifested as a giant birefringence in a slab setup. Our findings provide a promising route towards 3D optical Dirac physics in all-dielectric photonic crystals.

preprint2021arXiv

Shadowed triplet pairings in Hund's metals with spin-orbit coupling

Hund's coupling in multiorbital systems allows for the possibility of even-parity orbital-antisymmetric spin-triplet pairing, which can be stabilized by spin-orbit coupling (SOC). While this pairing expressed in the orbital basis is uniform and spin-triplet, it appears in the band basis as a pseudospin-singlet, with the momentum dependence determined by the SOC and the underlying triplet character remaining in the form of interband pairing active away from the Fermi energy. Here, we examine the role of momentum-dependent SOC in generating nontrivial pairing symmetries, as well as the hidden triplet nature associated with this interorbital pairing, which we dub a "shadowed triplet". Applying this concept to Sr$_{2}$RuO$_{4}$, we first derive several forms of SOC with $d$-wave form factors from a microscopic model, and subsequently we show that for a range of SOC parameters, a pairing state with $s+id_{xy}$ symmetry can be stabilized. Such a pairing state is distinct from pure spin-singlet and -triplet pairings due to its unique character of pseudospin-energy locking. We discuss experimental probes to differentiate the shadowed triplet pairing from conventional pseudospin-triplet and -singlet pairings.

preprint2020arXiv

Distinct reduction of Knight shift in superconducting state of Sr$_2$RuO$_4$ under uniaxial strain

Shortly after the discovery of superconductivity in Sr$_2$RuO$_4$, spin-triplet pairing was proposed and further corroborated by a constant Knight shift ($K$) across the transition temperature (T$_c$). However, a recent experiment observed a drop in $K$ at T$_c$ which becomes larger under uniaxial strain, ruling out several spin-triplet scenarios. Here we show that inter-orbital spin-triplet pairing features a d-vector that rotates when uniaxial strain is applied, leading to a larger drop in the spin polarization perpendicular to the strain direction, distinct from spin-singlet pairing. We propose that anisotropic spin polarization under strain will ultimately differentiate triplet vs.~singlet pairing.

preprint2020arXiv

Field-Driven Gapless Spin Liquid in the Spin-1 Kitaev Honeycomb Model

Recent proposals for spin-1 Kitaev materials, such as honeycomb Ni oxides with heavy elements of Bi and Sb, have shown that these compounds naturally give rise to antiferromagnetic (AFM) Kitaev couplings. Conceptual interest in such AFM Kitaev systems has been sparked by the observation of a transition to a gapless $U(1)$ spin liquid at intermediate field strengths in the AFM spin-1/2 Kitaev model. However, all hitherto known spin-1/2 Kitaev materials exhibit ferromagnetic bond-directional exchanges. Here we discuss the physics of the spin-1 Kitaev model in a magnetic field and show, by extensive numerical analysis, that for AFM couplings it exhibits an extended gapless quantum spin liquid at intermediate field strengths. The close analogy to its spin-1/2 counterpart suggests that this gapless spin liquid is a $U(1)$ spin liquid with a neutral Fermi surface, that gives rise to enhanced thermal transport signatures.

preprint2020arXiv

Honeycomb-Lattice Mott insulator on Tantalum Disulphide

Effects of electron many-body interactions amplify in an electronic system with a narrow bandwidth opening a way to exotic physics. A narrow band in a two-dimensional (2D) honeycomb lattice is particularly intriguing as combined with Dirac bands and topological properties but the material realization of a strongly interacting honeycomb lattice described by the Kane-Mele-Hubbard model has not been identified. Here we report a novel approach to realize a 2D honeycomb-lattice narrow-band system with strongly interacting 5$d$ electrons. We engineer a well-known triangular lattice 2D Mott insulator 1T-TaS$_2$ into a honeycomb lattice utilizing an adsorbate superstructure. Potassium (K) adatoms at an optimum coverage deplete one-third of the unpaired $d$ electrons and the remaining electrons form a honeycomb lattice with a very small hopping. Ab initio calculations show extremely narrow Z$_2$ topological bands mimicking the Kane-Mele model. Electron spectroscopy detects an order of magnitude bigger charge gap confirming the substantial electron correlation as confirmed by dynamical mean field theory. It could be the first artificial Mott insulator with a finite spin Chern number.

preprint2020arXiv

Phase diagram of the spin-1/2 Kitaev-Gamma chain and emergent SU(2) symmetry

We study the phase diagram of a one-dimensional version of the Kitaev spin-1/2 model with an extra ``$Γ$-term", using analytical, density matrix renormalization group and exact diagonalization methods. Two intriguing phases are found. In the gapless phase, although the exact symmetry group of the system is discrete, the low energy theory is described by an emergent SU(2)$_1$ Wess-Zumino-Witten (WZW) model. On the other hand, the spin-spin correlation functions exhibit SU(2) breaking prefactors, even though the exponents and the logarithmic corrections are consistent with the SU(2)$_1$ predictions. A modified nonabelian bosonization formula is proposed to capture such exotic emergent ``partial" SU(2) symmetry. In the ordered phase, there is numerical evidence for an $O_h\rightarrow D_4$ spontaneous symmetry breaking.

preprint2020arXiv

Testing Ising Topological Order in $α$-RuCl$_3$ Under In-Plane Magnetic Fields

Material realization of the non-Abelian Kitaev spin liquid phase - an example of Ising topological order (ITO) - has been the subject of intense research in recent years. The $4d$ honeycomb Mott insulator $α$-RuCl$_3$ has emerged as a leading candidate, as it enters a field-induced magnetically disordered state where a half-integer quantized thermal Hall conductivity $κ_{xy}$ was reported. Further, a recent report of a sign change in the quantized $κ_{xy}$ across a certain crystallographic direction is strong evidence for a topological phase transition between two ITOs with opposite Chern numbers. Although this is a fascinating result, independent verification remains elusive, and one may ask if there is a thermodynamic quantity sensitive to the phase transition. Here we propose that the magnetotropic coefficient $k$ under in-plane magnetic fields would serve such a purpose. We report a singular feature in $k$ that indicates a topological phase transition across the $\hat{b}$-axis where ITO is prohibited by a $C_2$ symmetry. If the transition in $α$-RuCl$_3$ is indeed a direct transition between ITOs, then this feature in $k$ should be observable.

preprint2020arXiv

The Heart of Entanglement: Chiral, Nematic, and Incommensurate Phases in the Kitaev-Gamma Ladder in a Field

The bond-dependent Kitaev model on the honeycomb lattice with anyonic excitations has recently attracted considerable attention. However, in solid state materials other spin interactions are present, and among such additional interactions, the off-diagonal symmetric Gamma interaction, another type of bond-dependent term, has been particularly challenging to fully understand. A minimal Kitaev-Gamma (KG) model has been investigated by various numerical techniques under a magnetic field, but definite conclusions about field-induced spin liquids remain elusive and one reason may lie in the limited sizes of the two-dimensional geometry it is possible to access numerically. We therefore focus on the KG model defined on a two-leg ladder which is much more amenable to a complete study, and determine the entire phase diagram in the presence of a magnetic field along [111]-direction. Due to the competition between the interactions and the field, an extremely rich phase diagram emerges with fifteen distinct phases. Focusing near the antiferromagnetic Kitaev region, we identify nine different phases solely within this region: several incommensurate magnetically ordered phases, spin nematic, and two chiral phases with enhanced entanglement. Of particular interest is a highly entangled phase with staggered chirality with zero-net flux occurring at intermediate field, which along with its companion phases outline a heart-shaped region of high entanglement, the heart of entanglement. We compare our results for the ladder with a C3 symmetric cluster of the two-dimensional honeycomb lattice, and offer insights into possible spin liquids in the two-dimensional limit.