Researcher profile

Federico Becca

Federico Becca contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

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

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.

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

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.

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.

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.

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$&#39;s, showing that the bosonization approach is not able to represent the time evolution of generic one-dimensional models after a quantum quench.

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&#39; 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&#39;/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&#39;/t < 0.75. We discuss these results in the context of organic charge-transfer salts.

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.

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.

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.

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.

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.

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.

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.

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