Federico Becca

Superconductivity in the two-dimensional t-J model

Using computational techniques, it is shown that pairing is a robust property of hole-doped antiferromagnetic insulators. In one dimension and for two-leg ladder systems, a BCS-like variational wave function with long-bond spin singlets and a Jastrow factor provides an accurate representation of the ground state of the t-J model, even though strong quantum fluctuations destroy the off-diagonal superconducting long-range order in this case. However, in two dimensions it is argued-and numerically confirmed using several techniques, especially quantum Monte Carlo-that quantum fluctuations are not strong enough to suppress superconductivity.

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 multireference projected Hartree-Fock methods. 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.

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.

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

Dynamical behavior across the Mott transition of two bands with different bandwidths

We investigate the role of the bandwidth difference in the Mott metal-insulator transition of a two-band Hubbard model in the limit of infinite dimensions by means of a Gutzwiller variational wave function as well as by dynamical mean-field theory. The variational calculation predicts a two-stage quenching of the charge degrees of freedom, in which the narrower band undergoes a Mott transition before the wider one, both in the presence and in the absence of a Hunds exchange coupling. However, this scenario is not fully confirmed by the dynamical mean-field theory calculation, which shows that, although the quasiparticle residue of the narrower band is zero within our numerical accuracy, low-energy spectral weight still exists inside the Mott-Hubbard gap, concentrated into two peaks symmetric around the chemical potential. This spectral weight vanishes only when the wider band ceases to conduct too. Although our results are compatible with several scenarios\char22{}e.g., a narrow-gap semiconductor or a semimetal\char22{}we argue that the most plausible one is that the two peaks coexist with a narrow resonance tied at the chemical potential, with a spectral weight below our numerical accuracy. This quasiparticle resonance is expected to vanish when the wider band undergoes the Mott transition.

Quantum compass model on the square lattice

Using exact diagonalizations, Greens function Monte Carlo simulations and high-order perturbation theory, we study the low-energy properties of the two-dimensional spin-1/2 compass model on the square lattice defined by the Hamiltonian

Direct evidence for a gaplessZ2spin liquid by frustrating Néel antiferromagnetism

H = - \sum_{\bm{r}} (J_x \sigma_{\bm{r}}^x \sigma_{\bm{r} + \bm{e}_x}^x + J_z \sigma_{\bm{r}}^z \sigma_{\bm{r} + \bm{e}_z}^z)

Charge-density waves and superconductivity as an alternative to phase separation in the infinite- U Hubbard-Holstein model

. When

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

J_x\ne J_z

Projected wave function study of Z

, we show that, on clusters of dimension