Perturbative many-body expansion for electrostatic energy and field for system of polarizable charged spherical ions in a dielectric medium

Recent Developments in the Methods and Applications of Electrostatic Theory

ConspectusThe review improves our understanding of how electrostatic interactions in the electrolyte, gas phase, and on surfaces can drive the fragmentation and assembly of particles. This is achieved through the overview of our advanced theoretical and computational modeling toolbox suitable for interpretation of experimental observations and discovery of novel, tunable assemblies and architectures. In the past decade, we have produced a significant, fundamental body of work on the development of comprehensive theories based on a rigorous mathematical foundation. These solutions are capable of accurate predictions of electrostatic interactions between dielectric particles of arbitrary size, anisotropy, composition, and charge, interacting in solvents, ionized medium, and on surfaces. We have applied the developed electrostatic approaches to describe physical and chemical phenomena in dusty plasma and planetary environments, in Coulomb fission and electrospray ionization processes, and in soft matter, including a counterintuitive but widespread attraction between like-charged particles.Despite its long history, the search for accurate methods to provide a deeper understanding of electrostatic interactions remains a subject of significant interest, as manifested by a constant stream of theoretical and experimental publications. While major international effort in this area has focused predominantly on the computational modeling of biocatalytic and biochemical performance, we have expanded the boundaries of accuracy, generality, and applicability of underlying theories. Simple solvation models, often used in calculating the electrostatic component of molecular solvation energy and polarization effects of solvent, rarely go beyond the induced dipole approximation because of computational costs. These approximations are generally adequate at larger separation distances; however, as particles approach the touching point, more advanced charged-induced multipolar descriptions of the electrostatic interactions are required to describe accurately a collective behavior of polarizable neutral and charged particles. At short separations, the electrostatic forces involving polarizable dielectric and conducting particles become nonadditive which necessitates further developments of quantitatively accurate many-body approaches. In applications, the electrostatic response of materials is commonly controlled by externally applied electric fields, an additional complex many-body problem that we have addressed most recently, both theoretically and numerically.This review reports on the most significant results and conclusions underpinning these recent advances in electrostatic theory and its applications. We first discuss the limitations of classical approaches to interpreting electrostatic phenomena in electrolytes and complex plasmas, leading to an extended analytical theory suitable for accurate estimation of the electrostatic forces in a dilute solution of a strong electrolyte. We then introduce the concept and numerical realization of many-body electrostatic theory focusing on its performance in selected experimental cases. These experiments underpin, among other applications, electrostatic self-assembly of two-dimensional lattice structures, melting of ionic colloidal crystals in an external electric field, and coalescence of charged clusters.

Manipulating Interactions between Dielectric Particles with Electric Fields: A General Electrostatic Many-Body Framework

We derive a rigorous analytical formalism and propose a numerical method for the quantitative evaluation of the electrostatic interactions between dielectric particles in an external electric field. This formalism also allows for inhomogeneous charge distributions, and, in particular, for the presence of pointlike charges on the particle surface. The theory is based on a boundary integral equation framework and yields analytical expressions for the interaction energy and net forces that can be computed in linear scaling cost, with respect to the number of interacting particles. We include numerical results that validate the proposed method and show the limitations of the fixed dipole approximation at small separation between interacting particles. The proposed method is also applied to study the stability and melting of ionic colloidal crystals in an external electric field.

From predictive modelling to machine learning and reverse engineering of colloidal self-assembly

An overwhelming diversity of colloidal building blocks with distinct sizes, materials and tunable interaction potentials are now available for colloidal self-assembly. The application space for materials composed of these building blocks is vast. To make progress in the rational design of new self-assembled materials, it is desirable to guide the experimental synthesis efforts by computational modelling. Here, we discuss computer simulation methods and strategies used for the design of soft materials created through bottom-up self-assembly of colloids and nanoparticles. We describe simulation techniques for investigating the self-assembly behaviour of colloidal suspensions, including crystal structure prediction methods, phase diagram calculations and enhanced sampling techniques, as well as their limitations. We also discuss the recent surge of interest in machine learning and reverse-engineering methods. Although their implementation in the colloidal realm is still in its infancy, we anticipate that these data-science tools offer new paradigms in understanding, predicting and (inverse) design of novel colloidal materials.

Dynamic simulations of many-body electrostatic self-assembly

Two experimental studies relating to electrostatic self-assembly have been the subject of dynamic computer simulations, where the consequences of changing the charge and the dielectric constant of the materials concerned have been explored. One series of calculations relates to experiments on the assembly of polymer particles that have been subjected to tribocharging and the simulations successfully reproduce many of the observed patterns of behaviour. A second study explores events observed following collisions between single particles and small clusters composed of charged particles derived from a metal oxide composite. As before, observations recorded during the course of the experiments are reproduced by the calculations. One study in particular reveals how particle polarizability can influence the assembly process.This article is part of the theme issue 'Modern theoretical chemistry'.

Free Energy Landscapes of Alanine Oligopeptides in Rigid-Body and Hybrid Water Models

Replica exchange molecular dynamics is used to study the effect of different rigid-body (mTIP3P, TIP4P, SPC/E) and hybrid (H1.56, H3.00) water models on the conformational free energy landscape of the alanine oligopeptides (acAnme and acA5nme), in conjunction with the CHARMM22 force field. The free energy landscape is mapped out as a function of the Ramachandran angles. In addition, various secondary structure metrics, solvation shell properties, and the number of peptide-solvent hydrogen bonds are monitored. Alanine dipeptide is found to have similar free energy landscapes in different solvent models, an insensitivity which may be due to the absence of possibilities for forming i-(i + 4) or i-(i + 3) intrapeptide hydrogen bonds. The pentapeptide, acA5nme, where there are three intrapeptide backbone hydrogen bonds, shows a conformational free energy landscape with a much greater degree of sensitivity to the choice of solvent model, though the three rigid-body water models differ only quantitatively. The pentapeptide prefers nonhelical, non-native PPII and β-sheet populations as the solvent is changed from SPC/E to the less tetrahedral liquid (H1.56) to an LJ-like liquid (H3.00). The pentapeptide conformational order metrics indicate a preference for open, solvent-exposed, non-native structures in hybrid solvent models at all temperatures of study. The possible correlations between the properties of solvent models and secondary structure preferences of alanine oligopeptides are discussed, and the competition between intrapeptide, peptide-solvent, and solvent-solvent hydrogen bonding is shown to be crucial in the relative free energies of different conformers.