“…This scaling will be used throughout this paper, and in particular in all the figures. To obtain numerical values for the overlap integrals γ ij , where i, j = ±, we solved (5) using the DVR method [33,34] with increasingly finer mesh, with iterations for each mesh to make the Φ ± functions and the the nonlinear term self-consisent. Values for the γ ij are shown as a function of barrier height, V b , for σ=1.5, for gN =1, 10 and 100, in Fig.…”