Related papers: Computing electrostatic potentials using regulariz…
We present a novel technique by which highly-segmented electrostatic configurations can be solved. The Robin Hood method is a matrix-inversion algorithm optimized for solving high density boundary element method (BEM) problems. We…
This paper introduces a novel numerical method for the inverse problem of electroencephalography(EEG). We pose the inverse EEG problem as an optimal control (OC) problem for Poisson's equation. The optimality conditions lead to a…
We study the computation of equilibrium points of electrostatic potentials: locations in space where the electrostatic force arising from a collection of charged particles vanishes. This is a novel scenario of optimization in which…
In this paper, we introduce and analyze a class of numerical schemes that demonstrate remarkable superiority in terms of efficiency, the preservation of positivity, energy stability, and high-order precision to solve the time-dependent…
This paper studies tensor-based Robust Principal Component Analysis (RPCA) using atomic-norm regularization. Given the superposition of a sparse and a low-rank tensor, we present conditions under which it is possible to exactly recover the…
Attaining a reliable complete basis set (CBS) limit remains a significant challenge in ab initio correlated electronic-structure calculations. Building on our previous work for atoms and diatomic molecules, we present a finite-element (FE)…
First-order energy dissipative schemes in time are available in literature for the Poisson-Nernst-Planck (PNP) equations, but second-order ones are still in lack. This work proposes novel second-order discretization in time and finite…
An accurate force calculation with the Poisson-Boltzmann equation is challenging, as it requires the electric field on the molecular surface. Here, we present a calculation of the electric field on the solute-solvent interface that is exact…
Over the past decade, Finite Element Method (FEM) has served as a foundational numerical framework for approximating the terms of Time Series Expansion (TSE) as solutions to transient Partial Differential Equation (PDE). However, the…
We propose a novel efficient algorithm to solve Poisson equation in irregular two dimensional domains for electrostatics. It can handle Dirichlet, Neumann or mixed boundary problems in which the filling media can be homogeneous or…
Constructing accurate, high dimensional molecular potential energy surfaces (PESs) for polyatomic molecules is challenging. Reproducing Kernel Hilbert space (RKHS) interpolation is an efficient way to construct such PESs. However, the…
In numerical simulations of many charged systems at the micro/nano scale, a common theme is the repeated solution of the Poisson-Boltzmann equation. This task proves challenging, if not entirely infeasible, largely due to the nonlinearity…
We consider the nonlinear Poisson-Boltzmann equation in the context of electrostatic models for a biological macromolecule, embedded in a bounded domain containing a solution of an arbitrary number of ionic species which is not necessarily…
This paper develops a probabilistic numerical method for solution of partial differential equations (PDEs) and studies application of that method to PDE-constrained inverse problems. This approach enables the solution of challenging inverse…
The three-dimensional electrostatic field configuration in gas ionization detectors has been simulated using an efficient and precise nearly exact boundary element method (NEBEM) solver set up to solve an integral equation of the first…
The recent development of spectral method has been praised for its high-order convergence in simulating complex physical problems. The combination of embedded boundary method and spectral method becomes a mainstream way to tackle…
We consider numerical methods for the Poisson-Nernst-Planck-Cahn-Hilliard (PNPCH) equations with steric interactions. We propose a novel energy stable numerical scheme that respects mass conservation and positivity at the discrete level.…
This work considers charged systems described by the modified Poisson--Nernst--Planck (PNP) equations, which incorporate ionic steric effects and the Born solvation energy for dielectric inhomogeneity. Solving the steady-state modified PNP…
Folding grid value vectors of size $2^L$ into $L$th order tensors of mode sizes $2\times \cdots\times 2$, combined with low-rank representation in the tensor train format, has been shown to lead to highly efficient approximations for…
Electrokinetic phenomena in nanopore sensors and microfluidic devices require accurate simulation of coupled fluid-electrostatic interactions in geometrically complex domains with irregular boundaries and adaptive mesh refinement. We…