Search

Water mediates electrostatic interactions via the orientation of its dipoles around ions, molecules, and interfaces. This induced water polarization consequently influences multiple phenomena. In particular, water polarization modulated by nanoconfinement affects ion adsorption and transport, biomolecular self-assembly, and surface chemical reactions. Therefore, it is of paramount importance to understand how water-mediated interactions change at the nanoscale. Here we show that near the graphene surface anion-cation interactions do not obey the translational and isotropic symmetries of Coulomb's law. We identify a new property, referred to as non-reciprocity, which describes the non-equivalent and directional interaction between two oppositely charged ions near the confining surface when their positions with respect to the interface are exchanged. Specifically, upon exchange of the two ions' positions along the surface normal direction the interaction energy changes by about 5$k_BT$. In both cases, confinement enhances the attraction between two oppositely charged ions near the graphene surface, while intercalation of one ion into the graphene layers shifts the interaction to repulsive. While the water permittivity in confinement is different from that in bulk, the effects observed here via molecular dynamics simulations and X-ray reflectivity experiments cannot be accounted for by current permittivity models. Our work shows that the water structure is not enough to infer electrostatic interactions near interfaces., Comment: 18 pages, 13 figures, Accepted in Physical Review Research

Solid polymer electrolytes are considered a promising alternative to traditional liquid electrolytes in energy storage applications because of their good mechanical properties, and excellent thermal and chemical stability. A gap, however, still exists in understanding ion transport mechanisms and improving ion transport in solid polymer electrolytes. Therefore, it is crucial to bridge composition--structure and structure--property relationships. Here we demonstrate that size asymmetry, $\lambda$, represented by the ratio of counterion to charged monomer size, plays a key role in both the nanostructure and in the ionic dynamics. More specifically, when the nanostructure is modified by the external electric field such that the mobility cannot be described by linear response theory, two situations arise. The ionic mobility increases as $\lambda$ decreases (small counterions) in the weak electrostatics (high dielectric constant) regime. Whereas in systems with strong electrostatic interactions, ionomers with higher size symmetry ($\lambda \approx 1$) display higher ionic mobility. Moreover, ion transport is found to be dominated by the hopping of the ions and not by moving ionic clusters (also known as "vehicular'' charge transport). These results serve as a guide for designing ion-containing polymers for ion transport related applications.

We demonstrate that a highly charged polyelectrolyte confined in a spherical cavity undergoes reversible transformations between amorphous conformations to a four-fold symmetry morphology as a function of dielectric mismatch between the media inside and outside the cavity. Surface polarization due to dielectric mismatch exhibits an extra "confinement" effect, which are most pronounced within a certain range of the cavity radius and the electrostatic strength between the monomers and counterions and multivalent counterions. For cavities with a charged surface, surface polarization leads to an increased amount of counterions adsorbed in the outer side, further compressing the confined polyelectrolyte into a four-fold symmetry morphology. The equilibrium conformation of the chain is dependent upon several key factors including the relative permittivities of the media inside and outside the cavity, multivalent counterion concentration, cavity radius relative to the chain length, and interface charge density. Our findings offer insights into the effects of dielectric mismatch in packaging and delivery of polyelectrolytes across media with different relative permittivities. Moreover, the reversible transformation of the polyelectrolyte conformations in response to environmental permittivity allows for potential applications in biosensing and medical monitoring.

Polyelectrolytes under confinement are crucial for energy storage and for understanding biomolecular functions. Using molecular dynamics simulations, we analyze a polyelectrolyte solution confined between two oppositely charged planar dielectric surfaces and include surface polarization effects due to dielectric mismatch at the two electrodes. Although the effect of polarization on the charge distribution seems minor, we find that surface polarization enhances energy storage and also leads to the emergence of negative differential capacitance in confined polyelectrolyte solutions.

We present a novel Ensemble Monte Carlo Growth method to sample the equilibrium thermodynamic properties of random chains. The method is based on the multicanonical technique of computing the density of states in the energy space. Such a quantity is temperature independent, and therefore microcanonical and canonical thermodynamic quantities, including the free energy, entropy, and thermal averages, can be obtained by re-weighting with a Boltzmann factor. The algorithm we present combines two approaches: the first is the Monte Carlo ensemble growth method, where a "population" of samples in the state space is considered, as opposed to traditional sampling by long random walks, or iterative single-chain growth. The second is the flat-histogram Monte Carlo, similar to the popular Wang-Landau sampling, or to multicanonical chain-growth sampling. We discuss the performance and relative simplicity of the proposed algorithm, and we apply it to known test cases., Comment: 7 pages, 6 figures; Added references and figures, corrected typos, improved notation

We describe a highly optimized implementation of MPI domain decomposition in a GPU-enabled, general-purpose molecular dynamics code, HOOMD-blue (Anderson and Glotzer, arXiv:1308.5587). Our approach is inspired by a traditional CPU-based code, LAMMPS (Plimpton, J. Comp. Phys. 117, 1995), but is implemented within a code that was designed for execution on GPUs from the start (Anderson et al., J. Comp. Phys. 227, 2008). The software supports short-ranged pair force and bond force fields and achieves optimal GPU performance using an autotuning algorithm. We are able to demonstrate equivalent or superior scaling on up to 3,375 GPUs in Lennard-Jones and dissipative particle dynamics (DPD) simulations of up to 108 million particles. GPUDirect RDMA capabilities in recent GPU generations provide better performance in full double precision calculations. For a representative polymer physics application, HOOMD-blue 1.0 provides an effective GPU vs. CPU node speed-up of 12.5x., Comment: 30 pages, 14 figures