Researcher profile

Elbio Dagotto

Elbio Dagotto contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
27works
0followers
5topics
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

27 published item(s)

preprint2023arXiv

Spinon continuum in the Heisenberg quantum chain compound Sr$_2$V$_3$O$_9$

Magnetic excitations in the spin chain candidate Sr$_2$V$_3$O$_9$ have been investigated by inelastic neutron scattering on a single crystal sample. A spinon continuum with a bandwidth of $\sim22$ meV is observed along the chain formed by alternating magnetic V$^{4+}$ and nonmagnetic V$^{5+}$ ions. Incipient magnetic Bragg peaks due to weak ferromagnetic interchain couplings emerge when approaching the magnetic transition at $T_N\sim 5.3$ K while the excitations remain gapless within the instrumental resolution. Comparisons to the Bethe ansatz, density matrix renormalization group (DMRG) calculations, and effective field theories confirm Sr$_2$V$_3$O$_9$ as a host of weakly coupled $S = 1/2$ chains dominated by antiferromagnetic intrachain interactions of $\sim7.1$(1) meV.

preprint2022arXiv

Electronic structure, magnetic properties and pairing tendencies of the copper-based honeycomb lattice Na$_2$Cu$_2$TeO$_6$

Spin-$1/2$ chains with alternating antiferromagnetic and ferromagnetic couplings have attracted considerable interest due to the topological character of their spin excitations. Here, using density functional theory and density matrix renormalization group methods, we have systematically studied the dimerized chain system Na$_2$Cu$_2$TeO$_6$. Near the Fermi level, the dominant states are mainly contributed by the Cu $3d_{x^2-y^2}$ orbitals highly hybridized with the O $2p$ orbitals in the nonmagnetic phase, leading to an "effective" single-orbital low-energy model. Furthermore, the bandwidth of the Cu $3d_{x^2-y^2}$ states is small ($\sim 0.8$ eV), suggesting that electronic correlations will strongly affect this system. By introducing such electronic correlations, we found this system is a Mott insulator. Moreover, by calculating the magnetic exchange interactions ($J_1$, $J_2$ and $J_3$), we explained the size and sign of the exchange interactions in Na$_2$Cu$_2$TeO$_6$, in agreement with neutron experiments. In addition, we constructed a single-orbital Hubbard model for this dimerized chain system, where the quantum fluctuations are taken into account. Both AFM and FM coupling ($\uparrow$-$\downarrow$-$\downarrow$-$\uparrow$) along the chain were found in our DMRG and Lanczos calculations, in agreement with DFT and neutron results. We also calculated the hole pairing binding energy $ΔE$ which becomes negative at Hubbard $U \sim 11$ eV, indicating incipient pairing tendencies. Finally, we also looked at various cases of hole doping that always exhibit tight pairs. Thus, we believe our results for Na$_2$Cu$_2$TeO$_6$ could provide guidance to experimentalists and theorists working on this dimerized chain system, such as short-range magnetic coupling, doping effects, and possible pairing tendencies.

preprint2022arXiv

Interfacial phase frustration stabilizes unconventional skyrmion crystals

Chiral magnetic phases with an unconventional topological twist in the magnetization are of huge interest due to their potential in spintronics applications. Here, we present a general method to induce such exotic magnetic phases using interfacial phase frustration within artificially grown superlattices. To demonstrate our method, we consider a multilayer with two different chiral magnetic phases as the competing orders at the top and bottom and show, using Monte Carlo calculations, that the interfacial phase frustration is realized at the central layer. In particular, we obtain three unconventional phases: a checkerboard skyrmion crystal, an incommensurate skyrmion stripe, and a ferrimagnetic skyrmion crystal. In these frustration-induced phases, the spin chirality driven topological Hall conductivity can be largely enhanced. This method provides a playground to realize unconventional magnetic phases in any family of materials that can be grown in superlattices.

preprint2022arXiv

Magnetic excitations, non-classicality and quantum wake spin dynamics in the Hubbard chain

Recent work has demonstrated that quantum Fisher information (QFI), a witness of multipartite entanglement, and magnetic Van Hove correlations $G(r,t)$, a probe of local real-space real-time spin dynamics, can be successfully extracted from inelastic neutron scattering on spin systems through accurate measurements of the dynamical spin structure factor $S(k,ω)$. Here we apply theoretically these ideas to the half-filled Hubbard chain with nearest-neighbor hopping, away from the strong-coupling limit. This model has nontrivial redistribution of spectral weight in $S(k,ω)$ going from the non-interacting limit ($U=0$) to strong coupling ($U\rightarrow \infty$), where it reduces to the Heisenberg quantum spin chain. We use the density matrix renormalization group (DMRG) to find $S(k,ω)$, from which QFI is then calculated. We find that QFI grows with $U$. With realistic energy resolution it becomes capable of witnessing bipartite entanglement above $U=2.5$ (in units of the hopping), where it also changes slope. This point is also proximate to slope changes of the bandwidth $W(U)$ and the half-chain von Neumann entanglement entropy. We compute $G(r,t)$ by Fourier-transforming $S(k,ω)$. The results indicate a crossover in the short-time short-distance dynamics at low $U$ characterized by ferromagnetic lightcone wavefronts, to a Heisenberg-like behavior at large $U$ featuring antiferromagnetic lightcones and spatially period-doubled antiferromagnetism. We find this crossover has largely been completed by $U=3$. Our results thus provide evidence that, in several aspects, the strong-coupling limit of the Hubbard chain is reached qualitatively already at a relatively modest interaction strength. We discuss experimental candidates for observing the $G(r,t)$ dynamics found at low $U$.

preprint2022arXiv

Magnetic ground states of honeycomb lattice Wigner crystals

Lattice Wigner crystal states stabilized by long-range Coulomb interactions have recently been realized in two-dimensional moiré materials. We employ large-scale unrestricted Hartree-Fock techniques to unveil the magnetic phase diagrams of honeycomb lattice Wigner crystals. For the three lattice filling factors with the largest charge gaps, n = 2/3, 1/2, 1/3, the magnetic phase diagrams contain multiple phases, including ones with non-collinear and non-coplanar spin arrangements. We discuss magnetization evolution with external magnetic field, which has potential as an experimental signature of exotic spin states. Our theoretical results could potentially be validated in moiré materials formed from group VI transition metal dichalcogenide twisted homobilayers.

preprint2022arXiv

Semi-Dirac and Weyl Fermions in Transition Metal Oxides

We show that a class of compounds with $I$4/$mcm$ crystalline symmetry hosts three-dimensional semi-Dirac fermions. Unlike the known two-dimensional semi-Dirac points, the degeneracy of these three-dimensional semi-Dirac points is not lifted by spin-orbit coupling due to the protection by a nonsymmorphic symmetry -- screw rotation in the $a-b$ plane and a translation along the $c$ axis. This crystalline symmetry is found in tetragonal perovskite oxides, realizable in thin films by epitaxial strain that results in a$^0$a$^0$c$^-$-type octahedral rotation. Interestingly, with broken time-reversal symmetry, two pairs of Weyl points emerge from the semi-Dirac points within the Brillouin zone, and an additional lattice distortion leads to enhanced intrinsic anomalous Hall effect. The ability to tune the Berry phase by epitaxial strain can be useful in novel oxide-based electronic devices.

preprint2022arXiv

Strongly anisotropic electronic and magnetic structures in oxide dichlorides RuOCl$_2$ and OsOCl$_2$

Here, using density functional theory and density matrix renormalization group methods, we investigate the electronic and magnetic properties of RuOCl$_2$ and OsOCl$_2$ with $d^4$ electronic configurations. Different from a previous study using VOI$_2$ with $d^1$ configuration, these systems with $4d^4$ or $5d^4$ do not exhibit a ferroelectric instability along the $a$-axis. Due to the fully-occupied $d_{xy}$ orbital in RuOCl$_2$ and OsOCl$_2$, the Peierls instability distortion disappears along the $b$-axis, leading to an undistorted I${\rm mmm}$ phase (No. 71). Furthermore, we observe strongly anisotropic electronic and magnetic structures along the $a$-axis. The large crystal-field splitting energy (between $d_{xz/yz}$ and $d_{xy}$ orbitals) and large hopping between nearest-neighbor Ru and Os atoms suppresses the spin-orbital effect in $M$OCl$_2$ ($M$ = Ru or Os) with electronic density $n = 4$, resulting in a spin-1 system instead of a $J = 0$ singlet ground state. Moreover, we find staggered antiferromagnetic order with $π$ wavevector along the $M$-O chain direction ($a$-axis) while the magnetic coupling along the $b$-axis is weak. Based on Wannier functions from first-principles calculations, we calculated the relevant hopping amplitudes and crystal-field splitting energies of the $t_{2g}$ orbitals for the Os atoms to construct a multi-orbital Hubbard model for the $M$-O chains. Staggered AFM with $\uparrow$-$\downarrow$-$\uparrow$-$\downarrow$ spin structure dominates in our DMRG calculations, in agreement with DFT calculations.

preprint2022arXiv

Theoretical study of the crystal and electronic properties of $α$-RuI$_3$

The material $α$-RuCl$_3$, with a two-dimensional Ru-honeycomb sublattice, has attracted considerable attention because it may be a realization of the Kitaev quantum spin liquid (QSL). Recently, a new honeycomb material, $α$-RuI$_3$, was prepared under moderate high-pressure and it is stable under ambient conditions. However, different from $α$-RuCl$_3$, $α$-RuI$_3$ was reported to be a paramagnetic metal without long-range magnetic order down to $0.35$ K. Here, the structural and electronic properties of the quasi-two-dimensional $α$-RuI$_3$ are theoretically studied. First, based on first-principles density functional theory (DFT) calculations, the ABC stacking honeycomb-layer $R\overline{3}$ (No. 148) structure is found to be the most likely stacking order for $α$-RuI$_3$ along the $c$-axis. Furthermore, both $R\overline{3}$ and $P\overline{3}1c$ are dynamically stable because no imaginary frequency modes were obtained in the phononic dispersion spectrum. Moreover, the different physical behavior of $α$-RuI$_3$ compared to $α$-RuCl$_3$ can be understood naturally. The strong hybridization between Ru $4d$ and I $5p$ orbitals decreases the effective atomic Hubbard repulsion $U$, leading the electrons of RuI$_3$ to be less localized than in RuCl$_3$. As a consequence, the effective repulsion $U$ is reduced from Cl to I, leading to the metallic nature of $α$-RuI$_3$. Based on the DFT+$U$ ($U_{\rm eff} = 2$ eV), plus spin-orbital coupling (SOC), we obtained a spin-orbit Mott insulating behavior for $α$-RuCl$_3$ and, by the same procedure, a metallic behavior for $α$-RuI$_3$, in good agreement with experimental results. Furthermore, when introducing a large (unrealistic) $U_{\rm eff} = 6$ eV, the spin-orbit Mott gap opens in $α$-RuI$_3$ as well, supporting the physical picture we are proposing.

preprint2022arXiv

Topological flat bands in a kagomé lattice multiorbital system

Flat bands and dispersive Dirac bands are known to coexist in the electronic bands in a two-dimensional kagome lattice. Including the relativistic spin-orbit coupling, such systems often exhibit nontrivial band topology, allowing for gapless edge modes between flat bands at several locations in the band structure, and dispersive bands or at the Dirac band crossing. Here, we theoretically demonstrate that a multiorbital system on a kagome lattice is a versatile platform to explore the interplay between nontrivial band topology and electronic interaction. Specifically, here we report that the multiorbital kagome model with the atomic spin-orbit coupling naturally supports topological bands characterized by nonzero Chern numbers $\cal C$, including a flat band with $|{\cal C}| =1$. When such a flat band is $1/3$ filled, the non-local repulsive interactions induce a fractional Chern insulating state. We also discuss the possible realization of our findings in real kagome materials.

preprint2021arXiv

Coupled Hubbard ladders at weak coupling: Pairing and spin excitations

The Hubbard model provides a simple framework in which one can study how certain aspects of the electronic structure of strongly interacting systems can be tuned to optimize the superconducting pairing correlations and how these changes affect the mechanisms giving rise to them. Here we use a weak-coupling random phase approximation to study a two-dimensional Hubbard model with a unidirectional modulation of the hopping amplitudes as the system evolves from the uniform square lattice to an array of weakly coupled two-leg ladders. We find that the pairing correlations retain their dominant $d_{x^2-y^2}$-wave like structure and that they are significantly enhanced for a slightly modulated lattice. This enhancement is traced backed to an increase in the strength of the spin-fluctuation pairing interacting due to favorable Fermi surface nesting in the modulated system. We then use a random-phase approximation BCS framework to examine the evolution of the neutron resonance in the superconducting state. We find that it changes only weakly for moderate modulations, but breaks up into two distinct resonances at incommensurate wave-vectors in the limit of weakly coupled ladders.

preprint2021arXiv

Damped Dirac magnon in a metallic kagome antiferromagnet FeSn

The kagome lattice is a fertile platform to explore topological excitations with both Fermi-Dirac and Bose-Einstein statistics. While relativistic Dirac Fermions and flat-bands have been discovered in the electronic structure of kagome metals, the spin excitations have received less attention. Here we report inelastic neutron scattering studies of the prototypical kagome magnetic metal FeSn. The spectra display well-defined spin waves extending up to 120 meV. Above this energy, the spin waves become progressively broadened, reflecting interactions with the Stoner continuum. Using linear spin wave theory, we determine an effective spin Hamiltonian that reproduces the measured dispersion. This analysis indicates that the Dirac magnon at the K-point remarkably occurs on the brink of a region where well-defined spin waves become unobservable. Our results emphasize the influential role of itinerant carriers on the topological spin excitations of metallic kagome magnets.

preprint2021arXiv

Magnetization dynamics fingerprints of an excitonic condensate $t_{2g}^{4}$ magnet

The competition between spin-orbit coupling $λ$ and electron-electron interaction $U$ leads to a plethora of novel states of matter, extensively studied in the context of $t_{2g}^4$ and $t_{2g}^5$ materials, such as ruthenates and iridates. Excitonic magnets -- the antiferromagnetic state of bounded electron-hole pairs -- is a prominent example of phenomena driven by those competing energy scales. Interestingly, recent theoretical studies predicted that excitonic magnets can be found in the ground-state of spin-orbit-coupled $t_{2g}^4$ Hubbard models. Here, we present a detailed computational study of the magnetic excitations in that excitonic magnet, employing one-dimensional chains (via density matrix renormalization group) and small two-dimensional clusters (via Lanczos). Specifically, first we show that the low-energy spectrum is dominated by a dispersive (acoustic) magnonic mode, with extra features arising from the $λ=0$ state in the phase diagram. Second, and more importantly, we found a novel magnetic excitation forming a high-energy optical mode with the highest intensity at wavevector $q\to 0$. In the excitonic condensation regime at large $U$, we also have found a novel high-energy $π$-mode composed solely of orbital excitations. These unique fingerprints of the $t_{2g}^4$ excitonic magnet are important in the analysis of neutron and RIXS experiments.

preprint2021arXiv

Orbital ordering in the layered perovskite material CsVF$_4$

In strongly correlated electronic systems, several novel physical properties are induced by the orbital degree of freedom. In particular, orbital degeneracy near the Fermi level leads to spontaneous symmetry breaking, such as the nematic state in FeSe and the orbital ordering in several perovskite systems. Here, the novel layered perovskite material CsVF$_4$, with a $3d^2$ electronic configuration, was systematically studied using density functional theory and a multiorbital Hubbard model within the Hatree-Fock approximation. Our results show that CsVF$_4$ should be magnetic, with a G-type antiferromagnetic arrangement in the $ab$ plane and weak antiferromagnetic exchange along the $c$-axis, in agreement with experimental results. Driven by the Jahn-Teller distortion in the VF$_6$ octahedra that shorten the $c$-axis, the system displays an interesting electron occupancy $d_{xy}^1(d_{xz}d_{yz})^1$ corresponding to the lower nondegenerate $d_{xy}$ orbital being half-filled and the other two degenerate $d_{yz}$ and $d_{xz}$ orbitals sharing one electron per site. We show that this degeneracy is broken and a novel $d_{yz}$/$d_{xz}$ staggered orbital pattern is here predicted by both the first-principles and Hubbard model calculations. This orbital ordering is driven by the electronic instability associated with degeneracy removal to lower the energy.

preprint2021arXiv

Range of biquadratic and triquadratic Heisenberg effective couplings deduced from multiorbital Hubbard models

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 $(\mathbf{S}_{i}\cdot\mathbf{S}_{j})$ and $(\mathbf{S}_{i}\cdot\mathbf{S}_{j})^{2}$ are allowed, with respective couplings $J_1$ and $J_2$. For the case of three orbitals per site, the mappping is into a $S=3/2$ Heisenberg model with $(\mathbf{S}_{i}\cdot\mathbf{S}_{j})$, $(\mathbf{S}_{i}\cdot\mathbf{S}_{j})^{2}$, and $(\mathbf{S}_{i}\cdot\mathbf{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 allows to establish bounds on how large the ratios $J_2/J_1$ and $J_3/J_1$ can be. We show that those ratios are severely limited and, as a crude 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.

preprint2020arXiv

BCS-BEC crossover in a $(t_{2g})^4$ Excitonic Magnet

The condensation of spin-orbit-induced excitons in $(t_{2g})^4$ electronic systems is attracting considerable attention. In the large Hubbard U limit, antiferromagnetism was proposed to emerge from the Bose-Einstein Condensation (BEC) of triplons ($J_{\textrm{eff}} = 1$). In this publication, we show that even for the weak and intermediate U regimes, the spin-orbit exciton condensation is possible leading also to staggered magnetic order. The canonical electron-hole excitations (excitons) transform into local triplon excitations at large U , and this BEC strong coupling regime is smoothly connected to the intermediate U excitonic insulator region. We solved the degenerate three-orbital Hubbard model with spin-orbit coupling ($λ$) in one-dimensional geometry using the Density Matrix Renormalization Group, while in two-dimensional square clusters we use the Hartree-Fock approximation (HFA). Employing these techniques, we provide the full $λ$ vs U phase diagrams for both one- and two- dimensional lattices. Our main result is that at the intermediate Hubbard U region of our focus, increasing $λ$ at fixed U the system transitions from an incommensurate spin-density-wave metal to a Bardeen-Cooper-Schrieffer (BCS) excitonic insulator, with coherence length r coh of O(a) and O(10a) in 1d and 2d, respectively, with a the lattice spacing. Further increasing $λ$, the system eventually crosses over to the BEC limit (with r coh << a).

preprint2020arXiv

Charge-density-wave melting in the one-dimensional Holstein model

We study the Holstein model of spinless fermions, which at half-filling exhibits a quantum phase transition from a metallic Tomonaga-Luttinger liquid phase to an insulating charge-density-wave (CDW) phase at a critical electron-phonon coupling strength. In our work, we focus on the real-time evolution starting from two different types of initial states that are CDW ordered: (i) ideal CDW states with and without additional phonons in the system and (ii) correlated ground states in the CDW phase. We identify the mechanism for CDW melting in the ensuing real-time dynamics and show that it strongly depends on the type of initial state. We focus on the far-from-equilibrium regime and emphasize the role of electron-phonon coupling rather than dominant electronic correlations, thus complementing a previous study of photo-induced CDW melting [H. Hashimoto and S. Ishihara, Phys. Rev. B 96, 035154 (2017)]. The numerical simulations are performed by means of matrix-product-state based methods with a local basis optimization (LBO). Within these techniques, one rotates the local (bosonic) Hilbert spaces adaptively into an optimized basis that can then be truncated while still maintaining a high precision. In this work, we extend the time-evolving block decimation (TEBD) algorithm with LBO, previously applied to single-polaron dynamics, to a half-filled system. We demonstrate that in some parameter regimes, a conventional TEBD method without LBO would fail. Furthermore, we introduce and use a ground-state density-matrix renormalization group method for electron-phonon systems using local basis optimization. In our examples, we account for up to $M_{\rm ph} = 40$ bare phonons per site by working with $O(10)$ optimal phonon modes.

preprint2020arXiv

Density Matrix Renormalization Group Study of Nematicity in Two Dimensions: Application to a Spin-$1$ Bilinear-Biquadratic Model on the Square Lattice

Nematic order is an exotic property observed in several strongly correlated systems, such as the iron-based superconductors. Using large-scale density matrix renormalization group (DMRG) techniques, we study at zero-temperature the nematic spin liquid that competes with spin dipolar and quadrupolar orders. We use these nematic orders to characterize different quantum phases and quantum phase transitions. More specifically, we study a spin-$1$ bilinear-biquadratic Heisenberg model on the square lattice with couplings beyond nearest neighbors. We focus on parameter regions around the highly symmetric $SU(3)$ point where the bilinear and biquadratic interactions are equal. With growing further-neighbor biquadratic interactions, we identify different spin dipolar and quadrupolar orders. We find that the DMRG results on cylindrical geometries correctly detect nematicity in different quantum states and accurately characterize the quantum phase transitions among them. Therefore, spin-driven nematicity -- here defined as the spontaneous breaking of the lattice invariance under a 90$^o$ rotation -- is an order parameter which can be studied directly in DMRG calculations in two dimensions in different quantum states.

preprint2020arXiv

Emergence of Superconductivity in Doped Multiorbital Hubbard Chains

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&#39;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&#39;s bandwidth -- could realize spin-singlet pairing if on-site anisotropies are small. If these anisotropies are robust, spin-triplet pairing emerges.

preprint2020arXiv

First-principles study of the low-temperature charge density wave phase in the quasi-one-dimensional Weyl chiral compound (TaSe$_4$)$_2$I

Using {\it ab initio} density functional theory, we study the lattice phase transition of quasi-one-dimensional (TaSe$_4$)$_2$I. In the undistorted state, the strongly anisotropic semimetal band structure presents two non-equivalent Weyl points. In previous efforts, two possible Ta-tetramerization patterns were proposed to be associated with the low-temperature structure. Our phonon calculations indicate that the orthorhombic $F222$ CDW-I phase is the most likely ground state for this quasi-one-dimensional system. In addition, the monoclinic $C2$ CDW-II phase may also be stable according to the phonon dispersion spectrum. Since these two phases have very similar energies in our DFT calculations, both these Ta-tetramerization distortions likely compete or coexist at low temperatures. The semimetal to insulator transition is induced by a Fermi-surface-driven instability that supports the Peierls scenario, which affects the Weyl physics developed above $T_{\rm CDW}$. Furthermore, the spin-orbit coupling generates Rashba-like band splittings in the insulating CDW phases.

preprint2020arXiv

Flat Bands and Ferrimagnetic Order in Electronically Correlated Dice-Lattice Ribbons

We study ribbons of the dice two-dimensional lattice (that we call ``dice ladders&#39;&#39;) known to have nontrivial topological properties, such as Chern numbers 2 [Wang and Y. Ran, Phys. Rev. B {\bf 84}, 241103 (2011)]. Our main results are two folded: (1) Analyzing the tight-binding model in the presence of Rashba spin-orbit coupling and an external magnetic field, we observed that dice ladders qualitatively display properties similar to their two-dimensional counterpart all the way to the limit of only two legs in the short direction. This includes flat bands near the Fermi level, edge currents and edge charge localization near zero energy when open boundary conditions are used, two chiral edge modes, and a nonzero Hall conductance. (2) We studied the effect of Hubbard correlation $U$ in the two-leg dice ladder using Lanczos and density matrix renormalization group techniques. We show that increasing $U$ the flat bands split without the need of introducing external fields. Moreover, robust ferrimagnetic order develops. Overall, our work establishes dice ladders as a promising playground to study the combined effect of topology and correlation effects, one of the frontiers in Quantum Materials.

preprint2020arXiv

Fractionalized Excitations Revealed by Entanglement Entropy

Fractionalized excitations develop in many unusual many-body states such as quantum spin liquids, disordered phases that cannot be described using any local order parameter. Because these exotic excitations correspond to emergent degrees of freedom, how to probe them and establish their existence is a long-standing challenge. We present a general procedure to reveal the fractionalized excitations using real-space entanglement entropy in critical spin liquids that are particularly relevant to experiments. Moreover, we show how to use the entanglement entropy to construct the corresponding spinon Fermi surface. Our work defines a new pathway to establish and characterize exotic excitations in novel quantum phases of matter.

preprint2020arXiv

Magnetic Switching in Weyl Semimetal-Superconductor Heterostructures

We present a new switching mechanism that utilizes the proximity coupling between the surface spin texture of a Weyl semimetal and a superconductor, in a Weyl semimetal/superconductor/Weyl semimetal trilayer heterostructure. We show that the superconductivity in the middle layer can be fully suppressed by the surface spin texture of the Weyl semimetals in the presence of an external magnetic field, but it can be recovered again by only changing the field direction. The restoration of the middle-layer superconductivity indicates a sharp transition to a low-resistance state. This sharp switching effect, realizable using a Weyl semimetal because of its strong spin-momentum locking and surface spin polarization, is a promising avenue for novel superconducting spin-valve applications.

preprint2020arXiv

Origin of the Magnetic and Orbital ordering in $α$-Sr$_2$CrO$_4$

Motivated by recent experimental progress in transition metal oxides with the K$_2$NiF$_4$ structure, we investigate the magnetic and orbital ordering in $α$-Sr$_2$CrO$_4$. Using first principles calculations, first we derive a three-orbital Hubbard model, which reproduces the {\it ab initio} band structure near the Fermi level. The unique reverse splitting of $t_{2g}$ orbitals in $α$-Sr$_2$CrO$_4$, with the $3d^2$ electronic configuration for the Cr$^{4+}$ oxidation state, opens up the possibility of orbital ordering in this material. Using real-space Hartree-Fock for multi-orbital systems, we constructed the ground-state phase diagram for the two-dimensional compound $α$-Sr$_2$CrO$_4$. We found stable ferromagnetic, antiferromagnetic, antiferro-orbital, and staggered orbital stripe ordering in robust regions of the phase diagram. Furthermore, using the density matrix renormalization group method for two-leg ladders with the realistic hopping parameters of $α$-Sr$_2$CrO$_4$, we explore magnetic and orbital ordering for experimentally relevant interaction parameters. Again, we find a clear signature of antiferromagnetic spin ordering along with antiferro-orbital ordering at moderate to large Hubbard interaction strength. We also explore the orbital-resolved density of states with Lanczos, predicting insulating behavior for the compound $α$-Sr$_2$CrO$_4$, in agreement with experiments. Finally, an intuitive understanding of the results is provided based on a hierarchy between orbitals, with $d_{xy}$ driving the spin order, while electronic repulsion and the effective one dimensionality of the movement within the $d_{xz}$ and $d_{yz}$ orbitals driving the orbital order.

preprint2020arXiv

Planar Topological Hall Effect from Conical Spin Spirals

Planar Topological Hall effect (PTHE) is known to arise from in-plane skyrmion tubes with an external magnetic field applied in the plane of the charge current. Here, we propose that the PTHE can robustly appear also from an unexpected source, the conical spin spiral, realizable in a variety of materials. We show that for both in-plane skyrmion tubes and conical spin spirals, the PTHE is two-fold symmetric with respect to the magnetic-field angle, while it is four-fold symmetric for 90$^{\circ}$-rotated domains of conical spin spirals. We predict that the symmetry and magnetic-field response of the PTHE can distinguish between the two scenarios, unambiguously probing the nature of the antisymmetric spin-exchange interaction and the resulting magnetic texture.

preprint2020arXiv

Prediction of exotic magnetic states in the alkali metal quasi-one-dimensional iron selenide compound Na$_2$FeSe$_2$

The magnetic and electronic phase diagram of a model for the quasi-one-dimensional alkali metal iron selenide compound Na$_2$FeSe$_2$ is presented. The novelty of this material is that the valence of iron is Fe$^{2+}$ contrary to most other iron-chain compounds with valence Fe$^{3+}$. Using first-principles techniques, we developed a three-orbital tight-binding model that reproduces the {\it ab initio} band structure near the Fermi level. Including Hubbard and Hund couplings and studying the model via the density matrix renormalization group and Lanczos methods, we constructed the ground state phase diagram. A robust region where the block state $\uparrow \uparrow \downarrow \downarrow \uparrow \uparrow \downarrow \downarrow$ is stabilized was unveiled. The analog state in iron ladders, employing 2$\times$2 ferromagnetic blocks, is by now well-established, but in chains a block magnetic order has not been observed yet in real materials. The phase diagram also contains a large region of canonical staggered spin order $\uparrow \downarrow \uparrow \downarrow \uparrow \downarrow \uparrow$ at very large Hubbard repulsion. At the block to staggered transition region, a novel phase is stabilized with a mixture of both states: an inhomogeneous orbital-selective charge density wave with the exotic spin configuration $\uparrow \uparrow \downarrow \uparrow \downarrow \downarrow \uparrow \downarrow$. Our predictions for Na$_2$FeSe$_2$ may guide crystal growers and neutron scattering experimentalists towards the realization of block states in one-dimensional iron-selenide chain materials.

preprint2020arXiv

Quantum Transitions of Nematic Phases in a Spin-$1$ Bilinear-Biquadratic Model and Their Implications for FeSe

Since its discovery, iron-based superconductivity has been known to develop near an antiferromagnetic order, but this paradigm fails in the iron chalcogenide FeSe, whose single-layer version holds the record for the highest superconducting transition temperature in the iron-based superconductors. The striking puzzle that FeSe displays nematic order (spontaneously broken lattice rotational symmetry) while being non-magnetic, has led to several competing proposals for its origin in terms of either the $3d$-electron&#39;s orbital degrees of freedom or spin physics in the form of frustrated magnetism. Here we argue that the phase diagram of FeSe under pressure could be qualitatively described by a quantum spin model with highly frustrated interactions. We implement both the site-factorized wave-function analysis and the large-scale density matrix renormalization group (DMRG) in cylinders to study the spin-$1$ bilinear-biquadratic model on the square lattice, and identify quantum transitions from the well-known $(π,0)$ antiferromagnetic state to an exotic $(π,0)$ antiferroquadrupolar order, either directly or through a $(π/2,π)$ antiferromagnetic state. These many phases, while distinct, are all nematic. We also discuss our theoretical ground-state phase diagram for the understanding of the experimental low-temperature phase diagram obtained by the NMR [P. S. Wang {\it et al.}, Phys. Rev. Lett. 117, 237001 (2016)] and X-ray scattering [K. Kothapalli {\it et al.}, Nature Communications 7, 12728 (2016)] measurements in pressurized FeSe. Our results suggest that superconductivity in a wide range of iron-based materials has a common origin in the antiferromagnetic correlations of strongly correlated electrons.

preprint2019arXiv

Iron telluride ladder compounds: Predicting the structural and magnetic properties of BaFe$_2$Te$_3$

Since the discovery of pressure-induced superconductivity in the two-leg ladder system BaFe$_2X_3$ ($X$=S, Se), with the 3$d$ iron electronic density $n = 6$, the quasi-one-dimensional iron-based ladders have attracted considerable attention. Here, we use Density Functional Theory (DFT) to predict that the novel $n = 6$ iron ladder BaFe$_2$Te$_3$ could be stable with a similar crystal structure as BaFe$_2$Se$_3$. Our results also indicate that BaFe$_2$Te$_3$ will display the complex 2$\times$2 Block-type magnetic order. Due to the magnetic striction effects of this Block order, BaFe$_2$Te$_3$ should be a magnetic noncollinear ferrielectric system with a net polarization $0.31$ $μ$C/cm$^2$. Compared with the S- or Se-based iron ladders, the electrons of the Te-based ladders are more localized, implying that the degree of electronic correlation is enhanced for the Te case which may induce additional interesting properties. The physical and structural similarity with BaFe$_2$Se$_3$ also suggests that BaFe$_2$Te$_3$ could become superconducting under high pressure.