No Arabic abstract
We use unbiased computational methods to elucidate the onset and properties of pair superfluidity in two-species fermionic and bosonic systems with onsite interspecies attraction loaded in one-dimensional optical lattice. We compare results from quantum Monte Carlo (QMC) and density matrix renormalization group (DMRG), emphasizing the one-to-one correspondence between the Drude weight tensor, calculated with DMRG, and the various winding numbers extracted from the QMC. Our results show that, for any nonvanishing attractive interaction, pairs form and are the sole contributors to superfluidity, there are no individual contributions due to the separate species. For weak attraction, the pair size diverges exponentially, i.e. Bardeen-Cooper-Schrieffer (BCS) pairing requiring huge systems to bring out the pair-only nature of the superfluid. This crucial property is largely overlooked in many studies, thereby misinterpreting the origin and nature of the superfluid. We compare and contrast this with the repulsive case and show that the behavior is very different, contradicting previous claims about drag superfluidity and the symmetry of properties for attractive and repulsive interactions. Finally, our results show that the situation is similar for soft core bosons: superfluidity is due only to pairs, even for the smallest attractive interaction strength compatible with the largest system sizes that we could attain.
We calculate the superfluid weight and the polarization amplitude for the one-dimensional bosonic Hubbard model focusing on the strong-coupling regime. Other than analytic calculations we apply two methods: variational Monte Carlo based on the Baeriswyl wave function and exact diagonalization. The former gives zero superfluid response at integer filling, while the latter gives a superfluid response at finite hopping. From the polarization amplitude we derive the variance and the associated size scaling exponent. Again, the variational study does not produce a finite superfluid weight at integer filling (size scaling exponent is near one), but the Fourier transform of the polarization amplitude behaves in a similar way to the result of exact diagonalization, with a peak at small hopping, and suddenly decreasing at the insulator-superfluid transition. On the other hand, exact diagonalization studies result in a finite spread of the total position which increases with the size of the system. In the superfluid phase the size scaling exponent is two as expected. Importantly, our work addresses the ambiguities that arise in the calculation of the superfluid weight in variational calculations, and we comment on the prediction of Anderson about the superfluid response of the model at integer filling.
We study the ground state phase diagram of a nonlinear two-photon Rabi-Hubbard (RH) model in one dimension using quantum Monte Carlo (QMC) simulations and density matrix renormalization group (DMRG) calculations. Our model includes a nonlinear photon-photon interaction term. Absent this term, the RH model has only one phase, the normal disordered phase, and suffers from spectral collapse at larger values of the photon-qubit interaction or inter-cavity photon hopping. The photon-photon interaction, no matter how small, stabilizes the system which now exhibits {it two} quantum phase transitions: Normal phase to {it photon pair} superfluid (PSF) transition and PSF to single particle superfluid (SPSF). The discrete $Z_4$ symmetry of the Hamiltonian spontaneously breaks in two stages: First it breaks partially as the system enters the PSF and then completely breaks when the system finally enters the SPSF phase. We show detailed numerical results supporting this, and map out the ground state phase diagram.
We address some open questions regarding the phase diagram of the one-dimensional Hubbard model with asymmetric hopping coefficients and balanced species. In the attractive regime we present a numerical study of the passage from on-site pairing dominant correlations at small asymmetries to charge-density waves in the region with markedly different hopping coefficients. In the repulsive regime we exploit two analytical treatments in the strong- and weak-coupling regimes in order to locate the onset of phase separation at small and large asymmetries respectively.
In ultracold atoms settings, inelastic light scattering is a preeminent technique to reveal static and dynamic properties at nonzero momentum. In this work, we investigate an array of one-dimensional trapped Bose gases, by measuring both the energy and the momentum imparted to the system via light scattering experiments. The measurements are performed in the weak perturbation regime, where these two quantities - the energy and momentum transferred - are expected to be related to the dynamical structure factor of the system. We discuss this relation, with special attention to the role of in-trap dynamics on the transferred momentum.
The density distribution of the one-dimensional Hubbard model in a harmonic trapping potential is investigated in order to study the effect of the confining trap. Strong superimposed oscillations are always present on top of a uniform density cloud, which show universal scaling behavior as a function of increasing interactions. An analytical formula is proposed on the basis of bosonization, which describes the density oscillations for all interaction strengths. The wavelength of the dominant oscillation changes with interaction, which indicates the crossover to a spin-incoherent regime. Using the Bethe ansatz the shape of the uniform fermion cloud is analyzed in detail, which can be described by a universal scaling form.