Related papers: A convergence analysis of a structure-preserving g…
We present for static density functional theory and time-dependent density functional theory calculations an all-electron method which employs high-order hierarchical finite element bases. Our mesh generation scheme, in which structured…
Elastic flow for closed curves can involve significant deformations. Mesh-based approximation schemes require tangentially redistributing vertices for long-time computations. We present and analyze a method that uses the Dirichlet energy…
In the present work, we introduce a Self-Consistent Density-Functional Embedding technique, which leaves the realm of standard energy-functional approaches in Density Functional Theory and targets directly the density-to-potential mapping…
In this paper, we present an efficient numerical method to address a thermodynamically consistent gas flow model in porous media involving compressible gas and deformable rock. The accurate modeling of gas flow in porous media often poses…
The analysis of structure-preserving numerical methods for the Poisson--Nernst--Planck (PNP) system has attracted growing interests in recent years. In this work, we provide an optimal rate convergence analysis and error estimate for finite…
Accurately modeling the electronic structure of water across scales, from individual molecules to bulk liquid, remains a grand challenge. Traditional computational methods face a critical trade-off between computational cost and efficiency.…
The ground state electron density -- obtainable using Kohn-Sham Density Functional Theory (KS-DFT) simulations -- contains a wealth of material information, making its prediction via machine learning (ML) models attractive. However, the…
We develop structure-preserving numerical methods for the Serre-Green-Naghdi equations, a model for weakly dispersive free-surface waves. We consider both the classical form, requiring the inversion of a non-linear elliptic operator, and a…
Kohn-Sham spin-density functional theory provides an efficient and accurate model to study electron-electron interaction effects in quantum dots, but its application to large systems is a challenge. An efficient algorithm for the…
Kohn-Sham density functional theory (KS-DFT) has found widespread application in accurate electronic structure calculations. However, it can be computationally demanding especially for large-scale simulations, motivating recent efforts…
We present a multigrid algorithm for self consistent solution of the Kohn-Sham equations in real space. The entire problem is discretized on a real space mesh with a high order finite difference representation. The resulting self consistent…
A model is developed, based on the density functional perturbation theory and the inverse Kohn-Sham method, that can be used to improve relativistic nuclear energy density functionals towards an exact but unknown Kohn-Sham…
We propose a novel Skew Gradient Embedding (SGE) framework for systematically reformulating thermodynamically consistent partial differential equation (PDE) models-capturing both reversible and irreversible processes-as generalized gradient…
As a counterpoint to recent numerical methods for crystal surface evolution, which agree well with microscopic dynamics but suffer from significant stiffness that prevents simulation on fine spatial grids, we develop a new numerical method…
In this paper, we consider numerical approximation of constrained gradient flows of planar closed curves, including the Willmore and the Helfrich flows. These equations have energy dissipation and the latter has conservation properties due…
In this paper, a novel adaptive finite element method is proposed to solve the Kohn-Sham equation based on the moving mesh (nonnested mesh) adaptive technique and the augmented subspace method. Different from the classical self-consistent…
We present a computationally efficient approach to perform large-scale all-electron density functional theory calculations by enriching the classical finite element basis with compactly supported atom-centered numerical basis functions that…
A finite difference numerical scheme is proposed and analyzed for the Cahn-Hilliard-Stokes system with Flory-Huggins energy functional. A convex splitting is applied to the chemical potential, which in turns leads to the implicit treatment…
In this paper, a novel $h$-adaptive isogeometric solver utilizing high-order hierarchical splines is proposed to solve the all-electron Kohn--Sham equation. In virtue of the smooth nature of Kohn--Sham wavefunctions across the domain,…
Kohn-Sham (KS) density functional theory (DFT) is a very efficient method for calculating various properties of solids as, for instance, the total energy, the electron density, or the electronic band structure. The KS-DFT method leads to…