The relative stability of two-dimensional soft quasicrystals in systems with two length scales is examined using a recently developed projection method, which provides a unified numerical framework to compute the free energy of periodic crystal and quasicrystals. Accurate free energies of numerous ordered phases, including dodecagonal, decagonal, and octagonal quasicrystals, are obtained for a simple model, i.e., the Lifshitz-Petrich free-energy functional, of soft quasicrystals with two length scales. The availability of the free energy allows us to construct phase diagrams of the system, demonstrating that, for the Lifshitz-Petrich model, the dodecagonal and decagonal quasicrystals can become stable phases, whereas the octagonal quasicrystal stays as a metastable phase.

The phase behaviour of a free energy functional with two length scales is examined by comparing the free energy of different candidate phases including three-dimensional icosahedral quasicrystals. Accurate free energy of the quasicrystals has been obtained using the recently developed projection method. The results reveal that the icosahedral quasicrystal and body-centred-cubic spherical phase are the stable ordered phases of the model. Furthermore, the difference between the results obtained from the projection method and the one-mode approximation has been analyzed in detail. The present study extends previous results on two-dimensional systems, demonstrating that the interactions between density waves at two length scales can stabilize two- and three-dimensional quasicrystals.

Quasiperiodic systems are important space-filling ordered structures, without decay and translational invariance. How to solve quasiperiodic systems accurately and efficiently is a great challenge. A useful approach, the projection method (PM) [ J. Comput. Phys., 256 (2014), pp. 428–440], has been proposed to compute quasiperiodic systems. Various studies have demonstrated that the PM is an accurate and efficient method to solve quasiperiodic systems. However, there is a lack of theoretical analysis of the PM. In this paper, we present a rigorous convergence analysis of the PM by establishing a mathematical framework of quasiperiodic functions and their high-dimensional periodic functions. We also give a theoretical analysis of the quasiperiodic spectral method (QSM) based on this framework. Results demonstrate that the PM and QSM both have exponential decay, and the QSM (PM) is a generalization of the periodic Fourier spectral (pseudospectral) method. Then, we analyze the computational complexity of the PM and QSM in calculating quasiperiodic systems. The PM can use a fast Fourier transform, while the QSM cannot. Moreover, we investigate the accuracy and efficiency of the PM, QSM, and periodic approximation method in solving the linear time-dependent quasiperiodic Schrödinger equation.

Large sparse linear algebraic systems can be found in a variety of scientific and engineering fields and many scientists strive to solve them in an efficient and robust manner. In this paper, we propose an interpretable neural solver, the Fourier neural solver (FNS), to address them. FNS is based on deep learning and a fast Fourier transform. Because the error between the iterative solution and the ground truth involves a wide range of frequency modes, the FNS combines a stationary iterative method and frequency space correction to eliminate different components of the error. Local Fourier analysis shows that the FNS can pick up on the error components in frequency space that are challenging to eliminate with stationary methods. Numerical experiments on the anisotropic diffusion equation, convection–diffusion equation, and Helmholtz equation show that the FNS is more efficient and more robust than the state-of-the-art neural solver.

We present an analysis of the approximation error for a [Formula: see text]-dimensional quasiperiodic function [Formula: see text] with Diophantine frequencies, approximated by a periodic function with the fundamental domain [Formula: see text]. When [Formula: see text] has a certain regularity, its global behavior can be described by a finite number of Fourier components and has a polynomial decay at infinity. The dominant part of the periodic approximation error is bounded by [Formula: see text], where [Formula: see text] belongs to the best simultaneous approximation sequence and [Formula: see text] is the number of different irrational elements in the [Formula: see text]th dimension component of Fourier frequencies, respectively. Meanwhile, we discuss the approximation rate. Finally, these analytical results are verified by some examples.