A Mass- and Energy-Conserving Numerical Model for a Fractional Gross–Pitaevskii System in Multiple Dimensions

This manuscript studies a double fractional extended p-dimensional coupled Gross–Pitaevskii-type system. This system consists of two parabolic partial differential equations with equal interaction constants, coupling terms, and spatial derivatives of the Riesz type. Associated with the mathematical model, there are energy and non-negative mass functions which are conserved throughout time. Motivated by this fact, we propose a finite-difference discretization of the double fractional Gross–Pitaevskii system which inherits the energy and mass conservation properties. As the continuous model, the mass is a non-negative constant and the solutions are bounded under suitable numerical parameter assumptions. We prove rigorously the existence of solutions for any set of initial conditions. As in the continuous system, the discretization has a discrete Hamiltonian associated. The method is implicit, multi-consistent, stable and quadratically convergent. Finally, we implemented the scheme computationally to confirm the validity of the mass and energy conservation properties, obtaining satisfactory results.