BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Aihui Zhou

Aihui Zhou appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical Analysismath-phmath.MPmath.OCphysics.comp-phcond-mat.mtrl-scicond-mat.quant-gasquant-ph

Catalog footprint

What is connected

18works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2026arXiv

A quasi-orthogonal iterative method for eigenvalue problems

For large-scale eigenvalue problems requiring many mutually orthogonal eigenvectors, traditional numerical methods suffer substantial computational and communication costs with limited parallel scalability, primarily due to explicit orthogonalization. To address these challenges, we propose a quasi-orthogonal iterative method that dispenses with explicit orthogonalization and orthogonal initial data. It inherently preserves quasi-orthogonality (the iterates asymptotically tend to be orthogonal) and enhances robustness against numerical perturbations. Rigorous analysis confirms its energy-decay property and convergence of energy, gradient, and iterate. Numerical experiments validate the theoretical results, demonstrate key advantages of strong robustness and high-precision numerical orthogonality preservation, and thereby position our iterative method as an efficient, stable alternative for large-scale eigenvalue computations.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

An extended plane wave framework for the electronic structure calculations of twisted bilayer material systems

In this paper, we propose an extended plane wave framework to make the electronic structure calculations of the twisted bilayer 2D material systems practically feasible. Based on the foundation in [Y. Zhou, H. Chen, A. Zhou, J. Comput. Phys. 384, 99 (2019)], following extensions take place: (1) an tensor-producted basis set, which adopts PWs in the incommensurate dimensions, and localized basis in the interlayer dimension, (2) a practical application of a novel cutoff techniques we have recently developed, and (3) a quasi-band structure picture under the small twisted angles and weak interlayer coupling limits. With (1) and (2) now the dimensions of Hamiltonian matrix are reduced by about 2 orders of magnitude compared with the original framework. And (3) enables us to better organize the calculations and understand the results. For numerical examples, we study the electronic structures of the linear bilayer graphene lattice system with the magic twisted angle ($\sim 1.05^{\circ}$). The famous flat bands have been reproduced with their features in quantitative agreement with those from experiments and other theoretical calculations. Moreover, the extended framework has much less computational cost compared to the commensurate cell approximations, and is more extendable compared to the traditional model hamiltonians and tight binding models. Finally this framework can readily accommodate nonlinear models thus will laid the foundations for more effective yet accurate Density Functional Theory (DFT) calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Mathematical Analysis and Numerical Approximations of Density Functional Theory Models for Metallic Systems

In this paper, we investigate the energy minimization model of the ensemble Kohn-Sham density functional theory for metallic systems, in which a pseudo-eigenvalue matrix and a general smearing approach are involved. We study the invariance and the existence of the minimizer of the energy functional. We propose an adaptive double step size strategy and the corresponding preconditioned conjugate gradient methods for solving the energy minimization model. Under some mild but reasonable assumptions, we prove the global convergence of our algorithms. Numerical experiments show that our algorithms are efficient, especially for large scale metallic systems. In particular, our algorithms produce convergent numerical approximations for some metallic systems, for which the traditional self-consistent field iterations fail to converge.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Symmetrized two-scale finite element discretizations for partial differential equations with symmetric solutions

In this paper, a symmetrized two-scale finite element method is proposed for a class of partial differential equations with symmetric solutions. With this method, the finite element approximation on a fine tensor product grid is reduced to the finite element approximations on a much coarse grid and a univariant fine grid. It is shown by both theory and numerics including electronic structure calculations that the resulting approximation still maintains an asymptotically optimal accuracy. Consequently the symmetrized two-scale finite element method reduces computational cost significantly.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A plane wave study on the localized-extended transitions in the one-dimensional incommensurate systems

Based on our recently proposed plane wave framework, we theoretically study the localized-extended transition in the one dimensional incommensurate systems with cosine type of potentials, which are in close connection to many recent experiments in the ultracold atom and photonic crystal. We formulate a propagator based scattering picture for the transition at the ground state and single particle mobility edge, in which the deeper connection between the incommensurate potentials, eigenstate compositions and transition mechanism is revealed. We further show that there exists a upper limit of localization length for all localized eigenstates, leading to an fundamental difference to the Anderson localization. Numerical calculations are presented alongside the analysis to justify our statements. The theoretical analysis and numerical methods can also be generalized to systems in higher dimensions, with different potentials or beyond the single particle regime, which would benefit the future studies in the related fields.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Adaptive Step Size Strategy for Orthogonality Constrained Line Search Methods

In this paper, we propose an adaptive step size strategy for a class of line search methods for orthogonality constrained minimization problems, which avoids the classic backtracking procedure. We prove the convergence of the line search methods equipped with our adaptive step size strategy under some mild assumptions.We then apply the adaptive algorithm to electronic structure calculations, which show that our strategy is efficient and recommended.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Pactical Newton Methods for Electronic Structure Calculations

In this paper, we propose and analyze some practical Newton methods for electronic structure calculations. We show the convergence and the local quadratic convergence rate for the Newton method when the Newton search directions are well-obtained. In particular, we investigate some basic implementation issues in determining the search directions and step sizes which ensures the convergence of the subproblem at each iteration and accelerates the algorithm, respectively. It is shown by our numerical experiments that our Newton methods perform better than the existing conjugate gradient method, and the Newton method with the adaptive step size strategy is even more efficient.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Two-Grid based Adaptive Proper Orthogonal Decomposition Algorithm for Time Dependent Partial Differential Equations

In this article, we propose a two-grid based adaptive proper orthogonal decomposition (POD) method to solve the time dependent partial differential equations. Based on the error obtained in the coarse grid, we propose an error indicator for the numerical solution obtained in the fine grid. Our new algorithm is cheap and easy to be implement. We apply our new method to the solution of time-dependent advection-diffusion equations with the Kolmogorov flow and the ABC flow. The numerical results show that our method is more efficient than the existing POD methods.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Parallel Orbital-Updating Approach for Electronic Structure Calculations

In this paper, we propose an orbital iteration based parallel approach for electronic structure calculations. This approach is based on our understanding of the single-particle equations of independent particles that move in an effective potential. With this new approach, the solution of the single-particle equation is reduced to some solutions of independent linear algebraic systems and a small scale algebraic problem. It is demonstrated by our numerical experiments that this new approach is quite efficient for full-potential calculations for a class of molecular systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Accuracy of generalized gradient approximation functionals for density functional perturbation theory calculations

We assess the validity of various exchange-correlation functionals for computing the structural, vibrational, dielectric, and thermodynamical properties of materials in the framework of density-functional perturbation theory (DFPT). We consider five generalized-gradient approximation (GGA) functionals (PBE, PBEsol, WC, AM05, and HTBS) as well as the local density approximation (LDA) functional. We investigate a wide variety of materials including a semiconductor (silicon), a metal (copper), and various insulators (SiO$_2$ $α$-quartz and stishovite, ZrSiO$_4$ zircon, and MgO periclase). For the structural properties, we find that PBEsol and WC are the closest to the experiments and AM05 performs only slightly worse. All three functionals actually improve over LDA and PBE in contrast with HTBS, which is shown to fail dramatically for $α$-quartz. For the vibrational and thermodynamical properties, LDA performs surprisingly very good. In the majority of the test cases, it outperforms PBE significantly and also the WC, PBEsol and AM05 functionals though by a smaller margin (and to the detriment of structural parameters). On the other hand, HTBS performs also poorly for vibrational quantities. For the dielectric properties, none of the functionals can be put forward. They all (i) fail to reproduce the electronic dielectric constant due to the well-known band gap problem and (ii) tend to overestimate the oscillator strengths (and hence the static dielectric constant).

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Adaptive Finite Element Approximations for Kohn-Sham Models

The Kohn-Sham equation is a powerful, widely used approach for computation of ground state electronic energies and densities in chemistry, materials science, biology, and nanosciences. In this paper, we study the adaptive finite element approximations for the Kohn-Sham model. Based on the residual type a posteriori error estimators proposed in this paper, we introduce an adaptive finite element algorithm with a quite general marking strategy and prove the convergence of the adaptive finite element approximations. Using D{\" o}rfler's marking strategy, we then get the convergence rate and quasi-optimal complexity. We also carry out several typical numerical experiments that not only support our theory,but also show the robustness and efficiency of the adaptive finite element computations in electronic structure calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Convergence Rate and Quasi-Optimal Complexity of Adaptive Finite Element Computations for Multiple Eigenvalues

In this paper, we study an adaptive finite element method for multiple eigenvalue problems of a class of second order elliptic equations. By using some eigenspace approximation technology and its crucial property which is also presented in this paper, we extend the results in \cite{dai-xu-zhou08} to multiple eigenvalue problems, we obtain both convergence rate and quasi-optimal complexity of the adaptive finite element eigenvalue approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Gradient type optimization methods for electronic structure calculations

The density functional theory (DFT) in electronic structure calculations can be formulated as either a nonlinear eigenvalue or direct minimization problem. The most widely used approach for solving the former is the so-called self-consistent field (SCF) iteration. A common observation is that the convergence of SCF is not clear theoretically while approaches with convergence guarantee for solving the latter are often not competitive to SCF numerically. In this paper, we study gradient type methods for solving the direct minimization problem by constructing new iterations along the gradient on the Stiefel manifold. Global convergence (i.e., convergence to a stationary point from any initial solution) as well as local convergence rate follows from the standard theory for optimization on manifold directly. A major computational advantage is that the computation of linear eigenvalue problems is no longer needed. The main costs of our approaches arise from the assembling of the total energy functional and its gradient and the projection onto the manifold. These tasks are cheaper than eigenvalue computation and they are often more suitable for parallelization as long as the evaluation of the total energy functional and its gradient is efficient. Numerical results show that they can outperform SCF consistently on many practically large systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Symmetry-based Decomposition Approach to Eigenvalue Problems: Formulation, Discretization, and Implementation

In this paper, we propose a decomposition approach for eigenvalue problems with spatial symmetries, including the formulation, discretization as well as implementation. This approach can handle eigenvalue problems with either Abelian or non-Abelian symmetries, and is friendly for grid-based discretizations such as finite difference, finite element or finite volume methods. With the formulation, we divide the original eigenvalue problem into a set of subproblems and require only a smaller number of eigenpairs for each subproblem. We implement the decomposition approach with finite elements and parallelize our code in two levels. We show that the decomposition approach can improve the efficiency and scalability of iterative diagonalization. In particular, we apply the approach to solving Kohn--Sham equations of symmetric molecules consisting of hundreds of atoms.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Hohenberg-Kohn Theorem for Coulomb Type Systems

Density functional theory (DFT) has become a basic tool for the study of electronic structure of matter, in which the Hohenberg-Kohn theorem plays a fundamental role in the development of DFT. Unfortunately, the existing proofs are incomplete even incorrect; besides, the statement of the Hohenberg-Kohn theorem for many-electron Coulomb systems is not perfect. In this paper, we shall restate the Hohenberg-Kohn theorem for Coulomb type systems and present a rigorous proof by using the Fundamental Theorem of Algebra.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Numerical Analysis of Finite Dimensional Approximations of Kohn-Sham Models

In this paper, we study finite dimensional approximations of Kohn-Sham models, which are widely used in electronic structure calculations. We prove the convergence of the finite dimensional approximations and derive the a priori error estimates for ground state energies and solutions. We also provide numerical simulations for several molecular systems that support our theory.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Convergence and Optimal Complexity of Adaptive Finite Element Methods

In this paper, we study adaptive finite element approximations in a perturbation framework, which makes use of the existing adaptive finite element analysis of a linear symmetric elliptic problem. We prove the convergence and complexity of adaptive finite element methods for a class of elliptic partial differential equations. For illustration, we apply the general approach to obtain the convergence and complexity of adaptive finite element methods for a nonsymmetric problem, a nonlinear problem as well as an unbounded coefficient eigenvalue problem.

0 citations0 reviewsNo code linkedOpen access