Search

We develop a multiscale coarse-grain model of the NIST Monoclonal Antibody Reference Material 8671 (NISTmAb) to enable systematic computational investigations of high-concentration physical instabilities such as phase separation, clustering, and aggregation. Our multiscale coarse-graining strategy captures atomic-resolution interactions with a computational approach that is orders of magnitude more efficient than atomistic models, assuming the biomolecule can be decomposed into one or more rigid bodies with known, fixed structures. This method reduces interactions between tens of thousands of atoms to a single anisotropic interaction site. The anisotropic interaction between unique pairs of rigid bodies is precomputed over a discrete set of relative orientations and stored, allowing interactions between arbitrarily oriented rigid bodies to be interpolated from the precomputed table during coarse-grained Monte Carlo simulations. We present this approach for lysozyme and lactoferrin as a single rigid body and for the NISTmAb as three rigid bodies bound by a flexible hinge with an implicit solvent model. This coarse-graining strategy predicts experimentally measured radius of gyration and second osmotic virial coefficient data, enabling routine Monte Carlo simulation of medically relevant concentrations of interacting proteins while retaining atomistic detail. All methodologies used in this work are available in the open-source software Free Energy and Advanced Sampling Simulation Toolkit. [ABSTRACT FROM AUTHOR]

FEASST is an open-source Monte Carlo software package for particle-based simulations. This software, which was released in 2017, has been used to study phase equilibrium, self-assembly, aggregation or gelation in biological materials, colloids, polymers, ionic liquids, and adsorption in porous networks. We highlight some of the unique features available in FEASST, such as flat-histogram grand canonical ensemble, Gibbs ensemble, and Mayer-sampling simulations with support for anisotropic models and parallelization with flat-histogram and prefetching. We also discuss how the challenges of supporting a variety of Monte Carlo algorithms were overcome by an object-oriented design. This also allows others to extend classes, which improves software interoperability, as inspired by LAMMPS classes and user packages. This article describes version 0.25.1 with benchmarks, compilation instructions, and introductory tutorials for running, restarting, and testing simulations, user guidelines, software design strategies, alternative interfaces, and the test-driven development strategy. [ABSTRACT FROM AUTHOR]

Theories of small systems play an important role in the fundamental understanding of finite size effects in statistical mechanics, as well as the validation of molecular simulation results as no computer can simulate fluids in the thermodynamic limit. Previously, a shell particle was included in the isothermal–isobaric ensemble in order to resolve an ambiguity in the resulting partition function. The shell particle removed either redundant volume states or redundant translational degrees of freedom of the system and yielded quantitative differences from traditional simulations in this ensemble. In this work, we investigate the effect of including a shell particle in the canonical, grand canonical, and Gibbs ensembles. For systems comprised of a pure component ideal gas, analytical expressions for various thermodynamic properties are obtained. We also derive the Metropolis Monte Carlo simulation acceptance criteria for these ensembles with shell particles, and the results of the simulations of an ideal gas are in excellent agreement with the theoretical predictions. The system size dependence of various important ensemble averages is also analyzed. [ABSTRACT FROM AUTHOR]

We introduce Gaussian Process Regression (GPR) as an enhanced method of thermodynamic extrapolation and interpolation. The heteroscedastic GPR models that we introduce automatically weight provided information by its estimated uncertainty, allowing for the incorporation of highly uncertain, high-order derivative information. By the linearity of the derivative operator, GPR models naturally handle derivative information and, with appropriate likelihood models that incorporate heterogeneous uncertainties, are able to identify estimates of functions for which the provided observations and derivatives are inconsistent due to the sampling bias that is common in molecular simulations. Since we utilize kernels that form complete bases on the function space to be learned, the estimated uncertainty in the model takes into account that of the functional form itself, in contrast to polynomial interpolation, which explicitly assumes the functional form to be fixed. We apply GPR models to a variety of data sources and assess various active learning strategies, identifying when specific options will be most useful. Our active-learning data collection based on GPR models incorporating derivative information is finally applied to tracing vapor–liquid equilibrium for a single-component Lennard-Jones fluid, which we show represents a powerful generalization to previous extrapolation strategies and Gibbs–Duhem integration. A suite of tools implementing these methods is provided at https://github.com/usnistgov/thermo-extrap. [ABSTRACT FROM AUTHOR]

We describe a method for deriving surface functionalization patterns for colloidal systems that can induce self-assembly into any chosen periodic symmetry at a planar interface. The result is a sequence of letters, s ∈ {A,T,C,G}, or a gene, that describes the perimeter of the colloidal object and programs its self-assembly. This represents a genome that is finite and can be exhaustively enumerated. These genes derive from symmetry, which may be topologically represented by two-dimensional parabolic orbifolds; since these orbifolds are surfaces that may be derived from first principles, this represents an ab initio route to colloid functionality. The genes are human readable and can be employed to easily design colloidal units. We employ a biological (genetic) analogy to demonstrate this and illustrate their connection to the designs of Maurits Cornelis (M. C.) Escher. [ABSTRACT FROM AUTHOR]

Variational autoencoders (VAEs) are rapidly gaining popularity within molecular simulation for discovering low-dimensional, or latent, representations, which are critical for both analyzing and accelerating simulations. However, it remains unclear how the information a VAE learns is connected to its probabilistic structure and, in turn, its loss function. Previous studies have focused on feature engineering, ad hoc modifications to loss functions, or adjustment of the prior to enforce desirable latent space properties. By applying effectively arbitrarily flexible priors via normalizing flows, we focus instead on how adjusting the structure of the decoding model impacts the learned latent coordinate. We systematically adjust the power and flexibility of the decoding distribution, observing that this has a significant impact on the structure of the latent space as measured by a suite of metrics developed in this work. By also varying weights on separate terms within each VAE loss function, we show that the level of detail encoded can be further tuned. This provides practical guidance for utilizing VAEs to extract varying resolutions of low-dimensional information from molecular dynamics and Monte Carlo simulations. [ABSTRACT FROM AUTHOR]

Thermodynamic extrapolation has previously been used to predict arbitrary structural observables in molecular simulations at temperatures (or relative chemical potentials in open-system mixtures) different from those at which the simulation was performed. This greatly reduces the computational cost in mapping out phase and structural transitions. In this work, we explore the limitations and accuracy of thermodynamic extrapolation applied to water, where qualitative shifts from anomalous to simple-fluid-like behavior are manifested through shifts in the liquid structure that occur as a function of both temperature and density. We present formulas for extrapolating in volume for canonical ensembles and demonstrate that linear extrapolations of water's structural properties are only accurate over a limited density range. On the other hand, linear extrapolation in temperature can be accurate across the entire liquid state. We contrast these extrapolations with classical perturbation theory techniques, which are more conservative and slowly converging. Indeed, we show that such behavior is expected by demonstrating exact relationships between extrapolation of free energies and well-known techniques to predict free energy differences. An ideal gas in an external field is also studied to more clearly explain these results for a toy system with fully analytical solutions. We also present a recursive interpolation strategy for predicting arbitrary structural properties of molecular fluids over a predefined range of state conditions, demonstrating its success in mapping qualitative shifts in water structure with density. [ABSTRACT FROM AUTHOR]

While ionic liquids have promising applications as industrial solvents, predicting their fluid phase properties and coexistence remains a challenge. Grand canonical Monte Carlo simulation is an effective method for such predictions, but equilibration is hampered by the apparent requirement to insert and delete neutral sets of ions simultaneously in order to maintain charge neutrality. For relatively high densities and low temperatures, previously developed methods have been shown to be essential in improving equilibration by gradual insertion and deletion of these neutral sets of ions. We introduce an expanded ensemble approach which may be used in conjunction with these existing methods to further improve efficiency. Individual ions are inserted or deleted in one Monte Carlo trial rather than simultaneous insertion/deletion of neutral sets. We show how charge neutrality is maintained and show rigorous quantitative agreement between the conventional and the proposed expanded ensemble approaches, but with up to an order of magnitude increase in efficiency at high densities. The expanded ensemble approach is also more straightforward to implement than simultaneous insertion/deletion of neutral sets, and its implementation is demonstrated within open source software. [ABSTRACT FROM AUTHOR]

We present analytic expressions for the second virial coefficient of hard particles in slit and cylindrical pores with arbitrary square-well fluid-fluid and fluid-solid interactions, number of fluid components, and pore sizes. We derive multiple methods to predict the isotherms in confined materials based on the virial expansion and find that the thermodynamic and structural properties calculated from the virial expansion to second order in the dilute limit display excellent qualitative agreement with previous simulation results. [ABSTRACT FROM AUTHOR]

The theoretical framework to evaluate small-angle scattering (SAS) profiles for multi-component macromolecular solutions is re-examined from the standpoint of molecular simulations in the grand-canonical ensemble, where the chemical potentials of all species in solution are fixed. This statistical mechanical ensemble resembles more closely scattering experiments, capturing concentration fluctuations that arise from the exchange of molecules between the scattering volume and the bulk solution. The resulting grand-canonical expression relates scattering intensities to the different intra- and intermolecular pair distribution functions, as well as to the distribution of molecular concentrations on the scattering volume. This formulation represents a generalized expression that encompasses most of the existing methods to evaluate SAS profiles from molecular simulations. The grand-canonical SAS methodology is probed for a series of different implicit-solvent, homogeneous systems at conditions ranging from dilute to concentrated. These systems consist of spherical colloids, dumbbell particles, and highly flexible polymer chains. Comparison of the resulting SAS curves against classical methodologies based on either theoretical approaches or canonical simulations (i.e., at a fixed number of molecules) shows equivalence between the different scattering intensities so long as interactions between molecules are net repulsive or weakly attractive. On the other hand, for strongly attractive interactions, grand-canonical SAS profiles deviate in the low- and intermediate-q range from those calculated in a canonical ensemble. Such differences are due to the distribution of molecules becoming asymmetric, which yields a higher contribution from configurations with molecular concentrations larger than the nominal value. Additionally, for flexible systems, explicit discrimination between intra- and inter-molecular SAS contributions permits the implementation of model-free, structural analysis such as Guinier’s plots at high molecular concentrations, beyond what the traditional limits are for such analysis. [ABSTRACT FROM AUTHOR]

We describe a methodology for extrapolating the structural properties of multicomponent fluids from one thermodynamic state to another. These properties generally include features of a system that may be computed from an individual configuration such as radial distribution functions, cluster size distributions, or a polymer’s radius of gyration. This approach is based on the principle of using fluctuations in a system’s extensive thermodynamic variables, such as energy, to construct an appropriate Taylor series expansion for these structural properties in terms of intensive conjugate variables, such as temperature. Thus, one may extrapolate these properties from one state to another when the series is truncated to some finite order. We demonstrate this extrapolation for simple and coarse-grained fluids in both the canonical and grand canonical ensembles, in terms of both temperatures and the chemical potentials of different components. The results show that this method is able to reasonably approximate structural properties of such fluids over a broad range of conditions. Consequently, this methodology may be employed to increase the computational efficiency of molecular simulations used to measure the structural properties of certain fluid systems, especially those used in high-throughput or data-driven investigations. [ABSTRACT FROM AUTHOR]

We derive an approach for extrapolating the free energy landscape of multicomponent systems in the grand canonical ensemble, obtained from flat-histogram Monte Carlo simulations, from one set of temperature and chemical potentials to another. This is accomplished by expanding the landscape in a Taylor series at each value of the order parameter which defines its macrostate phase space. The coefficients in each Taylor polynomial are known exactly from fluctuation formulas, which may be computed by measuring the appropriate moments of extensive variables that fluctuate in this ensemble. Here we derive the expressions necessary to define these coefficients up to arbitrary order. In principle, this enables a single flat-histogram simulation to provide complete thermodynamic information over a broad range of temperatures and chemical potentials. Using this, we also showhowto combine a small number of simulations, each performed at different conditions, in a thermodynamically consistent fashion to accurately compute properties at arbitrary temperatures and chemical potentials. This method may significantly increase the computational efficiency of biased grand canonical Monte Carlo simulations, especially for multicomponent mixtures. Although approximate, this approach is amenable to high-throughput and data-intensive investigations where it is preferable to have a large quantity of reasonably accurate simulation data, rather than a smaller amount with a higher accuracy. [ABSTRACT FROM AUTHOR]

Virial coefficients are predicted over a large range of both temperatures and model parameter values (i.e., alchemical transformation) from an individual Mayer-sampling Monte Carlo simulation by statistical mechanical extrapolation with minimal increase in computational cost.With this extrapolation method, a Mayer-sampling Monte Carlo simulation of the SPC/E (extended simple point charge) water model quantitatively predicted the second virial coefficient as a continuous function spanning over four orders of magnitude in value and over three orders of magnitude in temperature with less than a 2% deviation. In addition, the same simulation predicted the second virial coefficient if the site charges were scaled by a constant factor, from an increase of 40% down to zero charge. This method is also shown to perform well for the third virial coefficient and the exponential parameter [ABSTRACT FROM AUTHOR]

Proteins aggregate and precipitate from high concentration solutions in a wide variety of problems of natural and technological interest. Consequently, there is a broad interest in developing new ways to model the thermodynamic and kinetic aspects of protein stability in these crowded cellular or solution environments. We use a coarse-grained modeling approach to study the effects of different crowding agents on the conformational equilibria of proteins and the thermodynamic phase behavior of their solutions. At low to moderate protein concentrations, we find that crowding species can either stabilize or destabilize the native state, depending on the strength of their attractive interaction with the proteins. At high protein concentrations, crowders tend to stabilize the native state due to excluded volume effects, irrespective of the strength of the crowder-protein attraction. Crowding agents reduce the tendency of protein solutions to undergo a liquid-liquid phase separation driven by strong protein-protein attractions. The aforementioned equilibrium trends represent, to our knowledge, the first simulation predictions for how the properties of crowding species impact the global thermodynamic stability of proteins and their solutions. [DOI: 10.1115/1.3127259]

We derive a method for extrapolating the grand canonical free energy landscape of a multicomponent fluid system from one temperature to another. Previously, we introduced this statistical mechanical framework for the case where kinetic energy contributions to the classical partition function were neglected for simplicity [N. A. Mahynski et al., J. Chem. Phys. 146, 074101 (2017)]. Here, we generalize the derivation to admit these contributions in order to explicitly illustrate the differences that result. Specifically, we show how factoring out kinetic energy effects a priori, in order to consider only the configurational partition function, leads to simpler mathematical expressions that tend to produce more accurate extrapolations than when these effects are included. We demonstrate this by comparing and contrasting these two approaches for the simple cases of an ideal gas and a non-ideal, square-well fluid. [ABSTRACT FROM AUTHOR]

We present a method for predicting the free energy landscape of fluids at low temperatures from flat-histogram grand canonical Monte Carlo simulations performed at higher ones. We illustrate our approach for both pure and multicomponent systems using two different sampling methods as a demonstration. This allows us to predict the thermodynamic behavior of systems which undergo both first order and continuous phase transitions upon cooling using simulations performed only at higher temperatures. After surveying a variety of different systems, we identify a range of temperature differences over which the extrapolation of high temperature simulations tends to quantitatively predict the thermodynamic properties of fluids at lower ones. Beyond this range, extrapolation still provides a reasonably well-informed estimate of the free energy landscape; this prediction then requires less computational effort to refine with an additional simulation at the desired temperature than reconstruction of the surface without any initial estimate. In either case, this method significantly increases the computational efficiency of these flat-histogram methods when investigating thermodynamic properties of fluids over a wide range of temperatures. For example, we demonstrate how a binary fluid phase diagram may be quantitatively predicted for many temperatures using only information obtained from a single supercritical state. [ABSTRACT FROM AUTHOR]

We demonstrate an extensible flat-histogram Monte Carlo simulation methodology for studying the adsorption of multicomponent fluids in flexible porous solids. This methodology allows us to easily obtain the complete free energy landscape for the confined fluid-solid system in equilibrium with a bulk fluid of any arbitrary composition. We use this approach to study the adsorption of a prototypical coarse-grained binary fluid in "Hookean" solids, where the free energy of the solid may be described as a simple spring. However, our approach is fully extensible to solids with arbitrarily complex free energy profiles. We demonstrate that by tuning the fluid-solid interaction ranges, the inhomogeneous fluid structure inside the pore can give rise to enhanced selective capture of a larger species through cooperative adsorption with a smaller one. The maximum enhancement in selectivity is observed at low to intermediate pressures and is especially pronounced when the larger species is very dilute in the bulk. This suggest a mechanism by which the selective capture of a minor component from a bulk fluid may be enhanced. [ABSTRACT FROM AUTHOR]

We investigate the thermodynamic properties of various super-critical model adsorptive systems with different fluid–solid attractive strengths using the confined density virial expansion, with coefficients calculated using the Mayer-sampling Monte Carlo method up to fifth order. We find that the virial expansion converges for adsorptive systems over a density range corresponding approximately to the film-formation regime. Beyond this regime, higher-order effects become increasingly important. The virial expansion of the density profile is also investigated. It is determined that this expansion gives insight into the structure associated with adsorption. We also find that weakly attractive systems have a more negative second virial coefficient than strongly attractive systems. This runs counter to the usual interpretation of bulk fluid virial coefficients. This is due to the infinite-dilution limit being very different for adsorbed fluids compared to bulk fluids. [ABSTRACT FROM AUTHOR]