Nonlinear Schrodinger equation for a PT-symmetric delta-function double well

The time-independent nonlinear Schrodinger equation is solved for two attractive delta-function-shaped potential wells where an imaginary loss term is added in one well, and a gain term of the same size but with opposite sign in the other. We show that for vanishing nonlinearity the model captures all the features known from studies of parity-time-(PT-) symmetric optical waveguides, e. g., the coalescence of modes at an exceptional point at a critical value of the loss/gain parameter, and the breaking of PT symmetry beyond. With the nonlinearity present, the equation is a model for a Bose-Einstein condensate with loss and gain in a double-well potential. We find that the nonlinear Hamiltonian picks as stationary eigenstates exactly such solutions which render the nonlinear Hamiltonian itself PT symmetric, but observe coalescence and bifurcation scenarios different from those known from linear PT-symmetric Hamiltonians.