Non-Equilibrium Modeling of Concentration-Driven processes with Constant Chemical Potential Molecular Dynamics Simulations

Abstract

Conspectus Concentration-driven processes in solution, i.e., phenomena that are sustained by persistent concentration gradients, such as crystallization and surface adsorption, are fundamental chemical processes. Understanding such phenomena is crucial for countless applications, from pharmaceuticals to biotechnology. Molecular dynamics (MD), both in- and out-of-equilibrium, plays an essential role in the current understanding of concentration-driven processes. Computational costs, however, impose drastic limitations on the accessible scale of simulated systems, hampering the effective study of such phenomena. In particular, due to these size limitations, closed system MD of concentration-driven processes is affected by solution depletion/enrichment that unavoidably impacts the dynamics of the chemical phenomena under study. As a notable example, in simulations of crystallization from solution, the transfer of monomers between the liquid and crystal phases results in a gradual depletion/enrichment of solution concentration, altering the driving force for phase transition. In contrast, this effect is negligible in experiments, given the macroscopic size of the solution volume. Because of these limitations, accurate MD characterization of concentration-driven phenomena has proven to be a long-standing simulation challenge. While disparate equilibrium and nonequilibrium simulation strategies have been proposed to address the study of such processes, the methodologies are in continuous development. In this context, a novel simulation technique named constant chemical potential molecular dynamics (CμMD) was recently proposed. CμMD employs properly designed, concentration-dependent external forces that regulate the flux of solute species between selected subregions of the simulation volume. This enables simulations of systems under a constant chemical drive in an efficient and straightforward way. The CμMD scheme was originally applied to the case of crystal growth from solution and then extended to the simulation of various physicochemical processes, resulting in new variants of the method. This Account illustrates the CμMD method and the key advances enabled by it in the framework of in silico chemistry. We review results obtained in crystallization studies, where CμMD allows growth rate calculations and equilibrium shape predictions, and in adsorption studies, where adsorption thermodynamics on porous or solid surfaces was correctly characterized via CμMD. Furthermore, we will discuss the application of CμMD variants to simulate permeation through porous materials, solution separation, and nucleation upon fixed concentration gradients. While presenting the numerous applications of the method, we provide an original and comprehensive assessment of concentration-driven simulations using CμMD. To this end, we also shed light on the theoretical and technical foundations of CμMD, underlining the novelty and specificity of the method with respect to existing techniques while stressing its current limitations. Overall, the...

Graphene: From Solid Support for Nucleobase Assisted Self-Assemblies to Functional Material for DNA Sequencing

“…To this end, ionic conductivity and transport phenomena in nanopores and nanochannels is an active field of investigation with better representations required to capture the associated dynamics. Some recent developments to this end have been the use of nonequilibrium modeling to capture the concentration dependent dynamics. , Recent experimental studies have explored transverse detection of DNA strands using both graphene and MoS 2 nanopores, wherein nucleobase detection was observed in both the cases. Theoretically many groups have explored transverse current signatures of nucleobases trapped within a nanopore; however, a limitation with this work is that these DFT based calculations are performed using single layer nanopores. − A multilayer description must be adopted for a better description of nanopore interactions as evidenced by simulation studies. , We expect to see more investigations from both experimental and theoretical points of view aimed at investigating the transport of DNA molecules across nanopores.…”

Enhanced Sampling of Crystal Nucleation with Graph Representation Learnt Variables

“…For a more comprehensive summary of finite-size effects, refer to ref . To address this issue, various methods have been developed, including the introduction of constant chemical potential ensemble and analytical corrections to single-species , and solvated systems. − In future work, we hope to combine data-driven approaches with those from statistical mechanics to not just automate biasing variable discovery, as was done in this work, but also address finite-size effects.…”

Properties of aqueous electrolyte solutions at carbon electrodes: effects of concentration and surface charge on solution structure, ion clustering and thermodynamics in the electric double layer

“…In our simulations, the chemical potential of ions and water as a function of x in the steady state is constant. 31 Therefore, with knowledge of the ion molalities and electric potential in the EDL, eqn (11) can be rearranged to determine how the presence of the surface-where the density of ions and dielectric constant of the solution are changing compared with the bulk-affects the activity of ions as captured by g ion : −RT ln[g ion (m Na (x))g ion (m Cl (x))] = RT[ln m Na (x)m Cl (x)] + (2u − 1)Fj(x) (12) We label this quantity Dm E ion . Fig.…”

“…14,27,28 The asymmetric ordering of ions results in charge uctuations in this region-typically four-to-ve liquid layers deep-and a departure from descriptions of the EDL expected from the established mean-eld models described above. 9,29 Thanks to the adoption of the Constant Chemical Potential Molecular Dynamics (CmMD) method, 30,31 which maintains a constant thermodynamic driving force associated with ion adsorption, we were able to quantify the electric potential drop across the EDL and the excess chemical potential for ions at the solid-solution interface, 22,32 In CmMD, the use of an explicit molecular reservoir coupled to the model interface prevents any ion depletion in the bulk solution, which would otherwise occur in typical nite-sized MD simulations when ions adsorb at an interface.…”