Source author record

G. G. Batrouni

G. G. Batrouni appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

22works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

22 published item(s)

preprint2021arXiv

Pairing and superconductivity in quasi one-dimensional flat band systems: Creutz and sawtooth lattices

We study the pairing and superconducting properties of the attractive Hubbard model in two quasi one-dimensional topological lattices: the Creutz and sawtooth lattices. They share two peculiar properties: each of their band structures exhibits a flat band with a non-trivial winding number. The difference, however, is that only the Creutz lattice is genuinely topological, owing to a chiral (sub-lattice) symmetry, resulting in a quantized winding number and zero energy edge modes for open boundary conditions. We use mean field and exact density matrix renormalization group in our work. Our three main results are: (a) For both lattice systems, the superconducting weight, $D_s$, is linear in the coupling strength, $U$, for low values of $U$; (b) for small $U$, $D_s$ is proportional to the quantum metric for the Creutz system but not for the sawtooth system because its sublattices are not equivalent; (c) conventional BCS mean field is not appropriate for such systems with inequivalent sublattices. We show that, for a wide range of densities and coupling strengths, these systems are very well described by a full multi-band mean field method where the pairing parameters and the local particle densities on the inequivalent sublattices are variational mean field parameters.

preprint2014arXiv

Charge Gaps at Fractional Fillings in Boson Hubbard Ladders

The Bose-Hubbard Hamiltonian describes the competition between superfluidity and Mott insulating behavior at zero temperature and commensurate filling as the strength of the on-site repulsion is varied. Gapped insulating phases also occur at non-integer densities as a consequence of longer ranged repulsive interactions. In this paper we explore the formation of gapped phases in coupled chains due instead to anisotropies $t_x \neq t_y$ in the bosonic hopping, extending the work of Crepin {\it et al.} [Phys. Rev. B 84, 054517 (2011)] on two coupled chains, where a gap was shown to occur at half filling for arbitrarily small interchain hopping $t_y$. Our main result is that, unlike the two-leg chains, for three- and four-leg chains, a charge gap requires a finite nonzero critical $t_y$ to open. However, these finite values are surprisingly small, well below the analogous values required for a fermionic band gap to open.

preprint2014arXiv

Competing phases, phase separation and co-existence in the extended one-dimensional bosonic Hubbard model

We study the phase diagram of the one-dimensional bosonic Hubbard model with contact ($U$) and near neighbor ($V$) interactions focusing on the gapped Haldane insulating (HI) phase which is characterized by an exotic nonlocal order parameter. The parameter regime ($U$, $V$ and $μ$) where this phase exists and how it competes with other phases such as the supersolid (SS) phase, is incompletely understood. We use the Stochastic Green Function quantum Monte Carlo algorithm as well as the density matrix renormalization group to map out the phase diagram. Our main conclusions are that the HI exists only at $ρ=1$, the SS phase exists for a very wide range of parameters (including commensurate fillings) and displays power law decay in the one body Green function. In addition, we show that at fixed integer density, the system exhibits phase separation in the $(U,V)$ plane.

preprint2014arXiv

Superconducting Transitions in Flat Band Systems

The physics of strongly correlated quantum particles within a flat band was originally explored as a route to itinerant ferromagnetism and, indeed, a celebrated theorem by Lieb rigorously establishes that the ground state of the repulsive Hubbard model on a bipartite lattice with unequal number of sites in each sublattice must have nonzero spin S at half-filling. Recently, there has been interest in Lieb geometries due to the possibility of novel topological insulator, nematic, and Bose-Einstein condensed (BEC) phases. In this paper, we extend the understanding of the attractive Hubbard model on the Lieb lattice by using Determinant Quantum Monte Carlo to study real space charge and pair correlation functions not addressed by the Lieb theorems.

preprint2013arXiv

Competing Supersolid and Haldane Insulator phases in the extended one-dimensional bosonic Hubbard model

The Haldane Insulator is a gapped phase characterized by an exotic non-local order parameter. The parameter regimes at which it might exist, and how it competes with alternate types of order, such as supersolid order, are still incompletely understood. Using the Stochastic Green Function (SGF) quantum Monte Carlo (QMC) and the Density Matrix Renormalization Group (DMRG), we study numerically the ground state phase diagram of the one-dimensional bosonic Hubbard model (BHM) with contact and near neighbor repulsive interactions. We show that, depending on the ratio of the near neighbor to contact interactions, this model exhibits charge density waves (CDW), superfluid (SF), supersolid (SS) and the recently identified Haldane insulating (HI) phases. We show that the HI exists only at the tip of the unit filling CDW lobe and that there is a stable SS phase over a very wide range of parameters.

preprint2013arXiv

Exotic phases of interacting p-band bosons

We study a model of interacting bosons that occupy the first excited p-band states of a two-dimensional optical lattice. In contrast to the much studied single band Bose-Hubbard Hamiltonian, this more complex model allows for non-trivial superfluid phases associated with condensation at non-zero momentum and staggered order of the orbital angular momentum in addition to the superfluid-Mott insulator transition. More specifically, we observe staggered orbital angular momentum order in the Mott phase at commensurate filling and superfluidity at all densities. We also observe a transition between the staggered angular momentum superfluid phase and a striped superfluid, with an alternation of the phase of the superfluid along one direction. The transition between these two phases was observed in a recent experiment, which is then qualitatively well described by our model.

preprint2013arXiv

Interacting spin-1 bosons in a two-dimensional optical lattice

We study, using quantum Monte Carlo (QMC) simulations, the ground state properties of spin-1 bosons trapped in a square optical lattice. The phase diagram is characterized by the mobility of the particles (Mott insulating or superfluid phase) and by their magnetic properties. For ferromagnetic on-site interactions, the whole phase diagram is ferromagnetic and the Mott insulators-superfluid phase transitions are second order. For antiferromagnetic on-site interactions, spin nematic order is found in the odd Mott lobes and in the superfluid phase. Furthermore, the superfluid-insulator phase transition is first or second order depending on whether the density in the Mott is even or odd. Inside the even Mott lobes, we observe a singlet-to-nematic transition for certain values of the interactions. This transition appears to be first order.

preprint2013arXiv

Phase Stability in the Two dimensional Anisotropic Boson Hubbard Hamiltonian

The two dimensional square lattice hard-core boson Hubbard model with near neighbor interactions has a `checkerboard' charge density wave insulating phase at half-filling and sufficiently large intersite repulsion. When doped, rather than forming a supersolid phase in which long range charge density wave correlations coexist with a condensation of superfluid defects, the system instead phase separates. However, it is known that there are other lattice geometries and interaction patterns for which such coexistence takes place. In this paper we explore the possibility that anisotropic hopping or anisotropic near neighbor repulsion might similarly stabilize the square lattice supersolid. By considering the charge density wave structure factor and superfluid density for different ratios of interaction strength and hybridization in the $\hat x$ and $\hat y$ directions, we conclude that phase separation still occurs.

preprint2012arXiv

Density of States and Magnetic Correlations at a Metal-Mott Insulator Interface

The possibility of novel behavior at interfaces between strongly and weakly correlated materials has come under increased study recently. In this paper, we use determinant Quantum Monte Carlo to determine the inter-penetration of metallic and Mott insulator physics across an interface in the two dimensional Hubbard Hamiltonian. We quantify the behavior of the density of states at the Fermi level and the short and long range antiferromagnetism as functions of the distance from the interface and with different interaction strength, temperature and hopping across the interface. Induced metallic behavior into the insulator is evident over several lattice spacings, whereas antiferromagnetic correlations remain small on the metallic side. At large interface hopping, singlets form between the two boundary layers, shielding the two systems from each other.

preprint2012arXiv

Finite temperature phase diagram of spin-1/2 bosons in two-dimensional optical lattice

We study a two-species bosonic Hubbard model on a two-dimensional square lattice by means of quantum Monte Carlo simulations and focus on finite temperature effects. We show in two different cases, ferro- and antiferromagnetic spin-spin interactions, that the phase diagram is composed of solid Mott phases, liquid phases and superfluid phases. In the antiferromagnetic case, the superfluid (SF) is polarized while the Mott insulator (MI) and normal Bose liquid (NBL) phases are not. On the other hand, in the ferromagnetic case, none of the phases is polarized. The superfluid-liquid transition is of the Berezinsky-Kosterlitz-Thouless type whereas the solid-liquid passage is a crossover.

preprint2012arXiv

Gapless inhomogeneous superfluid phase with spin-dependent disorder

We show that the presence of a spin-dependent random potential in a superconductor or a superfluid atomic gas leads to distinct transitions at which the energy gap and average order parameter vanish, generating an intermediate gapless superfluid phase, in marked contrast to the case of spin-symmetric randomness where no such gapless superfluid phase is seen. By allowing the pairing amplitude to become inhomogeneous, the gapless superconducting phase persists to considerably higher disorder compared with the prediction of Abrikosov-Gorkov. The low-lying excited states are located predominantly in regions where the pairing amplitude vanishes and coexist with the superfluid regions with a finite pairing. Our results are based on inhomogeneous Bogoliubov-de Gennes mean field theory for a two dimensional attractive Hubbard model with spin-dependent disorder.

preprint2012arXiv

Ground state phase diagram of spin-1/2 bosons in a two-dimensional optical lattice

We study a two-species bosonic Hubbard model on a two-dimensional square lattice by means of quantum Monte Carlo simulations. In addition to the usual contact repulsive interactions between the particles, the Hamiltonian has an interconversion term which allows the transformation of two particles from one species to the other. The phases are characterized by their solid or superfluid properties and by their polarization, i.e. the difference in the populations. When inter-species interactions are smaller than the intra-species ones, the system is unpolarized, whereas in the opposite case the system is unpolarized in even Mott insulator lobes and polarized in odd Mott lobes and also in the superfluid phase. We show that in the latter case the transition between the Mott insulator of total density two and the superfluid can be either of second or first order depending on the relative values of the interactions, whereas the transitions are continuous in all other cases.

preprint2012arXiv

Optimized Confinement of Fermions in Two Dimensions

One of the challenging features of studying model Hamiltonians with cold atoms in optical lattices is the presence of spatial inhomogeneities induced by the confining potential, which results in the coexistence of different phases. This paper presents Quantum Monte Carlo results comparing meth- ods for confining fermions in two dimensions, including conventional diagonal confinement (DC), a recently proposed 'off-diagonal confinement' (ODC), as well as a trap which produces uniform den- sity in the lattice. At constant entropy and for currently accessible temperatures, we show that the current DC method results in the strongest magnetic signature, primarily because of its judicious use of entropy sinks at the lattice edge. For d-wave pairing, we show that a constant density trap has the more robust signal and that ODC can implement a constant density profile. This feature is important to any prospective search for superconductivity in optical lattices.

preprint2012arXiv

Pairing in a two-dimensional Fermi gas with population imbalance

Pairing in a population imbalanced Fermi system in a two-dimensional optical lattice is studied using Determinant Quantum Monte Carlo (DQMC) simulations and mean-field calculations. The approximation-free numerical results show a wide range of stability of the Fulde-Ferrell-Larkin-Ovshinnikov (FFLO) phase. Contrary to claims of fragility with increased dimensionality we find that this phase is stable across wide range of values for the polarization, temperature and interaction strength. Both homogeneous and harmonically trapped systems display pairing with finite center of mass momentum, with clear signatures either in momentum space or real space, which could be observed in cold atomic gases loaded in an optical lattice. We also use the harmonic level basis in the confined system and find that pairs can form between particles occupying different levels which can be seen as the analog of the finite center of mass momentum pairing in the translationally invariant case. Finally, we perform mean field calculations for the uniform and confined systems and show the results to be in good agreement with QMC. This leads to a simple picture of the different pairing mechanisms, depending on the filling and confining potential.

preprint2011arXiv

Pairing in population imbalanced Fermion systems

We use Quantum Monte Carlo (QMC) simulations to study the pairing mechanism in a one-dimensional fermionic system governed by the Hubbard model with attractive contact interaction and with imbalance between the two spin populations. This is done for the uniform system and also for the system confined in a harmonic trap to compare with experiments on confined ultra-cold atoms. In the uniform case we determine the phase diagram in the polarization-temperature plane and find that the "Fulde-Ferrell-Larkin-Ovchinnikov" (FFLO) phase is robust and persists to higher temperature for higher polarization. In the confined case, we also find that the FFLO phase is stabilized by higher polarization and that it is within the range of detection of experiments currently underway.

preprint2011arXiv

Superconductivity and charge order of confined Fermi systems

The low-temperature properties of the two-dimensional attractive Hubbard model are strongly influenced by the fermion density. Away from half-filling, there is a finite-temperature transition to a phase with s-wave pairing order. However, the critical temperature is suppressed to zero at half-filling, where long-range charge-density-wave order also appears, degenerate with superconductivity. This paper presents Determinant Quantum Monte Carlo simulations of the attractive Hubbard model in the presence of a confining potential V which makes the fermion density ρ inhomogeneous across the lattice. Pair correlations are shown to be large at low temperatures in regions of the trapped system with incommensurate filling, and to exhibit a minimum as the local density ρ(i) passes through one fermion per site. In this ring of ρ=1, charge order is enhanced. A comparison is made between treating V within the local-density approximation (LDA) and in an ab initio manner. It is argued that certain sharp features of the LDA result at integer filling do not survive the proximity of doped sites. The critical temperature of confined systems of fixed characteristic density is estimated.

preprint2010arXiv

Finite temperature QMC study of the one-dimensional polarized Fermi gas

Quantum Monte Carlo (QMC) techniques are used to provide an approximation-free investigation of the phases of the one-dimensional attractive Hubbard Hamiltonian in the presence of population imbalance. The temperature at which the "Fulde-Ferrell-Larkin-Ovchinnikov" (FFLO) phase is destroyed by thermal fluctuations is determined as a function of the polarization. It is shown that the presence of a confining potential does not dramatically alter the FFLO regime, and that recent experiments on trapped atomic gases likely lie just within the stable temperature range.

preprint2010arXiv

Isentropes of spin-1 bosons in an optical lattice

We analyze the effects of adiabatic ramping of optical lattices on the temperature of spin-1 bosons in a homogeneous lattice. Using mean-field theory, we present the isentropes in the temperature-interaction strength ($T,U_0$) plane for ferromagnetic, antiferromagnetic, and zero spin couplings. Following the isentropic lines, temperature changes can be determined during adiabatic loading of current experiments. We show that the heating-cooling separatrix lies on the superfluid-Mott phase boundary with cooling occuring within the superfluid and heating in the Mott insulator, and quantify the effects of spin coupling on the heating rate. We find that the mean-field isentropes for low initial entropy terminate at the superfluid-Mott insulator phase boundary.

preprint2010arXiv

Phase diagram of spin-1/2 bosons in one-dimensional optical lattice

Systems of two coupled bosonic species are studied using Mean Field Theory and Quantum Monte Carlo. The phase diagram is characterized both based on the mobility of the particles (Mott insulating or superfluid) and whether or not the system is magnetic (different populations for the two species). The phase diagram is shown to be population balanced for negative spin-dependent interactions, regardless of whether it is insulating or superfluid. For positive spin-dependent interactions, the superfluid phase is always polarized, the two populations are imbalanced. On the other hand, the Mott insulating phase with even commensurate filling has balanced populations while the odd commensurate filling Mott phase has balanced populations at very strong interaction and polarizes as the interaction gets weaker while still in the Mott phase.

preprint2010arXiv

Pure Mott phases in confined ultracold atomic systems

We propose a novel scheme for confining atoms to optical lattices by engineering a spatially-inhomogeneous hopping matrix element in the Hubbard-model (HM) description, a situation we term off-diagonal confinement (ODC). We show, via an exact numerical solution of the boson HM with ODC, that this scheme possesses distinct advantages over the conventional method of confining atoms using an additional trapping potential, including the presence of incompressible Mott phases at commensurate filling and a phase diagram that is similar to the uniform HM. The experimental implementation of ODC will thus allow a more faithful realization of correlated phases of interest in cold atom experiments.

preprint2008arXiv

Monte Carlo Study of an Inhomogeneous Blume-Capel Model

Systems of particles in a confining potential exhibit a spatially dependent density which fundamentally alters the nature of phase transitions that occur. A specific instance of this situation, which is being extensively explored currently, concerns the properties of ultra-cold, optically trapped atoms. Of interest is how the superfluid-insulator transition is modified by the inhomogeneity, and, indeed, the extent to which a sharp transition survives at all. This paper explores a classical analog of these systems, the Blume-Capel model with a spatially varying single ion anisotropy and/or temperature gradient. We present results both for the nature of the critical properties and for the validity of the "local density approximation" which is often used to model the inhomogeneous case. We compare situations when the underlying uniform transition is first and second order.

preprint2007arXiv

Quantum Monte Carlo study of a one-dimensional phase-fluctuating condensate

We study numerically the low temperature behavior of a one-dimensional Bose gas trapped in an optical lattice. For a sufficient number of particles and weak repulsive interactions, we find a clear regime of temperatures where density fluctuations are negligible but phase fluctuations are considerable, i.e., a quasicondensate. In the weakly interacting limit, our results are in very good agreement with those obtained using a mean-field approximation. In coupling regimes beyond the validity of mean-field approaches, a phase-fluctuating condensate also appears, but the phase-correlation properties are qualitatively different. It is shown that quantum depletion plays an important role.