Source author record

Elbio Dagotto

Elbio Dagotto 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

94works
7topics
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

94 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'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.

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'') 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'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.

preprint2016arXiv

Bicollinear Antiferromagnetic Order, Monoclinic Distortion, and Reversed Resistivity Anisotropy in FeTe as a Result of Spin-Lattice Coupling

The bicollinear antiferromagnetic order experimentally observed in FeTe is shown to be stabilized by the coupling $\tilde g_{12}$ between monoclinic lattice distortions and the spin-nematic order parameter with $B_{\rm 2g}$ symmetry, within a three-orbital spin-fermion model studied with Monte Carlo techniques. A finite but small value of $\tilde g_{12}$ is required, with a concomitant lattice distortion compatible with experiments, and a tetragonal-monoclinic transition strongly first order. Remarkably, the bicollinear state found here displays a planar resistivity with the "reversed" puzzling anisotropy discovered in transport experiments.Orthorhombic distortions are also incorporated and phase diagrams interpolating between pnictides and chalcogenides are presented. We conclude that the spin-lattice coupling discussed here is sufficient to explain the challenging properties of FeTe.

preprint2016arXiv

Inversion of ferrimagnetic magnetization by ferroelectric switching via a novel magnetoelectric coupling

Although several multiferroic materials/heterostructures have been extensively studied, finding strong magnetoelectric couplings for the electric field control of the magnetization remains challenging. Here, a novel interfacial magnetoelectric coupling based on three components (ferroelectric dipole, magnetic moment, and antiferromagnetic order) is analytically formulated. As an extension of carrier-mediated magnetoelectricity, the new coupling is shown to induce an electric-magnetic hysteresis loop. Realizations employing BiFeO$_3$ bilayers grown along the [$111$] axis are proposed. Without involving magnetic phase transitions, the magnetization orientation can be switched by the carrier modulation driven by the field effect, as confirmed using first-principles calculations.

preprint2016arXiv

Magnetic properties and pairing tendencies of the iron-based superconducting ladder BaFe$_2$S$_3$: combined ab initio and density matrix renormalization group study

The recent discovery of superconductivity under high pressure in the two-leg ladder compound BaFe$_2$S$_3$ [H. Takahashi et al., Nature Materials 14, 1008 (2015)] opens a broad avenue of research, because it represents the first report of pairing tendencies in a quasi one-dimensional iron-based high critical temperature superconductor. Similarly as in the case of the cuprates, ladders and chains can be far more accurately studied using many-body techniques and model Hamiltonians than their layered counterparts, particularly if several orbitals are active. In this publication, we derive a two-orbital Hubbard model from first principles that describes individual ladders of BaFe$_2$S$_3$. The model is studied with the density matrix renormalization group. These first reported results are exciting for two reasons: (i) at half-filling, ferromagnetic order emerges as the dominant magnetic pattern along the rungs of the ladder, and antiferromagnetic order along the legs, in excellent agreement with neutron experiments; (ii) pairs form in the strong coupling regime, as found by studying the binding energy of two holes doped on the half-filled system. In addition, Orbital Selective Mott Phase characteristics develop with doping, with only one Wannier orbital receiving the hole carriers while the other remains half-filled. These results suggest that the analysis of models for iron-based two-leg ladders could clarify the origin of pairing tendencies and other exotic properties of iron-based high critical temperature superconductors.

preprint2016arXiv

On-site Atractive Multiorbital Hamiltonian for $d$-Wave Superconductors

We introduce a two-orbital Hamiltonian on a square lattice that contains on-site attractive interactions involving the two $e_g$ orbitals. Via a canonical mean-field procedure similar to the one applied to the well-known negative-$U$ Hubbard model, it is shown that the new model develops $d$-wave ($B_{1g}$) superconductivity with nodes along the diagonal directions of the square Brillouin zone. This result is also supported by exact diagonalization of the model in a small cluster. The expectation is that this relatively simple attractive model could be used to address the properties of multiorbital $d$-wave superconductors in the same manner that the negative-$U$ Hubbard model is widely applied to the study of the properties of $s$-wave single-orbital superconductors. In particular, we show that by splitting the $e_g$ orbitals and working at three-quarters filling, such that the $x^2-y^2$ orbital dominates at the Fermi level but the $3z^2-r^2$ orbital contribution is nonzero, the $d$-wave pairing state found here phenomenologically reproduces several properties of the superconducting state of the high $T_c$ cuprates.

preprint2016arXiv

Orbital Selective Directional Conductor in the Two Orbital Hubbard Model

Employing a recently developed many-body technique that allows for the incorporation of thermal effects, the rich phase diagram of a two dimensional two orbital (degenerate $d_{xz}$ and $d_{yz}$) Hubbard model is presented varying temperature and the repulsion $U$. Our main result is the finding at intermediate $U$ of a novel antiferromagnetic orbital selective state where an effective dimensional reduction renders one direction insulating and the other metallic. Possible realizations of this state are discussed. In addition we also study nematicity above the Néel temperature. After a careful finite-size scaling analysis the nematicity temperature window appears to survive in the bulk limit although it is very narrow.

preprint2016arXiv

Orbital-selective Mott phases of a one-dimensional three-orbital Hubbard model studied using computational techniques

A recently introduced one-dimensional three-orbital Hubbard model displays orbital-selective Mott phases with exotic spin arrangements such as spin block states [J. Rincón {\em et al.}, Phys. Rev. Lett. \textbf{112}, 106405 (2014)]. In this publication we show that the Constrained-Path Quantum Monte Carlo (CPQMC) technique can accurately reproduce the phase diagram of this multiorbital one-dimensional model, paving the way to future CPQMC studies in systems with more challenging geometries, such as ladders and planes. The success of this approach relies on using the Hartree-Fock technique to prepare the trial states needed in CPQMC. We also study a simplified version of the model where the pair-hopping term is neglected and the Hund coupling is restricted to its Ising component. The corresponding phase diagrams are shown to be only mildly affected by the absence of these technically difficult-to-implement terms. This is confirmed by additional Density Matrix Renormalization Group and Determinant Quantum Monte Carlo calculations carried out for the same simplified model, with the latter displaying only mild Fermion sign problems. We conclude that these methods are able to capture quantitatively the rich physics of the several orbital-selective Mott phases (OSMP) displayed by this model, thus enabling computational studies of the OSMP regime in higher dimensions, beyond static or dynamic mean field approximations.

preprint2015arXiv

Isotropic Quenched Disorder Triggers a Robust Nematic State in Electron-Doped Pnictides

The phase diagram of electron-doped pnictides is studied varying the temperature, electronic density, and isotropic quenched disorder strength by means of computational techniques applied to a three-orbital ($xz$, $yz$, $xy$) spin-fermion model with lattice degrees of freedom. In experiments, chemical doping introduces disorder but in theoretical studies the relationship between electronic doping and the randomly located dopants, with their associated quenched disorder, is difficult to address. In this publication, the use of computational techniques allows us to study independently the effects of electronic doping, regulated by a global chemical potential, and impurity disorder at randomly selected sites. Surprisingly, our Monte Carlo simulations reveal that the fast reduction with doping of the Néel $T_N$ and the structural $T_S$ transition temperatures, and the concomitant stabilization of a robust nematic state, is primarily controlled by the magnetic dilution associated with the in-plane isotropic disorder introduced by Fe substitution. In the doping range studied, changes in the Fermi Surface produced by electron doping affect only slightly both critical temperatures.

preprint2015arXiv

Parallelized Traveling Cluster Approximation to Study Numerically Spin-Fermion Models on Large Lattices

Lattice spin-fermion models are important to study correlated systems where quantum dynamics allows for a separation between slow and fast degrees of freedom. The fast degrees of freedom are treated quantum mechanically while the slow variables, generically refereed to as the "spins", are treated classically. At present, exact diagonalization coupled with classical Monte Carlo (ED+MC) is extensively used to solve numerically a general class of lattice spin-fermion problems. In this common setup, the classical variables (spins) are treated via the standard MC method while the fermion problem is solved by exact diagonalization. The "Traveling Cluster Approximation" (TCA) is a real space variant of the ED+MC method that allows to solve spin-fermion problems on lattice sizes with up to $10^3$ sites. In this publication, we present a novel reorganization of the TCA algorithm in a manner that can be efficiently parallelized. This allows us to solve generic spin-fermion models easily on $10^4$ lattice sites and with some effort on $10^5$ lattice sites, representing the record lattice sizes studied for this family of models.

preprint2015arXiv

Strain Doping: Reversible Single-Axis Control of a Complex Oxide Lattice via Helium Implantation

We report on the use of helium ion implantation to independently control the out-of-plane lattice constant in epitaxial La0.7Sr0.3MnO3 thin films without changing the in-plane lattice constants. The process is reversible by a vacuum anneal. Resistance and magnetization measurements show that even a small increase in the out-of-plane lattice constant of less than 1% can shift the metal-insulator transition and Curie temperatures by more than 100 °C. Unlike conventional epitaxy-based strain tuning methods which are constrained not only by the Poisson effect but by the limited set of available substrates, the present study shows that strain can be independently and continuously controlled along a single axis. This permits novel control over orbital populations through Jahn-Teller effects, as shown by Monte Carlo simulations on a double-exchange model. The ability to reversibly control a single lattice parameter substantially broadens the phase space for experimental exploration of predictive models and leads to new possibilities for control over materials' functional properties.

preprint2014arXiv

Ab initio quantum Monte Carlo calculations of spin superexchange in cuprates: the benchmarking case of Ca$_2$CuO$_3$

In view of the continuous theoretical efforts aimed at an accurate microscopic description of the strongly correlated transition metal oxides and related materials, we show that with continuum quantum Monte Carlo (QMC) calculations it is possible to obtain the value of the spin superexchange coupling constant of a copper oxide in a quantitatively excellent agreement with experiment. The variational nature of the QMC total energy allows us to identify the best trial wave function out of the available pool of wave functions, which makes the approach essentially free from adjustable parameters and thus truly ab initio. The present results on magnetic interactions suggest that QMC is capable of accurately describing ground state properties of strongly correlated materials.

preprint2014arXiv

BaFe$_2$Se$_3$: a high $T_{\rm C}$ magnetic multiferroic with large ferrielectric polarization

The iron-selenides are important because of their superconducting properties. Here, an unexpected phenomenon is predicted to occur in an iron-selenide compound with a quasi-one-dimensional ladder geometry: BaFe$_2$Se$_3$ should be a magnetic ferrielectric system, driven by its magnetic block order via exchange striction. A robust performance (high $T_{\rm C}$ and large polarization) is expected. Different from most multiferroics, BaFe$_2$Se$_3$ is ferrielectric, with a polarization that mostly cancels between ladders. However, its strong magnetostriction still produces a net polarization that is large ($\sim$$0.1$ $μ$C/cm$^2$) as compared with most magnetic multiferroics. Its fully ferroelectric state, with energy only slightly higher than the ferrielectric, has a giant improper polarization $\sim$$2-3$ $μ$C/cm$^2$.

preprint2014arXiv

Competition between covalent bonding and charge transfer at complex-oxide interfaces

Here we study the electronic properties of cuprate/manganite interfaces. By means of atomic resolution electron microscopy and spectroscopy, we produce a subnanometer scale map of the transition metal oxidation state profile across the interface between the high $T_c$ superconductor YBa$_2$Cu$_3$O$_{7-δ}$ and the colossal magnetoresistance compound (La,Ca)MnO$_3$. A net transfer of electrons from manganite to cuprate with a peculiar non-monotonic charge profile is observed. Model calculations rationalize the profile in terms of the competition between standard charge transfer tendencies (due to band mismatch), strong chemical bonding effects across the interface, and Cu substitution into the Mn lattice, with different characteristic length scales.

preprint2014arXiv

Diverging Nematic Susceptibility, Physical Meaning of T* scale, and Pseudogap in the Spin Fermion Model for Pnictides

Using Monte Carlo simulations with a tunable uniaxial strain, for the first time the nematic susceptibility of a model for the pnictides is calculated. The results are in good agreement with the experiments by J-H. Chu {\it et al.}, Science {\bf 337}, 710 (2012). Via a Ginzburg-Landau analysis, our study suggests a nematicity in pnictides primarily originating on magnetism, but with the lattice/orbital boosting up critical temperatures and separating the structural $T_S$ and Néel $T_N$ transitions. At $T>T_S$, Curie-Weiss behavior is observed with the characteristic temperature $T^*$ unveiled by Chu {\it et al.} being the $T_N$ of the purely electronic system. In this temperature regime, short-range magnetic order with wavevectors $(π,0)-(0,π)$ induce local nematic fluctuations and a density-of-states pseudogap, compatible with several experiments.

preprint2014arXiv

Double magnetic resonance and spin anisotropy in Fe-based superconductors due to static and fluctuating antiferromagnetic orders

Motivated by recent neutron scattering experiments in Fe-based superconductors, we study how the magnetic resonance in the superconducting state is affected by the simultaneous presence of either static or fluctuating magnetic orders using the random phase approximation. We find that for the underdoped materials with coexisting superconducting and antiferromagnetic orders, spin rotational symmetry is explicitly broken at the ordering momentum $Q_1 = (π,0)$. Only the longitudinal susceptibility exhibits the resonance mode, whereas a spin-wave Goldstone mode develops in the transverse component. Meanwhile, at the frustrated momentum $Q_2 = (0,π)$, the susceptibility becomes isotropic in spin space and the magnetic resonance exists for both components. Furthermore, the resonance energies at $Q_1$ and $Q_2$ have distinct scales, which provides a natural explanation for the recently observed double resonance peaks. In addition, we show that near optimal doping the existence of strong magnetic fluctuations, which are modeled here via a Gaussian mode, can still induce the spin anisotropy in the magnetic susceptibility.

preprint2014arXiv

Exotic magnetic order in the orbital-selective Mott regime of multiorbital systems

The orbital-selective Mott phase (OSMP) of multiorbital Hubbard models has been extensively analyzed before using static and dynamical mean-field approximations. In parallel, the properties of Block states (antiferromagnetically coupled ferromagnetic spin clusters) in Fe-based superconductors have also been much discussed. The present effort uses numerically exact techniques in one-dimensional systems to report the observation of Block states within the OSMP regime, connecting two seemingly independent areas of research, and providing analogies with the physics of Double-Exchange models.

preprint2014arXiv

Magnetic states of the five-orbital Hubbard model for one-dimensional iron-based superconductors

The magnetic phase diagrams of models for quasi one-dimensional compounds belonging to the iron-based superconductors family are presented. The five-orbital Hubbard model and the real-space Hartree-Fock approximation are employed, supplemented by density functional theory to obtain the hopping amplitudes. Phase diagrams are constructed varying the Hubbard $U$ and Hund $J$ couplings and at zero temperature. The study is carried out at electronic density (electrons per iron) $n = 5.0$, which is of relevance for the already known material TlFeSe$_2$, and also at $n = 6.0$, where representative compounds still need to be synthesized. At $n = 5.0$ there is a clear dominance of staggered spin order along the chain direction. At $n = 6.0$ and the realistic Hund coupling $J/U = 0.25$, the phase diagram is far richer including a variety of ``block'' states involving ferromagnetic clusters that are antiferromagnetically coupled, in qualitative agreement with recent Density Matrix Renormalization Group calculations for the three-orbital Hubbard model in a different context. These block states arise from the competition between ferromagnetic order (induced by double exchange, and prevailing at large $J/U$) and antiferromagnetic order (dominating at small $J/U$). The density of states and orbital compositions of the many phases are also provided.

preprint2014arXiv

Orbital occupancy and charge doping in iron-based superconductors

Iron-based superconductors (FBS) comprise several families of compounds having the same atomic building blocks for superconductivity, but large discrepancies among their physical properties. A longstanding goal in the field has been to decipher the key underlying factors controlling TC and the various doping mechanisms. In FBS materials this is complicated immensely by the different crystal and magnetic structures exhibited by the different families. In this paper, using aberration-corrected scanning transmission electron microscopy (STEM) coupled with electron energy loss spectroscopy (EELS), we observe a universal behavior in the hole concentration and magnetic moment across different families. All the parent materials have the same total number of electrons in the Fe 3d bands; however, the local Fe magnetic moment varies due to different orbital occupancy. Although the common understanding has been that both long-range and local magnetic moments decrease with doping, we find that, near the onset of superconductivity, the local magnetic moment increases and shows a dome-like maximum near optimal doping, where no ordered magnetic moment is present. In addition, we address a longstanding debate concerning how Co substitutions induces superconductivity in the 122 arsenide family, showing that the 3d band filling increases a function of doping. These new microscopic insights into the properties of FBS demonstrate the importance of spin fluctuations for the superconducting state, reveal changes in orbital occupancy among different families of FBS, and confirm charge doping as one of the mechanisms responsible for superconductivity in 122 arsenides.

preprint2014arXiv

Photoexcitation of electronic instabilities in one-dimensional charge-transfer systems

We investigate the real-time dynamics of photoexcited electronic instabilities in a charge-transfer system model, using the time-dependent density matrix renormalization group method. The model of choice was the quarter-filled one-dimensional extended Peierls-Hubbard Hamiltonian interacting with classical few-cycle electromagnetic radiation. The results show that only one electronic instability drives the main features of the photogenerated time-dependent behavior. Indeed, the photoresponse of the system shows a large enhancement of the $4k_F$ (bond and charge) instability whereas the $2k_F$ state remains largely unaffected. This conclusion holds regardless of the nature of the optical excitations and whether the system is perturbed resonantly or not. Our results suggest potential applications of charge-transfer systems with slow phononic dynamics as optoelectronic switching devices.

preprint2014arXiv

Quantum phase transition between orbital-selective Mott states in Hund's metals

We report a quantum phase transition between orbital-selective Mott states, with different localized orbitals, in a Hund's metals model. Using the density matrix renormalization group, the phase diagram is constructed varying the electronic density and Hubbard $U$, at robust Hund's coupling. We demonstrate that this transition is preempted by charge fluctuations and the emergence of free spinless fermions, as opposed to the magnetically-driven Mott transition. The Luttinger correlation exponent is shown to have a universal value in the strong-coupling phase, whereas it is interaction dependent at intermediate couplings. At weak coupling we find a second transition from a normal metal to the intermediate-coupling phase.

preprint2014arXiv

Testing the Monte Carlo - Mean Field approximation in the one-band Hubbard model

The canonical one-band Hubbard model is studied using a computational method that mixes the Monte Carlo procedure with the mean field approximation. This technique allows us to incorporate thermal fluctuations and the development of short-range magnetic order above ordering temperatures, contrary to the crude finite-temperature Hartree-Fock approximation, which incorrectly predicts a Néel temperature $T_N$ that grows linearly with the Hubbard $U/t$. The effective model studied here contains quantum and classical degrees of freedom. It thus belongs to the "spin fermion" model family widely employed in other contexts. Using exact diagonalization, supplemented by the traveling cluster approximation, for the fermionic sector, and classical Monte Carlo for the classical fields, the Hubbard $U/t$ vs. temperature $T/t$ phase diagram is studied employing large three and two dimensional clusters. We demonstrate that the method is capable of capturing the formation of local moments in the normal state without long-range order, the non-monotonicity of $T_N$ with increasing $U/t$, the development of gaps and pseudogaps in the density of states, and the two-peak structure in the specific heat. Extensive comparisons with determinant quantum Monte Carlo results suggest that the present approach is qualitatively, and often quantitatively, accurate, particularly at intermediate and high temperatures. Finally, we study the Hubbard model including plaquette diagonal hopping (i.e. the $t-t^\prime$ Hubbard model) in two dimensions and show that our approach allows us to study low temperature properties where determinant quantum Monte Carlo fails due to the fermion sign problem. Future applications of this method include multi-orbital Hubbard models such as those needed for iron-based superconductors.

preprint2013arXiv

$B_{1g}$-like pairing states in two-leg ladder iron superconductors

Motivated by the recent report of superconductivity in Fe-based ladder materials, we study the pairing state of a multi-orbital t-J model defined on two-leg ladders using the standard mean-field theory. We find that the superconducting order parameters change sign between the $d_{xz}$ and $d_{yz}$ orbitals in most of the phase diagram. By analogy with the two-dimensional Fe planes, we conclude that the leading pairing channel of this state belongs to the $B_{1g}$ symmetry class, which is distinct from the common $s_{\pm}$ gap with the $A_{1g}$ symmetry. By smoothly interpolating from planes into ladders, we show that a first-order transition occurs between these two competing phases when the dimension of the system is reduced.

preprint2013arXiv

Electrophoretic-like gating used to control metal-insulator transitions in electronically phase separated manganite wires

Electronically phase separated manganite wires are found to exhibit controllable metal-insulator transitions under local electric fields. The switching characteristics are shown to be fully reversible, polarity independent, and highly resistant to thermal breakdown caused by repeated cycling. It is further demonstrated that multiple discrete resistive states can be accessed in a single wire. The results conform to a phenomenological model in which the inherent nanoscale insulating and metallic domains are rearranged through electrophoretic-like processes to open and close percolation channels.

preprint2013arXiv

First Principles Study of the Magnetic Properties of LaOMnAs

Recent experiments reported giant magnetoresistance at room temperature in LaOMnAs. Here a density functional theory calculation is performed to investigate magnetic properties of LaOMnAs. The ground state is found to be the G-type antiferromagnetic order within the $ab$ plane but coupled ferromagnetically between planes, in agreement with recent neutron investigations. The electronic band structures suggest an insulating state which is driven by the particular G-type magnetic order, while a metallic state accompanies the ferromagnetic order. This relation between magnetism and conductance may be helpful to qualitatively understand the giant magnetoresistance effects.

preprint2013arXiv

Full Control of Magnetism in Manganite Bilayer by Ferroelectric Polarization

An oxide heterostructure made of manganite bilayers and ferroelectric perovskites is predicted to lead to the full control of magnetism when switching the ferroelectric polarizations. By using asymmetric polar interfaces in the superlattices, more electrons occupy the Mn layer at the $n$-type interface side than at the $p$-type side. This charge disproportionation can be enhanced or suppressed by the ferroelectric polarization. Quantum model and density functional theory calculations reach the same conclusion: a ferromagnetic-ferrimagnetic phase transition with maximal change $>90%$ of the total magnetization can be achieved by switching the polarization's direction. This function is robust and provides full control of the magnetization's magnitude, not only its direction, via electrical methods.

preprint2013arXiv

Magnetic Phase Diagram of a Five-Orbital Hubbard Model in the Real-Space Hartree Fock Approximation Varying the Electronic Density

Using the real-space Hartree Fock approximation, the magnetic phase diagram of a five-orbital Hubbard model for the iron-based superconductors is studied varying the electronic density $n$ in the range from 5 to 7 electrons per transition metal atom. The Hubbard interaction $U$ is also varied, at a fixed Hund coupling $J/U=0.25$. Several qualitative trends and a variety of competing magnetic states are observed. At $n$=5, a robust G-type antiferromagnetic insulator is found, in agreement with experimental results for BaMn$_2$As$_2$. As $n$ increases away from 5, magnetic states with an increasing number of nearest-neighbors ferromagnetic links become energetically stable. This includes the well-known C-type antiferromagnetic state at $n$=6, the E-phase known to exist in FeTe, and also a variety of novel states not found yet experimentally, some of them involving blocks of ferromagnetically oriented spins. Regions of phase separation, as in Mn-oxides, have also been detected. Comparison with previous theoretical investigations indicate that these qualitative trends may be generic characteristics of phase diagrams of multiorbital Hubbard models.

preprint2013arXiv

Nematic State of the Pnictides Stabilized by the Interplay Between Spin, Orbital, and Lattice Degrees of Freedom

The nematic state of the iron-based superconductors is studied in the undoped limit of the three-orbital ($xz$, $yz$, $xy$) spin-fermion model via the introduction of lattice degrees of freedom. Monte Carlo simulations show that in order to stabilize the experimentally observed lattice distortion and nematic order, and to reproduce photoemission experiments, {\it both} the spin-lattice and orbital-lattice couplings are needed. The interplay between their respective coupling strengths regulates the separation between the structural and Néel transition temperatures. Experimental results for the temperature dependence of the resistivity anisotropy and the angle-resolved photoemission (ARPES) orbital spectral weight are reproduced by the present numerical simulations.

preprint2013arXiv

Quantum-Confinement-Induced Magnetism in LaNiO$_3$-LaMnO$_3$ Superlattices

The emergence of magnetic reconstructions at the interfaces of oxide heterostructures are often explained via subtle modifications in the electronic densities, exchange couplings, or strain. Here an additional possible route for induced magnetism is studied in the context of the (LaNiO$_3$)$_n$/(LaMnO$_3$)$_n$ superlattices using a hybrid tight-binding model. In the LaNiO$_3$ region, the induced magnetizations decouple from the intensity of charge leakage from Mn to Ni, but originate from the spin-filtered quantum confinement present in these nanostructures. In general, the induced magnetization is the largest for the (111)-stacking and the weakest for the (001)-stacking superlattices, results compatible with the exchange bias effects reported by Gibert et al. Nat. Mater. 11, 195 (2012).

preprint2013arXiv

Tunneling Electroresistance Induced by Interfacial Phase Transitions in Ultrathin Oxide Heterostructures

The ferroelectric (FE) control of electronic transport is one of the emerging technologies in oxide heterostructures. Many previous studies in FE tunnel junctions (FTJs) exploited solely the differences in the electrostatic potential across the FTJs that are induced by changes in the FE polarization direction. Here, we show that in practice the junction current ratios between the two polarization states can be further enhanced by the electrostatic modification in the correlated electron oxide electrodes, and that FTJs with nanometer thin layers can effectively produce a considerably large electroresistance ratio at room temperature. To understand these surprising results, we employed an additional control parameter, which is related to the crossing of electronic and magnetic phase boundaries of the correlated electron oxide. The FE-induced phase modulation at the heterointerface ultimately results in an enhanced electroresistance effect. Our study highlights that the strong coupling between degrees of freedom across heterointerfaces could yield versatile and novel applications in oxide electronics.

preprint2012arXiv

Anisotropy of Electrical Transport in Pnictide Superconductors Studied Using Monte Carlo Simulations of the Spin-Fermion Model

The undoped three-orbital spin fermion model for the Fe-based superconductors is studied via Monte Carlo techniques in two-dimensional clusters. At low temperatures, the magnetic and one-particle spectral properties are in good agreement with neutron and photoemission experiments. Our most important results are the resistance vs. temperature curves that display all the features experimentally observed in BaFe$_2$As$_2$ detwinned single crystals (under uniaxial stress), including a low-temperature anisotropy between the two directions followed by a peak at the magnetic ordering temperature, here induced by short-range spin order and concomitant Fermi Surface orbital order.

preprint2012arXiv

Complex State Found in the CMR Regime of Models for Colossal Magnetoresistive Manganites

In this report, we study in detail the competitor to the FM metallic state at electronic density $x=1/4$ in the CMR regime using the two-orbital double-exchange model with Jahn-Teller lattice distortions on two-dimensional clusters, employing a very careful large-scale cooling down process in the Monte Carlo simulations to avoid being trapped in metastable states. Our investigations show that this competing insulator has a very unexpected complex structure, involving diagonal stripes with alternating regions with FM and CE-like order. The level of complexity of this new state even surpasses that of the recently unveiled spin-orthogonal-stripe states and their associated high degeneracy. This new state complements the long-standing scenario of phase separation, since the alternating FM-CE pattern appears even in the present study which is carried out in the clean limit. The present and recent investigations are also in agreement with the many "glassy" characteristics of the CMR state found experimentally, due to the high degeneracy of the insulating states involved in the process. Results for the spin-structure factor of the new states are also here provided to facilitate the analysis of neutron scattering experiments for these materials.

preprint2012arXiv

Dynamics of doublon-holon pairs in Hubbard two-leg ladders

The dynamics of holon-doublon pairs is studied in Hubbard two-leg ladders using the time-dependent Density Matrix Renormalization Group method. We find that the geometry of the two-leg ladder, that is qualitatively different from a one-dimensional chain due to the presence of a spin-gap, strongly affects the propagation of a doublon-holon pair. Two distinct regimes are identified. For weak inter-leg coupling, the results are qualitatively similar to the case of the propagation previously reported in Hubbard chains, with only a renormalization of parameters. More interesting is the case of strong inter-leg coupling where substantial differences arise, particularly regarding the double occupancy and properties of the excitations such as the doublon speed. Our results suggest a connection between the presence of a spin gap and qualitative changes in the doublon speed, indicating a weak coupling between the doublon to magnetic excitations.

preprint2012arXiv

Magnetic and Orbital Order in ($R$MnO$_3$)$_n$/($A$MnO$_3$)$_{2n}$ Superlattices Studied via a Double-Exchange Model with Strain

The two-orbital double-exchange model is employed for the study of the magnetic and orbital orders in ($R$MnO$_3$)$_n$/($A$MnO$_3$)$_{2n}$ ($R$: rare earths; $A$: alkaline earths) superlattices. The A-type antiferromagnetic order is observed in a broad region of parameter space for the case of SrTiO$_3$ as substrate, in agreement with recent experiments and first-principles calculations using these superlattices. In addition, also a C-type antiferromagnetic state is predicted to be stabilized when using substrates like LaAlO$_3$ with smaller lattice constants than SrTiO$_3$, again in agreement with first principles results. The physical mechanism for the stabilization of the A- and C- magnetic transitions is driven by the orbital splitting of the $x^2-y^2$ and $3z^2-r^2$ orbitals. This splitting is induced by the $Q_3$ mode of Jahn-Teller distortions created by the strain induced by the substrates. In addition to the special example of (LaMnO$_3$)$_n$/(SrMnO$_3$)$_{2n}$, our phase diagrams can be valuable for the case where the superlattices are prepared employing narrow bandwidth manganites. In particular, several non-homogenous magnetic profiles are predicted to occur in narrow bandwidth superlattices, highlighting the importance of carrying out investigations in this mostly unexplored area of research.

preprint2012arXiv

Magnetic States of the Two-Leg Ladder Alkali Metal Iron Selenides $A$Fe$_2$Se$_3$

Recent neutron scattering experiments addressing the magnetic state of the two-leg ladder selenide compound BaFe$_2$Se$_3$ have unveiled a dominant spin arrangement involving ferromagnetically ordered 2$\times$2 iron-superblocks, that are antiferromagnetically coupled among them (the "block-AFM" state). Using the electronic five-orbital Hubbard model, first principles techniques to calculate the electronic hopping amplitudes between irons, and the real-space Hartree-Fock approximation to handle the many-body effects, here it is shown that the exotic block-AFM state is indeed stable at realistic electronic densities close to $n \sim 6.0$. Another state (the "CX" state) with parallel spins along the rungs and antiparallel along the legs of the ladders is close in energy. This state becomes stable in other portions of the phase diagrams, such as with hole doping, as also found experimentally via neutron scattering applied to KFe$_2$Se$_3$. In addition, the present study unveils other competing magnetic phases that could be experimentally stabilized varying either $n$ chemically or the electronic bandwidth by pressure. Similar results were obtained using two-orbital models, studied here via Lanczos and DMRG techniques. A comparison of the results obtained with the realistic selenides hoppings amplitudes for BaFe$_2$Se$_3$ against those found using the hopping amplitudes for pnictides reveals several qualitative similarities, particularly at intermediate and large Hubbard couplings.

preprint2012arXiv

Magnetism and its microscopic origin in iron-based high-temperature superconductors

High-temperature superconductivity in the iron-based materials emerges from, or sometimes coexists with, their metallic or insulating parent compound states. This is surprising since these undoped states display dramatically different antiferromagnetic (AF) spin arrangements and N$\rm \acute{e}$el temperatures. Although there is general consensus that magnetic interactions are important for superconductivity, much is still unknown concerning the microscopic origin of the magnetic states. In this review, progress in this area is summarized, focusing on recent experimental and theoretical results and discussing their microscopic implications. It is concluded that the parent compounds are in a state that is more complex than implied by a simple Fermi surface nesting scenario, and a dual description including both itinerant and localized degrees of freedom is needed to properly describe these fascinating materials.

preprint2012arXiv

RPA Analysis of a Two-orbital Model for the BiS2-based Superconductors

The random-phase approximation (RPA) is here applied to a two-orbital model for the BiS2-based superconductors that was recently proposed by Usui et al., arXiv:1207.3888. Varying the density of doped electrons per Bi site, n, in the range 0.46 < n < 1.0, the spin fluctuations promote competing A1g and B2g superconducting states with similar pairing strengths, in analogy with the A1g-B1g near degeneracy found also within RPA in models for pnictides. At these band fillings, two hole-pockets centered at (0,0) and (π,π) display nearly parallel Fermi Surface segments close to wavevector(π/2,π/2), whose distance increases with n. After introducing electronic interactions treated in the RPA, the inter-pocket nesting of these segments leads to pair scattering with a rather "local" character in k-space. The similarity between the A1g and B2g channels observed here should manifest in experiments on BiS2-based superconductors if the pairing is caused by spin fluctuations.

preprint2012arXiv

The Unexpected Properties of Alkali Metal Iron Selenide Superconductors

The iron-based superconductors that contain FeAs layers as the fundamental building block in the crystal structures have been rationalized in the past using ideas based on the Fermi Surface nesting of hole and electron pockets when in the presence of weak Hubbard $U$ interactions. This approach seemed appropriate considering the small values of the magnetic moments in the parent compounds and the clear evidence based on photoemission experiments of the required electron and hole pockets. However, recent results in the context of alkali metal iron selenides, with generic chemical composition $A_x$Fe$_{2-y}$Se$_2$ ($A$ = alkali element), have drastically challenged those previous ideas since at particular compositions $y$ the low-temperature ground states are insulating and display antiferromagnetic magnetic order with large iron magnetic moments. Moreover, angle resolved photoemission studies have revealed the absence of hole pockets at the Fermi level in these materials. The present status of this exciting area of research, with the potential to alter conceptually our understanding of the iron-based superconductors, is here reviewed, covering both experimental and theoretical investigations. Other recent related developments are also briefly reviewed, such as the study of selenide two-leg ladders and the discovery of superconductivity in a single layer of FeSe. The conceptual issues considered established for the alkali metal iron selenides, as well as the several issues that still require further work, are discussed in the text.

preprint2011arXiv

Ab Initio Study of the Intrinsic Exchange Bias at the SrRuO$_3$/SrMnO$_3$ Interface

In a recent publication (S. Dong et al., Phys. Rev. Lett.103, 127201 (2009)), two (related) mechanisms were proposed to understand the intrinsic exchange bias present in oxides heterostructures involving G-type antiferromagnetic perovskites. The first mechanism is driven by the Dzyaloshinskii-Moriya interaction, which is a spin-orbit coupling effect. The second is induced by the ferroelectric polarization, and it is only active in heterostructures involving multiferroics. Using the SrRuO$_3$/SrMnO$_3$ superlattice as a model system, density-functional calculations are here performed to verify the two proposals. This proof-of-principle calculation provides convincing evidence that qualitatively supports both proposals.

preprint2011arXiv

Charge Stripes in the Two-Orbital Hubbard Model for Pnictides

The two-orbital Hubbard model for the pnictides is studied numerically using the real-space Hartree-Fock approximation on finite clusters. Upon electron doping, states with a nonuniform distribution of charge are stabilized. The observed patterns correspond to charge stripes oriented perpendicular to the direction of the spin stripes of the undoped magnetic ground state. While these charge striped states are robust when the undoped state has a Hubbard gap, their existence when the intermediate-coupling magnetic metallic state of pnictides is doped was also observed for particular model parameters. Results for hole doping and implications for recent experiments that reported electronic nematic states and spin incommensurability are also briefly discussed.

preprint2011arXiv

Competing Pairing Symmetries in a Generalized Two-Orbital Model for the Pnictides

We introduce and study an extended "t-U-J" two-orbital model for the pnictides that includes Heisenberg terms deduced from the strong coupling expansion. Including these J terms explicitly allows us to enhance the strength of the (pi, 0)-(0, pi) spin order which favors the presence of tightly bound pairing states even in the small clusters that are here exactly diagonalized. The A1g and B2g pairing symmetries are found to compete in the realistic spin-ordered and metallic regime. The dynamical pairing susceptibility additionally unveils low-lying B1g states, suggesting that small changes in parameters may render any of the three channels stable.

preprint2011arXiv

Emergent dimensional reduction of the spin sector in a model for narrow-band manganites

The widely used Double-Exchange model for manganites is shown to support various "striped" phases at filling fractions 1/n (n=3, 4, 5...), in the previously unexplored regime of narrow bandwidth and small Jahn-Teller coupling. Working in two dimensions, our main result is that these stripes can be individually spin flipped without a physically relevant change in the energy, i.e., we find a large groundstate manifold with nearly degenerate energies. The two-dimensional spin system thus displays an unexpected dynamically generated dimensional reduction into decoupled one-dimensional stripes, even though the electronic states remain two-dimensional. Relations of our results with recent literature addressing compass models in quantum computing are discussed.

preprint2011arXiv

Microscopic model for the ferroelectric field effect in oxide heterostructures

A microscopic model Hamiltonian for the ferroelectric field effect is introduced for the study of oxide heterostructures with ferroelectric components. The long-range Coulomb interaction is incorporated as an electrostatic potential, solved self-consistently together with the charge distribution. A generic double-exchange system is used as the conducting channel, epitaxially attached to the ferroelectric gate. The observed ferroelectric screening effect, namely the charge accumulation/depletion near the interface, is shown to drive interfacial phase transitions that give rise to robust magnetoelectric responses and bipolar resistive switching, in qualitative agreement with previous density functional theory calculations. The model can be easily adapted to other materials by modifying the Hamiltonian of the conducting channel, and it is useful in simulating ferroelectric field effect devices particularly those involving strongly correlated electronic components where ab-initio techniques are difficult to apply.

preprint2011arXiv

Multiferroic Properties of CaMn$_7$O$_{12}$

We report that CaMn$_7$O$_{12}$ is a new magnetic multiferroic material. The appearance of a ferroelectric polarization coinciding with the magnetic phase transition ($\sim90$ K) suggests the presence of ferroelectricity induced by magnetism, further confirmed by its strong magnetoelectric response. With respect to other known magnetic multiferroics, CaMn$_7$O$_{12}$ displays attractive multiferroic properties, such as a high ferroelectric critical temperature and large polarization. More importantly, these results open a new avenue to search for magnetic multiferroics in the catalogue of doped oxides.

preprint2011arXiv

Optical Conductivity Anisotropy in the Undoped Three-Orbital Hubbard Model for the Pnictides

The resistivity anisotropy unveiled in the study of detwinned single crystals of the undoped 122 pnictides is here studied using the two-dimensional three-orbital Hubbard model in the mean-field approximation. Calculating the Drude weight in the x and y directions at zero temperature for a Q=(π,0) magnetically ordered state, the conductance along the antiferromagnetic direction is shown to be larger than along the ferromagnetic direction. This effect is caused by the suppression of the d_{yz} orbital at the Fermi surface, but additional insight based on the momentum dependence of the transitions induced by the current operator is provided. It is shown that the effective suppression of the inter-orbital hopping d_{xy} and d_{yz} along the y direction is the main cause of the anisotropy.

preprint2011arXiv

Pairing Symmetries of a Hole-Doped Extended Two-Orbital Model for the Pnictides

The hole-doped ground state of a recently introduced extended "t-U-J" two-orbital Hubbard model for the Fe-based superconductors is studied via exact diagonalization methods on small clusters. Similarly as in the previously studied case of electron doping, A. Nicholson et al., Phys. Rev. Lett. 106 21702 (2011), upon hole doping it is observed that there are several competing pairing symmetries including A_{1g}, B_{1g}, and B_{2g}. However, contrary to the electron-doped case, the ground state of the hole-doped state has pseudocrystal momentum k=(π,π) in the unfolded Brillouin zone. In the two Fe-atom per unit cell representation, this indicates that the ground state involves anti-bonding, rather than bonding, combinations of the orbitals of the two Fe atoms in the unit-cell. The lowest state with k=(0,0) has only a slightly higher energy. These results indicate that this simple two-orbital model may be useful to capture some subtle aspects of the hole-doped pnictides since calculations for the five-orbital model have unveiled a hole pocket centered at M (k=(π,π)) in the unfolded Brillouin zone.

preprint2011arXiv

Role of Degeneracy, Hybridization, and Nesting in the Properties of Multi-Orbital Systems

To understand the role that degeneracy, hybridization, and nesting play in the magnetic and pairing properties of multiorbital Hubbard models we here study numerically two types of two- orbital models, both with hole-like and electron-like Fermi surfaces (FS's) that are related by nesting vectors (π, 0) and (0, π). In one case the bands that determine the FS's arise from strongly hybridized degenerate dxz and dyz orbitals, while in the other the two bands are determined by non-degenerate and non-hybridized s-like orbitals. Using a variety of techniques, in the weak coupling regime it is shown that only the model with hybridized bands develops metallic magnetic order, while the other model exhibits an ordered excitonic orbital-transverse spin state that is insulating and does not have a local magnetization. However, both models display similar insulating magnetic stripe ordering in the strong coupling limit. These results indicate that nesting is a necessary but not sufficient condition for the development of ordered states with finite local magnetization in multiorbital Hubbard systems; the additional ingredient appears to be that the nested portions of the bands need to have the same orbital flavor. This condition can be achieved via strong hybridization of the orbitals in weak coupling or via the FS reconstruction induced by the Coulomb interactions in the strong coupling regime. This effect also impacts the pairing symmetry as demonstrated by the study of the dominant pairing channels for the two models.

preprint2011arXiv

Study of the magnetic state of K$_{0.8}$Fe$_{1.6}$Se$_2$ using the five-orbital Hubbard model in the Hartree-Fock approximation

Motivated by the recent discovery of Fe-based superconductors close to an antiferromagnetic insulator in the experimental phase diagram, here the five-orbital Hubbard model (without lattice distortions) is studied using the real-space Hartree-Fock approximation, employing a 10x10 Fe cluster with Fe vacancies in a $\sqrt{5}x\sqrt{5}$ pattern. Varying the Hubbard and Hund couplings, and at electronic density $n$=6.0, the phase diagram contains an insulating state with the same spin pattern as observed experimentally, involving 2x2 ferromagnetic plaquettes coupled with one another antiferromagnetically. The presence of local ferromagnetic tendencies is in qualitative agreement with Lanczos results for the three-orbital model also reported here. The magnetic moment ~ 3$μ_B$/Fe is in good agreement with experiments. Several other phases are also stabilized in the phase diagram, in agreement with recent calculations using phenomenological models.

preprint2010arXiv

Conducting Jahn-Teller Domain Walls in Undoped Manganites

We investigate the electronic properties of multi-domain configurations in models for undoped manganites by means of variational and Monte Carlo techniques. These materials display simultaneous Jahn-Teller distortions and magnetic ordering. We find that a band of electronic states appears associated with Jahn-Teller domain walls, and this band is localized in the direction perpendicular to the walls. The energy and width of this band depends on the conformational properties of the domain walls. At finite temperatures, the conductance along the domain walls, induced by the localized domain wall bands, is orders magnitude larger than in the bulk.

preprint2010arXiv

Constraints Imposed by Symmetry on Pairing Operators for the Pnictides

Considering model Hamiltonians that respect the symmetry properties of the pnictides, it is argued that pairing interactions that couple electrons at different orbitals with an orbital-dependent pairing strength inevitably lead to interband pairing matrix elements, at least in some regions of the Brillouin zone. Such interband pairing has not been considered of relevance in multiorbital systems in previous investigations. It is also observed that if, instead, a purely intraband pairing interaction is postulated, this requires that the pairing operator has the form $Δ^+(k)=f(k) \sum_α d^+_{k,α,\uparrow}d^+_{-k, α,\downarrow}$ where $α$ labels the orbitals considered in the model and f(k) arises from the spatial location of the coupled electrons or holes. This means that the gaps at two different Fermi surfaces involving momenta $k_F$ and $k'_F$ can only differ by the ratio $f(k_F)/ f(k'_F)$ and that electrons in different orbitals must be subject to the same pairing attraction, thus requiring fine tuning. These results suggest that previously neglected interband pairing tendencies could actually be of relevance in a microscopic description of the pairing mechanism in the pnictides.

preprint2010arXiv

Electron-lattice coupling and partial nesting as the origin of Fermi-Arcs in manganites

A tight-binding model for eg electrons coupled to Jahn-Teller lattice distortions is studied via unbiased Monte- Carlo simulations. By focusing on the periodicity of the Jahn Teller distortions, and the one-particle spectral function, our results clarify the physical origin of the Fermi-arcs phase observed in layered manganites. In a range of parameters where no broken symmetry phase exists, the nearly nested Fermi surface favors certain correlations between Jahn Teller distortions. The spectral weight near the Brillouin zone edge is suppressed, leading to the pseudogap in the density of states. We discuss the stability of this phase as a function of temperature and coupling strength for different hole dopings.

preprint2010arXiv

Highly anisotropic resistivities in the double-exchange model for strained manganites

The highly anisotropic resistivities in strained manganites are theoretically studied using the two-orbital double-exchange model. At the nanoscale, the anisotropic double-exchange and Jahn-Teller distortions are found to be responsible for the robust anisotropic resistivities observed here via Monte Carlo simulations. An unbalanced in the population of orbitals caused by strain is responsible for these effects. In contrast, the anisotropic superexchange is found to be irrelevant to explain our results. Our model study suggests that highly anisotropic resistivities could be present in a wide range of strained manganites, even without (sub)micrometer-scale phase separation. In addition, our calculations also confirm the formation of anisotropic clusters in phase-separated manganites, which magnifies the anisotropic resistivities.

preprint2010arXiv

Neutron and ARPES Constraints on the Couplings of the Multiorbital Hubbard Model for the Pnictides

The results of neutron scattering and angle-resolved photoemission experiments for the Fe-pnictide parent compounds, and their metallic nature, are shown to impose severe constraints on the range of values that can be considered "realistic" for the intraorbital Hubbard repulsion U and Hund coupling J in multiorbital Hubbard models treated in the mean-field approximation. Phase diagrams for three- and five-orbital models are here provided, and the physically realistic regime of couplings is highlighted, to guide future theoretical work into the proper region of parameters of Hubbard models. In addition, using the random phase approximation, the pairing tendencies in these realistic coupling regions are investigated. It is shown that the dominant spin-singlet pairing channels in these coupling regimes correspond to nodal superconductivity, with strong competition between several states that belong to different irreducible representations. This is compatible with experimental bulk measurements that have reported the existence of nodes in several Fe-pnictide compounds.

preprint2010arXiv

Orbital-weight redistribution triggered by spin order in the pnictides

The one-particle spectral function and its orbital composition are investigated in a three-orbital model for the undoped parent compounds of the iron-based superconductors. In the realistic parameter regime, where results best fit experimental data, it is observed that the magnetization in the xz and yz orbitals are markedly different and the Fermi surface presents mostly xz character, as recently observed in photoemission experiments [T. Shimojima et al., Phys. Rev. Lett. 104, 057002 (2010)]. Since the ferro-orbital order in this regime is at most a few percent, these results are mainly driven by the magnetic order. An analogous analysis for a five-orbital model leads to similar conclusions.

preprint2010arXiv

Real-time dynamics of particle-hole excitations in Mott insulator-metal junctions

Charge excitations in Mott insulators (MIs) are distinct from their band-insulator counterparts and they can provide a mechanism for energy harvesting in solar cells based on strongly correlated electronic materials. In this paper, we study the real-time dynamics of holon-doublon pairs in a MI connected to metallic leads using the time-dependent density matrix renormalization group method. The transfer of charge across the MI-metal interface is controlled by both the electron-electron interaction strength within the MI and the voltage difference between the leads. We find an overall enhancement of the charge transfer as compared to the case of a (noninteracting) band insulator-metal interface with a similar band gap. Moreover, the propagation of holon-doublon excitations within the MI dynamically changes the spin-spin correlations, introducing time-dependent phase shifts in the spin structure factor.

preprint2010arXiv

Three orbital model for the iron-based superconductors

The theoretical need to study the properties of the Fe-based high-T_c superconductors with reliable many-body techniques requires us to determine the minimum number of orbital degrees of freedom that will capture the physics of these materials. While the shape of the Fermi surface (FS) obtained with the local density approximation (LDA) can be reproduced by a two-orbital model, it has been argued that the bands that cross the chemical potential result from the strong hybridization of three of the Fe 3d orbitals. For this reason, a three-orbital Hamiltonian obtained with the Slater-Koster formalism by considering the hybridization of the As p orbitals with the Fe d_xz,d_yz, and d_xy orbitals is discussed here. This model reproduces qualitatively the FS shape and orbital composition obtained by LDA calculations for undoped pnictides when four electrons per Fe are considered. Within a mean-field approximation, its magnetic and orbital properties in the undoped case are described. With increasing Coulomb repulsion, four regimes are obtained: (1) paramagnetic, (2) magnetic (pi,0) spin order, (3) the same (pi,0) spin order but now including orbital order, and finally (4) a magnetic and orbital ordered insulator. The spin-singlet pairing operators allowed by the lattice and orbital symmetries are also constructed. It is found that for pairs of electrons involving up to diagonal nearest-neighbors sites, the only fully gapped and purely intraband spin-singlet pairing operator is given by Delta(k)=f(k)\sum_{alpha} d_{k,alpha,up}d_{-k,alpha,down} with f(k)=1 or f(k)=cos(k_x)cos(k_y) which would arise only if the electrons in all different orbitals couple with equal strength to the source of pairing.

preprint2010arXiv

Unveiling First Order CMR Transitions in the Two-Orbital Model for Manganites

Large scale Monte Carlo simulation results for the two-orbital model for manganites, including Jahn-Teller lattice distortions, are here presented. At hole density x=1/4 and in the vicinity of the region of competition between the ferromagnetic metallic and spin-charge-orbital ordered insulating phases, the colossal magnetoresistance (CMR) phenomenon is observed with a magnetoresistance ratio ~10,000% Our main result is that this CMR transition is found to be of first order in some portions of the phase diagram, in agreement with early results from neutron scattering, specific heat, and magnetization, thus solving a notorious discrepancy between experiments and previous theoretical studies. The first-order characteristics of the transition survive, and are actually enhanced, when weak quenched disorder is introduced.

preprint2009arXiv

Double-exchange model study of multiferroic $R$MnO$_3$ perovskites

In this proceeding, recent theoretical investigations by the authors on the multiferroic $R$MnO$_3$ perovskites are briefly reviewed at first. Using the double-exchange model, the realistic spiral spin order in undoped manganites such as TbMnO$_3$ and DyMnO$_3$ is well reproduced by incorporating a weak next-nearest neighbor superexchange ($\sim10%$ of nearest neighbor superexchange) and moderate Jahn-Teller distortion. The phase transitions from the A-type antiferromagnet (as in LaMnO$_3$), to the spiral phase (as in TbMnO$_3$), and finally to the E-type antiferromagnet (as in HoMnO$_3$), with decreasing size of the $R$ ions, were also explained. Moreover, new results of phase diagram of the three-dimensional lattice are also included. The ferromagnetic tendency recently discovered in the LaMnO$_3$ and TbMnO$_3$ thin films is explained by considering the substrate stress. Finally, the relationship between our double-exchange model and a previously used $J_1$-$J_2$-$J_3$ model is further discussed from the perspective of spin wave excitations.

preprint2009arXiv

Electronic and Magnetic Properties of $\rm{RMnO}_3/\rm{AMnO}_3$ Heterostructures

In this paper, we investigate the ground state properties of the $\rm{RMnO}_3/\rm{AMnO}_3$ ($\rm{RMO/AMO}$) heterostructures (R=trivalent cation, and A=divalent cation) by using a two-orbital double-exchange model supplemented by the Poisson's equation. We find that the state stabilized near the interface of the heterostructure is similar to the state of the bulk compound $\rm{(R,A)MO}$ at electronic density close to 0.5. Depending on the bandwidth it will be a charge and orbital ordered CE state or an A-AF state with $x^2-y^2$ orbital order. Our results can explain some properties of long-period superlattices. As another interesting result, we find exotic intermediate states stabilized in between the interface and the bulk-like regimes of the heterostructure. For instance, a spin "canted CE" state and others. They may not have an analog in experimentally known bulk phase diagrams, but provide a natural interpolation between magnetically-ordered states that are stable in the bulk at different electronic densities.

preprint1999arXiv

Hole-Density Evolution of the One-Particle Spectral Function in Doped Ladders

The spectral function $A({\bf q}, ω)$ of doped $t-J$ ladders is presented on clusters with up to $2 \times 20$ sites at zero temperature applying a recently developed technique that uses up to $\sim 6 \times 10^6$ rung-basis states. Similarities with photoemission results for the 2D cuprates are observed, such as the existence of a gap at $(π,0)$ near half-filling (caused by hole pair formation) and flat bands in its vicinity. These features should be observable in ARPES experiments on ladders. The main result of the paper is the nontrivial evolution of the spectral function from a narrow band at $x=0$, to a quasi-noninteracting band at $x \geq 0.5$. It was also observed that the low-energy peaks of a cluster spectra acquire finite line-widths as their energies move away from the chemical potential.

preprint1999arXiv

Indications of Spin-Charge Separation in the 2D Extended t-J Model

The 2D extended $\rm t-J$ model is studied computationally in a broad region of parameter space, motivated by recent photoemission experiments for the undoped cuprate $\rm Ca_2 Cu O_2 Cl_2$ (F. Ronning et al., Science {\bf 282}, 2067 (1998)). The one-hole ground state is shown to develop robust antiferromagnetic (AF) correlations between spins separated by the mobile hole (i.e. across the hole). This effect tends to decouple charge from spin, and the quasiparticle weight becomes negligible, particularly at momenta $(0,π)-(π,0)$. Studies with more holes show precursors of metallic stripe formation, with holes sharing their individual spin arrangements, and AF-correlations generated across the stripe.

preprint1997arXiv

Enhancement of Antiferromagnetic Correlations Induced by Nonmagnetic Impurities: Origin and Predictions for NMR Experiments

Spin models that have been proposed to describe dimerized chains, ladders, two dimensional antiferromagnets, and other compounds are here studied when some spins are replaced by spinless vacancies, such as it occurs by $Zn$ doping. A small percentage of vacancies rapidly destroys the spin gap, and their presence induces enhanced antiferromagnetic correlations near those vacancies. The study is performed with computational techniques which includes Lanczos, world-line Monte Carlo, and the Density Matrix Renormalization Group methods. Since the phenomenon of enhanced antiferromagnetism is found to occur in several models and cluster geometries, a common simple explanation for its presence may exist. It is argued that the resonating-valence-bond character of the spin correlations at short distances of a large variety of models is responsible for the presence of robust staggered spin correlations near vacancies and lattice edges. The phenomenon takes place regardless of the long distance properties of the ground state, and it is caused by a ``pruning'' of the available spin singlets in the vicinity of the vacancies. The effect produces a broadening of the low temperature NMR signal for the compounds analyzed here. This broadening should be experimentally observable in the structurally dimerized chain systems $Cu(NO_3)_2\cdot2.5H_2O$, $CuWO_4$, $(VO)_2P_2O_7$, and $Sr_{14}Cu_{24}O_{41}$, in ladder materials such as $Sr Cu_2 O_3$, in the spin-Peierls systems $CuGeO_3$ and $NaV_2 O_5$, and in several others since it is a universal effect common to a wide variety of models and compounds.

preprint1997arXiv

Hole Doping Evolution of the Quasiparticle Band in Models of Strongly Correlated Electrons for the High-T_c Cuprates

Quantum Monte Carlo (QMC) and Maximum Entropy (ME) techniques are used to study the spectral function $A({\bf p},ω)$ of the one band Hubbard model in strong coupling including a next-nearest-neighbor electronic hopping with amplitude $t'/t= -0.35$. These values of parameters are chosen to improve the comparison of the Hubbard model with angle-resolved photoemission (ARPES) data for $Sr_2 Cu O_2 Cl_2$. A narrow quasiparticle (q.p.) band is observed in the QMC analysis at the temperature of the simulation $T=t/3$, both at and away from half-filling. Such a narrow band produces a large accumulation of weight in the density of states at the top of the valence band. As the electronic density $< n >$ decreases further away from half-filling, the chemical potential travels through this energy window with a large number of states, and by $< n > \sim 0.70$ it has crossed it entirely. The region near momentum $(0,π)$ and $(π,0)$ in the spectral function is more sensitive to doping than momenta along the diagonal from $(0,0)$ to $(π,π)$. The evolution with hole density of the quasiparticle dispersion contains some of the features observed in recent ARPES data in the underdoped regime. For sufficiently large hole densities the ``flat'' bands at $(π,0)$ cross the Fermi energy, a prediction that could be tested with ARPES techniques applied to overdoped cuprates. The population of the q.p. band introduces a {\it hidden} density in the system which produces interesting consequences when the quasiparticles are assumed to interact through antiferromagnetic fluctuations and studied with the BCS gap equation formalism. In particular, a region of extended s-wave is found to compete with d-wave in the overdoped regime, i.e. when the chemical potential has almost entirely crossed the q.p.

preprint1997arXiv

Local Enhancement of Antiferromagnetic Correlations by Nonmagnetic Impurities

The local enhancement of antiferromagnetic correlations near vacancies observed in a variety of spin systems is analyzed in a single framework. Variational calculations suggest that the resonating-valence-bond character of the spin correlations at short distances is responsible for the enhancement. Numerical results for uniform spin chains, with and without frustration, dimerized chains, ladders, and two dimensional clusters are in agreement with our conjecture. This short distance phenomenon occurs independently of the long distance behavior of the spin correlations in the undoped system. Experimental predictions for a variety of compounds are briefly discussed.

preprint1996arXiv

Ferromagnetism in Electronic Models for Manganites

Ground state properties of the Kondo model for manganese oxides in one dimension are studied using numerical techniques. The large Hund coupling ($J_{H}$) limit is specially analyzed. A robust region of fully saturated ferromagnetism (FM) is identified at all densities. For open boundary conditions it is shown exactly that the ground state is FM at $J_{H} = \infty$. Hole-spin phase separation competing with FM was also observed when a large exchange $J$ between the $Mn^{3+}$ ions is used. As the spin of the transition metal ion grows, the hole mobility decreases providing a tentative explanation for the differences between Cu-oxides and Mn-oxides.

preprint1996arXiv

Photoemission Spectra in t-J Ladders with Two Legs

Photoemission spectra for the isotropic two-leg t-J ladder are calculated at various hole-doping levels using exact diagonalization techniques. Low-energy sharp features caused by short-range antiferromagnetic correlations are observed at finite doping levels close to half-filling, above the naive Fermi momentum. These features should be observable in angle-resolved photoemission experiments. In addition, the formation of a d-wave pairing condensate as the ratio J/t is increased leads to dynamically generated spectral weight for momenta close to $k_F$ where the $d_{x^2-y^2 }$-order parameter is large.

preprint1996arXiv

Rapid Suppression of the Spin Gap in Zn-doped CuGeO_3 and SrCu_2O_3

The influence of non-magnetic impurities on the spectrum and dynamical spin structure factor of a model for CuGeO$_3$ is studied. A simple extension to Zn-doped ${\rm Sr Cu_2 O_3}$ is also discussed. Using Exact Diagonalization techniques and intuitive arguments we show that Zn-doping introduces states in the Spin-Peierls gap of CuGeO$_3$. This effect can beunderstood easily in the large dimerization limit where doping by Zn creates ``loose'' S=1/2 spins, which interact with each other through very weak effective antiferromagnetic couplings. When the dimerization is small, a similar effect is observed but now with the free S=1/2 spins being the resulting S=1/2 ground state of severed chains with an odd number of sites. Experimental consequences of these results are discussed. It is interesting to observe that the spin correlations along the chains are enhanced by Zn-doping according to the numerical data presented here. As recent numerical calculations have shown, similar arguments apply to ladders with non-magnetic impurities simply replacing the tendency to dimerization in CuGeO$_3$ by the tendency to form spin-singlets along the rungs in SrCu$_2$O$_3$.

preprint1995arXiv

Magnetic Raman Scattering in Two-Dimensional Spin-1/2 Heisenberg Antiferromagnets: Spectral Shape Anomaly and Magnetostrictive Effects

We calculate the Raman spectrum of the two-dimensional (2D) spin-1/2 Heisenberg antiferromagnet by exact diagonalization and quantum Monte Carlo techniques on clusters of up to 144 sites and, on a 16-site cluster, by considering the phonon-magnon interaction which leads to random fluctuations of the exchange integral. Results are in good agreement with experiments on various high-T_c precursors, such as La_2CuO_4 and YBa_2Cu_3O_{6.2}. In particular, our calculations reproduce the broad lineshape of the two-magnon peak, the asymmetry about its maximum, the existence of spectral weight at high energies, and the observation of nominally forbidden A_{1g} scattering.

preprint1993arXiv

Random Exchange Disorder in the Spin-1/2 XXZ Chain

The one-dimensional XXZ model is studied in the presence of disorder in the Heisenberg Exchange Integral. Recent predictions obtained from renormalization group calculations are investigated numerically using a Lanczos algorithm on chains of up to 18 sites. It is found that in the presence of strong X-Y-symmetric random exchange couplings, a ``random singlet'' phase with quasi-long-range order in the spin-spin correlations persists. As the planar anisotropy is varied, the full zero-temperature phase diagram is obtained and compared with predictions of Doty and Fisher [Phys. Rev. B {\bf 45 }, 2167 (1992)].