Search

We propose an unconditionally energy-stable, orthonormality-preserving, component-wise splitting iterative scheme for the Kohn-Sham gradient flow based model in the electronic structure calculation. We first study the scheme discretized in time but still continuous in space. The component-wise splitting iterative scheme changes one wave function at a time, similar to the Gauss-Seidel iteration for solving a linear equation system. At the time step n , the orthogonality of the wave function being updated to other wave functions is preserved by projecting the gradient of the Kohn-Sham energy onto the subspace orthogonal to all other wave functions known at the current time, while the normalization of this wave function is preserved by projecting the gradient of the Kohn-Sham energy onto the subspace orthogonal to this wave function at t n + 1 / 2. The unconditional energy stability is nontrivial, and it comes from a subtle treatment of the two-electron integral as well as a consistent treatment of the two projections. Rigorous mathematical derivations are presented to show our proposed scheme indeed satisfies the desired properties. We then study the fully-discretized scheme, where the space is further approximated by a conforming finite element subspace. For the fully-discretized scheme, not only the preservation of orthogonality and normalization (together we called orthonormalization) can be quickly shown using the same idea as for the semi-discretized scheme, but also the highlight property of the scheme, i.e., the unconditional energy stability can be rigorously proven. The scheme allows us to use large time step sizes and deal with small systems involving only a single wave function during each iteration step. Several numerical experiments are performed to verify the theoretical analysis, where the number of iterations is indeed greatly reduced as compared to similar examples solved by the Kohn-Sham gradient flow based model in the literature. • This paper proposes a novel and efficient numerical scheme for the Kohn-Sham gradient flow based model. • The scheme is an unconditionally energy-stable, orthonormality-preserving, component-wise splitting iterative scheme. • The scheme does not modify the original energy, allows large time step sizes, and solves small systems at each time step. • The rigorous proof is presented in the paper. • Several numerical examples are illustrated to verify the theoretical analysis. [ABSTRACT FROM AUTHOR]

• The real part of the star-genvalue problem corresponds to an eigenvalue problem in the Wigner formalism. • The imaginary part of the star-genvalue problem corresponds to the stationary Wigner equation. • The eigenvalue problem can be solved via the method of imaginary time propagation (ITP). • The stationary Wigner equation can be used to accelerate the convergence of ITP. • The Hermite spectral method turns out to be efficient in the numerical simulations. In this paper, the ground state Wigner function of a many-body system is explored theoretically and numerically. First, an eigenvalue problem for Wigner function is derived based on the energy operator of the system. The validity of finding the ground state through solving this eigenvalue problem is obtained by building a correspondence between its solution and the solution of stationary Schrödinger equation. Then, a numerical method is designed for solving proposed eigenvalue problem in one dimensional case, which can be briefly described by i) a simplified model is derived based on a quantum hydrodynamic model Cai et al. (2013) [7] to reduce the dimension of the problem, ii) an imaginary time propagation method is designed for solving the model, and numerical techniques such as solution reconstruction are proposed for the feasibility of the method. Results of several numerical experiments verify our method, in which the potential application of the method for large scale system is demonstrated by examples with density functional theory. [ABSTRACT FROM AUTHOR]

In this paper, we propose a parallel orbital-updating based optimization method for electronic structure calculations. With our method, the solution of the minimization problem for the Kohn-Sham energy functional with respect to N orbitals is replaced by the solution of N independent minimization problems for the energy functional with respect to one orbital and the orthogonalization of the N updated orbitals. This new method allows a two-level parallelization. This feature makes our approach has a great advantage in large scale parallel computing. The numerical experiments show that our new method is reliable and efficient. Hence, our new method has a great potential for large scale electronic structure calculations on modern supercomputers. • Propose a parallel orbital-updating based optimization approach for KS model. • The new approach has a feature of two-level parallelization. • Provide two practical parallel orbital-updating based optimization algorithms. • Numerical experiments show that these new algorithms are reliable and efficient. [ABSTRACT FROM AUTHOR]