Long-range Electrostatic Interactions Research Articles

Manipulation of Hierarchical Chiral Self-assembly and Anion Recognition by Supramolecular Systems of β-Glucopyranoside, Pillar[5]arenes, and Polyoxometalates.

Hierarchical chiral structures have broad applications in optical devices, asymmetric catalysis, and biological systems. The delicate balance of various interactions are key to the self-assembly of chiral structures. Herein, a ternary co-assembly consisting of cationic pillar[5]arenes (P5As), anionic β-glucopyranoside (βGlcD/βGlcL), and Anderson-type polyoxometalates (POMs) were constructed. Through adjusting the stoichiometry of βGlcD, the assemblies were effectively controlled to form hierarchical nano-leaf assemblies with twisted nanoribbons in a homochiral direction. The co-assemblies exhibit strong Cotton effects, and successfully induced the chirality of Anderson-type POMs. More interestingly, by changing the central metal in Anderson-type POMs (XMo6 O24 3- (X=Cr, Al, and Ga)), even though the three clusters have the same numbers of charge and size, the hierarchical chirality of the related assemblies varied in the morphology of the assemblies and the Cotton effect in the CD spectra. Results in theoretical calculations and ITC titration indicates that the tiny difference in long-range electrostatic interaction would result in the anion recognition of POMs, modulated by βGlcD through host-guest inclusion and hydrogen bonding in the assembly process.

Visualizing interfacial collective reaction behaviour of Li-S batteries.

Benefiting from high energy density (2,600 Wh kg-1) and low cost, lithium-sulfur (Li-S) batteries are considered promising candidates for advanced energy-storage systems1-4. Despite tremendous efforts in suppressing the long-standing shuttle effect of lithium polysulfides5-7, understanding of the interfacial reactions of lithium polysulfides at the nanoscale remains elusive. This is mainly because of the limitations of in situ characterization tools in tracing theliquid-solid conversion of unstable lithium polysulfides at high temporal-spatial resolution8-10. There is an urgent need to understand the coupled phenomena inside Li-S batteries, specifically, the dynamic distribution, aggregation, deposition and dissolution of lithium polysulfides. Here, by using in situ liquid-cell electrochemical transmission electron microscopy, we directly visualized the transformation of lithium polysulfides over electrode surfaces at the atomic scale. Notably, an unexpected gathering-induced collective charge transfer of lithium polysulfides was captured onthe nanocluster active-centre-immobilized surface. It further induced an instantaneous deposition of nonequilibrium Li2S nanocrystals from the dense liquid phase of lithium polysulfides. Without mediation of active centres, the reactions followed a classical single-molecule pathway, lithium polysulfides transforming into Li2S2 and Li2S step by step. Molecular dynamics simulations indicated that the long-range electrostatic interaction between active centres and lithium polysulfides promoted the formation of adense phase consisting of Li+ and Sn2- (2 < n ≤ 6), and the collective charge transfer in thedense phase was further verified by ab initio molecular dynamics simulations. The collective interfacial reaction pathway unveils a new transformation mechanism and deepens the fundamental understanding of Li-S batteries.

Charge-driven liquid-crystalline behavior of ligand-functionalized nanorods in apolar solvent.

Concentrated colloidal suspensions of nanorods often exhibit liquid-crystalline (LC) behavior. The transition to a nematic LC phase, with long-range orientational order of the particles, is usually well-captured by Onsager's theory for hard rods, at least qualitatively. The theory shows how the volume fraction at the transition decreases with increasing aspect ratio of the rods. It also explains that the long-range electrostatic repulsive interaction occurring between rods stabilized by their surface charge can significantly increase their effective diameter, resulting in a decrease in the volume fraction at the transition, as compared to sterically stabilized rods. Here, we report on a system of ligand-stabilized LaPO4 nanorods, of aspect ratio ≈11, dispersed in apolar medium exhibiting the counter-intuitive observation that the onset of nematic self-assembly occurs at an extremely low volume fraction of ≈0.25%, which is lower than observed (≈3%) with the same particles when charged-stabilized in polar solvent. Furthermore, the nanorod volume fraction at the transition increases with increasing concentration of ligands, in a similar way as in polar media where increasing the ionic strength leads to surface charge screening. This peculiar system was investigated by dynamic light scattering, Fourier-transform infrared spectroscopy, zetametry, electron microscopy, polarized light microscopy, photoluminescence measurements, and X-ray scattering. Based on these experimental data, we formulate several tentative scenarios that might explain this unexpected phase behavior. However, at this stage, its full understanding remains a pending theoretical challenge. Nevertheless, this study shows that dispersing anisotropic nanoparticles in an apolar solvent may sometimes lead to spontaneous ordering events that defy our intuitive ideas about colloidal systems.

PANNA 2.0: Efficient neural network interatomic potentials and new architectures.

We present the latest release of PANNA 2.0 (Properties from Artificial Neural Network Architectures), a code for the generation of neural network interatomic potentials based on local atomic descriptors and multilayer perceptrons. Built on a new back end, this new release of PANNA features improved tools for customizing and monitoring network training, better graphics processing unit support including a fast descriptor calculator, new plugins for external codes, and a new architecture for the inclusion of long-range electrostatic interactions through a variational charge equilibration scheme. We present an overview of the main features of the new code, and several benchmarks comparing the accuracy of PANNA models to the state of the art, on commonly used benchmarks as well as richer datasets.

Nucleoside Analog Reverse-Transcriptase Inhibitors in Membrane Environment: Molecular Dynamics Simulations.

The behavior of four drugs from the family of nucleoside analog reverse-transcriptase inhibitors (zalcitabine, stavudine, didanosine, and apricitabine) in a membrane environment was traced using molecular dynamics simulations. The simulation models included bilayers and monolayers composed of POPC and POPG phospholipids. It was demonstrated that the drugs have a higher affinity towards POPG membranes than POPC membranes due to attractive long-range electrostatic interactions. The results obtained for monolayers were consistent with those obtained for bilayers. The drugs accumulated in the phospholipid polar headgroup region. Two adsorption modes were distinguished. They differed in the degree of penetration of the hydrophilic headgroup region. Hydrogen bonds between drug molecules and phospholipid heads were responsible for adsorption. It was shown that apricitabine penetrated the hydrophilic part of the POPC and POPG membranes more effectively than the other drugs. Van der Waals interactions between S atoms and lipids were responsible for this.

EE-ONIOM-CT Method to Efficiently Account for the Missing Interactions in ONIOM: Energies and Analytic Gradients.

Hybrid methods such as ONIOM (QM:QM) are widely used for the study of local processes in large systems. However, the intrinsic need for system partitioning often leads to a less-than-desirable performance for some important chemical processes. This is due to the missing interactions in the chemically important model region (i.e., active site) of the high-level theory. The missing interactions can be categorized into two classes, viz. charge transfer (i.e., charge redistribution) and long-range electrostatic interactions. Our group presented two entirely different methods to treat these deficiencies individually. ONIOM-CT and ONIOM-EE methods have been demonstrated to improve the performance of ONIOM by incorporating charge transfer and missing electrostatic interactions, respectively. In general, the inclusion of the missing interactions separately in two different calculations may not be sufficient to reach a high accuracy. Thus, it is highly desirable to develop a method to correct both deficiencies simultaneously. In this work, we aim to connect the methods ONIOM-CT and ONIOM-EE for a more comprehensive treatment. A "stepwise" model was found to be necessary for a robust performance. This model employs a stepwise procedure by first satisfying the ONIOM-CT condition for charge balance before accounting for the electrostatic interactions from the rest of the system perturbatively. This has the advantage of easy interpretation due to the clear separation of the two effects. We demonstrate the performance of our method using embedding charges determined from a Mulliken population analysis. An efficient analytic gradient expression for this method is derived and implemented by requiring three sets of z-vector self-consistent equations. The performance of our method is assessed against full system calculations in high-level theory for a set of three proton transfer reactions representing different degrees of electrostatic embedding.

Q-pac: A Python package for machine learned charge equilibration models.

Many state-of-the art machine learning (ML) interatomic potentials are based on a local or semi-local (message-passing) representation of chemical environments. They, therefore, lack a description of long-range electrostatic interactions and non-local charge transfer. In this context, there has been much interest in developing ML-based charge equilibration models, which allow the rigorous calculation of long-range electrostatic interactions and the energetic response of molecules and materials to external fields. The recently reported kQEq method achieves this by predicting local atomic electronegativities using Kernel ML. This paper describes the q-pac Python package, which implements several algorithmic and methodological advances to kQEq and provides an extendable framework for the development of ML charge equilibration models.

Discovery of a Series of Indane-Containing NBTIs with Activity against Multidrug-Resistant Gram-Negative Pathogens.

The rise of multidrug-resistant (MDR) Gram-negative bacteria is a major global health problem necessitating the discovery of new classes of antibiotics. Novel bacterial topoisomerase inhibitors (NBTIs) target the clinically validated bacterial type II topoisomerases with a distinct binding site and mechanism of action to fluoroquinolone antibiotics, thus avoiding cross-resistance to this drug class. Here we report the discovery of a series of NBTIs incorporating a novel indane DNA binding moiety. X-ray cocrystal structures of compounds 2 and 17a bound to Staphylococcus aureus DNA gyrase-DNA were determined, revealing specific interactions with the enzyme binding pocket at the GyrA dimer interface and a long-range electrostatic interaction between the basic amine in the linker and the carboxylate of Asp83. Exploration of the structure-activity relationship within the series led to the identification of lead compound 18c, which showed potent broad-spectrum activity against a panel of MDR Gram-negative bacteria.

Hierarchical Model of Weak Polyelectrolytes with Ionization and Conformation Consistency

Conventional models of weak polyelectrolytes are formulated in terms of either an average degree of ionization or a rigid polymer conformation. While the decoupling of segment-level ionization and polymer structure simplifies the theoretical procedure and is beneficial from a computational perspective, it misses intrachain correlations (i.e., interactions beyond adjacent monomers) that are important for systems at a low polymer concentration with strong electrostatic interactions. In this work, we propose a hierarchical method that incorporates the liquid-state theories with the site-binding model to account for the effects of the local ionic environment on both inter- and intrachain correlations. A wormlike-chain model is introduced to describe the coupling of polymer conformation and monomeric ionization with the long-range electrostatic interactions accounted for by using the Gibbs–Bogoliubov variational principle. Through an extensive comparison of the theoretical predictions with experimental titration curves for poly(acrylic acid) in different alkali chloride solutions, we demonstrate that the hierarchical model is able to quantify the charging behavior of weak polyelectrolytes over a broad range of solution conditions.

Accurate Fourth-Generation Machine Learning Potentials by Electrostatic Embedding.

In recent years, significant progress has been made in the development of machine learning potentials (MLPs) for atomistic simulations with applications in many fields from chemistry to materials science. While most current MLPs are based on environment-dependent atomic energies, the limitations of this locality approximation can be overcome, e.g., in fourth-generation MLPs, which incorporate long-range electrostatic interactions based on an equilibrated global charge distribution. Apart from the considered interactions, the quality of MLPs crucially depends on the information available about the system, i.e., the descriptors. In this work we show that including─in addition to structural information─the electrostatic potential arising from the charge distribution in the atomic environments significantly improves the quality and transferability of the potentials. Moreover, the extended descriptor allows current limitations of two- and three-body based feature vectors to be overcome regarding artificially degenerate atomic environments. The capabilities of such an electrostatically embedded fourth-generation high-dimensional neural network potential (ee4G-HDNNP), which is further augmented by pairwise interactions, are demonstrated for NaCl as a benchmark system. Employing a data set containing only neutral and negatively charged NaCl clusters, even small energy differences between different cluster geometries can be resolved, and the potential shows an impressive transferability to positively charged clusters as well as the melt.

Dielectric Saturation in Water from a Long-Range Machine Learning Model.

Machine learning-based neural network potentials have the ability to provide ab initio-level predictions while reaching large length and time scales often limited to empirical force fields. Traditionally, neural network potentials rely on a local description of atomic environments to achieve this scalability. These local descriptions result in short-range models that neglect long-range interactions necessary for processes like dielectric screening in polar liquids. Several approaches to including long-range electrostatic interactions within neural network models have appeared recently, and here we investigate the transferability of one such model, the self-consistent field neural network (SCFNN), which focuses on learning the physics associated with long-range response. By learning the essential physics, one can expect that such a neural network model should exhibit at least partial transferability. We illustrate this transferability by modeling dielectric saturation in a SCFNN model of water. We show that the SCFNN model can predict nonlinear response at high electric fields, including saturation of the dielectric constant, without training the model on these high field strengths and the resulting liquid configurations. We then use these simulations to examine the nuclear and electronic structure changes underlying dielectric saturation. Our results suggest that neural network models can exhibit transferability beyond the linear response regime and make genuine predictions when the relevant physics is properly learned.

Thermodynamics of Light-Induced Nanoscale Polar Structures in Ferroelectric Superlattices.

We study the thermodynamics of nanoscale polar structures in PbTiO3/SrTiO3 ferroelectric superlattices induced by above-bandgap optical excitation using a phase-field model explicitly considering both structural and electronic processes. We demonstrate that the light-excited carriers provide the charge compensation of polarization bound charges and the lattice thermal energy, both of which are key to the thermodynamic stabilization of a previously observed supercrystal, a three-dimensionally periodic nanostructure, within a window of substrate strains, while different mechanical and electrical boundary conditions can stabilize a number of other nanoscale polar structures by balancing the competing short-range exchange interactions responsible for the domain wall energy and long-range electrostatic and elastic interactions. The insights into the light-induced formation and richness of nanoscale structures from this work offer theoretical guidance for exploring and manipulating the thermodynamic stability of nanoscale polar structures employing a combination of thermal, mechanical, and electrical stimuli as well as light.

Machine Learning Assisted Simulations of Electrochemical Interfaces: Recent Progress and Challenges.

The electrochemical interface, where the adsorption of reactants and electrocatalytic reactions take place, has long been a focus of attention. Some of the important processes on it tend to possess relatively slow kinetic characteristics, which are usually beyond the scope of ab initio molecular dynamics. The newly emerging technique, machine learning methods, provides an alternative approach to achieve thousands of atoms and nanosecond time scale while ensuring precision and efficiency. In this Perspective, we summarize in detail the recent progress and achievements made by the introduction of machine learning to simulate electrochemical interfaces, and focus on the limitations of current machine learning models, such as accurate descriptions of long-range electrostatic interactions and the kinetics of the electrochemical reactions occurring at the interface. Finally, we further point out the future directions for machine learning to expand in the field of electrochemical interfaces.

The Open Catalyst 2022 (OC22) Dataset and Challenges for Oxide Electrocatalysts

The development of machine learning models for electrocatalysts requires a broad set of training data to enable their use across a wide variety of materials. One class of materials that currently lacks sufficient training data is oxides, which are critical for the development of OER catalysts. To address this, we developed the OC22 dataset, consisting of 62,331 DFT relaxations (~9,854,504 single point calculations) across a range of oxide materials, coverages, and adsorbates. We define generalized total energy tasks that enable property prediction beyond adsorption energies; we test baseline performance of several graph neural networks; and we provide pre-defined dataset splits to establish clear benchmarks for future efforts. In the most general task, GemNet-OC sees a ~36% improvement in energy predictions when combining the chemically dissimilar OC20 and OC22 datasets via fine-tuning. Similarly, we achieved a ~19% improvement in total energy predictions on OC20 and a ~9% improvement in force predictions in OC22 when using joint training. We demonstrate the practical utility of a top performing model by capturing literature adsorption energies and important OER scaling relationships. We expect OC22 to provide an important benchmark for models seeking to incorporate intricate long-range electrostatic and magnetic interactions in oxide surfaces. Dataset and baseline models are open sourced, and a public leaderboard is available to encourage continued community developments on the total energy tasks and data.

SH3 binding pathway and kinetics in the presence of salt.

Protein-protein interactions are involved in a wide range of cellular processes, like cell signaling, in which these interactions often bind to intrinsically disordered proteins. SH3 domains are common protein interaction domains which are typically negatively charged and typically bind proline-rich disordered segments, which are positively charged. This project focuses on the role of electrostatics in the binding pathway of the yeast protein Abp1 SH3 domain (AbpSH3) and the intrinsically disordered peptide ArkA. We are investigating the binding pathway by simulating ArkA-AbpSH3 binding using molecular dynamics in the presence of sodium chloride, which serves to screen electrostatic interactions. Previous simulations in the absence of salt showed that ArkA initially forms a disordered encounter complex with AbpSH3 before reaching the fully bound state. In the absence of salt, the encounter complex, characterized by many transient non-specific electrostatic interactions, forms rapidly followed by a slower transition to the fully bound state. The addition of 800 mM sodium chloride leads to a slower encounter complex formation and a decrease in electrostatic contacts. Screening of long-range electrostatic interactions could also disrupt electrostatic steering to the binding surface during the initial encounter complex formation. Simulating the binding between ApbSH3 and ArkA in the presence of salt can provide insight into the role of electrostatics in SH3 binding in general and can be compared to experimental data on AbpSH3 binding at different salt concentrations.

Impact of Mineralizing Agents on Aluminum Distribution and Acidity of ZSM-5 Zeolites.

Intensive research on improving the catalytic properties of zeolites is focused on modulating their acidity and the distribution of associated Al sites. Herein, by studying a series of ZSM-5 zeolites over a broad range of Al content, we demonstrate how the nature of the mineralizing agent (F- or OH- ) used in hydrothermal syntheses directly impacts Al sites distribution. The proportions of Al sites, probed by 27 Al NMR, depend on the Si/Al ratio for F- , but remain identical for OH- (from Si/Al=30 to 760). This leads to contrasting variations in weak and strong acidities. Such opposite effect of mineralizers is explained by the spatial location of negative charges and the resulting balance between short- and long-range electrostatic interactions. This understanding paves the way for additional and simple opportunities to control zeolites' acidity.

DNA-protamine condensates under low salt conditions: molecular dynamics simulation with a simple coarse-grained model focusing on electrostatic interactions.

Protamine, a small, strongly positively-charged protein, plays a key role in achieving chromatin condensation inside sperm cells and is also involved in the formulation of nanoparticles for gene therapy and packaging of mRNA-based vaccines against viral infection and cancer. The detailed mechanisms of such condensations are still poorly understood especially under low salt conditions where electrostatic interaction predominates. Our previous study, with a refined coarse-grained model in full consideration of the long-range electrostatic interactions, has demonstrated the crucial role of electrostatic interaction in protamine-controlled reversible DNA condensation. Therefore, we herein pay our attention only to the electrostatic interaction and devise a coarser-grained bead-spring model representing the right linear charge density on protamine and DNA chains but treating other short-range interactions as simply as possible, which would be suitable for real-scale simulations. Effective pair potential calculations and large-scale molecular dynamics simulations using this extremely simple model reproduce the phase behaviour of DNA in a wide range of protamine concentrations under low salt conditions, again revealing the importance of the electrostatic interaction in this process and providing a detailed nanoscale picture of bundle formation mediated by a charge disproportionation mechanism. Our simulations also show that protamine length alters DNA overcharging and in turn redissolution thresholds of DNA condensates, revealing the important role played by entropies and correlated fluctuations of condensing agents and thus offering an additional opportunity to design tailored nanoparticles for gene therapy. The control mechanism of DNA-protamine condensates will also provide a better microscopic picture of biomolecular condensates, i.e., membraneless organelles arising from liquid-liquid phase separation, that are emerging as key principles of intracellular organization. Such condensates controlled by post-translational modification of protamine, in particular phosphorylation, or by variations in protamine length from species to species may also be responsible for the chromatin-nucleoplasm patterning observed during spermatogenesis in several vertebrate and invertebrate species.

A step-wise ion hydration model of aqueous electrolyte solution: The 2:2, 2:1 and 1:2 punches

In prior work we showed that the current step-wise hydration framework shows promise for modeling the concentration-activity relationship in aqueous solutions of ions representing 1:1 electrolytes, alcohols, and sugars, and accurately predicting the partitioning of the water between the bound (hydrated) pool and an unbound (free) pool. In the current work, we continue developing this framework by applying it to ions representing 2:2, 2:1 and 1:2 electrolyte in aqueous solutions. Properly treating ion hydration as an ideal solution along with employing mean spherical approximation of the long-range electrostatic interactions yields an accurate model for most electrolytes considered. This simple model of aqueous solutions is robust for electrolytes and organic solutes.

How Does Gecko Keratin Stick to Hydrophilic and Hydrophobic Surfaces in the Presence and Absence of Water? An Atomistic Molecular Dynamics Investigation.

We developed a united-atom model of gecko keratin to investigate the influence of electrostatic and van der Waals contributions to gecko adhesion in scenarios representing gecko's natural habitats. The keratin model assumes that only intrinsically disordered regions directly contact the surface. Contact angles of two generic substrate surfaces that we created match those previously used in AFM experiments on gecko adhesion. Force probe molecular dynamics simulations pulling the keratin from the surface show that the pull-off force increases with increased water content and is inversely related to the water contact angle of the surface. This matches experimental trends. We investigated the number and charge density of keratin and water at the surface, confirming a water-mediating effect and are able to show that keratin folds polar groups to the hydrophilic surface. We decomposed energetic contributions during pull-off, and our computational model shows that in contrast to popular hypotheses, long-range electrostatic interactions determine much of the pull-off process. The contribution of electrostatics to adhesion may be in the order of the van der Waals contributions.

Cluster expansions of multicomponent ionic materials: Formalism and methodology

The cluster expansion (CE) method has seen continuous and increasing use in the study of configuration-dependent properties of crystalline materials. The original development of the CE method along with the underlying mathematical formalism and assumptions was focused on the study of metallic alloys. Since then the methodology has been actively and successfully used in the study of ionic materials as well. In this work, we present a cohesive reformulation of the mathematical formalism underlying the CE method based on a synthesis of its original formulation, several additions and extensions that have been proposed since, and a revised representation of its constituent mathematical objects. We then proceed to describe some of the formal implications of using the methodology for charge-neutral configurations in ionic systems. In particular, we discuss the reduction of the size of configuration spaces and the resulting linear dependencies that arise among correlation functions that span the larger unconstrained configuration space. Additionally, we explore the effects of long-range electrostatic interactions. We also demonstrate how the previously proposed use of a point electrostatic term successfully accounts for the majority of the longer-range electrostatic interactions, and leaves the cluster expansion terms to capture mostly short-range interactions. Finally, we present and discuss a variety of recently developed methodologies, including training structure selection, oxidation state assignment, structure mapping, and regression algorithms, that are necessary to address these formal mathematical notions for a practical implementation of the CE method in the study of multicomponent ionic materials.