“…Thirdly, in biased cells nematicons propagate in a wide linear index well and possess walk-off in the vertical plane (x, z): the model used in [8] does not account for this dynamics, which was addressed later in [11,15]. Reality is far more complex than a single Schrödinger-Poisson equation: (i) in real samples there are interfaces at finite distances, as we discussed in [1]; ii) several nonlinear effects act together, thus models considering only the reorientational response are approximations [50,51]; iii) actual reorientation is driven by a more complicated equation than a single Poisson equation, even in the perturbative regime [38,52]; iv) NLC described by a molecular director field are a simplification of an underlying many body system which is subject to continuous temporal fluctuations [38,53]. The elegant SMM, despite its mathematical simplicity, explains qualitatively nonlinear light propagation, predicting an oscillatory (breathing) behavior qualitatively different from the breather dynamics in a standard NLSE [54][55][56], stability in (2 + 1)D [57] and interaction between solitons [17,58,59].…”