Asaad R. Sakhel,Robert J. Ragan,William J. Mullin
Abstract:The Gross-Pitaevskii equation (GPE) in a double well potential produces solutions that break the symmetry of the underlying non-interacting Hamiltonian, i.e., asymmetric solutions. The GPE is derived from the more general second-quantized Fock Schroedinger equation (FSE). We investigate whether such solutions appear in the more general case or are artifacts of the GPE. We use two-mode analyses for a variational treatment of the GPE and to treat the Fock equation. An exact diagonalization of the FSE in dual-condensates yields degenerate ground states that are very accurately fitted by phase-state representations of the degenerate asymmetric states found in the GPE. The superposition of degenerate asymmetrical states forms a cat state. An alternative form of cat state results from a change of the two-mode basis set.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to evaluate whether the solutions of the Gross - Pitaevskii equation (GPE) in the double - well potential can accurately describe the asymmetric quantum condensates resulting from nonlinear bifurcations. Specifically, the authors explored whether these asymmetric solutions generated by GPE also appear in the more general Fock Schrödinger equation (FSE), or are merely a by - product of GPE.
### Background and Problem
The Gross - Pitaevskii equation (GPE) is a single - particle functional equation for describing Bose - Einstein condensates (BEC), and it has been very successful in dealing with problems in atomic traps. Especially in the double - well potential, GPE can exhibit various phenomena, such as symmetry breaking, self - trapping, etc. However, GPE is derived from the more general second - quantized Fock Hamiltonian. By replacing the field operator with a classical function, it ignores the small effect of the non - commutativity of the field operator. Therefore, a natural question is: Do the various effects observed in GPE also appear in the original Fock Schrödinger equation (FSE), or are they merely by - products of GPE?
### Research Methods
To answer this question, the authors used two mode - analysis methods:
1. **Variational treatment of GPE**: By choosing the states of the two lowest non - interacting single - particle Hamiltonians as variational trial wave functions, the solutions of GPE are studied.
2. **Variational treatment of FSE**: Also use the two - mode approximation method to deal with the Fock Hamiltonian to study the solutions of FSE.
### Main Findings
1. **Results of GPE**:
- When the interaction parameter \(\eta N\) is less than a certain critical value, the energy of the symmetric state is higher than that of the asymmetric state, resulting in the formation of the asymmetric state.
- The energy of the asymmetric state obtained by variational treatment is very close to the exact solution.
2. **Results of FSE**:
- In FSE, when the interaction parameter \(\eta N\) is less than a certain critical value, double condensates will also appear, and these states are asymmetric.
- Through the exact diagonalization including the tunneling term, the solutions of FSE show asymmetric states similar to those of GPE.
### Conclusion
The main conclusion of the paper is that GPE has high accuracy in describing the asymmetric quantum condensates resulting from nonlinear bifurcations in the double - well potential. These asymmetric states not only appear in GPE, but also in the more general Fock Schrödinger equation, indicating that these phenomena are inherent in the physical system, rather than specific by - products of GPE.
### Formula Examples
- **Time - dependent form of GPE**:
\[
i\frac{\partial \Phi}{\partial t}=\left(-\frac{\partial^{2}}{\partial x^{2}}+V_{\text{ex}}(x)+\eta|\Phi|^{2}\right)\Phi
\]
- **One - dimensional form of Fock Hamiltonian**:
\[
H = \int dx\,\Psi^{\dagger}(x)\left(-\frac{\partial^{2}}{\partial x^{2}}+V_{\text{ex}}(x)\right)\Psi(x)+\frac{1}{2}\int dxdx'\,\Psi^{\dagger}(x)\Psi^{\dagger}(x')W(x - x')\Psi(x')\Psi(x)
\]
- **Diagonal energy part**:
\[
H_0=\epsilon_0\frac{n_0}{N}+\epsilon_1\frac{n_1}{N}+\frac{\eta N}{2}\left(\chi_{40}\frac{n_0(n_0 - 1)}{N^2}+\chi_{04}\frac{n_1(n_1 - 1)}{N^2}+4\ch