2,055 research outputs found
Bayesian selection for coarse-grained models of liquid water
The necessity for accurate and computationally efficient representations of
water in atomistic simulations that can span biologically relevant timescales
has born the necessity of coarse-grained (CG) modeling. Despite numerous
advances, CG water models rely mostly on a-priori specified assumptions. How
these assumptions affect the model accuracy, efficiency, and in particular
transferability, has not been systematically investigated. Here we propose a
data driven, comparison and selection for CG water models through a
Hierarchical Bayesian framework. We examine CG water models that differ in
their level of coarse-graining, structure, and number of interaction sites. We
find that the importance of electrostatic interactions for the physical system
under consideration is a dominant criterion for the model selection. Multi-site
models are favored, unless the effects of water in electrostatic screening are
not relevant, in which case the single site model is preferred due to its
computational savings. The charge distribution is found to play an important
role in the multi-site model's accuracy while the flexibility of the
bonds/angles may only slightly improve the models. Furthermore, we find
significant variations in the computational cost of these models. We present a
data informed rationale for the selection of CG water models and provide
guidance for future water model designs
Molecular Structure and Modeling of Water-Air and Ice-Air Interfaces Monitored by Sum-Frequency Generation.
From a glass of water to glaciers in Antarctica, water-air and ice-air interfaces are abundant on Earth. Molecular-level structure and dynamics at these interfaces are key for understanding many chemical/physical/atmospheric processes including the slipperiness of ice surfaces, the surface tension of water, and evaporation/sublimation of water. Sum-frequency generation (SFG) spectroscopy is a powerful tool to probe the molecular-level structure of these interfaces because SFG can specifically probe the topmost interfacial water molecules separately from the bulk and is sensitive to molecular conformation. Nevertheless, experimental SFG has several limitations. For example, SFG cannot provide information on the depth of the interface and how the orientation of the molecules varies with distance from the surface. By combining the SFG spectroscopy with simulation techniques, one can directly compare the experimental data with the simulated SFG spectra, allowing us to unveil the molecular-level structure of water-air and ice-air interfaces. Here, we present an overview of the different simulation protocols available for SFG spectra calculations. We systematically compare the SFG spectra computed with different approaches, revealing the advantages and disadvantages of the different methods. Furthermore, we account for the findings through combined SFG experiments and simulations and provide future challenges for SFG experiments and simulations at different aqueous interfaces
Efficient minimization of multipole electrostatic potentials in torsion space
The development of models of macromolecular electrostatics capable of delivering improved fidelity to quantum mechanical calculations is an active field of research in computational chemistry. Most molecular force field development takes place in the context of models with full Cartesian coordinate degrees of freedom. Nevertheless, a number of macromolecular modeling programs use a reduced set of conformational variables limited to rotatable bonds. Efficient algorithms for minimizing the energies of macromolecular systems with torsional degrees of freedom have been developed with the assumption that all atom-atom interaction potentials are isotropic. We describe novel modifications to address the anisotropy of higher order multipole terms while retaining the efficiency of these approaches. In addition, we present a treatment for obtaining derivatives of atom-centered tensors with respect to torsional degrees of freedom. We apply these results to enable minimization of the Amoeba multipole electrostatics potential in a system with torsional degrees of freedom, and validate the correctness of the gradients by comparison to finite difference approximations. In the interest of enabling a complete model of electrostatics with implicit treatment of solvent-mediated effects, we also derive expressions for the derivative of solvent accessible surface area with respect to torsional degrees of freedom
Minimal Basis Iterative Stockholder: Atoms in Molecules for Force-Field Development
Atomic partial charges appear in the Coulomb term of many force-field models
and can be derived from electronic structure calculations with a myriad of
atoms-in-molecules (AIM) methods. More advanced models have also been proposed,
using the distributed nature of the electron cloud and atomic multipoles. In
this work, an electrostatic force field is defined through a concise
approximation of the electron density, for which the Coulomb interaction is
trivially evaluated. This approximate "pro-density" is expanded in a minimal
basis of atom-centered s-type Slater density functions, whose parameters are
optimized by minimizing the Kullback-Leibler divergence of the pro-density from
a reference electron density, e.g. obtained from an electronic structure
calculation. The proposed method, Minimal Basis Iterative Stockholder (MBIS),
is a variant of the Hirshfeld AIM method but it can also be used as a
density-fitting technique. An iterative algorithm to refine the pro-density is
easily implemented with a linear-scaling computational cost, enabling
applications to supramolecular systems. The benefits of the MBIS method are
demonstrated with systematic applications to molecular databases and extended
models of condensed phases. A comparison to 14 other AIM methods shows its
effectiveness when modeling electrostatic interactions. MBIS is also suitable
for rescaling atomic polarizabilities in the Tkatchenko-Sheffler scheme for
dispersion interactions.Comment: 61 pages, 12 figures, 2 table
Direct Wolf summation of a polarizable force field for silica
We extend the Wolf direct, pairwise r^(-1) summation method with spherical
truncation to dipolar interactions in silica. The Tangney-Scandolo interatomic
force field for silica takes regard of polarizable oxygen atoms whose dipole
moments are determined by iteration to a self-consistent solution. With Wolf
summation, the computational effort scales linearly in the system size and can
easily be distributed among many processors, thus making large-scale
simulations of dipoles possible. The details of the implementation are
explained. The approach is validated by estimations of the error term and
simulations of microstructural and thermodynamic properties of silica.Comment: See http://link.aip.org/link/?JCP/132/194109 - 8 pages, 6 figures.
Changes in v3: Copyright notice added, minor typographical changes. Changes
in v2: 1. Inserted Paragraph in Sec. IV B describing the limitations of the
TS potential. 2. We corrected transcription errors in Tab. II, and adjusted
the deviation percentages mentioned in Sec. IV B, first paragraph,
accordingl
- …