We have examined the behavior of the compressibility, the dc-conductivity, the single-particle gap, and the Drude weight as probes of the density-driven metal-insulator transition in the Hubbard model on a square lattice. These quantities have been o
btained through determinantal quantum Monte Carlo simulations at finite temperatures on lattices up to 16 X 16 sites. While the compressibility, the dc-conductivity, and the gap are known to suffer from `closed-shell effects due to the presence of artificial gaps in the spectrum (caused by the finiteness of the lattices), we have established that the former tracks the average sign of the fermionic determinant (<sign>), and that a shortcut often used to calculate the conductivity may neglect important corrections. Our systematic analyses also show that, by contrast, the Drude weight is not too sensitive to finite-size effects, being much more reliable as a probe to the insulating state. We have also investigated the influence of the discrete imaginary-time interval (Deltatau) on <sign>, on the average density (rho), and on the double occupancy (d): we have found that <sign> and rho are more strongly dependent on Delta tau away from closed-shell configurations, but d follows the Deltatau^2 dependence in both closed- and open-shell cases.
We study the attractive fermionic Hubbard model on a honeycomb lattice using determinantal quantum Monte Carlo simulations. By increasing the interaction strength U (relative to the hopping parameter t) at half-filling and zero temperature, the syste
m undergoes a quantum phase transition at 5.0 < U_c/t < 5.1 from a semi-metal to a phase displaying simultaneously superfluid behavior and density order. Doping away from half-filling, and increasing the interaction strength at finite but low temperature T, the system always appears to be a superfluid exhibiting a crossover between a BCS and a molecular regime. These different regimes are analyzed by studying the spectral function. The formation of pairs and the emergence of phase coherence throughout the sample are studied as U is increased and T is lowered.
The single band, two dimensional Hubbard Hamiltonian has been extensively studied as a model for high temperature superconductivity. While Quantum Monte Carlo simulations within the dynamic cluster approximation are now providing considerable evidenc
e for a d-wave superconducting state at low temperature, such a transition remains well out of reach of finite lattice simulations because of the sign problem. We show here that a bilayer Hubbard model, in which one layer is electron doped and one layer is hole doped, can be studied to lower temperatures and exhibits an interesting signal of d-wave pairing. The results of our simulations bear resemblance to a recent report on the magnetic and superconducting properties of Ba$_2$Ca$_3$Cu$_4$O$_8$F$_2$ which contains both electron and hole doped CuO$_2$ planes. We also explore the phase diagram of bilayer models in which each sheet is at half-filling.
The `dynamic Hubbard Hamiltonian describes interacting fermions on a lattice whose on-site repulsion is modulated by a coupling to a fluctuating bosonic field. We investigate one such model, introduced by Hirsch, using the determinant Quantum Monte C
arlo method. Our key result is that the extended s-wave pairing vertex, repulsive in the usual static Hubbard model, becomes attractive as the coupling to the fluctuating Bose field increases. The sign problem prevents us from exploring a low enough temperature to see if a superconducting transition occurs. We also observe a stabilization of antiferromagnetic correlations and the Mott gap near half-filling, and a near linear behavior of the energy as a function of particle density which indicates a tendency toward phase separation.
We investigate the phases of the ionic Hubbard model in a two-dimensional square lattice using determinant quantum Monte Carlo (DQMC). At half-filling, when the interaction strength or the staggered potential dominate we find Mott and band insulators
, respectively. When these two energies are of the same order we find a metallic region. Charge and magnetic structure factors demonstrate the presence of antiferromagnetism only in the Mott region, although the externally imposed density modulation is present everywhere in the phase diagram. Away from half-filling, other insulating phases are found. Kinetic energy correlations do not give clear signals for the existence of a bond-ordered phase.
