Integrating Charge Equilibration with Equivariant Machine-Learning Interatomic Potentials
Martin Vondrák1,2, William J Baldwin3, Gábor Csányi3,4
1University of Bayreuth, Bavarian Center for Battery Technology (BayBatt), Bayreuth 95447, Germany.
Journal of Chemical Theory and Computation
|June 20, 2026
Summary
Machine-learning interatomic potentials enhanced with charge equilibration accurately model charged defects in ZnO. However, limitations arise in larger systems, highlighting the need for advanced architectures in materials science.
Area of Science:
- Materials Science
- Computational Chemistry
- Physics
Background:
- Machine-learning interatomic potentials (MLIPs) excel in modeling atomic interactions but struggle with long-range electrostatic effects and charge transfer.
- Existing MLIPs often fail to accurately capture phenomena crucial for materials properties, such as nonlocal electronic behavior.
Purpose of the Study:
- To augment the Multi-Atomic Cluster Expansion (MACE) potential with a charge equilibration (QEq) framework for self-consistent charge redistribution.
- To evaluate the capabilities and limitations of ML-enhanced QEq in modeling charged defects and transferable potentials.
Main Methods:
- Integration of a charge equilibration (QEq) framework into the equivariant Multi-Atomic Cluster Expansion (MACE) potential.
- Application of the ML-enhanced QEq model to charged oxygen vacancies in wurtzite ZnO and a transferable water potential.
Main Results:
- The ML-enhanced QEq model accurately reproduced charge state-dependent relaxations and migration pathways for ZnO defects in small to medium supercells.
- Limitations of the quadratic QEq formalism were observed in larger ZnO systems, leading to charge delocalization and loss of distinct charge states.
- Initializing long-range water models from pretrained short-range representations improved data efficiency and transferability from gas-phase clusters to bulk liquid.
Conclusions:
- ML-enhanced QEq shows promise for capturing complex defect physics and enabling transferable potentials.
- The study highlights the fundamental constraints of QEq formalisms in large systems and the importance of physically informed architectures.
- Pretrained representations are crucial for extending ML potentials to systems dominated by long-range electrostatics and nonlocal charge response.
Related Concept Videos
Equilibrium Conditions for a Particle
When an object is in equilibrium, it is either at rest or moving with a constant velocity. There are two types of equilibrium: static and dynamic. Static equilibrium occurs when an object is at rest, while dynamic equilibrium occurs when an object is moving with a constant velocity. In both cases, there must be a balance of forces acting on the object.
To understand the concept of equilibrium, let us first consider the forces acting on an object. When different forces act on an object, they can...
To understand the concept of equilibrium, let us first consider the forces acting on an object. When different forces act on an object, they can...
The Electrical Double Layer
In the region where two bulk phases meet, an intricate electric charge distribution arises due to charge transfer, ion adsorption, molecular orientation, and charge distortion. This complex distribution is commonly referred to as the electrical double layer.When a solid electrode interfaces with ions in an electrolyte solution, the speed of electron transfer dictates the rates of oxidation and reduction. The electrode acquires a charge through the escape of atoms into the solution as cations or...
Energy Associated With a Charge Distribution
The work done to bring a charge through a distance r is given by the potential difference between the initial and the final position. To assemble a collection of point charges, the total work done can be expressed in terms of the product of each pair of charges divided by their separation distance, defined with respect to a suitable origin. Solving this expression gives the energy stored in a point charge distribution.
Potential-Energy Criterion for Equilibrium
Potential energy or potential function plays an essential role in determining the stability of a mechanical system. If a system is subjected to both gravitational and elastic forces, the potential function of the system can be expressed as the algebraic sum of gravitational and elastic potential energy. If the system is in equilibrium and is displaced by a small amount, then the work done on the system equals the negative of the change in the system's potential energy from the initial to the...
Calculations of Electric Potential I
Consider a ring of radius R with a uniform charge density λ. What will the electric potential be at point M, which is located on the axis of the ring at a distance x from the center of the ring?
The ring is divided into infinitesimal small arcs such that point M is equidistant from all the arcs. Here, the cylindrical coordinate system is used to calculate the electric potential at point M. A general element of the arc between angles θ and θ + dθ is of the length Rdθ and has a charge of λRdθ.
The ring is divided into infinitesimal small arcs such that point M is equidistant from all the arcs. Here, the cylindrical coordinate system is used to calculate the electric potential at point M. A general element of the arc between angles θ and θ + dθ is of the length Rdθ and has a charge of λRdθ.
Calculations of Electric Potential II
An electric dipole is a system of two equal but opposite charges, separated by a fixed distance. This system is used to model many real-world systems, including atomic and molecular interactions. One of these systems is the water molecule, but only under certain circumstances. These circumstances are met inside a microwave oven, where electric fields with alternating directions make the water molecules change orientation. This vibration is equivalent to heat at the molecular level.
Consider a...
Consider a...

