Simulation of inhomogeneous distributions of ultracold atoms in an optical lattice via a massively parallel implementation of nonequilibrium strong-coupling perturbation theory
arXiv:1309.5994 · doi:10.1103/PhysRevE.89.023306
Abstract
We present a nonequilibrium strong-coupling approach to inhomogeneous systems of ultracold atoms in optical lattices. We demonstrate its application to the Mott-insulating phase of a two-dimensional Fermi-Hubbard model in the presence of a trap potential. Since the theory is formulated self-consistently, the numerical implementation relies on a massively parallel evaluation of the self-energy and the Green's function at each lattice site, employing thousands of CPUs. While the computation of the self-energy is straightforward to parallelize, the evaluation of the Green's function requires the inversion of a large sparse matrix, with . As a crucial ingredient, our solution heavily relies on the smallness of the hopping as compared to the interaction strength and yields a widely scalable realization of a rapidly converging iterative algorithm which evaluates all elements of the Green's function. Results are validated by comparing with the homogeneous case via the local-density approximation. These calculations also show that the local-density approximation is valid in non-equilibrium setups without mass transport.
14 pages, 9 figures
References in corpus (7)
- Many-Body Physics with Ultracold Gases
- Repulsively bound atom pairs in an optical lattice
- Atomic three-body loss as a dynamical three-body interaction
- Spin-charge separation in cold Fermi-gases: a real time analysis
- Spectroscopy of ultracold atoms by periodic lattice modulations
- Probing nearest-neighbor correlations of ultracold fermions in an optical lattice
- Andreev-like reflections with cold atoms