This paper presents a deep analysis of a time-dependent Schrödinger equation with fractional time derivative. After the discretization of the spatial operator the equation is reformulated in terms of a special system of fractional differential equations; it occurs that the eigenvalues of the coefficient matrix lay on the boundary of the stability region of fractional differential equations. The main difficulties in solving this system are hence related to the simultaneous presence of persisting oscillations (possibly with high frequency as it is typical with Schrödinger equations) and a persisting memory (as a consequence of the fractional order); moreover, an accurate spatial discretization gives rise to systems of large to very large size, involving a noteworthy computational complexity. By means of a theoretical analysis the exact solution is split into two or three terms (depending on the order of the fractional derivative), thus to face the numerical computation by different and suitably selected methods: direct evaluation of matrix functions for the terms characterized by smooth behaviour but with persistent memory and a step-by-step strategy, in conjunction with matrix function, for the oscillating term. In both cases, Krylov subspace methods are employed for the computation of matrix functions and convergence results are presented

On the time-fractional Schrödinger equation: Theoretical analysis and numerical solution by matrix Mittag-Leffler functions / Garrappa, Roberto; Moret, Igor; Popolizio, Marina. - In: COMPUTERS & MATHEMATICS WITH APPLICATIONS. - ISSN 0898-1221. - STAMPA. - 74:5(2017), pp. 977-992. [10.1016/j.camwa.2016.11.028]

On the time-fractional Schrödinger equation: Theoretical analysis and numerical solution by matrix Mittag-Leffler functions

Marina Popolizio
2017-01-01

Abstract

This paper presents a deep analysis of a time-dependent Schrödinger equation with fractional time derivative. After the discretization of the spatial operator the equation is reformulated in terms of a special system of fractional differential equations; it occurs that the eigenvalues of the coefficient matrix lay on the boundary of the stability region of fractional differential equations. The main difficulties in solving this system are hence related to the simultaneous presence of persisting oscillations (possibly with high frequency as it is typical with Schrödinger equations) and a persisting memory (as a consequence of the fractional order); moreover, an accurate spatial discretization gives rise to systems of large to very large size, involving a noteworthy computational complexity. By means of a theoretical analysis the exact solution is split into two or three terms (depending on the order of the fractional derivative), thus to face the numerical computation by different and suitably selected methods: direct evaluation of matrix functions for the terms characterized by smooth behaviour but with persistent memory and a step-by-step strategy, in conjunction with matrix function, for the oscillating term. In both cases, Krylov subspace methods are employed for the computation of matrix functions and convergence results are presented
2017
On the time-fractional Schrödinger equation: Theoretical analysis and numerical solution by matrix Mittag-Leffler functions / Garrappa, Roberto; Moret, Igor; Popolizio, Marina. - In: COMPUTERS & MATHEMATICS WITH APPLICATIONS. - ISSN 0898-1221. - STAMPA. - 74:5(2017), pp. 977-992. [10.1016/j.camwa.2016.11.028]
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11589/169254
Citazioni
  • Scopus 22
  • ???jsp.display-item.citation.isi??? 21
social impact