Abstract
Although the quantum theory of the optical response of individual atoms to coherent light with frequencies close to electronic transitions and the fluid equations for a gas are well known and understood from first principles, they are developed independently of each other and therefore cannot be applied directly to describe many of the quantum collective and transport phenomena that occur in cold atomic gases, especially in what regards their interaction with optical pulses and beams. Few attempts have been made to derive a consistent formalism and theory that are capable to model this type of systems, and those which exist rely on the adaptation of several ad-hoc hypothesis and simplifications, such as space and time dependent density operators. In this paper we provide the theoretical foundations and establish a formalism capable of paving the way for the development of new simulation tools and to explore new problems in nonlinear optics out of equilibrium.

Abstract
This paper reports the development of a numerical module for the HiLight simulation platform dedicated to the propagation of light in nonlocal nonlinear optical media and the adaptions implemented for it to be used as a numerical test-bed to evaluate the impact of new extensions of the Theory of General Relativity in the dynamics of a N-body system. The phenomenology of light in nonlocal and nonlinear media is very rich and can be described by a multitude of effective models, with different levels of detail and approximations, which coincide with few or no differences with those found in many other fields of physics. In particular, nonlocal extensions of the Generalized Nonlinear Schrodinger equation (also known as the generalized Schrodinger-Newton system) constitute a wide class of physical models that can be found in both optics and in the studies of alternative theories for gravity. Therefore, numerical solvers developed for the former can be adapted to address the later. Indeed, this paper discusses the adaptation of a numerical solver of the generalized Schrodinger-Newton system based on GPGPU supercomputing, initially developed to investigate the properties of light in exotic nonlocal media, to tackle the dynamics of large distributions of matter whose interaction is governed by extensions of the Theory of General Relativity, namely those based on non-minimal coupling between curvature and matter. This paper analysis the structure of the resulting simulation module, its performance and validation tests.

Abstract
This paper proposes the use of nematic liquid crystals as tunable setups to implement optical analogues of physical systems and phenomena that are hard or even impossible to study experimentally under controlled conditions. Optical analogues share the same physical model with the systems that they emulate and can be understood as a form of physical simulations or optical computation. However, their success relies not only on the existence of media with optical properties capable of emulating the models associated with the original system as they interact with light, but also on the possibility of being able to tune those properties in order to cover the multitude of conditions or range of parameters. In particular, the Schrodinger-Newton model is a good target for this kind of studies as it can describe a plethora of different phenomena in physics and can be implemented in the laboratory using optical analogues, usually using thermo-optical materials. However, such materials have limitations, and in this work we propose nematic liquid crystals as a more advantageous alternative. We discuss how nematic liquid crystals can be used as a tunable support medium for optical analogues of superfluids by analyzing the dispersion relation of light under specific conditions and using numerical simulations based on GPGPU supercomputing to verify our findings. Extending on this, we explore more direct manifestations of superfluid effects in nematic liquid crystals, such as drag-force cancellation in the superfluid regime and the possibility of creating a roton-minimum in the dispersion relation.

Abstract
This paper explores and compares three different plasmonic optical fibre sensor configurations, based on D-type and suspended core fibres combined with metallic nanowires, and investigates how their different geometrical parameters can affect the coupling between the guided optical mode supported by fibres and the localized plasmonic modes, and how that ultimately results in improved sensor performance. Fibre optical sensors based on plasmonic resonances with metallic nanostructures have revolutionized the field of optical sensing because they have permitted to obtain sharper and fine-tuned resonances with higher sensitivity. The essence for exploring the properties of localized plasmonic modes and their coupling with the optical guided mode depends not only on the choice of the materials employed in the device, but also on the geometry of the different components and their relative position, which ultimately determines the spatial distributions of optical power of the different modes and consequently their overlap and coupling. In this work, we use numerical simulations based on finite element methods to demonstrate the importance of shaping the features of the guided optical mode to promote the coupling with the localized modes, in the two types of fibres considered. The results clarify some of the fundamental aspects behind the operation of these devices and provide novel proposals for enhanced refractive index sensors.

Abstract
The introduction of metallic nanostructures in optical fibers has revolutionized the field of plasmonic sensors since they produce sharper and fine-tuned resonances resulting in higher sensitivities and resolutions. This article evaluates the performance of three different plasmonic optical fiber sensors based on D-type and suspended core fibers with metallic nanowires. It addresses how their different materials, geometry of the components, and their relative position can influence the coupling between the localized plasmonic modes and the guided optical mode. It also evaluates how that affects the spatial distributions of optical power of the different modes and consequently their overlap and coupling, which ultimately impacts the sensor performance. In this work, we use numerical simulations based on finite element methods to validate the importance of tailoring the features of the guided optical mode to promote an enhanced coupling with the localized modes. The results in terms of sensitivity and resolution demonstrate the advantages of using suspended core fibers with metallic nanowires.

Abstract
We report on the development of numerical module for the HiLight simulation platform based on GPGPU supercomputing to solve a system of coupled fields governed by the Generalized Nonlinear Schrodinger Equation with local and/or nonlocal nonlinearities. This models plays an important role in describing a plethora of different phenomena in various areas of physics. In optics, this model was initially used to describe the propagation of light through local and/or nonlocal systems under the paraxial approximation, but more recently it has been extensively used as a support model to develop optical analogues. However, establishing the relation between the original system and the analogue, as well as, between their model and the actual experimental setup is not an easy task. First and foremost because in most cases the governing equations are nonintegrable, preventing from obtaining analytical solutions and hindering the optimization of the experiments. Alternatively, despite numerical methods not providing exact solutions, they allow to test different experimental scenarios and provide a better insight to what to expect in an actual experiment, while giving access to all the variables of the optical system being simulated. However, the numerical solution of a system of N-coupled Schrodinger fields in systems with two or three spatial dimensions requires massive computation resources, and must employ advanced supercomputing and parallelization techniques, such as GPGPU. This paper focuses on the numerical aspects behind this challenge, describing the structure of our simulation module, its performance and the tests performed.