Source author record

Adriana Moreo

Adriana Moreo 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

37works
3topics
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

37 published item(s)

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

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.

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

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.

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

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.

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

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

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.

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

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

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

Spectral density in a nematic state of iron pnictides

Using cluster-perturbation theory, we calculate the spectral density A(k,w) for a nematic phase of models describing pnictide superconductors, where very short-range magnetic correlations choose the ordering vector (pi,0) over the equivalent (0,pi) and thus break the fourfold rotation symmetry of the underlying lattice without inducing long-range magnetic order. In excellent agreement with angle resolved photo-emission spectroscopy (ARPES), we find that the yz bands at X move to higher energies. When onsite Coulomb repulsion brings the system close to a spin--density-wave (SDW) and renormalizes the band width by a factor of approx. 2, even small anisotropic couplings of 10 to 15 meV strongly distort the bands, splitting the formerly degenerate states at X and Y by approx. 70 meV and shifting the yz states at X above the chemical potential. This similarity to the SDW bands is in excellent agreement with ARPES. An important difference to the SDW bands is that the yz bands still cross the Fermi level, again in agreement with experiment. We find that orbital weights near the Fermi surface provide a better characterization than overall orbital densities and orbital polarization.

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

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

Comment on "Nonmagnetic Impurity Resonances as a Signature of Sign-Reversal Pairing in Fe-As-Based Superconductors"

In a recent Letter [PRL 103, 186402 (2009)], the energy band structure of Fe-As-based superconductors is fitted with a tight-binding model with two Fe ions per unit cell and two degenerate xz and yz orbitals per Fe ion. The author claims that the proposed model, which differs markedly from a model previously used by other authors for the same two orbitals in the same compounds, possesses the symmetry required to describe the Fe-As planes in iron-pnictide superconductors. In this comment we argue that this is not the case.

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

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

Spin-Polarized Semiconductor Induced by Magnetic Impurities in Graphene

Magnetic impurities adsorbed on graphene are coupled magnetically via the itinerant electrons. This interaction opens a gap in the band structure of graphene. The result strongly depends on how the magnetic impurities are distributed. While random doping produces a semiconductor, if all or most impurities are located in the same sublattice, the spin degeneracy is removed and a spin-polarized semiconductor arises.

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.

preprint2005arXiv

Effect of Adiabatic Phonons on Striped and Homogeneous Ground States

The effects of adiabatic phonons on a spin-fermion model for high T_c cuprates are studied using numerical simulations. In the absence of electron-phonon interactions (EPI), stripes in the ground state are observed for certain dopings while homogeneous states are stabilized in other regions of parameter space. Different modes of adiabatic phonons are added to the Hamiltonian:breathing, shear and half-breathing modes. Diagonal and off-diagonal electron-phonon couplings are considered. It is observed that strong diagonal EPI generate stripes in previously homogeneous states, while in striped ground states an increase in the diagonal couplings tends to stabilize the stripes, inducing a gap in the density of states (DOS) and rendering the ground state insulating. The off-diagonal terms, on the other hand, destabilize the stripes creating inhomogeneous ground states with a pseudogap at the chemical potential in the DOS. The breathing mode stabilizes static diagonal stripes; while the half-breathing (shear) modes stabilize dynamical (localized) vertical and horizontal stripes. The EPI induces decoherence of the quasi-particle peaks in the spectral functions.

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.