Related Experiment Video
Updated: Jun 22, 2026

Rapid in-silico Battery Electrolyte Electrochemical Reaction Generation using 3T-VASP Multi-Scale Energy Minimization
Published on: August 22, 2025
Rigorous treatment of electrostatics for spatially varying dielectrics based on energy minimization
O I Obolensky1, T P Doerr, R Ray
1National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, Maryland 20894, USA.
Abstract:
An energy minimization formulation of electrostatics that allows computation of the electrostatic energy and forces to any desired accuracy in a system with arbitrary dielectric properties is presented. An integral equation for the scalar charge density is derived from an energy functional of the polarization vector field. This energy functional represents the true energy of the system even in nonequilibrium states. Arbitrary accuracy is achieved by solving the integral equation for the charge density via a series expansion in terms of the equation's kernel, which depends only on the geometry of the dielectrics. The streamlined formalism operates with volume charge distributions only, not resorting to introducing surface charges by hand. Therefore, it can be applied to any spatial variation of the dielectric susceptibility, which is of particular importance in applications to biomolecular systems. The simplicity of application of the formalism to real problems is shown with analytical and numerical examples.
Related Concept Videos
Electrostatic Boundary Conditions in Dielectrics
Consider a case where both the mediums across a boundary are two different dielectric materials. Recall that the electric field and electric displacement are proportional and related through the material's permittivity.
Electrostatic Boundary Conditions
The surface integral of an electric field is given by Gauss's law in integral form and is related to...
Gauss's Law in Dielectrics
Dielectric Polarization in a Capacitor
The Electrical Double Layer
Electric Field of a Charged Disk
The system's symmetry is in the cylindrical directions across the plane of the charge. As a result, the electric fields created by various surface charge elements nullify each other in the direction parallel to the surface. Thereby, the resulting electric field is perpendicular to the plane. Since the disk is...

