Abstract:As a typical quantum many body problem, we consider the time evolution of density matrix elements in the Bose-Hubbard model. For an arbitrary initial state, these quantities can be obtained from an SDE or stochastic differential equation system. To this SDE system, a Girsanov transformation can be applied. This has the effect that all the information from the initial state moves into the drift part, into the mean field part, of the transformed system. In the large N limit with fixed, the diffusive part of the transformed system vanishes and as a result, the exact quantum dynamics is given by an ODE system which turns out to be the time dependent discrete Gross Pitaevskii equation. For the two site Bose-Hubbard model, the GP equation reduces to the mathematical pendulum and the difference of expected number of particles at the two lattice sites is equal to the velocity of that pendulum which is either oscillatory or it can have rollovers which then corresponds to the self trapping or insulating phase. As a by-product, we also find an equivalence of the mathematical pendulum with a quartic double well potential. Collapse and revivals are a more subtle phenomenom, in order to see these the diffusive part of the SDE system or quantum corrections have to be taken into account. This can be done with an approximation and collapse and revivals can be reproduced, numerically and also through an analytic calculation. Since expectation values of Fresnel or Wiener diffusion processes, we write the density matrix elements exactly in this way, can be obtained from parabolic second order PDEs, we also obtain various exact PDE representations. The paper has been written with the goal to come up with an efficient calculation scheme for quantum many body systems and as such the formalism is generic and applies to arbitrary dimension, arbitrary hopping matrices and, with suitable adjustments, to fermionic models.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is about the time - evolution problem of density matrix elements in the Bose - Hubbard model. Specifically, the author studies the dynamic behavior of this quantum many - body problem through stochastic calculus and Girsanov transformation.
### Main problems of the paper:
1. **Description of time - evolution**: The paper considers the time - evolution of density matrix elements of any initial state in the Bose - Hubbard model. These quantities can be obtained through a system of stochastic differential equations (SDEs).
2. **Application of Girsanov transformation**: By applying the Girsanov transformation to this SDE system, the information of the initial state is transferred to the drift part (i.e., the mean - field part) of the transformed system. This makes the diffusion part disappear in the large - N limit (when \( g=\frac{U}{N} \) is fixed), thus obtaining an exact quantum - dynamic description.
3. **Discrete Gross - Pitaevskii equation**: In the large - N limit, the exact quantum dynamics of the system is given by a system of ordinary differential equations (ODEs), and this ODE system is actually the time - dependent discrete Gross - Pitaevskii equation.
4. **Explanation of physical phenomena**: For the two - site Bose - Hubbard model, the GP equation can be simplified to the mathematical pendulum model, and the particle - number difference corresponds to the speed of the pendulum. In addition, the collapse and revival phenomena are also explored, and these phenomena are reproduced by approximate methods.
### Formula summary:
- Hamiltonian of the Bose - Hubbard model:
\[
H=-J\sum_{\langle i,j\rangle}a_{i}^{\dagger}a_{j}+\frac{U}{2}\sum_{j}a_{j}^{\dagger}a_{j}^{\dagger}a_{j}a_{j}+\sum_{j}\epsilon_{j}a_{j}^{\dagger}a_{j}
\]
- Time - evolution operator:
\[
e^{-itH}=e^{-it(h_{0}+h_{\text{int}})}
\]
- Girsanov transformation in the SDE system:
\[
d\tilde{x}_{j,t}=dx_{j,t}-\sqrt{2u}dt(v_{j}\bar{v}_{j})
\]
\[
d\tilde{y}_{j,t}=dy_{j,t}-\sqrt{2u}dt(v_{j}\bar{v}_{j})
\]
- Time - evolution of density matrix elements:
\[
(\psi_{t},a_{i}^{\dagger}a_{j}\psi_{t})_{F}=E[\tilde{v}_{i}\bar{\tilde{v}}_{j}]
\]
### Summary:
By introducing stochastic calculus and Girsanov transformation, the paper provides an effective method for calculating the time - evolution of quantum many - body systems. This method is not only applicable to the Bose - Hubbard model, but also can be generalized to any dimension, any transition matrix, and fermion models.