3,635 research outputs found

    Solitary waves in the Nonlinear Dirac Equation

    Get PDF
    In the present work, we consider the existence, stability, and dynamics of solitary waves in the nonlinear Dirac equation. We start by introducing the Soler model of self-interacting spinors, and discuss its localized waveforms in one, two, and three spatial dimensions and the equations they satisfy. We present the associated explicit solutions in one dimension and numerically obtain their analogues in higher dimensions. The stability is subsequently discussed from a theoretical perspective and then complemented with numerical computations. Finally, the dynamics of the solutions is explored and compared to its non-relativistic analogue, which is the nonlinear Schr{\"o}dinger equation. A few special topics are also explored, including the discrete variant of the nonlinear Dirac equation and its solitary wave properties, as well as the PT-symmetric variant of the model

    Uncertainty quantification for kinetic models in socio-economic and life sciences

    Full text link
    Kinetic equations play a major rule in modeling large systems of interacting particles. Recently the legacy of classical kinetic theory found novel applications in socio-economic and life sciences, where processes characterized by large groups of agents exhibit spontaneous emergence of social structures. Well-known examples are the formation of clusters in opinion dynamics, the appearance of inequalities in wealth distributions, flocking and milling behaviors in swarming models, synchronization phenomena in biological systems and lane formation in pedestrian traffic. The construction of kinetic models describing the above processes, however, has to face the difficulty of the lack of fundamental principles since physical forces are replaced by empirical social forces. These empirical forces are typically constructed with the aim to reproduce qualitatively the observed system behaviors, like the emergence of social structures, and are at best known in terms of statistical information of the modeling parameters. For this reason the presence of random inputs characterizing the parameters uncertainty should be considered as an essential feature in the modeling process. In this survey we introduce several examples of such kinetic models, that are mathematically described by nonlinear Vlasov and Fokker--Planck equations, and present different numerical approaches for uncertainty quantification which preserve the main features of the kinetic solution.Comment: To appear in "Uncertainty Quantification for Hyperbolic and Kinetic Equations

    The low dimensional dynamical system approach in General Relativity: an example

    Get PDF
    In this paper we explore one of the most important features of the Galerkin method, which is to achieve high accuracy with a relatively modest computational effort, in the dynamics of Robinson-Trautman spacetimes.Comment: 7 pages, 5 figure

    A system of ODEs for a Perturbation of a Minimal Mass Soliton

    Full text link
    We study soliton solutions to a nonlinear Schrodinger equation with a saturated nonlinearity. Such nonlinearities are known to possess minimal mass soliton solutions. We consider a small perturbation of a minimal mass soliton, and identify a system of ODEs similar to those from Comech and Pelinovsky (2003), which model the behavior of the perturbation for short times. We then provide numerical evidence that under this system of ODEs there are two possible dynamical outcomes, which is in accord with the conclusions of Pelinovsky, Afanasjev, and Kivshar (1996). For initial data which supports a soliton structure, a generic initial perturbation oscillates around the stable family of solitons. For initial data which is expected to disperse, the finite dimensional dynamics follow the unstable portion of the soliton curve.Comment: Minor edit

    A time-splitting pseudospectral method for the solution of the Gross-Pitaevskii equations using spherical harmonics with generalised-Laguerre basis functions

    Get PDF
    We present a method for numerically solving a Gross-Pitaevskii system of equations with a harmonic and a toroidal external potential that governs the dynamics of one- and two-component Bose-Einstein condensates. The method we develop maintains spectral accuracy by employing Fourier or spherical harmonics in the angular coordinates combined with generalised-Laguerre basis functions in the radial direction. Using an error analysis, we show that the method presented leads to more accurate results than one based on a sine transform in the radial direction when combined with a time-splitting method for integrating the equations forward in time. In contrast to a number of previous studies, no assumptions of radial or cylindrical symmetry is assumed allowing the method to be applied to 2D and 3D time-dependent simulations. This is accomplished by developing an efficient algorithm that accurately performs the generalised-Laguerre transforms of rotating Bose-Einstein condensates for different orders of the Laguerre polynomials. Using this spatial discretisation together with a second order Strang time-splitting method, we illustrate the scheme on a number of 2D and 3D computations of the ground state of a non-rotating and rotating condensate. Comparisons between previously derived theoretical results for these ground state solutions and our numerical computations show excellent agreement for these benchmark problems. The method is further applied to simulate a number of time-dependent problems including the Kelvin-Helmholtz instability in a two-component rotating condensate and the motion of quantised vortices in a 3D condensate

    Structure preserving schemes for the continuum Kuramoto model: phase transitions

    Get PDF
    The construction of numerical schemes for the Kuramoto model is challenging due to the structural properties of the system which are essential in order to capture the correct physical behavior, like the description of stationary states and phase transitions. Additional difficulties are represented by the high dimensionality of the problem in presence of multiple frequencies. In this paper, we develop numerical methods which are capable to preserve these structural properties of the Kuramoto equation in the presence of diffusion and to solve efficiently the multiple frequencies case. The novel schemes are then used to numerically investigate the phase transitions in the case of identical and non identical oscillators
    • …
    corecore