Abstract

Interacting one-dimensional quantum systems play a pivotal role in physics. Exact solutions can be obtained for the homogeneous case using the Bethe ansatz and bosonisation techniques. However, these approaches are not applicable when external confinement is present. Recent theoretical advances beyond the Bethe ansatz and bosonisation allow us to predict the behaviour of one-dimensional confined systems with strong short-range interactions, and new experiments with cold atomic Fermi gases have already confirmed these theories. Here we demonstrate that a simple linear combination of the strongly interacting solution with the well-known solution in the limit of vanishing interactions provides a simple and accurate description of the system for all values of the interaction strength. This indicates that one can indeed capture the physics of confined one-dimensional systems by knowledge of the limits using wave functions that are much easier to handle than the output of typical numerical approaches. We demonstrate our scheme for experimentally relevant systems with up to six particles. Moreover, we show that our method works also in the case of mixed systems of particles with different masses. This is an important feature because these systems are known to be non-integrable and thus not solvable by the Bethe ansatz technique.

Highlights

  • Compared to entirely numerical approaches, which typically have the solution represented on a large basis set with many non-zero contributions

  • To describe the main idea we will focus on two-component Fermi systems of N↑ particles with spin projection up and N↓ particles with spin projection down

  • While we only discuss the two-component Fermi system below, the case of bosons or Bose-Fermi mixtures is similar in spirit and in formalism and we return briefly to this extension in the outlook

Read more

Summary

Results and Discussion

Let |γ0〉be an energy eigenstate in the non-interaction limit (that is, g = 0) with eigenenergy E0. We use a variational approach to solving the Schrödinger equation by identifying stationary points of the trial state energy functional Eq (5). We select the values of α0 and α∞ that optimise the energy of the trial state for a given value of g. This will yield energies and eigenstates that, approximate, turn out to be extremely accurate as discussed below. The coefficients that yield stationary points of Eq (5) are given by (see the Supplementary Materials for details) α0 α∞

This gives the energy
We can then select λ such that
Methods
Author Contributions
Additional Information
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call