BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Federico Becca

Federico Becca appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.str-elcond-mat.quant-gascond-mat.supr-conphysics.comp-phcond-mat.othercond-mat.stat-mechquant-ph

Catalog footprint

What is connected

46works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

46 published item(s)

preprint2022arXiv

Accuracy of Restricted Boltzmann Machines for the one-dimensional $J_1-J_2$ Heisenberg model

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Charge-density waves in kagome-lattice extended Hubbard models at the van Hove filling

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Stripes in the extended $t-t^\prime$ Hubbard model: A Variational Monte Carlo analysis

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$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Investigation of the Néel phase of the frustrated Heisenberg antiferromagnet by differentiable symmetric tensor networks

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Dynamical properties of Néel and valence-bond phases in the $J_1-J_2$ model on the honeycomb lattice

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Assessing the orbital selective Mott transition with variational wave functions

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$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Comment on "$Z_{2}$ spin liquid phase on the kagome lattice: a new saddle point", by Tao Li [arXiv:1601.02165 (2016)]

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Hidden Mott transition and large-$U$ superconductivity in the two-dimensional Hubbard model

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Spin liquid nature in the Heisenberg $J_{1}$-$J_{2}$ triangular antiferromagnet

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Spontaneous symmetry breaking in correlated wave functions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Variational wave functions for the $S=1/2$ Heisenberg model on the anisotropic triangular lattice: Spin liquids and spiral orders

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Solutions of the Two Dimensional Hubbard Model: Benchmarks and Results from a Wide Range of Numerical Algorithms

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Spin-$\frac{1}{2}$ Heisenberg $J_1$-$J_2$ antiferromagnet on the kagome lattice

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Variational Monte Carlo study of chiral spin liquid in the extended Heisenberg model on the Kagome lattice

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Variational Monte Carlo study of gapless spin liquid in the spin-$1/2$ XXZ antiferromagnetic model on the kagome lattice

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Lanczos steps to improve variational wave functions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Light-cone effect and supersonic correlations in one- and two-dimensional bosonic superfluids

We study the spreading of density-density correlations in Bose-Hubbard models after a quench of the interaction strength, using time-dependent variational Monte Carlo simulations. It gives access to unprecedented long propagation times and to dimensions higher than one. In both one and two dimensions, we find ballistic light-cone spreading of correlations and extract accurate values of the light-cone velocity in the superfluid regime. We show that the spreading of correlations is generally supersonic, with a light-cone propagating faster than sound modes but slower than the maximum group velocity of density excitations, except at the Mott transition, where all the characteristic velocities are equal. Further, we show that in two dimensions the correlation spreading is highly anisotropic and presents nontrivial interference effects.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nature of the phases in the frustrated XY model on the honeycomb lattice

We study the phase diagram of the frustrated XY model on the honeycomb lattice by using accurate correlated wave functions and variational Monte Carlo simulations. Our results suggest that a spin-liquid state is energetically favorable in the region of intermediate frustration, intervening between two magnetically ordered phases. The latter ones are represented by classically ordered states supplemented with a long-range Jastrow factor, which includes relevant correlations and dramatically improves the description provided by the purely classical solution of the model. The construction of the spin-liquid state is based on a decomposition of the underlying bosonic particles in terms of spin-1/2 fermions (partons), with a Gutzwiller projection enforcing no single occupancy, as well as a long-range Jastrow factor.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

One-dimensional spin liquid, collinear, and spiral phases from uncoupled chains to the triangular lattice

We investigate the Hubbard model on the anisotropic triangular lattice with two hopping parameters $t$ and $t^\prime$ in different spatial directions, interpolating between decoupled chains ($t=0$) and the isotropic triangular lattice ($t=t^\prime$). Variational wave functions that include both Jastrow and backflow terms are used to compare spin-liquid and magnetic phases with different pitch vectors describing both collinear and coplanar (spiral) order. For relatively large values of the on-site interaction $U/t^\prime \gtrsim 10$ and substantial frustration, i.e., $0.3\lesssim t/t^\prime \lesssim 0.8$, the spin-liquid state is clearly favored over magnetic states. Spiral magnetic order is only stable in the vicinity of the isotropic point, while collinear order is obtained in a wide range of inter-chain hoppings from small to intermediate frustration.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Spiral antiferromagnets beyond the spin-wave approximation: frustrated $XY$ and Heisenberg models in the honeycomb lattice

We examine the stability of classical states with a generic incommensurate spiral order against quantum fluctuations. Specifically, we focus on the frustrated spin-1/2 $XY$ and Heisenberg models on the honeycomb lattice with nearest-neighbor $J_1$ and next-nearest-neighbor $J_2$ antiferromagnetic couplings. Our variational approach is based on the Jastrow wave functions, which include quantum correlations on top of classical spin waves. We perform a systematic optimization of wave vectors and Jastrow pseudo-potentials within this class of variational states and find that quantum fluctuations favor collinear states over generic coplanar spirals. The Néel state with ${\bf Q}=(0,0)$ extends its stability well beyond the classical value $J_2/J_1=1/6$. Most importantly, the collinear states with ${\bf Q}=(0,2π/\sqrt{3})$ (and the two symmetry-related states) are found to be stable in a large regime with intermediate frustration, while at the classical level they are limited to the point $J_2/J_1=0.5$. For large frustration, the $120^\circ$ state is stabilized for finite values of $J_2/J_1$ in both models.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Superconductivity from spoiling magnetism in the Kondo lattice model

We find evidence that superconductivity intrudes into the paramagnetic-to-magnetic transition of the Kondo lattice model if magnetic frustration is added. Specifically, we study by variational method the model on a square lattice in the presence of both nearest- ($t$) and next-nearest-neighbor ($t^\prime$) hopping of the conduction electrons. We find that when $t^\prime/t$ is finite a $d$-wave superconducting dome emerges between the magnetic and paramagnetic metal phases and close to the compensated regime, i.e., the number of conduction electrons equal to the number of localized spin-1/2 moments. Superconductivity is further strengthened by a direct antiferromagnetic exchange between the localized moments, to such an extent that we even observe coexistence with magnetic order.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Vanishing spin gap in a competing spin-liquid phase in the kagome Heisenberg antiferromagnet

We provide strong numerical evidence, using improved variational wave functions, for a ground state with vanishing spin gap in the spin-$1/2$ quantum Heisenberg model on the kagome lattice. Starting from the algebraic $U(1)$ Dirac spin liquid state proposed by Y. Ran $et al.$ [Phys. Rev. Lett. $98$, $117205$ ($2007$)] and iteratively applying a few Lanczos steps, we compute the lowest $S=2$ excitation constructed by exciting spinons close to the Dirac nodes. Our results are compatible with a vanishing spin gap in the thermodynamic limit and in consonance with a power-law decay of long distance spin-spin correlations in real space. The competition with a gapped (topological) spin liquid is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Direct evidence for a gapless $Z_2$ spin liquid by frustrating Néel antiferromagnetism

By direct calculations of the spin gap in the frustrated Heisenberg model on the square lattice, with nearest- ($J_1$) and next-nearest-neighbor ($J_2$) super-exchange couplings, we provide a solid evidence that the spin-liquid phase in the frustrated regime $0.45 \lesssim J_2/J_1 \lesssim 0.6$ is gapless. Our numerical method is based on a variational wave function that is {\it systematically} improved by the application of few Lanczos steps and allows us to obtain reliable extrapolations in the thermodynamic limit. The peculiar nature of the non-magnetic state is unveiled by the existence of S=1 gapless excitations at $k=(π,0)$ and $(0,π)$. The magnetic transition can be described and interpreted by a variational state that is built from Abrikosov fermions having a $Z_2$ gauge structure and four Dirac points in the spinon spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Gapless spin-liquid phase in the kagome spin-1/2 Heisenberg antiferromagnet

We study the energy and the static spin structure factor of the ground state of the spin-1/2 quantum Heisenberg antiferromagnetic model on the kagome lattice. By the iterative application of a few Lanczos steps on accurate projected fermionic wave functions and the Green's function Monte Carlo technique, we find that a gapless (algebraic) U(1) Dirac spin liquid is competitive with previously proposed gapped (topological) Z2 spin liquids. By performing a finite-size extrapolation of the ground-state energy, we obtain an energy per site E/J=-0.4365(2), which is equal, within three error bars, to the estimates given by the density-matrix renormalization group (DMRG). Our estimate is obtained for a translationally invariant system, and, therefore, does not suffer from boundary effects, like in DMRG. Moreover, on finite toric clusters at the pure variational level, our energies are lower compared to those from DMRG calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Quantum quenches in one-dimensional gapless systems: Does bosonization work?

We present a comparison between the bosonization results for quantum quenches and exact diagonalizations in microscopic models of interacting spinless fermions in a one-dimensional lattice. We show that important features are missed by the bosonization technique, which predicts the persistence of long-wavelength critical properties in the long-time evolution. Instead, numerical analysis provides puzzling evidences: while the momentum distribution appears to be consistent with the presence of a singularity at $k_F$, density-density correlations at small momenta clearly display a thermal-like behavior, namely ${\bar{N}}_q \simeq {\rm const}$ (where the overbar indicates the long-time average). This feature at small momenta is preserved in presence of an interaction term that breaks integrability, together with a rounding of the singularities at finite $q$'s, showing that the bosonization approach is not able to represent the time evolution of generic one-dimensional models after a quantum quench.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Spin-liquid versus spiral-order phases in the anisotropic triangular lattice

We study the competition between magnetic and spin-liquid phases in the Hubbard model on the anisotropic triangular lattice, which is described by two hopping parameters t and t' in different spatial directions and is relevant for layered organic charge-transfer salts. By using a variational approach that includes spiral magnetic order, we provide solid evidence that a spin-liquid phase is stabilized in the strongly-correlated regime and close to the isotropic limit t'/t=1. Otherwise, a magnetically ordered spiral state is found, connecting the (collinear) Néel and the (coplanar) 120° phases. The pitch vector of the spiral phase obtained from the unrestricted Hartree-Fock approximation is substantially renormalized in presence of electronic correlations, and the Néel phase is stabilized in a wide regime of the phase diagram, i.e., for t'/t < 0.75. We discuss these results in the context of organic charge-transfer salts.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Variational Monte Carlo approach to the two-dimensional Kondo lattice model

We study the phase diagram of the Kondo-lattice model with nearest-neighbor hopping in the square lattice by means of the variational Monte Carlo technique. Specifically, we analyze a wide class of variational wave functions that allow magnetic and superconducting order parameters, so to assess the possibility that superconductivity might emerge close to the magnetic instability, as often observed in heavy fermion systems. Indeed, we do find evidence of d-wave superconductivity in the paramagnetic sector, i.e., when magnetic order is not allowed in the variational wave function. However, when magnetism is allowed, it completely covers the superconducting region, which thus disappears from the phase diagram.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Absence of static stripes in the two-dimensional $t{-}J$ model by an accurate and systematic quantum Monte Carlo approach

We examine the two-dimensional $t{-}J$ model by using variational approach combined with well established quantum Monte Carlo techniques [S. Sorella {\it et al.}, \prl {\bf 88}, 117002 (2002)] that are used to improve systematically the accuracy of the variational ansatz. Contrary to recent density-matrix renormalization group and projected entangled-pair state calculations [P. Corboz {\it et al.}, \prb {\bf 84}, 041108(R) (2011)], a uniform phase is found for $J/t=0.4$, even when the calculation is biased with an ansatz that explicitly contains stripe order. Moreover, in the small hole doping regime, i.e., $δ\lesssim 0.1$, our results support the coexistence of antiferromagnetism and superconductivity.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Localization and Glassy Dynamics Of Many-Body Quantum Systems

When classical systems fail to explore their entire configurational space, intriguing macroscopic phenomena like aging and glass formation may emerge. Also closed quanto-mechanical systems may stop wandering freely around the whole Hilbert space, even if they are initially prepared into a macroscopically large combination of eigenstates. Here, we report numerical evidences that the dynamics of strongly interacting lattice bosons driven sufficiently far from equilibrium can be trapped into extremely long-lived inhomogeneous metastable states. The slowing down of incoherent density excitations above a threshold energy, much reminiscent of a dynamical arrest on the verge of a glass transition, is identified as the key feature of this phenomenon. We argue that the resulting long-lived inhomogeneities are responsible for the lack of thermalization observed in large systems. Such a rich phenomenology could be experimentally uncovered upon probing the out-of-equilibrium dynamics of conveniently prepared quantum states of trapped cold atoms which we hereby suggest.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Strong renormalization of the Fermi-surface topology close to the Mott transition

The underlying Fermi surface is a key concept for strongly-interacting electron models and has been introduced to generalize the usual notion of the Fermi surface to generic (superconducting or insulating) systems. By using improved correlated wave functions that contain backflow and Jastrow terms, we examine the two-dimensional $t{-}t^\prime$ Hubbard model and find a non-trivial renormalization of the topology of the underlying Fermi surface close to the Mott insulator. Moreover, we observe a sharp crossover region, which arises from the metal-insulator transition, from a weakly interacting metal at small coupling to a resonating valence-bond superconductor at intermediate coupling. A violation of the Luttinger theorem is detected at low hole dopings.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Valence-bond crystals in the kagome spin-1/2 Heisenberg antiferromagnet: a symmetry classification and projected wave function study

In this paper, we do a complete classification of valence-bond crystals (VBCs) on the kagome lattice based on general arguments of symmetry only and thus identify many new VBCs for different unit cell sizes. For the spin-1/2 Heisenberg antiferromagnet, we study the relative energetics of competing gapless spin liquids (SLs) and VBC phases within the class of Gutzwiller-projected fermionic wave functions using variational Monte Carlo techniques, hence implementing exactly the constraint of one fermion per site. By using a state-of-the-art optimization method, we conclusively show that the U(1) Dirac SL is remarkably stable towards dimerizing into all 6-, 12- and 36-site unit cell VBCs. This stability is also preserved on addition of a next-nearest-neighbor super-exchange coupling of both antiferromagnetic and ferromagnetic (FM) type. However, we find that a 36-site unit cell VBC is stabilized on addition of a very small next-nearest-neighbor FM super-exchange coupling, i.e. |J2|~0.045, and this VBC is the same in terms of space-group symmetry as that obtained in an effective quantum dimer model study. It breaks reflection symmetry, has a nontrivial flux pattern and is a strong dimerization of the uniform RVB SL.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Backflow correlations in the Hubbard model: an efficient tool for the metal-insulator transition and the large-U limit

We show that backflow correlations in the variational wave function for the Hubbard model greatly improve the previous results given by the Slater-Jastrow state, usually considered in this context. We provide evidence that, within this approach, it is possible to have a satisfactory connection with the strong-coupling regime. Moreover, we show that, for the Hubbard model on the lattice, backflow correlations are essentially short range, inducing an effective attraction between empty (holons) and doubly occupied sites (doublons). In presence of frustration, we report the evidence that the metal to Mott-insulator transition is marked by a discontinuity of the double occupancy, together with a similar discontinuity of the kinetic term that does not change the number of holons and doublons, while the other kinetic terms are continuous across the transition. Finally, we show the estimation of the charge gap, obtained by particle-hole excitations {\it à la Feynman} over the ground-state wave function.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Bose-glass, superfluid, and rung-Mott phases of hard-core bosons in disordered two-leg ladders

By means of Monte Carlo techniques, we study the role of disorder on a system of hard-core bosons in a two-leg ladder with both intra-chain ($t$) and inter-chain ($t^\prime$) hoppings. We find that the phase diagram as a function of the boson density, disorder strength, and $t^\prime/t$ is far from being trivial. This contrasts the case of spin-less fermions where standard localization arguments apply and an Anderson-localized phase pervades the whole phase diagram. A compressible Bose-glass phase always intrudes between the Mott insulator with zero (or one) bosons per site and the superfluid that is stabilized for weak disorder. At half filling, there is a direct transition between a (gapped) rung-Mott insulator and a Bose glass, which is driven by exponentially rare regions where disorder is suppressed. Finally, by doping the rung-Mott insulator, a direct transition to the superfluid is possible only in the clean system, whereas the Mott phase is always surrounded by the a Bose glass when disorder is present. The phase diagram based on our numerical evidence is finally reported.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Phase diagram of hard-core bosons on a frustrated zig-zag ladder

We study hard-core bosons with unfrustrated nearest-neighbor hopping $t$ and repulsive interaction $V$ on a zig-zag ladder. As a function of the boson density $ρ$ and $V/t$, the ground state displays different quantum phases. A standard one-component Tomonaga-Luttinger liquid is stable for $ρ<1/3$ (and $ρ>2/3$) at any value of $V/t$. At commensurate densities $ρ=1/3$, $1/2$, and $2/3$ insulating (crystalline) phases are stabilized for a sufficiently large interaction $V$. For intermediate densities $1/3<ρ<2/3$ and large $V/t$, the ground state shows a clear evidence of a bound state of two bosons, implying gapped single-particle excitations but gapless excitations of boson pairs. Here, the low-energy properties may be described by a two-component Tomonaga-Luttinger liquid with a finite gap in the antisymmetric sector. Finally, for the same range of boson densities and weak interactions, the system is again a one-component Tomonaga-Luttinger liquid with no evidence of any breaking of discrete symmetries, in contrast to the frustrated case, where a ${\cal Z}_2$ symmetry breaking has been predicted.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Projected wave function study of Z2 spin liquids on the kagome lattice for the spin-1/2 quantum Heisenberg antiferromagnet

Motivated by recent density-matrix renormalization group (DMRG) calculations [Yan, Huse, and White, Science 332, 1173 (2011)], which claimed that the ground state of the nearest-neighbor spin-1/2 Heisenberg antiferromagnet on the kagome lattice geometry is a fully gapped spin liquid with numerical signatures of Z2 gauge structure, and a further theoretical work [Lu, Ran, and Lee, Phys. Rev. B 83, 224413 (2011)], which gave a classification of all Schwinger-fermion mean-field fully symmetric Z2 spin liquids on the kagome lattice, we have thoroughly studied Gutzwiller-projected fermionic wave functions by using quantum variational Monte Carlo techniques, hence implementing exactly the constraint of one fermion per site. In particular, we investigated the energetics of all Z2 candidates (gapped and gapless) that lie in the neighborhood of the energetically competitive U(1) gapless spin liquids. By using a state-of-the-art optimization method, we were able to conclusively show that the U(1) Dirac state is remarkably stable with respect to all Z2 spin liquids in its neighborhood, and in particular for opening a gap toward the so-called Z2[0,π]β state, which was conjectured to describe the ground state obtained by the DMRG method. Finally, we also considered the addition of a small second nearest-neighbor exchange coupling of both antiferromagnetic and ferromagnetic type, and obtained similar results, namely, a U(1) Dirac spin-liquid ground state.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Valence-bond crystal in the extended kagome spin-1/2 quantum Heisenberg antiferromagnet: A variational Monte Carlo approach

The highly-frustrated spin-1/2 quantum Heisenberg model with both nearest ($J_1$) and next-nearest ($J_2$) neighbor exchange interactions is revisited by using an extended variational space of projected wave functions that are optimized with state-of-the-art methods. Competition between modulated valence-bond crystals (VBCs) proposed in the literature and the Dirac spin liquid (DSL) is investigated. We find that the addition of a {\it small} ferromagnetic next-nearest-neighbor exchange coupling $|J_2|>0.09 J_1$ leads to stabilization of a 36-site unit cell VBC, although the DSL remains a local minimum of the variational parameter landscape. This implies that the VBC is not trivially connected to the DSL: instead it possesses a non-trivial flux pattern and large dimerization.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Characterization of the Bose-glass phase in low-dimensional lattices

We study by numerical simulation a disordered Bose-Hubbard model in low-dimensional lattices. We show that a proper characterization of the phase diagram on finite disordered clusters requires the knowledge of probability distributions of physical quantities rather than their averages. This holds in particular for determining the stability region of the Bose-glass phase, the compressible but not superfluid phase that exists whenever disorder is present. This result suggests that a similar statistical analysis should be performed also to interpret experiments on cold gases trapped in disordered lattices, limited as they are to finite sizes.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Extracting the Mott gap from energy measurements in trapped atomic gases

We show that the measure of the so-called {\it release-energy}, which is an experimentally accessible quantity, makes it possible to assess the value of the Mott gap in the presence of the confinement potential that is unavoidable in the actual experimental setup. Indeed, the curve of the release-energy as a function of the total number of particles shows kinks that are directly related to the existence of excitation gaps. Calculations are presented within the Gutzwiller approach, but the final results go beyond this simple approximation and represent a genuine feature of the real system. In the case of harmonic confinement, the Mott gaps may be renormalized with respect to the uniform case. On the other hand, in the case of the recently proposed off-diagonal confinement, our results show an almost perfect agreement with the homogeneous case.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Ground-state properties of the disordered Hubbard model in two dimensions

We study the interplay between electron correlation and disorder in the two-dimensional Hubbard model at half-filling by means of a variational wave function that can interpolate between Anderson and Mott insulators. We give a detailed description of our improved variational state and explain how the physics of the Anderson-Mott transition can be inferred from equal-time correlations functions, which can be easily computed within the variational Monte Carlo scheme. The ground-state phase diagram is worked out in both the paramagnetic and the magnetic sector. Whereas in the former a direct second-order Anderson-Mott transition is obtained, when magnetism is allowed variationally, we find evidence for the formation of local magnetic moments that order before the Mott transition. Although the localization length increases before the Mott transition, we have no evidence for the stabilization of a true metallic phase. The effect of a frustrating next-nearest-neighbor hopping $t^\prime$ is also studied in some detail. In particular, we show that $t^\prime$ has two primary effects. The first one is the narrowing of the stability region of the magnetic Anderson insulator, also leading to a first-order magnetic transition. The second and most important effect of a frustrating hopping term is the development of a ``glassy'' phase at strong couplings, where many paramagnetic states, with disordered local moments, may be stabilized.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Interaction induced Fermi-surface renormalization in the $t_1{-}t_2$ Hubbard model close to the Mott-Hubbard transition

We investigate the nature of the interaction-driven Mott-Hubbard transition of the half-filled $t_1{-}t_2$ Hubbard model in one dimension, using a full-fledged variational Monte Carlo approach including a distance-dependent Jastrow factor and backflow correlations. We present data for the evolution of the magnetic properties across the Mott-Hubbard transition and on the commensurate to incommensurate transition in the insulating state. Analyzing renormalized excitation spectra, we find that the Fermi surface renormalizes to perfect nesting right at the Mott-Hubbard transition in the insulating state, with a first-order reorganization when crossing into the conducting state.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Reptation quantum Monte Carlo for lattice Hamiltonians with a directed-update scheme

We provide an extension to lattice systems of the reptation quantum Monte Carlo algorithm, originally devised for continuous Hamiltonians. For systems affected by the sign problem, a method to systematically improve upon the so-called fixed-node approximation is also proposed. The generality of the method, which also takes advantage of a canonical worm algorithm scheme to measure off-diagonal observables, makes it applicable to a vast variety of quantum systems and eases the study of their ground-state and excited-states properties. As a case study, we investigate the quantum dynamics of the one-dimensional Heisenberg model and we provide accurate estimates of the ground-state energy of the two-dimensional fermionic Hubbard model.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The itinerant ferromagnetic phase of the Hubbard model

Using a newly developed quantum Monte Carlo technique, we provide strong evidence for the stability of a saturated ferromagnetic phase in the high-density regime of the two-dimensional infinite-U Hubbard model. By decreasing the electron density, a discontinuous transition to a paramagnetic phase is observed, accompanied by a divergence of the susceptibility on the paramagnetic side. This behavior, resulting from a high degeneracy among different spin sectors, is consistent with an infinite-order phase transition. The remarkable stability of itinerant ferromagnetism renews the hope to describe this phenomenon within a purely kinetic mechanism and will facilitate the validation of experimental quantum simulators with cold atoms loaded in optical lattices.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Spin-liquid and magnetic phases in the anisotropic triangular lattice: the case of $κ$-(ET)$_2$X

The two-dimensional Hubbard model on the anisotropic triangular lattice, with two different hopping amplitudes $t$ and $t^\prime$, is relevant to describe the low-energy physics of $κ$-(ET)$_2$X, a family of organic salts. The ground-state properties of this model are studied by using Monte Carlo techniques, on the basis of a recent definition of backflow correlations for strongly-correlated lattice systems. The results show that there is no magnetic order for reasonably large values of the electron-electron interaction $U$ and frustrating ratio $t^\prime/t = 0.85$, suitable to describe the non-magnetic compound with X=Cu$_2$(CN)$_3$. On the contrary, Néel order takes place for weaker frustrations, i.e., $t^\prime/t \sim 0.4 ÷0.6$, suitable for materials with X=Cu$_2$(SCN)$_2$, Cu[N(CN)$_2$]Cl, or Cu[N(CN)$_2$]Br.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Metal-insulator transition and strong-coupling spin liquid in the $t{-}t^\prime$ Hubbard model

We study the phase diagram of the frustrated $t{-}t^\prime$ Hubbard model on the square lattice by using a novel variational wave function. Taking the clue from the backflow correlations that have been introduced long-time ago by Feynman and Cohen and have been used for describing various interacting systems on the continuum (like liquid $^3$He, the electron jellium, and metallic Hydrogen), we consider many-body correlations to construct a suitable approximation for the ground state of this correlated model on the lattice. In this way, a very accurate {\it ansatz} can be achieved both at weak and strong coupling. We present the evidence that an insulating and non-magnetic phase can be stabilized at strong coupling and sufficiently large frustrating ratio $t^\prime/t$.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Stripes and spin-incommensurabilities are favored by lattice anisotropies

Structural distortions in cuprate materials give a natural origin for anisotropies in electron properties. We study a modified one-band t-J model in which we allow for different hoppings and antiferromagnetic couplings in the two spatial directions ($t_x \ne t_y$ and $J_x \ne J_y$). Incommensurate peaks in the spin structure factor show up only in the presence of a lattice anisotropy, whereas charge correlations, indicating enhanced fluctuations at incommensurate wave vectors, are almost unaffected with respect to the isotropic case.

0 citations0 reviewsNo code linkedOpen access