“…Light propagating in nonlinear and nonlocal systems, under the paraxial approximation, is described by a Schrödinger-Newton system, which consists in a Schrödinger equation coupled with a Newtonian-like potential (or, more generally, a Poisson potential), where the former is responsible for describing the evolution of the envelope of the beam through the medium, and the last one for describing the distribution of the refractive index, and examples of materials where this model is used are nematic liquid crystals [8,9,12], thermo-optical materials [10,11] and quantum gases [13,14]. Indeed, the Schrödinger-Newton model is capable of describing a much wider set of systems, ranging from boson stars [10] to dark matter [11] and superfluidity [9,15], to name a few, and due to the similarities between these mathematical descriptions, this system has been proposed and used for implementing optical analogues [9][10][11]15], where systems that are hard or even impossible to study are emulated in the laboratory under controlled conditions. Over the years, many numerical models have been developed and improved to simulate this class of systems and in the last years, at our research group, we have developed a set of high-performance solvers based on GPGPU supercomputing to numerically study this class of systems, and successfully applied them in the study of superfluidity in nematic liquid crystals [9] and persistent currents in atomic gases [16].…”