The work required to solve for the fully interacting N boson wavefunction, which is widely believed to scale exponentially with N, has been previously shown to scale as N0 when the problem is rearranged using analytic building blocks. The exponential complexity reappears in an exponential scaling with the order of our perturbation series allowing exact analytical calculations for very large N systems through low order. In this paper, we analyse the growth in complexity with order when a normal mode basis is used.
The resources required to solve the general interacting quantum
N-body problem scale exponentially with
N, making the solution of this problem very difficult when
N is large. In a previous series of papers we develop an approach for a fully interacting wave function for a confined system of identical bosons with a general two-body interaction. This method tames the
N scaling by developing a perturbation series that is order-by-order invariant under a point group isomorphic with
SN. Group theory and graphical techniques are then used to solve for the wave function exactly and analytically at each order, yielding a solution for the general
N-body problem which scales as
N0 at any given order. Recently this formalism has been used to obtain the first-order fully interacting wave function for a system of harmonically confined bosons interacting harmonically. In this paper, we report the application of this
N-body wave function to a system of
N fully interacting bosons in three dimensions. We derive an expression for the density profile for a confined system of harmonically interacting bosons. Choosing this simple interaction is not necessary or even advantageous for our method, however this choice allows a direct comparison of our exact results through first order with exact results obtained in an independent solution. Our density profile for the wave function through first order in three dimensions is indistinguishable from the first-order exact result obtained independently and shows strong convergence to the exact result to all orders.
The resources required to solve the general interacting quantum N-body problem scale exponentially with N, making the solution of this problem very difficult when N is large. In a previous series of papers we develop an approach for a fully-interacting wave function with a general two-body interaction which tames the N-scaling by developing a perturbation series that is order-by-order invariant under a point group isomorphic with S_N . Group theory and graphical techniques are then used to solve for the wave function exactly and analytically at each order. Recently this formalism has been used to obtain the first-order, fully-interacting wave function for a system of harmonically-confined bosons interacting harmonically. In this paper, we report the first application of this N-body wave function to a system of N fully-interacting bosons in three dimensions. We determine the density profile for a confined system of harmonically-interacting bosons. Choosing this simple interaction is not necessary or even advantageous for our method, however this choice allows a direct comparison of our exact results through first order with exact results obtained in an independent solution. Our density profile through first-order in three dimensions is indistinguishable from the first-order exact result obtained independently and shows strong convergence to the exact result to all orders.
Group theoretic and graphical techniques are used to derive the N-body wave function for a system of identical bosons with general interactions through first-order in a perturbation approach. This method is based on the maximal symmetry present at lowest order in a perturbation series in inverse spatial dimensions. The symmetric structure at lowest order has a point group isomorphic with the SN group, the symmetric group of N particles, and the resulting perturbation expansion of the Hamiltonian is order-by-order invariant under the permutations of the SN group. This invariance under SN imposes severe symmetry requirements on the tensor blocks needed at each order in the perturbation series. We show here that these blocks can be decomposed into a basis of binary tensors invariant under SN. This basis is small (25 terms at first order in the wave function), independent of N, and is derived using graphical techniques. This checks the N6 scaling of these terms at first order by effectively separating the N scaling problem away from the rest of the physics. The transformation of each binary tensor to the final normal coordinate basis requires the derivation of Clebsch–Gordon coefficients of SN for arbitrary N. This has been accomplished using the group theory of the symmetric group. This achievement results in an analytic solution for the wave function, exact through first order, that scales as N0, effectively circumventing intensive numerical work. This solution can be systematically improved with further analytic work by going to yet higher orders in the perturbation series.