-
Insulating and metallic phases in the one-dimensional Hubbard-Su-Schrieffer-Heeger model: Insights from a backflow-inspired variational wave function
Authors:
Davide Piccioni,
Francesco Ferrari,
Michele Fabrizio,
Federico Becca
Abstract:
The interplay between electron-electron and electron-phonon interactions is studied in a one-dimensional lattice model, by means of a variational Monte Carlo method based on generalized Jastrow-Slater wave functions. Here, the fermionic part is constructed by a pair-product state, which explicitly depends on the phonon configuration, thus including the electron-phonon coupling in a backflow-inspir…
▽ More
The interplay between electron-electron and electron-phonon interactions is studied in a one-dimensional lattice model, by means of a variational Monte Carlo method based on generalized Jastrow-Slater wave functions. Here, the fermionic part is constructed by a pair-product state, which explicitly depends on the phonon configuration, thus including the electron-phonon coupling in a backflow-inspired way. We report the results for the Hubbard model in presence of the Su-Schrieffer-Heeger coupling to optical phonons, both at half-filling and upon hole doping. At half-filling, the ground state is either a translationally invariant Mott insulator, with gapless spin excitations, or a Peierls insulator, which breaks translations and has fully gapped excitations. Away from half-filling, the charge gap closes in both Mott and Peierls insulators, turning the former into a conventional Luttinger liquid (gapless in all excitation channels). The latter, instead, retains a finite spin gap that closes only above a threshold value of the doping. Even though consistent with the general theory of interacting electrons in one dimension, the existence of such a phase (with gapless charge but gapped spin excitations) has never been demonstrated in a model with repulsive interaction and with only two Fermi points. Since the spin-gapped metal represents the one-dimensional counterpart of a superconductor, our results furnish evidence that a true off-diagonal long-range order may exist in the two-dimensional case.
△ Less
Submitted 3 July, 2024;
originally announced July 2024.
-
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.
-
Fine-tuning Neural Network Quantum States
Authors:
Riccardo Rende,
Sebastian Goldt,
Federico Becca,
Luciano Loris Viteritti
Abstract:
Recent progress in the design and optimization of Neural-Network Quantum States (NQS) have made them an effective method to investigate ground-state properties of quantum many-body systems. In contrast to the standard approach of training a separate NQS from scratch at every point of the phase diagram, we demonstrate that the optimization of a NQS at a highly expressive point of the phase diagram…
▽ More
Recent progress in the design and optimization of Neural-Network Quantum States (NQS) have made them an effective method to investigate ground-state properties of quantum many-body systems. In contrast to the standard approach of training a separate NQS from scratch at every point of the phase diagram, we demonstrate that the optimization of a NQS at a highly expressive point of the phase diagram (i.e., close to a phase transition) yields interpretable features that can be reused to accurately describe a wide region across the transition. We demonstrate the feasibility of our approach on different systems in one and two dimensions by initially pretraining a NQS at a given point of the phase diagram, followed by fine-tuning only the output layer for all other points. Notably, the computational cost of the fine-tuning step is very low compared to the pretraining stage. We argue that the reduced cost of this paradigm has significant potential to advance the exploration of condensed matter systems using NQS, mirroring the success of fine-tuning in machine learning and natural language processing.
△ Less
Submitted 18 March, 2024; v1 submitted 12 March, 2024;
originally announced March 2024.
-
Transformer Wave Function for the Shastry-Sutherland Model: emergence of a Spin-Liquid Phase
Authors:
Luciano Loris Viteritti,
Riccardo Rende,
Alberto Parola,
Sebastian Goldt,
Federico Becca
Abstract:
Quantum magnetism in two-dimensional systems represents a lively branch of modern condensed-matter physics. In the presence of competing super-exchange couplings, magnetic order is frustrated and can be suppressed down to zero temperature, leading to exotic ground states. The Shastry-Sutherland model, describing $S=1/2$ degrees of freedom interacting in a two-dimensional lattice, portrays a simple…
▽ More
Quantum magnetism in two-dimensional systems represents a lively branch of modern condensed-matter physics. In the presence of competing super-exchange couplings, magnetic order is frustrated and can be suppressed down to zero temperature, leading to exotic ground states. The Shastry-Sutherland model, describing $S=1/2$ degrees of freedom interacting in a two-dimensional lattice, portrays a simple example of highly-frustrated magnetism, capturing the low-temperature behavior of SrCu$_2$(BO$_3$)$_2$ with its intriguing properties. Here, we investigate this problem by using a Vision Transformer to define an extremely accurate variational wave function. From a technical side, a pivotal achievement relies on using a deep neural network with real-valued parameters, parametrized with a Transformer, to map physical spin configurations into a high-dimensional feature space. Within this abstract space, the determination of the ground-state properties is simplified, requiring only a single output layer with complex-valued parameters. From the physical side, we supply strong evidence for the stabilization of a spin-liquid between the plaquette and antiferromagnetic phases. Our findings underscore the potential of Neural-Network Quantum States as a valuable tool for probing uncharted phases of matter, opening opportunities to establish the properties of many-body systems.
△ Less
Submitted 12 February, 2024; v1 submitted 28 November, 2023;
originally announced November 2023.
-
Spin-phonon interactions on the kagome lattice: Dirac spin liquid versus valence-bond solids
Authors:
Francesco Ferrari,
Federico Becca,
Roser Valenti
Abstract:
We investigate the impact of the spin-phonon coupling on the S=1/2 Heisenberg model on the kagome lattice. For the pure spin model, there is increasing evidence that the low-energy properties can be correctly described by a Dirac spin liquid, in which spinons with a conical dispersion are coupled to emergent gauge fields. Within this scenario, the ground-state wave function is well approximated by…
▽ More
We investigate the impact of the spin-phonon coupling on the S=1/2 Heisenberg model on the kagome lattice. For the pure spin model, there is increasing evidence that the low-energy properties can be correctly described by a Dirac spin liquid, in which spinons with a conical dispersion are coupled to emergent gauge fields. Within this scenario, the ground-state wave function is well approximated by a Gutzwiller-projected fermionic state [Y. Ran, M. Hermele, P.A. Lee, and X.-G. Wen, Phys. Rev. Lett. 98, 117205 (2007)]. However, the existence of U(1) gauge fields may naturally lead to instabilities when small perturbations are included. Since phonons are ubiquitous in real materials, they may play a relevant role in the determination of the actual physical properties of the kagome antiferromagnet. We perform a step forward in this direction, including phonon degrees of freedom (at the quantum level) and applying a variational approach based upon Gutzwiller-projected fermionic Ansätze. Our results suggest that the Dirac spin liquid is stable for small spin-phonon couplings, while valence-bond solids are obtained at large couplings. Even though different distortions can be induced by the spin-phonon interaction, the general aspect is that the energy is lowered by maximizing the density of perfect hexagons in the dimerization pattern.
△ Less
Submitted 18 April, 2024; v1 submitted 16 November, 2023;
originally announced November 2023.
-
A simple linear algebra identity to optimize Large-Scale Neural Network Quantum States
Authors:
Riccardo Rende,
Luciano Loris Viteritti,
Lorenzo Bardone,
Federico Becca,
Sebastian Goldt
Abstract:
Neural-network architectures have been increasingly used to represent quantum many-body wave functions. These networks require a large number of variational parameters and are challenging to optimize using traditional methods, as gradient descent. Stochastic Reconfiguration (SR) has been effective with a limited number of parameters, but becomes impractical beyond a few thousand parameters. Here,…
▽ More
Neural-network architectures have been increasingly used to represent quantum many-body wave functions. These networks require a large number of variational parameters and are challenging to optimize using traditional methods, as gradient descent. Stochastic Reconfiguration (SR) has been effective with a limited number of parameters, but becomes impractical beyond a few thousand parameters. Here, we leverage a simple linear algebra identity to show that SR can be employed even in the deep learning scenario. We demonstrate the effectiveness of our method by optimizing a Deep Transformer architecture with $3 \times 10^5$ parameters, achieving state-of-the-art ground-state energy in the $J_1$-$J_2$ Heisenberg model at $J_2/J_1=0.5$ on the $10\times10$ square lattice, a challenging benchmark in highly-frustrated magnetism. This work marks a significant step forward in the scalability and efficiency of SR for Neural-Network Quantum States, making them a promising method to investigate unknown quantum phases of matter, where other methods struggle.
△ Less
Submitted 2 August, 2024; v1 submitted 9 October, 2023;
originally announced October 2023.
-
Piercing the Dirac spin liquid: From a single monopole to chiral states
Authors:
Sasank Budaraju,
Yasir Iqbal,
Federico Becca,
Didier Poilblanc
Abstract:
The parton approach for quantum spin liquids gives a transparent description of low-energy elementary excitations, e.g., spinons and emergent gauge-field fluctuations. The latter ones are directly coupled to the hopping/pairing of spinons. By using the fermionic representation of the $U(1)$ Dirac state on the kagome lattice and variational Monte Carlo techniques to include the Gutzwiller projectio…
▽ More
The parton approach for quantum spin liquids gives a transparent description of low-energy elementary excitations, e.g., spinons and emergent gauge-field fluctuations. The latter ones are directly coupled to the hopping/pairing of spinons. By using the fermionic representation of the $U(1)$ Dirac state on the kagome lattice and variational Monte Carlo techniques to include the Gutzwiller projection, we analyse the effect of modifying the gauge fields in the spinon kinematics. In particular, we construct low-energy monopole excitations, which are shown to be gapless in the thermodynamic limit. States with a finite number of monopoles or with a finite density of them are also considered, with different patterns of the gauge fluxes. We show that these chiral states are not stabilized in the Heisenberg model with nearest-neighbor super-exchange couplings, and the Dirac state corresponds to the lowest-energy Ansatz within this family of variational wave functions. Our results support the idea that spinons with a gapless conical spectrum coexist with gapless monopole excitations, even for the spin-1/2 case.
△ Less
Submitted 24 November, 2023; v1 submitted 3 July, 2023;
originally announced July 2023.
-
A Jastrow wave function for the spin-1 Heisenberg chain: the string order revealed by the mapping to the classical Coulomb gas
Authors:
Davide Piccioni,
Christian Apostoli,
Federico Becca,
Guglielmo Mazzola,
Alberto Parola,
Sandro Sorella,
Giuseppe E. Santoro
Abstract:
We show that a two-body Jastrow wave function is able to capture the ground-state properties of the $S=1$ antiferromagnetic Heisenberg chain with the single-ion anisotropy term, in both the topological and trivial phases. Here, the optimized Jastrow pseudo potential assumes a very simple form in Fourier space, i.e., $v_{q} \approx 1/q^2$, which is able to give rise to a finite string-order paramet…
▽ More
We show that a two-body Jastrow wave function is able to capture the ground-state properties of the $S=1$ antiferromagnetic Heisenberg chain with the single-ion anisotropy term, in both the topological and trivial phases. Here, the optimized Jastrow pseudo potential assumes a very simple form in Fourier space, i.e., $v_{q} \approx 1/q^2$, which is able to give rise to a finite string-order parameter in the topological regime. The results are analysed by using an exact mapping from the quantum expectation values over the variational state to the classical partition function of the one-dimensional Coulomb gas of particles with charge $q=\pm 1$. Here, two phases are present at low temperatures: the first one is a diluted gas of dipoles (bound states of particles with opposite charges), which are randomly oriented (describing the trivial phase); the other one is a dense liquid of dipoles, which are aligned thanks to the residual dipole-dipole interactions (describing the topological phase, with the finite string order being related to the dipole alignment). Our results provide an insightful interpretation of the ground-state nature of the spin-1 antiferromagnetic Heisenberg model.
△ Less
Submitted 23 May, 2023;
originally announced May 2023.
-
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.
-
Transformer variational wave functions for frustrated quantum spin systems
Authors:
Luciano Loris Viteritti,
Riccardo Rende,
Federico Becca
Abstract:
The Transformer architecture has become the state-of-art model for natural language processing tasks and, more recently, also for computer vision tasks, thus defining the Vision Transformer (ViT) architecture. The key feature is the ability to describe long-range correlations among the elements of the input sequences, through the so-called self-attention mechanism. Here, we propose an adaptation o…
▽ More
The Transformer architecture has become the state-of-art model for natural language processing tasks and, more recently, also for computer vision tasks, thus defining the Vision Transformer (ViT) architecture. The key feature is the ability to describe long-range correlations among the elements of the input sequences, through the so-called self-attention mechanism. Here, we propose an adaptation of the ViT architecture with complex parameters to define a new class of variational neural-network states for quantum many-body systems, the ViT wave function. We apply this idea to the one-dimensional $J_1$-$J_2$ Heisenberg model, demonstrating that a relatively simple parametrization gets excellent results for both gapped and gapless phases. In this case, excellent accuracies are obtained by a relatively shallow architecture, with a single layer of self-attention, thus largely simplifying the original architecture. Still, the optimization of a deeper structure is possible and can be used for more challenging models, most notably highly-frustrated systems in two dimensions. The success of the ViT wave function relies on mixing both local and global operations, thus enabling the study of large systems with high accuracy.
△ Less
Submitted 11 June, 2023; v1 submitted 10 November, 2022;
originally announced November 2022.
-
Static and dynamical signatures of Dzyaloshinskii-Moriya interactions in the Heisenberg model on the kagome lattice
Authors:
Francesco Ferrari,
Sen Niu,
Juraj Hasik,
Yasir Iqbal,
Didier Poilblanc,
Federico Becca
Abstract:
Motivated by recent experiments on Cs$_2$Cu$_3$SnF$_{12}$ and YCu$_{3}$(OH)$_{6}$Cl$_{3}$, we consider the ${S=1/2}$ Heisenberg model on the kagome lattice with nearest-neighbor super-exchange $J$ and (out-of-plane) Dzyaloshinskii-Moriya interaction $J_D$, which favors (in-plane) ${\bf Q}=(0,0)$ magnetic order. By using both variational Monte Carlo (based upon Gutzwiller-projected fermionic wave f…
▽ More
Motivated by recent experiments on Cs$_2$Cu$_3$SnF$_{12}$ and YCu$_{3}$(OH)$_{6}$Cl$_{3}$, we consider the ${S=1/2}$ Heisenberg model on the kagome lattice with nearest-neighbor super-exchange $J$ and (out-of-plane) Dzyaloshinskii-Moriya interaction $J_D$, which favors (in-plane) ${\bf Q}=(0,0)$ magnetic order. By using both variational Monte Carlo (based upon Gutzwiller-projected fermionic wave functions) and tensor-network approaches (built from infinite projected-entangled pair/simplex states), we show that the ground state develops a finite magnetization for $J_D/J \gtrsim 0.03 - 0.04$, while the gapless spin liquid remains stable for smaller values of the Dzyaloshinskii-Moriya interaction. The relatively small value of $J_D/J$ for which magnetic order sets in is particularly relevant for the interpretation of low-temperature behaviors of kagome antiferromagnets, including ZnCu$_{3}$(OH)$_{6}$Cl$_{2}$. In addition, we assess the spin dynamical structure factors and the corresponding low-energy spectrum, by using the variational Monte Carlo technique. The existence of a continuum of excitations above the magnon modes is reported within the magnetically ordered phase, similarly to what has been detected by inelastic neutron scattering on Cs$_{2}$Cu$_{3}$SnF$_{12}$.
△ Less
Submitted 20 February, 2023; v1 submitted 19 July, 2022;
originally announced July 2022.
-
Symmetric projected entangled-pair states analysis of a phase transition in coupled spin-1/2 ladders
Authors:
Juraj Hasik,
Glen B. Mbeng,
Sylvain Capponi,
Federico Becca,
Andreas M. Läuchli
Abstract:
Infinite projected entangled-pair states (iPEPS) have been introduced to accurately describe many-body wave functions on two-dimensional lattices. In this context, two aspects are crucial: the systematic improvement of the {\it Ansatz} by the optimization of its building blocks, i.e., tensors characterized by bond dimension $D$, and the extrapolation scheme to reach the "thermodynamic" limit…
▽ More
Infinite projected entangled-pair states (iPEPS) have been introduced to accurately describe many-body wave functions on two-dimensional lattices. In this context, two aspects are crucial: the systematic improvement of the {\it Ansatz} by the optimization of its building blocks, i.e., tensors characterized by bond dimension $D$, and the extrapolation scheme to reach the "thermodynamic" limit $D \to \infty$. Recent advances in variational optimization and scaling based on correlation lengths demonstrated the ability of iPEPS to capture the spontaneous breaking of a continuous symmetry in phases such as the antiferromagnetic (Néel) phase with high fidelity, in addition to valence-bond solids which are already well described by finite-$D$ iPEPS. In contrast, systems in the vicinity of continuous quantum phase transitions still present a challenge for iPEPS, especially when non-abelian symmetries are involved. Here, we consider the iPEPS Ansatz to describe the continuous transition between the (gapless) antiferromagnet and the (gapped) paramagnet that exists in the $S=1/2$ Heisenberg model on coupled two-leg ladders. In particular, we show how accurate iPEPS results can be obtained down to a narrow interval around criticality and analyze the scaling of the order parameter in the Néel phase in a spatially anisotropic situation.
△ Less
Submitted 13 June, 2022;
originally announced June 2022.
-
Charge-density waves in kagome-lattice extended Hubbard models at the van Hove filling
Authors:
Francesco Ferrari,
Federico Becca,
Roser Valentí
Abstract:
The Hubbard model on the kagome lattice is presently often considered as a minimal model to describe the rich low-temperature behavior of AV$_{3}$Sb$_{5}$ compounds (with A=K, Rb, Cs), including charge-density waves (CDWs), superconductivity, and possibly broken time-reversal symmetry. Here, we investigate, via variational Jastrow-Slater wave functions, the properties of its ground state when both…
▽ More
The Hubbard model on the kagome lattice is presently often considered as a minimal model to describe the rich low-temperature behavior of AV$_{3}$Sb$_{5}$ compounds (with A=K, Rb, Cs), including charge-density waves (CDWs), superconductivity, and possibly broken time-reversal symmetry. Here, we investigate, via variational Jastrow-Slater wave functions, the properties of its ground state when both onsite $U$ and nearest-neighbor $V$ Coulomb repulsions are considered at the van Hove filling. Our calculations reveal the presence of different interaction-driven CDWs and, contrary to previous renormalization-group studies, the absence of ferromagnetism and charge- or spin-bond order. No signatures of chiral phases are detected. Remarkably, the CDWs triggered by the nearest-neighbor repulsion possess charge disproportionations that are not compatible with the ones observed in AV$_{3}$Sb$_{5}$. As an alternative mechanism to stabilize charge-bond order, we consider the electron-phonon interaction, modeled by coupling the hopping amplitudes to quantum phonons, as in the Su-Schrieffer-Heeger model. Our results show the instability towards a tri-hexagonal distortion with $2\times 2$ periodicity, in a closer agreement with experimental findings.
△ Less
Submitted 15 August, 2022; v1 submitted 21 April, 2022;
originally announced April 2022.
-
A real-space many-body marker for correlated ${\mathbb Z}_2$ topological insulators
Authors:
Ivan Gilardoni,
Federico Becca,
Antimo Marrazzo,
Alberto Parola
Abstract:
Taking the clue from the modern theory of polarization [R. Resta, Rev. Mod. Phys. {\bf 66}, 899 (1994)], we identify an operator to distinguish between ${\mathbb Z}_2$-even (trivial) and ${\mathbb Z}_2$-odd (topological) insulators in two spatial dimensions. Its definition extends the position operator [R. Resta and S. Sorella, Phys. Rev. Lett. {\bf 82}, 370 (1999)], which was introduced in one-di…
▽ More
Taking the clue from the modern theory of polarization [R. Resta, Rev. Mod. Phys. {\bf 66}, 899 (1994)], we identify an operator to distinguish between ${\mathbb Z}_2$-even (trivial) and ${\mathbb Z}_2$-odd (topological) insulators in two spatial dimensions. Its definition extends the position operator [R. Resta and S. Sorella, Phys. Rev. Lett. {\bf 82}, 370 (1999)], which was introduced in one-dimensional systems. We first show a few examples of non-interacting models, where single-particle wave functions are defined and allow for a direct comparison with standard techniques on large system sizes. Then, we illustrate its applicability for an interacting model on a small cluster, where exact diagonalizations are available. Its formulation in the Fock space allows a direct computation of expectation values over the ground-state wave function (or any approximation of it), thus allowing us to investigate generic interacting systems, such as strongly-correlated topological insulators.
△ Less
Submitted 30 November, 2022; v1 submitted 24 February, 2022;
originally announced February 2022.
-
Accuracy of Restricted Boltzmann Machines for the one-dimensional $J_1-J_2$ Heisenberg model
Authors:
Luciano Loris Viteritti,
Francesco Ferrari,
Federico Becca
Abstract:
Neural networks have been recently proposed as variational wave functions for quantum many-body systems [G. Carleo and M. Troyer, Science 355, 602 (2017)]. In this work, we focus on a specific architecture, known as Restricted Boltzmann Machine (RBM), and analyse its accuracy for the spin-1/2 $J_1-J_2$ antiferromagnetic Heisenberg model in one spatial dimension. The ground state of this model has…
▽ More
Neural networks have been recently proposed as variational wave functions for quantum many-body systems [G. Carleo and M. Troyer, Science 355, 602 (2017)]. In this work, we focus on a specific architecture, known as Restricted Boltzmann Machine (RBM), and analyse its accuracy for the spin-1/2 $J_1-J_2$ antiferromagnetic Heisenberg model in one spatial dimension. The ground state of this model has a non-trivial sign structure, especially for $J_2/J_1>0.5$, forcing us to work with complex-valued RBMs. Two variational Ansätze are discussed: one defined through a fully complex RBM, and one in which two different real-valued networks are used to approximate modulus and phase of the wave function. In both cases, translational invariance is imposed by considering linear combinations of RBMs, giving access also to the lowest-energy excitations at fixed momentum $k$. We perform a systematic study on small clusters to evaluate the accuracy of these wave functions in comparison to exact results, providing evidence for the supremacy of the fully complex RBM. Our calculations show that this kind of Ansätze is very flexible and describes both gapless and gapped ground states, also capturing the incommensurate spin-spin correlations and low-energy spectrum for $J_2/J_1>0.5$. The RBM results are also compared to the ones obtained with Gutzwiller-projected fermionic states, often employed to describe quantum spin models [F. Ferrari, A. Parola, S. Sorella and F. Becca, Phys. Rev. B 97, 235103 (2018)]. Contrary to the latter class of variational states, the fully-connected structure of RBMs hampers the transferability of the wave function from small to large clusters, implying an increase of the computational cost with the system size.
△ Less
Submitted 14 April, 2022; v1 submitted 15 February, 2022;
originally announced February 2022.
-
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.
-
Gutzwiller-projected states for the $J_1$-$J_2$ Heisenberg model on the kagome lattice: achievements and pitfalls
Authors:
Yasir Iqbal,
Francesco Ferrari,
Aishwarya Chauhan,
Alberto Parola,
Didier Poilblanc,
Federico Becca
Abstract:
We assess the ground-state phase diagram of the $J_1$-$J_2$ Heisenberg model on the kagome lattice by employing Gutzwiller-projected fermionic wave functions. Within this framework, different states can be represented, defined by distinct unprojected fermionic Hamiltonians that include hopping and pairing terms, as well as a coupling to local Zeeman fields to generate magnetic order. For $J_2=0$,…
▽ More
We assess the ground-state phase diagram of the $J_1$-$J_2$ Heisenberg model on the kagome lattice by employing Gutzwiller-projected fermionic wave functions. Within this framework, different states can be represented, defined by distinct unprojected fermionic Hamiltonians that include hopping and pairing terms, as well as a coupling to local Zeeman fields to generate magnetic order. For $J_2=0$, the so-called U(1) Dirac state, in which only hopping is present (such as to generate a $π$-flux in the hexagons), has been shown to accurately describe the exact ground state [Y. Iqbal, F. Becca, S. Sorella, and D. Poilblanc, Phys. Rev. B 87, 060405 (2013); Y.-C. He, M. P. Zaletel, M. Oshikawa, and F. Pollmann, Phys. Rev. X 7, 031020 (2017)]. Here, we show that its accuracy improves in presence of a small antiferromagnetic super-exchange $J_2$, leading to a finite region where the gapless spin liquid is stable; then, for $J_2/J_1=0.11(1)$, a first-order transition to a magnetic phase with pitch vector ${\bf q}=(0,0)$ is detected, by allowing magnetic order within the fermionic Hamiltonian. Instead, for small ferromagnetic values of $|J_2|/J_1$, the situation is more contradictory. While the U(1) Dirac state remains stable against several perturbations in the fermionic part (i.e., dimerization patterns or chiral terms), its accuracy clearly deteriorates on small systems, most notably on $36$ sites where exact diagonalization is possible. Then, upon increasing the ratio $|J_2|/J_1$, a magnetically ordered state with $\sqrt{3} \times \sqrt{3}$ periodicity eventually overcomes the U(1) Dirac spin liquid. Within the ferromagnetic $J_{2}$ regime, evidence is shown in favor of a first-order transition at $J_2/J_1=-0.065(5)$.
△ Less
Submitted 19 September, 2021; v1 submitted 5 August, 2021;
originally announced August 2021.
-
Effects of spin-phonon coupling in frustrated Heisenberg models
Authors:
Francesco Ferrari,
Roser Valenti,
Federico Becca
Abstract:
The existence and stability of spin-liquid phases represent a central topic in the field of frustrated magnetism. While a few examples of spin-liquid ground states are well established in specific models (e.g. the Kitaev model on the honeycomb lattice), recent investigations have suggested the possibility of their appearance in several Heisenberg-like models on frustrated lattices. An important re…
▽ More
The existence and stability of spin-liquid phases represent a central topic in the field of frustrated magnetism. While a few examples of spin-liquid ground states are well established in specific models (e.g. the Kitaev model on the honeycomb lattice), recent investigations have suggested the possibility of their appearance in several Heisenberg-like models on frustrated lattices. An important related question concerns the stability of spin liquids in the presence of small perturbations in the Hamiltonian. In this respect, the magnetoelastic interaction between spins and phonons represents a relevant and physically motivated perturbation, which has been scarcely investigated so far. In this work, we study the effect of the spin-phonon coupling on prototypical models of frustrated magnetism. We adopt a variational framework based upon Gutzwiller-projected wave functions implemented with a spin-phonon Jastrow factor, providing a full quantum treatment of both spin and phonon degrees of freedom. The results on the frustrated $J_1-J_2$ Heisenberg model on one- and two-dimensional (square) lattices show that, while a valence-bond crystal is prone to lattice distortions, a gapless spin liquid is stable for small spin-phonon couplings. In view of the ubiquitous presence of lattice vibrations, our results are particularly important to demonstrate the possibility that gapless spin liquids may be realized in real materials.
△ Less
Submitted 14 July, 2021; v1 submitted 5 May, 2021;
originally announced May 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.
-
Gapless spin liquids in disguise
Authors:
Francesco Ferrari,
Alberto Parola,
Federico Becca
Abstract:
We show that gapless spin liquids, which are potential candidates to describe the ground state of frustrated Heisenberg models in two dimensions, become trivial insulators on cylindrical geometries with an even number of legs. In particular, we report calculations for Gutzwiller-projected fermionic states on strips of square and kagome lattices. By choosing different boundary conditions for the fe…
▽ More
We show that gapless spin liquids, which are potential candidates to describe the ground state of frustrated Heisenberg models in two dimensions, become trivial insulators on cylindrical geometries with an even number of legs. In particular, we report calculations for Gutzwiller-projected fermionic states on strips of square and kagome lattices. By choosing different boundary conditions for the fermionic degrees of freedom, both gapless and gapped states may be realized, the latter ones having a lower variational energy. The direct evaluation of static and dynamical correlation functions, as well as overlaps between different states, allows us to demonstrate the sharp difference between the ground-state properties obtained within cylinders or directly in the two-dimensional lattice. Our results shed light on the difficulty to detect bona fide gapless spin liquids in such cylindrical geometries.
△ Less
Submitted 18 May, 2021; v1 submitted 9 March, 2021;
originally announced March 2021.
-
Investigation of the Néel phase of the frustrated Heisenberg antiferromagnet by differentiable symmetric tensor networks
Authors:
Juraj Hasik,
Didier Poilblanc,
Federico Becca
Abstract:
The recent progress in the optimization of two-dimensional tensor networks [H.-J. Liao, J.-G. Liu, L. Wang, and T. Xiang, Phys. Rev. X ${\bf 9}$, 031041 (2019)] based on automatic differentiation opened the way towards precise and fast optimization of such states and, in particular, infinite projected entangled-pair states (iPEPS) that constitute a generic-purpose ${\it Ansatz}$ for lattice proble…
▽ More
The recent progress in the optimization of two-dimensional tensor networks [H.-J. Liao, J.-G. Liu, L. Wang, and T. Xiang, Phys. Rev. X ${\bf 9}$, 031041 (2019)] based on automatic differentiation opened the way towards precise and fast optimization of such states and, in particular, infinite projected entangled-pair states (iPEPS) that constitute a generic-purpose ${\it Ansatz}$ for lattice problems governed by local Hamiltonians. In this work, we perform an extensive study of a paradigmatic model of frustrated magnetism, the $J_1-J_2$ Heisenberg antiferromagnet on the square lattice. By using advances in both optimization and subsequent data analysis, through finite correlation-length scaling, we report accurate estimations of the magnetization curve in the Néel phase for $J_2/J_1 \le 0.45$. The unrestricted iPEPS simulations reveal an $U(1)$ symmetric structure, which we identify and impose on tensors, resulting in a clean and consistent picture of antiferromagnetic order vanishing at the phase transition with a quantum paramagnet at $J_2/J_1 \approx 0.46(1)$. The present methodology can be extended beyond this model to study generic order-to-disorder transitions in magnetic systems.
△ Less
Submitted 10 September, 2020; v1 submitted 4 September, 2020;
originally announced September 2020.
-
Variational wave functions for the spin-Peierls transition in the Su-Schrieffer-Heeger model with quantum phonons
Authors:
Francesco Ferrari,
Roser Valenti,
Federico Becca
Abstract:
We introduce variational wave functions to evaluate the ground-state properties of spin-phonon coupled systems described by the Su-Schrieffer-Heeger model. Quantum spins and phonons are treated on equal footing within a Monte Carlo sampling, and different regimes are investigated. We show that the proposed variational Ansatz yields good agreement with previous density-matrix renormalization group…
▽ More
We introduce variational wave functions to evaluate the ground-state properties of spin-phonon coupled systems described by the Su-Schrieffer-Heeger model. Quantum spins and phonons are treated on equal footing within a Monte Carlo sampling, and different regimes are investigated. We show that the proposed variational Ansatz yields good agreement with previous density-matrix renormalization group results in one dimension and is able to accurately describe the spin-Peierls transition. This variational approach is neither constrained by the magnetoelastic-coupling strength nor by the dimensionality of the systems considered, thus allowing future investigations in more general cases, which are relevant to spin-liquid and topological phases in two spatial dimensions.
△ Less
Submitted 3 October, 2020; v1 submitted 28 July, 2020;
originally announced July 2020.
-
Gapless spin liquid and valence-bond solid in the $J_1-J_2$ Heisenberg model on the square lattice: insights from singlet and triplet excitations
Authors:
Francesco Ferrari,
Federico Becca
Abstract:
The spin-$1/2$ $J_1-J_2$ Heisenberg model on the square lattice represents one of the simplest examples in which the effets of magnetic interactions may suppress magnetic order, eventually leading to a pure quantum phase with no local order parameters. This model has been extensively studied in the last three decades, with conflicting results. Here, by using Gutzwiller-projected wave functions and…
▽ More
The spin-$1/2$ $J_1-J_2$ Heisenberg model on the square lattice represents one of the simplest examples in which the effets of magnetic interactions may suppress magnetic order, eventually leading to a pure quantum phase with no local order parameters. This model has been extensively studied in the last three decades, with conflicting results. Here, by using Gutzwiller-projected wave functions and recently developed methods to assess the low-energy spectrum, we show the existence of a level crossing between the lowest-energy triplet and singlet excitations for $J_2/J_1 \approx 0.54$. This fact supports the existence of a phase transition between a gapless spin liquid (which is stable for $0.48 \lesssim J_2/J_1 \lesssim 0.54$) and a valence-bond solid (for $0.54 \lesssim J_2/J_1 \lesssim 0.6$), even though no clear sign of dimer order is visible in the correlations functions. These results, which confirm recent density-matrix renormalization calculations on cylindrical clusters [L. Wang and A.W. Sandvik, Phys. Rev. Lett. 121, 107202 (2018)] reconcile the contraddicting results obtained within different approaches over the years.
△ Less
Submitted 26 May, 2020;
originally announced May 2020.
-
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.
-
Dynamical properties of Néel and valence-bond phases in the $J_1-J_2$ model on the honeycomb lattice
Authors:
Francesco Ferrari,
Federico Becca
Abstract:
By using a variational Monte Carlo technique based upon Gutzwiller-projected fermionic states, we investigate the dynamical structure factor of the antiferromagnetic $S=1/2$ Heisenberg model on the honeycomb lattice, in presence of first-neighbor ($J_1$) and second-neighbor ($J_2$) couplings, for ${J_2 < 0.5 J_1}$. The ground state of the system shows long-range antiferromagnetic order for…
▽ More
By using a variational Monte Carlo technique based upon Gutzwiller-projected fermionic states, we investigate the dynamical structure factor of the antiferromagnetic $S=1/2$ Heisenberg model on the honeycomb lattice, in presence of first-neighbor ($J_1$) and second-neighbor ($J_2$) couplings, for ${J_2 < 0.5 J_1}$. The ground state of the system shows long-range antiferromagnetic order for ${J_2/J_1 \lesssim 0.23}$, plaquette valence-bond order for ${0.23 \lesssim J_2/J_1 \lesssim 0.36}$, and columnar dimer order for ${J_2/J_1 \gtrsim 0.36}$. Within the antiferromagnetic state, a well-defined magnon mode is observed, whose dispersion is in relatively good agreement with linear spin-wave approximation for $J_2=0$. When a nonzero second-neighbor super-exchange is included, a roton-like mode develops around the $K$ point (i.e., the corner of the Brillouin zone). This mode softens when $J_2/J_1$ is increased and becomes gapless at the transition point, $J_2/J_1 \approx 0.23$. Here, a broad continuum of states is clearly visible in the dynamical spectrum, suggesting that nearly-deconfined spinon excitations could exist, at least at relatively high energies. For larger values of $J_2/J_1$, valence-bond order is detected and the spectrum of the system becomes clearly gapped, with a triplon mode at low energies. This is particularly evident for the spectrum of the dimer valence-bond phase, in which the triplon mode is rather well separated from the continuum of excitations that appears at higher energies.
△ Less
Submitted 19 December, 2019;
originally announced December 2019.
-
Neural Gutzwiller-projected variational wave functions
Authors:
Francesco Ferrari,
Federico Becca,
Juan Carrasquilla
Abstract:
Variational wave functions have enabled exceptional scientific breakthroughs related to the understanding of novel phases of matter. Examples include the Bardeen-Cooper-Schrieffer theory of superconductivity, the description of the fractional quantum Hall effect through the Laughlin state, and Feynman's variational understanding of large-scale quantum effects in liquid Helium. More recently, Gutzw…
▽ More
Variational wave functions have enabled exceptional scientific breakthroughs related to the understanding of novel phases of matter. Examples include the Bardeen-Cooper-Schrieffer theory of superconductivity, the description of the fractional quantum Hall effect through the Laughlin state, and Feynman's variational understanding of large-scale quantum effects in liquid Helium. More recently, Gutzwiller-projected wave functions, typically constructed from fermionic degrees of freedom, have been employed to examine quantum spin models in the presence of competing interactions, where exotic phases with no spontaneous symmetry breaking and fractional excitations may exist. In this work, we investigate the aforementioned fermionic wave functions supplemented with neural networks, specifically with the so-called restricted Boltzmann machine (RBM), to boost their accuracy and obtain reliable approximations to the ground state of generic spin models. In particular, we apply our neural augmented fermionic construction to the description of both magnetically ordered and disordered phases of increasing complexity, including cases where the ground state displays a non-trivial sign structure. Even though the RBM state is by far more effective for Néel states endowed with a particularly simple sign structure, it provides a significant improvement over the original fermionic state in highly frustrated regimes where a complex sign structure is anticipated, thus marking the path to an understanding of strongly-correlated spin models on the lattice via neural Gutzwiller-projected variational wave functions.
△ Less
Submitted 13 September, 2019; v1 submitted 2 June, 2019;
originally announced June 2019.
-
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.
-
Optimization of infinite projected entangled pair states: The role of multiplets and their breaking
Authors:
Juraj Hasik,
Federico Becca
Abstract:
The infinite projected entangled pair states (iPEPS) technique [J. Jordan {\it et al.}, Phys. Rev. Lett. {\bf 101}, 250602 (2008)] has been widely used in the recent years to assess the properties of two-dimensional quantum systems, working directly in the thermodynamic limit. This formalism, which is based upon a tensor-network representation of the ground-state wave function, has several appeali…
▽ More
The infinite projected entangled pair states (iPEPS) technique [J. Jordan {\it et al.}, Phys. Rev. Lett. {\bf 101}, 250602 (2008)] has been widely used in the recent years to assess the properties of two-dimensional quantum systems, working directly in the thermodynamic limit. This formalism, which is based upon a tensor-network representation of the ground-state wave function, has several appealing features, e.g., encoding the so-called area law of entanglement entropy by construction; still, the method presents critical issues when dealing with the optimization of tensors, in order to find the best possible approximation to the exact ground state of a given Hamiltonian. Here, we discuss the obstacles that arise in the optimization by imaginary-time evolution within the so-called simple and full updates and connect them to the emergence of a sharp multiplet structure in the "virtual" indices of tensors. In this case, a generic choice of the bond dimension $D$ is not compatible with the multiplets and leads to a symmetry breaking (e.g., generating a finite magnetic order). In addition, varying the initial guess, different final states may be reached, with very large deviations in the magnetization value. In order to exemplify this behavior, we show the results of the $S=1/2$ Heisenberg model on an array of coupled ladders, for which a vanishing magnetization below the critical interladder coupling is recovered only for selected values of $D$, while a blind optimization with a generic $D$ gives rise to a finite magnetization down to the limit of decoupled ladders.
△ Less
Submitted 20 August, 2019; v1 submitted 6 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.
-
Dynamical structure factor of the $J_1-J_2$ Heisenberg model on the triangular lattice: magnons, spinons, and gauge fields
Authors:
Francesco Ferrari,
Federico Becca
Abstract:
Understanding the nature of the excitation spectrum in quantum spin liquids is of fundamental importance, in particular for the experimental detection of candidate materials. However, current theoretical and numerical techniques have limited capabilities, especially in obtaining the dynamical structure factor, which gives a crucial characterization of the ultimate nature of the quantum state and m…
▽ More
Understanding the nature of the excitation spectrum in quantum spin liquids is of fundamental importance, in particular for the experimental detection of candidate materials. However, current theoretical and numerical techniques have limited capabilities, especially in obtaining the dynamical structure factor, which gives a crucial characterization of the ultimate nature of the quantum state and may be directly assessed by inelastic neutron scattering. In this work, we investigate the low-energy properties of the $S=1/2$ Heisenberg model on the triangular lattice, including both nearest-neighbor $J_1$ and next-nearest-neighbor $J_2$ super-exchanges, by a dynamical variational Monte Carlo approach that allows accurate results on spin models. For $J_2=0$, our calculations are compatible with the existence of a well-defined magnon in the whole Brillouin zone, with gapless excitations at $K$ points (i.e., at the corners of the Brillouin zone). The strong renormalization of the magnon branch (also including roton-like minima around the $M$ points, i.e., midpoints of the border zone) is described by our Gutzwiller-projected state, where Abrikosov fermions are subject to a non-trivial magnetic $π$-flux threading half of the triangular plaquettes. When increasing the frustrating ratio $J_2/J_1$, we detect a progessive softening of the magnon branch at $M$, which eventually becomes gapless within the spin-liquid phase. This feature is captured by the band structure of the unprojected wave function (with $2$ Dirac points for each spin component). In addition, we observe an intense signal at low energies around the $K$ points, which cannot be understood within the unprojected picture and emerges only when the Gutzwiller projection is considered, suggesting the relevance of gauge fields for the low-energy physics of spin liquids.
△ Less
Submitted 19 July, 2019; v1 submitted 13 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.
-
Spectral signatures of fractionalization in the frustrated Heisenberg model on the square lattice
Authors:
Francesco Ferrari,
Federico Becca
Abstract:
We employ a variational Monte Carlo approach to efficiently obtain the dynamical structure factor for the spin-1/2 $J_1-J_2$ Heisenberg model on the square lattice. Upon increasing the frustrating ratio $J_2/J_1$, the ground state undergoes a continuous transition from a Néel antiferromagnet to a $\mathbb{Z}_{2}$ gapless spin liquid. We identify the characteristic spectral features in both phases…
▽ More
We employ a variational Monte Carlo approach to efficiently obtain the dynamical structure factor for the spin-1/2 $J_1-J_2$ Heisenberg model on the square lattice. Upon increasing the frustrating ratio $J_2/J_1$, the ground state undergoes a continuous transition from a Néel antiferromagnet to a $\mathbb{Z}_{2}$ gapless spin liquid. We identify the characteristic spectral features in both phases and highlight the existence of a broad continuum of excitations in the proximity of the spin-liquid phase. The magnon branch, which dominates the spectrum of the unfrustrated Heisenberg model, gradually loses its spectral weight, thus releasing nearly-deconfined spinons, whose signatures are visible even in the magnetically ordered state. Our results show how free spinons emerge across a quantum critical point, providing evidence for the fractionalization of magnons into deconfined spinons.
△ Less
Submitted 17 September, 2018; v1 submitted 23 May, 2018;
originally announced May 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.
-
Dynamical structure factor of the $J_1-J_2$ Heisenberg model in one dimension: the variational Monte Carlo approach
Authors:
Francesco Ferrari,
Alberto Parola,
Sandro Sorella,
Federico Becca
Abstract:
The dynamical spin structure factor is computed within a variational framework to study the one-dimensional $J_1-J_2$ Heisenberg model. Starting from Gutzwiller-projected fermionic wave functions, the low-energy spectrum is constructed from two-spinon excitations. The direct comparison with Lanczos calculations on small clusters demonstrates the excellent description of both gapless and gapped (di…
▽ More
The dynamical spin structure factor is computed within a variational framework to study the one-dimensional $J_1-J_2$ Heisenberg model. Starting from Gutzwiller-projected fermionic wave functions, the low-energy spectrum is constructed from two-spinon excitations. The direct comparison with Lanczos calculations on small clusters demonstrates the excellent description of both gapless and gapped (dimerized) phases, also including incommensurate structures for $J_2/J_1>0.5$. Calculations on large clusters show how the intensity evolves when increasing the frustrating ratio and give an unprecedented accurate characterization of the dynamical properties of (non-integrable) frustrated spin models.
△ Less
Submitted 7 June, 2018; v1 submitted 6 March, 2018;
originally announced March 2018.
-
Persistence of the gapless spin liquid in the breathing kagome Heisenberg antiferromagnet
Authors:
Yasir Iqbal,
Didier Poilblanc,
Ronny Thomale,
Federico Becca
Abstract:
The nature of the ground state of the spin $S=1/2$ Heisenberg antiferromagnet on the kagome lattice with breathing anisotropy (i.e., with different superexchange couplings $J_{\vartriangle}$ and $J_{\triangledown}$ within elementary up- and down-pointing triangles) is investigated within the framework of Gutzwiller projected fermionic wave functions and Monte Carlo methods. We analyze the stabilit…
▽ More
The nature of the ground state of the spin $S=1/2$ Heisenberg antiferromagnet on the kagome lattice with breathing anisotropy (i.e., with different superexchange couplings $J_{\vartriangle}$ and $J_{\triangledown}$ within elementary up- and down-pointing triangles) is investigated within the framework of Gutzwiller projected fermionic wave functions and Monte Carlo methods. We analyze the stability of the U(1) Dirac spin liquid with respect to the presence of fermionic pairing that leads to a gapped $\mathbb{Z}_{2}$ spin liquid. For several values of the ratio $J_{\triangledown}/J_{\vartriangle}$, the size scaling of the energy gain due to the pairing fields and the variational parameters are reported. Our results show that the energy gain of the gapped spin liquid with respect to the gapless state either vanishes for large enough system size or scales to zero in the thermodynamic limit. Similarly, the optimized pairing amplitudes (responsible for opening the spin gap) are shown to vanish in the thermodynamic limit. Our outcome is corroborated by the application of one and two Lanczos steps to the gapless and gapped wave functions, for which no energy gain of the gapped state is detected when improving the quality of the variational states. Finally, we discuss the competition with the "simplex" $\mathbb{Z}_{2}$ resonating-valence-bond spin liquid, valence-bond crystal, and nematic states in the strongly anisotropic regime, i.e., $J_{\triangledown} \ll J_{\vartriangle}$.
△ Less
Submitted 15 March, 2018; v1 submitted 12 December, 2017;
originally announced December 2017.
-
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.
-
Competition between spin liquids and valence-bond order in the frustrated spin-$1/2$ Heisenberg model on the honeycomb lattice
Authors:
Francesco Ferrari,
Samuel Bieri,
Federico Becca
Abstract:
Using variational wave functions and Monte Carlo techniques, we study the antiferromagnetic Heisenberg model with first-neighbor $J_1$ and second-neighbor $J_2$ antiferromagnetic couplings on the honeycomb lattice. We perform a systematic comparison of magnetically ordered and nonmagnetic states (spin liquids and valence-bond solids) to obtain the ground-state phase diagram. Néel order is stabiliz…
▽ More
Using variational wave functions and Monte Carlo techniques, we study the antiferromagnetic Heisenberg model with first-neighbor $J_1$ and second-neighbor $J_2$ antiferromagnetic couplings on the honeycomb lattice. We perform a systematic comparison of magnetically ordered and nonmagnetic states (spin liquids and valence-bond solids) to obtain the ground-state phase diagram. Néel order is stabilized for small values of the frustrating second-neighbor coupling. Increasing the ratio $J_2/J_1$, we find strong evidence for a continuous transition to a nonmagnetic phase at $J_2/J_1 \approx 0.23$. Close to the transition point, the Gutzwiller-projected uniform resonating valence bond state gives an excellent approximation to the exact ground-state energy. For $0.23 \lesssim J_2/J_1 \lesssim 0.4$, a gapless $Z_2$ spin liquid with Dirac nodes competes with a plaquette valence-bond solid. In contrast, the gapped spin liquid considered in previous works has significantly higher variational energy. Although the plaquette valence-bond order is expected to be present as soon as the Néel order melts, this ordered state becomes clearly favored only for $J_2/J_1 \gtrsim 0.3$. Finally, for $0.36 \lesssim J_2/J_1 \le 0.5$, a valence-bond solid with columnar order takes over as the ground state, being also lower in energy than the magnetic state with collinear order. We perform a detailed finite-size scaling and standard data collapse analysis, and we discuss the possibility of a deconfined quantum critical point separating the Néel antiferromagnet from the plaquette valence-bond solid.
△ Less
Submitted 5 October, 2017; v1 submitted 23 June, 2017;
originally announced June 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.
-
Non local parity order in the two-dimensional Mott insulator
Authors:
Serena Fazzini,
Federico Becca,
Arianna Montorsi
Abstract:
The Mott insulator is characterized by having small deviations around the (integer) average particle density n, with pairs with n-1 and n+1 particles forming bound states. In one dimension, the effect is captured by a non-zero value of a non-local "string" of parities, which instead vanishes in the superfluid phase where density fluctuations are large. Here, we investigate the interaction induced…
▽ More
The Mott insulator is characterized by having small deviations around the (integer) average particle density n, with pairs with n-1 and n+1 particles forming bound states. In one dimension, the effect is captured by a non-zero value of a non-local "string" of parities, which instead vanishes in the superfluid phase where density fluctuations are large. Here, we investigate the interaction induced transition from the superfluid to the Mott insulator, in the paradigmatic Bose Hubbard model at n=1. By means of quantum Monte Carlo simulations and finite size scaling analysis on LxM ladders, we explore the behavior of "brane" parity operators for L going to infinity from one dimension (i.e., M=1) to two dimensions (i.e., M going to infinity). We confirm the conjecture that, adopting a standard definition, their average value decays to zero in two dimensions also in the insulating phase, evaluating the scaling factor of the "perimeter law" [S.P. Rath et al., Ann. Phys. (N.Y.) 334, 256 (2013)]. Upon introducing a further phase in the brane parity, we show that its expectation value becomes non-zero in the insulator, while still vanishing at the transition to the superfluid phase. These quantities are directly accessible to experimental measures, thus providing an insightful signature of the Mott insulator.
△ Less
Submitted 2 May, 2017; v1 submitted 31 December, 2016;
originally announced January 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.
-
Comment on "$Z_{2}$ spin liquid phase on the kagome lattice: a new saddle point", by Tao Li [arXiv:1601.02165 (2016)]
Authors:
Yasir Iqbal,
Didier Poilblanc,
Federico Becca
Abstract:
In a recent paper [arXiv:1601.02165], Tao Li claimed that a gapped $\mathbb{Z}_2$ spin liquid, obtained using projected Gutzwiller fermionic wave functions, can be stabilized in the Heisenberg model on the kagome lattice. We perform very accurate variational calculations and confirm with unprecedented accuracy the fact, that the $\mathbb{Z}_{2}$ spin liquid is a local energy minimum that goes away…
▽ More
In a recent paper [arXiv:1601.02165], Tao Li claimed that a gapped $\mathbb{Z}_2$ spin liquid, obtained using projected Gutzwiller fermionic wave functions, can be stabilized in the Heisenberg model on the kagome lattice. We perform very accurate variational calculations and confirm with unprecedented accuracy the fact, that the $\mathbb{Z}_{2}$ spin liquid is a local energy minimum that goes away with system size, thus reaffirming the scenario of a gapless $U(1)$ Dirac spin liquid ground state.
△ Less
Submitted 5 June, 2016;
originally announced June 2016.
-
Spin liquid nature in the Heisenberg $J_{1}$-$J_{2}$ triangular antiferromagnet
Authors:
Yasir Iqbal,
Wen-Jun Hu,
Ronny Thomale,
Didier Poilblanc,
Federico Becca
Abstract:
We investigate the spin-$\frac{1}{2}$ Heisenberg model on the triangular lattice in the presence of nearest-neighbor $J_1$ and next-nearest-neighbor $J_2$ antiferromagnetic couplings. Motivated by recent findings from density-matrix renormalization group (DMRG) claiming the existence of a gapped spin liquid with signatures of spontaneously broken lattice point group symmetry [Zhu and White, Phys.…
▽ More
We investigate the spin-$\frac{1}{2}$ Heisenberg model on the triangular lattice in the presence of nearest-neighbor $J_1$ and next-nearest-neighbor $J_2$ antiferromagnetic couplings. Motivated by recent findings from density-matrix renormalization group (DMRG) claiming the existence of a gapped spin liquid with signatures of spontaneously broken lattice point group symmetry [Zhu and White, Phys. Rev. B 92, 041105 (2015); Hu, Gong, Zhu, and Sheng, Phys. Rev. B 92, 140403 (2015)], we employ the variational Monte Carlo (VMC) approach to analyze the model from an alternative perspective that considers both magnetically ordered and paramagnetic trial states. We find a quantum paramagnet in the regime $0.08\lesssim J_2/J_1\lesssim 0.16$, framed by $120^{\circ}$ coplanar (stripe collinear) antiferromagnetic order for smaller (larger) $J_2/J_1$. By considering the optimization of spin-liquid wave functions of a different gauge group and lattice point group content as derived from Abrikosov mean-field theory, we obtain the gapless $U(1)$ Dirac spin liquid as the energetically most preferable state in comparison to all symmetric or nematic gapped $\mathbb{Z}_{2}$ spin liquids so far advocated by DMRG. Moreover, by the application of few Lanczos iterations, we find the energy to be the same as the DMRG result within error-bars. To further resolve the intriguing disagreement between VMC and DMRG, we complement our methodological approach by the pseudofermion functional renormalization group (PFFRG) to compare the spin structure factors for the paramagnetic regime calculated by VMC, DMRG, and PFFRG. This model promises to be an ideal test-bed for future numerical refinements in tracking the long-range correlations in frustrated magnets.
△ Less
Submitted 11 April, 2016; v1 submitted 22 January, 2016;
originally announced January 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.
-
Variational Monte Carlo study of gapless spin liquid in the spin-$1/2$ XXZ antiferromagnetic model on the kagome lattice
Authors:
Wen-Jun Hu,
Shou-Shu Gong,
Federico Becca,
D. N. Sheng
Abstract:
By using the variational Monte Carlo technique, we study the spin-$1/2$ XXZ antiferromagnetic model (with easy-plane anisotropy) on the kagome lattice. A class of Gutzwiller projected fermionic states with a spin Jastrow factor is considered to describe either spin liquids (with $U(1)$ or $Z_2$ symmetry) or magnetically ordered phases (with ${\bf q}=(0,0)$ or ${\bf q}=(4π/3,0)$). We find that the…
▽ More
By using the variational Monte Carlo technique, we study the spin-$1/2$ XXZ antiferromagnetic model (with easy-plane anisotropy) on the kagome lattice. A class of Gutzwiller projected fermionic states with a spin Jastrow factor is considered to describe either spin liquids (with $U(1)$ or $Z_2$ symmetry) or magnetically ordered phases (with ${\bf q}=(0,0)$ or ${\bf q}=(4π/3,0)$). We find that the magnetic states are not stable in the thermodynamic limit. Moreover, there is no energy gain to break the gauge symmetry from $U(1)$ to $Z_2$ within the spin-liquid states, as previously found in the Heisenberg model. The best variational wave function is therefore the $U(1)$ Dirac state, supplemented by the spin Jastrow factor. Furthermore, a vanishing $S=2$ spin gap is obtained at the variational level, in the whole regime from the $XY$ to the Heisenberg model.
△ Less
Submitted 17 November, 2015; v1 submitted 23 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.
-
Lanczos steps to improve variational wave functions
Authors:
Federico Becca,
Wen-Jun Hu,
Yasir Iqbal,
Alberto Parola,
Didier Poilblanc,
Sandro Sorella
Abstract:
Gutzwiller-projected fermionic states can be efficiently implemented within quantum Monte Carlo calculations to define extremely accurate variational wave functions for Heisenberg models on frustrated two-dimensional lattices, not only for the ground state but also for low-energy excitations. The application of few Lanczos steps on top of these states further improves their accuracy, allowing calc…
▽ More
Gutzwiller-projected fermionic states can be efficiently implemented within quantum Monte Carlo calculations to define extremely accurate variational wave functions for Heisenberg models on frustrated two-dimensional lattices, not only for the ground state but also for low-energy excitations. The application of few Lanczos steps on top of these states further improves their accuracy, allowing calculations on large clusters. In addition, by computing both the energy and its variance, it is possible to obtain reliable estimations of exact results. Here, we report the cases of the frustrated Heisenberg models on square and Kagome lattices.
△ Less
Submitted 8 December, 2014;
originally announced December 2014.
-
Variational Monte Carlo study of chiral spin liquid in the extended Heisenberg model on the Kagome lattice
Authors:
Wen-Jun Hu,
Wei Zhu,
Yi Zhang,
Shoushu Gong,
Federico Becca,
D. N. Sheng
Abstract:
We investigate the extended Heisenberg model on the Kagome lattice by using Gutzwiller projected fermionic states and the variational Monte Carlo technique. In particular, when both second- and third-neighbor super-exchanges are considered, we find that a gapped spin liquid described by non-trivial magnetic fluxes and long-range chiral-chiral correlations is energetically favored compared to the g…
▽ More
We investigate the extended Heisenberg model on the Kagome lattice by using Gutzwiller projected fermionic states and the variational Monte Carlo technique. In particular, when both second- and third-neighbor super-exchanges are considered, we find that a gapped spin liquid described by non-trivial magnetic fluxes and long-range chiral-chiral correlations is energetically favored compared to the gapless U(1) Dirac state. Furthermore, the topological Chern number, obtained by integrating the Berry curvature, and the degeneracy of the ground state, by constructing linearly independent states, lead us to identify this flux state as the chiral spin liquid with $C=1/2$ fractionalized Chern number.
△ Less
Submitted 22 May, 2015; v1 submitted 5 November, 2014;
originally announced November 2014.
-
Spin-$\frac{1}{2}$ Heisenberg $J_1$-$J_2$ antiferromagnet on the kagome lattice
Authors:
Yasir Iqbal,
Didier Poilblanc,
Federico Becca
Abstract:
We report variational Monte Carlo calculations for the spin-$\frac{1}{2}$ Heisenberg model on the kagome lattice in the presence of both nearest-neighbor $J_1$ and next-nearest-neighbor $J_2$ antiferromagnetic superexchange couplings. Our approach is based upon Gutzwiller projected fermionic states that represent a flexible tool to describe quantum spin liquids with different properties (e.g., gap…
▽ More
We report variational Monte Carlo calculations for the spin-$\frac{1}{2}$ Heisenberg model on the kagome lattice in the presence of both nearest-neighbor $J_1$ and next-nearest-neighbor $J_2$ antiferromagnetic superexchange couplings. Our approach is based upon Gutzwiller projected fermionic states that represent a flexible tool to describe quantum spin liquids with different properties (e.g., gapless and gapped). We show that, on finite clusters, a gapped $\mathbb{Z}_{2}$ spin liquid can be stabilized in the presence of a finite $J_2$ superexchange, with a substantial energy gain with respect to the gapless $U(1)$ Dirac spin liquid. However, this energy gain vanishes in the thermodynamic limit, implying that, at least within this approach, the $U(1)$ Dirac spin liquid remains stable in a relatively large region of the phase diagram. For $J_2/J_1 \gtrsim 0.3$, we find that a magnetically ordered state with ${\bf q}={\bf 0}$ overcomes the magnetically disordered wave functions, suggesting the end of the putative gapless spin-liquid phase.
△ Less
Submitted 9 January, 2015; v1 submitted 27 October, 2014;
originally announced October 2014.