|Home | About | Journals | Submit | Contact Us | Français|
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.
Understanding the properties of low-dimensional systems is not merely an academic pursuit. Technologically promising systems such as nanotubes, nanowires, and organic conductors have one-dimensional nature1,2,3, while there is much evidence that high-temperature superconductors owe their spectacular properties to an effective two-dimensional structure4,5. However, in the case of interacting particles in one dimension there are still many outstanding fundamental issues in their quantum mechanical description. An avenue within which these problems can be studied is that of cold atomic gases6,7 where experiments in one-dimensional (1D) confinement can be performed with tunable interactions for systems of bosons8,9,10,11,12,13,14 or fermions15. Most recently, 1D Fermi systems have been constructed with full control over the particle number16 and thus engineered few-body systems are now available. This provides opportunities to study pairing, impurity physics, magnetism, and strongly interacting particles from the bottom up17,18,19,20,21.
The role of strong interactions in important quantum phenomena such as superconductivity and magnetism drives research into the regime of strong interactions also for 1D systems. In this respect, new theoretical approaches to confined quantum systems with strong short-range interactions have been proposed in the last couple of years22,23,24,25,26,27,28,29. Within these new developments it has become clear that the strongly interacting limit has an emergent Heisenberg spin model description. While this was realized some time ago for the case of a homogeneous system30, in the presence of confinement the Heisenberg model obtained has non-trivial nearest-neighbour interactions that depend on the confining potential. This can be exploited for tailoring systems to have desired static and dynamic properties25. In the opposite limit where we have vanishing interactions, confined 1D systems are trivially solved as the single-particle Schrödinger equation can be solved numerically to any desired level of accuracy. The natural question to ask is how to describe 1D systems with interaction strengths that are somewhere in the region between the two extreme limits?
In this paper, we propose a deceptively simple ansatz that linearly combines our knowledge of the weakly and strongly interacting limits. As an ansatz containing just two wave functions it is much easier to work with as 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. The Hamiltonian for the N=N↑+N↓ system is
where m is the mass and the sums run over the number of particles, i=1, …, N, xi is the coordinate, and pi the momentum of the ith particle. The external confinement, Vext(x), is assumed to be the same for all particles. Furthermore, we will assume that Vext(x) is parity invariant and has at least N bound states. The interaction strength g is positive for repulsive interactions and negative for attractive interactions. We note that since the parity operator commutes with the Hamiltonian, parity is conserved when g changes. 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. For the majority of the discussion we assume all particles have equal mass, m, but we also discuss the important extension to mass-imbalanced systems. Notice that in the case of two-component Fermi systems, particles of the same spin projection will not interact due to the antisymmetry required by the Pauli principle which implies that the zero-range interaction will vanish for similar spin projections. For this reason we may use the general form of the Hamiltonian in Eq. (1) with the sum over i<j also for fermionic systems.
Consider now an N-body system and assume that we have solved the problem for g=0 with energy eigenstate |γ0 and for 1/g=0 with eigenstate |γ∞. Now form the linear combination
This interpolatory ansatz is motivated by the intuitive idea that the wave function of a system with intermediate-strength interactions contains a mixture of qualities from the wave functions with weak and strong interactions. We may compute γ|H|γ as a function of α0 and α∞ and look for optimum values. As we shall demonstrate in this paper, |γ provides a simple yet accurate description of the system for any value of g. The only exception is deeply bound states to which we return below.
We will demonstrate that the interpolatory ansatz in Eq. (2) can capture the qualitative features of the eigenstates of the Hamiltonian in Eq. (1) and is quantitatively accurate at the level of a few percent. Furthermore, we show that the expression for the optimum energy of |γ may be modified slightly to make it perturbatively correct in both limits of the interaction strength. With this modification, the interpolatory ansatz provides an approximation to the eigenenergy with an accuracy that is comparable to state-of-the-art numerical methods, though the ansatz is far simpler.
We provide a proof of principle by considering some important examples from the few-body limit that are experimentally relevant at the moment. For this purpose we restrict to a harmonic confinement, i.e. , throughout (here ω is the angular frequency of the oscillator).
The results we present here are:
The ansatz we propose can be used to get very simple expressions for different observables as one needs to compute only a few matrix elements between the g=0 and 1/g=0 states of interest. Our method is directly extendable to bosonic systems or mixed systems, as long as one has access to the two limiting wave functions, and while we have focused on harmonically confined systems it can be straightforwardly extended to any other form of confinement. Furthermore, one may systematically improve the ansatz by adding more states.
Let |γ0 be an energy eigenstate in the non-interaction limit (that is, g=0) with eigenenergy E0. We now adiabatically change the interaction strength from g=0 to |g|→∞ and in turn we adiabatically change |γ0 into a new state denoted |γ∞ with energy E∞. The wave function of the state |γ∞ vanishes whenever the position coordinates of any two particles coincide, and |γ∞ is thus unaffected by the interaction potential
As our ansatz we now construct the trial state
where α0 and α∞ are real parameters. Assuming |γ0 and |γ∞ are normalised, the energy of the trial state is
where we let ΔEE∞−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 thus 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, although approximate, turn out to be extremely accurate as discussed below.
Before presenting the results of the variational calculation, we shall briefly examine Eq. (5) in the limiting cases of the interaction strength: If we require that α∞/α0→0 for g→0 such that the trial state approaches |γ0 in the non-interacting limit, Eq. (5) gives a first-order term in g of the form γ0|V|γ0, in agreement with first order perturbation theory. We caution, however, that requiring α0/α∞→0 (i.e. |γ→|γ∞) for 1/g→0, does not automatically ensure that the first order expansion of Eq. (5) in 1/g is equal to that of the exact eigenstate. We will come back to this point later on.
This gives the energy
The denominator in Eq. (7) is positive because |γ0|γ∞|≤1 by the Cauchy-Schwarz inequality. Hence, is the energy maximum while is the energy minimum. To ensure the correct energies in the limits g=0 and 1/g=0, we use the energy minimum, , to approximate the eigenenergy whenever g>0, while for g<0 we use the energy maximum, .
While the interpolatory ansatz is extremely simple, it has a shortcoming in the 1/g→0 limit where it does not reproduce the slope of the energy. As is shown in the following, we may, however, modify the ansatz slightly to correct this behaviour. Letting q−1/g, the first-order expansion has the slope
where K0=γ0|V|γ0/g is the corresponding slope of the energy curve in the limit of vanishing interactions. This demonstrates that the important quantity in the slope is the overlap γ0|γ∞. In the original philosophy of the ansatz, we exploit that we know both wave functions in this overlap exactly and thus also the overlap itself; this leaves no unfixed parameters. Realising that this yields a discrepancy we have explored how to modify this assumption in order to improve the approximation.
To this end, we note that the derivation of Eq. (7) actually does not depend on |γ∞ being an energy eigenstate. It must be a state with energy E∞ (with respect to the non-interacting Hamiltonian), but beyond that the only requirement of |γ∞ is that V|γ∞=0. Furthermore, Eopt only depends on |γ∞ through the squared wave-function overlap γ0|γ∞2. Thus, if we substitute γ0|γ∞2 in Eqs (7 and 8) with some parameter λ, we may regard Eopt and as functions of λ. We can then select λ such that Eopt(λ) becomes perturbatively correct, that is , or,
where is the slope of the true eigenenergy curve at q=0 (which is known exactly using the formalism of A.G.V. et al.22). We shall refer to this perturbatively correct modification, Eopt(λ), as the modified ansatz. Note that this modification breaks variational bounds since we have no a priori knowledge of any trial state whose energy is Eopt(λ). By going from the original interpolatory ansatz to the modified ansatz, we lose information about the wave function, but gain the correct slope of the energy at 1/g→0. Below we will see that the modified ansatz increases the accuracy significantly as compared to the original ansatz.
We now proceed to discuss examples of the ansatz. Throughout our discussion, the external potential is taken to be harmonic and the same for all particles. This is the most widely studied and experimentally relevant case at the moment so it will be our focus here. Henceforth, we use natural units =m=ω=1 such that energies are given in units of ω, lengths in units of , and interaction strengths g in units of σω.
The N=2 case is special as analytical results for any g are available due to the seminal work of Busch et al.31. It is therefore an important benchmark case for our approach. First we note that we will only be interested in the relative energy as the center of mass decouples in the harmonic trap. Note that this decoupling is not an essential assumption of our method and is merely a convenience.
The details on how to construct the ansatz states for the two-particle case are given in the Methods section below. The energy spectrum using the interpolatory ansatz of Eq. (7) is shown in Fig. 1. Only the even-parity solutions are shown on the plot as odd-parity states are unchanged by the zero-range interaction. The figure also includes experimental measurements of the ground-state energy. The experiment has been conducted by S. Jochim’s ultra-cold atoms group in Heidelberg, and the data originates from Wenz et al.19. This data has been corrected for the imperfections of the trap as described in Wenz et al.19. As seen, the experimental data agrees with the interpolatory ansatz well within the experimental uncertainties.
If we expand the energy in terms of g, the first-order term agrees with the result of ordinary non-degenerate perturbation theory. Similarly, in the limit 1/g=0, the energy curve (as a function of q) of the true eigenstate has the same slope as that of our interpolatory ansatz.
The inset in Fig. 1 shows a zoom of the energy spectrum on the ground state and compares the energy predicted by Eq. (7) with the exact solution. For q>0 the energy of the interpolatory ansatz is within 0.05 of the exact energy in the vicinity of q~0.4 and even less elsewhere. The error decreases as we move up in the spectrum. For q<0 the deviation is less than 0.006 for the ground state; again greatest around q~−0.4. On this side, the error increases for the excited states, but we find that it is bounded by about 0.03.
We conclude that the ansatz is extremely accurate for the N=2 case where we can compare to analytical results31. One may also compare the wave functions and again find extremely good agreement (see Supplementary Fig. S1).
For attractive interactions, a deeply bound molecular state exists that we have so far ignored. However, it turns out that the ansatz of Eq. (7) can be extended to also give extremely accurate results for the deeply bound state. As is shown in the Supplementary Materials, this can be done by including an additional state in the ansatz that has the correct asymptotic behaviour as g→−∞. We stress that this is in fact in complete agreement with the universal philosophy of the ansatz method, i.e. interpolation between (known) extremes. Thus to address deeply bound states one needs a state in the extreme limit of large negative energy. This yields a very precise approximation also for the deeply bound state. This highlights the universal nature of our approach. We will not pursue the deeply bound states any further in this paper.
The simplest non-trivial example of a three-body two-component Fermi system has N↑=2 and N↓=1. The interaction potential is
and we let x1 be the position of the spin-down fermion, while x2 and x3 denote the position of spin-up fermions. Again the third interaction term will vanish for identical fermions but we keep it for generality.
Eigenstates of the harmonic Hamiltonian are described by two quantum numbers, ν≥0 and μ≥1, which we shall call the radial quantum number and the angular quantum number, respectively. The eigenenergies are
The quantum numbers ν and μ can be used to describe the energy eigenstates in both the non-interacting limit and the infinite-interaction limit.
Recall that |γ0 and |γ∞ denote eigenstates in the non-interacting limit and the infinite-interaction limit, respectively, and that we are looking for states that are adiabatically connected. As states with different radial quantum number ν are orthogonal, we assume that the adiabatically connected states have the same ν values. Parity p is exactly conserved and thus also the same for two states in the ansatz. The quantity that changes with g is therefore the μ quantum number, and we call the limiting values μ0 and μ∞, respectively. The angular quantum numbers are related by
Figure 2 shows the energy spectrum for E<8 (deeply bound states are not considered) both as calculated using the interpolatory ansatz and by exact numerical diagonalisation. The interpolatory ansatz clearly describes the qualitative features of the spectrum well; both for repulsive and attractive interactions. In this energy-range, the ν=0 states also give a good quantitative match to the numerical results. However, the error seems to grow with ν.
We see in Fig. 2 that trial states with ν=0 and Δμ=1 offer a particularly good approximation to the corresponding eigenstates in the repulsive region (see e.g., the first excited state). This is because the slope of the energy curve, that is , at q=0 is the same as that of the exact eigenenergy. For the trial states with Δμ=2, this is generally not the case. This result suggests that is an important quantity to reproduce correctly in an attempt such as the present to describe the physics of our problem through a simple ansatz. Note also, that the slope of the energy curve has a discontinuity at 1/g=0, which contradicts the expectation that the states go smoothly through this region.
If we now enforce the correct slope by a modification of the interpolatory ansatz as proposed in the discussion following Eq. (8), we arrive at the spectrum shown in Fig. 3. We see that the modified ansatz agrees better with the numerical results than the original ansatz; especially for states with ν>0. There are, however, still deviations on the attractive side of the spectrum, and for high energies, also on the repulsive side. We shall give a quantitative discussion of the quality of the approximation when discussing the impurity system below. We have included the experimental measurements of the ground-state energy19 in the figure and see that both ansatz and modified ansatz agrees very nicely with experiment, although for large g the modified ansatz naturally does better.
In the case of different masses, one typically uses another length scale given by where , and m1 is the mass of the spin-down fermion while m2 and m3 are the masses of the spin-up fermions. We now consider the case where m1=M and m2=m3=m.
In Fig. 4 we show the energy spectrum obtained by the interpolatory ansatz for M/m=1/2 and M/m=2 and compare this with numerically calculated results using the correlated Gaussian approach (see the Supplementary Materials for details on the numerical methods). The agreement between the ansatz and the numerical results is striking for the low-energy part of the spectrum considered here, and we see that the ansatz can be extended also to mass-imbalanced systems.
We now consider a system of N fermions among which one particle (x1) is spin-down and the N↑=N−1 remaining particles (x2, …, xN) are spin-up. Taking into account that interactions between identical fermions vanish, we can write the interaction term as
We restrict the discussion to the ground state with repulsive interactions as this has been a focus of recent experimental attention19. The considerations can be extended to obtain more states in the spectrum, to the attractive side, and to deeply bound states in the same manner as in the previous examples.
We denote by |γA the Slater determinant of the single-particle harmonic eigenstates xk|n=ψn(xk) for k=1, …, N and n=0, …, N−1. Here ψn(x) is the single-particle eigenstate of the harmonic oscillator Hamiltonian in coordinate space with quantum number n and argument x. The state |γA is antisymmetric with respect to interchange of any two coordinates, so γA|V|γA=0. We can write the antisymmetric state as
if we define
where x=(x1, …, xN), SN is the symmetric group of order N, and Π(σ) indicates the integration region xσ(1)<…<xσ(N). In each region, Π(σ), the wave function of the ground state in the infinite-interaction limit is proportional to x|γA. The (normalised) ground state in the infinite-interaction limit is then
for some coefficients a=(a1, …, aN) to be determined by the method of A.G.V. et al.22 (see the Supplementary Materials). The ground state in the non-interaction limit is the single-particle harmonic eigenstate ψ0(x1) multiplied by the Slater determinant for the remaining particles in the states ψn(xk) for k=2, …, N and n=0, …, N−2.
For systems of four and five particles, the interpolatory ansatz with these |γ0 and |γ∞ yields integrals that can be evaluated analytically. For larger systems, the integrals are readily evaluated numerically. The resulting energies for N=4–6 are plotted as dashed lines in Fig. 5. Here we see a very good agreement between the ansatz and the numerically exact results, and in turn excellent agreement with experimental measurements19. This indicates that the energetics of the system is captured by the interpolatory ansatz with high accuracy.
However, as we discussed briefly above, the ansatz does not generally reproduce the correct first-order energy term at 1/g=0. Assuming that the non-interacting state should remain untouched, this prompts us to investigate whether another state can be found in the 1/g=0 limit that can replace |γ∞ in Eq. (4) and in turn give the exact result for the slope of the energy at 1/g=0. As noted above, Eq. (7) remains valid if we substitute the eigenstate |γ∞ with any other state with energy E∞ that obeys V|γ∞=0. In particular, any choice of a in Eq. (16) would work, and Eopt only depends on a through the wave-function overlap γ0|γ∞. Hence, because Eopt is monotonic in γ0|γ∞2, we may optimise the energy with respect to a for all values of g by maximising γ0|γ∞(a)2 with respect to a. Leaving the details to the Supplementary Materials, the optimum is
For the ground state of the N=3–6 systems, however, γ0|γ∞(amax)2 is very close to the known exact value of γ0|γ∞2, and is not large enough to make the slope of the energy correct in the strongly-interacting limit, that is, γ0|γ∞(amax)2<λ with λ given by Eq. (9). This indicates that we cannot find a state in the infinite-interaction limit that both has the correct zeroth-order energy, E∞, and satisfies the delta-boundary conditions, V|γ∞=0. However, we caution that this is under the assumptions that the wave function of |γ∞ is continuous and has discontinuous derivatives that satisfy the boundary conditions imposed by the zero-range interaction.
We see that the interpolatory ansatz–with whatever choice of |γ∞–does not reproduce the correct energy slope. The modified ansatz, however, does. It is worthwhile to note that the inability to find a state, |γ=α0|γ0+α∞|γ∞, whose expectation value is that predicted by the modified ansatz, does not imply that the modified ansatz cannot be used. The two methods rely on different information about the system: The interpolatory ansatz requires the knowledge of the wave function at g=0 and 1/g=0, whereas the modified ansatz uses the knowledge of the energy behaviour. Therefore, one can use the modified anzatz to estimate the energies, and the interpolatory ansatz to approximate the wave functions.
Figure 5 compares the ground-state energy of the modified ansatz with results from an exact numerical diagonalisation as well as experimental data19. Here we see an improved agreement for large g. The deviations of the modified ansatz from exact numerical diagonalisation are shown in greater detail in Fig. 6. The errors are more than an order of magnitude smaller than those of the unmodified ansatz (shown in Supplementary Fig. S4).
We plot the error scaled against N ↑ instead of E0 or E, because E0 scales as while ΔE=N↑. As seen in Fig. 5, the maximum in the error moves to larger interaction strength for higher N↑, i.e. with system size. There seems to be only very little increase in the magnitude of the error for larger N↑ for the system sizes we have studied. Notice that the error is slightly negative for 1/g0 in the case N=3–6. Since we know the exact energies and slopes around 1/g=0, the likely cause of this is that we are pushing the accuracy of the exact numerical diagonalisation method here.
Our first observation is that the approximation of the ground-state energy offered by the modified ansatz is so good that it can compete and even in some cases beat state-of-the-art numerical exact methods. The drawback seems quite evident, i.e. we cannot be sure that any state exists in the infinite-interaction limit that when used in Eq. (5) would reproduce the energy of the modified ansatz. Thus, the modified ansatz does not immediately give us any information about the wave function in spite of its near perfect approximation of the energy. We also note that the modified ansatz breaks the variational bound on the ground state and can in principle have lower energy as compared to the exact result. However, we have clearly demonstrated that the slope of the energy at 1/g=0 is an extremely important quantity for these systems as it appears crucial to reproduce in order to capture the energetics. This highlights the important role played by the recently developed theory of strongly interacting confined systems22.
Finally, we discuss the so-called Anderson overlap, which is the wave-function overlap between the non-interacting eigenstate, |γ0, and the interacting state, |γ, for some value of the interaction strength, g. This quantity is related to the Anderson orthogonality catastrophe32, which states that the Anderson overlap is zero in the thermodynamic limit; that is γ0|γ→0 for N→∞, and in particular γ0|γ∞→0 for N→∞.
In Fig. 7 we illustrate how the overlap γ0|γ2–with |γ given by the interpolatory ansatz–decreases as a function of g. Due to the fact that we only consider a finite-sized system, the overlap does not approach zero as g→∞, but the plot clearly shows that the overlap from our interpolatory ansatz tends to decrease as expected. The fact that the ansatz gives a very accurate approximation for the energy of the system does not immediately imply that this is also the case for the wave function. We leave this question for future studies, and the overlaps presented here are thus predictions based on the ansatz. It should be compared either to elaborate exact numerical calculations or to experimental measurements.
We have proposed a simple interpolatory ansatz for approximating the energy eigenstate of a confined, one-dimensional system of interacting particles. The ansatz is a linear combination of known eigenstates in the extreme limits of the interaction strength, g→0 and 1/g→0, respectively. Thanks to recent advances in the description of the eigenstates in the 1/g=0 limit, both these wave functions are now available. An analytical expression for the optimum energy of this ansatz is presented which is an elementary function of only two matrix elements; the interaction energy of the eigenstate at g=0, and the wave-function overlap between the eigenstates in the two limits of g→0 and 1/g→0. By focusing on harmonically trapped impurity systems of fermions, we have demonstrated that the ansatz is able to capture the physics of such a system. It gives us a highly accurate approximation for the energy and it also gives us a very simple expression for the wave function.
For both the two- and three-particle systems, we have been able to reproduce the entire energy spectrum with the interpolatory ansatz, save for deeply bound molecular states. The ansatz can be extended to describe deeply bound states as well; this has been shown specifically for the two-particle system. Taking the three-particle system as an example, we have also demonstrated that the ansatz works equally well for mass-imbalanced systems. A future extension of this study might investigate mass-imbalanced systems of four or more particles. It should be noted that the bottleneck here is that there are generally quite few known results about mass-imbalanced systems in the 1/g→0 limit33. It is an open problem to find a general method that yields exact eigenstates in the strongly interacting regime for mass-imbalanced systems.
A drawback of our ansatz is that it is in general not perturbatively correct in the strongly-interacting regime. More precisely, if we take the first order derivative of the energy with respect to 1/g it deviates from the known exact result. We may, however, modify the expression for the energy of the interpolatory ansatz slightly such that it is perturbatively correct to linear order in g for g→0 and to linear order in 1/g for 1/g→0. The modified ansatz has great simplicity and accuracy at a level that is competitive with state-of-the-art numerical methods for obtaining the energy of the ground state for arbitrary g. Due to its simplicity it should provide a very useful tool.
We note that although the results presented here assume that the particles are trapped in an external harmonic trap, the formalism is completely general and can be applied for arbitrary external potentials with at least N bound single-particle states for N-body systems. As we have shown, the relative deviation of the energy obtained from the ansatz grows only very slowly with N, and there should be no problem in extending the technique to even larger systems than considered here. The decisive quantity is the overlap of the g=0 and 1/g=0 wave functions. This overlap may be computed using the same methods that have recently been used to compute spin chain models for strongly interacting fermions34 and thus scaling to larger N of order 30 or 40 is certainly within reach.
A future direction would be to consider a generalised version of the interpolatory ansatz where one systematically adds more states at g=0 and 1/g=0 in order to gradually improve the comparison. In addition, it is relatively straightforward to apply the interpolatory ansatz in the case of strongly interacting bosons35,36, or mixed systems22,28. The requirements are knowledge of states in the two limits and their overlaps so that the interpolation can be performed. The formulas for the interpolated energy given here still apply. An example could be an impurity interacting strongly with a Tonks-Girardeau gas of hard-core bosons, which is a topic of great recent interest37,38,39.
In the following, we provide the details of the methods used in applying the interpolatory ansatz to two-, three- and many-particle systems.
We consider a system of two distinguishable fermions and we define . In the absence of interactions, the energy eigenstates are the harmonic eigenstates denoted |n with integer n≥0 (we are only concerned with the motion relative to the center of mass).
and the wave function of the state is
where 1F1 is the confluent hypergeometric function of the first kind.
Let |γ0=|n1 be an energy eigenstate in the non-interaction limit. Here |n1 is a single-particle eigenstate of a harmonic oscillator Hamiltonian (in the relative coordinate x) with quantum number n1. Furthermore, let |γ∞ be the corresponding eigenstate in the infinite-interaction limit–that is, |γ0→|γ∞ as the interaction strength is changed adiabatically from g=0 to 1/g=0. By conservation of parity, the states |γ0 and |γ∞ have the same parity.
If |γ0 is odd, x=0|γ0=0 and thus γ0|V|γ0=0. Hence, odd harmonic eigenstates do not change as we introduce a non-zero interaction. The even harmonic eigenstates do, however, change; the correct even eigenstate in the infinite-interaction limit is
with n2=n1+1 for repulsive interactions (g>0) and n2=n1−1 for attractive interactions (g<0).
Before we employ the interpolatory ansatz, we first separate out the center-of-mass motion using hyperspherical coordinates35,41. This is done merely for convenience and is not in any way essential for the approach. Defining and , the hyperradius is given by and the hyperangle is defined by tanϕ=y/x. The Hamiltonian of the relative motion becomes
where 2 is the Laplacian in polar coordinates (ρ, ϕ) and
is the interaction.
The quantum numbers ν and μ can be used to describe the energy eigenstates in both the non-interacting limit and the infinite-interaction limit. The eigenstate wave function in both limits has the general form41
where is a generalised Laguerre polynomial.
In the non-interacting limit, the (normalised) angular part of the wave function is
where p is the parity of the wave function. In this limit p=(−1)μ.
In the infinite-interaction limit, the wave function vanishes at the lines ϕ=−5π/6, −π/2, −π/6, π/6, π/2, 5π/6. In the regions between these lines, the wave function solves the Schrödinger wave equation for the harmonic oscillator without interactions.
The eigenenergies in the limit 1/g=0 is given by Eq. (11) with μ=3, 6, 9, 12, …. Each allowed eigenenergy is three-fold degenerate (not counting states with a different radial quantum number ν): One of the eigenstates in each energy triplet is a harmonic eigenstate with μ=3, 6, 9, 12, … which is unaffected by the interactions. Between the two remaining eigenstates in the triplet, one is odd and the other is even.
When μ is even (that is, μ=6, 12, 18, …) the non-trivial eigenstates in the infinite-interaction limit have the angular wave function
if p is their parity. The coefficients in the four remaining regions follow by symmetry considerations.
When μ is odd (μ=3, 9, 15, …),
The transformation into hyperspherical coordinates proceeds along the same lines with a modified Jacobian. Generally, the coordinates x=(x1, x2, x3) are transformed to Jacobi coordinates, x′=Jx, through the transformation matrix
where M23=m2+m3, M123=m1+m2+m3 and the ‘reduced’ masses are defined as and . This transformation allows us to separate the center-of-mass motion from the relative motion, the solutions of the former being the well-known harmonic eigenstates. Afterwards, we transform the remaining relative coordinates into hyperspherical coordinates, ρ and ϕ, by and tan(ϕ)=.
From now on, we assume that m1=M and m1=m2=m. The interaction potential can then be written as
The μ eigenvalue can be found by using parity symmetry, the Pauli principle and the delta-boundary conditions of the interaction potential. Using these conditions, one can setup an equation from which μ can be obtained in the limits g=0 and 1/g=0: For solutions with odd parity, μ solves the equation
For even-parity solutions, the equation is
The wave function in the infinite-interaction limit now vanishes along the delta-boundary lines ϕ=±θ0, ±π/2, π±θ0. Only when M=m, is the wave function non-zero in all six regions separated by the delta-boundary lines. In order to illustrate this, we look at two explicit examples where M/m=1/2 and M/m=2, respectively. For the case of M/m=1/2, we have θ00.421 (or 24.1°), and the lowest-energy solutions of Eqs (30 and 31) at 1/g=0 have μ0.731 (both odd and even) while the second-lowest has μ3.75 (odd). The angular part of the wave function for the ground state and the first excited state in the infinite-interaction limit are
where is the parity of the state. The angular part for the non-degenerate second excited state is
When M/m=2, the roles are reversed. The ground state is now non-degenerate with μ2.552 at g=∞ and θ00.615 (or 35.2°). In addition, the wave function has the same form as Eq. (33), only with different θ0 and μ. For the first and second excited states, the wave function has the form of Eq. (32).
One might think that there would be some continuous crossover from M<m to M=m and then to M>m, but this is not the case. Indeed, M=m is a singular case where the wave function is non-zero in all six regions.
How to cite this article: Andersen, M. E. S. et al. An interpolatory ansatz captures the physics of one-dimensional confined Fermi systems. Sci. Rep. 6, 28362; doi: 10.1038/srep28362 (2016).
Part of this work is based on the Bachelor thesis of M.E.S.A. The authors thank C. Forssén, D.V. Fedorov, and A.S. Jensen for discussions, as well as C. Forssén and X. Cui for feedback on the manuscript. We thank the experimental team of the S. Jochim group for extended discussions and for sharing their results. This work was funded by the Danish Council for Independent Research DFF Natural Sciences and the DFF Sapere Aude program.
Author Contributions A.G.V., E.J.L. and N.T.Z. devised the project. M.E.S.A. and A.S.D. developed the formalism under the supervision of A.G.V. and N.T.Z. The numerical calculations were carried out by A.S.D. and M.E.S.A. with assitance from A.G.V. and E.J.L. The initial draft of the paper was written by M.E.S.A., A.S.D. and N.T.Z. All authors contributed to the revisions that led to the final version.