Search papers, labs, and topics across Lattice.
The paper introduces Soft-FQEq, a differentiable layer for charge-aware machine learning potentials (MLIPs) that enforces fragment-resolved charge conservation in reactive systems. This method addresses the limitations of global charge equilibration (QEq) in capturing interfacial gradients at electrochemical interfaces, which leads to spurious charge transfer. By making fragment identification a differentiable function of atomic geometry, Soft-FQEq enables MLIPs to accurately model charge distributions and electrochemical potentials at interfaces, as demonstrated on IrO2/H2O/Na+/ClO4- interfaces.
Reactive MLIPs can now accurately model electrochemical interfaces thanks to a differentiable layer that dynamically identifies molecular fragments and conserves charge within them.
Predictive simulation of electrochemical interfaces requires atomistic models that capture reactive bond rearrangements, long-range electrostatics, and charge distributions reflecting the electronic distinctness of electrode and electrolyte. Existing charge-aware machine-learned interatomic potentials (MLIPs) built on global charge equilibration (QEq) settle electrode and electrolyte at a common electrochemical potential, leaving no room for the interfacial gradient that the double layer requires and admitting spurious charge transfer between electronically disconnected regions. Per-fragment charge equilibration is the established remedy in classical molecular dynamics, but reliance on predefined molecular topology has confined it to non-reactive systems. We lift this restriction by making fragment identification itself a differentiable function of atomic geometry, yielding soft fragment-constrained charge equilibration (Soft-FQEq) -- a solver layer that restores fragment-resolved charge conservation in reactive MLIPs. The layer consumes four scalar MLP readouts from a shared atomic-feature network -- per-atom electronegativity, source charge, short-range energy, and a soft bond connectivity -- and returns equilibrated charges together with per-fragment chemical potentials. We implement Soft-FQEq as an extension of the hippynn framework on a HIP-NN feature network and train it on DFT energies, forces, and DDEC6 charges for IrO2/H2O/Na+/ClO4- interfaces. The trained model recovers a clear electrode-to-electrolyte gradient in the per-atom electrochemical potential. With the same trained weights but the fragment-constrained solver replaced by global QEq at inference, this gradient collapses to an essentially uniform profile, directly showing that the gradient cannot be sustained within global QEq while the fragment formulation recovers it.
