Finite domain simulations with adaptive boundaries: Accurate potentials and nonequilibrium movesets

Abstract

We extend the theory of hybrid explicit/implicit solvent models to include an explicit domain that grows and shrinks in response to a solute's evolving configuration. The goal of this model is to provide an appropriate but not excessive amount of solvent detail, and the inclusion of an adjustable boundary provides a significant computational advantage for solutes that explore a range of configurations. In addition to the theoretical development, a successful implementation of this method requires (1) an efficient moveset that propagates the boundary as a new coordinate of the system, and (2) an accurate continuum solvent model with parameters that are transferable to an explicit domain of any size. We address these challenges and develop boundary updates using Monte Carlo moves biased by nonequilibrium paths. We obtain the desired level of accuracy using a "decoupling interface" that we have previously shown to remove boundary artifacts common to hybrid solvent models. Using an uncharged, coarse-grained solvent model, we then study the efficiency of nonequilibrium paths that a simulation takes by quantifying the dissipation. In the spirit of optimization, we study this quantity over a range of simulation parameters.