Macroscopic surface charges from microscopic simulations

A theory for the stabilization of polar crystal surfaces by a liquid environment

Polar crystal surfaces play an important role in the functionality of many materials and have been studied extensively over many decades. In this article, a theoretical framework is presented that extends existing theories by placing the surrounding solution environment on an equal footing with the crystal itself; this is advantageous, e.g., when considering processes such as crystal growth from solution. By considering the polar crystal as a stack of parallel plate capacitors immersed in a solution environment, the equilibrium adsorbed surface charge density is derived by minimizing the free energy of the system. In analogy to the well-known diverging surface energy of a polar crystal surface at zero temperature, for a crystal in solution it is shown that the "polar catastrophe" manifests as a diverging free energy cost to perturb the system from equilibrium. Going further than existing theories, the present formulation predicts that fluctuations in the adsorbed surface charge density become increasingly suppressed with increasing crystal thickness. We also show how, in the slab geometry often employed in both theoretical and computational studies of interfaces, an electric displacement field emerges as an electrostatic boundary condition, the origins of which are rooted in the slab geometry itself, rather than the use of periodic boundary conditions. This aspect of the work provides a firmer theoretical basis for the recent observation that standard "slab corrections" fail to correctly describe, even qualitatively, polar crystal surfaces in solution.

Chemical Nature of Heterogeneous Electrofreezing of Supercooled Water Revealed on Polar (Pyroelectric) Surfaces

The ability to control the icing temperature of supercooled water (SCW) is of supreme importance in subfields of pure and applied sciences. The ice freezing of SCW can be influenced heterogeneously by electric effects, a process known as electrofreezing. This effect was first discovered during the 19th century; however, its mechanism is still under debate. In this Account we demonstrate, by capitalizing on the properties of polar crystals, that heterogeneous electrofreezing of SCW is a chemical process influenced by an electric field and specific ions. Polar crystals possess a net dipole moment. In addition, they are pyroelectric, displaying short-lived surface charges at their hemihedral faces at the two poles of the crystals as a result of temperature changes. Accordingly, during cooling or heating, an electric field is created, which is negated by the attraction of compensating charges from the environment. This process had an impact in the following experiments. The icing temperatures of SCW within crevices of polar crystals are higher in comparison to icing temperatures within crevices of nonpolar analogs. The role played by the electric effect was extricated from other effects by the performance of icing experiments on the surfaces of pyroelectric quasi-amorphous SrTiO₃. During those studies it was found that on positively charged surfaces the icing temperature of SCW is elevated, whereas on negatively charged surfaces it is reduced. Following investigations discovered that the icing temperature of SCW is impacted by an ionic current created within a hydrated layer on top of hydrophilic faces residing parallel to the polar axes of the crystals. In the absence of such current on analogous hydrophobic surfaces, the pyroelectric effect does not influence the icing temperature of SCW. Those results implied that electrofreezing of SCW is a process influenced by specific compensating ions attracted by the pyroelectric field from the aqueous solution. When freezing experiments are performed in an open atmosphere, bicarbonate and hydronium ions, created by the dissolution of atmospheric CO₂ in water, influence the icing temperature. The bicarbonate ions, when attracted by positively charged pyroelectric surfaces, elevate the icing temperature, whereas their counterparts, hydronium ions, when attracted by the negatively charged surfaces reduce the icing temperature. Molecular dynamic simulations suggested that bicarbonate ions, concentrated within the near positively charged interfacial layer, self-assemble with water molecules to create stabilized slightly distorted “ice-like” hexagonal assemblies which mimic the hexagons of the crystals of ice. This occurs by replacing, within those ice-like hexagons, two hydrogen bonds of water by C–O bonds of the HCO₃^– ion. On the basis of these simulations, it was predicted and experimentally confirmed that other trigonal planar ions such as NO₃^–, guanidinium⁺, and the quasi-hexagonal biguanidinium⁺ ion elevate the icing temperature. These ions were coined as “ice makers”. Other ions including hydronium, Cl^–, and SO₄^–2 interfere with the formation of ice-like assemblies and operate as “ice breakers”. The higher icing temperatures induced within the crevices of the hydrophobic polar crystals in comparison to the nonpolar analogs can be attributed to the proton ordering of the water molecules. In contrast, the icing temperatures on related hydrophilic surfaces are influenced both by compensating charges and by proton ordering.

Long-Range Ionic and Short-Range Hydration Effects Govern Strongly Anisotropic Clay Nanoparticle Interactions

The aggregation of clay particles in aqueous solution is a ubiquitous everyday process of broad environmental and technological importance. However, it is poorly understood at the all-important atomistic level since it depends on a complex and dynamic interplay of solvent-mediated electrostatic, hydrogen bonding, and dispersion interactions. With this in mind, we have performed an extensive set of classical molecular dynamics simulations (included enhanced sampling simulations) on the interactions between model kaolinite nanoparticles in pure and salty water. Our simulations reveal highly anisotropic behavior, in which the interaction between the nanoparticles varies from attractive to repulsive depending on the relative orientation of the nanoparticles. Detailed analysis reveals that at large separation (>1.5 nm), this interaction is dominated by electrostatic effects, whereas at smaller separations, the nature of the water hydration structure becomes critical. This study highlights an incredible richness in how clay nanoparticles interact, which should be accounted for in, for example, coarse-grained models of clay nanoparticle aggregation.

Self-consistent determination of long-range electrostatics in neural network potentials

Machine learning has the potential to revolutionize the field of molecular simulation through the development of efficient and accurate models of interatomic interactions. Neural networks can model interactions with the accuracy of quantum mechanics-based calculations, but with a fraction of the cost, enabling simulations of large systems over long timescales. However, implicit in the construction of neural network potentials is an assumption of locality, wherein atomic arrangements on the nanometer-scale are used to learn interatomic interactions. Because of this assumption, the resulting neural network models cannot describe long-range interactions that play critical roles in dielectric screening and chemical reactivity. Here, we address this issue by introducing the self-consistent field neural network - a general approach for learning the long-range response of molecular systems in neural network potentials that relies on a physically meaningful separation of the interatomic interactions - and demonstrate its utility by modeling liquid water with and without applied fields.