Alejandro Muramatsu

Quantum spin liquid emerging in two-dimensional correlated Dirac fermions

At sufficiently low temperatures, condensed-matter systems tend to develop order. A notable exception to this behaviour is the case of quantum spin liquids, in which quantum fluctuations prevent a transition to an ordered state down to the lowest temperatures. There have now been tentative observations of such states in some two-dimensional organic compounds, yet quantum spin liquids remain elusive in microscopic two-dimensional models that are relevant to experiments. Here we show, by means of large-scale quantum Monte Carlo simulations of correlated fermions on a honeycomb lattice (a structure realized in, for example, graphene), that a quantum spin liquid emerges between the state described by massless Dirac fermions and an antiferromagnetically ordered Mott insulator. This unexpected quantum-disordered state is found to be a short-range resonating valence-bond liquid, akin to the one proposed for high-temperature superconductors: the possibility of unconventional superconductivity through doping therefore arises in our system. We foresee the experimental realization of this model system using ultra-cold atoms, or group IV elements arranged in honeycomb lattices.

Strongly Correlated Fermions after a Quantum Quench

Using the adaptive time-dependent density-matrix renormalization group method, we study the time evolution of strongly correlated spinless fermions on a one-dimensional lattice after a sudden change of the interaction strength. For certain parameter values, two different initial states (e.g., metallic and insulating) lead to observables which become indistinguishable after relaxation. We find that the resulting quasistationary state is nonthermal. This result holds for both integrable and nonintegrable variants of the system.

Mott Domains of Bosons Confined on Optical Lattices

In the absence of a confining potential, the boson-Hubbard model exhibits a superfluid to Mott insulator quantum phase transition at commensurate fillings and strong coupling. We use quantum Monte Carlo simulations to study the ground state of the one-dimensional bosonic Hubbard model in a trap. Some, but not all, aspects of the Mott insulating phase persist. Mott behavior occurs for a continuous range of incommensurate fillings, very different from the unconfined case, and the establishment of the Mott phase does not proceed via a traditional quantum phase transition. These results have important implications for interpreting experiments on ultracold atoms on optical lattices.

Local quantum criticality in confined fermions on optical lattices.

Using quantum Monte Carlo simulations, we show that the one-dimensional fermionic Hubbard model in a harmonic potential displays quantum critical behavior at the boundaries of a Mott-insulating region. A local compressibility defined to characterize the Mott-insulating phase has a nontrivial critical exponent. Both the local compressibility and the variance of the local density show universality with respect to the confining potential. We determine a generic phase diagram, which allows the prediction of the phases to be observed in experiments with ultracold fermionic atoms trapped on optical lattices.

Hard-core bosons on optical superlattices: Dynamics and relaxation in the superfluid and insulating regimes

We study the ground-state properties and nonequilibrium dynamics of hard-core bosons confined in one-dimensional lattices in the presence of an additional periodic potential (superlattice) and a harmonic trap. The dynamics is analyzed after a sudden switch-on or switch-off of the superlattice potential, which can bring the system into insulating or superfluid phases, respectively. A collapse and revival of the zero-momentum peak can be seen in the first case. We study in detail the relaxation of these integrable systems towards equilibrium. We show how after relaxation time averages of physical observables, like the momentum distribution function, can be predicted by means of a generalization of the Gibbs distribution.

Spectral properties of the one-dimensional Hubbard model.

The spectral properties of the 1-D Hubbard model are obtained from quantum Monte Carlo simulations using the maximum entropy method. The one-particle excitations are characterized by dispersive cosine-like bands. Velocities for spin- and charge excitations are obtained that lead to a conformal charge c=0.98 +/- 0.05 for the largest system simulated (N=84). An exact sum-rule for the spin-excitations is fulfilled accurately with deviations of at most 10% only around 2 kF.

Quantum Monte Carlo study of confined fermions in one-dimensional optical lattices

Using quantum Monte Carlo (QMC) simulations we study the ground-state properties of the one-dimensional fermionic Hubbard model in traps with an underlying lattice. Since due to the confining potential the density is space dependent, Mott-insulating domains always coexist with metallic regions, such that global quantities are not appropriate to describe the system. We define a local compressibility that characterizes the Mott-insulating regions and analyze other local quantities. It is shown that the momentum distribution function, a quantity that is commonly considered in experiments, fails in giving a clear signal of the Mott-insulating transition. Furthermore, we analyze a mean-field approach to these systems and compare it with the numerically exact QMC results. Finally, we determine a generic form for the phase diagram that allows us to predict the phases to be observed in the experiments.

Phase diagram of the three-dimensional Hubbard model at half filling

Abstract:We investigate the phase diagram of the three-dimensional Hubbard model at half filling using quantum Monte Carlo (QMC) simulations. The antiferromagnetic Néel temperature is determined from the specific heat maximum in combination with finite-size scaling of the magnetic structure factor. Our results interpolate smoothly between the asymptotic solutions for weak and strong coupling, respectively, in contrast to previous QMC simulations. The location of the metal-insulator transition in the paramagnetic phase above is determined using the electronic compressibility as criterion.

Time evolution of correlations in strongly interacting fermions after a quantum quench

Using the adaptive time-dependent density-matrix renormalization group, we study the time evolution of density correlations of interacting spinless fermions on a one-dimensional lattice after a sudden change in the interaction strength. Over a broad range of model parameters, the correlation function exhibits a characteristic light-cone-like time evolution representative of a ballistic transport of information. Such behavior is observed both when quenching an insulator into the metallic region and also when quenching within the insulating region. However, when a metallic state beyond the quantum critical point is quenched deep into the insulating regime, no indication for ballistic transport is observed. Instead, stable domain walls in the density correlations emerge during the time evolution, consistent with the predictions of the Kibble-Zurek mechanism.

Single-hole dynamics in thet−Jmodel on a square lattice

We present quantum Monte Carlo (QMC) simulations for a single hole in a