No Arabic abstract
We use the density-matrix renormalization group method to investigate ground-state and dynamic properties of the one-dimensional Bose-Hubbard model, the effective model of ultracold bosonic atoms in an optical lattice. For fixed maximum site occupancy $n_b=5$, we calculate the phase boundaries between the Mott insulator and the `superfluid phase for the lowest two Mott lobes. We extract the Tomonaga-Luttinger parameter from the density-density correlation function and determine accurately the critical interaction strength for the Mott transition. For both phases, we study the momentum distribution function in the homogeneous system, and the particle distribution and quasi-momentum distribution functions in a parabolic trap. With our zero-temperature method we determine the photoemission spectra in the Mott insulator and in the `superfluid phase of the one-dimensional Bose-Hubbard model. In the insulator, the Mott gap separates the quasi-particle and quasi-hole dispersions. In the `superfluid phase the spectral weight is concentrated around zero momentum.
We study finite-temperature transport properties of the one-dimensional Hubbard model using the density matrix renormalization group. Our aim is two-fold: First, we compute both the charge and the spin current correlation function of the integrable model at half filling. The former decays rapidly, implying that the corresponding Drude weight is either zero or very small. Second, we calculate the optical charge conductivity sigma(omega) in presence of small integrability-breaking next-nearest neighbor interactions (the extended Hubbard model). The DC conductivity is finite and diverges as the temperature is decreased below the gap. Our results thus suggest that the half-filled, gapped Hubbard model is a normal charge conductor at finite temperatures. As a testbed for our numerics, we compute sigma(omega) for the integrable XXZ spin chain in its gapped phase.
We compute the phase diagram of the one-dimensional Bose-Hubbard model with a quasi-periodic potential by means of the density-matrix renormalization group technique. This model describes the physics of cold atoms loaded in an optical lattice in the presence of a superlattice potential whose wave length is incommensurate with the main lattice wave length. After discussing the conditions under which the model can be realized experimentally, the study of the density vs. the chemical potential curves for a non-trapped system unveils the existence of gapped phases at incommensurate densities interpreted as incommensurate charge-density wave phases. Furthermore, a localization transition is known to occur above a critical value of the potential depth V_2 in the case of free and hard-core bosons. We extend these results to soft-core bosons for which the phase diagrams at fixed densities display new features compared with the phase diagrams known for random box distribution disorder. In particular, a direct transition from the superfluid phase to the Mott insulating phase is found at finite V_2. Evidence for reentrances of the superfluid phase upon increasing interactions is presented. We finally comment on different ways to probe the emergent quantum phases and most importantly, the existence of a critical value for the localization transition. The later feature can be investigated by looking at the expansion of the cloud after releasing the trap.
We study the real-time and real-space dynamics of charge in the one-dimensional Hubbard model in the limit of high temperatures. To this end, we prepare pure initial states with sharply peaked density profiles and calculate the time evolution of these nonequilibrium states, by using numerical forward-propagation approaches to chains as long as 20 sites. For a class of typical states, we find excellent agreement with linear-response theory and unveil the existence of remarkably clean charge diffusion in the regime of strong particle-particle interactions. Moreover, we demonstrate that this diffusive behavior does not depend on certain details of our initial conditions, i.e., it occurs for five different realizations with random and nonrandom internal degrees of freedom, single and double occupation of the central site, and displacement of spin-up and spin-down particles.
Non-Hermitian phenomena offer a novel approach to analyze and interpret spectra in the presence of interactions. Using the density-matrix renormalization group (DMRG), we demonstrate the existence of exceptional points for the one-particle Greens function of the 1D alternating Hubbard chain with chiral symmetry, with a corresponding Fermi arc at zero frequency in the spectrum. They result from the non-Hermiticity of the effective Hamiltonian describing the Greens function and only appear at finite temperature. They are robust and can be topologically characterized by the zeroth Chern number. This effect illustrates a case where temperature has a strong effect in 1D beyond the simple broadening of spectral features. Finally, we demonstrate that exceptional points appear even in the two-particle Greens function (charge structure factor) where an effective Hamiltonian is difficult to establish, but move away from zero frequency due to a distinct symmetry constraint.
We study the zero-temperature phase diagram of the half-filled one-dimensional ionic Hubbard model. This model is governed by the interplay of the on-site Coulomb repulsion and an alternating one-particle potential. Various many-body energy gaps, the charge-density-wave and bond-order parameters, the electric as well as the bond-order susceptibilities, and the density-density correlation function are calculated using the density-matrix renormalization group method. In order to obtain a comprehensive picture, we investigate systems with open as well as periodic boundary conditions and study the physical properties in different sectors of the phase diagram. A careful finite-size scaling analysis leads to results which give strong evidence in favor of a scenario with two quantum critical points and an intermediate spontaneously dimerized phase. Our results indicate that the phase transitions are continuous. Using a scaling ansatz we are able to read off critical exponents at the first critical point. In contrast to a bosonization approach, we do not find Ising critical exponents. We show that the low-energy physics of the strong coupling phase can only partly be understood in terms of the strong coupling behavior of the ordinary Hubbard model.