Abstract. We study a non-Hermitian PT −symmetric generalization of an Nparticle, two-mode Bose-Hubbard system, modeling for example a Bose-Einstein condensate in a double well potential coupled to a continuum via a sink in one of the wells and a source in the other. The effect of the interplay between the particle interaction and the non-Hermiticity on characteristic features of the spectrum is analyzed drawing special attention to the occurrence and unfolding of exceptional points (EPs). We find that for vanishing particle interaction there are only two EPs of order N + 1 which under perturbation unfold either into [(N + 1)/2] eigenvalue pairs (and in case of N + 1 odd, into an additional zero-eigenvalue) or into eigenvalue triplets (third-order eigenvalue rings) and (N + 1) mod 3 single eigenvalues, depending on the direction of the perturbation in parameter space. This behavior is described analytically using perturbational techniques. More general EP unfoldings into eigenvalue rings up to (N + 1)th order are indicated.

We investigate an N-particle Bose-Hubbard dimer with an additional effective decay term in one of the sites. A mean-field approximation for this non-Hermitian many-particle system is derived, based on a coherent state approximation. The resulting nonlinear, non-Hermitian two-level dynamics, in particular, the fixed point structures showing characteristic modifications of the self-trapping transition, are analyzed. The mean-field dynamics is found to be in reasonable agreement with the full many-particle evolution.

We investigate the many-particle and mean-field correspondence for a non-Hermitian N-particle BoseHubbard dimer where a complex onsite energy describes an effective decay from one of the modes. Recently a generalized mean-field approximation for this non-Hermitian many-particle system yielding an alternative complex nonlinear Schrödinger equation was introduced. Here we give details of this mean-field approximation and show that the resulting dynamics can be expressed in a generalized canonical form that includes a metric gradient flow. The interplay of nonlinearity and non-Hermiticity introduces a qualitatively new behavior to the mean-field dynamics: The presence of the non-Hermiticity promotes the self-trapping transition, while damping the self-trapping oscillations, and the nonlinearity introduces a strong sensitivity to the initial conditions in the decay of the normalization. Here we present a complete characterization of the mean-field dynamics and the fixed point structure. We also investigate the full many-particle dynamics, which shows a rich variety of breakdown and revival as well as tunneling phenomena on top of the mean-field structure.