Skin modes in a nonlinear Hatano-Nelson model

Non-Hermitian lattices with nonreciprocal couplings under open boundary conditions are known to possess linear modes exponentially localized on one edge of the chain. This phenomenon, dubbed non-Hermitian skin effect, induces all input waves in the linearized limit of the system to unidirectionally propagate toward the system's preferred boundary. Here we investigate the fate of the non-Hermitian skin effect in the presence of Kerr-type nonlinearity within the well-established Hatano-Nelson lattice model. Our method is to probe the presence of nonlinear stationary modes which are localized at the favored edge, when the Hatano-Nelson model deviates from the linear regime. Based on perturbation theory, we show that families of nonlinear skin modes emerge from the linear ones at any nonreciprocal strength. Our findings reveal that, in the case of focusing nonlinearity, these families of nonlinear skin modes tend to exhibit enhanced localization, bridging the gap between weakly nonlinear modes and the highly nonlinear states (discrete solitons) when approaching the anti-continuum limit with vanishing coupling. Conversely, for defocusing nonlinearity, these nonlinear skin modes tend to become more extended than their linear counterpart. To assess the stability of these solutions, we conduct a linear stability analysis across the entire spectrum of obtained nonlinear modes and also explore representative examples of their evolution dynamics.