We investigate how the spin-1 Affleck-Kennedy-Lieb-Tasaki (AKLT) Hamiltonian can emerge from a microscopic fermionic model based on half-filled Hubbard tripods. We first show that a single tripod hosts a robust threefold-degenerate low-energy manifold corresponding to an effective $S = 1$ degree of freedom. This manifold prevails over a broad range of interactions and remains stable against moderate disorder. We then combine exact diagonalization with fourth-order quasi-degenerate perturbation theory to derive an effective bilinear-biquadratic spin model for a pair of coupled tripods and identify coupling regimes where the target ratio is approached. In particular, tuning leg-center hopping together with two symmetry-inequivalent leg-leg hoppings yields the characteristic singlet-triplet degeneracy associated with a biquadratic-to-bilinear ratio close to 1/3. Extending the analysis to three tripods, we compare nonequivalent coupling geometries and find a strategy that suppresses unwanted longer-range and multispin terms while preserving the target nearest-neighbor couplings in the weak-coupling regime. These results establish a concrete bottom-up route from Hubbard clusters to valence-bond-solid spin physics in tunable quantum-dot arrays.
We investigate how the spin-1 Affleck-Kennedy-Lieb-Tasaki (AKLT) Hamiltonian can emerge from a microscopic fermionic model based on half-filled Hubbard tripods. We first show that a single tripod hosts a robust threefold-degenerate low-energy manifold corresponding to an effective S=1 degree of freedom. This manifold prevails over a broad range of interactions and remains stable against moderate disorder. We then combine exact diagonalization with fourth-order quasi-degenerate perturbation theory to derive an effective bilinear-biquadratic spin model for a pair of coupled tripods and identify coupling regimes where the target ratio is approached. In particular, tuning leg-center hopping together with two symmetry-inequivalent leg-leg hoppings yields the characteristic singlet-triplet degeneracy associated with a biquadratic-to-bilinear ratio close to 1/3. Extending the analysis to three tripods, we compare nonequivalent coupling geometries and find a strategy that suppresses unwanted longer-range and multispin terms while preserving the target nearest-neighbor couplings in the weak-coupling regime. These results establish a concrete bottom-up route from Hubbard clusters to valence-bond-solid spin physics in tunable quantum-dot arrays.
Abstract We introduce a variational state for one-dimensional two-orbital Hubbard models that intuitively explains the recent computational discovery of pairing in these systems when hole doped. Our ansatz is an optimized linear superposition of Affleck–Kennedy–Lieb–Tasaki valence-bond states, rendering the combination a valence-bond liquid dubbed orbital resonant valence bond. We show that the undoped (one-electron/orbital) quantum state of two sites coupled into a global spin singlet is exactly written employing only spin-1/2 singlets linking orbitals at nearest-neighbor sites. Generalizing to longer chains defines our variational state visualized geometrically expressing our chain as a two-leg ladder, with one orbital per leg. As in Anderson’s resonating valence-bond state, our undoped variational state contains preformed singlet pairs that via doping become mobile, leading to superconductivity. Doped real materials with one-dimensional substructures, two near-degenerate orbitals, and intermediate Hubbard U/W strengths—W the carrier’s bandwidth—could realize spin-singlet pairing if on-site anisotropies are small. If these anisotropies are robust, spin-triplet pairing emerges.
The Affleck-Kennedy-Lieb-Tasaki (AKLT) Hamiltonian is a particular instance of a general class of model Hamiltonians defined in lattices with coordination z where each site hosts a spins S=z/2, interacting both with linear and nonlinear exchange couplings. In two dimensions, the AKLT model features a gap in the spectrum, and its ground state is a valence bond solid state; that is a universal resource for measurement based quantum computing, motivating the quest of physical systems that realize this Hamiltonian. Given a finite-size system described with a specific instance of this general class of models, we address the question of how to assess if such system is a realization of the AKLT model using inelastic tunnel spectroscopy implemented with scanning tunnel microscopy (IETS-STM). We propose two approaches. First, in the case of a dimer, we show how to leverage nonequilibrium IETS-STM to obtain the energies of all excited states, and determine thereby the magnitude of both linear and nonlinear exchange interactions. Second, we explore how IETS can probe the in-gap excitations associated with edge spins in a hexamer. In the AKLT limit, spins S=3/2 at the edge of the lattice have coordination 2, giving rise to S=1/2 dangling spins that can be probed with IETS. We propose a S=1/2 effective Hamiltonian to describe the interactions between these dangling spins in the neighborhood of the AKLT point, where their degeneracy lifted.
Low-dimensional quantum systems host a variety of exotic states, such as symmetry-protected topological ground states in spin-1 Haldane chains. Real-world realizations of such states could serve as practical quantum simulators for quantum phases if the interactions can be controlled. However, many proposed models, such as the Affleck-Kennedy-Lieb-Tasaki (AKLT) state, require unconventional forms of spin interactions beyond standard Heisenberg terms, which do not naturally emerge from microscopic (Coulomb) interactions. Here, we demonstrate a general strategy to induce a biquadratic term between two spin-1 sites and to tune its strength β by placing pairs of spin-1/2 spacers in between them. β is controlled by the ratio of the Heisenberg couplings between the spin-1 sites and the spacer spins, and between the spacer spins themselves. Increasing this ratio increases the magnitude of ∣β∣ and decreases the correlation length of edge states. Detailed atomistic calculations reveal that chains of nanographene flakes with 22 and 13 atoms, respectively, which could be realized by state-of-the-art bottom-up growth technology, yield precisely the couplings required to approach the AKLT state. These findings deliver a blueprint for engineering unconventional interactions in bottom-up synthesized quantum simulators.
We construct an exact spin-1 parent Hamiltonian for constrained valence-bond loop manifolds on the checkerboard and pyrochlore lattices. The Hamiltonian is local, SU(2)- and time-reversal-invariant, and built from positive-semidefinite projectors acting on triangular faces. Each projector removes only the maximally polarized state of a triangle, so the model is frustration-free. Its zero-energy states are generated by an AKLT-like construction in which each spin-1 moment is resolved into two virtual spin-21 degrees of freedom, singlets are formed inside every crossed plaquette or tetrahedron, and the physical spin-1 Hilbert space is recovered by projection. The resulting ground states are fully packed singlet-loop states on the corner-sharing lattice. Thus, a loop or dimer constraint, usually introduced as part of an effective Rokhsar--Kivelson description, appears here as the exact zero-energy manifold of a microscopic spin Hamiltonian. We analyze spin correlations within this manifold and show that, for a fixed loop covering, they are determined by loop connectivity. We also project symmetry-allowed perturbations into the ground-state manifold and derive the resulting low-energy pseudospin dynamics. The checkerboard and pyrochlore cases differ sharply. On the pyrochlore lattice, tetrahedral symmetry removes simple local bias terms, and the leading nontrivial next-nearest-neighbour Heisenberg perturbation gives an emergent spin-21 XY model on the diamond lattice of tetrahedron centers. These results give an exact spin-1 microscopic starting point for constrained valence-bond physics in two and three dimensions, and show how loop, dimer, and gauge-theoretic descriptions can be approached from a small-spin, SU(2)-invariant frustrated magnet.
Inspired by the Affleck-Kennedy-Lieb-Tasaki (AKLT) model, we present exact solutions for a spin-1 chain with Kitaev-like couplings. We consider an expanded Kitaev model with bilinear and biquadratic terms. At an exactly solvable point, the Hamiltonian can be reexpressed as a sum of projection operators. Unlike the AKLT model where projectors act on total spin, we project the component of spin along the bond direction. This leads to exponential ground-state degeneracy, expressed in terms of fractionalized spin-1/2 objects. Each ground state can be expressed concisely as a matrix product state. We construct a phase diagram by varying the relative strength of bilinear and biquadratic terms. The fractionalized states provide a qualitative picture for the spin-1 Kitaev model, yielding approximate forms for the ground state and low-lying excitations.
Abstract We studied a multi-orbital Hubbard model at half-filling for two and three orbitals per site on a two-site cluster via full exact diagonalization, in a wide range for the onsite repulsion U, from weak to strong coupling, and multiple ratios of the Hund coupling J H to U. The hopping matrix elements among the orbitals were also varied extensively. At intermediate and large U, we mapped the results into a Heisenberg model. For two orbitals per site, the mapping is into a S = 1 Heisenberg model where by symmetry both nearest-neighbor (S i ⋅ S j ) and ( S i ⋅ S j ) 2 are allowed, with respective couplings J 1 and J 2. For the case of three orbitals per site, the mapping is into a S = 3/2 Heisenberg model with (S i ⋅ S j ), ( S i ⋅ S j ) 2 , and ( S i ⋅ S j ) 3 terms, and respective couplings J 1, J 2, and J 3. The strength of these coupling constants in the Heisenberg models depend on the U, J H, and hopping amplitudes of the underlying Hubbard model. Our study provides a first crude estimate to establish bounds on how large the ratios J 2/J 1 and J 3/J 1 can be. We show that those ratios appear rather limited and, as a qualitative guidance, we conclude that J 2/J 1 is less than 0.4 and J 3/J 1 is less than 0.2, establishing bounds on effective models for strongly correlated Hubbard systems. Moreover, the intermediate Hubbard U regime was found to be the most promising to enhance J 2/J 1 and J 3/J 1.