Current Organizations / Institutes

The organizations/institutes from the scholar's latest article


Athanassios S. Fokas has a total of 99 co-authors. Here are the top co-authors by number of shared publications:


Latest Publications

  • Open Access

    From the Radon Transform to Deep Learning in Medical Image Reconstruction

    by Nicholas E. Protonotarios, Sotirios Messinis, George A. Kastis, Nikolaos Dikaios, Athanassios S. Fokas

    Medical image reconstruction began with an exact mathematical result. In 1917, Radon showed how to recover a function from its line integrals, providing the mathematical basis for computed tomography, positron emission tomography, and single-photon emission computed tomography. This review follows the development of reconstruction from analytical inversion to modern deep learning. We first present the Radon transform, its attenuated and three-dimensional generalizations, and the corresponding inversion formulae, including filtered backprojection. We then review iterative and optimization-based methods, including maximum-likelihood expectation-maximization, variational regularization, and compressed sensing. The discussion proceeds to deep learning, covering convolutional and transformer architectures, algorithm unrolling, physics-informed networks, neural operators, and generative diffusion priors. Our central argument is that analytical, iterative, and learned reconstruction are different formulations of the same inverse problem rather than competing paradigms. Many learned components derive from classical formulations, and many classical operators admit learnable counterparts. We conclude by discussing open challenges related to theoretical guarantees, trustworthiness, data efficiency, federated learning, and regulatory validation. The most reliable learned systems are likely to retain the mathematical and physical structure of established reconstruction methods.

  • Open Access

    A Preliminary Analysis of a Physics-Informed Neural Network for the Forward Problem in EEG

    by Athanassios S. Fokas, Alireza Afzal Aghaei, Parham Hashemzadeh

    The distributed inverse source problem in electroencephalography (EEG) requires the determination of a current-independent, geometry-dependent auxiliary function, which is defined by a Poisson partial differential equation (PDE), where its solution is referred to as the forward problem. In this study, we investigate the feasibility of employing a mesh-free Physics-Informed Neural Network (PINN) for obtaining this auxiliary function. The proposed architecture integrates Kolmogorov–Arnold Networks (KANs) into an extended PINN (XPINN) framework augmented with Multi-scale Fourier feature mappings to capture potential field discontinuities across piecewise-homogeneous tissue interfaces. The PINN loss functional incorporates the governing PDE, Neumann boundary conditions, flux continuity and reference data for specific neuronal source and electrode configurations. Numerical experiments on a three-layer spherical head model demonstrate that the XPIKAN surrogate achieves a relative L2 error below 1% on unseen sensor coordinates. Factorial sensitivity analyses confirm stable model generalization across varying source-sensor configurations without the need for dense volumetric meshes. As a result, XPIKAN provides a meshless, continuous, and differentiable solution that offers faster inference time compared to classical solvers like finite element or boundary element methods and enables exact gradient computation for inverse source localization.

  • Open Access

    A Novel Asymptotic Technique for Integrals Involving the Hankel Contour and the Bleistein Asymptotic Formula

    by Athanassios S. Fokas, Jonatan Lenells

    Several important functions, including the gamma function, as well as several infinite sums, admit integral representations involving the Hankel contour. In addition, the large t asymptotic analysis of several recently derived identities satisfied by the Riemann zeta function requires the computation of the asymptotic form of certain integrals which also involve the Hankel contour; these integrals depend on a real parameter, α. A rigorous asymptotic technique is presented here for computing such integrals to all orders. For certain values of α, the relevant formula, in addition to an asymptotic series of explicit terms, also contains a specific integral. It is shown that, remarkably, the leading behavior of this integral can be written in the form of the leading order of the Bleistein integral. The latter integral arises in the implementation of the classical steepest descent method in the case that the stationary point coincides with one of the boundary points of the integral under consideration.

  • Open Access

    An explicit relationship between the error terms in Atkinson’s formula and in the formula for the Dirichlet divisor problem

    by Athanassios Fokas

    A simple formula is derived for the difference of the functions defining the error terms in Atkinson’s formula and in the formula for the Dirichlet divisor problem. Namely, this difference equals π/2 plus a function that is simply related to the square of the Riemann zeta function. Specifically, the large-t leading-order asymptotics of this function is expressed in terms of a double sum that differs from the large-t asymptotics of (ζ(1/2+it))2 only in the factor of 1/ln⁡(2πmn/t). Interestingly, to the leading order for large t, the derivative of the above function coincides with the square of the absolute value of the Riemann zeta function. This remarkable fact implies a very simple formula for the derivative of the function defining the error term in the Dirichlet divisor problem.

  • Higher-order diffusion and Cahn–Hilliard-type models revisited on the half-line

    by Andreas Chatziafratis, Alain Miranville, Georgia Karali, Athanasios S. Fokas, Elias C. Aifantis

    In this paper, we solve explicitly and analyze rigorously inhomogeneous initial-boundary-value problems (IBVP) for several fourth-order variations of the traditional diffusion equation and the associated linearized Cahn–Hilliard (C-H) model (also Kuramoto–Sivashinsky equation), formulated in the spatiotemporal quarter-plane. Such models are of relevance to heat-mass transfer phenomena, solid-fluid dynamics and the applied sciences. In particular, we derive formally effective solution representations, justifying a posteriori their validity. This includes the reconstruction of the prescribed initial and boundary data, which requires careful analysis of the various integral terms appearing in the formulae, proving that they converge in a strictly defined sense. In each IBVP, the novel formula is utilized to rigorously deduce the solution’s regularity and asymptotic properties near the boundaries of the domain, including uniform convergence, eventual (long-time) periodicity under (eventually) periodic boundary conditions, and null noncontrollability. Importantly, this analysis is indispensable for exploring the (non)uniqueness of the problem’s solution and a new counter-example is constructed. Our work is based on the synergy between: (i) the well-known Fokas unified transform method and (ii) a new approach recently introduced for the rigorous analysis of the Fokas method and for investigating qualitative properties of linear evolution partial differential equations (PDE) on semi-infinite strips. Since only up to third-order evolution PDE have been investigated within this novel framework to date, we present our analysis and results in an illustrative manner and in order of progressively greater complexity, for the convenience of readers. The solution formulae established herein are expected to find utility in well-posedness and asymptotics studies for nonlinear counterparts too.


  • Open Access

    Conscious and unconscious processes in vision and homeostasis

    by Athanassios S. Fokas, Nikos K. Logothetis

    The central and autonomic communications affecting cognition, emotions, and visceral functions, are determined by the interaction of unconscious and conscious processes. In this regard, we discuss two basic hypotheses. First, unconscious and conscious processes form a dynamic organizational continuum. Second, emotions, which are unconscious forms of feeling, characterize the response of an organism to any internal change disturbing homeostasis or to any change in the exterior of the organism detected via specialized sensory probes. The former hypothesis is illustrated by discussing aspects of visual perception. The validity of the second hypothesis is supported by discussing interactions between unconscious and conscious process necessary for maintaining energy and water balance.


  • Open Access

    Towards a unified approach to sciences and the arts

    by Athanassios S. Fokas

    The process of ‘comprehending how we comprehend’, requires elucidating fundamental cognitive mechanisms used by the brain. In these regards, this article discusses two basic hypotheses. First, unconscious and conscious processes form a continuum. Second, humans have a predisposition to construct real versions of mental images and of unconscious structures, as well as to assign to them specific symbols. In addition, it is claimed that our process of comprehending is aided by invoking certain concepts which are important because, apparently, they reflect fundamental cognitive mechanisms. Among them are the notions of continuity, associations, abstraction, reduction, unification, and generalization, as well as the dialectic pairs of local versus global processes and simplicity versus complexity.

  • Top Cited Papers
    Open Access

    Can artificial intelligence reach human thought?

    by Athanassios S Fokas

    The transformative achievements of deep learning have led several scholars to raise the question of whether artificial intelligence (AI) can reach and then surpass the level of human thought. Here, after addressing methodological problems regarding the possible answer to this question, it is argued that the definition of intelligence proposed by proponents of the AI as “the ability to accomplish complex goals,” is appropriate for machines but does not capture the essence of human thought. After discussing the differences regarding understanding between machines and the brain, as well as the importance of subjective experiences, it is emphasized that most proponents of the eventual superiority of AI ignore the importance of the body proper on the brain, the laterization of the brain, and the vital role of the glia cells. By appealing to the incompleteness theorem of Gödel’s and to the analogous result of Turing regarding computations, it is noted that consciousness is much richer than both mathematics and computations. Finally, and perhaps most importantly, it is stressed that artificial algorithms attempt to mimic only the conscious function of parts of the cerebral cortex, ignoring the fact that, not only every conscious experience is preceded by an unconscious process but also that the passage from the unconscious to consciousness is accompanied by loss of information.


  • Open Access

    Chebyshev Interpolation Using Almost Equally Spaced Points and Applications in Emission Tomography

    by Vangelis Marinakis, Athanassios S. Fokas, George A. Kastis, Nicholas E. Protonotarios

    Since their introduction, Chebyshev polynomials of the first kind have been extensively investigated, especially in the context of approximation and interpolation. Although standard interpolation methods usually employ equally spaced points, this is not the case in Chebyshev interpolation. Instead of equally spaced points along a line, Chebyshev interpolation involves the roots of Chebyshev polynomials, known as Chebyshev nodes, corresponding to equally spaced points along the unit semicircle. By reviewing prior research on the applications of Chebyshev interpolation, it becomes apparent that this interpolation is rather impractical for medical imaging. Especially in clinical positron emission tomography (PET) and in single-photon emission computerized tomography (SPECT), the so-called sinogram is always calculated at equally spaced points, since the detectors are almost always uniformly distributed. We have been able to overcome this difficulty as follows. Suppose that the function to be interpolated has compact support and is known at q equally spaced points in −1,1. We extend the domain to −a,a, a>1, and select a sufficiently large value of a, such that exactlyq Chebyshev nodes are included in −1,1, which are almost equally spaced. This construction provides a generalization of the concept of standard Chebyshev interpolation to almost equally spaced points. Our preliminary results indicate that our modification of the Chebyshev method provides comparable, or, in several cases including Runge’s phenomenon, superior interpolation over the standard Chebyshev interpolation. In terms of the L∞ norm of the interpolation error, a decrease of up to 75% was observed. Furthermore, our approach opens the way for using Chebyshev polynomials in the solution of the inverse problems arising in PET and SPECT image reconstruction.