RESUMO
Cassandra is an open source atomistic Monte Carlo software package that is effective in simulating the thermodynamic properties of fluids and solids. The different features and algorithms used in Cassandra are described, along with implementation details and theoretical underpinnings to various methods used. Benchmark and example calculations are shown, and information on how users can obtain the package and contribute to it are provided. © 2017 Wiley Periodicals, Inc.
RESUMO
Despite more than 40 years of research in condensed-matter physics, state-of-the-art approaches for simulating the radial distribution function (RDF) g(r) still rely on binning pair-separations into a histogram. Such methods suffer from undesirable properties, including subjectivity, high uncertainty, and slow rates of convergence. Moreover, such problems go undetected by the metrics often used to assess RDFs. To address these issues, we propose (I) a spectral Monte Carlo (SMC) quadrature method that yields g(r) as an analytical series expansion; and (II) a Sobolev norm that assesses the quality of RDFs by quantifying their fluctuations. Using the latter, we show that, relative to histogram-based approaches, SMC reduces by orders of magnitude both the noise in g(r) and the number of pair separations needed for acceptable convergence. Moreover, SMC reduces subjectivity and yields simple, differentiable formulas for the RDF, which are useful for tasks such as coarse-grained force-field calibration via iterative Boltzmann inversion.
RESUMO
Acceptance rules for reaction ensemble Monte Carlo (RxMC) simulations containing classically modeled atomistic degrees of freedom are derived for complex molecular systems where insertions and deletions are achieved gradually by utilizing the continuous fractional component (CFC) method. A self-consistent manner in which to utilize statistical mechanical data contained in ideal gas free energy parameters during RxMC moves is presented. The method is tested by applying it to two previously studied systems containing intramolecular degrees of freedom: the propene metathesis reaction and methyl-tert-butyl-ether (MTBE) synthesis. Quantitative agreement is found between the current results and those of Keil et al. (J. Chem. Phys. 2005, 122, 164705) for the propene metathesis reaction. Differences are observed between the equilibrium concentrations of the present study and those of Lísal et al. (AIChE J. 2000, 46, 866-875) for the MTBE reaction. It is shown that most of this difference can be attributed to an incorrect formulation of the Monte Carlo acceptance rule. Efficiency gains using CFC MC as opposed to single stage molecule insertions are presented.
RESUMO
The phase behavior of an embedded-charge model for lysozyme developed by Carlsson and co-workers (J. Phys. Chem. B 2001, 105, 9040) is investigated using grand canonical transition matrix Monte Carlo simulation. Within this model, protein-protein interactions are approximated through a combination of hard-sphere repulsion, isotropic hydrophobic attraction, and screened electrostatic interactions through a series of embedded point charges located at the positions of charged amino acid groups within lysozyme. Liquid-liquid phase diagrams are constructed for a wide range of solution conditions and compared with experimental data. Our results indicate that the model is generally capable of describing qualitative trends in the evolution of protein phase behavior with variation of pH and ionic strength. From a quantitative perspective, model estimates for both the change in critical temperature with variation of the solution conditions and the critical concentration do not agree with experimental results. We find the width of model coexistence curves to be independent of solution conditions and narrow relative to experimentally obtained phase envelopes. Connections between the value of the second virial coefficient evaluated at the critical temperature and the location of the liquid-liquid phase envelope are also examined.
