We use the Bose-Hubbard Hamiltonian to study quantum fluctuations in canonical equilibrium ensembles of bosonic Josephson junctions at relatively high temperatures, comparing the results for finite particle numbers to the classical limit that is atta
ined as $N$ approaches infinity. We consider both attractive and repulsive atom-atom interactions, with especial focus on the behavior near the T=0 quantum phase transition that occurs, for large enough $N$, when attractive interactions surpass a critical level. Differences between Bose-Hubbard results for small $N$ and those of the classical limit are quite small even when $N sim 100$, with deviations from the limit diminishing as 1/N.
This Mathematica 7.0/8.0 package upgrades and extends the quantum computer simulation code called QDENSITY. Use of the density matrix was emphasized in QDENSITY, although that code was also applicable to a quantum state description. In the present ve
rsion, the quantum state version is stressed and made amenable to future extensions to parallel computer simulations. The add-on QCWAVE extends QDENSITY in several ways. The first way is to describe the action of one, two and three- qubit quantum gates as a set of small ($2 times 2, 4times 4$ or $8times 8$) matrices acting on the $2^{n_q}$ amplitudes for a system of $n_q$ qubits. This procedure was described in our parallel computer simulation QCMPI and is reviewed here. The advantage is that smaller storage demands are made, without loss of speed, and that the procedure can take advantage of message passing interface (MPI) techniques, which will hopefully be generally available in future Mathemati
