-
Intertwined superconductivity and orbital selectivity in a three-orbital Hubbard model for the iron pnictides
Authors:
Vito Marino,
Alberto Scazzola,
Federico Becca,
Massimo Capone,
Luca F. Tocchio
Abstract:
We study a three-orbital Hubbard-Kanamori model relevant for iron-based superconductors using variational wave functions, which explicitly include spatial correlations and electron pairing. We span the nonmagnetic sector ranging from a filling $n=4$, which is representative of undoped iron-based superconductors, to $n=3$. In the latter case, a Mott insulating state is found, with each orbital at h…
▽ More
We study a three-orbital Hubbard-Kanamori model relevant for iron-based superconductors using variational wave functions, which explicitly include spatial correlations and electron pairing. We span the nonmagnetic sector ranging from a filling $n=4$, which is representative of undoped iron-based superconductors, to $n=3$. In the latter case, a Mott insulating state is found, with each orbital at half filling. In the strong-coupling regime, when the electron density is increased, we find a spontaneous differentiation between the occupation of $d_{xz}$ and $d_{yz}$ orbitals, which leads to an orbital-selective state with a nematic character that becomes stronger at increasing density. One of these orbitals stays half filled for all densities while the other one hosts (together with the $d_{xy}$ orbital) the excess of electron density. Most importantly, in this regime, long-range pairing correlations appear in the orbital with the largest occupation. Our results highlight a strong link between orbital-selective correlations, nematicity, and superconductivity, which requires the presence of a significant Hund's coupling.
△ Less
Submitted 19 June, 2024;
originally announced June 2024.
-
Variational Monte Carlo study of stripes as a function of doping in the $t-t'$ Hubbard model
Authors:
Antonio Lechiara,
Vito Marino,
Luca F. Tocchio
Abstract:
We perform variational Monte Carlo simulations of the single-band Hubbard model on the square lattice with both nearest ($t$) and next-nearest ($t'$) neighbor hoppings. Our work investigates the consequences of increasing hole doping on the instauration of stripes and the behavior of the superconducting order parameter, with a discussion on how the two phenomena affect each other. We consider two…
▽ More
We perform variational Monte Carlo simulations of the single-band Hubbard model on the square lattice with both nearest ($t$) and next-nearest ($t'$) neighbor hoppings. Our work investigates the consequences of increasing hole doping on the instauration of stripes and the behavior of the superconducting order parameter, with a discussion on how the two phenomena affect each other. We consider two different values of the next-nearest neighbor hopping parameter, that are appropriate for describing cuprate superconductors. We observe that stripes are the optimal state in a wide doping range; the stripe wavelength reduces at increasing doping, until stripes melt into a uniform state for large values of doping. Superconducting pair-pair correlations, indicating the presence of superconductivity, are always suppressed in the presence of stripes. Our results suggest that the phase diagram for the single-band Hubbard model is dominated by stripes, with superconductivity being possible only in a narrow doping range between striped states and a nonsuperconducting metal.
△ Less
Submitted 8 July, 2024; v1 submitted 4 March, 2024;
originally announced March 2024.
-
Variational Benchmarks for Quantum Many-Body Problems
Authors:
Dian Wu,
Riccardo Rossi,
Filippo Vicentini,
Nikita Astrakhantsev,
Federico Becca,
Xiaodong Cao,
Juan Carrasquilla,
Francesco Ferrari,
Antoine Georges,
Mohamed Hibat-Allah,
Masatoshi Imada,
Andreas M. Läuchli,
Guglielmo Mazzola,
Antonio Mezzacapo,
Andrew Millis,
Javier Robledo Moreno,
Titus Neupert,
Yusuke Nomura,
Jannes Nys,
Olivier Parcollet,
Rico Pohle,
Imelda Romero,
Michael Schmid,
J. Maxwell Silvester,
Sandro Sorella
, et al. (8 additional authors not shown)
Abstract:
The continued development of novel many-body approaches to ground-state problems in physics and chemistry calls for a consistent way to assess its overall progress. Here we introduce a metric of variational accuracy, the V-score, obtained from the variational energy and its variance. We provide the most extensive curated dataset of variational calculations of many-body quantum systems to date, ide…
▽ More
The continued development of novel many-body approaches to ground-state problems in physics and chemistry calls for a consistent way to assess its overall progress. Here we introduce a metric of variational accuracy, the V-score, obtained from the variational energy and its variance. We provide the most extensive curated dataset of variational calculations of many-body quantum systems to date, identifying cases where state-of-the-art numerical approaches show limited accuracy, and novel algorithms or computational platforms, such as quantum computing, could provide improved accuracy. The V-score can be used as a metric to assess the progress of quantum variational methods towards quantum advantage for ground-state problems, especially in regimes where classical verifiability is impossible.
△ Less
Submitted 9 February, 2023;
originally announced February 2023.
-
Stripes in the extended $t-t^\prime$ Hubbard model: A Variational Monte Carlo analysis
Authors:
Vito Marino,
Federico Becca,
Luca F. Tocchio
Abstract:
By using variational quantum Monte Carlo techniques, we investigate the instauration of stripes (i.e., charge and spin inhomogeneities) in the Hubbard model on the square lattice at hole doping $δ=1/8$, with both nearest- ($t$) and next-nearest-neighbor hopping ($t^\prime$). Stripes with different wavelengths $λ$ (denoting the periodicity of the charge inhomogeneity) and character (bond- or site-c…
▽ More
By using variational quantum Monte Carlo techniques, we investigate the instauration of stripes (i.e., charge and spin inhomogeneities) in the Hubbard model on the square lattice at hole doping $δ=1/8$, with both nearest- ($t$) and next-nearest-neighbor hopping ($t^\prime$). Stripes with different wavelengths $λ$ (denoting the periodicity of the charge inhomogeneity) and character (bond- or site-centered) are stabilized for sufficiently large values of the electron-electron interaction $U/t$. The general trend is that $λ$ increases going from negative to positive values of $t^\prime/t$ and decreases by increasing $U/t$. In particular, the $λ=8$ stripe obtained for $t^\prime=0$ and $U/t=8$ [L.F. Tocchio, A. Montorsi, and F. Becca, SciPost Phys. {\bf 7}, 21 (2019)] shrinks to $λ=6$ for $U/t\gtrsim 10$. For $t^\prime/t<0$, the stripe with $λ=5$ is found to be remarkably stable, while for $t^\prime/t>0$, stripes with wavelength $λ=12$ and $λ=16$ are also obtained. In all these cases, pair-pair correlations are highly suppressed with respect to the uniform state (obtained for large values of $|t^\prime/t|$), suggesting that striped states are not superconducting at $δ=1/8$.
△ Less
Submitted 3 June, 2022; v1 submitted 8 November, 2021;
originally announced November 2021.
-
The Hubbard model on triangular $N$-leg cylinders: chiral and non-chiral spin liquids
Authors:
Luca F. Tocchio,
Arianna Montorsi,
Federico Becca
Abstract:
The existence of a gapped chiral spin liquid has been recently suggested in the vicinity of the metal-insulator transition of the Hubbard model on the triangular lattice, by intensive density-matrix renormalization group (DMRG) simulations [A. Szasz, J. Motruk, M.P. Zaletel, and J.E. Moore, Phys. Rev. X {\bf 10}, 021042 (2020)]. Here, we report the results obtained within the variational Monte Car…
▽ More
The existence of a gapped chiral spin liquid has been recently suggested in the vicinity of the metal-insulator transition of the Hubbard model on the triangular lattice, by intensive density-matrix renormalization group (DMRG) simulations [A. Szasz, J. Motruk, M.P. Zaletel, and J.E. Moore, Phys. Rev. X {\bf 10}, 021042 (2020)]. Here, we report the results obtained within the variational Monte Carlo technique based upon Jastrow-Slater wave functions, implemented with backflow correlations. As in DMRG calculations, we consider $N$-leg cylinders. For $N=4$ and in the presence of a next-nearest neighbor hopping, a chiral spin liquid emerges between the metal and the insulator with magnetic quasi-long-range order. Within our approach, the chiral state is gapped and breaks the reflection symmetry. By contrast, for both $N=5$ and $N=6$, the chiral spin liquid is not the state with the lowest variational energy: in the former case a nematic spin liquid is found in the entire insulating regime, while for the less frustrated case with $N=6$ the results are very similar to the one obtained on two-dimensional clusters [L.F. Tocchio, A. Montorsi, and F. Becca, Phys. Rev. B {\bf 102}, 115150 (2020)], with an antiferromagnetic phase close to the metal-insulator transition and a nematic spin liquid in the strong-coupling regime.
△ Less
Submitted 15 October, 2021; v1 submitted 3 May, 2021;
originally announced May 2021.
-
Magnetic and spin-liquid phases in the frustrated $t-t^\prime$ Hubbard model on the triangular lattice
Authors:
Luca F. Tocchio,
Arianna Montorsi,
Federico Becca
Abstract:
The Hubbard model and its strong-coupling version, the Heisenberg one, have been widely studied on the triangular lattice to capture the essential low-temperature properties of different materials. One example is given by transition metal dichalcogenides, as 1T$-$TaS$_2$, where a large unit cell with $13$ Ta atom forms weakly-coupled layers with an isotropic triangular lattice. By using accurate v…
▽ More
The Hubbard model and its strong-coupling version, the Heisenberg one, have been widely studied on the triangular lattice to capture the essential low-temperature properties of different materials. One example is given by transition metal dichalcogenides, as 1T$-$TaS$_2$, where a large unit cell with $13$ Ta atom forms weakly-coupled layers with an isotropic triangular lattice. By using accurate variational Monte Carlo calculations, we report the phase diagram of the $t-t^\prime$ Hubbard model on the triangular lattice, highlighting the differences between positive and negative values of $t^\prime/t$; this result can be captured only by including the charge fluctuations that are always present for a finite electron-electron repulsion. Two spin-liquid regions are detected: one for $t^\prime/t<0$, which persists down to intermediate values of the electron-electron repulsion, and a narrower one for $t^\prime/t>0$. The spin-liquid phase appears to be gapless, though the variational wave function has a nematic character, in contrast to the Heisenberg limit. We do not find any evidence for non-magnetic Mott phases in the proximity of the metal-insulator transition, at variance with the predictions (mainly based upon strong-coupling expansions in $t/U$) that suggest the existence of a weak-Mott phase that intrudes between the metal and the magnetically ordered insulator.
△ Less
Submitted 17 September, 2020; v1 submitted 19 May, 2020;
originally announced May 2020.
-
Metallic and insulating stripes and their relation with superconductivity in the doped Hubbard model
Authors:
Luca F. Tocchio,
Arianna Montorsi,
Federico Becca
Abstract:
The dualism between superconductivity and charge/spin modulations (the so-called stripes) dominates the phase diagram of many strongly-correlated systems. A prominent example is given by the Hubbard model, where these phases compete and possibly coexist in a wide regime of electron dopings for both weak and strong couplings. Here, we investigate this antagonism within a variational approach that i…
▽ More
The dualism between superconductivity and charge/spin modulations (the so-called stripes) dominates the phase diagram of many strongly-correlated systems. A prominent example is given by the Hubbard model, where these phases compete and possibly coexist in a wide regime of electron dopings for both weak and strong couplings. Here, we investigate this antagonism within a variational approach that is based upon Jastrow-Slater wave functions, including backflow correlations, which can be treated within a quantum Monte Carlo procedure. We focus on clusters having a ladder geometry with $M$ legs (with $M$ ranging from $2$ to $10$) and a relatively large number of rungs, thus allowing us a detailed analysis in terms of the stripe length. We find that stripe order with periodicity $λ=8$ in the charge and $2λ=16$ in the spin can be stabilized at doping $δ=1/8$. Here, there are no sizable superconducting correlations and the ground state has an insulating character. A similar situation, with $λ=6$, appears at $δ=1/6$. Instead, for smaller values of dopings, stripes can be still stabilized, but they are weakly metallic at $δ=1/12$ and metallic with strong superconducting correlations at $δ=1/10$, as well as for intermediate (incommensurate) dopings. Remarkably, we observe that spin modulation plays a major role in stripe formation, since it is crucial to obtain a stable striped state upon optimization. The relevance of our calculations for previous density-matrix renormalization group results and for the two-dimensional case is also discussed.
△ Less
Submitted 15 July, 2019; v1 submitted 7 May, 2019;
originally announced May 2019.
-
Nematic and Antiferromagnetic Quantum Criticality in a Multi-Orbital Hubbard Model for Iron Pnictides
Authors:
Wen-Jun Hu,
Haoyu Hu,
Rong Yu,
Hsin-Hua Lai,
Luca F. Tocchio,
Federico Becca,
Qimiao Si
Abstract:
The extent to which quantum criticality drives the physics of iron pnictides is a central question in the field. Earlier theoretical considerations were based on an effective field theory, and the proposed realization in P-doped iron arsenides has received extensive experimental evidence. To connect the quantum critical behavior with the underlying electronic physics, it is important to analyze it…
▽ More
The extent to which quantum criticality drives the physics of iron pnictides is a central question in the field. Earlier theoretical considerations were based on an effective field theory, and the proposed realization in P-doped iron arsenides has received extensive experimental evidence. To connect the quantum critical behavior with the underlying electronic physics, it is important to analyze it within microscopic models. Here, we do so for a multi-orbital model containing both Hubbard and Hund's interactions, by a variational Monte Carlo method based on Jastrow-Slater wave functions that allow for a non-perturbative treatment of electron correlations. We find strong evidence for the existence of a unique quantum critical point, where both nematic and $(π,0)$ antiferromagnetic orders develop together, in the bad-metal regime of the phase diagram. Implications of our results for the iron-based superconductivity are discussed.
△ Less
Submitted 29 March, 2019;
originally announced March 2019.
-
Superconductivity in the Hubbard model: a hidden-order diagnostics from the Luther-Emery phase on ladders
Authors:
Luca F. Tocchio,
Federico Becca,
Arianna Montorsi
Abstract:
Short-range antiferromagnetic correlations are known to open a spin gap in the repulsive Hubbard model on ladders with $M$ legs, when $M$ is even. We show that the spin gap originates from the formation of correlated pairs of electrons with opposite spin, captured by the hidden ordering of a spin-parity operator. Since both spin gap and parity vanish in the two-dimensional limit, we introduce the…
▽ More
Short-range antiferromagnetic correlations are known to open a spin gap in the repulsive Hubbard model on ladders with $M$ legs, when $M$ is even. We show that the spin gap originates from the formation of correlated pairs of electrons with opposite spin, captured by the hidden ordering of a spin-parity operator. Since both spin gap and parity vanish in the two-dimensional limit, we introduce the fractional generalization of spin parity and prove that it remains finite in the thermodynamic limit. Our results are based upon variational wave functions and Monte Carlo calculations: performing a finite size-scaling analysis with growing $M$, we show that the doping region where the parity is finite coincides with the range in which superconductivity is observed in two spatial dimensions. Our observations support the idea that superconductivity emerges out of spin gapped phases on ladders, driven by a spin-pairing mechanism, in which the ordering is conveniently captured by the finiteness of the fractional spin-parity operator.
△ Less
Submitted 22 January, 2019; v1 submitted 29 October, 2018;
originally announced October 2018.
-
Metal-insulator transitions, superconductivity, and magnetism in the two-band Hubbard model
Authors:
Caterina De Franco,
Luca F. Tocchio,
Federico Becca
Abstract:
We explore the ground-state properties of the two-band Hubbard model with degenerate electronic bands, parametrized by nearest-neighbor hopping $t$, intra- and inter-orbital on-site Coulomb repulsions $U$ and $U^\prime$, and Hund coupling $J$, focusing on the case with $J>0$. Using Jastrow-Slater wave functions, we consider both states with and without magnetic/orbital order. Electron pairing can…
▽ More
We explore the ground-state properties of the two-band Hubbard model with degenerate electronic bands, parametrized by nearest-neighbor hopping $t$, intra- and inter-orbital on-site Coulomb repulsions $U$ and $U^\prime$, and Hund coupling $J$, focusing on the case with $J>0$. Using Jastrow-Slater wave functions, we consider both states with and without magnetic/orbital order. Electron pairing can also be included in the wave function, in order to detect the occurrence of superconductivity for generic electron densities $n$. When no magnetic/orbital order is considered, the Mott transition is continuous for $n=1$ (quarter filling); instead, at $n=2$ (half filling), it is first order for small values of $J/U$, while it turns out to be continuous when the ratio $J/U$ is increased. A significant triplet pairing is present in a broad region around $n=2$. By contrast, singlet superconductivity (with $d$-wave symmetry) is detected only for small values of the Hund coupling and very close to half filling. When including magnetic and orbital order, the Mott insulator acquires antiferromagnetic order for $n=2$; instead, for $n=1$ the insulator has ferromagnetic and antiferro-orbital orders. In the latter case, a metallic phase is present for small values of $U/t$ and the metal-insulator transition becomes first order. In the region with $1<n<2$, we observe that ferromagnetism (with no orbital order) is particularly robust for large values of the Coulomb repulsion and that triplet superconductivity is strongly suppressed by the presence of antiferromagnetism. The case with $J=0$, which has an enlarged SU(4) symmetry due to the interplay between spin and orbital degrees of freedom, is also analyzed.
△ Less
Submitted 29 August, 2018; v1 submitted 4 May, 2018;
originally announced May 2018.
-
Superconductivity, charge-density waves, antiferromagnetism, and phase separation in the Hubbard-Holstein model
Authors:
Seher Karakuzu,
Luca F. Tocchio,
Sandro Sorella,
Federico Becca
Abstract:
By using variational wave functions and quantum Monte Carlo techniques, we investigate the interplay between electron-electron and electron-phonon interactions in the two-dimensional Hubbard-Holstein model. Here, the ground-state phase diagram is triggered by several energy scales, i.e., the electron hopping $t$, the on-site electron-electron interaction $U$, the phonon energy $ω_0$, and the elect…
▽ More
By using variational wave functions and quantum Monte Carlo techniques, we investigate the interplay between electron-electron and electron-phonon interactions in the two-dimensional Hubbard-Holstein model. Here, the ground-state phase diagram is triggered by several energy scales, i.e., the electron hopping $t$, the on-site electron-electron interaction $U$, the phonon energy $ω_0$, and the electron-phonon coupling $g$. At half filling, the ground state is an antiferromagnetic insulator for $U \gtrsim 2g^2/ω_0$, while it is a charge-density-wave (or bi-polaronic) insulator for $U \lesssim 2g^2/ω_0$. In addition to these phases, we find a superconducting phase that intrudes between them. For $ω_0/t=1$, superconductivity emerges when both $U/t$ and $2g^2/tω_0$ are small; then, by increasing the value of the phonon energy $ω_0$, it extends along the transition line between antiferromagnetic and charge-density-wave insulators. Away from half filling, phase separation occurs when doping the charge-density-wave insulator, while a uniform (superconducting) ground state is found when doping the superconducting phase. In the analysis of finite-size effects, it is extremely important to average over twisted boundary conditions, especially in the weak-coupling limit and in the doped case.
△ Less
Submitted 27 November, 2017; v1 submitted 1 September, 2017;
originally announced September 2017.
-
Charge orders in organic charge-transfer salts
Authors:
Ryui Kaneko,
Luca F. Tocchio,
Roser Valentí,
Federico Becca
Abstract:
Motivated by recent experimental suggestions of charge-order-driven ferroelectricity in organic charge-transfer salts, such as $κ$-(BEDT-TTF)$_2$Cu[N(CN)$_2$]Cl, we investigate magnetic and charge-ordered phases that emerge in an extended two-orbital Hubbard model on the anisotropic triangular lattice at $3/4$ filling. This model takes into account the presence of two organic BEDT-TTF molecules, w…
▽ More
Motivated by recent experimental suggestions of charge-order-driven ferroelectricity in organic charge-transfer salts, such as $κ$-(BEDT-TTF)$_2$Cu[N(CN)$_2$]Cl, we investigate magnetic and charge-ordered phases that emerge in an extended two-orbital Hubbard model on the anisotropic triangular lattice at $3/4$ filling. This model takes into account the presence of two organic BEDT-TTF molecules, which form a dimer on each site of the lattice, and includes short-range intramolecular and intermolecular interactions and hoppings. By using variational wave functions and quantum Monte Carlo techniques, we find two polar states with charge disproportionation inside the dimer, hinting to ferroelectricity. These charge-ordered insulating phases are stabilized in the strongly correlated limit and their actual charge pattern is determined by the relative strength of intradimer to interdimer couplings. Our results suggest that ferroelectricity is not driven by magnetism, since these polar phases can be stabilized also without antiferromagnetic order and provide a possible microscopic explanation of the experimental observations. In addition, a conventional dimer-Mott state (with uniform density and antiferromagnetic order) and a nonpolar charge-ordered state (with charge-rich and charge-poor dimers forming a checkerboard pattern) can be stabilized in the strong-coupling regime. Finally, when electron-electron interactions are weak, metallic states appear, with either uniform charge distribution or a peculiar $12$-site periodicity that generates honeycomb-like charge order.
△ Less
Submitted 27 October, 2017; v1 submitted 24 May, 2017;
originally announced May 2017.
-
Hidden Mott transition and large-$U$ superconductivity in the two-dimensional Hubbard model
Authors:
Luca F. Tocchio,
Federico Becca,
Sandro Sorella
Abstract:
We consider the one-band Hubbard model on the square lattice by using variational and Green's function Monte Carlo methods, where the variational states contain Jastrow and backflow correlations on top of an uncorrelated wave function that includes BCS pairing and magnetic order. At half filling, where the ground state is antiferromagnetically ordered for any value of the on-site interaction $U$,…
▽ More
We consider the one-band Hubbard model on the square lattice by using variational and Green's function Monte Carlo methods, where the variational states contain Jastrow and backflow correlations on top of an uncorrelated wave function that includes BCS pairing and magnetic order. At half filling, where the ground state is antiferromagnetically ordered for any value of the on-site interaction $U$, we can identify a hidden critical point $U_{\rm Mott}$, above which a finite BCS pairing is stabilized in the wave function. The existence of this point is reminiscent of the Mott transition in the paramagnetic sector and determines a separation between a Slater insulator (at small values of $U$), where magnetism induces a potential energy gain, and a Mott insulator (at large values of $U$), where magnetic correlations drive a kinetic energy gain. Most importantly, the existence of $U_{\rm Mott}$ has crucial consequences when doping the system: We observe a tendency to phase separation into a hole-rich and a hole-poor region only when doping the Slater insulator, while the system is uniform by doping the Mott insulator. Superconducting correlations are clearly observed above $U_{\rm Mott}$, leading to the characteristic dome structure in doping. Furthermore, we show that the energy gain due to the presence of a finite BCS pairing above $U_{\rm Mott}$ shifts from the potential to the kinetic sector by increasing the value of the Coulomb repulsion.
△ Less
Submitted 16 November, 2016; v1 submitted 22 July, 2016;
originally announced July 2016.
-
Emergent lattices with geometrical frustration in doped extended Hubbard models
Authors:
Ryui Kaneko,
Luca F. Tocchio,
Roser Valentí,
Claudius Gros
Abstract:
Spontaneous charge ordering occurring in correlated systems may be considered as a possible route to generate effective lattice structures with unconventional couplings. For this purpose we investigate the phase diagram of doped extended Hubbard models on two lattices: (i) the honeycomb lattice with on-site $U$ and nearest-neighbor $V$ Coulomb interactions at $3/4$ filling ($n=3/2$) and (ii) the t…
▽ More
Spontaneous charge ordering occurring in correlated systems may be considered as a possible route to generate effective lattice structures with unconventional couplings. For this purpose we investigate the phase diagram of doped extended Hubbard models on two lattices: (i) the honeycomb lattice with on-site $U$ and nearest-neighbor $V$ Coulomb interactions at $3/4$ filling ($n=3/2$) and (ii) the triangular lattice with on-site $U$, nearest-neighbor $V$, and next-nearest-neighbor $V'$ Coulomb interactions at $3/8$ filling ($n=3/4$). We consider various approaches including mean-field approximations, perturbation theory, and variational Monte Carlo. For the honeycomb case (i), charge order induces an effective triangular lattice at large values of $U/t$ and $V/t$, where $t$ is the nearest-neighbor hopping integral. The nearest-neighbor spin exchange interactions on this effective triangular lattice are antiferromagnetic in most of the phase diagram, while they become ferromagnetic when $U$ is much larger than $V$. At $U/t\sim (V/t)^3$, ferromagnetic and antiferromagnetic exchange interactions nearly cancel out, leading to a system with four-spin ring-exchange interactions. On the other hand, for the triangular case (ii) at large $U$ and finite $V'$, we find no charge order for small $V$, an effective kagome lattice for intermediate $V$, and one-dimensional charge order for large $V$. These results indicate that Coulomb interactions induce [case (i)] or enhance [case(ii)] emergent geometrical frustration of the spin degrees of freedom in the system, by forming charge order.
△ Less
Submitted 14 November, 2016; v1 submitted 21 June, 2016;
originally announced June 2016.
-
Variational wave functions for the $S=1/2$ Heisenberg model on the anisotropic triangular lattice: Spin liquids and spiral orders
Authors:
Elaheh Ghorbani,
Luca F. Tocchio,
Federico Becca
Abstract:
By using variational wave functions and quantum Monte Carlo techniques, we investigate the complete phase diagram of the Heisenberg model on the anisotropic triangular lattice, where two out of three bonds have super-exchange couplings $J$ and the third one has instead $J^\prime$. This model interpolates between the square lattice and the isotropic triangular one, for $J^\prime/J \le 1$, and betwe…
▽ More
By using variational wave functions and quantum Monte Carlo techniques, we investigate the complete phase diagram of the Heisenberg model on the anisotropic triangular lattice, where two out of three bonds have super-exchange couplings $J$ and the third one has instead $J^\prime$. This model interpolates between the square lattice and the isotropic triangular one, for $J^\prime/J \le 1$, and between the isotropic triangular lattice and a set of decoupled chains, for $J/J^\prime \le 1$. We consider all the fully-symmetric spin liquids that can be constructed with the fermionic projective-symmetry group classification [Y. Zhou and X.-G. Wen, arXiv:cond-mat/0210662] and we compare them with the spiral magnetic orders that can be accommodated on finite clusters. Our results show that, for $J^\prime/J \le 1$, the phase diagram is dominated by magnetic orderings, even though a spin-liquid state may be possible in a small parameter window, i.e., $0.7 \lesssim J^\prime/J \lesssim 0.8$. In contrast, for $J/J^\prime \le 1$, a large spin-liquid region appears close to the limit of decoupled chains, i.e., for $J/J^\prime \lesssim 0.6$, while magnetically ordered phases with spiral order are stabilized close to the isotropic point.
△ Less
Submitted 5 February, 2016; v1 submitted 10 December, 2015;
originally announced December 2015.
-
Spontaneous symmetry breaking in correlated wave functions
Authors:
Ryui Kaneko,
Luca F. Tocchio,
Roser Valentí,
Federico Becca,
Claudius Gros
Abstract:
We show that Jastrow-Slater wave functions, in which a density-density Jastrow factor is applied onto an uncorrelated fermionic state, may possess long-range order even when all symmetries are preserved in the wave function. This fact is mainly related to the presence of a sufficiently strong Jastrow term (also including the case of full Gutzwiller projection, suitable for describing spin models).…
▽ More
We show that Jastrow-Slater wave functions, in which a density-density Jastrow factor is applied onto an uncorrelated fermionic state, may possess long-range order even when all symmetries are preserved in the wave function. This fact is mainly related to the presence of a sufficiently strong Jastrow term (also including the case of full Gutzwiller projection, suitable for describing spin models). Selected examples are reported, including the spawning of Néel order and dimerization in spin systems, and the stabilization of charge and orbital order in itinerant electronic systems.
△ Less
Submitted 20 March, 2016; v1 submitted 29 October, 2015;
originally announced October 2015.
-
Assessing the orbital selective Mott transition with variational wave functions
Authors:
Luca F. Tocchio,
Federico Arrigoni,
Sandro Sorella,
Federico Becca
Abstract:
We study the Mott metal-insulator transition in the two-band Hubbard model with different hopping amplitudes $t_1$ and $t_2$ for the two orbitals on the two-dimensional square lattice by using {\it non-magnetic} variational wave functions, similarly to what has been considered in the limit of infinite dimensions by dynamical mean-field theory. We work out the phase diagram at half filling (i.e., t…
▽ More
We study the Mott metal-insulator transition in the two-band Hubbard model with different hopping amplitudes $t_1$ and $t_2$ for the two orbitals on the two-dimensional square lattice by using {\it non-magnetic} variational wave functions, similarly to what has been considered in the limit of infinite dimensions by dynamical mean-field theory. We work out the phase diagram at half filling (i.e., two electrons per site) as a function of $R=t_2/t_1$ and the on-site Coulomb repulsion $U$, for two values of the Hund's coupling $J=0$ and $J/U=0.1$. Our results are in good agreement with previous dynamical mean-field theory calculations, demonstrating that the non-magnetic phase diagram is only slightly modified from infinite to two spatial dimensions. Three phases are present: a metallic one, for small values of $U$, where both orbitals are itinerant; a Mott insulator, for large values of $U$, where both orbitals are localized because of the Coulomb repulsion; and the so-called orbital-selective Mott insulator (OSMI), for small values of $R$ and intermediate $U$'s, where one orbital is localized while the other one is still itinerant. The effect of the Hund's coupling is two-fold: on one side, it favors the full Mott phase over the OSMI; on the other side, it stabilizes the OSMI at larger values of $R$.
△ Less
Submitted 17 February, 2016; v1 submitted 26 May, 2015;
originally announced May 2015.
-
Solutions of the Two Dimensional Hubbard Model: Benchmarks and Results from a Wide Range of Numerical Algorithms
Authors:
J. P. F. LeBlanc,
Andrey E. Antipov,
Federico Becca,
Ireneusz W. Bulik,
Garnet Kin-Lic Chan,
Chia-Min Chung,
Youjin Deng,
Michel Ferrero,
Thomas M. Henderson,
Carlos A. Jiménez-Hoyos,
E. Kozik,
Xuan-Wen Liu,
Andrew J. Millis,
N. V. Prokof'ev,
Mingpu Qin,
Gustavo E. Scuseria,
Hao Shi,
B. V. Svistunov,
Luca F. Tocchio,
I. S. Tupitsyn,
Steven R. White,
Shiwei Zhang,
Bo-Xiao Zheng,
Zhenyue Zhu,
Emanuel Gull
Abstract:
Numerical results for ground state and excited state properties (energies, double occupancies, and Matsubara-axis self energies) of the single-orbital Hubbard model on a two-dimensional square lattice are presented, in order to provide an assessment of our ability to compute accurate results in the thermodynamic limit. Many methods are employed, including auxiliary field quantum Monte Carlo, bare…
▽ More
Numerical results for ground state and excited state properties (energies, double occupancies, and Matsubara-axis self energies) of the single-orbital Hubbard model on a two-dimensional square lattice are presented, in order to provide an assessment of our ability to compute accurate results in the thermodynamic limit. Many methods are employed, including auxiliary field quantum Monte Carlo, bare and bold-line diagrammatic Monte Carlo, method of dual fermions, density matrix embedding theory, density matrix renormalization group, dynamical cluster approximation, diffusion Monte Carlo within a fixed node approximation, unrestricted coupled cluster theory, and multi-reference projected Hartree-Fock. Comparison of results obtained by different methods allows for the identification of uncertainties and systematic errors. The importance of extrapolation to converged thermodynamic limit values is emphasized. Cases where agreement between different methods is obtained establish benchmark results that may be useful in the validation of new approaches and the improvement of existing methods.
△ Less
Submitted 15 December, 2015; v1 submitted 9 May, 2015;
originally announced May 2015.
-
One-dimensional spin liquid, collinear, and spiral phases from uncoupled chains to the triangular lattice
Authors:
Luca F. Tocchio,
Claudius Gros,
Roser Valentí,
Federico Becca
Abstract:
We investigate the Hubbard model on the anisotropic triangular lattice with two hopping parameters $t$ and $t^\prime$ in different spatial directions, interpolating between decoupled chains ($t=0$) and the isotropic triangular lattice ($t=t^\prime$). Variational wave functions that include both Jastrow and backflow terms are used to compare spin-liquid and magnetic phases with different pitch vect…
▽ More
We investigate the Hubbard model on the anisotropic triangular lattice with two hopping parameters $t$ and $t^\prime$ in different spatial directions, interpolating between decoupled chains ($t=0$) and the isotropic triangular lattice ($t=t^\prime$). Variational wave functions that include both Jastrow and backflow terms are used to compare spin-liquid and magnetic phases with different pitch vectors describing both collinear and coplanar (spiral) order. For relatively large values of the on-site interaction $U/t^\prime \gtrsim 10$ and substantial frustration, i.e., $0.3\lesssim t/t^\prime \lesssim 0.8$, the spin-liquid state is clearly favored over magnetic states. Spiral magnetic order is only stable in the vicinity of the isotropic point, while collinear order is obtained in a wide range of inter-chain hoppings from small to intermediate frustration.
△ Less
Submitted 23 June, 2014; v1 submitted 18 March, 2014;
originally announced March 2014.
-
Phase diagram of the triangular extended Hubbard model
Authors:
Luca F. Tocchio,
Claudius Gros,
Xue-Feng Zhang,
Sebastian Eggert
Abstract:
We study the extended Hubbard model on the triangular lattice as a function of filling and interaction strength. The complex interplay of kinetic frustration and strong interactions on the triangular lattice leads to exotic phases where long-range charge order, antiferromagnetic order, and metallic conductivity can coexist. Variational Monte Carlo simulations show that three kinds of ordered metal…
▽ More
We study the extended Hubbard model on the triangular lattice as a function of filling and interaction strength. The complex interplay of kinetic frustration and strong interactions on the triangular lattice leads to exotic phases where long-range charge order, antiferromagnetic order, and metallic conductivity can coexist. Variational Monte Carlo simulations show that three kinds of ordered metallic states are stable as a function of nearest neighbor interaction and filling. The coexistence of conductivity and order is explained by a separation into two functional classes of particles: part of them contributes to the stable order, while the other part forms a partially filled band on the remaining substructure. The relation to charge ordering in charge transfer salts is discussed.
△ Less
Submitted 13 February, 2014;
originally announced February 2014.
-
Phase diagram of the square lattice bilayer Hubbard model: a variational Monte Carlo study
Authors:
Robert Rüger,
Luca F. Tocchio,
Roser Valentí,
Claudius Gros
Abstract:
We investigate the phase diagram of the square lattice bilayer Hubbard model at half filling with the variational Monte Carlo method for both the magnetic and the paramagnetic case as a function of interlayer hopping t_perp and on-site Coulomb repulsion U. With this study we resolve some discrepancies in previous calculations based on the dynamical mean field theory, and we are able to determine t…
▽ More
We investigate the phase diagram of the square lattice bilayer Hubbard model at half filling with the variational Monte Carlo method for both the magnetic and the paramagnetic case as a function of interlayer hopping t_perp and on-site Coulomb repulsion U. With this study we resolve some discrepancies in previous calculations based on the dynamical mean field theory, and we are able to determine the nature of the phase transitions between metal, Mott insulator and band insulator. In the magnetic case we find only two phases: An antiferromagnetic Mott insulator at small t_perp for any value of U and a band insulator at large t_perp. At large U values we approach the Heisenberg limit. The paramagnetic phase diagram shows at small t_perp a metal to Mott insulator transition at moderate U values and a Mott to band insulator transition at larger U values. We also observe a reentrant Mott insulator to metal transition and metal to band insulator transition for increasing t_perp in the range of 5.5t < U < 7.5t. Finally, we discuss the obtained phase diagrams in relation to previous studies based on different many-body approaches.
△ Less
Submitted 23 February, 2014; v1 submitted 25 November, 2013;
originally announced November 2013.
-
Importance of anisotropy in the spin-liquid candidate Me3EtSb[Pd(dmit)2]2
Authors:
A. C. Jacko,
Luca F. Tocchio,
Harald O. Jeschke,
Roser Valenti
Abstract:
Organic charge transfer salts based on the molecule Pd(dmit)2 display strong electronic correlations and geometrical frustration, leading to spin liquid, valence bond solid, and superconducting states, amongst other interesting phases. The low energy electronic degrees of freedom of these materials are often described by a single band model; a triangular lattice with a molecular orbital representi…
▽ More
Organic charge transfer salts based on the molecule Pd(dmit)2 display strong electronic correlations and geometrical frustration, leading to spin liquid, valence bond solid, and superconducting states, amongst other interesting phases. The low energy electronic degrees of freedom of these materials are often described by a single band model; a triangular lattice with a molecular orbital representing a Pd(dmit)2 dimer on each site. We use ab initio electronic structure calculations to construct and parametrize low energy effective model Hamiltonians for a class of Me(4-n) EtnX[Pd(dmit)2]2 (X=As,P,N,Sb) salts and investigate how best to model these systems by using variational Monte Carlo (VMC) simulations. Our findings suggest that the prevailing model of these systems as a t-t' triangular lattice is incomplete, and that a fully anisotropic triangular lattice (FATL) description produces importantly different results, including a significant lowering of the critical U of the spin-liquid phase.
△ Less
Submitted 14 November, 2013; v1 submitted 21 August, 2013;
originally announced August 2013.
-
Spin-liquid versus spiral-order phases in the anisotropic triangular lattice
Authors:
Luca F. Tocchio,
Hélène Feldner,
Federico Becca,
Roser Valentí,
Claudius Gros
Abstract:
We study the competition between magnetic and spin-liquid phases in the Hubbard model on the anisotropic triangular lattice, which is described by two hopping parameters t and t' in different spatial directions and is relevant for layered organic charge-transfer salts. By using a variational approach that includes spiral magnetic order, we provide solid evidence that a spin-liquid phase is stabili…
▽ More
We study the competition between magnetic and spin-liquid phases in the Hubbard model on the anisotropic triangular lattice, which is described by two hopping parameters t and t' in different spatial directions and is relevant for layered organic charge-transfer salts. By using a variational approach that includes spiral magnetic order, we provide solid evidence that a spin-liquid phase is stabilized in the strongly-correlated regime and close to the isotropic limit t'/t=1. Otherwise, a magnetically ordered spiral state is found, connecting the (collinear) Néel and the (coplanar) 120° phases. The pitch vector of the spiral phase obtained from the unrestricted Hartree-Fock approximation is substantially renormalized in presence of electronic correlations, and the Néel phase is stabilized in a wide regime of the phase diagram, i.e., for t'/t < 0.75. We discuss these results in the context of organic charge-transfer salts.
△ Less
Submitted 31 January, 2013; v1 submitted 10 September, 2012;
originally announced September 2012.
-
Mott correlated states in the underdoped two-dimensional Hubbard model: variational Monte Carlo versus a dynamical cluster approximation
Authors:
Luca F. Tocchio,
Hunpyo Lee,
Harald O. Jeschke,
Roser Valentí,
Claudius Gros
Abstract:
We investigate the properties of the frustrated underdoped Hubbard model on the square lattice using two complementary approaches, the dynamical cluster extension of dynamical mean field theory, and variational Monte Carlo simulations of Gutzwiller-Jastrow wavefunctions with backflow corrections. We compare and discuss data for the energy and the double occupancies, as obtained from both approache…
▽ More
We investigate the properties of the frustrated underdoped Hubbard model on the square lattice using two complementary approaches, the dynamical cluster extension of dynamical mean field theory, and variational Monte Carlo simulations of Gutzwiller-Jastrow wavefunctions with backflow corrections. We compare and discuss data for the energy and the double occupancies, as obtained from both approaches. At small dopings, we observe a rapid crossover from a weakly correlated metal at low interaction strength U to a non-Fermi liquid correlated state with strong local spin correlations. Furthermore, we investigate the stability of the correlated state against phase separation. We observe phase separation only for large values of U or very large frustration. No phase separation is present for the parameter range relevant for the cuprates.
△ Less
Submitted 20 December, 2012; v1 submitted 26 July, 2012;
originally announced July 2012.
-
Strong renormalization of the Fermi-surface topology close to the Mott transition
Authors:
Luca F. Tocchio,
Federico Becca,
Claudius Gros
Abstract:
The underlying Fermi surface is a key concept for strongly-interacting electron models and has been introduced to generalize the usual notion of the Fermi surface to generic (superconducting or insulating) systems. By using improved correlated wave functions that contain backflow and Jastrow terms, we examine the two-dimensional $t{-}t^\prime$ Hubbard model and find a non-trivial renormalization o…
▽ More
The underlying Fermi surface is a key concept for strongly-interacting electron models and has been introduced to generalize the usual notion of the Fermi surface to generic (superconducting or insulating) systems. By using improved correlated wave functions that contain backflow and Jastrow terms, we examine the two-dimensional $t{-}t^\prime$ Hubbard model and find a non-trivial renormalization of the topology of the underlying Fermi surface close to the Mott insulator. Moreover, we observe a sharp crossover region, which arises from the metal-insulator transition, from a weakly interacting metal at small coupling to a resonating valence-bond superconductor at intermediate coupling. A violation of the Luttinger theorem is detected at low hole dopings.
△ Less
Submitted 14 May, 2012; v1 submitted 11 January, 2012;
originally announced January 2012.
-
Backflow correlations in the Hubbard model: an efficient tool for the metal-insulator transition and the large-U limit
Authors:
Luca F. Tocchio,
Federico Becca,
Claudius Gros
Abstract:
We show that backflow correlations in the variational wave function for the Hubbard model greatly improve the previous results given by the Slater-Jastrow state, usually considered in this context. We provide evidence that, within this approach, it is possible to have a satisfactory connection with the strong-coupling regime. Moreover, we show that, for the Hubbard model on the lattice, backflow c…
▽ More
We show that backflow correlations in the variational wave function for the Hubbard model greatly improve the previous results given by the Slater-Jastrow state, usually considered in this context. We provide evidence that, within this approach, it is possible to have a satisfactory connection with the strong-coupling regime. Moreover, we show that, for the Hubbard model on the lattice, backflow correlations are essentially short range, inducing an effective attraction between empty (holons) and doubly occupied sites (doublons). In presence of frustration, we report the evidence that the metal to Mott-insulator transition is marked by a discontinuity of the double occupancy, together with a similar discontinuity of the kinetic term that does not change the number of holons and doublons, while the other kinetic terms are continuous across the transition. Finally, we show the estimation of the charge gap, obtained by particle-hole excitations {\it à la Feynman} over the ground-state wave function.
△ Less
Submitted 15 February, 2011;
originally announced February 2011.
-
Tunnelling matrix elements with antiferromagnetic Gutzwiller wave functions
Authors:
Andrea Di Ciolo,
Luca F. Tocchio,
Claudius Gros
Abstract:
We use a generalized Gutzwiller Approximation (GA) elaborated to evaluate matrix elements with partially projected wave functions and formerly applied to homogeneous systems.
In the present paper we consider projected single-particle (hole) excitations for electronic systems with antiferromagnetic (AFM) order and obtain the corresponding tunnelling probabilities. The accuracy and the reliability…
▽ More
We use a generalized Gutzwiller Approximation (GA) elaborated to evaluate matrix elements with partially projected wave functions and formerly applied to homogeneous systems.
In the present paper we consider projected single-particle (hole) excitations for electronic systems with antiferromagnetic (AFM) order and obtain the corresponding tunnelling probabilities. The accuracy and the reliability of our analytical approximation is tested using the Variational Monte Carlo (VMC). Possible comparisons with experimental results are also discussed.
△ Less
Submitted 25 October, 2010;
originally announced October 2010.
-
Interaction induced Fermi-surface renormalization in the $t_1{-}t_2$ Hubbard model close to the Mott-Hubbard transition
Authors:
Luca F. Tocchio,
Federico Becca,
Claudius Gros
Abstract:
We investigate the nature of the interaction-driven Mott-Hubbard transition of the half-filled $t_1{-}t_2$ Hubbard model in one dimension, using a full-fledged variational Monte Carlo approach including a distance-dependent Jastrow factor and backflow correlations. We present data for the evolution of the magnetic properties across the Mott-Hubbard transition and on the commensurate to incommens…
▽ More
We investigate the nature of the interaction-driven Mott-Hubbard transition of the half-filled $t_1{-}t_2$ Hubbard model in one dimension, using a full-fledged variational Monte Carlo approach including a distance-dependent Jastrow factor and backflow correlations. We present data for the evolution of the magnetic properties across the Mott-Hubbard transition and on the commensurate to incommensurate transition in the insulating state. Analyzing renormalized excitation spectra, we find that the Fermi surface renormalizes to perfect nesting right at the Mott-Hubbard transition in the insulating state, with a first-order reorganization when crossing into the conducting state.
△ Less
Submitted 11 May, 2010; v1 submitted 21 December, 2009;
originally announced December 2009.
-
Spin-liquid and magnetic phases in the anisotropic triangular lattice: the case of $κ$-(ET)$_2$X
Authors:
Luca F. Tocchio,
Alberto Parola,
Claudius Gros,
Federico Becca
Abstract:
The two-dimensional Hubbard model on the anisotropic triangular lattice, with two different hopping amplitudes $t$ and $t^\prime$, is relevant to describe the low-energy physics of $κ$-(ET)$_2$X, a family of organic salts. The ground-state properties of this model are studied by using Monte Carlo techniques, on the basis of a recent definition of backflow correlations for strongly-correlated lat…
▽ More
The two-dimensional Hubbard model on the anisotropic triangular lattice, with two different hopping amplitudes $t$ and $t^\prime$, is relevant to describe the low-energy physics of $κ$-(ET)$_2$X, a family of organic salts. The ground-state properties of this model are studied by using Monte Carlo techniques, on the basis of a recent definition of backflow correlations for strongly-correlated lattice systems. The results show that there is no magnetic order for reasonably large values of the electron-electron interaction $U$ and frustrating ratio $t^\prime/t = 0.85$, suitable to describe the non-magnetic compound with X=Cu$_2$(CN)$_3$. On the contrary, Néel order takes place for weaker frustrations, i.e., $t^\prime/t \sim 0.4 ÷0.6$, suitable for materials with X=Cu$_2$(SCN)$_2$, Cu[N(CN)$_2$]Cl, or Cu[N(CN)$_2$]Br.
△ Less
Submitted 4 September, 2009; v1 submitted 12 June, 2009;
originally announced June 2009.
-
Metal-insulator transition and strong-coupling spin liquid in the $t{-}t^\prime$ Hubbard model
Authors:
Federico Becca,
Luca F. Tocchio,
Sandro Sorella
Abstract:
We study the phase diagram of the frustrated $t{-}t^\prime$ Hubbard model on the square lattice by using a novel variational wave function. Taking the clue from the backflow correlations that have been introduced long-time ago by Feynman and Cohen and have been used for describing various interacting systems on the continuum (like liquid $^3$He, the electron jellium, and metallic Hydrogen), we c…
▽ More
We study the phase diagram of the frustrated $t{-}t^\prime$ Hubbard model on the square lattice by using a novel variational wave function. Taking the clue from the backflow correlations that have been introduced long-time ago by Feynman and Cohen and have been used for describing various interacting systems on the continuum (like liquid $^3$He, the electron jellium, and metallic Hydrogen), we consider many-body correlations to construct a suitable approximation for the ground state of this correlated model on the lattice. In this way, a very accurate {\it ansatz} can be achieved both at weak and strong coupling. We present the evidence that an insulating and non-magnetic phase can be stabilized at strong coupling and sufficiently large frustrating ratio $t^\prime/t$.
△ Less
Submitted 3 October, 2008;
originally announced October 2008.
-
Role of backflow correlations for the non-magnetic phase of the t-t' Hubbard model
Authors:
Luca F. Tocchio,
Federico Becca,
Alberto Parola,
Sandro Sorella
Abstract:
We introduce an efficient way to improve the accuracy of projected wave functions, widely used to study the two-dimensional Hubbard model. Taking the clue from the backflow contribution, whose relevance has been emphasized for various interacting systems on the continuum, we consider many-body correlations to construct a suitable approximation for the ground state at intermediate and strong coup…
▽ More
We introduce an efficient way to improve the accuracy of projected wave functions, widely used to study the two-dimensional Hubbard model. Taking the clue from the backflow contribution, whose relevance has been emphasized for various interacting systems on the continuum, we consider many-body correlations to construct a suitable approximation for the ground state at intermediate and strong couplings. In particular, we study the phase diagram of the frustrated $t{-}t^\prime$ Hubbard model on the square lattice and show that, thanks to backflow correlations, an insulating and non-magnetic phase can be stabilized at strong coupling and sufficiently large frustrating ratio $t^\prime/t$.
△ Less
Submitted 10 May, 2008;
originally announced May 2008.