Stationary transmission through lattices with asymmetric nonlinear quadratic-cubic defect (Q2166507)

The authors write down the transmission coefficients for a discrete stationary nonlinear Schrödinger equation and look at mechanisms that lead to different transmission coefficient behavior in different directions. A starting point for their work is a time dependent discrete nonlinear Schrödinger equation (Equation (1) in [the first author et al., Phys. Lett., A 447, Article ID 128301, 8 p. (2022; Zbl 07574940)]). From this equation, a discrete stationary nonlinear Schrödinger equation is derived (Equation (2) in [the first author et al., Phys. Lett., A 447, Article ID 128301, 8 p. (2022; Zbl 07574940)]) \[ \psi_{n-1}=-\psi_{n+1}+(-\omega+V_n + \gamma_n|\psi_n|^2-\rho_n|\psi_n|)\psi_n \] $n\in\mathbb{Z}$, $\omega=-2\cos(k)$, $0\leq k\leq\pi$. Boundary conditions are imposed (Equation (3) in [\textit{M. A. Wasay} et al., Phys. Lett., A 447, Article ID 128301, 8 p. (2022; Zbl 07574940)]). Recursion relations can be set up (Equation (3) in [\textit{S. Lepri} and \textit{G. Casati}, ``Asymmetric wave propagation in nonlinear systems'', Phys. Rev. Lett. 106, No. 16, Article ID 164101, 4 p. (2011; \url{doi:10.1103/PhysRevLett.106.164101})]). Two primary cases that they consider are when $V_n$, $\gamma_n$, $\rho_n$ are nonzero for $1\leq n\leq N$ and zero elsewhere. $N$ equals either $2$ or $3$. For these two cases, the boundary conditions and the recursion relations are used to find exact expressions for the transmission coefficients $t(k,|T|^2)$, $t(-k,|T|^2)$ (Equations (15) and (16) in [\textit{M. A. Wasay} et al., Phys. Lett., A 447, Article ID 128301, 8 p. (2022; Zbl 07574940)] and the paragraph after Equation (16) in [\textit{M. A. Wasay} et al., Phys. Lett., A 447, Article ID 128301, 8 p. (2022; Zbl 07574940)]). The complex number $T$ appears as part of the boundary conditions (Equation (3) in [\textit{M. A. Wasay} et al., Phys. Lett., A 447, Article ID 128301, 8 p. (2022; Zbl 07574940)]). One of the main findings of their work is that at special values of $k$, $V_n$, $\gamma_n$, $\rho_n$, the transmission coefficient $t(-k,|T|^2)$ as a function of $|T|^2$ has an additional local maximum compared to $t(k,|T|^2)$ for the $N=2$ case (Figure 1, left panel, in [\textit{M. A. Wasay} et al., Phys. Lett., A 447, Article ID 128301, 8 p. (2022; Zbl 07574940)]). Other differing properties of the two transmission coefficients are noted (Figures 3, 4, 5 and 6 in [\textit{M. A. Wasay} et al., Phys. Lett., A 447, Article ID 128301, 8 p. (2022; Zbl 07574940)]). They also perform numerical studies of the long-time behavior of solutions to the time dependent discrete nonlinear Schrödinger equation with Gaussian initial conditions (Equation (1), Figures 7, 8, Equation (18) and Table 1 in [\textit{M. A. Wasay} et al., Phys. Lett., A 447, Article ID 128301, 8 p. (2022; Zbl 07574940)]).