Size and shape of Mott regions for fermionic atoms in a two-dimensional optical lattice

We investigate the harmonic-trap control of size and shape of Mott regions in the Fermi Hubbard model on a square optical lattice. The use of Lanczos diagonalization on clusters with twisted boundary conditions, followed by an average over 50-80 samples, drastically reduce finite-size effects in some ground state properties; calculations in the grand canonical ensemble together with a local-density approximation (LDA) allow us to simulate the radial density distribution. We have found that as the trap closes, the atomic cloud goes from a metallic state, to a Mott core, and to a Mott ring; the coverage of Mott atoms reaches a maximum at the core-ring transition. A `phase diagram' in terms of an effective density and the on-site repulsion is proposed, as a guide to maximize the Mott coverage. We also predict that the usual experimentally accessible quantities, the global compressibility and the average double occupancy (rather, its density derivative) display detectable signatures of the core-ring transition. Some spin correlation functions are also calculated, and predict the existence N\'eel ordering within Mott cores and rings.Comment: 5 pages, 6 figure

Multiperiodic magnetic structures in Hubbard superlattices

We consider fermions in one-dimensional superlattices (SL's), modeled by site-dependent Hubbard-U couplings arranged in a repeated pattern of repulsive (i.e., U>0) and free (U=0) sites. Density Matrix Renormalization Group (DMRG) diagonalization of finite systems is used to calculate the local moment and the magnetic structure factor in the ground state. We have found four regimes for magnetic behavior: uniform local moments forming a spin-density wave (SDW), `floppy' local moments with short-ranged correlations, local moments on repulsive sites forming long-period SDW's superimposed with short-ranged correlations, and local moments on repulsive sites solely with long-period SDW's; the boundaries between these regimes depend on the range of electronic densities, rho, and on the SL aspect ratio. Above a critical electronic density, rho_{uparrow downarrow}, the SDW period oscillates both with rho and with the spacer thickness. The former oscillation allows one to reproduce all SDW wave vectors within a small range of electronic densities, unlike the homogeneous system. The latter oscillation is related to the exchange oscillation observed in magnetic multilayers. A crossover between regimes of `thin' to `thick' layers has also been observed.Comment: 9 two-column pages, 10 figure

Luttinger liquid superlattices: realization of gapless insulating phases

We investigate Luttinger Liquid superlattices, a periodic structure composed of two kinds of one-dimensional systems of interacting electrons. We calculate several properties of the low-energy sector: the effective charge and spin velocities, the compressibility, various correlation functions, the Landauer conductance and the Drude weight. The low-energy properties are subsumed into effective parameters, much like homogeneous one-dimensional systems. A generic result is the weighted average nature of these parameters, in proportion to the spatial extent of the underlying subunits, pointing to the possibility of ``engineered'' structures. As a specific realization, we consider a one-dimensional Hubbard superlattice, which consists of a periodic arrangement of two long Hubbard chains with different coupling constants and different hopping amplitudes. This system exhibits a rich phase diagram with several phases, both metallic and insulating. We have found that gapless insulating phases are present over a wide range of parameters.Comment: 16 pages, 15 figures, RevTeX

Disordered two-dimensional superconductors: roles of temperature and interaction strength

We have considered the half-filled disordered attractive Hubbard model on a square lattice, in which the on-site attraction is switched off on a fraction $f$ of sites, while keeping a finite $U$ on the remaining ones. Through Quantum Monte Carlo (QMC) simulations for several values of $f$ and $U$, and for system sizes ranging from $8\times 8$ to $16\times 16$, we have calculated the configurational averages of the equal-time pair structure factor $P_s$, and, for a more restricted set of variables, the helicity modulus, $\rho_s$, as functions of temperature. Two finite-size scaling {\it ansatze} for $P_s$ have been used, one for zero-temperature and the other for finite temperatures. We have found that the system sustains superconductivity in the ground state up to a critical impurity concentration, $f_c$, which increases with $U$, at least up to U=4 (in units of the hopping energy). Also, the normalized zero-temperature gap as a function of $f$ shows a maximum near $f\sim 0.07$, for $2\lesssim U\lesssim 6$. Analyses of the helicity modulus and of the pair structure factor led to the determination of the critical temperature as a function of $f$, for $U=3,$ 4 and 6: they also show maxima near $f\sim 0.07$, with the highest $T_c$ increasing with $U$ in this range. We argue that, overall, the observed behavior results from both the breakdown of CDW-superconductivity degeneracy and the fact that free sites tend to "push" electrons towards attractive sites, the latter effect being more drastic at weak couplings.Comment: 9 two-column pages, 14 figures, RevTe