Dynamically Polarizable Force Fields for Surface Simulations via Multi-output Classification Neural Networks

Simulation of polymeric mixed ionic and electronic conductors with a combined classical and quantum mechanical model

In organic polymeric materials with mixed ionic and electronic conduction (OMIEC), the excess charge in doped polymers is very mobile and the dynamics of the polymer chain cannot be accurately described with a model including only fixed point charges. Ions and polymer are comparatively slower and a methodology to capture the correlated motions of excess charge and ions is currently unavailable. Considering a prototypical interface encountered in this type of materials, we constructed a scheme based on the combination of MD and QM/MM to evaluate the classical dynamics of polymer, water and ions, while allowing the excess charge of the polymer chains to rearrange following the external electrostatic potential. We find that the location of the excess charge varies substantially between chains. The excess charge changes across multiple timescales as a result of fast structural fluctuations and slow rearrangement of the polymeric chains. Our results indicate that such effects are likely important to describe the phenomenology of OMIEC, but additional features should be added to the model to enable the study of processes such as electrochemical doping.

Constant chemical potential-quantum mechanical-molecular dynamics simulations of the graphene-electrolyte double layer

We present the coupling of two frameworks-the pseudo-open boundary simulation method known as constant potential molecular dynamics simulations (CμMD), combined with quantum mechanics/molecular dynamics (QMMD) calculations-to describe the properties of graphene electrodes in contact with electrolytes. The resulting CμQMMD model was then applied to three ionic solutions (LiCl, NaCl, and KCl in water) at bulk solution concentrations ranging from 0.5 M to 6 M in contact with a charged graphene electrode. The new approach we are describing here provides a simulation protocol to control the concentration of electrolyte solutions while including the effects of a fully polarizable electrode surface. Thanks to this coupling, we are able to accurately model both the electrode and solution side of the double layer and provide a thorough analysis of the properties of electrolytes at charged interfaces, such as the screening ability of the electrolyte and the electrostatic potential profile. We also report the calculation of the integral electrochemical double layer capacitance in the whole range of concentrations analyzed for each ionic species, while the quantum mechanical simulations provide access to the differential and integral quantum capacitance. We highlight how subtle features, such as the adsorption of potassium graphene or the tendency of the ions to form clusters contribute to the ability of graphene to store charge, and suggest implications for desalination.

Fluids and Electrolytes under Confinement in Single-Digit Nanopores

Confined fluids and electrolyte solutions in nanopores exhibit rich and surprising physics and chemistry that impact the mass transport and energy efficiency in many important natural systems and industrial applications. Existing theories often fail to predict the exotic effects observed in the narrowest of such pores, called single-digit nanopores (SDNs), which have diameters or conduit widths of less than 10 nm, and have only recently become accessible for experimental measurements. What SDNs reveal has been surprising, including a rapidly increasing number of examples such as extraordinarily fast water transport, distorted fluid-phase boundaries, strong ion-correlation and quantum effects, and dielectric anomalies that are not observed in larger pores. Exploiting these effects presents myriad opportunities in both basic and applied research that stand to impact a host of new technologies at the water-energy nexus, from new membranes for precise separations and water purification to new gas permeable materials for water electrolyzers and energy-storage devices. SDNs also present unique opportunities to achieve ultrasensitive and selective chemical sensing at the single-ion and single-molecule limit. In this review article, we summarize the progress on nanofluidics of SDNs, with a focus on the confinement effects that arise in these extremely narrow nanopores. The recent development of precision model systems, transformative experimental tools, and multiscale theories that have played enabling roles in advancing this frontier are reviewed. We also identify new knowledge gaps in our understanding of nanofluidic transport and provide an outlook for the future challenges and opportunities at this rapidly advancing frontier.

Capturing charge and size effects of ions at the graphene-electrolyte interface using polarizable force field simulations

We present a systematic investigation capturing the charge and size effects of ions interacting with a graphene surface using polarizable simulations. Our results utilizing the Drude polarizable force field (FF) for ions, water and graphene surfaces, show that the graphene parameters previously developed by us are able to accurately capture the dynamics at the electrolyte-graphene interface. For monovalent ions, with increasing size, the solvation shell plays a crucial role in controlling the ion-graphene interactions. Smaller monovalent ions directly interact with the graphene surface, while larger ions interact with the graphene surface via a well-formed solvation shell. For divalent ions, both interaction modes are observed. For the anion Cl^-, we observe direct interaction between the ions and the graphene surface. The anion-graphene interactions are strongly driven by the polarizability of the graphene surface. These effects are not captured in the absence of polarization by additive FF simulations. The present study underlines the importance of polarizability in capturing the interfacial phenomenon at the solid-solute interface.