Abstract We obtain ab initio the Hubbard parameters for Rydberg-dressed atoms in a one-dimensional (1D) sinusoidal optical lattice on the basis of maximally-localized Wannier states. Finite range, soft-core interatomic interactions become the trait of Rydberg admixed atoms, which can be extended over many neighboring lattice sites. In contrast to dipolar gases, where the interactions follow an inverse cubic law, the key feature of Rydberg-dressed interactions is the possibility of making neighboring couplings to the same magnitude as that of the onsite ones. The maximally-localized Wannier functions (MLWFs) are typically calculated via a spread-minimization procedure (Marzari N and Vanderbilt D 1997 Phys. Rev. B 56 12847 ) and are always found to be real functions apart from a trivial global phase when an isolated set of Bloch bands are considered. For an isolated single Bloch band, the above procedure reduces to a simple quasi-momentum-dependent unitary phase transformation. Here, instead of minimizing the spread, we employ a diagonal phase transformation which eliminates the imaginary part of the Wannier functions. The resulting Wannier states are found to be maximally localized and in exact agreement with those obtained via a spread-minimization procedure. Using these findings, we calculate the Hubbard couplings from the Rydberg admixed interactions, including dominant density-assisted tunneling (DAT) coefficients. Finally, we provide realistic lattice parameters for the state-of-the-art experimental Rydberg-dressed rubidium setup.
Abstract We obtain ab initio the Hubbard parameters for Rydberg-dressed atoms in a one-dimensional (1D) sinusoidal optical lattice on the basis of maximally-localized Wannier states. Finite range, soft-core interatomic interactions become the trait of Rydberg admixed atoms, which can be extended over many neighboring lattice sites. In contrast to dipolar gases, where the interactions follow an inverse cubic law, the key feature of Rydberg-dressed interactions is the possibility of making neighboring couplings to the same magnitude as that of the onsite ones. The maximally-localized Wannier functions (MLWFs) are typically calculated via a spread-minimization procedure (Marzari N and Vanderbilt D 1997 Phys. Rev. B 56 12847 ) and are always found to be real functions apart from a trivial global phase when an isolated set of Bloch bands are considered. For an isolated single Bloch band, the above procedure reduces to a simple quasi-momentum-dependent unitary phase transformation. Here, instead of minimizing the spread, we employ a diagonal phase transformation which eliminates the imaginary part of the Wannier functions. The resulting Wannier states are found to be maximally localized and in exact agreement with those obtained via a spread-minimization procedure. Using these findings, we calculate the Hubbard couplings from the Rydberg admixed interactions, including dominant density-assisted tunneling (DAT) coefficients. Finally, we provide realistic lattice parameters for the state-of-the-art experimental Rydberg-dressed rubidium setup.
We derive ab initio local Hubbard models for several optical-lattice potentials of current interest, including the honeycomb and kagome lattices, verifying their accuracy on each occasion by comparing the interpolated band structures against the originals. To achieve this, we calculate the maximally localized generalized Wannier basis by implementing the steepest-descent algorithm of Marzari and Vanderbilt [Phys. Rev. B 56, 12847 (1997)] directly in one and two dimensions. To avoid local minima we develop an initialization procedure that is both robust and requires no prior knowledge of the optimal Wannier basis.
Optical lattices provide ideal experimental tools for simulating a wide variety of physical systems. They exhibit excellent coherence properties, their microscopic dynamics are well-understood, and they can be dynamically controlled with great precision. Recent advances in cooling and trapping gases of dipolar molecules and Rydberg atoms have further increased the versatility of optical lattices as quantum simulators, particularly for solid state systems possessing long-range interactions. In this thesis we further examine the microscopic properties of cold atoms trapped in an optical lattice, going beyond the case where the lattice axes are orthogonal. We present a transformation of the single-particle energy eigenstates to a complete, orthonormal basis of so-called 'generalised Wannier functions' that are localised to each lattice site. Generalised Wannier functions, as opposed to ordinary Wannier functions, are formed by mixing states between bands and are required when the bands are degenerate. We modify an algorithm devised by Marzari and Vanderbilt (Phys. Rev. B 56(20) 12847) for calculating maximally localised generalised Wannier functions, specifically to determine these functions for optical lattice systems of arbitrary geometry, including those with degenerate lower bands. Moreover, we overcome an issue with the initialisation of the original algorithm to find the maximally localised set without fail. We then present results for a variety of optical lattice systems in one- and two-dimensions, including hexagonal and Kagomé geometries. In each case we use the generalised Wannier functions to calculate the hopping and interaction parameters, and thus determine a reduced form of the system's Hamiltonian that contains only local interactions. We also study two spin-1/2 systems with long-range interactions as paradigm examples of the models that can be simulated using optical lattices loaded with dipolar molecules. We first examine the critical properties of a 1D spin-1/2 Ising chain possessing exponentially decaying interactions as the range of the interactions is increased. Using the infinite density matrix renormalisation group method, we calculate the location of the critical point and observe a smooth approach to the classical mean-field result in the limit of infinite-ranged interactions. We go on to consider the quantum dynamics of a single spin-1/2 coupled to a bath of interacting spins as a model for decoherence in solid state quantum memories. The large degree of symmetry of the bath, described by the infinite-ranged Lipkin-Meshkov-Glick model, allows us to find parameter regimes where the initial qubit state is revived at well-defined times after the qubit preparation. We observe that these times may become independent of the bath size for large baths and thus enable faithful qubit storage even in the presence of strong coupling to the bath.
Ultracold atomic gases in optical lattices provide an unique framework to study quantum phenomena in strongly correlated systems. In addition to the precise control over all relevant parameters in the experiment, these experiments can be mapped to a fundamental model from solid-state physics. For the bosonic version of this model, the so-called Bose-Hubbard model, a phase transition from a superfluid to a Mott insulator was theoretically predicted and later experimentally observed in an ultracold gas of Rb-87 atoms in three-dimensional as well as in one-dimensional optical lattices. Apart from homogeneous optical lattices one can introduce more complex lattice topologies such as two-color superlattices which give rise to a rich phase diagram including more exotic phases like the Bose-glass. We employ various powerful many-body techniques like exact diagonalization in complete and truncated Hilbert spaces and the Density-Matrix Renormalization Group (DMRG) algorithm to study the phase diagrams of the one-dimensional Bose-Hubbard and the one-dimensional Bose-Fermi-Hubbard Hamiltonian. Most theoretical studies of these systems discuss the phase diagrams with respect to the generic parameters of the Hubbard model. These Hubbard parameters, however, depend non-trivially on the control parameters used in experiments. The focus of this work is on the ab-initio calculation of the phase diagram of ultracold Rb-87 in one-dimensional optical superlattices starting directly from the experimental setup. To this end, we first employ band-structure calculations to extract the Hubbard parameters from the experimental parameters. Then, we use state-of-the-art DMRG techniques to solve the many-body problem for realistic particle numbers and lattice sizes that occur in experiments. Our results show that by using the intensities of the two laser fields forming the two-color superlattice as control parameters while keeping all other experimental parameters fixed, it is possible to access all relevant quantum phases of the system. Furthermore, we found out that the critical values of the laser intensities for the different phase transitions depend strongly on a third parameter that has to be included for a realistic description of the experiment. This third parameter is the strength of a harmonic trapping potential which accounts for the Gaussian shape of the laser fields and an additional magnetic potential used to confine the atoms in the center of the trap.
Originally, the Hubbard model was derived for describing the behavior of strongly correlated electrons in solids. However, for over a decade now, variations of it have also routinely been implemented with ultracold atoms in optical lattices, allowing their study in a clean, essentially defect-free environment. Here, we review some of the vast literature on this subject, with a focus on more recent non-standard forms of the Hubbard model. After giving an introduction to standard (fermionic and bosonic) Hubbard models, we discuss briefly common models for mixtures, as well as the so-called extended Bose-Hubbard models, that include interactions between neighboring sites, next-neighbor sites, and so on. The main part of the review discusses the importance of additional terms appearing when refining the tight-binding approximation for the original physical Hamiltonian. Even when restricting the models to the lowest Bloch band is justified, the standard approach neglects the density-induced tunneling (which has the same origin as the usual on-site interaction). The importance of these contributions is discussed for both contact and dipolar interactions. For sufficiently strong interactions, the effects related to higher Bloch bands also become important even for deep optical lattices. Different approaches that aim at incorporating these effects, mainly via dressing the basis, Wannier functions with interactions, leading to effective, density-dependent Hubbard-type models, are reviewed. We discuss also examples of Hubbard-like models that explicitly involve higher p orbitals, as well as models that dynamically couple spin and orbital degrees of freedom. Finally, we review mean-field nonlinear Schrödinger models of the Salerno type that share with the non-standard Hubbard models nonlinear coupling between the adjacent sites. In that part, discrete solitons are the main subject of consideration. We conclude by listing some open problems, to be addressed in the future.
We propose a numerical method using the discrete variable representation (DVR) for constructing real-valued Wannier functions localized in a unit cell for both symmetric and asymmetric periodic potentials. We apply these results to finding Wannier functions for ultracold atoms trapped in laser-generated optical lattices. Following S. Kivelson [Phys. Rev. B 26, 4269 (1982)], for a symmetric lattice with inversion symmetry, we construct Wannier functions as eigenstates of the position operators $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{x}, \stackrel{\ifmmode \hat{}\else \^{}\fi{}}{y}$, and $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{z}$ restricted to single-particle Bloch functions belonging to one or more bands. To ensure that the Wannier functions are real-valued, we numerically obtain the band structure and real-valued eigenstates using a uniform Fourier grid DVR. We then show, by a comparison of tunneling energies, that the Wannier functions are accurate for both inversion-symmetric and asymmetric potentials to better than 10 significant digits when using double-precision arithmetic. The calculations are performed for an optical lattice with double-wells per unit cell with tunable asymmetry along the x axis and a single sinusoidal potential along the perpendicular directions. Localized functions at the two potential minima within each unit cell are similarly constructed, but using a superposition of single-particle solutions from the two lowest bands. We finally use these localized basis functions to determine the two-body interaction energies in the Bose-Hubbard model and show the dependence of these energies on lattice asymmetry.
In systems with periodic potential fields, building relatively local Wannier functions can significantly simplify the Hamiltonian and enhance our understanding of the system's ground state and dynamic properties. In this work, we improve the current method of building the Wannier functions of ultracold atomic systems, including the case in the presence or absence of interactions. In noninteracting systems, we propose a method to directly obtain the real-valued maximally localized Wannier functions (MLWFs) by using real-valued eigenstates, and verify the effectiveness of this method in a two-dimensional (2D) degenerate system. In interacting systems, we obtain the effect of high-energy bands on the lowest-energy band by using the accurate calculation results of the two-particle system. In the two-particle system, we consider the effect of the entanglement between the particles and obtain the optimal two-particle Wannier functions. These Wannier functions are then further utilized to obtain the parameters of the extended Bose-Hubbard model. The effectiveness of the method is verified by taking a one-dimensional (1D) system with contact interaction as an example. In the three-particle and four-particle systems, compared calculation results with the original system and the unmodified two-band Bose-Hubbard model, we find that the effective Hamiltonian is more accurate than the unmodified two-band model. This verifies the effectiveness of our method, and the parameters obtained can reflect the original system well, which provides an effective method for accurate modeling of interacting systems.