Many-body attractions do not stabilize gas-liquid phase separation in aqueous dispersions of charged colloids within the Poisson-Boltzmann framework
Abstract
Attractive three-body interactions have been reported for like-charged colloids in low-salt suspensions, based on both finite-element Poisson-Boltzmann calculations and direct experimental measurements, and have been proposed as a mechanism to drive colloidal clustering. However, these Poisson-Boltzmann calculations typically neglect charge regulation and higher-order many-body effects. Here, we construct machine-learned (ML) many-body interaction potentials for charge-regulating colloids, trained on finite-element Poisson-Boltzmann calculations, to accurately capture three-body and higher-order contributions. We find that the three-body contribution to the many-body potential as obtained from Poisson-Boltzmann calculations on isolated colloid triplets is strongly attractive, consistent with previous work, whereas the four-body contribution for an equilateral pyramid configuration of four colloids is repulsive. We then construct ML many-body potentials for charged colloids using finite-element Poisson-Boltzmann calculations on clusters of 13 colloids, and find that the incorporation of higher-body interactions weakens the cohesive nature of the interactions. We identify a parameter regime exhibiting gas-liquid or gas-solid phase separation using the ML potentials in molecular dynamics simulations. However, when we include clusters of 48 colloids in the training data, the cohesion diminishes further, and molecular dynamics simulations using these potentials no longer include broad phase separation in aqueous dispersions of charged colloids. Finally, we compute the potential of mean force of pairs and triplets of colloids using primitive model simulations. We find that the resulting potentials are in good agreement with those obtained from the Poisson-Boltzmann calculations, thereby supporting the validity of the Poisson-Boltzmann approach for determining many-body interactions.
Versions
Related
- A Unified Framework for Bell-Shaped Gaussian Quadratures in Deterministic Boltzmann Solvers
- Near-real-time, meter-scale 3D urban wind modeling for low-altitude micrometeorology: numerical verification of a GPU-accelerated lattice Boltzmann framework
- ShARK: A Stochastic Transport Framework for the Relativistic Relaxation Time Approximation Boltzmann Equation
- A unified multirate lattice Boltzmann framework for thermosolutal dendritic solidification
- Is there an accelerating nonspreading wave packet in the Schrödinger equation with higher even-order dispersions?
BibTeX
@misc{rele2026manybody,
title = {{Many-body attractions do not stabilize gas-liquid phase separation in aqueous dispersions of charged colloids within the Poisson-Boltzmann framework}},
author = {Thijs ter Rele and René van Roij and Marjolein Dijkstra},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2606.21184},
eprint = {2606.21184},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2606.21184},
note = {CC BY-NC-ND 4.0},
}
RIS
TY - UNPB
TI - Many-body attractions do not stabilize gas-liquid phase separation in aqueous dispersions of charged colloids within the Poisson-Boltzmann framework
AU - Thijs ter Rele
AU - René van Roij
AU - Marjolein Dijkstra
PY - 2026
DA - 2026/09/18
PB - arXiv
DO - 10.48550/arxiv.2606.21184
UR - https://arxiv.org/abs/2606.21184
AB - Attractive three-body interactions have been reported for like-charged colloids in low-salt suspensions, based on both finite-element Poisson-Boltzmann calculations and direct experimental measurements, and have been proposed as a mechanism to drive colloidal clustering. However, these Poisson-Boltzmann calculations typically neglect charge regulation and higher-order many-body effects. Here, we construct machine-learned (ML) many-body interaction potentials for charge-regulating colloids, trained on finite-element Poisson-Boltzmann calculations, to accurately capture three-body and higher-order contributions. We find that the three-body contribution to the many-body potential as obtained from Poisson-Boltzmann calculations on isolated colloid triplets is strongly attractive, consistent with previous work, whereas the four-body contribution for an equilateral pyramid configuration of four colloids is repulsive. We then construct ML many-body potentials for charged colloids using finite-element Poisson-Boltzmann calculations on clusters of 13 colloids, and find that the incorporation of higher-body interactions weakens the cohesive nature of the interactions. We identify a parameter regime exhibiting gas-liquid or gas-solid phase separation using the ML potentials in molecular dynamics simulations. However, when we include clusters of 48 colloids in the training data, the cohesion diminishes further, and molecular dynamics simulations using these potentials no longer include broad phase separation in aqueous dispersions of charged colloids. Finally, we compute the potential of mean force of pairs and triplets of colloids using primitive model simulations. We find that the resulting potentials are in good agreement with those obtained from the Poisson-Boltzmann calculations, thereby supporting the validity of the Poisson-Boltzmann approach for determining many-body interactions.
ER -